Cette thèse a pour objet l'étude de l'incertitude attachée à l'estimation de la géometrie d'une scène à partir d'un couple stéréoscopique d'images. La mise en correspondance des points homologues d'un couple suppose la similarité locale des deux images et nécessite une information radiométrique discriminante. Dans de nombreuses situations cependant (déformations géométriques, bruit d'acquisition, manque de contraste, ....), ces hypothèses sont mises en défaut et les erreurs d'appariemment qui en résultent dépendent fortement de l'information contenue dans le couple et non du sytème stéréoscopique lui-meme. <br />Afin d'aborder ce problème, nous proposons un cadre bayésien et l'application de méthodes de Monte Carlo par chaînes de Markov. Celles-ci consistent à simuler la distribution conditionnelle du champ de disparité connaissant le couple stéréoscopique et permettent de déterminer les zones où des erreurs importantes peuvent apparaitre avec une probabilité éventuellement faible. Différents modèles stochastiques sont comparés et testés a partir de scènes stéréoscopiques SPOT, et nous donnons quelques pistes pour étendre ces modèles à d'autres types d'images. Nous nous intéressons également au probleme de l'estimation des paramètres de ces modèles et proposons un certain nombre d'algorithmes permettant une estimation automatique. Enfin, une part importante du travail est consacrée à l'étude d'algorithmes de simulation reposant sur la théorie des chaînes de Markov. L'apport essentiel réside dans l'extension de l'algorithme de Metropolis-Hastings dans une perspective multi-dimensionnelle. Une application performante reposant sur l'utilisation de la loi gaussienne est donnée. De plus, nous montrons comment le recours à des techniques d'échantillonnage d'importance permet de diminuer efficacement le temps de calcul.

Les approches markoviennes en imagerie et vision par ordinateur offrent un cadre mathématique élégant pour résoudre certains problèmes complexes. Le plus souvent, la fonction d'énergie globale modélisant le problème demeure difficile à minimiser. La première contribution de cette thèse consiste alors à proposer des algorithmes d'optimisation efficace de la classe d'énergies markoviennes multi-étiquettes de premier ordre ayant une attache aux données quelconque et un a priori convexe. Les algorithmes proposés reposent sur la technique de coupe-minimum sur graphe et des stratégies d'optimisation itérative par des nouveaux mouvements de partitions larges et multi-étiquettes, qui permettent d'avoir un compromis entre la qualité de l'optimum atteint et la complexité algorithmique. Le cadre applicatif principal de cette thèse est la reconstruction du relief par interférométrie radar à synthèse d'ouverture. Cette méthode de calcul de modèles numériques de terrain est confrontée le plus souvent à la nature très bruitée des données radar et aussi à la complexité des scènes naturelles et urbaines à reconstruire et à leur grande dimension. Ainsi, la seconde contribution de ces travaux de thèse consiste à proposer des modèles markoviens robustes face à la diversité des scènes à reconstruire dans le cas général et des algorithmes d'optimisation qui leur sont appropriés. L'approche générale est testée et validée sur un jeu de données radar synthétiques et réelles (ERS et ESAR).

La description de la géométrie et de la distribution des structures fragiles présentes dans les rifts contribuent à la compréhension des mécanismes de la déformation extensive à l'échelle de la lithosphère. Une analyse quantitative est envisagée sur une portion du rift Est-africain (Gregory Rift au Kenya) et sur la marge Est de la Mer Rouge (Yémen), à partir de données structurales, d'images satellitaires et de données topographiques. Les études de terrain permettent de préciser la géométrie et la chronologie des déformations pour chaque objet étudié (Kenya ou Yémen). Le traitement puis l'interprétation des images numériques SPOT ou Landsat TM permet (i) de compléter les études structurales de terrain et de les intégrer à petite échelle, (ii) de quantifier la distribution en deux dimensions (2D) des failles sur le Gregory Rift. Ces études en 2D s'avèrent insuffisantes pour la quantification de l'extension. L'utilisation de modèles numériques de terrain obtenus par auto-corrélation de couples d'images stéréoscopiques SPOT conduit à la reconstitution de l'état initial, ante-déformation, des systèmes de blocs basculés qui se développent au Yémen. La restauration de l'état initial des blocs basculés comporte une étape "débasculement" ou de dépliage selon le modèle réalisé, et une étape d'ajustement manuel des blocs entre eux. L'ajustement est réalisé en tenant compte soit de la forme des blocs, soit des trajectoires de déformation données par le dépliage des structures. La comparaison de l'état initial restauré et la connaissance de l'état final déformé permet une approche cinématique de l'extension, en terme de déplacements finis. La quantité d'extension est déterminée. Les champs de déplacements 2D (en carte) sont tracés pour chaque secteur étudié au Yémen. Les déplacements 3D sont calculés et projetés sur les plans de failles séparant les blocs basculés, afin de calculer le tenseur de contraintes associé. L'interprétation des résultats obtenus, combinée à l'étude structurale réalisée sur le terrain ou sur les images, fournit (i) les limites de validité de la méthode de quantification employée, (ii) un modèle de découpage en blocs crustaux du sud de la plaque arabique, (iii) la géométrie des phases précoces lors de l'ouverture du Golfe d'Aden.

