Linear Gradient Characterization Using A Rapid Thin-slice Measurement
Ryan K Robison1 and James G Pipe1

1Imaging Research, Barrow Neurological Institute, Phoenix, AZ, United States

Synopsis

Gradient impulse response functions can predict gradient behavior in the presence of eddy currents and other sources of error. This work presents a rapid method for generating the 0th and 1st order gradient impulse response functions using existing software measurement techniques. The proposed method is applied to spiral imaging but could be similarly applied to EPI, projection reconstruction, or any other imaging sequence.

Introduction

This work describes a method for measuring the 0th and 1st order gradient impulse response functions (GIRF) using existing intrinsic MR methods and for applying these response functions to accurately estimate the k-space trajectory and B0 phase errors in the presence of eddy currents and channel delays. Measurement of the GIRF is rapid (< 3 minutes) and needs only be done once for a given gradient system. This method will be shown to reduce gradient related artifacts in the context of 2D spiral imaging.

Methods

Two-dimensional multi-slice Spiral and Cartesian data were acquired on a healthy volunteer. For the spiral acquisition, previously determined channel delays and gradient corrections were disabled to highlight the proposed method (eddy current correction was not disabled). The process illustrated in figure 1 was used to measure the B0 eddy current response, H0, and the linear gradient field response, H1. An impulse response measurement framework was used as set forth by Vannesjo et al1, but the actual measurement was accomplished with a birdcage coil and a uniform spherical phantom using the B0 eddy current measurement method of Brodsky et al2.

The phase produced by B0 eddy currents can be predicted according to

$$\hat \theta_{B0}(\omega) = G_x(\omega) H_{0x}(\omega) + G_y(\omega) H_{0y}(\omega) + G_z(\omega) H_{0z}(\omega),$$

and the actual k-space trajectory can be predicted according to

$$\hat k_l(\omega) = k_l(\omega) H_{1l}(\omega),$$

where $$$\hat \theta_{B0}$$$ is the predicted B0 eddy current induced phase, $$$H_{0l}$$$ and $$$H_{1l}$$$ are the measured 0th and 1st order GIRFs and $$$l$$$ is one of the three physical gradient axes $$$l \in \{x, y, z\}$$$, $$$G_l$$$ are the theoretical gradient waveforms, and $$$\hat k_l$$$ and $$$k_l$$$ are the predicted and theoretical k-space coordinates. Artifact reduction was achieved by applying the conjugate of the predicted B0 phase to the spiral data and by reconstructing this data using the predicted k-space coordinates.

To validate the predicted k-space coordinates, the actual k-space trajectory was independently measured using the self-encoding gradient method of Alley et al3. This method requires several minutes to measure all interleaves of a 2D spiral trajectory but is very accurate. Similarly, the B0 eddy current induced phase from the spiral trajectory was measured independently using the method of Brodsky et al2. All data processing and reconstruction was performed in GPI4.

Results

Spiral data reconstructed with and without correction are shown in figure 2. Cartesian data are also shown as a reference. As is evident from the figure, the B0 eddy current phase and k-space trajectory predicted using the measured H0 and H1 GIRFs provide substantial improvements in image quality. The predicted k-space trajectory is compared in figure 3 to the designed theoretical spiral trajectory and the independently measured spiral trajectory. As can be seen from the figure, the predicted trajectory very closely follows the measured trajectory, and both differ from the theoretical trajectory. Figure 4 compares the independently measured B0 eddy current induced phase to the predicted phase. The predicted phase follows the shape of the measured phase.

Discussion and Conclusion

This work has presented a framework for predicting the effects of eddy currents, channel delays, and other gradient errors using a fast, one-time measurement. It should be noted that the k-space trajectory measurements extracted using this method also benefit in accuracy from the removal of the phase associated with B0 eddy currents. To our knowledge, the effects of B0 eddy current induced phase have not been taken into account in any of the previously published methods that measure k-space using a thin-slice excitation.

The predicted B0 eddy current phase and k-space trajectory follow the measured trajectories very closely. The match is less good for the B0 eddy current phase; however, the match is close enough that image quality is similar whether the measured or predicted phase is used to correct the spiral data.

While the application here was spiral MRI, the measured GIRFs can be applied with equal efficacy to projection reconstruction MRI, EPI, or any other acquisition as long as the theoretical gradient waveforms are known.

Acknowledgements

This work was funded in part by Philips Healthcare

References

1. Vannesjo SJ, Haeberlin M, Kasper L, Pavan M, Wilm BJ, Barmet C, and Preussmann KP, "Gradient System Characterization by Impulse Response Measurements with a Dynamic Field Camera", MRM (2013), 69:583-593.

2. Brodsky EK, Klaers JL, Samsonov AA, Kijowski R, and Block WF, "Rapid Measurement and Correction of Phase Errors From B0 Eddy Currents: Impact on Image Quality for Non-Cartesian Imaging", MRM (2013), 69:509-515.

3. Alley MT, Glover GH, and Pelc NJ, "Gradient Characterization Using a Fourier-Transform Technique", MRM (1998), 39:581-587.

4. Zwart NR and Pipe JG, "Graphical Programming Interface: A Development Environment for MRI Methods", MRM (2015), 74:1449-1460.

Figures

Figure 1: Proposed method for measuring the gradient field responses in the presence of eddy currents and gradient imperfections. Triangular gradient waveforms are played on each physical gradient. Measurement data are processed to yield B0 eddy current and linear gradient field responses, which subsequently yield gradient impulse response functions.

Figure 2: Spiral data (scan time = 29”) were reconstructed without (a) and with (b) B0 eddy current and trajectory corrections as estimated using the proposed method. Cartesian data (c) (scan time = 1’14”) are also shown as a reference.

Figure 3: Comparison of the kx component of one interleave of a spiral k-space trajectory produced by the theoretical spiral design, the measured spiral k-space trajectory, and the predicted response found by applying the measured H1x response to the theoretical spiral trajectory.

Figure 4: Comparison of measured and predicted B0 eddy current phase. The predicted phase was produced by applying the measured H0 response to the designed spiral gradient waveforms.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
4283