Équipe IMAGeS : Images, Modélisation, Apprentissage, Géométrie et Statistique

Différences entre les versions de « Analyse statistique et problèmes inverses »

De Équipe IMAGeS : Images, Modélisation, Apprentissage, Géométrie et Statistique
Aller à la navigation Aller à la recherche
m (Admin a déplacé la page ASPI vers Analyse statistique et problèmes inverses)
 
(102 versions intermédiaires par 4 utilisateurs non affichées)
Ligne 1 : Ligne 1 :
= Analyse Statistique et Problèmes Inverses =
+
''Responsable'' : [[Erik-André Sauleau]]
  
Ce thème s'intéresse à l’ensemble des approches méthodologiques basées sur l’utilisation des statistiques et des outils associés à l’inférence bayésienne pour l’analyse, le recalage, la segmentation, la détection, la classification, la fusion et la visualisation des signaux et des images au sens large, incluant les observations multivariées ou multimodales. L’approche bayésienne que nous privilégions s’appuie sur la théorie des modèles graphiques et de la marginalisation, elle a pour but non seulement l’estimation des variables d’intérêt mais également l’évaluation des incertitudes liées aux paramètres estimés.
+
<br />
  
[[File:image_156.jpeg|frameless|center|907px]]
+
''Participants permanents'' : C. Collet, S. Faisan, C. Heinrich, F. Heitz, V. Mazet, N. Meyer, V. Noblet, F. Rousseau, E.-A. Sauleau, L. Thoraval
 +
<br />
 +
''Doctorants (thèses en cours)'' : A. Bouchon, F. Séverac.
  
===== Bilan des activités sur la période 2007-2011 =====
 
  
Le thème a travaillé sur les familles de méthodes markoviennes où les processus de segmentation en jeu sont des opérations complexes lorsque le bruit prédomine sur le signal utile et nécessite la prise en compte de l'incertitude sur les données observées. On peut tenter de contourner cette difficulté en utilisant, quand elle est disponible, une information multimodale ou multispectrale pour renforcer la robustesse du processus de décision, même si cette approche multivariée a ses limites, liées en particulier à la complexité croissante du modèle (phénomène de Hughes) et à la prise en compte adéquate du bruit multivarié (théorie des copules). Une approche complémentaire a été explorée, qui consiste à introduire dans nos modèles markoviens, la notion d'imprécision intrinsèque des étiquettes dans le processus d’étiquetage. Nous avons montré que la modélisation floue des états cachés sur chaînes ou sur champs markoviens pouvait améliorer le processus de segmentation, en prenant en compte la variation locale de la densité de probabilité du bruit (non stationnarité). Autour de la même problématique de segmentation, et dans le cadre de notre collaboration avec le département d’informatique de l’université de Ioannina (Grèce), nous avons contribué à l’élaboration de modèles contextuels gaussiens et d’algorithmes probabilistes innovants permettant d’autoriser des ruptures dans le vecteur décrivant les paramètres du mélange et de proposer des algorithmes d’estimation des paramètres et hyperparamètres impliqués. L’estimation a été réalisée par algorithme Expectation-Maximization et par approche bayésienne variationnelle. Ces modèles ont été mis en œuvre en imagerie physique (traitement de données polarimétriques) et pour des problèmes bidimensionnels en imagerie cérébrale.
+
Ce thème s'intéresse à l’ensemble des approches méthodologiques basées sur l’utilisation des statistiques et des outils associés à l’inférence bayésienne pour l’analyse d'images et la résolution de problèmes inverses. L'objectif est le développement de modèles, de méthodes d'inférence et d'algorithmes originaux, dans un contexte générique, ou dans le cadre d'applications particulières. Les problèmes traités incluent la reconstruction, le recalage, la segmentation, la détection, la classification, la fusion, et la visualisation des signaux et des images au sens large, incluant les observations multivariées ou multimodales.
  
