Portable Imaging with Rotating Inhomogeneous Magnetic Fields
Clarissa Zimmerman Cooley1

1A. A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital, Charlestown, MA, United States

Synopsis

As the premiere modality for brain imaging, MRI could find wider applicability if lightweight, portable systems were available for siting in unconventional locations. We construct and validate a truly portable (<100kg) and silent proof-of-concept scanner which replaces conventional gradient encoding with a rotating inhomogeneous low-field magnet. When rotated about the object, the inhomogeneous field pattern is used to create generalized projections. The system is validated with experimental 2D images, and extended to 3D imaging with the addition of Transmit Array Spatial Encoding (TRASE). This new scanner architecture demonstrates the potential for portability by simultaneously relaxing the magnet homogeneity criteria and eliminating the gradient coil.

Target Audience

MRI researchers interested in methods for portable imaging

Methods

A small rotating permanent magnet array with an inhomogeneous field is used to create a portable MRI scanner. The field variation of the inhomogeneous magnet is leveraged as the encoding field in the Y-Z plane instead of using additional gradient coils. This eliminates cooling requirements, and vastly reduces the weight and power requirements of the scanner.

The previously described 45 kg (36 cm dia. 35.6 cm length) Halbach cylinder magnet [1,2] (Fig. 1) uses 20 rungs of NdFeB magnets to create a 77.3 mT field (3.29 MHz proton freq.) with 32 kHz variation in the center slice 16 cm dia. FOV (mapped using small field probes). As the magnet is physically rotated around the object using a stepper motor, a RARE-type spin echo sequence (40 kHz, 256pt readout) records the generalized projections onto the nonlinear field pattern. The individual echoes are averaged to improve SNR.

The non-bijectivity of the approximately quadrupolar spatial encoding magnetic field (SEM) causes image aliasing that we disambiguate using parallel imaging as described by Schultz et al.[3]. We use an 8-channel receive array (Fig. 2a) with 8 cm circular loops for this purpose. The coil sensitivity map (Fig. 2b) is modeled using the Biot-Savart law and is different for each magnet rotation angle. The coil sensitivity profiles and field maps at each rotation angle are used to form the encoding matrix. We solve for 2D images (in the Y-Z plane) using the Algebraic Reconstruction Technique (ART) and Preconditioned Conjugate Gradients [2,4,5,6].

For 3D imaging, encoding is needed along the axis of the cylinder (X direction) as well. Transmit Array Spatial Encoding (TRASE) is a complementary, B1 imaging method [7,8], that can be used for 1D encoding in the X direction of the described scanner. To enable TRASE encoding in the inhomogeneous magnetic field, broadband WURST pulses were used [9]. The combination of the rotating SEM method for 2D encoding and the TRASE method for 3rd axis encoding, enables a 3D imaging method in the portable scanner [10].

Results

Figure 3 shows experimental images of a “MIT/MGH” phantom (CuSO4-doped water, 1.7 cm thick, 13 cm dia.) acquired with 7 Rx coils, 32 averages of a 6 spin-echo train (TR = 550 ms, echo spacing = 8 ms) using 91 magnet rotations of 2o. Data is also acquired from a single field probe at each rotation to track field drift for use during reconstruction (Fig. 3b). The acquisition time of 66 minutes results from serially acquiring data from the 8 coils into a single channel on the console and would be reduced to 7.3 minutes with true parallel acquisition. Figure 3c shows a 1cm thick lemon slice imaged using 5 Rx coils, 128 echo train (TR = 4500 ms, echo spacing = 8 ms), and 181 magnet rotations of 1°. The total acquisition time was 93 minutes (15.5 minutes w/ parallel Rx).

Figure 4 shows experimental 3D imaging results using a phantom with 3 water-filled compartments spaced 2 cm apart in X. A 1D projection along X shows the water-filled compartments. Three slices are reconstructed from data at X = -2cm, 0cm, and 2cm using the corresponding 2D B0 field maps.

Discussion

The nonlinear SEM of the Halbach magnet results in variable resolution over the FOV, including notable blurring in the center where the encoding field is shallow (Fig. 3c). If a sufficiently-strong linear term was added to the SEM, the spatially uniform encoding field region would not coincide with the axis of rotation. In this case, the “blind-spot” moves around the object resulting in less severe blurring [11] (an idea that guides the design of future magnets for this system). Field map errors from temperature drifts are significant but can be measured and included in the encoding matrix and mitigated in the reconstructed image (Fig. 3b).

