Spectroscopie Séquence de RMN
180°
90°
temps
t
5 mm
500 µL
Tube d’eau
RMN
haute résolution fréquence
Spectre 1H 42
Représentation dans les domaines
temporel : fréquentiel :
gaussienne
t
T période en (s), f (Hz)= 1/T = 2 f, pulsation en rd/s
?
43
Représentation dans les domaines
temporel : fréquentiel :
gaussienne
T période en (s), f (Hz)= 1/T = 2 f, pulsation en rd/s
43b
Représentation dans les domaines
temporel : fréquentiel :
?
T période en (s), f (Hz)= 1/T = 2 f, pulsation en rd/s
?
44
Représentation dans les domaines
temporel : fréquentiel :
T période en (s), f (Hz)= 1/T = 2 f, pulsation en rd/s
44b
Représentation dans les domaines
temporel : fréquentiel :
?
T période en (s), f (Hz)= 1/T f, fréquence
45
Représentation dans les domaines
temporel : fréquentiel :
+ f, fréquence
T période en (s), f (Hz)= 1/T 45b
Spectroscopie
4.0 3.5 3.0 2.5 2.0 [ppm]
10 ppm ( 5000 Hz @ 500 MHz)
5 mm
500 µL Spectre 1H
RMN Tube ?
haute résolution
46
Pourquoi tous les protons ne résonnent-ils
pas exactement à la même fréquence ?
Deux origines :
Déplacement chimique
Couplage scalaire
4.0 3.5 3.0 2.5 2.0 [ppm]
f1 f2 …
47
Déplacement chimique
B0
B0 b0
b0
Le noyau observé e
(proton par exemple)
ne « voit » pas un champ B0,
mais un champ plus faible
car la rotation du nuage électronique autour du noyau
(= un courant)
vient créer un champ magnétique b0 qui s’oppose à B0.
48
Déplacement chimique : table des décalages
TMS
51
Applications médicales : spectroscopie localisée in vivo
Spectres avec suppression
d’eau (foie droit de 2 patients)
à 1,5 T (HEH, Lyon)
Voxel: 2 x 2 x 2 cm3
stéatose
55
De la RMN à l’IRM
Imagerie par Résonance Magnétique
56
Espace fréquence x
S(t) S(t)
t t
S(f)
f f
f1 f2 57
Transformée de Fourier x
S(t)
Temps t
S(f) Fourier (1768-1830)
Fréquence
f (ou x !)
f1 f2 58
Imagerie procédé de localisation spatiale ?
idée 1: gradient de champ magnétique superposé à B0
B (T)
B0 + Gx x
B0
Gradient G = B/x
variation linéaire
qq 10 mT/m
x 59
Imagerie procédé de localisation spatiale ?
idée 2: F B0 [RMN] (F = 42,5 MHz , B0 = 1 tesla)
B (T) F (Hz)
B0 + Gx x1 F1=F0 + Gx x1
B0 F0
Gx (Hz/T)
Fréquence de résonance :
variation linéaire avec la coordonnée spatiale
x
x1 60
Idée 3: 3 directions (x, y, z) 3 gradients
z
x
x
Exemple : Gx
Bo
y 61
Imagerie procédé de localisation spatiale ?
idée 3: 3 directions (x, y, z) 3 gradients
Gradients
Gx = B/x,
Gy = B/y,
Gz = B/z,
produits par
3 jeux de bobines
parcourues par des
courants. 62
Codage par la fréquence :
lecture du signal en présence d’un gradient
B (T) objet B0 + G y
B0 B0
y L y
Signal à la fréquence F0 Dispersion en fréquence
sur F = G L
(Hz/T) 63
Codage par la fréquence :
lecture du signal en présence d’un gradient
= 90° = 90°
B1 B1 Gy
t t
Mxy Mxy
t t
64
Codage par la fréquence :
Mxy Mxy
t t
Analyse en fréquence du signal transformation de Fourier
Mxy
f ou y
Profil (projection) de l’objet suivant la au gradient
67
Codage par la fréquence : illustration.
= 90°
B1
Gy= 10 mT/m Gy
t
y
Mxy Fourier
f
t
ou
y= 0,1m f = ? y
Pour le proton 42,5 106 Hz/T, ici 10 mT/m x 0,1 m 1 mT f = 42,5 kHz
74
Codage par la fréquence :
Objet de largeur L=0,2 m, gradient G = 15 mT/m
Objet
de largeur L B0 + G y
B0
f
L y
F ?
Calculer la dispersion en fréquence F
Données: rapport gyromagnétique du proton 42,5 MHz/T
75
Codage par la fréquence :
Objet de largeur L=0,2 m, gradient G = 15 mT/m
Données: rapport gyromagnétique du proton 42,5 MHz/T
75b
Codages dans le plan image (2D)
Transformation
de Fourier 2D
Ny . Nz points
76
Codages dans le plan image (2D)
Ny . Nz points
Transformation
de Fourier
Inverse 2D
77
Codages dans le plan image (2D)
centre
(24x24 points)
Basses fréquences
spatiales
périphérie
Hautes fréquences
spatiales
Espace k Espace réel 78
L’affection des gradients (x, y, z)
est totalement arbitraire (même en oblique).
79
L’affection des gradients (x, y, z)
est totalement arbitraire (même en 3D).
80
Propriétés magnétiques du noyau : 1H principalement
1 Polarisation dans un champ magnétique field (teslas)
3 Codage (3 gradients)
4 Relaxation
5 signal
6 Image !
2 Excitation radiofréquence
(T, MHz)
@ 1.5 T : 64 MHz 81
IRM : le signal et le bruit #110
Signal en RMN/IRM
•Loi de Boltzman:
n-/n+ = e-E/kT = e- h B0 /kT ~ 1 - h B0 /kT
Et polarisation M0 ~ n/n ~ + h B0 /kT
S ~ M0 donc ~ (B0 ou 0)
Et capteur de flux e=-d/dt ~ 0
Donc, au total :
Signal de RMN, S ~ 02 ou 02
82
IRM : le signal et le bruit
Le rapport signal sur bruit dans l’image:
1. Proportionnel à la quantité de noyaux présents
dans le volume élémentaire observé Vv (voxel)
SB = k Vv h e
Vv = l h e l
83
IRM : le signal et le bruit
Le rapport signal sur bruit dans l’image:
2. Racine du nombre d’observations du signal
SB = k No et le temps d’acquisition Tacq = No TR
Exemple : signal=1, écart type du bruit B=0,1
S/B~10 S/B~20 S/B~40
No :1 4 16
Le signal linéairement
Le bruit (incohérent) avec la racine de No. 84
IRM : le signal et le bruit
Le rapport signal sur bruit dans l’image:
SB = k Vv No
et le temps d’acquisition de l’image
Tacq = No TR
No : le nombre de fois que l’on répète l’expérience
No = Ncodage1 . Ncodage2 . Nacc # accumulations
TR : temps de répétition de la séquence 85
Deux images acquises dans des conditions différentes
Déterminer : le gain (ou la perte) en signal/bruit de
l’image 2 par rapport à l’image 1
Données: matrice 256x256,
N accumulations 2 1
épaisseur de coupe: 3 mm 6 mm
Champ de vue: 300 x 300 mm2 250 x 250 mm2
1 2
86
Deux images acquises dans des conditions différentes
Gain (ou perte) en signal/bruit de 2 / 1 SB = k Vv No
86b
Une acquisition est réalisée avec un volume de voxel de
1,25x1,25x1,25 mm3 en 5 minutes.
Le rapport signal sur bruit est jugé correct, mais
la résolution est jugée médiocre et un volume de voxel de
1x1x1 mm3 est choisi pour une seconde acquisition.
Déterminer la durée d’acquisition désormais
nécessaire pour conserver le signal sur bruit de la
première acquisition.
87
Une acquisition est réalisée avec un volume de voxel de
1,25x1,25x1,25 mm3 en 5 minutes.
Le rapport signal sur bruit est jugé correct, mais
la résolution est jugée médiocre et un volume de voxel de
1x1x1 mm3 est choisi pour une seconde acquisition.
Durée d’acquisition nécessaire à signal sur bruit constant ?
87b
Notion de contraste sur bruit
image Signal1 =20
profil Signal2 =15
1
2
Bruit =1
C/B = |S1-S2|/B= |S/B1-S/B2| S/B1= 20
S/B2= 15
SB = k Vv No
88
Contraste sur bruit
Signal1 =20
S/B1= 20 Signal2 =15
S/B2= 15
Bruit =1
C/B =|S1-S2|/B= |S/B1-S/B2 |= 5
SB = k Vv No
88b
Contraste sur bruit
Deux tissus de même densité de proton et de temps de relaxation respectifs
{T1a=1s, T2a=0,1s} et {T1b=0,5s, T2b=0,08s}
sont examinés à l’aide de deux séquences d’écho de spin de paramètres
{TRx=1s, TEx=0,02s} puis {TRy=4s, TEy=0,1s}.
Le bruit est mesuré égal à 0,005 pour les deux séquences.
1. Calculer les contrastes sur bruit Cx et Cy obtenus entre les deux tissus a et b :
|Sa-Sb|/bruit pour les séquences x et y.
2. En prenant en compte le temps d’obtention des deux images x et y, comparer
les deux séquences en terme de contraste sur bruit par unité de temps. Justifier
précisément la réponse.
C/B = |S1-S2|/B= |S/B1-S/B2| SB = k Vv No
89
Contraste sur bruit
89b
IRM : l’instrument en quelques vues
Aimant
RF Gradients
Table
Calculateur
2
IRM, les composants : AIMANT, technologie
I : courant, n tours de fil
Rayon: a
B0
z
T
8 km de fil !
Exemple: 1.5
1.25
a= 0,5 m 1
I=500 A
0.75
0.5
B0/B0 !!!
n= 2500 tours 0.25
0
-1 - 0.5 0 0.5 1 m 3
IRM, les composants : AIMANT, technologie
Aimants résistifs et supraconducteurs : uniformité
Paire I, n I, n
d’Helmholtz,
ou solénoïde z
B0
T
Exemple: 1.4
Pas à
1.2
a= 0,5 m 1 l ’échelle
I=500 A 0.8
0.6
n= 1640 tours 0.4
0.2
zc=a/2 4
- 0.75 - 0.5- 0.25 0 0.25 0.5 0.75 1
m
IRM, les composants : AIMANT, technologie
Du fil supraconducteur (NbTi) à l’aimant...
5
IRM, les composants : AIMANT, technologie
fil supraconducteur (NbTi) + He liquide : 1,5 à 3 T
6
IRM, les composants : AIMANT, technologie
Du cryostat au … fil.
7
IRM, les composants : gradients
IRM standards :
G
G maximum : 10 à 60 mT/m
T
Temps de commutation T :
moins d’une ms (G Max.)
Pente (slew rate) : G/T en T/m/s
20 à 200
Linearité :
qq % sur le volume d’imagerie
8
IRM, les composants : gradients
Système de gradient : jeu complet
I ~ 50 A G ~ 10 mT/m 9
IRM, les composants : gradients
Bruit acoustique ?
Haut parleur:
Couplage
Aimant mécanique, bobine, aimant
Bobines
de
Gradient
Courant I(t)
Structure G(t)
i(t) Temps 10
IRM, les composants : RF émission
Fonction B0
Basculement de M0 grâce
à un champ magnétique = B1
radiofréquence B1.
(B1 B0)
B1
Contraintes:
Accès du patient et sécurité
11
IRM, les composants : RF émission
Bobine
Emetteur RF :
radiofréquence :
V
t
Puissance RF : 1 à 15 kW
12
IRM, les composants : RF réception
Fonction B0
mesure du signal de RMN
(bobine radiofréquence)
F0
+ M0
amplification, démodulation
(spectromètre)
fem
Contrainte force électromotrice (tension)
maximiser le rapport signal sur bruit (S/B) !
13
IRM, les composants : RF réception
Bobine de surface ?
14
IRM, les composants : RF réception
Illustration : utilisation de 8 canaux de réception
Eau dopée (Gd)
Eau pure
2 canaux « spine ».
réseau 6 canaux BA
15
IRM, les composants : RF réception #140
Exemple
Siemens 1.5 T Sonata, avec « 6 channel Body Anterior » 6BA.
16
IRM, les composants : RF réception #141
Exemple
Siemens 1.5 T Sonata, avec « 6 channel Body Anterior » 6BA.
17
IRM, les composants : RF réception
Images canal 1 à 8 :
Combinaison :
18
IRM : applications
Parc d’appareils IRM cliniques :
France : > 800 (~1100 CT)
Monde : ~ 50 000
Coûts approximatifs:
IRM (acquisition) : 1 à 2 M€
Antenne externe : 10 à 40 k€
Antenne usage unique : 150 €
Produit de contraste : 100 €
Examen (technique + médical): 300 €
19
IRM : anatomie, diagnostic
Corps entier rapide
20
IRM : anatomie, diagnostic
Rachis cervical
TSE 3 mm TSE 3 mm
21
IRM : anatomie, diagnostic
Angiographie des membres inférieurs
22
IRM : imagerie fonctionnelle cérébrale
Stimulation
activation neuronale
débit sanguin
HbO2 / Hb + paramagnetique
T2*
Signal 23
1%-10%FMRIB, Oxford University, Department of Clinical Neurology
IRM : anatomie, diagnostic
Foie : surcharge en fer => T2*
TE = 4ms TE = 4ms
TE = 14ms
TE = 14ms
Normal Forte surcharge
< 40 µmol/g 300 µmol/g
24
• RMN des tissus biologiques : + flux (sang),
diffusion,
contraste très important, eau/graisse,
méthode ~ inoffensive. température, ...
• IRM : localisation spatiale de l’information très efficace
et souple (gradients de champ magnétique).
Fin 161
36
