La régénération naturelle des peuplements de chênes verts (Quercus ilex) des systèmes agroforestiers de « dehesas » pose problème dans leur aire de distribution principale. La cause invoquée est l’absence de strate arbustive protégeant les semis de l’abroutissement et de la chaleur. La thèse, à travers la méthode régressive, propose d’explorer en détail cette piste. Le Parc Naturel de la Sierra de Grazalema (Andalousie, Espagne), en périphérie de la distribution des dehesas et inexploré sur ces questions, a été le lieu de nos investigations. La régénération naturelle a été analysée comme un concept évoluant dans le temps et les champs disciplinaires : c’est une technique naturalisée. A l’analyse de documents historiques, allant de 1765 à nos jours, s’ajoute un travail de terrain comprenant 71 transects mesurant les strates arbustive et arborée sur cinq terres publiques et trois terres privées. Une répartition en classes de diamètre pour la strate arborée a été effectuée et des variogrammes ont été utilisés pour analyser le patron de semis à plusieurs échelles. La distribution des classes de diamètre, suivant une distribution en « J inversé », indique une réussite de la régénération avec cependant une forte hétérogénéité entre les transect et le statut foncier. Un nombre important de jeunes arbres a été trouvés : 4516,65 semis/ha, 986 fourrés/ha et 543 gaules/ha. Les variogrammes indiquent des patrons de régénération à plusieurs échelles, une dizaine de mètres, 45 à 90 mètres, 100 à 200 mètres et proche du kilomètre. Ces résultats, très différents des études précédentes et venant compléter des trajectoires historiques, indiquent qu’il n’y a pas une forme de régénération mais plusieurs en fonction des époques mais aussi des milieux. Il y a des cas de régénération sans un passage par une « matorralisation », des cas s’inscrivant dans des dynamiques plus proches de séries (régressives ou progressives) et enfin des cas de régénération sensu stricto. / The regeneration of Holm oak (Quercus ilex) populations remains problematic in most dehesa agroforestry systems, allegedly due to the absence of a shrub layer protecting seedlings from grazing and heat. This thesis explores the latter claim using the regressive method. Our study focuses on the Sierra de Grazalema Natural Park, a site located at the periphery of dehesa regions in Andalusia, Spain and not previously analyzed in this manner. Regeneration was interpreted as a naturalized technique that evolves over time and within disciplinary fields. Our study features an analysis of historical documents, from 1765 to the present day, as well as fieldwork including 71 transects measuring shrub and tree layers on five public lands and three private lands. The tree layer was examined using a distribution into diameter classes and variograms analyzed seeding patterns on several scales. The reverse J-shaped diameter distribution indicates successful regeneration, but with a strong heterogeneity both among transects and due to land status. A significant number of young trees were found: 4516.65 seedlings / ha, 986 thickets / ha and 543 saplings / ha. Variograms indicate regeneration patterns present at several scales; at about ten meters, at 45 to 90 meters, at 100 to 200 meters and near one kilometer. These results, very different from previous studies, complement historical trajectories, and indicate that there is not one form of regeneration but several, related to both historical eras and environments. Patterns of regeneration include cases without a matorral shrublands and woodlands stage, cases belonging to a regressive or progressive series, and finally, cases of regeneration sensu stricto.

