Ce travail de thèse traite de la relation entre les propriétés thermiques effectives de matériaux fibreux à base de bois d'une part et leur microstructure d'autre part. La technique de prise de moyenne est utilisée, dans le cadre général du non-équilibre thermique local, pour écrire les équations de transfert macroscopique de matériaux pouvant présenter une anisotropie locale de conductivité thermique (liée ici aux propriétés intrinsèques des fibres de bois) et exprimer le tenseur de conductivité thermique effective. Le modèle est renseigné par l'analyse quantitative d'images 3D de matériaux réels acquises par microtomographie X. Des outils issus de la morphologie mathématique sont mis en oeuvre pour caractériser finement la microstructure du réseau fibreux et en quantifier l'anisotropie locale et globale. Un volume élémentaire représentatif (VER) qui satisfait aux contraintes de la méthode de prise de moyenne peut être défini grâce à ces informations. Une méthode de segmentation basée sur un algorithme de squelettisation par amincissement homotopique est également développée afin d'identifier chaque fibre individuellement, ce qui permet d'accéder à de nombreux paramètres comme la longueur, la tortuosité ou encore le nombre de contacts. Le modèle thermique macroscopique développé ainsi que son implémentation numérique sont d'abord validés par une comparaison avec les prédictions théoriques dans des cas simples puis les valeurs des tenseurs de conductivité thermique effective, obtenus à partir des images de différents matériaux, sont finalement confrontés aux résultats expérimentaux.

La tomographie par rayons X permet un accès direct aux caractéristiques internes des billons scannés par l'intermédiaire des variations de densité et de teneur en eau. L'objectif de ce travail est de démontrer la faisabilité de la détection automatique de caractéristiques dans des billons à des fins d'analyses scientifiques. Nous disposons pour cela d'une base d'images tomographiques de 24 épicéas communs (Picea abies (L.) Karst.) obtenues avec un scanner médical. Les arbres sont représentatifs de différents statuts sociaux et proviennent de quatre peuplements du nord-est de la France, eux-mêmes sont représentatifs de plusieurs classes d'âge, de densité et de fertilité. Les procédures de détection développées sont les suivantes : - détection de la moelle dans les billons, y compris en présence de nœuds et/ou d'excentricité des cernes. La précision de la localisation est inférieure au millimètre ; - détection de la limite aubier/duramen dans les billons, y compris en présence des nœuds (principale source de difficulté). L'erreur sur le diamètre du duramen est de 1.8mm soit une erreur relative de 1.3%. ; - détection de la localisation des verticilles et une comparaison à une méthode optique ; - détection des nœuds individualisés. Cette procédure permet de compter les nœuds dans un billon et de les localiser (hauteur dans le billon et azimut); cependant, la validation de la méthode et l'extraction du diamètre et de l'inclinaison des branches reste à effectuer. Une application de ce travail a été l'analyse de la variabilité de la quantité d'aubier dans le tronc: en intra-arbre, la largeur d'aubier était constante sous la base du houppier; en inter-arbre, une forte corrélation avec la quantité de branches vivantes a notamment été mise en évidence. De nombreuses analyses sont envisageables à partir des résultats de notre travail, parmi lesquelles: l'étude architecturale des arbres avec le suivi de la moelle dans les billons et l'occurrence des prises de relais sur l'axe principal, l'analyse des variations radiales de la forme du duramen, l'analyse de la distribution des nœuds dans les billons.

