“Physics of magnetic resonance imaging”的意思、由来-开放百科全书

The physics of magnetic resonance imaging (MRI) concerns fundamental physical considerations of MRI techniques and technological aspects of MRI devices. MRI is a medical imaging technique mostly used in radiology and nuclear medicine in order to investigate the anatomy and physiology of the body, and to detect pathologies including tumors, inflammation, neurological conditions such as stroke, disorders of muscles and joints, and abnormalities in the heart and blood vessels among others. Contrast agents may be injected intravenously or into a joint to enhance the image and facilitate diagnosis. Unlike CT and X-ray, MRI uses no ionizing radiation and is therefore a safe procedure suitable for diagnosis in children and repeated runs. Patients with specific non-ferromagnetic metal implants, cochlear implants, and cardiac pacemakers nowadays may also have an MRI in spite of effects of the strong magnetic fields. This does not apply on older devices, details for medical professionals are provided by the device's manufacturer.

Certain atomic nuclei are able to absorb and emit radio frequency energy when placed in an external magnetic field. In clinical and research MRI, hydrogen atoms are most often used to generate a detectable radio-frequency signal that is received by antennas in close proximity to the anatomy being examined. Hydrogen atoms are naturally abundant in people and other biological organisms, particularly in water and fat. For this reason, most MRI scans essentially map the location of water and fat in the body. Pulses of radio waves excite the nuclear spin energy transition, and magnetic field gradients localize the signal in space. By varying the parameters of the pulse sequence, different contrasts may be generated between tissues based on the relaxation properties of the hydrogen atoms therein.

When inside the magnetic field (B₀) of the scanner, the magnetic moments of the protons align with the direction of the field. They may align in two configurations parallel (in the direction of B₀), or anti-parallel opposing B₀. While each individual proton can only have one of two alignments, the collection of protons appear to behave as though they can have any alignment. Most protons align parallel to B₀ as this is a lower energy state. A radio frequency pulse is then applied, which can excite protons from parallel to anti-parallel alignment, only the latter are relevant to the rest of the discussion. In response to the force bringing them back to their equilibrium orientation, the protons undergo a rotating motion (precession), much like a spin wheel under the effect of gravity. The protons will return to the low energy state by the process of spin-lattice relaxation. This appears as a magnetic flux, which yields a changing voltage in receiver coils to give the signal. The frequency at which a proton or group of protons in a voxel resonates depends on the strength of the local magnetic field around the proton or group of protons, a stronger field corresponds to a larger energy difference and higher frequency photons. By applying additional magnetic fields (gradients) that vary linearly over space, specific slices to be imaged can be selected, and an image is obtained by taking the 2-D Fourier transform of the spatial frequencies of the signal (k-space). Due to the magnetic Lorentz force from B₀ on the current flowing in the gradient coils, the gradient coils will try to move producing loud knocking sounds, for which patients require hearing protection.

Nuclear magnetism

Subatomic particles have the quantum mechanical property of spin.^[1] Certain nuclei such as ¹H (protons), ²H, ³He, ²³Na or ³¹P, have a non–zero spin and therefore a magnetic moment. In the case of the so-called spin//Precession">precess around an axis along the direction of the field. Protons align in two energy eigenstates (the Zeeman effect): one low-energy and one high-energy, which are separated by a very small splitting energy.

Resonance and relaxation

Quantum mechanics is required to accurately model the behaviour of a single proton, however classical mechanics can be used to describe the behaviour of the ensemble of protons adequately. As with other spin

particles, whenever the spin of a single proton is measured it can only have one of two results commonly called parallel and anti-parallel. When we discuss the state of a proton or protons we are referring to the wavefunction of that proton which is a linear combination of the parallel and anti-parallel states.

In the presence of the magnetic field, B₀, the protons will appear to precess at the Larmor frequency determined by the particle's gyro-magnetic ratio and the strength of the field. The static fields used most commonly in MRI cause precession which corresponds to a radiowave photon.