Le fonctionnement hydrologique des petits bassins versants ruraux de quelques km² à régime intermittent est complexe car de nombreux processus affectent le cheminement des eaux de surface. Il en résulte une variabilité de la densité de drainage et de la dynamique de l'écoulement au sein du réseau hydrographique. Cette dynamique de fonctionnement est peu prise en considération du fait des difficultés d'ordres techniques et économiques pour suffisamment observer les mouvements de l'eau à la surface. Or, c'est une information essentielle pour décrire la connectivité hydrologique du bassin qui représente la distribution spatiale et temporelle des surfaces contributives à l'écoulement par leur connexion au réseau de drainage. De telles connaissances seraient utiles pour traiter la question de la séparation entre production et transfert effectuée dans les modèles hydrologiques, avec la perspective de proposer des simulations plus justes physiquement. L'objectif de la thèse est de proposer des approches spatiales pour mieux intégrer la dynamique de fonctionnement du réseau de drainage en lien étroit avec la réponse hydrologique du bassin. Le Mercier (7 km²) est le site expérimental situé en tête du bassin versant de l'Yzeron (142 km²) localisé à l'ouest de l'agglomération lyonnaise. Ce bassin sur socle cristallin est composé essentiellement de surfaces agricoles et de forêts. Son fonctionnement hydrologique est affecté par l'existence de zones humides contributives. Un réseau de routes et de fossés anthropiques s'ajoute aux talwegs naturels. Les méthodes développées relèvent de deux approches : (1) la microtopographie issue d'un MNT LiDAR (Light Detection and Ranging) permet d'identifier et décrire des extensions fonctionnelles du réseau de drainage d'une part au niveau des linéaires artificiels avec un apport minimal de données externes, et d'autre part au niveau des talwegs naturels en distinguant la présence ou l'absence d'un chenal de drainage, signe d'un potentiel d'écoulement concentré. (2) Un dispositif de 18 capteurs limnimétriques est mis en œuvre pour assurer un suivi permanent de la réponse hydrologique par emboitement de stations au sein du réseau hydrographique. Cette observation « multi-locale » permet de mesurer l'évolution de la densité de drainage, d'identifier localement la hiérarchie des facteurs qui affectent la réponse et de distinguer différentes dynamiques de transfert dans le réseau de drainage. Les résultats des approches par la microtopographie et par observations multi-locales aident à identifier des régions du bassin au comportement différent. Ils permettent notamment de mieux comprendre les interrelations entre occupation du sol et processus hydrologiques, voire géomorphologie et processus. Ces résultats valident donc l'intérêt du capteur LiDAR et celui d'un dispositif in situ souple et adaptable pour proposer un « pattern de drainage » réaliste en limitant le recours au terrain. Enfin, ce pattern décrivant la tendance d'organisation spatiale des écoulements, est paramétré dans une fonction de transfert géomorphologique calculée sur la base des cheminements fournis par un MNT. L'usage de cet outil constitue une ébauche mais conforte l'idée d'un couplage entre le pattern de drainage et la réponse du bassin versant dans des conditions d'intensité pluvieuse soutenue et d'humidité modérée pour expliquer la réponse rapide du bassin. L'ensemble des résultats justifie donc la mise en avant de la nature transitoire du réseau de drainage pour paramétrer des modèles spatialisés avec la perspective d'améliorer leurs capacités prédictives. / Hydrological functioning of small temporary catchments depends on several processes governing transfer from surface water paths. As a result flow dynamics and drainage density are highly variable in space and time. But this complex dynamic is not enough taken into account because of technical and economical limitations. However, it is essential to describe hydrological connectivity as a spatial and temporal pattern of contributive areas to the drainage network. Get this pattern may facilitate the conceptual distinction between production and transfer functions to improve spatially distributed models. The aim of the study is to develop methods to describe spatial and temporal patterns of the drainage network in relation to catchment hydrological responses. The Mercier headwater catchment (7 km²) is located into the Yzeron catchment near Lyon (France). The land use is principally composed of agricultural plots and forested areas. The hydrographical network consists of natural thalwegs and many roadside ditches and agrarian ditches. Both approaches are developed for this purpose: first, the micro-topography from a LiDAR DEM helps to describe potential flow lengths from drainage network extensions during rainfall events. On the one hand, main artificial ditches are mapped from the DEM with minimal corrections from ancillary data. On the other hand, channelized or unchannelized reaches are located from the DEM into the natural thalwegs. Second, a water level sensor device is set up to record hydrological response from 18 stations located in nested sub-catchments into the hydrological network. These synoptic measurements are used to estimate temporal changes in drainage density, to analyze local hydrological functioning, or to describe flood propagation to the outlet. Results from both approaches lead to the identification of specific behaviors inside the hydrological network controlled by functional thresholds. These patterns help to better understand the relationship between land use and hydrological processes. The results also show the interest of LiDAR DEM and the suitability of adaptable distributed measurements as a substitute to heavy in situ studies for the identification of drainage patterns. Finally, to test the assumption of a gradual drainage network extension during a rainfall event, a simplified drainage pattern is processed into a geomorphological transfer function. This tool is fairly easy to set and is used to initiate an association between local hydrological knowledge and global catchment response. Both approaches lead to conclude that various drainage pattern must be integrated into spatially distributed models according to hydrological conditions, rather than a single hydrological network.

