Spatial encoding in Magnetic Resonance Imaging Jean-Marie BONNY
What s Qu est an image ce qu une? image? «a reproduction of a material object by a camera or a related technique» Multi-dimensional signal I(r) r = vector of spatial coordinates (x,y,z) T Signal intensity of the image I is scalar = Image with a single frame I is a vector = Multi-channel image
The first Qu est MR image ce qu une image? Made by Paul LAUTERBUR, Nature (1973) Zeumatography comes from the reek word zeugma, or yoke, to signify the fact that the technique links chemical and spatial information
There Qu est are 3 spatial ce qu une dimensions image? In real world, images should be in 3D Human brain V = 3 x 3 x 3 mm 3
1D or 2D Qu est images ce qu une image? The concept of «profile» x Ix, y, zdydz I Mulkern et al, NMR in Biomed icine (1999) PMID: 119533
Multi-channel Qu est ce images qu une? image? MRI provides different contrasts Clark et al, Journal of the Science of Food and Agriculture (1998) DOI: 1.12/(SICI)197-1(199811)78:3<349::AID-JSFA125>3..CO;2-X
The MRI Qu est challenges ce qu une image? More sensitivity Improvement of spatial and temporal resolutions More specificity Tissue discrimination through the signature given by the signal I
How to Qu est encode ce spatial qu une information image?? The target Amplitude M T of transverse magnetization over space M T r r : Spatial Cartesian coordinates r x y z Assumptions No relaxation Uncoupled spins
Signal Qu est without ce encoding qu une image? Free induction decay RF / ACQ Flip angle t
Signal Qu est without ce encoding qu une image? FID S( r, t) MT r, r expi t S Complex signal because of quadrature detection because of reception gain S r, t 2 f r, t' t f dt' Phase accumulated during t period S( r, t) dr M T r dr reception field, amplifiers ( ) On-resonance signal coming from the whole sample
Quadrature Qu est detection ce qu une image? Components of M T in the rotating frame Phase/off-resonance information ReS ImS Hoult, Concepts in Magnetic Resonance (2)
What is Qu est a gradient ce qu une? image? Mathematical operator Vector of partial derivatives with respect to each coordinate For a 3D function fun r funx, y, z gradfun r fun r fun fun x y fun z r r r
Qu est ce qu une image? fun = magnetic field (MF) B Linearly varying magnetic fields x,y,z in each direction is superposed to the static field (B ) B B r r r.... B z y x B B z y x z y x
y Some Qu est examples ce qu une of linearly image varying? MF B r B x. x y. y z. z 3 z x y 1.8 3 2 1-1 -2-3 1.5.6.4.2 -.2 -.4 -.6 1.5 -.8 -.5 -.5-1 y -1-1 -1 -.8 -.6 -.4 -.2.2.4.6.8 1 x x
y Some Qu est examples ce qu une of linearly image varying? MF B r 3 1 z x B x. x y. y z. z y 1.8 4 2-2.6.4.2 -.2-4 1.5 -.4 -.6 1.5 -.8 -.5 -.5-1 y -1-1 -1 -.8 -.6 -.4 -.2.2.4.6.8 1 x x
Qu est ce qu une image? radient of the magnetic field r z y x z y x z y x z y x z z y x B y z y x B x z y x B B......... corresponds to the gradient of the MF The gradient coils generate
radient Qu est coils ce / Principles qu une image? Hidalgo-Tobon, Concepts in Magnetic Resonance (21)
radient Qu est coils ce / How qu une it looks image like??
Effect Qu est of currents ce qu une image? The Biot-Savart law (182) MF generated by a current that circulates through wires db dlr I 3 4 r I = current intensity B I
Time varying Qu est ce fields qu une image? Currents = independant functions of time I t I I I x y z t t t t x y z t t t x y z t
Effects Qu est of MF ce gradients qu une image on NMR? signal Time AND spatially dependent frequency B r B. r Larmor equation B r, t B t r. f 2 2 r, t B t r. Phase of the transverse magnetization t r t 2 f r, t' f dt' t ',. rdt' t
Qu est ce qu une image? The k variable radient antiderivatives t dt t dt t t t t r.k r r r ' '. '. ', t z t y t x z y x t dt t dt t dt t t k t k t k dt t t ' ' ' ' ' ' ' ' k k looks like a «wave vector»
Effects Qu est of MF ce gradients qu une image on NMR? signal Signal coming from the whole sample S T, t S( r, t) dr M r expi r t Using the «k» formulation S T. S t M r expi r kt dr dr 1 k M r expi r. kd r FT M r T T
The k-space Qu est ce qu une image? Rearranged acquired NMR signal Time-domain NMR signal is represented as function of k Reparameterization of acquired signal S k FT FT 1 M T r 1 S k FT FT r M T M r T S k FT M T r k-space FT -1 Image
S k M T r k-space Image k y FT FT -1 y k x x
S k M T r k-space Image k y FT FT -1 y k x x Remember S is a complex number! z is orthogonal to the slide!
How riding Qu est about ce qu une the k-space image? (1) Change the gradients with time t x y z t t t (2) Calculate the wave-vector = integration t t t ' k dt' 3D trajectory
How filling Qu est the ce k-space qu une? image? (1) Change the gradients = displacement along the 3D trajectory (2) Acquire the signal = fill the k-space with the (complex) value of NMR signal
How filling Qu est the ce k-space qu une? image? Many ways Multi-steps / Cartesian TR n.tr and many others!
Pure «phase» encoding
RF ACQ 9 k y x k x y z
RF ACQ 9 k y x k x y z
RF ACQ 9 k y x k x y z
RF ACQ 9 k y x k x y z
RF ACQ 9 k y x k x y z
Pure phase Qu est encoded ce qu une imaging image? Cartesian trajectory over the k-space radients before the acquisition «Phase» encoding gradient pulses Encode the 3 directions Whole FID available All the chemical shift information is available! Acquisition time TA = Nx Ny Nz TR = 64.64.64. 1s = 73h
Spin Warp
RF ACQ 9 k y x k x y z
RF ACQ 9 k y x k x y z
RF ACQ 9 k y x k x y z
Spin warp Qu est ce qu une image? Cartesian trajectory over the k-space radients during the acquisition «Read» gradient pulses Encode one direction «Phase» encoding in the 2 other directions Acquisition time TA = Nx Ny Nz TR = 64.64. 1s = 1h8
Radial encoding
RF ACQ 9 k y x y tanφ= y / x k x z
RF ACQ 9 k y x k x y z
Radial Qu est encoding ce qu une image? Radial trajectory over the k-space Non-uniform density radients during the acquisition Only read-gradient pulses No phase encoding Acquisition time TA = N N TR = 64.64. 1s = 1h8
Echo planar
RF ACQ 9 k y x k x y z
Echo planar Qu est ce qu une image? (Quasi) Cartesian trajectory over the k-space Alternated read gradients during the acquisition Short phase encoding in-between the read gradients Blip Acquisition time TA = TR = 1s
Design of any trajectories t kt t ' dt' t ' t k t Lee, Neuroimage (21)
«Exotic» trajectories Khrapitchev, J Magn Reson (26)
(Some) Limits
Spatial Qu est resolution ce qu une image? Signal-to-noise ratio V = voxel volume Voxel size SNR V x. y. z x kmax x k 1 max High resolution = high = Strong gradients k max
Sampling Qu est ce qu une image? -1-8 -6-4 -2 2 4 6 8 1-1 -8-6 -4-2 2 4 6 8 1 Object 2 4 6 8 1 12 14 Spin warp encoding FOV = 128 mm 16 x 16 16 2 4 6 8 1 12 14 16 Spin warp encoding FOV = 128 mm 64 x 64
Aliasing Qu est ce qu une image? Field of view (FOV) should contain the object Solutions exist for reducing the object size Saturation band Selective pulse Surface coil
Off-resonance Qu est ce / radient qu une image imperfections? Off-resonance Frequency shift r Trajectory perturbation Differences between experimental and theoretical gradient waveforms S T * k t M r exp i r. k t 2 r t dr Trajectory perturbation All the off-resonance perturbations contribute
Off-resonance Qu est ce qu une image? -1-8 -6-4 -2 2 4 6 8 1-1 -8-6 -4-2 2 4 6 8 1 Object 2 4 6 8 1 12 14 Spin warp encoding FOV = 128 mm 16 x 16 No perturbations -5-4 -3-2 16 2 4 6 8 1 12 14 16-1 1 2 3 4 2 5-5 -4-3 -2-1 1 2 3 4 5 Off-resonance 4 6 8 1 With perturbations 12 14 16 2 4 6 8 1 12 14 16
Inhomogeneities Qu est ce qu une of B1 fields image? Quadrature surface coil 7 T Transmission Reception Resulting image Homogeneous phantom Wang et al, Magn Reson Med (22) 48:362-369
Inhomogeneities Qu est ce qu une of B1 fields image? 1.5 Tesla 8 Tesla Novak et al., Magn. Reson. Imaging (25)
Selective Qu est excitation ce qu une image? For reducing the number of directions to encode Band-limited pulse Pulse applied WITH a gradient pulse M T Off-resonance (khz)