The net longitudinal magnetization in thermodynamical equilibrium is due to a tiny excess of protons in the lower energy state. This gives a net polarization that is parallel to the external field. Application of a radio frequency (RF) pulse can tip this net polarization vector sideways (with, i.e., a so-called 90° pulse), or even reverse it (with a so-called 180° pulse). The protons will come into phase with the RF pulse and therefore each other.

The recovery of longitudinal magnetization is called longitudinal or T₁ relaxation and occurs exponentially with a time constant T₁. The loss of phase coherence in the transverse plane is called transverse or T₂ relaxation. T₁ is thus associated with the enthalpy of the spin system, or the number of nuclei with parallel versus anti-parallel spin. T₂ on the other hand is associated with the entropy of the system, or the number of nuclei in phase.

When the radio frequency pulse is turned off, the transverse vector component produces an oscillating magnetic field which induces a small current in the receiver coil. This signal is called the free induction decay (FID). In an idealized nuclear magnetic resonance experiment, the FID decays approximately exponentially with a time constant T₂. However, in practical MRI there are small differences in the static magnetic field at different spatial locations ("inhomogeneities") that cause the Larmor frequency to vary across the body. This creates destructive interference, which shortens the FID. The time constant for the observed decay of the FID is called the T{{su|p=*|b=2}} relaxation time, and is always shorter than T₂. At the same time, the longitudinal magnetization starts to recover exponentially with a time constant T₁ which is much larger than T₂ (see below).

