Page 1
The Magnetic Resonance Phenomenon
Mark CohenUCLA Brain Mapping Division
Departments of Psychiatry, Neurology, Radiology, Psychology, Biomedical Physics. Biomedical Engineering
©2007 Mark Cohen, all rights reserved
Nuclear Spin
! “Spin” is a property of many particles. It is a type of angular momentum.
! Angular momentum is a vector quantity.
! Quantum properties prohibit knowing both magnitude and three-dimensional orientation. We can know both the z-component and magnitude.
! = 1.0546 X 10-34 J-s
S = ! s(s +1), where s = {0, 1
2,1, 3
2,2,…}
©2007 Mark Cohen, all rights reserved
Spin States
! A Spin 1/2 particle has two states (“up/
down”, “1 and 2”, !/")
! In a magnetic field, B0, the two states have
different energies
©2007 Mark Cohen, all rights reserved
Proton Precession
Applied Magnetic Field: B
Precession: #
Spin
# = $ X B
$H " 267.52 Rad/sec/Tesla
" 42.577 MHz/Tesla
Page 2
©2007 Mark Cohen, all rights reserved
Protons in Applied Field
Applied Magnetic Field
©2007 Mark Cohen, all rights reserved
Transition to Equilibrium
Zero Field Field Applied
Up/Down state transitions require quantized energy input
Energy
©2007 Mark Cohen, all rights reserved
T1
CSF
Brain
Fat
Time (seconds)
Magnetization(signal)
1
0.5
01 2 3 40
M (t) = M0(1! e
! trT1)
©2007 Mark Cohen, all rights reserved
Protons in Applied Field
Applied MagneticField
Due to their angular momentum, Protons precess in the magnetic field.
Page 3
©2007 Mark Cohen, all rights reserved
Proton Responses to Applied Magnetic Field
! Spin Alignment Along Net Applied Field spins align parallel or anti-parallel to the applied field! Precession About the Magnetic Field at a precession frequency of: $ % B, known as the Larmor frequency! Spin Alignment Occurs at the Rate, T1
©2007 Mark Cohen, all rights reserved
T1:
The Characteristic Time forLongitudinal Relaxation
©2007 Mark Cohen, all rights reserved
the Resonance Phenomenon
When a second magnetic field (B1) is applied, rotating at the Larmor rate, the proton will precess about it.
B1 Field Axis
Precession Angle About RF Field
Static Magnetic Field: B0
©2007 Mark Cohen, all rights reserved
An RF Pulse Converts LongitudinalMagnetization to Signal
90° RF Pulse
LongitudinalMagnetization
MR Signal
Precession
Page 4
©2007 Mark Cohen, all rights reserved
In-Phase Precession
N
S
Receiver
NMR Signal
N
S
©2007 Mark Cohen, all rights reserved
Out of Phase Precession
N
S
Receiver
NMR Signal
©2007 Mark Cohen, all rights reserved
T2
The Characteristic Time forTransverse Decay
©2007 Mark Cohen, all rights reserved
Summary Animation
Page 5
©2007 Mark Cohen, all rights reserved
T2 and TE
0
1
0.5
060 120
Time (milliseconds)
Signal
CSF
Brain
Fat
S(t) = Mxy(t) = M
0e
! teT 2
©2007 Mark Cohen, all rights reserved
Partial Saturation Sequence
NMR Signal
Sequence of 90° Pulses
©2007 Mark Cohen, all rights reserved
Effects of TE at long TR
©2007 Mark Cohen, all rights reserved
Effects of TR (density weighted)
TE=17! NEX=1! Thick=3mmMatrix=256x256 BW=16kHz
Page 6
©2007 Mark Cohen, all rights reserved
Contrast, TR and TE
TR
Long
Short
Short Long
TE
T2-Weighted
T1-Weighted
ProtonDensity
©2007 Mark Cohen, all rights reserved
Contrast, TR and TE
TR
Long
Short
Short LongTE
Density
T1
T2
©2007 Mark Cohen, all rights reserved
Contrast, TR and TE
TR
Long
Short
Short Long
TE
T2-Weighted
T1-Weighted
ProtonDensity
S = k!M0(1" e
" trT1)e
" teT 2
©2007 Mark Cohen, all rights reserved
Contrast, TR and TE
TR
Long
Short
Short LongTE
Density
T1
T2
Page 7
©2007 Mark Cohen, all rights reserved
T2-weighted EPI scan
Metastatic (cancer) lesions, and many others, typically appear bright on T2-weighted MR images
©2007 Mark Cohen, all rights reserved
Liver Mets
TE = 26 TE = 50 TE = 100
©2007 Mark Cohen, all rights reserved
Observed TransverseRelaxation Rate, T2*,is the sum of several components:
the
1
T2*
1
T2
1
T2
1
T2= + +
D"
MolecularFieldInhomogeneity
Diffusion
©2007 Mark Cohen, all rights reserved
MR Formulæ
Contrast Summary:
Spin Echo Signal = k&'0(1 - e-tr/T1)e-te/T2
& is the proton densityk represents instrument effects
The “Bloch” Equation:
dM/dt = $M X B1 + (M0 - Mz)/T1 - (Mx + My)/
T2
Page 8
©2007 Mark Cohen, all rights reserved
Hahn Spin Echo
2. 90° Pulse
3. T2* Relaxation
4. 180° Pulse
5. Spin Rephasing
6. Spin Echo
1. Equilibrium
1. Equilibrium
©2007 Mark Cohen, all rights reserved
Summary Animation
©2007 Mark Cohen, all rights reserved
Inversion Recovery
©2007 Mark Cohen, all rights reserved
Spin Echo
Page 9
©2007 Mark Cohen, all rights reserved
3D T1 Images
3D T1 Images
UCLA BrainMapping Division
TE = 3.2
TR = 14.4
124 slices
1.25 mm thick
1 NEX
Flip Angle 20°
TI = 500
©2007 Mark Cohen, all rights reserved
Sample Data Set (normal)
Sample Data Set (normal)
Fast Spin Echo
3 mm Slices
3D IR-SPGR
TE = 3.2, TI = 700
©2007 Mark Cohen, all rights reserved
Contrast to Noise Ratio
trte
0
0.2
0.13
6
tr, te in seconds
-5% +3%0%
Contrast = [(1! e! tr /1.2
)e! te /.08
], gray matter
![(1! e! tr /1.0
)e! te /.07
], white matter
Gray – White
©2007 Mark Cohen, all rights reserved
Apodization from Long Readouts
Phantom Readout = 2T2* Readout = 4T2*
Page 10
©2007 Mark Cohen, all rights reserved
The Larmor Relation
84
42
840.0 1.0 2.0
Frequency
(MHz)
Magnetic Field (Tesla)
# = $ X B = 2#f
©2007 Mark Cohen, all rights reserved
Magnetic Field Gradients
Position
Field Strength
©2007 Mark Cohen, all rights reserved
Slice
Selection Gradient Axis
©2007 Mark Cohen, all rights reserved
Frequency Encoding
Frequency Gradient Axis
Page 11
©2007 Mark Cohen, all rights reserved
Fourier Projectionss
MR ImageRaw Data
FFT of Raw Data
©2007 Mark Cohen, all rights reserved
realimaginary
Back Projection
Image Domain
Fourier Domain
GradientEncoding
2D FourierTransform
©2007 Mark Cohen, all rights reserved
Equivalent Strategies in k-space*
Gradient
Samples
Gradient
Samples
Gradient
Samples
Time
Gradient
Time
*Ignoring effects of signal decay and sample motion©2007 Mark Cohen, all rights reserved
GradientPre-encoding
Gradient
Samples
Time
Signal
Page 12
©2007 Mark Cohen, all rights reserved
Interleaved Spatial Encoding
Gradient 1
Samples
A2
A1
Gradient 2
Points indicated in
black are affected
by Gradient 1 but
NOT by Gradient 2
Points indicated in Blue are
affected by both gradients
…
…
…
©2007 Mark Cohen, all rights reserved
In Conventional MRISampling is ApportionedInto Brief EpisodesSeparated by a “TR” Period
Scan Time
©2007 Mark Cohen, all rights reserved
Spatial Encoding in a Pulse Sequence
Grad 1
Samples
A2
A1
Grad 2A2
A1
A2
A1
RF
Grad 0
tr
…
…
…
…
…
©2007 Mark Cohen, all rights reserved
k-space
!G(x,y,t)dtk(x,y,t) = $t = 0
T
0 T
GX
GY
ky
kx
t = 0
t = T
Page 13
©2007 Mark Cohen, all rights reserved
Conventional K-Space Trajectory
+K
-Kfrequency
phase
TR
©2007 Mark Cohen, all rights reserved
Echo-Planar k-space Trajectoryk-phase
k-frequency
©2007 Mark Cohen, all rights reserved
Spiral
Gx
Gy
kx
ky
©2007 Mark Cohen, all rights reserved
Equivalent Strategies in k-space*
Gradient
Samples
*Ignoring effects of signal decay and sample motion
Gradient
Samples
These all have equal areas.
Page 14
©2007 Mark Cohen, all rights reserved
A Pulse Sequence Controls
• Slice Location• Slice Orientation• Slice Thickness• Number of Slices• Resolution (FOV and
Matrix)
• Contrast TR, TE, TI, Flip Angle, Diffusion, etc…• Artifact Correction Saturation Pulses, Flow
Comp, Fat Suppression, etc…
©2007 Mark Cohen, all rights reserved
Phase Maps
readout
RF
slice select
te
• Time shift in data collection amounts to a phase offset
• Spins precessing at different rates (different magnetic fields) will acquire different phase shifts
©2007 Mark Cohen, all rights reserved
Imaging System Components
X Y Z
Gradient PowerSystems
RF Transmitter
Magnet RF Receiver
Viewing Console
Scan Controller
©2007 Mark Cohen, all rights reserved
System Components
Page 15
©2007 Mark Cohen, all rights reserved
MR Field Gradient Coil
©2007 Mark Cohen, all rights reserved
Why is MRI So Noisy?
Typical GuitarAmpli!er: 100 Watts
Loudspeaker
~0.2 Tesla
Ultra-fast (echo-planar) gradient Ampli!er: 865,000 Watts
MR Field Gradient!1.5 Tesla
©2007 Mark Cohen, all rights reserved
Gradient Echo Sequence
RF
SliceSelect
PhaseEncode
Readout
©2007 Mark Cohen, all rights reserved
Spin Echo Sequence
RF
SliceSelect
PhaseEncode
Readout
TE
Page 16
©2007 Mark Cohen, all rights reserved
Properties of K-Space
+K
-K
+K
-Kfrequency
phase
frequency
phase
• Increasing K values Represent Higher Resolution
• Finer Grain Sampling Results in Wider FOV
• Reflections Across K=0 are approximate Complex Conjugates
©2007 Mark Cohen, all rights reserved
Raw Data Symmetry
k[t] = k[-t]
k[t] = -k[-t]
Detector 1
Detector 2
©2007 Mark Cohen, all rights reserved
k-space conjugate symmetry
For a Stationary Object,in a Homogeneous Field:
S(kx,ky) = S(kx,-ky)
where S(kx, ky) is the signal at (kx,ky).
Example: if , then .
S(kx,ky) = a + ibS(kx,-ky) = a – ib
©2007 Mark Cohen, all rights reserved
Gradient Coil Characteristics
L # 1 mH
Gradient Coil
i
Gradient Strength = k i
k # 1 Gauss/cm / 100 Amps
Page 17
©2007 Mark Cohen, all rights reserved
Rise Time, Current and Voltage
250 Amps
0.1 msec
For: 250 Amps in 100 µsec VL = 2500 Volts
Power = 2500 Volts x 250 Amps = 6.25 x 105 Watts
di
dt=
VL
1 mHL " 1 mH
i
©2007 Mark Cohen, all rights reserved
Resonant Gradient
R
LC
+
-
VC
iL
VC i
L
+
–(+ –
©2007 Mark Cohen, all rights reserved
Contrast to Noise Ratio
©2007 Mark Cohen, all rights reserved
CNR vs. Resolution
Minimum Imaging Time
256 X 256
128 X 128
64 X 64
Noise free
Page 18
©2007 Mark Cohen, all rights reserved
Imaging Time Held Constant
CNR vs. Resolution
256 X 256
128 X 128
64 X 64
Noise free
256 X 256
128 X 128
64 X 64
©2007 Mark Cohen, all rights reserved
CNR vs. Resolution
Signal/Noise Ratio Held Constant
©2007 Mark Cohen, all rights reserved
An “Equation” in Resolution
Because MR is an emission modality the temporal resolution, spatial resolution and contrast are inter-dependent:
where B0 is the field strength.
Signal = k B voxel size !imaging time
– contrast
0
©2007 Mark Cohen, all rights reserved
Motion Artifact
http://porkpie.loni.ucla.edu/BMD_HTML/SharedCode/Motion/motion.html
Page 19
©2007 Mark Cohen, all rights reserved
Bandwidth and Readout
• Position is encoded by FREQUENCY
• Bandwidth refers to the Frequency Difference
from the center of the image to its edge:
#Frequency per pixel = =
• Bandwidth decreases with readout duration:
Bandwidth =
2* Bandwidth
number of pixels
number of pixels
2 * readout duration
1
readout duration
©2007 Mark Cohen, all rights reserved
Narrow BandwidthWide Bandwidth
Bandwidth and SNR
Decreasing the Bandwidth Improves SNR:
Imaging Time is INCREASED and high frequency noise is excluded
SignalIntensity
FrequencyNoise
©2007 Mark Cohen, all rights reserved
Bandwidth
TE=11-14! NEX=1! Thick=3mmTR=500! Matrix=256x256
BW=4kHz BW=8kHz BW=16kHz
©2007 Mark Cohen, all rights reserved
Bandwidth
TE=11-14NEX=1!ick=3mmTR=500Matrix=256x256
BW=4kHz BW=8kHz BW=16kHz
Page 20
©2007 Mark Cohen, all rights reserved
Shape and Bandwidth
©2007 Mark Cohen, all rights reserved
The Origin of Chemical Shift
In water, electrons move from Hydrogen towards Oxygen.
This exposes the Proton to a slightly higher magnetic field.
Electrons in lipid are shared equally between Hydrogen and Oxygen
Resonance Frequencies
Water Lipid
Higher Frequency
H
OH
H H H
H H H
C C CH
©2007 Mark Cohen, all rights reserved
Chemical Shift Artifact
If the frequency width of each pixel is less than the frequency difference between water and lipid,
then water and lipid will appear in separate pixels
Higher Frequency
©2007 Mark Cohen, all rights reserved
Chemical Shift
d
Water Fat
The Fat-Water chemical shift is about 3.5 ppm or:
Which is: with a 32 kHz readout 75 Hz @ 0.5 Tesla < 1 pixel150 Hz @ 1.0 Tesla ! 1 pixel220 Hz @ 1.5 Tesla > 1 pixel440 Hz @ 3.0 Tesla ! 3.5 pixels
Amplitude
frequency
Lowering the Bandwidth/pixel increases the Chemical Shift in pixels
Page 21
©2007 Mark Cohen, all rights reserved
Aliasing
©2007 Mark Cohen, all rights reserved
Saturation
©2007 Mark Cohen, all rights reserved
Spikes
©2007 Mark Cohen, all rights reserved
k-Space Truncation (Gibbs)
Page 22
©2007 Mark Cohen, all rights reserved
What is the actual resolution of MRI?
Original Data MR Image
Single pixel “activation”©2007 Mark Cohen, all rights reserved
The Actual Resolution of ƒMRI
http://porkpie.loni.ucla.edu/BMD_HTML/SharedCode/MRArtifacts/MRArtifacts.html
©2007 Mark Cohen, all rights reserved
Distortions are More Severe at High Magnetic Field Strength
B0
z
x, y
Variation in sample magnetization of is proportional to !eld strength.
High Field images lose more signal from !eld inhomogeneity
Mid Field High Field
Research supported under DA13054©2007 Mark Cohen, all rights reserved
Crucial Tradeoffs
! Magneto-stimulation, spatial resolution, chemical shift artifact, gradient power and signal to noise ratio are ALL interdependent
Page 23
©2007 Mark Cohen, all rights reserved
EPI Readout Durations
0
0.5
1
0 20 40 60 80 100
T2* signal decay(T2* ~ 45 msec)UCLA 64x128
GE Product 64x128
UCLA 128x128
GE Product 128x128
Stanford Spiral 128x128
MR Signal
©2007 Mark Cohen, all rights reserved
Apodization from Long Readouts
Phantom Readout = 2T2* Readout = 4T2*