Thème de recherche : Modélisation
mathématique et simulation numérique
en IRM
Le cadre général
de recherche concerne la modélisation mathématique et la simulation
numérique de différents phénomènes intervenant en Imagerie par
Résonance Magnétique (IRM). Ce travail a été mené en
collaboration d’une part avec G. Caloz (IRMAR, UMR 6625, Université de
Rennes I), J. de Certaines et G. Cathelineau (Laboratoire de Résonance
Magnétique en Biologie et Médecine, Université de Rennes I), d’autre
part avec H. Benoit-Cattin et C. Odet (CREATIS, UMR 5515, INSA de Lyon)
et enfin avec L. Chupin (Institut
Camille Jordan, UMR 5208, INSA de Lyon).
Le thème principal de recherche est l’étude des
perturbations du processus d’acquisition d’une image IRM et ces
diverses conséquences. L’obtention d’une image IRM résulte de
l’utilisation de champs magnétiques (statique et radio-fréquence) dont
toute perturbation, qui peut être liée à l’appareillage ou au patient
lui même, se traduit par une déformation de l’image (on parle
d’artefact). L’étude mathématique des artefacts de susceptibilité
magnétique liés à la présence d’un champ magnétique statique
perturbateur induit par un implant médical porté par le patient
(problème majeur en imagerie IRM) qui a commencé dans le cadre de ma
thèse est maintenant achevée. Deux aspects complémentaires liés à ce
problème ont été traités :
- proposer une méthode de calcul de la perturbation du champ
magnétique statique créée par la présence d’un implant médical
métallique dans le cadre d’une expérience IRM;
- modéliser la façon dont les inhomogénéités de champ magnétique
statique engendrent les distorsions de l’image IRM.
La correction des images présentant des artefacts de
susceptibilité magnétique a été traité par B. Belaroussi dans sa thèse
[a] soutenue en octobre 2005 au laboratoire
CREATIS. Les différents aspects de cette partie de notre travail ont
fait l’objet des plusieurs publications [3], [4], [5], [6], [7]
et [8].
Un logiciel (SIMRI) permettant de simuler de manière
complète le fonctionnement d’un appareil IRM est développé au
laboratoire CREATIS. Je collabore à ce projet pour les aspects liés à
la modélisation mathématique des phénomènes physiques intervenant en
IRM, étape préliminaire à toute simulation. À terme, ce simulateur
devra permettre la génération d’images suffisamment réalistes afin
qu’elles soient utilisables pour calibrer et optimiser les techniques
de traitement automatique d’images IRM et ceci aussi bien en IRM bas
champ (utilisée en agro-alimentaire), champ moyen (utilisée en
médecine) et haut champ (utilisé en biologie). Ces traitements
automatique sont aujourd’hui indispensables pour permettre à l’IRM
d’évoluer d’un outil de visualisation à un outil de mesure. Cette
collaboration a donné lieu à la publication [2].
Un autre partie de ce travail concerne les
perturbations liées au champ radio-fréquence. D’une part, les
perturbations du champ RF sont à l’origine d’un type particulier de
distorsions de l’image IRM appelées artefacts radio-fréquence. La
modélisation et la simulation numérique de ce type d’artefact est plus
délicate que celle des artefacts de susceptibilité magnétique car cela
nécessite de reconsidérer l’équation de Bloch qui décrit l’évolution du
moment magnétique macroscopique au cours du temps. Une méthode
appropriée de résolution numérique de l’équation de Bloch a été
proposée et étudiée ; elle fait l’objet de la publication [1].
Outre la perturbation de l’image en cas de champ RF
parasite, un autre effet lié au champ radiofréquence est une possible
élévation de la température des tissus biologiques au voisinage
d’objets métalliques longs comme les cathéters ou les câble de
monitoring utilisés en IRM interventionnelle (micro-chirurgie sous
IRM). La modélisation et la simulation de ce phénomène est relativement
complexe. Elle nécessite de prendre en compte le couplage entre les
phénomènes magnétiques (du type courants de Foucault) et thermiques.
Une première approche de ce problème a été réalisée par H. Bouk’hil
dans sa thèse
[b] soutenue en 2003 au
Laboratoire de Résonance Magnétique en Biologie et Médecine de
l’Université de Rennes 1 à laquelle nous avons collaboré.
[a] Boubakeur Belaroussi, Correction par
traitement
d’images de l’artefact de susceptibilité magnétique dans les images
IRM, thèse de l’INSA de Lyon, 2005.
[b] Hind Bouk’hil, Contribution à la caractérisation
des effets thermiques liés aux biomatériaux métalliques en imagerie par
résonance magnétique, Thèse de l’Université de Rennes 1, Faculté de
Médecine, 2003.
Publications associées à ce thème de recherche
1- S. Balac and L. Chupin. Fast approximate solution of Bloch equation
for simulation of RF artifacts in Magnetic Resonance Imaging.
Mathematical
and
Computer
Modelling (48 : 1901-1913 (2008))
2- S. Balac, H. Benoit-Cattin, T. Lamotte and C. Odet. Analytic
solution to boundary integral computation of susceptibility induced
magnetic field inhomogeneities. Mathematical and Computer Modelling,
39(4-5):437–455, 2004.
3- S. Balac and G. Caloz. Induced magnetic field computations using a
boundary integral formulation. Applied Numerical Mathematics,
41(3):345–367, 2002.
4- S. Balac, G. Caloz, G. Cathelineau, B. Chauvel and J.D. De
Certaines. An integral representation method for numerical simulation
of MRI artifacts induced by metallic implants. Journal of Magnetic
Resonance in Medicine, 45(4):724–727, 2001.
5- S. Balac and G. Caloz. Mathematical modeling and numerical
simulation of magnetic susceptibility artifacts in Magnetic Resonance
Imaging. Computer Methods in Biomechanics and Biomedical Engineering,
3:335–349, 2000.
6- S. Balac. Simulation numérique des artefacts de susceptibilité
magnétique en IRM. Innovation et Technologie en Biologie et Médecine
(ITBM), 19(5):369–379, 1998.
7- S. Balac and G. Caloz. Magnetic susceptibility artifacts in
Magnetic Resonance Imaging : calculation of the magnetic field
disturbances. IEEE Trans. on Magnetics, 32(3):1645–1648, 1996.
8- B. Chauvel, G. Cathelineau, S. Balac, J. Lecerf and J.D. De
Certaines. Cancellation of metalinduced MRI artifacts with
dual-component and diamagnetic material : mathematical modelization and
experimental verification. Journal of Magnetic Resonance Imaging,
6(6)::936-938 (1996)
Codes de simulation numérique associés
à ce thème de recherche
(Suivre les liens pour télécharger ces codes de simulation)
CMAGXX:
une
application Modulef pour le calcul du champ
magnétique induit par un objet uniformément
aimanté.
ARAMIS
: un
programme Matlab de simulation
numérique des artefacts radio-fréquences en IRM pour une géométrie
sphérique.
SIMART :
un
programme Matlab de simulation
numérique des artefacts de susceptibilité
magnétique en IRM pour une géométrie
sphérique.
ATHOS
: un programme Matlab
de simulation numérique des artefacts de
susceptibilité
magnétique en IRM pour une géométrie
cylindrique.