MRI Physics 1 Image Acquisition l Chris Rorden

MRI Physics 1: Image Acquisition l Chris Rorden & Paul Morgan – – – l Magnetic Resonance Radio frequency absorption Relaxation: Radio frequency emission Gradients Increases signal/noise: antenna selection, field strength. Excellent references: – www. cis. rit. edu/htbooks/mri/ – www. easymeasure. co. uk/principlesmri. aspx 1

Units of magnetic strength l Tesla and Gauss are measures of magnetic field strength – – – Earth’s magnetic field ~0. 5 Gauss. 1 Tesla = 10, 000 Gauss. Our f. MRI system is 3 T. ~x 60, 000 earth’s field strength 2

Anatomy of an atom l l Atoms are the building blocks of our world. Composed of 3 components – – – l l Electrons: tiny negatively charged particles. Protons: heavy positively charged particles. Neutrons: heavy particles without charge. Electrons are often thought of like planets: tiny objects distantly orbiting a massive core. Neutrons and protons form the dense nucleus of an atom. 3

Atomic Nuclei l l l Atomic nuclei are composed of protons and neutrons. Protons determine the element, e. g. hydrogen has one proton, helium has two. Neutrons determines the isotope. – Helium: l l – Usually, two neutrons: 4 He. Rarely, only have one neutron: 3 He. Hydrogen l l l Usually, no neutrons: 1 H. Rarely, Deuterium (2 H) is found in ‘heavy’ water Tritium (3 H) is radioactive (spontaneous decay). 4 He 3 He 1 H 2 H 4

Nuclear Magnetic Resonance l Felix Block and Edward Purcell – – l 1946: the nuclei of some elements absorb and re-emit radio frequency energy when in magnetic field 1952: Nobel prize in physics Atoms with odd number of protons/neutrons spin in a magnetic field – – – Nuclear: properties of nuclei of atoms Magnetic: magnetic field required Resonance: interaction between magnetic field and radio frequency Bloch Purcell 5

How fast do hydrogen atoms spin? 180 Resonance Frequency (MHz) The rate of rotation is determined by the strength of the magnetic field. 40 Larmor equation f = B 0 For Hydrogen, = 42. 58 MHz/T At 1. 5 T, f = 63. 76 MHz At 3 T, f = 127. 7 MHz 1. 0 Field Strength (Tesla) 4. 0 6

Radiofrequency Pulses A radiofrequency (RF) pulse at the Larmor frequency will be absorbed. l This higher energy state tips the spin, so it is no longer aligned to the field. l An RF pulse at any other frequency will not influence the nuclei. l This is resonance. l 7

Source of MR signal l At rest, atoms align parallel to magnetic field. – l Lowest energy state. After RF absorption, net magnetization is orthogonal to magnetic field – – RF emission at larmor frequency RF emission decays with time 8

Absorption and Relaxation Our RF transmission is absorbed by atoms at Larmor frequency. l After the RF pulse, atoms will begin to realign with the magnetic field: relaxation l During this period, an RF signal is emitted. l This signal will be at the Larmor frequency. l An antenna can measure this signal. l 9

Hydrogen is the mainstay for MRI l We will focus on Hydrogen – – – l Larmor frequency varies for elements: – – l Hydrogen abundant in body (63% of atoms). Elements with even numbers of neutrons and protons have no spin, so we can not image them (4 He, 12 C). 23 Na and 31 P are relatively abundant, so can be imaged. 1 H = 42. 58 Mhz/T 13 C = 10. 7 Mhz/T 19 F = 40. 1 Mhz/T 31 P = 17. 7 Mhz/T Therefore, by sending in a RF pulse at a specific frequency we can selectively energize hydrogen. 10

l l l MRI signals are in the same range as FM radio and TV (30 -300 MHz). MRI frequency is non-ionizing radiation, unlike X-rays. Absorbed RF will cause heating. Specific absorption rate (SAR): measure of the energy absorbed by tissue. – – – Increases ~ with square of field strength. Higher SAR = more energy = more signal = more heating FDA limits SAR, and is a limiting factor for some protocols (3 W/kg averaged over 10 minutes). Non-Ionizing Radiation Heating l Ionizing Radiation Breaks Bonds Electromagnetic Spectrum Excites Electrons Excites Nuclei 11