Conclusion

Using an inhomogeneous rotating magnet for spatial encoding in lieu of gradient coils, we have constructed and demonstrated a lightweight scanner for 2D MR imaging with minimal power requirements. The 2D proof-of-concept images from this nearly head-sized imager show the ability of this encoding scheme to produce sufficient spatial resolution and sensitivity for the detection and characterization of many common neurological disorders such as hydrocephalus and traumatic space-occupying hemorrhages. Future work perfecting the calibration methods is likely to bring experimental image quality closer to the theoretical limit, but the resolution of the current system is sufficient for identifying gross pathologies.

Acknowledgements

Lawrence Wald, Jason Stockmann, Matthew Rosen, Bastien Guerin, Matthieu Sarracanie, Cristen LaPierre, Brandon Armstrong, Steven Cauley, Melissa Haskell, Charlotte Sappo, James Blau. Support by NIH R01EB018976.

References

[1] Zimmerman C, Blau J, Rosen MS, Wald LL. Design and construction of a Halbach array magnet for portable brain MRI. In Proc. of the ISMRM, Melbourne, Australia, 2012, p. 2575

[2] Cooley CZ, Stockmann JP, Armstrong BD, Sarracanie M, Lev MH, Rosen MS, Wald LL. Two-dimensional imaging in a lightweight portable MRI scanner without gradient coils. Magn Reson Med. 2015 Feb;73(2):872-83.

[3] Gerrit Schultz, Daniel Gallichan, Marco Reisert, Jürgen Hennig, and MaximZaitsev. MR image reconstruction from generalized projections. Magn. Reson.Med., 2013.

[4] R Gordon, R Bender, and G T Herman. Algebraic reconstruction techniques(ART) for three-dimensional electron microscopy and x-ray photography. J.Theor. Biol., 29(3):471–481, 1970.

[5] Magnus Hestenes and Eduard Stiefel. Methods of conjugate gradients for solving linear systems. Journal of Research of the National Bureau of Standards,49(6):409–436, 1952.

[6] Jason P Stockmann, Pelin Aksit Ciris, Gigi Galiana, Leo Tam, and R Todd Constable. O-space imaging: Highly efficient parallel imaging using second order nonlinear fields as encoding gradients with no phase encoding. Magn.Reson. Med., 64(2):447–456, 2010.

[7] Scott B King, D Yin, S Thingvold, Jonathan C Sharp, and Boguslaw Tomanek.Transmit array spatial encoding (TRASE): A new data acquisition method in MRI. in Proc. of the ISMRM, Seattle, Washington, USA, 2006, page 2628,2006.

[8] Jonathan C Sharp and Scott B King. MRI using radiofrequency magnetic fieldphase gradients. Magn. Reson. Med., 63(1):151–161, 2010.

[9] Stockmann JP, Cooley CZ, Sarracanie M, Rosen MS, Wald LL. Transmit Array Spatial Encoding (TRASE) with broadband WURST pulses for robust spatial encoding in inhomogeneous B0 fields. In Proc. of the ISMRM, Toronto, Canada, 2015, p. 0703.

[10] Cooley CZ, Stockmann JP, Rosen MS, Wald LL. 3D Imaging in a Portable MRI Scanner using Rotating Spatial Encoding Magnetic Fields and Transmit array spatial encoding (TRASE). In Proc. of the ISMRM, Toronto, Canada, 2015, p. 0703.

[11] Cooley CZ, Stockmann JP, Armstrong BD, Rosen MS, Wald LL. Spatial resolution in rotating spatial encoding magnetic field MRI (rSEM-MRI. In Proc. of the ISMRM, Milan, Italy, 2014, p. 0333.

Figures

Figure 1: (a) Halbach cylinder magnet on high friction rollers, (b) field map estimated as a polynomial fit to the measured field probe points (black dots).

Figure 2: (a) 8-channel RX array coil with disk phantom. (b) Calculated B1- profiles in the center plane of a surface coil located at the right side of the FOV. The arrows show three representative orientations of B0.

Figure 3: Experimental 16 cm FOV images. (a,b) 1.7cm thick phantom filled with CuSO4-doped water. (a) Temperature drift not corrected, (b) temperature drift corrected. (c) 1 cm thick lemon slice.

Figure 4: Image of a 3D phantom consisting of three water-filled compartments spaced cm apart in X. A 1D image along X for 1 rotation angle is shown (lower left). 3 image slices were reconstructed using the appropriate 2D field maps.



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