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.