# User:GyroMagician/Imaging

### Imaging

Signal Localization using gradients, initially PR but Fourier better (why?)

Creating an image of a sample requires a process to localize the NMR signal. This can be done by superimposing a small, spatially varying magnetic field onto the main field, making the resonant frequency a function of position. These applied fields are called field gradients.

Consider a magnetic field ${\displaystyle B_{grad}}$ aligned to the z-axis, the magnitude of which varies linearly in the x direction

${\displaystyle B_{grad}=x\cdot G_{x}{\hat {z}}\,}$

where ${\displaystyle G_{x}}$ is a constant. Applying the gradient field to the sample causes the Larmor frequency vary linearly as a function of position along the x-axis

${\displaystyle \omega (x)=\gamma (B_{0}+x.G_{x})\,}$

An MRI system typically has three independent orthogonal gradient coils, each designed to produce a field gradient ${\displaystyle G_{x}}$, ${\displaystyle G_{y}}$ and ${\displaystyle G_{z}}$ varying linearly along its respective axis, producing a total gradient field

${\displaystyle G=G_{x}{\hat {z}}+G_{y}{\hat {z}}+G_{z}{\hat {z}}\,}$

The magnetic field at point r is then

${\displaystyle B(r)=(B_{0}+G\cdot r){\hat {z}}\,}$

and the resonant frequency is

${\displaystyle \omega (r)=\omega _{0}+\gamma G\cdot r}$

#### Slice Selection

Slice selection is a technique used to excite a slab of spins in a sample and forms the first stage of most imaging sequences. Applying a linear field gradient across the sample makes the spin resonant frequency a function of distance along the gradient axis. If the slice select gradient is applied along the z-axis,

${\displaystyle \omega (z)=\gamma (B_{0}+z.G_{sl})\,}$