Des travaux ont été menés sur la décomposition conjointe de spectres sur des cubes radio-astronomiques. Afin de distinguer des structures cinématiques (gaz) dans ces cubes, nous avons proposé de décomposer les spectres indépendamment les uns des autres en une somme de raies dont le nombre et les paramètres (position, amplitude et largeur) sont estimés : c'est un problème d'extraction d'objets établi ici dans le cadre bayésien où l'algorithme RJMCMC permet ensuite de calculer une estimation. Nous avons développé une méthode basée sur la décomposition conjointe de spectres. L'approche privilégiée consiste à décomposer tous les spectres en même temps, et non plus séparément. Le modèle est établi dans un cadre bayésien : l'a priori sur les paramètres des raies (positions, amplitudes, largeurs) est un champ de Markov gaussien qui permet à la fois de contraindre les décompositions à être semblables et d'adoucir l'évolution des paramètres au sein de la séquence de spectres.
 
  
Des développements originaux utilisant des modèles graphiques pour la fusion d'images hétérogènes (prises avec différents capteurs, résolutions et bandes spectrales) ont été proposés à travers l'estimation d'un seul modèle, de résolutions spatiale et spectrale arbitraires. Il s'agit de construire un objet géométrique et radiométrique qui préserve l'information utile présente dans les données initiales. L'originalité de l'approche réside dans l'intégration des processus d'acquisition d'images et d'un modèle des objets observés (images 2D ou terrains 3D). Pour cela nous utilisons une approche probabiliste qui permet la modélisation des processus de formation des données images (géométrie, flou, échantillonnage, etc.), et qui s'appuie sur la théorie des modèles graphiques (modèle direct) et l'inversion par marginalisation et optimisation fonctionnelle. Le but visé est l'estimation des variables d'intérêt mais aussi l'évaluation des incertitudes liées aux paramètres estimés. La notion de fusion de données permet ainsi de pallier la redondance d'information et de profiter de la complémentarité des données. Grâce aux incertitudes, la fusion peut être effectuée de manière récursive au moyen de mises à jour successives du modèle.
+
== Modèles et décision statistique ==
  
D’autres travaux ont porté sur le développement de méthodes de modélisation statistique en imagerie tensorielle (avec des applications en détection de changements en IRM de diffusion, dont les données sont intrinsèquement multivariées). Les méthodes développées s’appuient sur les approches élaborées récemment pour définir des modèles et des méthodes d’analyse statistique sur des espaces qui ne sont pas des espaces vectoriels (variétés riemanniennes). Des modèles vectoriels gaussiens, associés à des tests d’hypothèses ont dans un premier temps été évalués. Des modèles log-normaux ont ensuite été considérés pour des tenseurs d’ordre 2 (espace des matrices définies positives). Les recherches en cours portent sur des tenseurs d’ordre supérieur à 2, qui offrent davantage de puissance de modélisation dans le contexte de l’IRMd.
+
=== Modèles markoviens ===
  
En complément des travaux précédents, nous nous intéressons à la représentation et à la modélisation probabiliste de jeux de données de grande dimension. Il s’agit d’abord de capturer les degrés de liberté des données, supposées appartenir à des variétés, c’est-à-dire de représenter le jeu de données de référence dans un espace de faible dimension. Il s’agit ensuite d’élaborer des modèles probabilistes dans cet espace réduit et de mettre en oeuvre des méthodes à noyaux pour réaliser partitionnement et classification, toujours dans l’espace réduit. Ces modèles sont ensuite utilisés pour l’analyse de nouvelles données, à des fins de diagnostic et de suivi thérapeutique.
+
[[File:chaine_Markov.png|frameless|left|400px]]
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
  
Dans le contexte de la modélisation des formes et de leur variabilité, des contributions ont également été apportées à l’introduction decontraintes statistiques de type géométrique dans les contours actifs orientés région et les modèles géodésiques. Des modèles compactsde formes (bases polynomiales de Legendre) associés à des représentations invariantes par transformation affines ont été développés et appliqués à la segmentation d’images.
+
* Modèles markoviens cachés de séquences d'événements. Application à la cartographie IRMf cérébrale. [http://scd-theses.u-strasbg.fr/886/01/Faisan2004.pdf Sylvain Faisan]. [http://dx.doi.org/10.1016/j.media.2006.09.003 Article 1]. [http://dx.doi.org/10.1016/j.acra.2004.09.012 Article 2.]
  