Making a spatial image l l To create spatial images, we need a way to cause different locations in the scanner to generate different signals. To do this, we apply gradients. Gradients make the magnetic field slightly stronger at one location compared to another. Lauterbur: first MRI: 2003 Nobel Prize. Lauterbur 12

127 Mhz RF pulse Field Strength l 128 Mhz l Gradients make field stronger at one location compared to another. Larmor frequency different along this dimension. RF pulse only energizes slice where field strength matches Larmor frequency. 127 Mhz l 126 Mhz Slice Selection Gradient Larmor Freq Z Position 13

Slice Selection Gradient Gradual slice selection gradients will select thick slices, while steep gradients select thinner slices. – – Z Position Field Strength Position of gradient determines which 2 D slice is selected. Field Strength l The strength of your scanner’s gradients can limit minimum slice thickness. FDA limits speed of gradient shift (d. B/dt) and some of our protocols can elicit slight tingling sensation or brief muscle twitches. Field Strength l Z Position 14

Phase encoding gradient l l Orthogonal gradient applied between RF pulse and readout This adjusts the phase along this dimension. Analogy: Phase encoding is like timezones. Clocks in different zones will have different phases. Y Position l Field Strength 15

Frequency encoding gradient l l Apply final orthogonal gradient when we wish to acquire image. Slice will emit signal at Lamour frequency, e. g. lines at higher fields will have higher frequency signals. Aka ‘Readout gradient’. X Position l Field Strength 16

Raw MRI image: K-Space Raw MRI is 2 D Frequency Map. Apply Fourier Transform to reconstruct image. Source: Traveler’s Guide to K-space (C. A. Mistretta) 17

Reconstruction l l Medical scanners automatically reconstruct your data. You can manually reconstruct data using Tokarczuk’s free Intel Reconstruction Tool: www. mricro. com/import. html Fourier Transforms are slow: 1021 -sample data requires >2 million multiplications (2*N 2) Fast Fourier Transform: 1024 -sample data requires 20, 000 multiplications. (2(N log N)) – – Optimal when data is power of two (64, 128, 256, 512), reverts to traditional Fourier for prime numbers This is why most image matrices are a power of 2. 18

MRI scanner anatomy l A helium-cooled superconducting magnet generates the static field. – – l Always on: only quench field in emergency. niobium titanium wire. Coils allow us to – – – Make static field homogenous (shims: solenoid coils) Briefly adjust magnetic field (gradients: solenoid coils) Transmit, record RF signal (RF coils: antennas) 19

MRI scanner anatomy Antennas Magnet Gradient RF Transmit RF Receive 20

EPI l l l In conventional MRI, we collect one line of our matrix with each RF pulse. So a 64 x 64 matrix with a TR of 2 sec will be generated in 128 seconds. Problem: this is unacceptable if the object changes rapidly: – – l Images of the heart: shape changes quickly, conventional image would be blurred. Imaging brain activity: could not see changes that occur within a few seconds. Echo Planar Imaging: By rapidly applying the frequency gradient, we can collect a 2 D slab with a single RF pulse. 21

Echo Planar Imaging l EPI was devised in 1977 by Sir Mansfield. – l EPI remains the principle technique for f. MRI acquisition. – l 2003 Nobel Prize. Alternative is ‘spiral’ imaging. Mansfield EPI Multi-shot EPI does have a cost: – Image warping due to slow encoding. 22 www. fmrib. ox. ac. uk/~karla/

MRI terminology l l Orientation: typically coronal, sagittal or axial, can be inbetween these (oblique) Matrix Size: – l Voxels in each dimension Field of view: – l Axial Orientation 64 x 64 Matrix 192 mm FOV 3 x 3 mm Resolution Spatial extent of each dimension. Resolution: – FOV/Matrix size. Sagittal Orientation 256 x 256 Matrix 256 mm FOV 1 x 1 mm Resolution 23