If an RF pulse is now applied to the sample, only those spins resonant at or close to the frequency of the RF pulse will be affected (such that ${\displaystyle B_{0} [16]), so that only a slice of the sample is excited (fig. 2.4). The position of the excited slice along ${\displaystyle G_{sl}}$ can be changed by changing the frequency of the RF pulse. The thickness of the slice is determined by the gradient amplitude.

#### Frequency Encoding

Frequency encoding uses the same principle as slice selection, but the field gradient is applied during image readout rather than excitation and so only influences the precession frequency of spins that have already been excited. Applying a linear field gradient ${\displaystyle G_{fr}}$ across the sample makes the precession frequency a function of position along the gradient axis

${\displaystyle \omega (x)=-\gamma (B_{0}+G_{\textit {fr}}x)=\omega _{0}+\gamma G_{\textit {fr}}x\,}$

The observed signal is then

${\displaystyle s(t)=\int _{sample}\rho (x)\exp(-i(\omega _{0}+\gamma G_{\textit {fr}}x)t)dx\,}$

where ${\displaystyle \rho (x)}$ is the projection of the spin density function along the x-axis

${\displaystyle \rho (x)=\int \int \rho (\mathbf {r} )dydz\,}$

The signal is then demodulated to remove the RF carrier frequency ${\displaystyle \omega _{0}}$, leaving the envelope

${\displaystyle s(t)=\int _{sample}\rho (x)\exp(-i(\gamma G_{\textit {fr}}t)x)dx\,}$

which is the Fourier transform of ${\displaystyle \rho (x)}$. Fourier transforming the acquired signal produces a projection of the sample magnetisation along the gradient axis (fig.1d_localisation). If the frequency encoding gradient is orthogonal to the slice select gradient, the two techniques can be combined to image the spin density of a column through the sample (fig.select_pencil).

#### Phase Encoding

Frequency encoding acquires a series of time points during a single echo (fig.phase_encode_a). The same result may be achieved by exciting the sample a number of times and applying a constant gradient for a different duration ${\displaystyle T_{\textit {ph}}}$ each time before sampling a single data point (fig.phase_encode_b)~\cite{kumar75}. If the ${\displaystyle y}$ gradient is used, the acquired signal is

${\displaystyle s(T_{\textit {ph}})=\int _{sample}\rho (y)\exp(-i\gamma G_{\textit {ph}}T_{\textit {ph}}y)dy\,}$

Repeating the experiment at many different values of ${\displaystyle T_{\textit {ph}}}$ samples the complete signal. The Fourier transform of the acquired signal is then the projection of signal intensity along the gradient axis as before.

While requiring more acquisitions than frequency encoding (one per sample point), phase encoding can be separately applied to all three orthogonal axes, localising the signal to a point and making 3D imaging possible. Phase encoding can also be used in combination with slice selection and frequency encoding, forming the basis for most common imaging sequences.

Varying ${\displaystyle T_{\textit {ph}}}$ allows the complete echo signal to be acquired. However, it also means that the time between the excitation pulse and the sampled point varies. During this time the echo signal undergoes ${\displaystyle T_{1}}$ and ${\displaystyle T_{2}^{*}}$ decay.

Returning to eqn.~\ref{eqn:phase_encode} it can be seen that the location of the sampled point is determined by the product $G_\textit{ph} T_\textit{ph}$. However, as the gradient is switched off by the time the signal is sampled, it is the sum effect rather than the exact waveform that is important. Eqn.~\ref{eqn:phase_encode} can be written more generally as

${\displaystyle s(T_{\textit {ph}})=\int _{sample}\rho (y)\exp(-i\phi (T_{\textit {ph}}))dy\,}$

where

${\displaystyle \phi (T_{\textit {ph}})=\gamma \int _{0}^{T_{\textit {ph}}}G_{\textit {ph}}\,d\tau }$

For a constant gradient amplitude ${\displaystyle \phi (T_{\textit {ph}})}$ reduces to ${\displaystyle G_{\textit {ph}}T_{\textit {ph}}}$. By varying the amplitude of the gradient, rather than the duration, exactly the same sampling of the echo is achieved but each point has undergone the same degree of ${\displaystyle T_{1}}$ and ${\displaystyle T_{2}^{*}}$ relaxation (fig.phase_encode_c)~\cite{edelstein80}.

#### Previous Text

A number of schemes have been devised for combining field gradients and radiofrequency excitation to create an image. One involves 2D or 3D reconstruction from projections, much as in Computed Tomography. Others involve building the image point-by-point or line-by-line. One even uses gradients in the rf field rather than the static field. Although each of these schemes is occasionally used in specialist applications, the majority of MR Images today are created either by the Two-Dimensional Fourier Transform (2DFT) technique with slice selection, or by the Three-Dimensional Fourier Transform (3DFT) technique. Another name for 2DFT is spin-warp. What follows here is a description of the 2DFT technique with slice selection.

Slice selection is achieved by applying a magnetic gradient in addition to the external magnetic field during the radio frequency pulse. Only one plane within the object will have protons that are on–resonance and contribute to the signal.

A real image can be considered as being composed of a number of spatial frequencies at different orientations. A two–dimensional Fourier transformation of a real image will express these waves as a matrix of spatial frequencies known as k–space. Low spatial frequencies are represented at the center of k–space and high spatial frequencies at the periphery. Frequency and phase encoding are used to measure the amplitudes of a range of spatial frequencies within the object being imaged.

The frequency encoding gradient is applied during readout of the signal and is orthogonal to the slice selection gradient. During application of the gradient the frequency differences in the readout direction progressively change. At the midpoint of the readout these differences are small and the low spatial frequencies in the image are sampled filling the center of k-space. Higher spatial frequencies will be sampled towards the beginning and end of the readout filling the periphery of k-space.

Phase encoding is applied in the remaining orthogonal plane and uses the same principle of sampling the object for different spatial frequencies. However, it is applied for a brief period before the readout and the strength of the gradient is changed incrementally between each radio frequency pulse. For each phase encoding step a line of k–space is filled.

Either a spin echo or a gradient echo can be used to refocus the magnetisation.

The 3DFT technique is rather similar except that there is no slice selection and phase-encoding is performed two separate directions.

Another scheme which is sometimes used, especially in brain scanning or where images are needed very rapidly, is called echo-planar imaging (EPI): in this case each rf excitation is followed by a whole train of gradient echoes with different spatial encoding.

#### Example of a Pulse Sequence

Simplified timing diagram for two-dimensional-Fourier-transform (2DFT) pulse sequence