10. MR Physics & Imaging
MRI is non-ionizing
MRI can have both anatomical and functional imaging
MRI uses precession imaging
NMR is nuclear magnetic resonance, but is not related to radioactive decay (the nucleus is spinning, but there’s no radioactivity)
MRI offers images of the same object in different contrasts to bring attention to different areas
Protons and neutrons are spinning inside the atom, which act like magnets
Spinning nuclei in an atom must have an odd atomic number or odd mass number to exhibit this behavior
If a nuclei has an even atomic number, the magnetic fields cancel out
Spinning nuclei possess angular momentum (J) called spin
Placed in a magnetic field, the spin of the protons and neutrons assume certain orientations
Each spin produces a microscopic magnetization vector (𝜇)
𝜇 = 𝛾 * J where gamma is the gyromagnetic ratio (specific to each material)
Gyromagnetic ratio gamma is radians per second*Teslas
The gyromagnetic ratio can also be gamma divided by 2pi
The body is mostly water, which has a lot of hydrogen and is a great source for MRI signal
Hydrogen has the highest gyromagnetic ratio
The net (macroscopic) magnetization of a magnetic field is zero in its natural state
A single particle will have its own nonzero magnetic field
Once an external magnetic field (B0) is applied…
The spins of each particle align in parallel or anti-parallel (a semi majority in parallel)
The net (bulk) magnetization vector with an external magnetic field is nonzero
The bulk magnetization (𝜇 vector) vector is the source of the MR signal (also called M vector)
B0 is always applied along the z direction (up-down on the body)
Field strengths in MRI are commonly 1.5 T or 3 T (Tesla)
A 1.5T magnet is about 30,000 times stronger than Earth’s magnetic field
RF cage has to be in the MRI room to absorb any signals from outside, otherwise FM radio frequencies would be measured
The vector 𝜇 precesses (spins) around the z axis st an angular frequency 𝜔 (radians per second)
𝜔0 = 𝛾 * B0
Larmor Frequency: The angular frequency at which the bulk magnetization factor precesses around the z-axis
f = 𝛾 / (2pi) * B0
Larmor frequency depends only on the material (𝛾) and the main magnetic field strength (B0)
Spin Equilibrium & spin excitation
At equilibrium, the net magnetization vector M precesses about the z-axis
When M is entirely along the z-axis, only the Mz (longitudinal) component exists, and the transverse component M(xy) = 0
If M tips away from the z-axis (during excitation), a component in the x-y plane is generated
The signal is always measured in the x-y plane
The signal component is always in the transverse x-y plane
RF excitation
An RF pulse moves M away from the equilibrium state (tips away from z-axis), by using energy from the RF pulse onto the magnetic field B1
B0 is the magnet by itself, and the magnet is always on
The coil with current running through it acts as the magnet, and is surrounded by liquid Helium to keep the coil cool (decreasing resistance and energy lost to heat)
Why do we want the maximum signal? To decrease noise
B1 is applied from external antennas
B1 is along on the x-axis, so B1x is the only non-zero value for B1
Resonance comes from B1 (the additional force) being at the same frequency as B0 (𝜔0) to make sure there is maximum energy being delivery
As M returns to equilibrium, an RF signal is produced by energy releasing from the absorbed energy that was given to the material from B1
RF Pulse Shape
While oscilating at 𝜔0, the RF pulse can have three shapes, rectangular, sinc, and triangular
Rotating frame of reference
Since the spins and B1 field are both rotating, you can image from a spinning reference
The transverse plane (xy) is rotating at 𝜔, which is equal to 𝜔0, so it can be considered a rotating plane (x’-y’)
Flip angle/tip angle
The flip angle depends on the (1) shape of the RF pulse, (2) field strength B1 and (3) duration of the RF pulse 𝜏p (tau p)
For a rectangular pulse after 𝜏 seconds, the vector M has rotated at an angle α
α = 𝛾 * *B1 ** 𝜏
For any B1(t)
α (always in radians) is the integration from 0 to 𝜏 over B1(t) with respect to t (in seconds)
For the max transverse component Mxy, α is set to pi/2 radians
α = 𝛾 * *B1 ** 𝜏 = pi/2
Relaxation: The process by which M returns to steady state configuration after a B1 is done being applied
The longitudinal Mz and transverse Mxy components vary as time changes, and are descibed by Bloch Equations
Where T1 is the longitudinal relaxation time (spin minus lattice relaxation)
T2 is the transvese relaxation time (spin minus spin relaxation)
The solutions become
Mz^0 is the steady state of Mz
Mz(0+) is transverse magnetization Mx’y’ immediately after the RF pulse
We always measure the x-y plane and do everything else after
T1 is usually a LOT larger that T2
The longitudinal time is usually a lot longer than the transverse time (around 10 times more)
Mx’y’ decays at a faster rate (T2* instead of T2) because of signal desync
Spin Echo
To correct for dephasing and loss of Mx’y’ signal, another RF pulse (180 degree pulse) is applied to re-phase the spins
Spin echo makes the slower spins catch up with faster oens to sync the signal
Relaxation times T1 and T2 are dependant on materials (very high for water, really low for muscle, fat, tendons)
Contrast in MRI is the difference between signal intensities generated by different tissues (B0, which is constant, vs gamma, which is material specific)
The difference in signal intensity is used to discriminate different tissues by representing them as brightness
Pulse sequence: Carefully timed set of scanner operations used to generate images.
Gradient-echo based pulse sequences are
Spin-echo based pulse sequences
90 degree pulse
180 degree pulse
Wait for relaxation echo
Record relaxation echo
Repeat
MR signal is always recorded in the xy plane
Signal intensity S of a SE sequence is
S = K x [H] x (1 - exp(-1 x TR / T1) x exp(-1 x TE/T2)
K is the scaling factor
[H] is the proton/spin density
Twice the number of spins means you have twice as large a signal
TR is the repetition time
TE is the echo time
T1 is the longitudinal relaxation time
T2 is the transverse relaxation time
Weighting
Proton density weighting: We want a long TR and a short TE
T1 weighting: We want an intermediate TR and a short TE (becomes a constant scaling factor, less noise)
T2 weighting: We want ta long TR (minimize T1 differences) and an intermediate TE
Gradients
Adding a small gradient field (extra magnetic field) allows the same material to be differentiated if it is at different distances
Gradient affects Larmor frequency
Gradients added using more coils
Gradients are localized in x, y, and z directions
Gz is a slice selection
z-directional coils are called Maxwell pairs
Gx is frequency encoding
x-directional coils are called Golay coils
Gy is phase encoding
y-directional coils are called Golay-type coils
Gradients vary the magnetic field strength along the z-direction, which affects the frequency and phase of the spins
Fourier Transform is used to break a signal into its freuqency components (sum of sinusoids) based on their amplitude and phase
Slice thickness
Centered at z=0 , B(z) = B0 + Gz(z)
Gz is gradient strength
The magnet is used to create polarization (creates the bulk of the magnetic field)
The RF oil is used to promote excitation (sends RF energy at resonance conditions)
The gradients are used to generate spatial localization (used to form the images)
The receiver RF coils are there to target specific organs
MRI is non-ionizing
MRI can have both anatomical and functional imaging
MRI uses precession imaging
NMR is nuclear magnetic resonance, but is not related to radioactive decay (the nucleus is spinning, but there’s no radioactivity)
MRI offers images of the same object in different contrasts to bring attention to different areas
Protons and neutrons are spinning inside the atom, which act like magnets
Spinning nuclei in an atom must have an odd atomic number or odd mass number to exhibit this behavior
If a nuclei has an even atomic number, the magnetic fields cancel out
Spinning nuclei possess angular momentum (J) called spin
Placed in a magnetic field, the spin of the protons and neutrons assume certain orientations
Each spin produces a microscopic magnetization vector (𝜇)
𝜇 = 𝛾 * J where gamma is the gyromagnetic ratio (specific to each material)
Gyromagnetic ratio gamma is radians per second*Teslas
The gyromagnetic ratio can also be gamma divided by 2pi
The body is mostly water, which has a lot of hydrogen and is a great source for MRI signal
Hydrogen has the highest gyromagnetic ratio
The net (macroscopic) magnetization of a magnetic field is zero in its natural state
A single particle will have its own nonzero magnetic field
Once an external magnetic field (B0) is applied…
The spins of each particle align in parallel or anti-parallel (a semi majority in parallel)
The net (bulk) magnetization vector with an external magnetic field is nonzero
The bulk magnetization (𝜇 vector) vector is the source of the MR signal (also called M vector)
B0 is always applied along the z direction (up-down on the body)
Field strengths in MRI are commonly 1.5 T or 3 T (Tesla)
A 1.5T magnet is about 30,000 times stronger than Earth’s magnetic field
RF cage has to be in the MRI room to absorb any signals from outside, otherwise FM radio frequencies would be measured
The vector 𝜇 precesses (spins) around the z axis st an angular frequency 𝜔 (radians per second)
𝜔0 = 𝛾 * B0
Larmor Frequency: The angular frequency at which the bulk magnetization factor precesses around the z-axis
f = 𝛾 / (2pi) * B0
Larmor frequency depends only on the material (𝛾) and the main magnetic field strength (B0)
Spin Equilibrium & spin excitation
At equilibrium, the net magnetization vector M precesses about the z-axis
When M is entirely along the z-axis, only the Mz (longitudinal) component exists, and the transverse component M(xy) = 0
If M tips away from the z-axis (during excitation), a component in the x-y plane is generated
The signal is always measured in the x-y plane
The signal component is always in the transverse x-y plane
RF excitation
An RF pulse moves M away from the equilibrium state (tips away from z-axis), by using energy from the RF pulse onto the magnetic field B1
B0 is the magnet by itself, and the magnet is always on
The coil with current running through it acts as the magnet, and is surrounded by liquid Helium to keep the coil cool (decreasing resistance and energy lost to heat)
Why do we want the maximum signal? To decrease noise
B1 is applied from external antennas
B1 is along on the x-axis, so B1x is the only non-zero value for B1
Resonance comes from B1 (the additional force) being at the same frequency as B0 (𝜔0) to make sure there is maximum energy being delivery
As M returns to equilibrium, an RF signal is produced by energy releasing from the absorbed energy that was given to the material from B1
RF Pulse Shape
While oscilating at 𝜔0, the RF pulse can have three shapes, rectangular, sinc, and triangular
Rotating frame of reference
Since the spins and B1 field are both rotating, you can image from a spinning reference
The transverse plane (xy) is rotating at 𝜔, which is equal to 𝜔0, so it can be considered a rotating plane (x’-y’)
Flip angle/tip angle
The flip angle depends on the (1) shape of the RF pulse, (2) field strength B1 and (3) duration of the RF pulse 𝜏p (tau p)
For a rectangular pulse after 𝜏 seconds, the vector M has rotated at an angle α
α = 𝛾 * *B1 ** 𝜏
For any B1(t)
α (always in radians) is the integration from 0 to 𝜏 over B1(t) with respect to t (in seconds)
For the max transverse component Mxy, α is set to pi/2 radians
α = 𝛾 * *B1 ** 𝜏 = pi/2
Relaxation: The process by which M returns to steady state configuration after a B1 is done being applied
The longitudinal Mz and transverse Mxy components vary as time changes, and are descibed by Bloch Equations
Where T1 is the longitudinal relaxation time (spin minus lattice relaxation)
T2 is the transvese relaxation time (spin minus spin relaxation)
The solutions become
Mz^0 is the steady state of Mz
Mz(0+) is transverse magnetization Mx’y’ immediately after the RF pulse
We always measure the x-y plane and do everything else after
T1 is usually a LOT larger that T2
The longitudinal time is usually a lot longer than the transverse time (around 10 times more)
Mx’y’ decays at a faster rate (T2* instead of T2) because of signal desync
Spin Echo
To correct for dephasing and loss of Mx’y’ signal, another RF pulse (180 degree pulse) is applied to re-phase the spins
Spin echo makes the slower spins catch up with faster oens to sync the signal
Relaxation times T1 and T2 are dependant on materials (very high for water, really low for muscle, fat, tendons)
Contrast in MRI is the difference between signal intensities generated by different tissues (B0, which is constant, vs gamma, which is material specific)
The difference in signal intensity is used to discriminate different tissues by representing them as brightness
Pulse sequence: Carefully timed set of scanner operations used to generate images.
Gradient-echo based pulse sequences are
Spin-echo based pulse sequences
90 degree pulse
180 degree pulse
Wait for relaxation echo
Record relaxation echo
Repeat
MR signal is always recorded in the xy plane
Signal intensity S of a SE sequence is
S = K x [H] x (1 - exp(-1 x TR / T1) x exp(-1 x TE/T2)
K is the scaling factor
[H] is the proton/spin density
Twice the number of spins means you have twice as large a signal
TR is the repetition time
TE is the echo time
T1 is the longitudinal relaxation time
T2 is the transverse relaxation time
Weighting
Proton density weighting: We want a long TR and a short TE
T1 weighting: We want an intermediate TR and a short TE (becomes a constant scaling factor, less noise)
T2 weighting: We want ta long TR (minimize T1 differences) and an intermediate TE
Gradients
Adding a small gradient field (extra magnetic field) allows the same material to be differentiated if it is at different distances
Gradient affects Larmor frequency
Gradients added using more coils
Gradients are localized in x, y, and z directions
Gz is a slice selection
z-directional coils are called Maxwell pairs
Gx is frequency encoding
x-directional coils are called Golay coils
Gy is phase encoding
y-directional coils are called Golay-type coils
Gradients vary the magnetic field strength along the z-direction, which affects the frequency and phase of the spins
Fourier Transform is used to break a signal into its freuqency components (sum of sinusoids) based on their amplitude and phase
Slice thickness
Centered at z=0 , B(z) = B0 + Gz(z)
Gz is gradient strength
The magnet is used to create polarization (creates the bulk of the magnetic field)
The RF oil is used to promote excitation (sends RF energy at resonance conditions)
The gradients are used to generate spatial localization (used to form the images)
The receiver RF coils are there to target specific organs