L'étude conjointe des parties immergée et émergée de l'Ile de La Réunion (Océan Indien), basée sur la valorisation des données marines multifaisceaux acquises ces 20 dernières années et confrontée à l'ensemble des connaissances géologiques et géophysiques, démontre l'importance considérable des processus de déstabilisations de flanc dans l'évolution des édifices de ce système volcanique océanique. Au total, une cinquantaine d'événements multi-échelles, primaires ou secondaires, d'origine aérienne, côtière ou sous-marine, ont démantelé de façon récurrente les volcans du Piton des Neiges, des Alizés et du Piton de La Fournaise au cours des 2 derniers millions d'années au moins. La plupart correspondent à des phénomènes instantanés et catastrophiques de type avalanche de débris. Leurs dépôts sont principalement accumulés en mer, au sein de quatre gigantesques excroissances topographiques, les Reliefs Sous-Marins, localisés à l'est, au nord, à l'ouest et au sud de l'île. Ces éperons triangulaires, de 20-30 km de large au niveau de la côte et 100-150 km à leur extrémité, à 70-80 km du littoral, se sont construits progressivement par la superposition et la juxtaposition des produits déplacés dans des épisodes de déstabilisation pouvant impliquer plusieurs centaines de km3 de matériaux. Une très faible proportion est retrouvée à terre, au cœur et en bordure littorale du massif du Piton des Neiges et constitue le faciès proximal des dépôts immergés. Les principales structures de la morphologie sub-aérienne de l'île sont interprétées comme (1) les zones sources de déstabilisations relativement récentes, souvent remodelées par l'érosion et/ou comblées par l'activité volcanique ultérieure, ou (2) les traces rémanentes d'amphithéâtres d'avalanches de débris anciens dont les failles bordières ont pu être réactivées dans des processus volcano-tectoniques (effondrements caldériques, étalement gravitaire, subsidence du complexe hypovolcanique du Piton des Neiges,...) et/ou lors d'épisodes ultérieurs de déstabilisation. Ces failles ont probablement guidé le démantèlement par érosion des massifs et ainsi facilité le creusement des dépressions centrales du Piton des Neiges (les trois cirques de Mafate, Salazie, Cilaos et le paléo-cirque des Marsouins) et des principales vallées de l'île. L'étude du volcan Mount Haddington (Ile de James Ross, Antarctique) et les expériences analogiques montrent que la déformation ductile d'horizons de faible résistance inclus dans la structure des volcans boucliers, et en particulier de deltas volcanoclastiques essentiellement constitués de dépôts d'avalanches de débris, est un élément clé dans la génération de glissements rotationnels et d'avalanches de débris ultérieures. Ce modèle est privilégié à La Réunion, en considérant un rôle prépondérant des processus gravitaires dans le déclenchement des déstabilisations des volcans de l'île.