knowt logo

10. MR Physics & Imaging

Introduction

  • 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

Physics

  • 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

Bulk Magnetization Vector

  • 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)

Generating Data

  • 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)

Acquisition & Contrast

  • 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 Intensity

  • 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

Overview

  • 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

GV

10. MR Physics & Imaging

Introduction

  • 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

Physics

  • 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

Bulk Magnetization Vector

  • 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)

Generating Data

  • 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)

Acquisition & Contrast

  • 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 Intensity

  • 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

Overview

  • 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