Résumé: Tel que le notent si bien Oliver et al. (1990), la plupart des méthodes d’interpolation classiques sont basées sur les modèles des méthodes numériques différents des modèles probabilistes de variation spatiale. Les distributions des variables des données spatiales se comportent plus comme des variables aléatoires et la théorie des variables régionalisées comporte un ensemble de méthodes stochastiques pour les analyser. Le krigeage et le cokrigeage, basés sur cette théorie, expriment les variations spatiales de variables en termes de variogrammes et ils minimisent les erreurs de prédiction qui sont elles-mêmes estimées. Appliquées aux données des variables pédologiques des champs de la Station de recherches de Lennoxville, les méthodes ou techniques d’interpolation de la géostatistique: les krigeages et cokrigeages ordinaires, universels et par bloc, ainsi dénommées, sont explorées et comparées en vue d’en identifier la méthode appropriée intégrable dans un système d’information géographique, pour obtenir l’optimum de la modélisation, la présentation et celui de l’analyse spatiale. Les valeurs krigées et cokrigées des techniques choisies sont obtenues à l’aide de l’algorithme Cokri (Marcotte, 1991, 1993) et sont comparées aux valeurs initiales dont elles en sont les estimations. Des tests de similarité appliqués à ces données et résultats attestent la similarité de leurs distributions respectives. En d’autres termes, les types de krigeage et de cokrigeage examinés bénéficient de la même homogénéité des données qui est vérifiée par des outils de comparaison visuels et des traitements géostatistiques. De plus l’utilisation de ces méthodes donne des erreurs minimes proches de la réalité. Ainsi, selon la nature spécifique des données spatiales, en l’occurence les données homogènes sans structure spatiale particulière cachée, un estimateur linéaire simple ou la moyenne arithmétique de ces données, suffit pour représenter les variables étudiées.||Abstract: As Oliver et al. (1990) pointed it out so well, most of iraditional methods of interpolation are based on numerical method models as distinct from stochastic models of spatial variation. Spatially distributed data behave more like random variables, however, and regionalized variable theory provides a set of stochastic methods for analyzing them. Kriging and cokriging, based on that theory, depend on expressing spatial variation of the property in terms of variograms, and they minimize the prediction erros which are themselves estimated. Used on pedological parameters data, oridinary, universai and block krigings and cokrigings, geostatistical interpolation methods or techniques, are expiored and compared in order to identify the appropriate method which can be integrated into a Geographic Information System to provide optimum modelting, presentation and spatial analysis. Kriged and cokriged values are obtained by Cokri algorithm (Marcotte, 1991, 1993), and are compared to estimated initial data. Similarity tests applied to these data and values show that their respective distributions are similar. In othe words, the examined kriging and cokriging types perfom equally well. They ail take advantage of the data homogeneity, verified by visual comparison and geostatiatical tools or techniques, for providing doser minimum error terms. Thus, according to the special nature of spatial data, namely the homogeneous ones without a hiden particular spatial structure, a simple linear estimator, such as the data arithmetical mean, is sufficient for representing die variables under consideration.