Volumes l l 3 D volumes are composed of stacks of 2 D slices, like a loaf of bread. Each slice has a thickness. – Thicker slices have more hydrogen, so more signal l l – l Volume of 1 x 1 x 1 mm voxel is 1 mm 3 Volume of 1 x 1 x 2 mm voxel is 2 mm 3 Thinner slices provide higher resolution. Optional: gap between slices. – – Reduces RF interference Allows fewer slices to cover whole brain. 1 mm Gap 3 mm 2 mm Thick 24

Creating 3 D volumes from 2 D EPI l Slices can be either collected sequentially (e. g. ascending axial slices) or in interleaved order (e. g. collect all odd slices, then all even slices). Sequential Interleaved … … Time 25

Sequential vs Interleaved Volumes l l l Interleaving reduces interference between slices. Useful if thin gap between slices. Unfortunately, any head movements will cause worse ‘spin history’ effects for interleaved slices than sequential. (TR will be shorter for regions that were previously in another slice). You must know slice order if you want to ‘slice time correct’ data (temporal processing lecture). Slice 4 Slice 3 Slice 2 Slice 1 lice 44 SSlice 3 S Slice 22 Slice 11 26

Signal to Noise Signal To noise is given by the formula Signal = V N l Where V is the volume and N is the number of samples averaged (referred to as ‘Nex’, as in ‘number of excitations’). l For example, to get the same SNR as a single 3 x 3 x 3 mm scan (27 mm 3) we would need to collect 12 2 x 2 x 2 mm (8 mm 3) scans or 768 1 x 1 x 1 mm (1 mm 3) scans. l 27

www. fmrib. ox. ac. uk/~karla/ Signal to Noise: Antennas l l The MRI antenna is called a coil. We use different coils for different body parts. For brains, the most common antenna is the head coil, which is a volume coil: it shows the whole brain. We can also use a surface coil: it gives great signal for a small field of view. Head coil Surface coil Volume coil Surface coil 28

Parallel Imaging (SENSE, i. Pat) l l l Parallel imaging uses multiple surface coils to generate a volume image. Dramatically reduces spatial distortion and increases signal. Optionally, you can acquire images more rapidly by only collecting a portion of kspace. – – – SENSE R=2 collects half of the lines. SENSE R=3 collects one third Reduces spatial distortion and increases speed of acquisition. Some loss in signal. 8 -channel array 29

Parallel Imaging (SENSE, i. Pat) Increasing SENSE reduction factor decreases acquisition time and spatial distortion, but high values lead to reduced signal. Effects of SENSE factor (R) on EPI l R=1 R=2 R=3 30

Signal and Field Strength Outside magnetic field • Spins randomly oriented In magnetic field: • Spins tend to align parallel or anti-parallel to magnetic field. • At room temperature, ~4 parts per million more protons per Tesla align with versus against field. • As field strength increases, there is a bigger energy difference between parallel and anti-parallel alignment (faster rotation = more energy). • A larger proportion will align parallel to field. • More energy will be released as nuclei align. • Therefore, MR signal increases with square of field strength. 31

Signal and Field Strength l l l Most clinical MRI: 1. 5 T Our f. MRI systems: 3. 0 T Maximum for Nb. Ti MRI ~11. 7 T World record for superconductive MRI ~22 T (Nb 3 Sn) Field strength influences: – Faster Larmor frequency : l Bigger energy difference between parallel and anti-parallel alignment Larger ratio of nuclei aligned = more signal • More signal as nuclei realign. • l Reduced TR and TE: less time to take images (next week). 32

Signal and Field Strength l l In theory: – Signal increases with square of field strength – Noise increases linearly with field strength – A 3 T scanner should have twice SNR of 1. 5 T scanner; 7 T should have ~4. 7 times SNR of 1. 5 T. Unfortunately, physiological artifacts also increase, so advantage is less in practice. Benefits: speed, resolution Costs: Artifacts, Money, SAR, wavelength effects, auditory noise 33

Magnetic attraction l Force in magnetic field (1. 8 T, unit dynes) – – l l l Water: -22 Copper: -2. 6 Copper Chloride: +280 Iron: +400, 000 Diamagnetic material repelled (e. g. H 2 O). Paramagnetic material attracted (e. g. Cu. Cl, Gd). Ferromagnetic material strongly attracted (Fe). – – Even without magnetic field, magnetic moments aligned Dangerous near MRI scanner 34