K-space in magnetic resonance imaging explained

In magnetic resonance imaging (MRI), the k-space or reciprocal space (a mathematical space of spatial frequencies) is obtained as the 2D or 3D Fourier transform of the image measured.It was introduced in 1979 by Likes^[1] and in 1983 by Ljunggren^[2] and Twieg.^[3]

In MRI physics, complex values are sampled in k-space during an MR measurement in a premeditated scheme controlled by a pulse sequence, i.e. an accurately timed sequence of radiofrequency and gradient pulses. In practice, k-space often refers to the temporary image space, usually a matrix, in which data from digitized MR signals are stored during data acquisition. When k-space is full (at the end of the scan) the data are mathematically processed to produce a final image. Thus k-space holds raw data before reconstruction.

It can be formulated by defining wave vectors

k _{FE} and k _{PE} for "frequency encoding" (FE) and "phase encoding" (PE):k _{FE=\bar{\gamma}} G _{FEm\Delta} t k _{PE=\bar{\gamma}} n\Delta G _{PE} \tau where \Delta t is the sampling time (the reciprocal of sampling frequency), \tau is the duration of G_PE, \bar{\gamma} (gamma bar) is the gyromagnetic ratio, m is the sample number in the FE direction and n is the sample number in the PE direction (also known as partition number). Then, the 2D-Fourier Transform of this encoded signal results in a representation of the spin density distribution in two dimensions. Thus position (x,y) and spatial frequency (k _{FE} , k _{PE} ) constitute a Fourier transform pair.

Typically, k-space has the same number of rows and columns as the final image and is filled with raw data during the scan, usually one line per TR (Repetition Time).

An MR image is a complex-valued map of the spatial distribution of the transverse magnetization M_xy in the sample at a specific time point after an excitation. Conventional qualitative interpretation of Fourier Analysis asserts that low spatial frequencies (near the center of k-space) contain the signal to noise and contrast information of the image, whereas high spatial frequencies (outer peripheral regions of k-space) contain the information determining the image resolution. This is the basis for advanced scanning techniques, such as the keyhole acquisition, in which a first complete k-space is acquired, and subsequent scans are performed for acquiring just the central part of the k-space; in this way, different contrast images can be acquired without the need of running full scans.

A nice symmetry property exists in k-space if the image magnetization M_xy is prepared to be proportional simply to a contrast-weighted proton density and thus is a real quantity. In such a case, the signal at two opposite locations in k-space is:

S( - k _{FE,-k} _{PE)} = S ^{*(k} _{FE ,k} _{PE)}

where the star (

^{*} ) denotes complex conjugation.Thus k-space information is somewhat redundant then, and an image can be reconstructed using only one half of the k-space, either in the PE (Phase Encode) direction saving scan time (such a technique is known as half Fourier or half scan) or in the FE (Frequency Encode) direction, allowing for lower sampling frequencies and/or shorter echo times (such a technique is known as half echo). However, these techniques are approximate due to phase errors in the MRI data which can rarely be completely controlled (due to imperfect static field shim, effects of spatially selective excitation, signal detection coil properties, motion etc.) or nonzero phase due to just physical reasons (such as the different chemical shift of fat and water in gradient echo techniques).

MRI k-space is related to NMR time-domain^[4] in all aspects, both being used for raw data storage. The only difference between the MRI k-space and the NMR time domain is that a gradient G is present in MRI data acquisition, but is absent in NMR data acquisition. As a result of this difference, the NMR FID signal and the MRI spin-echo signal take different mathematical forms:

FID = M _{0} cos(\omega _{0t)} exp( - t/T _{2)} and Spin-Echo = M _{0} sin(\omega _{rt)/(\omega} _{rt)} where\omega _{r=\omega} _{0} + \bar{\gamma} rG Due to the presence of the gradient G, the spatial information r (not the spatial frequency information k) is encoded onto the frequency \omega , and at the same time the time-domain is renamed as k-space.

References

1. US . 4307343 . patent . Moving Gradient Zeugmatography . 1981-12-22 . 1979-08-20 . Richard S. Likes . General Electric Company.
2. Ljunggren S. Journal of Magnetic Resonance 1983; 54:338.
3. Twieg D . The k-trajectory formulation of the NMR imaging process with applications in analysis and synthesis of imaging methods. . Medical Physics . 10 . 5 . 610–21 . 1983 . 6646065 . 10.1118/1.595331. 1983MedPh..10..610T .
4. Ernst RR, Bodenhausen G and Wokaun A (1987), Principles of nuclear magnetic resonance in one and two dimensions, Oxford University Press.

Further reading

• McRobbie D., et al. MRI, From picture to proton. 2003
• Hashemi Ray, et al. MRI, The Basics 2ED. 2004.