La partie centrale de la dépression Afar (Sud Mer Rouge) est bordée de blocs sans déformation récente (Danakil et d'Aysha). En raison du climat désertique, la morphologie des structures récentes et/ou actives est très bien conservée. Dans le cadre du programme TECTOSCOPE Afar financé par l'INSU/CNRS et pour la partie sud de la dépression, l'analyse détaillée de 8 couples d'images stéréoscopiques SPOT a permis non seulement de cartographier et d'interpréter le champ de failles mais aussi de mesurer le rejet vertical des failles. Cette interprétation a été réalisée d'une part sur des tirages photographiques Vizir réalisés à l'INSU, d'autre part sur station de travail par images "anaglyphes". Il a ainsi été possible de définir la géométrie et la cinématique du champ de déformation finie post-stratoïde (4Ma). Les caractéristiques du champ de failles ont été déterminées en utilisant des critères structuraux (intersection de failles, segmentation de coulées volcaniques, etc ... ) et l'analyse statistique de la distribution spatiale des orientations, rejets, et densités de la fracturation. Deux domaines de déformations (blocs basculés et horst-grabens) sont caractérisés et des zones de transfert entre ces domaines et les zones non déformées accommodent les variations latérales de la déformation. Cette étude met en évidence une direction/ moyenne d'extension N020° dans le centre de la dépression, s'orientant N040° en bordure du bloc danakil. Cette configuration structurale résulte d'une déformation progressive du domaine central de l'Afar et les hétérogénéités du champ de déformations sont interprétables en termes d'effets de bordure le long des blocs non déformés Dankaly et Aysha. Les implications cinématiques en sont déduites et comparées avec des expériences sur modèles réduits fragilesductiles inspirés du modèle de la "Biellette" Danakil (Sichler 1980). La rotation senestre de 10° du bloc Danakil offre une explication cohérente de la geométrie du champ de failles, des rotations dextres des blocs basculés du Sud de la dépression et de l'extension oblique du golfe de Tadjoura.

Ce document d'habilitation est consacré à l'étude de séries temporelles d'images de télédétection à basse (LSR) et à très haute résolution spatiale (VHSR). Les phénomènes étudiés concernent la météorologie et l'océanographie (données LSR) et l'agriculture et le milieu urbain (données VHSR). La fréquence d'acquisition des données satellites est inversement proportionnelle à la résolution spatiale. Ainsi, pour des données LSR, la cadence d'acquisition est élevée (une images pour $15min$ avec le satellite MSG --Météosat Seconde Génération) et cela autorise l'étude des mouvement atmosphériques turbulents observés à travers le mouvement des nuages, la circulation océanique, ... Les problèmes d'analyse d'images associés concernent l'estimation du mouvement, le suivi de courbe ou encore l'interpolation de données manquantes. En ce qui concerne les données VHSR, la période séparant deux images peut varier de quelques semaines à quelques mois. Les études associées sont alors dédiées à la détection de changements structurels entre deux images. Du point de vue méthodologique, l'analyse de données LSR est principalement réalisée en introduisant des connaissances physiques a priori dans les outils classiques d'analyse d'images. Une part importante est dédiée à l'utilisation de techniques d'assimilation variationnelle de données. Pour l'analyse de données VHSR, nous proposons des descripteurs spécifiques permettant de caractériser les motifs texturés que l'on a à manipuler. Ces descripteurs sont ensuite utilisés pour résoudre différents problèmes d'analyse de données VHSR tels que la segmentation, la classification, la détection de fronts texturés ou encore l'estimation de l'orientation. Enfin, un chapitre est consacré à la détection de changements où nous proposons des techniques pour le détection binaire et multi-labels.

Cette thèse a pour cadre l'étude de l'organisation spatiale des fibres du muscle cardiaque à partir de séries d'images tridimensionnelles acquises par IRM du Tenseur de Diffusion (IRMTD). Cette organisation constitue une propriété fondamentale du coeur sous-tendant la fonction contractile. Néanmoins elle est très complexe à obtenir au vu des difficultés inhérentes au mouvements cardiaque et respiratoire. Notre objectif consiste à développer de nouvelles approches, basées sur la prise en compte du mouvement du coeur et de la sensibilité au bruit de l'acquisition, pour l'estimation, l'analyse et la visualisation des fibres du myocarde. Dans ce cadre, mes travaux se déclinent selon trois axes principaux. Le premier compare, dans le contexte d'études cliniques ex vivo, les principales approches de régularisation opérant soit sur les images pondérées en diffusion soit sur les champs de tenseurs de diffusion. Les différences sont suffisamment faibles pour conclure que la qualité de nos données IRMTD est suffisante pour considérer toutes les méthodes de régularisation comme équivalentes. Partant de ce constat, une méthode de régularisation simple et rapide apparaî satisfaisante. Le second concerne la mise en place d'une méthode de tractographie spécialement conçue pour la spécificité cardiaque. Celle-ci est guidée par une fonctionnelle de coût globale qui permet l'estimation automatique des fibres cardiaques en une seule fois pour l'ensemble des données, et ce sans l'utilisation de points d'initialisation. Le dernier axe consiste en la distinction d'une population de fibres cardiaques en sous-groupes. Celle-ci s'appuie sur la comparaison de méthodes de classification de type géométrique et de type topologique exploitant toutes trois modes différents de représentation des fibres. Les résultats établissent que la classification pourrait permettre l'identification automatique de régions spécialisées dans le myocarde, ce qui pourrait grandement faciliter l'analyse et la comparaison des données IRMTD cardiaques pour la conception de thérapies patient-spécifiques.

Cette thèse se concentre sur l'étude de l'extraction et de la caractérisation des structures locales, dans le contexte de l'analyse d'images et des textures. S'appuyant sur les lignes de niveau des images ou sur la notion duale et moins structuré e d'orientation du gradient, les contributions de cette thèse se concentrent sur trois thèmes suivants: La première partie présente une nouvelle méthode pour l'analyse de texture qui dans l'esprit est similaire a la granulométrie morphologique, tout en permettant un haut degré d'invariance géométrique et radiométrique. Avec l'aide de la représentation par carte topographique, la deuxième partie de cette thèse développe une approche générale pour l'abstraction d'images, dont le but est de générer automatiquement des images abstraites à partir de photographies réalistes. Le sujet de la dernière partie de cette thèse est la détection des jonctions dans les images naturelles. L'approche s'appuie sur les directions locales de lignes de niveau à travers l'orientation du gradient de l'image. Nous introduisons un système générique d'analyse de jonction. Le premier avantage de la procédure proposée est un critère pour la détection automatique de jonctions. Celui-ci permet de traiter des parties texturées de l'image dans lesquelles aucune détection n'est attendue. Deuxièmement, la méthode donne une caractérisation des jonction en L-, Y et en X, y compris un calcul précis de leur type, de leur localisation et de leur échelle. Contrairement aux approches classiques, la caractérisation de l'échelle ne repose pas sur un espace-échelle linéaire, et permet donc d'obtenir une bonne précision géométrique.

Le périmysium a été analysé à l'aide de deux modalités d'imagerie, l'histologie et l'imagerie par résonance magnétique (IRM). L'histologie a permis l'observation des faisceaux primaires et l'IRM des faisceaux d'ordre supérieur. Dans les deux cas, une méthode de seuillage du périmysium a été mise au point et validée avec la construction d'une banque de cartes de référence probabilistes, résultat de la segmentation manuelle du réseau par un jury d'experts. Le périmysium segmenté a été caractérisé à l'aide de variable morpho-anatomiques et granulométriques. Celles-ci ont été ensuite reliées à la résistance mécanique et à la tendreté sensorielle de la viande cuite. Ce sont les structures macroscopiques, visualisées en IRM, qui jouent un rôle majeur sur le déterminisme de la tendreté de la viande bovine cuite. Ces résultats ouvrent des perspectives pour le développement de méthodes d'imagerie simples adaptées à la prédiction et au contrôle de la tendreté de la viande bovine.

L'objectif du travail présenté dans ce mémoire est de développer une méthode pour mettre en évidence des structures non observables directement en surface à partir de données structurales utilisées comme données de champ. Le paramètre choisi est la déformation finie. La fiabilité de la méthode est testée en comparant les restaurations de la déformation (d'après les données de déformation finie) avec les modèles cinématiques issus d'une approche classique de terrain. La région d'application est la Zone Houillère Briançonnaise. Elle est située juste au dessus du Front Briançonnais (FB) qui est une des structures majeures des Alpes occidentales associée au Chevauchement Pennique Frontal (limite entre les zones internes et externes). Le modèle cinématique issu de l'approche classique de terrain a permis de mettre en évidence trois événements principaux suivis d'un épisode extensif. L'événement D1 traduit très probablement la subduction de l'océan piémontais au cours de laquelle, restant à un niveau structural élevé dans le prisme d'accrétion tectonique, la Zone Houillère subit un écaillage vers l'Ouest alors qu'une partie de la Vanoise, située plus à l'Est, est entraînée dans la subduction. Au cours de D2, la subduction de la Vanoise se bloque et cette zone est alors charriée sur la Zone Houillère puis elle l'emboutit. L'événement D3 est associé à l'indentation de la croûte européenne par le poinçon adriatique. Cette indentation est, entre autre, la cause de la surrection des Massifs cristallins externes qui réactive une surface correspondant au FB originel (décollement précoce D1 et/ou structure D2) en faille normale. Le rebroussement résultant de la surrection des Massifs cristallins externes produit une structure de type roll-over. L'événement extensif D4, continuum plus superficiel de D3, provoque alors le basculement de tout l'édifice vers l'Ouest. Les données de déformation finie ont été acquises par l'utilisation d'algorithmes d'analyse d'images, semi-automatisées au cours de ce travail. Les méthodes utilisées sont fondées soit sur l'analyse de forme des grains (diamètres de Féret) soit sur l'analyse de leur position (méthode de Fry normée). La comparaison et l'analyse des résultats obtenus montrent que la déformation enregistrée au cours de l'événement 02 (surface régionale de référence ) est de très faible intensité. Le processus de déformation est l'aplatissement accompagné d'une variation de volume inférieure à 10%. Les données ont été ensuite interpolées et visualisées à l'aide de GOCAO (modeleur 3D développé par l'équipe de J.L. Mallet, ENSG-LIAD). L'inversion de la déformation en plan a été effectuée à l'aide du logiciel FaultPack en collaboration avec l'équipe de P.R. Cobbold (Univ. Rennes 1). Les restaurations, bien qu'incomplètes, confirment, au cours de l'événement D2, l'absence d'une déformation cisaillante, la faible intensité de la déformation et la prédominance de l'aplatissement à l'échelle de la Zone Houillère entre l'Isère et l'Arc. Cependant, ce genre d'approche devra être à l'avenir associé très étroitement à une restauration géométrique 3D préalable, le problème majeur, non encore résolu. étant de faire la part de chaque événement de déformation dans les fabriques mesurées.

