Imagerie par résonance magnétique
Pierre Maurel ([Link]@[Link]) Cours IM 32 / 88
Principes physiques
Magnétisme nucléaire
Spin : "Un petit aimant tournant sur lui-même"
Pierre Maurel ([Link]@[Link]) Cours IM 33 / 88
Principes physiques
On s’intéresse aux protons d’hydrogène pour l’imagerie
présent en grande quantité dans le corps (eau)
Pierre Maurel ([Link]@[Link]) Cours IM 34 / 88
Principes physiques
On s’intéresse aux protons d’hydrogène pour l’imagerie
présent en grande quantité dans le corps (eau)
Aimant très intense pour aligner les moments magnétiques
Pierre Maurel ([Link]@[Link]) Cours IM 34 / 88
Principes physiques
Ordre de grandeur
à température ambiante et à 1,5T (Tesla) environ 9,6 spins sont en excès
dans la position parallèle par million de protons
Pierre Maurel ([Link]@[Link]) Cours IM 35 / 88
Principes physiques
Ordre de grandeur
à température ambiante et à 1,5T (Tesla) environ 9,6 spins sont en excès
dans la position parallèle par million de protons
6 × 1022 protons par ml de tissu
Pierre Maurel ([Link]@[Link]) Cours IM 35 / 88
Principes physiques
Ordre de grandeur
à température ambiante et à 1,5T (Tesla) environ 9,6 spins sont en excès
dans la position parallèle par million de protons
6 × 1022 protons par ml de tissu
le nombre de spins contribuant au signal RMN à 1,5T est
9.6 × 10−6 · 6 × 1022 = 5 × 1017 par ml de tissu
Pierre Maurel ([Link]@[Link]) Cours IM 35 / 88
Principes physiques
Perturbation de l’aimantation
Pierre Maurel ([Link]@[Link]) Cours IM 36 / 88
Principes physiques
Perturbation de l’aimantation
Pour écarter M de l’équilibre, il faut appliquer une OEM à la fréquence de
Larmor :
γ
ν0 = B0
2π
où γ est le rapport gyromagnétique du noyau (constante physique
caractéristique de chaque noyau possédant un spin non nul).
fréquence de Larmor = fréquence de précession
"fréquence de résonance du système"
Pierre Maurel ([Link]@[Link]) Cours IM 36 / 88
Principes physiques
Pierre Maurel ([Link]@[Link]) Cours IM 37 / 88
Principes physiques
Pierre Maurel ([Link]@[Link]) Cours IM 38 / 88
Principes physiques
Retour à l’équilibre le long de B0 : Équations de Bloch
dMx (t) Mx (t)
= γ(M(t) × B(t))x −
dt T2
dMy (t) My (t)
= γ(M(t) × B(t))y −
dt T2
dMz (t) Mz (t) − M0
= γ(M(t) × B(t))z −
dt T1
Pierre Maurel ([Link]@[Link]) Cours IM 39 / 88
Principes physiques
Retour à l’équilibre le long de B0 : Équations de Bloch
dMx (t) Mx (t)
= γ(M(t) × B(t))x −
dt T2
dMy (t) My (t)
= γ(M(t) × B(t))y −
dt T2
dMz (t) Mz (t) − M0
= γ(M(t) × B(t))z −
dt T1
Relaxation = retour à l’équilibre de l’aimantation
Relaxation transversale en T2 Relaxation longitudinale en T1
Pierre Maurel ([Link]@[Link]) Cours IM 39 / 88
Principes physiques
Mesure du signal RMN
Mxy tend vers 0 tandis que Mz croît.
Pierre Maurel ([Link]@[Link]) Cours IM 40 / 88
Principes physiques
Mesure du signal RMN
Mxy tend vers 0 tandis que Mz croît.
la décroissance de Mxy se fait de manière exponentielle et engendre un
courant induit dans une bobine (ou antenne de réception) située sur l’axe
Oy.
Pierre Maurel ([Link]@[Link]) Cours IM 40 / 88
Principes physiques
Mesure du signal RMN
Mxy tend vers 0 tandis que Mz croît.
la décroissance de Mxy se fait de manière exponentielle et engendre un
courant induit dans une bobine (ou antenne de réception) située sur l’axe
Oy.
C’est ce signal qui est mesuré !
T2 qui caractérise la décroissance dépend du tissu → on obtient une
image avec des intensités différentes selon les tissus observés
Pierre Maurel ([Link]@[Link]) Cours IM 40 / 88
Principes physiques
Problème : d’où vient le signal enregistré par la bobine ?
Pierre Maurel ([Link]@[Link]) Cours IM 41 / 88
Principes physiques
Problème : d’où vient le signal enregistré par la bobine ?
Par défaut : de partout !
Pierre Maurel ([Link]@[Link]) Cours IM 41 / 88
Principes physiques
Problème : d’où vient le signal enregistré par la bobine ?
Par défaut : de partout !
Sélection de coupe
on ajoute un "gradient" au champ magnétique principal B0 qui dépend donc
maintenant de z
B0 (z) = B0 + Gz .z
Pierre Maurel ([Link]@[Link]) Cours IM 41 / 88
Principes physiques
Problème : d’où vient le signal enregistré par la bobine ?
Par défaut : de partout !
Sélection de coupe
on ajoute un "gradient" au champ magnétique principal B0 qui dépend donc
maintenant de z
B0 (z) = B0 + Gz .z
la fréquence de Larmor (fréquence résonance) est donc modifié :
γ
ν( z) = B0 (z)
2π
Pierre Maurel ([Link]@[Link]) Cours IM 41 / 88
Principes physiques
Problème : d’où vient le signal enregistré par la bobine ?
Par défaut : de partout !
Sélection de coupe
on ajoute un "gradient" au champ magnétique principal B0 qui dépend donc
maintenant de z
B0 (z) = B0 + Gz .z
la fréquence de Larmor (fréquence résonance) est donc modifié :
γ
ν( z) = B0 (z)
2π
en présence d’un gradient Gz la fréquence de résonance dépend de z
Pierre Maurel ([Link]@[Link]) Cours IM 41 / 88
Principes physiques
Problème : d’où vient le signal enregistré par la bobine ?
Par défaut : de partout !
Sélection de coupe
on ajoute un "gradient" au champ magnétique principal B0 qui dépend donc
maintenant de z
B0 (z) = B0 + Gz .z
la fréquence de Larmor (fréquence résonance) est donc modifié :
γ
ν( z) = B0 (z)
2π
en présence d’un gradient Gz la fréquence de résonance dépend de z
Le choix de la fréquence d’excitation RF permet de sélectionner une coupe
Pierre Maurel ([Link]@[Link]) Cours IM 41 / 88
Principes physiques
grâce au gradient de coupe (B0 n’est plus constant), on ne "bascule" que
les spins dans une coupe choisie
Pierre Maurel ([Link]@[Link]) Cours IM 42 / 88
Principes physiques
grâce au gradient de coupe (B0 n’est plus constant), on ne "bascule" que
les spins dans une coupe choisie
mais on enregistre alors un signal électrique (1D) donné par la relaxation
de tous les spins de cette coupe (2D)
Pierre Maurel ([Link]@[Link]) Cours IM 42 / 88
Principes physiques
grâce au gradient de coupe (B0 n’est plus constant), on ne "bascule" que
les spins dans une coupe choisie
mais on enregistre alors un signal électrique (1D) donné par la relaxation
de tous les spins de cette coupe (2D)
...
Pierre Maurel ([Link]@[Link]) Cours IM 42 / 88
Principes physiques
grâce au gradient de coupe (B0 n’est plus constant), on ne "bascule" que
les spins dans une coupe choisie
mais on enregistre alors un signal électrique (1D) donné par la relaxation
de tous les spins de cette coupe (2D)
...
Codage en phase (gradient Gy )
Codage en fréquence (gradient Gx )
Pierre Maurel ([Link]@[Link]) Cours IM 42 / 88
Principes physiques
grâce au gradient de coupe (B0 n’est plus constant), on ne "bascule" que
les spins dans une coupe choisie
mais on enregistre alors un signal électrique (1D) donné par la relaxation
de tous les spins de cette coupe (2D)
...
Codage en phase (gradient Gy )
Codage en fréquence (gradient Gx )
Espace k (pour une coupe)
Pierre Maurel ([Link]@[Link]) Cours IM 42 / 88
Rappel : Transformée de Fourier
1 2 3 4 5
A B C D E
Pierre Maurel ([Link]@[Link]) Cours IM 43 / 88
Rappel : Transformée de Fourier
1 2 3 4 5
A B C D E
1→B 2→D 3→A 4→C 5→E
Pierre Maurel ([Link]@[Link]) Cours IM 43 / 88
Dispositif d’acquisition
Pierre Maurel ([Link]@[Link]) Cours IM 44 / 88
IRM anatomique
variété de séquence/contraste
bonne résolution (< 1 mm)
Pierre Maurel ([Link]@[Link]) Cours IM 45 / 88
IRM anatomique
Contraste T1
substance blanche → blanche
substance grise → grise
LCR (ventricules) → noir
Comprendre l’IRM : Manuel d’auto-apprentissage (Bruno Kastler, Daniel Vetter, Zoltan Patay, Philippe Germain)
Pierre Maurel ([Link]@[Link]) Cours IM 46 / 88
IRM anatomique
Contraste T2
Contraste densité de protons ρ
Comprendre l’IRM : Manuel d’auto-apprentissage (Bruno Kastler, Daniel Vetter, Zoltan Patay, Philippe Germain)
Pierre Maurel ([Link]@[Link]) Cours IM 47 / 88
IRM anatomique
pathologies
Tumeur cérébrale Sclérose en plaques
[Clatz, TMI 2005]
Pierre Maurel ([Link]@[Link]) Cours IM 48 / 88
IRM anatomique
pathologies
Lésion cérébrale : hyposignal en contraste T1 et hypersignal en contraste T2
Comprendre l’IRM : Manuel d’auto-apprentissage (Bruno Kastler, Daniel Vetter, Zoltan Patay, Philippe Germain)
Pierre Maurel ([Link]@[Link]) Cours IM 49 / 88
IRM anatomique
imagerie foetale et néonatale
Dubois et al. Cerebral Cortex, 2007
Pierre Maurel ([Link]@[Link]) Cours IM 50 / 88
IRM fonctionnelle
signal BOLD (blood-oxygen-level dependent)
zones "activées" (dont le cerveau se sert) → afflux de sang oxygénés
Concentration d’oxyhémoglobine plus forte près des régions activées
Pierre Maurel ([Link]@[Link]) Cours IM 51 / 88
IRM fonctionnelle
Création de carte statistique et fusion avec l’IRM anatomique (pour
visualisation)
Aires plus activées lors de l’écoute de phrases en langue française que par
des phrases en langue étrangère
Pallier et al., Cerebral Cortex, 2003
Pierre Maurel ([Link]@[Link]) Cours IM 52 / 88
IRM de diffusion
Diffusion des molécules d’eau (et donc protons d’hydrogène)
Diffusion restreinte Diffusion restreinte
Diffusion libre
isotrope anisotrope
Pierre Maurel ([Link]@[Link]) Cours IM 53 / 88
IRM de diffusion
en chaque point estimation de la probabilité de diffusion en chaque
direction → Tenseur de diffusion
Modèle gaussien de diffusion : un tenseur représenté par une matrice
symétrique 3 × 3
Pierre Maurel ([Link]@[Link]) Cours IM 54 / 88
Problématiques classiques en IRM
IRM anatomique
Segmentation des différents tissus
Pierre Maurel ([Link]@[Link]) Cours IM 55 / 88
Problématiques classiques en IRM
IRM anatomique
Segmentation des différents tissus
Recalage inter-modalités
Pierre Maurel ([Link]@[Link]) Cours IM 55 / 88
Problématiques classiques en IRM
IRM fonctionnel
détection des zones d’activation cérébrale
Pierre Maurel ([Link]@[Link]) Cours IM 56 / 88
Problématiques classiques en IRM
IRM fonctionnel
détection des zones d’activation cérébrale
IRM de diffusion
tractographie des faisceaux de fibres de la matière blanche
Pierre Maurel ([Link]@[Link]) Cours IM 56 / 88
Problématiques classiques en IRM
d’autres variantes : ASL (Arterial Spin Labelling) pour étudier la perfusion
Pierre Maurel ([Link]@[Link]) Cours IM 57 / 88
IRM : récapitulatif
Résolution spatiale
anatomique : 1 mm
fonctionnelle : 3 mm
diffusion : 1 à 3 mm
Pierre Maurel ([Link]@[Link]) Cours IM 58 / 88
IRM : récapitulatif
Résolution spatiale
anatomique : 1 mm
fonctionnelle : 3 mm
diffusion : 1 à 3 mm
Résolution temporelle
anatomique : 15 min
fonctionnelle : 1s
diffusion : 10 min à 1h
Pierre Maurel ([Link]@[Link]) Cours IM 58 / 88
IRM : récapitulatif
Résolution spatiale
anatomique : 1 mm
fonctionnelle : 3 mm
diffusion : 1 à 3 mm
Résolution temporelle
anatomique : 15 min
fonctionnelle : 1s
diffusion : 10 min à 1h
Améliorations permanente de ces résolutions
nécessité d’aimants de plus en plus puissants (nocivité pour l’homme ?)
Pierre Maurel ([Link]@[Link]) Cours IM 58 / 88