Implémentation d'un algorithme de reconstruction d'images en élastographie par résonance magnétique

Tuteurs :
  • Jonathan Vappou
  • Nadège Corbin
Groupe :
  • Myriam ABDELKAFI (Responsable projet)
  • Margaux BEAUVALLET
  • Nicolas CASTILLON
  • Valentin GUEVEL

Contexte :


L'élastographie par résonance magnétique (E.R.M.) est une technique d'imagerie récente, utilisant la technologie de l'IRM.
L'onde émise se propage dans les tissus mous de l'organisme afin de différencier les tissus malades des tissus sains, en particulier dans le cas de tumeurs, cancers…


Ci dessus: ERM d'un foie (avec à gauche l'image d'amplitude obtenu grâce à l'IRM, au milieu l'image de phase où on voit l'onde se propager; et à droite l'élastogramme du foie)

Ci dessous: Elastographie d'un sein (la partie en bleu indique ici une lésion plus dure, et donc la présence possible d'une tumeur)



L'ERM permet au radiologue de “palper” l'image, de connaître l'élasticité du tissu.
Pour le moment, cette technique d'imagerie ne peut être utilisée que dans le cas de diagnostics. Cependant si l'on parvenait à l'améliorer, elle pourrait informer un radiologue interventionnel en temps réel lors d'une opération.


Objectif :


Implémentation en langage C d'un code Matlab (permettant la représentation de l'élastogramme d'un objet) pour un souci de portabilité et de rapidité.


Principe :


L'ERM nécessite de passer par 3 étapes:

  • Excitation mécanique
  • Encodage des déplacements de l'onde
  • Application de l'algorithme LFE (Local Frequency Estimation)



(Cliquez pour agrandir)

C'est ce troisième et dernier point qui constitue le coeur de notre projet. Il se base sur la résolution d'un problème inverse (c'est à dire obtenir un élastogramme à partir de d'une image de phase).
On applique une transformée de Fourier (FFT) spatiale en 2 dimensions à l'image de phase obtenue à partir du signal complexe S de l'IRM.
–> Φ = arctan(Im(S)/Re(S))

Dans le domaine de Fourier, on applique l'algorithme, qui permet d'estimer la fréquence spatiale k, avant d'appliquer une Transformée de Fourier inverse (IFFT).

Détails sur l'algorithme LFE


Avancement du projet :


Nous vous proposons de consulter sur cette page l'avancement de l'implémentation de notre programme en langage C.


Carnet de Bord
Diagramme de Gantt
Présentations


Présentation n°1 - Présentation : erm_-_wiki_-_presentation.pdf
Présentation n°1 - Codage Matlab : erm_-_wiki_-_codage_sous_matlab.pdf

Présentation n°2 - Anglais : erm_-_wiki_-_presentation_anglais.pdf

Présentation finale : presentation_finale.pdf

Rapport final : projet_ingenieur_tic1a_-_rapport_groupe_4_bis.pdf