Nous avons enfin travaillé sur l'analyse de sons pour la détection de signatures précoces associées à des pathologies pulmonaires. Nous avons proposé une méthode de détection des phases respiratoires (inspiration, expiration, apnée) et le recalage permettant la localisation des descripteurs dans le cycle basée sur l'utilisation de chaînes de Markov permettant de prendre en compte la structure imprécise mais régulière du cycle respiratoire. L'approche par théorie de l'évidence a également été abordée pour modéliser l'imprécision intrinsèque du modèle. Nous avons ensuite proposé et validé un nouveau descripteur d'anormalités respiratoires basé sur la distribution des coefficients du signal décomposé en paquet d'ondelettes, en modélisant l'attache aux données multivariées corrélés issues de diverses observations inter et intra patient. Cette approche a permis la définition
+
* [http://dx.doi.org/10.1016/j.patrec.2012.03.009 Comparaison entre un sujet et un modèle statistique.] Application à la détection d'anomalies en chirurgie maxillo-faciale. Sylvain Faisan.
 +
 
 +
* Modèles de Markov cachés. [http://lsiit-miv.u-strasbg.fr/lsiit/perso/collet/ftp/PhD/PhD_LeCam.pdf Application à l'analyse temps-fréquence des sons respiratoires pour l'identification de signatures pulmonaires.] Thèse Steven Le Cam.
 +
 
 +
* Modèles de Markov cachés. Application à la segmentation d'IRM 3D et à la détection de lésions. [http://scd-theses.u-strasbg.fr/1651/01/BRICQ_Stephanie_2008-protege.pdf]. Thèse Stéphanie Bricq [http://dx.doi.org/10.1016/j.media.2008.03.001 Article].
 +
 
 +
* Chaînes de Markov 3D+t. Application à la segmentation statistique 3d+t d’images IRM multimodales pour la quantification de la charge lésionnelle. <span style="color:red">Thèse en cours : François Lavigne.</span>
 +
 
 +
=== Modèles graphiques, modèles bayésiens hiérarchiques ===
 +
 
 +
[[File:modele_graphique.png|frameless|left|250px]]
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
 
 +
* Modèles graphiques. Application à la détection de changements du bâti en imagerie de télédétection multibande et hyperspectrale. <span style="color:red">Thèse en cours : Thanh Ngo.</span>
 +
 
 +
=== Modèles non linéaires compacts. Apprentissage de variétés ===
 +
 
 +
[[File:swissroll.png|frameless|left|380px]]
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
 
 +
* Modèles à variables latentes : modèles d'observations non gaussiens et modèles non linéaires dans l'espace réduit : [http://scd-theses.u-strasbg.fr/881/01/VIK2004.pdf Thèse Vik Torbjorn]. [http://dx.doi.org/10.1007/s10044-004-0230-5 Article 1], [http://dx.doi.org/10.1109/TPAMI.2007.1028 Article 2].
 +
 
 +
* Apprentissage statistique de variétés non linéaires : une approche par MDS. [http://scd-theses.u-strasbg.fr/1565/01/BRUCHER_Matthieu_2008.pdf Thèse Matthieu Brucher]. [http://asp.eurasipjournals.com/content/pdf/1687-6180-2008-862015.pdf Article].
 +
 
 +
* [http://scd-theses.u-strasbg.fr/2461/01/RENARD_Felix_2011.pdf Modèles non linéaires statistiques compacts pour des données de grande dimension. Application en neuro-imagerie : étude des formes, des structures, de la connectivité cérébrale. Thèse Félix Renard]
 +
 
 +
* Modèles non linéaires sur des espaces tensoriels. Application à la détection de changements en IRM de diffusion.] [https://publication-theses.unistra.fr/public/theses_doctorat/grigis_antoine_2012_ED269.pdf Thèse Antoine Grigis]. [http://dx.doi.org/10.1016/j.neuroimage.2012.02.049 Article 1]. [http://dx.doi.org/10.1016/j.media.2013.01.004 Article 2].
 +
 
 +
* Modélisation de tenseurs d'ordre supérieur en IRM de diffusion. Application à la comparaison de populations et à l'extraction de biomarqueurs. <span style="color:red">Thèse en cours : Theodosios Gkamas.</span>
 +
 
 +
* Modèles statistiques linéaires / non linéaires pour les études de groupe en imagerie ITD. <span style="color:red">Thèse en cours : Alix Bouchon.</span>
 +
 
 +
* Modélisation non linéaire compacte d'images multimodales TEP / ITD. Application à la comparaison patient-atlas de normalité. <span style="color:red">Thèse en cours : Florent Savary.</span>
 +
 
 +
== Problèmes inverses ==
 +
 
 +
=== Imagerie médicale ===
 +
 
 +
[[File:motion_compensation.png|frameless|left|300px]]
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
 
 +
* [http://dx.doi.org/10.1016/j.media.2010.04.005 Super-résolution en imagerie IRM, intégrant un a priori multimodal. François Rousseau.]
 +
 
 +
* [http://dx.doi.org/10.1109/TMI.2009.2030679 Reconstruction 3D en imagerie foetale in utero à partir de coupes IRM. François Rousseau.]
 +
 
 +
* [http://dx.doi.org/10.1109/ICIP.2009.5413470 A priori de formes en reconstruction tomographique orientée objet. Thèse Gil Gaullier.]
 +
 
 +
=== Imagerie polarimétrique ===
 +
[[File:poincare.png|frameless|left|300px]]
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
 
 +
* Inversion de données en imagerie polarimétrique : imagerie de Stockes, imagerie de Mueller, Non Local Means, inversion sous contrainte d'admissibilité. S. Faisan, C. Heinrich, J. Zallat. [http://dx.doi.org/10.1364/JOSAA.28.000465 Article 1]. [http://dx.doi.org/10.1364/JOSAA.29.002028 Article 2]. [http://dx.doi.org/10.1364/OE.21.004424 Article 3].
 +
 
 +
=== Signaux spectroscopiques ===
 +
[[File:spectrum.png|frameless|left|230px]]
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
* [http://scd-theses.u-strasbg.fr/1548/01/GENDRIN_Christelle_2008.pdf Spectroscopie proche infra rouge et imagerie chimique non destructive. Application au contrôle de qualité dans l'industrie pharmaceutique. Thèse Christelle Gendrin.]
 +
 
 +
* [https://publication-theses.unistra.fr/public/theses_doctorat/Belghith_Akram_2012_ED269.pdf Indexation de spectres HSQC et d'images IRMF appliquée à la détection de bio-marqueurs. Thèse Akram Belghith.]
 +
 
 +
* [http://lsiit-miv.u-strasbg.fr/mazet/spectrodec/index.html Décomposition conjointe de séquences de signaux spectroscopiques] : modèles bayésiens hiérarchiques, champ markoviens gaussiens, RJMCMC, spectroscopie de photoélectrons, radioastronomie. Vincent Mazet, Sylvain Faisan.
 +
 
 +
=== Imagerie hyperspectrale ===
 +
 
 +
[[File:AstroHyperspectral.jpg|frameless|left|200px]]
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
<br />
 +
 
 +
* [http://lsiit-miv.u-strasbg.fr/paseo/cubefusion.php Fusion bayésienne optimale de données hyperspectrales multi-observations  (2d+lambda+t).] Christophe Collet.
 +
 
 +
* [http://lsiit-miv.u-strasbg.fr/paseo/cubedeconvolution.php Déconvolution de données astronomiques hyperspectrales.] Vincent Mazet.
 +
 
 +
[[en:Statistical_analysis_and_inverse_problems]]

Version actuelle datée du 26 août 2016 à 10:44

Responsable : Erik-André Sauleau


Participants permanents : C. Collet, S. Faisan, C. Heinrich, F. Heitz, V. Mazet, N. Meyer, V. Noblet, F. Rousseau, E.-A. Sauleau, L. Thoraval
Doctorants (thèses en cours) : A. Bouchon, F. Séverac.


Ce thème s'intéresse à l’ensemble des approches méthodologiques basées sur l’utilisation des statistiques et des outils associés à l’inférence bayésienne pour l’analyse d'images et la résolution de problèmes inverses. L'objectif est le développement de modèles, de méthodes d'inférence et d'algorithmes originaux, dans un contexte générique, ou dans le cadre d'applications particulières. Les problèmes traités incluent la reconstruction, le recalage, la segmentation, la détection, la classification, la fusion, et la visualisation des signaux et des images au sens large, incluant les observations multivariées ou multimodales.


Modèles et décision statistique

Modèles markoviens

Chaine Markov.png








  • Modèles de Markov cachés. Application à la segmentation d'IRM 3D et à la détection de lésions. [1]. Thèse Stéphanie Bricq Article.
  • Chaînes de Markov 3D+t. Application à la segmentation statistique 3d+t d’images IRM multimodales pour la quantification de la charge lésionnelle. Thèse en cours : François Lavigne.

Modèles graphiques, modèles bayésiens hiérarchiques

Modele graphique.png










  • Modèles graphiques. Application à la détection de changements du bâti en imagerie de télédétection multibande et hyperspectrale. Thèse en cours : Thanh Ngo.

Modèles non linéaires compacts. Apprentissage de variétés

Swissroll.png









  • Modélisation de tenseurs d'ordre supérieur en IRM de diffusion. Application à la comparaison de populations et à l'extraction de biomarqueurs. Thèse en cours : Theodosios Gkamas.
  • Modèles statistiques linéaires / non linéaires pour les études de groupe en imagerie ITD. Thèse en cours : Alix Bouchon.
  • Modélisation non linéaire compacte d'images multimodales TEP / ITD. Application à la comparaison patient-atlas de normalité. Thèse en cours : Florent Savary.

Problèmes inverses

Imagerie médicale

Motion compensation.png









Imagerie polarimétrique

Poincare.png










  • Inversion de données en imagerie polarimétrique : imagerie de Stockes, imagerie de Mueller, Non Local Means, inversion sous contrainte d'admissibilité. S. Faisan, C. Heinrich, J. Zallat. Article 1. Article 2. Article 3.

Signaux spectroscopiques

Spectrum.png










Imagerie hyperspectrale

AstroHyperspectral.jpg