In MRI, the static magnetic field is caused to vary across the body (by using a field gradient), so that different spatial locations become associated with different precession frequencies. Usually these field gradients are pulsed, and it is the almost infinite variety of RF and gradient pulse sequences that gives MRI its versatility. Application of field gradient destroys the FID signal, but this can be recovered and measured by a refocusing gradient (to create a so-called "gradient echo"), or by a radio frequency pulse (to create a so-called "spin-echo"). The whole process can be repeated when some T₁-relaxation has occurred and the thermal equilibrium of the spins has been more or less restored. The repetition time (TR) is the time between two successive excitations of the same slice.^{[weighted images and may be used for liver imaging, as normal liver tissue retains the agent, but abnormal areas (e.g., scars, tumors) do not. They can also be taken orally, to improve visualization of the gastrointestinal tract, and to prevent water in the gastrointestinal tract from obscuring other organs (e.g., the pancreas). Diamagnetic agents such as barium sulfate have also been studied for potential use in the gastrointestinal tract, but are less frequently used.}

k-space

In 1983, Ljunggren^[9] and Twieg^[10] independently introduced the k-space formalism, a technique that proved invaluable in unifying different MR imaging techniques. They showed that the demodulated MR signal S(t) generated by freely precessing nuclear spins in the presence of a linear magnetic field gradient G equals the Fourier transform of the effective spin density. Mathematically:

In other words, as time progresses the signal traces out a trajectory in k-space with the velocity vector of the trajectory proportional to the vector of the applied magnetic field gradient.

By the term effective spin density we mean the true spin density

corrected for the effects of T₁ preparation, T₂ decay, dephasing due to field inhomogeneity, flow, diffusion, etc. and any other phenomena that affect that amount of transverse magnetization available to induce signal in the RF probe or its phase with respect to the receiving coil' s electromagnetic field.

From the basic k-space formula, it follows immediately that we reconstruct an image

simply by taking the inverse Fourier transform of the sampled data, viz.

Using the k-space formalism, a number of seemingly complex ideas became simple. For example, it becomes very easy to understand the role of phase encoding (the so-called spin-warp method). In a standard spin echo or gradient echo scan, where the readout (or view) gradient is constant (e.g., G), a single line of k-space is scanned per RF excitation. When the phase encoding gradient is zero, the line scanned is the k_x axis. When a non-zero phase-encoding pulse is added in between the RF excitation and the commencement of the readout gradient, this line moves up or down in k-space, i.e., we scan the line k_y = constant.

The k-space formalism also makes it very easy to compare different scanning techniques. In single-shot EPI, all of k-space is scanned in a single shot, following either a sinusoidal or zig-zag trajectory. Since alternating lines of k-space are scanned in opposite directions, this must be taken into account in the reconstruction. Multi-shot EPI and fast spin echo techniques acquire only part of k-space per excitation. In each shot, a different interleaved segment is acquired, and the shots are repeated until k-space is sufficiently well-covered. Since the data at the center of k-space represent lower spatial frequencies than the data at the edges of k-space, the T_E value for the center of k-space determines the image's T₂ contrast.

The importance of the center of k-space in determining image contrast can be exploited in more advanced imaging techniques. One such technique is spiral acquisition—a rotating magnetic field gradient is applied, causing the trajectory in k-space to spiral out from the center to the edge. Due to T₂ and T{{su|p=*|b=2}} decay the signal is greatest at the start of the acquisition, hence acquiring the center of k-space first improves

contrast to noise ratio (CNR) when compared to conventional zig-zag acquisitions, especially in the presence of rapid movement.

Since

and

are conjugate variables (with respect to the Fourier transform) we can use the Nyquist theorem to show that the step in k-space determines the field of view of the image (maximum frequency that is correctly sampled) and the maximum value of k sampled determines the resolution; i.e.,

Example of a pulse sequence

In the timing diagram, the horizontal axis represents time. The vertical axis represents: (top row) amplitude of radio frequency pulses; (middle rows) amplitudes of the three orthogonal magnetic field gradient pulses; and (bottom row) receiver analog-to-digital converter (ADC). Radio frequencies are transmitted at the Larmor frequency of the nuclide to be imaged. For example, for ¹H in a magnetic field of 1 T, a frequency of 42.5781 MHz would be employed. The three field gradients are labeled G_X (typically corresponding to a patient's left-to-right direction and colored red in diagram), G_Y (typically corresponding to a patient's front-to-back direction and colored green in diagram), and G_Z (typically corresponding to a patient's head-to-toe direction and colored blue in diagram). Where negative-going gradient pulses are shown, they represent reversal of the gradient direction, i.e., right-to-left, back-to-front or toe-to-head. For human scanning, gradient strengths of 1–100 mT/m are employed: Higher gradient strengths permit better resolution and faster imaging. The pulse sequence shown here would produce a transverse (axial) image.

The first part of the pulse sequence, SS, achieves "slice selection". A shaped pulse (shown here with a sinc modulation) causes a 90° nutation of longitudinal nuclear magnetization within a slab, or slice, creating transverse magnetization. The second part of the pulse sequence, PE, imparts a phase shift upon the slice-selected nuclear magnetization, varying with its location in the Y direction. The third part of the pulse sequence, another slice selection (of the same slice) uses another shaped pulse to cause a 180° rotation of transverse nuclear magnetization within the slice. This transverse magnetisation refocuses to form a spin echo at a time T_E. During the spin echo, a frequency-encoding (FE) or readout gradient is applied, making the resonant frequency of the nuclear magnetization vary with its location in the X direction. The signal is sampled n_FE times by the ADC during this period, as represented by the vertical lines. Typically n_FE of between 128 and 512 samples are taken.

The longitudinal magnetisation is then allowed to recover somewhat and after a time T_R the whole sequence is repeated n_PE times, but with the phase-encoding gradient incremented (indicated by the horizontal hatching in the green gradient block). Typically n_PE of between 128 and 512 repetitions are made.

The negative-going lobes in G_X and G_Z are imposed to ensure that, at time T_E (the spin echo maximum), phase only encodes spatial location in the Y direction.

Typically T_E is between 5 ms and 100 ms, while T_R is between 100 ms and 2000 ms.

After the two-dimensional matrix (typical dimension between 128 × 128 and 512 × 512) has been acquired, producing the so-called k-space data, a two-dimensional inverse Fourier transform is performed to provide the familiar MR image. Either the magnitude or phase of the Fourier transform can be taken, the former being far more common.

Overview of main sequences

MRI scanner

Construction and operation

The major components of an MRI scanner are: the main magnet, which polarizes the sample, the shim coils for correcting inhomogeneities in the main magnetic field, the gradient system which is used to localize the MR signal and the RF system, which excites the sample and detects the resulting NMR signal. The whole system is controlled by one or more computers.

Magnet

The magnet is the largest and most expensive component of the scanner, and the remainder of the scanner is built around it. The strength of the magnet is measured in teslas (T). Clinical magnets generally have a field strength in the range 0.1–3.0 T, with research systems available up to 9.4 T for human use and 21 T for animal systems.^[11]

In the United States, field strengths up to 4 T have been approved by the FDA for clinical use.^[12]

Just as important as the strength of the main magnet is its precision. The straightness of the magnetic lines within the center (or, as it is technically known, the iso-center) of the magnet needs to be near-perfect. This is known as homogeneity. Fluctuations (inhomogeneities in the field strength) within the scan region should be less than three parts per million (3 ppm). Three types of magnets have been used:

Most superconducting magnets have their coils of superconductive wire immersed in liquid helium, inside a vessel called a cryostat. Despite thermal insulation, sometimes including a second cryostat containing liquid nitrogen, ambient heat causes the helium to slowly boil off. Such magnets, therefore, require regular topping-up with liquid helium. Generally a cryocooler, also known as a coldhead, is used to recondense some helium vapor back into the liquid helium bath. Several manufacturers now offer 'cryogenless' scanners, where instead of being immersed in liquid helium the magnet wire is cooled directly by a cryocooler.^[13]

Magnets are available in a variety of shapes. However, permanent magnets are most frequently 'C' shaped, and superconducting magnets most frequently cylindrical. However, C-shaped superconducting magnets and box-shaped permanent magnets have also been used.

Magnetic field strength is an important factor in determining image quality. Higher magnetic fields increase signal-to-noise ratio, permitting higher resolution or faster scanning. However, higher field strengths require more costly magnets with higher maintenance costs, and have increased safety concerns. A field strength of 1.0–1.5 T is a good compromise between cost and performance for general medical use. However, for certain specialist uses (e.g., brain imaging) higher field strengths are desirable, with some hospitals now using 3.0 T scanners.

Shims

When the MR scanner is placed in the hospital or clinic, its main magnetic field is far from being homogeneous enough to be used for scanning. That is why before doing fine tuning of the field using a sample, the magnetic field of the magnet must be measured and shimmed.

After a sample is placed into the scanner, the main magnetic field is distorted by susceptibility boundaries within that sample, causing signal dropout (regions showing no signal) and spatial distortions in acquired images. For humans or animals the effect is particularly pronounced at air-tissue boundaries such as the sinuses (due to paramagnetic oxygen in air) making, for example, the frontal lobes of the brain difficult to image. To restore field homogeneity a set of shim coils is included in the scanner. These are resistive coils, usually at room temperature, capable of producing field corrections distributed as several orders of spherical harmonics.^[14]

After placing the sample in the scanner, the B₀ field is 'shimmed' by adjusting currents in the shim coils. Field homogeneity is measured by examining an FID signal in the absence of field gradients. The FID from a poorly shimmed sample will show a complex decay envelope, often with many humps. Shim currents are then adjusted to produce a large amplitude exponentially decaying FID, indicating a homogeneous B₀ field. The process is usually automated.^[15]

Gradients

Gradient coils are used to spatially encode the positions of protons by varying the magnetic field linearly across the imaging volume. The Larmor frequency will then vary as a function of position in the x, y and z-axes.

Gradient coils are usually resistive electromagnets powered by sophisticated amplifiers which permit rapid and precise adjustments to their field strength and direction. Typical gradient systems are capable of producing gradients from 20–100 mT/m (i.e., in a 1.5 T magnet, when a maximal z-axis gradient is applied, the field strength may be 1.45 T at one end of a 1 m long bore and 1.55 T at the other^[16]). It is the magnetic gradients that determine the plane of imaging—because the orthogonal gradients can be combined freely, any plane can be selected for imaging.

Scan speed is dependent on performance of the gradient system. Stronger gradients allow for faster imaging, or for higher resolution; similarly, gradient systems capable of faster switching can also permit faster scanning. However, gradient performance is limited by safety concerns over nerve stimulation.

Some important characteristics of gradient amplifiers and gradient coils are slew rate and gradient strength. As mentioned earlier, a gradient coil will create an additional, linearly varying magnetic field that adds or subtracts from the main magnetic field. This additional magnetic field will have components in all 3 directions, viz. x, y and z; however, only the component along the magnetic field (usually called the z-axis, hence denoted G_z) is useful for imaging. Along any given axis, the gradient will add to the magnetic field on one side of the zero position and subtract from it on the other side. Since the additional field is a gradient, it has units of gauss per centimeter or millitesla per meter (mT/m). High performance gradient coils used in MRI are typically capable of producing a gradient magnetic field of approximate 30 mT/m or higher for a 1.5 T MRI. The slew rate of a gradient system is a measure of how quickly the gradients can be ramped on or off. Typical higher performance gradients have a slew rate of up to 100–200 T·m⁻¹·s⁻¹. The slew rate depends both on the gradient coil (it takes more time to ramp up or down a large coil than a small coil) and on the performance of the gradient amplifier (it takes a lot of voltage to overcome the inductance of the coil) and has significant influence on image quality.

Radio frequency system

The radio frequency (RF) transmission system consists of an RF synthesizer, power amplifier and transmitting coil. That coil is usually built into the body of the scanner. The power of the transmitter is variable, but high-end whole-body scanners may have a peak output power of up to 35 kW,^[17] and be capable of sustaining average power of 1 kW. Although these electromagnetic fields are in the RF range of tens of megahertz (often in the shortwave radio portion of the electromagnetic spectrum) at powers usually exceeding the highest powers used by amateur radio, there is very little RF interference produced by the MRI machine. The reason for this, is that the MRI is not a radio transmitter. The RF frequency electromagnetic field produced in the "transmitting coil" is a magnetic near-field with very little associated changing electric field component (such as all conventional radio wave transmissions have). Thus, the high-powered electromagnetic field produced in the MRI transmitter coil does not produce much electromagnetic radiation at its RF frequency, and the power is confined to the coil space and not radiated as "radio waves." Thus, the transmitting coil is a good EM field transmitter at radio frequency, but a poor EM radiation transmitter at radio frequency.

The receiver consists of the coil, pre-amplifier and signal processing system. The RF electromagnetic radiation produced by nuclear relaxation inside the subject is true EM radiation (radio waves), and these leave the subject as RF radiation, but they are of such low power as to also not cause appreciable RF interference that can be picked up by nearby radio tuners (in addition, MRI scanners are generally situated in metal mesh lined rooms which act as Faraday cages.)

While it is possible to scan using the integrated coil for RF transmission and MR signal reception, if a small region is being imaged, then better image quality (i.e., higher signal-to-noise ratio) is obtained by using a close-fitting smaller coil. A variety of coils are available which fit closely around parts of the body such as the head, knee, wrist, breast, or internally, e.g., the rectum.

A recent development in MRI technology has been the development of sophisticated multi-element phased array^[18] coils which are capable of acquiring multiple channels of data in parallel. This 'parallel imaging' technique uses unique acquisition schemes that allow for accelerated imaging, by replacing some of the spatial coding originating from the magnetic gradients with the spatial sensitivity of the different coil elements. However, the increased acceleration also reduces the signal-to-noise ratio and can create residual artifacts in the image reconstruction. Two frequently used parallel acquisition and reconstruction schemes are known as SENSE^[19] and GRAPPA.^[20] A detailed review of parallel imaging techniques can be found here:^[21]

References

1. ^{{cite book | vauthors = Callaghan P |date = 1994 |title = Principles of Nuclear Magnetic Resonance Microscopy |pages = |publisher = Oxford University Press |isbn = 978-0-19-853997-1}}
2. ^[https://books.google.se/books?id=T2YcBQAAQBAJ&pg=PT26 Page 26] in: {{cite book|title=How does MRI work?: An Introduction to the Physics and Function of Magnetic Resonance Imaging| first1 = Dominik | last1 = Weishaupt | first2 = Victor D. | last2 = Koechli | first3 = Borut | last3 = Marincek | name-list-format = vanc |publisher=Springer Science & Business Media |year=2013 |isbn=978-3-662-07805-1}}
3. ^{{cite journal | vauthors = Poustchi-Amin M, Mirowitz SA, Brown JJ, McKinstry RC, Li T | title = Principles and applications of echo-planar imaging: a review for the general radiologist | journal = Radiographics | volume = 21 | issue = 3 | pages = 767–79 | date = 2000 | pmid = 11353123 | doi = 10.1148/radiographics.21.3.g01ma23767 }}
4. ^{{cite journal | vauthors = Feinberg DA, Moeller S, Smith SM, Auerbach E, Ramanna S, Gunther M, Glasser MF, Miller KL, Ugurbil K, Yacoub E | title = Multiplexed echo planar imaging for sub-second whole brain FMRI and fast diffusion imaging | journal = PLOS One | volume = 5 | issue = 12 | pages = e15710 | date = December 2010 | pmid = 21187930 | pmc = 3004955 | doi = 10.1371/journal.pone.0015710 | bibcode = 2010PLoSO...515710F }}
5. ^{{cite journal | vauthors = Weinmann HJ, Brasch RC, Press WR, Wesbey GE | title = Characteristics of gadolinium-DTPA complex: a potential NMR contrast agent | journal = AJR. American Journal of Roentgenology | volume = 142 | issue = 3 | pages = 619–24 | date = March 1984 | pmid = 6607655 | doi = 10.2214/ajr.142.3.619 | url = http://www.ajronline.org/cgi/pmidlookup?view=long&pmid=6607655 }}
6. ^{{cite journal | vauthors = Laniado M, Weinmann HJ, Schörner W, Felix R, Speck U | title = First use of GdDTPA/dimeglumine in man | journal = Physiological Chemistry and Physics and Medical NMR | volume = 16 | issue = 2 | pages = 157–65 | date = 1984 | pmid = 6505042 | doi = }}
7. ^{{cite journal | vauthors = Widder DJ, Greif WL, Widder KJ, Edelman RR, Brady TJ | title = Magnetite albumin microspheres: a new MR contrast material | journal = AJR. American Journal of Roentgenology | volume = 148 | issue = 2 | pages = 399–404 | date = February 1987 | pmid = 3492120 | doi = 10.2214/ajr.148.2.399 | url = http://www.ajronline.org/cgi/pmidlookup?view=long&pmid=3492120 }}
8. ^{{cite journal | vauthors = Weissleder R, Elizondo G, Wittenberg J, Rabito CA, Bengele HH, Josephson L | title = Ultrasmall superparamagnetic iron oxide: characterization of a new class of contrast agents for MR imaging | journal = Radiology | volume = 175 | issue = 2 | pages = 489–93 | date = May 1990 | pmid = 2326474 | doi = 10.1148/radiology.175.2.2326474 | url = http://radiology.rsnajnls.org/cgi/pmidlookup?view=long&pmid=2326474 }}
9. ^{{cite journal | vauthors = Ljunggren S | date = 1983 |title=A simple graphical representation of Fourier-based imaging methods |journal=Journal of Magnetic Resonance |volume=54 |issue=2 |pages=338–343 |doi=10.1016/0022-2364(83)90060-4|bibcode=1983JMagR..54..338L}}
10. ^{{cite journal | vauthors = Twieg DB | title = The k-trajectory formulation of the NMR imaging process with applications in analysis and synthesis of imaging methods | journal = Medical Physics | volume = 10 | issue = 5 | pages = 610–21 | date = 1983 | pmid = 6646065 | doi = 10.1118/1.595331 | bibcode = 1983MedPh..10..610T }}
11. ^{{cite web |title=In vivo MR Imaging at 21.1 T |url=http://www.magnet.fsu.edu/mediacenter/factsheets/documents/invivomri-maglabreports-2007-v14-i4.pdf}}
12. ^{{cite web |last=Duggan-Jahns |first=Terry |title=The Evolution of Magnetic Resonance Imaging: 3T MRI in Clinical Applications |url=http://www.eradimaging.com/site/article.cfm?ID=426 |work=eRADIMAGING.com |publisher=eRADIMAGING.com |access-date=24 June 2013}}
13. ^{{Cite book | vauthors = Obasih KM, Mruzek | veditors = Timmerhaus KD | title = Proceedings of the 1995 cryogenic engineering conference | chapter = Thermal design and analysis of a cryogenless superconducting magnet for interventional MRI therapy | publisher = Plenum Press | date = 1996 | location = New York | pages = 305–312 | isbn = 978-0-306-45300-7 }}
14. ^{{cite book | vauthors = Chen CN, Hoult DH |date = 1989 |title = Biomedical Magnetic Resonance Technology |publisher = Taylor & Francis |series = Medical Sciences |isbn = 978-0-85274-118-4 }}
15. ^{{cite journal | vauthors = Gruetter R | title = Automatic, localized in vivo adjustment of all first- and second-order shim coils | journal = Magnetic Resonance in Medicine | volume = 29 | issue = 6 | pages = 804–11 | date = June 1993 | pmid = 8350724 | doi = 10.1002/mrm.1910290613 }}
16. ^This unrealistically assumes that the gradient is linear out to the end of the magnet bore. While this assumption is fine for pedagogical purposes, in most commercial MRI systems the gradient droops significantly after a much smaller distance; indeed, the decrease in the gradient field is the main delimiter of the useful field of view of a modern commercial MRI system.
17. ^{{cite book | vauthors = Oppelt A | date = 2006 | title = Imaging Systems for Medical Diagnostics: Fundamentals, Technical Solutions and Applications for Systems Applying Ionizing Radiation, Nuclear Magnetic Resonance and Ultrasound | page = 566 | publisher = Wiley-VCH | isbn = 978-3-89578-226-8 }}
18. ^{{cite journal | vauthors = Roemer PB, Edelstein WA, Hayes CE, Souza SP, Mueller OM | title = The NMR phased array | journal = Magnetic Resonance in Medicine | volume = 16 | issue = 2 | pages = 192–225 | date = November 1990 | pmid = 2266841 | doi = 10.1002/mrm.1910160203 }}
19. ^{{cite journal | vauthors = Pruessmann KP, Weiger M, Scheidegger MB, Boesiger P | title = SENSE: sensitivity encoding for fast MRI | journal = Magnetic Resonance in Medicine | volume = 42 | issue = 5 | pages = 952–62 | date = November 1999 | pmid = 10542355 | doi = 10.1002/(SICI)1522-2594(199911)42:5<952::AID-MRM16>3.0.CO;2-S | citeseerx = 10.1.1.139.3032 }}
20. ^{{cite journal | vauthors = Griswold MA, Jakob PM, Heidemann RM, Nittka M, Jellus V, Wang J, Kiefer B, Haase A | title = Generalized autocalibrating partially parallel acquisitions (GRAPPA) | journal = Magnetic Resonance in Medicine | volume = 47 | issue = 6 | pages = 1202–10 | date = June 2002 | pmid = 12111967 | doi = 10.1002/mrm.10171 | citeseerx = 10.1.1.462.3159 }}
21. ^{{cite journal | vauthors = Blaimer M, Breuer F, Mueller M, Heidemann RM, Griswold MA, Jakob PM |date=2004 |title=SMASH, SENSE, PILS, GRAPPA: How to Choose the Optimal Method | url = http://cfmriweb.ucsd.edu/ttliu/be280a_05/blaimer05.pdf |journal= Topics in Magnetic Resonance Imaging |volume=15 |issue=4 |pages=223–236 |doi=10.1097/01.rmr.0000136558.09801.dd }}