L'étude des empilements de grains est importante pour la compréhension des propriétés des milieux granulaires. En particulier, la connaissance des aires de contact entre grains est utile pour décrire la cinétique d'une réaction solide-solide. Nous avons mis au point une méthode d'étude originale des caractéristiques d'un empilement de grains, fondée sur l'analyse d'images de sections de l'empilement enrobe dans une résine. A l'aide de cette méthode, nous avons développé deux protocoles d'analyse spécifiques. Le premier permet de définir complètement l'état d'homogénéité d'un mélange, en évaluant à la fois l'intensité et l'échelle de ségrégation. Le second permet d'estimer les nombres de contacts et les aires élémentaires de contact entre grains dans un mélange. La méthode a été appliquée à l'étude de mélanges de grains de WO₃ et de NiO. L'analyse de l'homogénéité a montré que l'agitation à l'aide d'un mélangeur " Turbula " conduisait a des mélanges peu homogènes. Nous avons pu atteindre un degré d'homogénéité plus avance d'une part en appliquant un mouvement d'agitation moins régulier, d'autre part en utilisant un ajout de liquide. Cependant, le mélangeage par agitation présente toujours une tendance à la démixtion, caractérisée par l'observation d'une ségrégation à large échelle après quelques heures d'agitation. L'analyse des contacts nous a permis de comparer les valeurs expérimentales des nombres de contacts dans nos mélanges avec des valeurs théoriques données pour des empilements de sphères. Certains écarts au modèle ont pu être expliques par les caractéristiques particulières du système granulaire utilise. Enfin, l'analyse de l'influence du mode de construction sur les caractéristiques de l'empilement a été effectuée, et comparée aux résultats d'une simulation numérique. Nous avons ainsi pu montrer les effets d'un mode de construction collectif sur les caractéristiques de l'empilement.