Depuis près de quinze ans, l'Andra mène un important programme scientifique pour acquérir des connaissances en vue d'un éventuel stockage de déchets radioactifs en milieu géologique profond. L'outil de recherche le plus important est le laboratoire souterrain de Meuse/Haute-Marne, implanté à 490 mètres de profondeur au cœur d'une formation d'argilites datées du Callovo-Oxfordien. Cette thèse a pour objectif de tirer le meilleur parti possible des données diagraphiques haute définition et d'optimiser leur étude pour préciser et caractériser les variations et la variabilité spatio-temporelle des argilites du site de Meuse/Haute-Marne dans une zone de 250 km2 autour du laboratoire actuel. Cette étude revient à répondre à la question " quel est le temps réellement enregistré par des sédiments ? ", une problématique qui peut avoir des applications importantes dans les mondes industriel et scientifique. L'étude des variations spatiales peut être menée à différentes échelles. Une première échelle de variabilité est donnée par les diagraphies conventionnelles, qui mettent en évidence des variations latérales d'épaisseur entre les forages Andra. Une deuxième échelle de variabilité est donnée par l'analyse détaillée d'un outil à haute résolution : le FMI®, sur lequel des alternances régulières de niveaux argileux et plus carbonatés sont clairement identifiées. Des traitements géostatistiques ont été effectués pour affiner l'étude de cette cyclicité. Trois niveaux de cyclicité ont pu être quantifiés par une analyse variographique de l'ensemble des données diagraphiques et des données FMI®. Ces dernières ont permis de distinguer les périodes les plus courtes. Une correspondance entre les périodes mises en évidence et les cycles orbitaux de Milankovitch (excentricité, obliquité et précession, dont les durées sont connues, a pu être établie. Elle constitue la première étape vers le passage à un repère géochronologique. La caractérisation des variations fines du taux de sédimentation a nécessité l'utilisation d'outils géostatistiques adaptés tels que l'analyse krigeante factorielle. L'estimation locale des périodes verticales a montré les variations fines du taux de sédimentation, repéré les niveaux condensés, et ainsi permis de quantifier ces variations et de compter le nombre exact de cycles dans un intervalle donné.<br/> Le passage à un repère géochronologique a d'abord permis d'estimer des durées et en particulier la durée d'unités biostratigraphiques dans des étages (Callovien et Oxfordien) où les durées absolues ne sont quasiment pas connues. Enfin, l'étude des variations du taux moyen de sédimentation dans un repère géochronologique a montré les évolutions verticales et latérales du taux de sédimentation à l'échelle du site de Meuse/Haute-Marne et le synchronisme de la plupart des niveaux condensés. Elle a infirmé certaines hypothèses concernant la présence de hiatus et en a validé d'autres, et dans ce cas a abouti à l'estimation de la durée de ces hiatus.

Actuellement, les technologies sans fil sont omniprésentes dans les lieux de vie, de travail ou detransport de la population. Les recherches en dosimétrie numérique non-ionisante ont connu desavancées considérables notamment grâce à la capacité de calcul croissante des ordinateurs. Cesdernières années, la prise en compte de la variabilité dans la dosimétrie est devenue un enjeumajeur. Les sources de variabilité sont nombreuses, parmi elles, les conditions d'exposition aurayonnement électromagnétique peuvent induire des doses absorbées très différentes.La modélisation de canal de propagation en environnement intérieur a permis d'avoir uneconnaissance précise des conditions d'exposition d'un corps humain plongé dans cet environnement.Les modèles stochastiques montrent le caractère statistique de l'exposition des corpshumains. Cette thèse a pour but de développer une méthode statistique du calcul de la doseabsorbée par le corps et d'adapter le modèle de canal stochastique à la dosimétrie.L'étude statistique de l'exposition a révélé la nécessité d'obtenir les valeurs de Débit d'AbsorptionSpécifique corps entier dans le cas d'expositions à une onde plane pour tous les anglesd'incidence possibles. Compte tenu des temps de calcul particulièrement long en dosimétrie, uneméthode d'interpolation efficace, le krigeage, des valeurs de Débit d'Absorption Spécifique a étémise en œuvre. La méthode du krigeage permet alors d'obtenir les valeurs de Débit d'AbsorptionSpécifique pour tous les angles d'incidences et donc de calculer les deux premiers momentsdu Débit d'Absorption Spécifique.L'analyse de sensibilité aux paramètres du canal des moments du Débit d'Absorption Spécifiquea permis de connaître l'impact de chacun de ces paramètres. Le modèle de canal a pu être simplifié et donc adapté à la dosimétrie tout en quantifiant l'erreur d'approximation qu'implique cette simplification.Cette thèse répond à l'enjeu de la prise en compte de la variabilité en dosimétrie dans un environnementcomplexe et les outils mis en œuvre ouvrent la voie de l'étude de l'exposition danstous les environnements complexes. / info:eu-repo/semantics/nonPublished

Cette thèse vise à développer des algorithmes de simulation conditionnelle de fonctions aléatoires à lois bivariables isofactorielles, c'est-à-dire telles qu'il existe une base de fonctions (facteurs) sans corrélations spatiales croisées. Les lois isofactorielles réalisent un compromis entre deux exigences souvent difficiles à concilier: d'une part, elles forment une vaste classe de modèles et permettent de s'adapter à la description de nombreux phénomènes; d'autre part, l'inférence statistique du modèle repose sur un faible nombre de paramètres. La première partie de la thèse met en relief les limitations et les approximations commises par une technique réputée "passe-partout": l'algorithme séquentiel, qui consiste à simuler les sites de l'espace les uns après les autres dans un ordre aléatoire, en estimant en chaque site la distribution de probabilité conditionnelle de la valeur inconnue par un krigeage d'indicatrices ou par un krigeage disjonctif. La seconde partie propose des modèles nouveaux et des algorithmes adaptés à leur simulation et au conditionnement à des données expérimentales. Plusieurs outils structuraux sont introduits et étudiés pour mener l'inférence des lois bivariables, notamment le madogramme, les variogrammes d'ordre inférieur à deux et les variogrammes d'indicatrices. Les concepts et méthodes sont finalement appliqués à un ensemble de données minières (accumulations en or et argent dans un gisement chilien en veine) caractérisées par une très forte dissymétrie des histogrammes.

Dans le cadre de recherches menées par l'ANDRA, qui visent à évaluer la faisabilité de créer un stockage de déchets radioactifs dans la formation argileuse du Callovo Oxfordien, cette thèse s'intéresse à deux étapes indispensables pour obtenir une image 3D du sous sol. La première étape concerne la cartographie de la géométrie de la formation. Celle-ci est déterminée à partir de temps pointés le long de profils sismiques, orientés suivant deux directions principales. Aux intersections, les valeurs enregistrées par deux profils sont différentes. Le faible nombre d'intersections ne permettant pas une étude détaillée du comportement spatial de ces différences, elles sont analysées par le biais du pseudo variogramme croisé. L'analyse variographique permet de choisir un modèle cohérent pour l'ensemble des profils. À partir de ce modèle une estimation du temps corrigée est présentée. La deuxième étape concerne l'étude de l'impédance, variable liée aux propriétés intrinsèques de la roche. Deux types de mesure sont disponibles. D'une part l'impédance sismique résultant de l'inversion géophysique de l'amplitude enregistrée le long de profils de sismique réflexion est largement étendue dans la zone d'intérêt. D'autre part, l'impédance log, résultant du produit de la densité de la roche et de la vitesse des ondes P, n'est connue uniquement dans quelques puits. Le comportement spatial bivariable, le long de la direction verticale, a permis de formuler l'hypothèse suivante : chaque mesure d'impédance est décomposée en une somme d'une variable Z (commune aux deux mesures) et d'un résidu différent pour chaque mesure. Sous certaines hypothèses supplémentaires, il est possible d'estimer la composante commune sur l'ensemble des sections sismiques. En complément à ces études, d'autres points ont été évoqués, par exemple la correspondance entre procédés géophysiques et géostatistiques, ou la prise en compte d'incertitude sur les amplitudes.

L’analyse du signal de parole offre un moyen privilégié pour l’évaluation clinique de la qualité de la voix en vue d’un diagnostique et d’une documentation quantitative des pathologies du larynx. Une analyse acoustique présente plusieurs avantages. En plus du coût peu élevé du système d’évaluation, elle est simple à mettre en œuvre, présente un caractère non invasif et documente quantitativement le degré d’enrouement perçu par le clinicien. L’objectif de l’analyse du signal acoustique est d’extraire des indices pertinents permettant de déterminer les caractéristiques de la voix afin de renseigner sur l’état du larynx du locuteur. Dans ce contexte, plusieurs indices acoustiques sont utilisés pour caractériser la parole produite par des locuteurs dysphoniques. Un nombre de ces indices reflète la déviation du signal de parole voisée par rapport à la périodicité parfaite. Les causes de ces dyspériodicités sont diverses : vibrations non modales des cordes vocales, bruit de modulation comprenant les variations cycle à cycle de la durée de cycle (jitter) et de l’amplitude (shimmer) dues aux perturbations externes, bruit additif dû à une turbulence excessive, etc. Les indices acoustiques des dyspériodicités vocales sont souvent obtenus à partir de fragments stationnaires de voyelles soutenues. La raison en est que les voyelles sans attaques et déclins sont faciles à analyser parce que les hypothèses de cyclicité et de stationnarité utilisées par les méthodes d’analyse sont valables pour beaucoup de locuteurs. En effet, les voyelles soutenues peuvent être supposées avec une bonne précision comme produites en maintenant invariant dans le temps les caractéristiques de la source vocale, du conduit vocal et des articulateurs et, donc, les paramètres des perturbations et du bruit sont facilement calculés pour les voyelles soutenues. La plupart des cliniciens considèrent la parole continue plus informative que les voyelles soutenues. Plusieurs arguments en faveur de l’analyse de la parole continue peuvent être avancés. La vibration des cordes vocales doit commuter continuellement pour donner lieu à l’apparition ou l’extinction du voisement, le voisement doit être maintenu alors que l’impédance supra-glottique change constamment, plus particulièrement durant les obstruents, et le larynx doit descendre et monter continuellement. Le larynx fonctionne donc dans des conditions non stationnaires et très variables. La parole continue contient donc les caractéristiques dynamiques de la source de la voix et du conduit vocal tels que l’attaque et le déclin et les variations dans la fréquence fondamentale et l’amplitude. De même, il semble que les locuteurs compensent moins pour leurs problèmes de voix lors de la production de la parole continue que lorsqu’ils produisent des voyelles soutenues. Les méthodes de traitement du signal de parole pour la caractérisation des troubles de la voix sont basées sur les hypothèses de stationnarité et de périodicité locales qui conduisent à des approches heuristiques permettant de détecter et d’isoler les périodes fondamentales ou les harmoniques spectrales. Comme conséquence, des erreurs d’insertion ou d’omission se produisent souvent lors de l’analyse des signaux fortement perturbés. Ces erreurs biaisent numériquement les indices acoustiques. Les mesures objectives ne sont donc fiables que lorsque l’analyse est effectuée sur des voyelles soutenues produites par des locuteurs faiblement ou modérément enroués. La précision des méthodes d’analyse dans le cadre de l’estimation des dyspériodicités vocales est une caractéristique fondamentale. En effet, les perturbations cycle à cycle des durées de cycle peuvent être inférieures à 1 % alors que les perturbations cycle à cycle de l’amplitude peuvent être inférieures à 10 %. Donc, les méthodes de traitement doivent être appliquées avec précaution vis à vis de la précision de mesure pour ne pas biaiser les valeurs numériques. Le bruit de quantification par exemple peut affecter la précision d’estimation s’il n’est pas pris en considération. Cette thèse se focalise sur le développement de méthodes robustes d’analyse acoustiques du signal de parole continue en vue de la caractérisation des troubles de la voix, la définition et l’évaluation d’indices acoustiques pour quantifier les dyspériodicités vocales et le développement d’une méthode d’évaluation perceptive fiable permettant de mesurer et comparer les performances des différentes méthodes d’analyse développées. Les méthodes développées sont de deux types : des méthodes d’analyse par bloc qui opèrent sur des trames de courte durée du signal et des méthodes adaptatives qui permettent d’analyser le signal à chaque échantillon en tenant compte de son caractère non stationnaire. Deux indices acoustiques sont utilisés pour quantifier les dyspériodicités vocales dans le signal de parole. Le premier indice, utilisé conventionnellement dans le cadre de l’évaluation objective de la qualité de la voix, est le rapport signal à dyspériodicité global. La valeur numérique de l’indice global est principalement déterminée par les segments vocaliques et donc il peut masquer certains évènements locaux, notamment dans le cas de l’analyse de la parole continue. Le second indice proposé comme alternative à l’indice global est le rapport signal à dyspériodicité segmental. Il a pour objectif de donner une plus forte pondération aux segments bruités de faibles niveaux qui sont peu pondérés dans le calcul de l’indice global. La méthode d’évaluation perceptive développée est basée sur la comparaison de paires de signaux. Elle permet d’obtenir des résultats fiables même lorsque l’évaluation est réalisée par des auditeurs naïfs n’ayant pas d’expérience dans l’évaluation de la qualité de la voix et permet une grande concordance inter-juges et intra-juges. Les performances de la méthode d’évaluation perceptive basée sur la comparaison de paires de signaux sont comparées à celles de la méthode conventionnelle utilisée en milieu clinique. Les performances des différentes méthodes d’analyse et des indices acoustiques sont mesurées en les testant sur des corpus comprenant des voyelles soutenues et de la parole continue. Les stimuli sont produits par des locuteurs normophoniques et dysphoniques et comprennent une large gamme de pathologies.

Les recherches en dosimétrie numérique non-ionisante ont connu des avancées considérables notamment grâce à la capacité de calcul croissante des ordinateurs. Ces dernières années, la prise en compte de la variabilité dans la dosimétrie est devenue un enjeu majeur. Les sources de variabilité sont nombreuses, parmi elles, les conditions d'exposition au rayonnement électromagnétique peuvent induire des doses absorbées très différentes. La modélisation de canal de propagation en environnement intérieur a permis d'avoir une connaissance précise des conditions d'exposition d'un corps humain plongé dans cet environnement. Cette thèse a pour but de développer une méthode statistique du calcul de la dose absorbée par le corps et d'adapter le modèle de canal stochastique à la dosimétrie. L'étude statistique de l'exposition a révélé la nécessité d'obtenir les valeurs de Débit d'Absorption Spécifique corps entier dans le cas d'expositions à une onde plane pour tous les angles d'incidence possibles. Compte tenu des temps de calcul particulièrement long en dosimétrie, une méthode d'interpolation efficace, le krigeage, des valeurs de Débit d'Absorption Spécifique a été mise en oeuvre. L'analyse de sensibilité aux paramètres du canal des moments du Débit d'Absorption Spécifique a permis de connaître l'impact de chacun de ces paramètres. Le modèle de canal a pu être simplifié et donc adapté à la dosimétrie tout en quantifiant l'erreur d'approximation qu'implique cette simplification. Cette thèse répond à l'enjeu de la prise en compte de la variabilité en dosimétrie dans un environnement complexe. / Research in non-ionizing numerical dosimetry has been improved thanks to high calculation capacity of computers. These years, integrating variability in the field of dosimetry has become a major issue. Sources of variability are numerous; among them, there are the exposure conditions to electromagnetic radiation which can lead to very different absorbed doses. Indoor channel modeling enables to have a deep knowledge of the exposure conditions of a human body located inside this indoor environment. The aim of this thesis is to develop a statistical method of calculation of the absorbed dose by the human body and to adapt the stochastic channel model to dosimetry. The statistical study of exposure reveals the need to obtain Specific Absorption Rate values for a plane wave exposure for all possible angles of incidence. Taking into account that computation in dosimetry is time consuming, an efficient interpolation method, kriging method, is implemented in order to get whole body Specific Absorption Rate values. Kriging method enables to obtain Specific Absorption Rate for all possible angles of incidence and then to calculate expectation and variance of Specific Absorption Rate. Sensitivity Analysis of expectation and variance to the statistical channel parameters reveals the impact of each parameter. The channel model has been simplified and then adapted to dosimetry by estimating the approximation error induced by this reduction. This thesis answers to the issue of integrating variability in dosimetry in a complex environment and develop the tools that open a new path in studying exposure in any complex environment.

L'AMÉLIORATION DES TECHNIQUES D'OBSERVATION DES PRÉCIPITATIONS À L'ÉCHELLE RÉGIONALE EST L'UN DES ENJEUX FORTS DE L'OBSERVATOIRE HYDROMÉTÉOROLOGIQUE MÉDITERRANÉEN CÉVENNES-VIVARAIS. SUITE À UNE PREMIÈRE EXPÉRIMENTATION RÉALISÉE À BOLLÈNE À L'AUTOMNE 2002, LES DONNÉES RECUEILLIES PAR LE RADAR MÉTÉOROLOGIQUE OPÉRATIONNEL DE BOLLÈNE DU RÉSEAU ARAMIS ONT SERVI DE SUPPORT À UNE CHAÎNE D'ALGORITHMES MISE AU POINT POUR L'ESTIMATION QUANTITATIVE DES PRÉCIPITATIONS AU SOL. <br />NOTRE CONTRIBUTION À CES TRAITEMENTS CONCERNE PARTICULIÈREMENT UNE SOURCE D'ERREUR POUR LA MESURE DES PRÉCIPITATIONS PAR RADAR MÉTÉOROLOGIQUE : L'HÉTÉROGÉNÉITÉ VERTICALE DES PRÉCIPITATIONS. NOUS NOUS PLAÇONS DANS LA CONTINUITÉ DE LA MÉTHODE À BASE STATISTIQUE D'IDENTIFICATION DU PROFIL VERTICAL DE RÉFLECTIVITÉ PAR INVERSION PROPOSÉE PAR ANDRIEU ET CREUTIN (1995). NOUS AVONS ÉTENDU CETTE MÉTHODE AU CAS DE SUPPORTS GÉOGRAPHIQUES ÉVOLUTIFS DANS LE TEMPS AFIN DE CONSIDÉRER LE TYPE DE PRÉCIPITATIONS (CONVECTIF, STRATIFORME). <br />L'ÉVALUATION DES TRAITEMENTS, PAR LA CARACTÉRISATION DE L'INCERTITUDE DE L'ESTIMATION DE PLUIE PAR RADAR, EST UN ASPECT IMPORTANT DU TRAVAIL. NOUS AVONS DÉVELOPPÉ UNE TECHNIQUE BASÉE SUR LA GÉOSTATISTIQUE POUR ÉTABLIR UNE PLUIE DE RÉFÉRENCE À PARTIR DE DONNÉES PLUVIOMÉTRIQUES. NOUS AVONS ÉGALEMENT DÉVELOPPÉ UN MODÈLE POUR CARACTÉRISER LA STRUCTURE STATISTIQUE DES ERREURS AINSI QUE LEUR STRUCTURE SPATIO-TEMPORELLE. <br />NOUS AVONS ENFIN DÉVELOPPÉ UNE APPROCHE D'IDENTIFICATION DU PVR À BASE PHYSIQUE. EN DIMINUANT LE NOMBRE DE PARAMÈTRE PAR RAPPORT À LA MÉTHODE INITIALE ET EN INTRODUISANT DES CONTRAINTES PLUS PHYSIQUES, CETTE APPROCHE VISE À AMÉLIORER LA ROBUSTESSE D'IDENTIFICATION DU PVR.

Face à la contradiction consistant à traiter les variables régionalisées écologiques sans tenir compte de leurs propriétés spatiales, nous développons des méthodes géomatiques, utilisant des techniques informatiques, et géostatistiques, appliquant la théorie des fonctions aléatoires. Après avoir introduit des éléments de géomatique et de géostatistique, et avoir précisé la nature spécifique de l'autocorrélation spatiale, nous introduisons les inférences design-based et model-based. A l'aide de fonctions aléatoires, nous étudions l'efficacité de l'échantillonnage, optimisons des prédicteurs, et calculons des intervalles de prédiction. Nous proposons une procédure d'optimisation des classes de distances lors du calcul du variogramme. Nous examinons également l'utilisation de l'intégrale du variogramme, et justifions la modélisation du variogramme par ajustement aux moindres carrés pondérés. Nous discutons de la précision du variogramme dans les cadres design-based et model-based, et au sens du jackknife. L'optimisation de l'échantillonnage d'une population finie en vue de l'estimation de la moyenne spatiale ou du variogramme est examinée à l'aide de plusieurs heuristiques d'optimisation combinatoire et de simulations de fonctions aléatoires. Le problème du test de la corrélation ou de l'association entre deux variables régionalisées est étudié, à nouveau en utilisant des simulations de fonctions aléatoires. Nous passons en revue plusieurs méthodes et recommandons les tests qui font explicitement référence à l'autocorrélation spatiale des variables régionalisées. Dans le cadre de la définition de l'association spatiale entre variables régionalisées, nous proposons une méthode hybride utilisant des quadtrees et une distance d'édition entre arborescences récursives. Enfin, nous étudions des mesures de la complexité spatiale, critiquons l'analyse fractale et proposons des méthodes alternatives, notamment une mesure de complexité topologique d'une carte en isolignes.