Measuring Magnetic Field Inhomogeneity From Spatial Distortion of Echo Planar Images
Peter Andrew Hardy1 and Erfan Akbari2

1Radiology, University of Kentucky, Lexington, KY, United States, 2Radiation Medicine, University of Kentucky, Lexington, KY, United States

Synopsis

We developed an imaging and analysis method to use the spatial distortion in echo planar images to estimate the magnetic field inhomogeneity. The method requires the acquisition of two images from an unmodified echo planar sequence with the direction of the phase encode reversed between the two. After combining the images appropriately an undistorted and a displacement image are produced. The displacement image is a map of the local magnetic field inhomogeneity.

Purpose

Our purpose was to develop an imaging method which could quickly estimate the homogeneity of the magnetic field using readily available, echo planar imaging techniques. This method is valuable both clinically and for routine quality assurance of the imager even when double echo gradient phase map imaging is not available.

Method

On a 1.5T Siemens Aera MR imager we acquired echo planar and gradient echo images of the ACR quality control phantom. We deliberately adjusted the linear and quadratic shim gradients to create known magnetic field inhomogeneities. For each condition, we acquired images with the phase encode gradient oriented in both positive and negative directions. We combined these two images by integrating the image intensity along each column in the phase encoding direction to produce one combined, distortion-free image and a second image corresponding to the amount of shift required to remove the distortion.3 Knowing the bandwidth in the phase encoding direction of the echo planar imaging technique we easily converted this image into a measurement of magnetic field inhomogeneity to within approximately 5%.

Introduction

Inhomogeneous magnetic fields leads to several problems in MRI including image geometric distortion, non-uniformity in image intensity and inhomogeneous fat suppression.2 As a result the American College of Radiology asks for regular measurements of field homogeneity. Additionally, measuring field homogeneity is valuable in clinical situations to optimize image performance. When available the best method to measure field homogeneity is a phase contrast technique where the phase in a gradient echo sequence is measured at two echo times. In many cases, however, a phase contrast technique is not readily available except to the Field Service Engineer. We developed a method using echo planar images whose image fidelity is very sensitive to field homogeneity. We also developed image analysis methods to automatically estimate the distortion and from that infer the field inhomogeneity.

Results

Figure 1 shows the measured inhomogeneity derived from the measurement of the dimension of the phantom in the phase encode direction versus the nominal field inhomogeneity derived from the deliberate misadjustment of the Gy linear gradient. Figure 2 shows in each row, the composite collection of echo planar images acquired with the phase encode gradient oriented down or up, as well as an image derived from the distorted images and, a an image representing the displacement necessary to transform the distorted image into the undistorted-combined image. This set of four images is repeated for the three settings of the quadratic shim gradient, Gxy=X2-Y2 of -50 µT/m2, 0, or +50 µT/m2. Figure 3 is a plot of the field inhomogeneity expressed as frequency offset in Hz. Note the expected quadratic shape of the field inhomogeneity arising from the X2-Y2 gradient. Figure 4 plots the vertical profiles taken through the displacement images shown in figure 2 taken through the phantom. Plots are given for the conditions of Gxy= -50 µT/m2, 0 , or +50 µT/m2.

Discussion

The EPI method was successful in estimating linear field inhomogeneities accurately. By swapping the direction of the phase encoding gradient and acquiring two EPI images in rapid succession one can estimate the in-plane distortion in the phase encoding direction and hence, the field homogeneity over a large volume. The total acquisition time is on the order of 2 TR which can be as low as 5 sec. The method can also be applied clinically where it may be used to estimate the magnetic field homogeneity where movement from cardiac or respiratory motion would normally corrupt the measurement.

Conclusion

The echo planar method for measuring field homogeneity can produce accurate estimates of in-plane magnetic field homogeneity throughout a large volume. The measurement time is short and the measurements can be performed on most MR imagers without special pulse sequences.

Acknowledgements

This work was supported by the department of Radiology at the University of Kentucky

References

1. P. Jezzard, R.S. Balaban, "Correction for geometric distortion in echo planar images from B0 field variations," Magnetic resonance in medicine : official journal of the Society of Magnetic Resonance in Medicine / Society of Magnetic Resonance in Medicine 34, 65-73 (1995).

2. E.M. Haacke, R. Brown, M. Thompson, R. Venkatesan, Magnetic resonance imaging: physical principles and sequence design, 1st ed. (John Wiley & Sons, New York, 1999).

3. H. Chang, J.M. Fitzpatrick, "A technique for accurate magnetic resonance imaging in the presence of field inhomogeneities," IEEE transactions on medical imaging 11, 319-329 (1992).

Figures

Figure 1. Magnetic field inhomogeneity derived from measuring the dimension of the ACR phantom in the phase encoding direction as it varied when the linear gradient Gy was misadjusted.

Figure3. Vertical profile across a phase contrast gradient echo sequence acquired on the ACR phantom. Blue profile acquired with quadratic gradient X2-Y2 offset by – 50 µT/m2. Red profile acquired with quadratic gradient X2-Y2 offset by + 50 µT/m2 .

Figure 4. Vertical profile across a reconstructed echo planar image acquired on the ACR phantom. Blue profile acquired with quadratic gradient X2-Y2 offset by – 50 µT/m2. Red profile acquired with quadratic gradient X2-Y2 offset by + 50 µT/m2.

Figure 2. Images of the ACR QC phantom acquired with an echo planar imaging technique in the presence of varying magnetic field inhomogeneities and with different directions of the phase encode gradient.



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