Cette thèse présente les travaux réalisés en vue de l'élaboration d'un modèle cardiaque associant croissance, mouvement et circulation sanguine pour permettre ensuite la construction d'un modèle patient à partir d'un modèle de population. Le premier axe de ce travail est la simulation de la croissance bi-ventriculaire. Un modèle existant de surface unique, calculé à l'aide de méthodes statistiques, a été généralisé à un modèle bi-ventriculaire puis appliqué à la tétralogie de Fallot (ToF). Le deuxième axe concerne la modélisation du mouvement cardiaque au niveau de la population. Un modèle d'ordre réduit basé sur un modèle Polyaffine et LogDemons a été proposé. Il simule la dynamique cardiaque avec peu de paramètres. Les paramètres de transformation sont analysés par des méthodes statistiques. Un modèle de mouvement moyen a été calculé pour représenter le mouvement standard de la population. Le troisième axe s'intéresse à la simulation de l'écoulement sanguin à l'échelle de la population. La complexité des simulations spécifiques à un patient a été réduite grâce à l'utilisation de méthodes d'analyse d'image, de dynamique des fluides numérique et de réduction d'ordre de modèle. La simulation du flux sanguin dans l'artère pulmonaire pour des patients ToF a permis de mieux comprendre l'impact du sang régurgité sur la pression et la vitesse. Étant donné nos contributions sur ces trois axes, nous sommes maintenant en bonne position pour élaborer le modèle couplé des contributions interdépendantes de la croissance, du mouvement et de l'écoulement sanguin. Ce modèle pourrait être utilisé afin d'aider la planification de la thérapie chez les patients atteints de maladies cardiaques.