Very-high order shimming in the human spinal cord using a dedicated 24-channel array coil
Ryan Topfer1, Grégoire Germain1, Jason P. Stockmann2, Karl Metzemaekers3, Hoby Hetherington4, Raphaël Paquin5, Piotr Starewicz3, Nikola Stikov1,6, and Julien Cohen-Adad1,7

1NeuroPoly Lab, Institute of Biomedical Engineering, École Polytechnique de Montréal, Montreal, QC, Canada, 2Athinoula A. Martinos Center for Biomedical Imaging, Department of Radiology, Massachusetts General Hospital, Charlestown, MA, United States, 3Resonance Research Inc., Billerica, MA, United States, 4Department of Radiology, University of Pittsburgh, Pittsburgh, PA, United States, 5Siemens Healthcare Ltd., Montreal, QC, Canada, 6Montreal Heart Institute, Université de Montréal, Montreal, QC, Canada, 7Functional Neuroimaging Unit, CRIUGM, Université de Montréal, Montreal, QC, Canada

Synopsis

Pathologies of the spinal cord are a primary cause of functional disability and chronic pain. Although MRI already plays a role in the evaluation of these pathologies, it continues to be hampered by artifacts due to magnetic field inhomogeneity. This study reports the first results applying a specially designed 24-channel shim array to correct magnetic field inhomogeneity in the human spinal cord. Shimming using the custom array improved field homogeneity in the thoracic spinal cord of the two initial subjects by 53.6 % and 31.4 % respectively.

Purpose

Pathologies of the spinal cord, such as those due to trauma, neurodegenerative diseases, and cancers, are a primary cause for functional disability (paralysis) and chronic pain. Though MRI already plays a role in the evaluation of spinal cord pathology, it continues to be hampered by image artifacts owing to magnetic field inhomogeneity, underscored in a recent review as “the greatest challenge for acquiring MR images in the spinal cord.”1 Underlying the field inhomogeneity are a number of distinct local susceptibility inclusions such as cartilage, vertebral bones, fat, and the oxygen of the nearby airways and lungs, the latter of which acquires a time-varying component by virtue of the subject’s respiration. This study reports the first results applying a specially designed 24-channel shim array to image the human spinal cord.

Methods

Hardware: The shim system consists of 24 independently driven, rectangular-planar electrical coils, atop which sits the custom-built 8-channel transceiver array. These two coil arrays, first described in an earlier work,2 insert into the patient bed-table (Fig. 1) so as to lie in close proximity to the subject’s spine. Acquisition: For this initial proof-of-concept experiment, two healthy subjects (S1 & S2) were scanned on a 3T system (Magnetom Tim Trio, Siemens Healthcare, Erlangen, Germany). Standard 2nd-order shimming was first performed using the Siemens shims. High resolution T2-weighted anatomical images were obtained of the spinal cord with parameters: TE/TR=4.92/30 ms (S1), 3.88/15 ms (S2); flip angle 10°; isotropic voxel size 1.1 mm3; FOV=55x220x220 mm3 (S1), 46x280x280 mm3 (S2). For field mapping, a multi-echo gradient-echo sequence was acquired with parameters: TE=[4.92, 7.38, 11.13, 14.88] ms, TR=200 ms; flip angle 50°; spatial resolution=2.2x2.2 mm2 in-plane, with 8 sagittal slices of thickness 3.0 mm, for an effective FOV=24x176x282 mm3. To assess the static field correction specifically, Subject 2 was asked to inspire and maintain apnea during field map acquisition (approx. 10 s). Processing: To create binary masks for 3d phase unwrapping,3 the corresponding magnitude images were set to zero for intensities below 1 percent of the maximum. Unwrapped phase differences between the first two echoes were normalized by the echo time difference and scaled to Hz to obtain the field maps. The high resolution anatomical images were passed to the sct_propseg tool of the Spinal Cord Toolbox (http://www.neuro.polymtl.ca/downloads) to delineate the shim volumes of interest (VOI) automatically in about 30s.4 Spinal cord VOIs as well as the shim reference maps2 were interpolated to the grid spacing of the field maps prior to shim optimization. Shim optimization was rapidly performed using the fmincon function of the MATLAB optimization toolbox, minimizing the sum of the shim fields and the measured field map over the VOI while accounting for non-linear system constraints (e.g. max current per channel, and max total current).5 Optimal shim currents were set and the field maps were reacquired, during which Subject 2 was asked to once again maintain an inspired apnea.

Results

Results are shown in Fig. 2. The segment of the spinal cord targeted for shimming was approximately 7 cm3 and 15 cm3 in subjects 1 and 2 respectively. Shimming over these volumes using the the 24-channel array reduced the standard deviation of the field from 49.19 Hz to 22.85 Hz for Subject 1 (an improvement of 53.6 %), and from 13.65 Hz to 9.36 Hz for Subject 2 (an improvement of 31.4 %).

Discussion

The greater degree of field inhomogeneity observed for Subject 1 was likely due in part to less effective initial shimming using the standard shims in combination with (and perhaps more importantly) the presence of respiration-induced field variation. This respiration-induced variation was reported to be on the order of 100 Hz in the spinal cord across levels T1-T8 at 3.0 T.6 Work is underway to synchronize shim updates to the subject respiration to correct for these distortions in real-time. Though the results remain preliminary and limited to two subjects, the very-high order shim insert offers a promising technique to improve EPI-based quantitative MRI and spectroscopy of the spinal cord.

Acknowledgements

Funded by the Canadian Institute of Health Research [CIHR FDN-143263], the Fonds de Recherche du Québec - Santé [28826], the Fonds de Recherche du Québec - Nature et Technologies [2015-PR-182754], the Natural Sciences and Engineering Research Council of Canada [435897-2013], the Quebec BioImaging Network, Polytechnique MEDITIS, and NIH R21 EB017338.

References

1. P. W. Stroman, C. Wheeler-Kingshott, M. Bacon, J. M. Schwab, R. Bosma, J. Brooks, D. Cadotte, T. Carlstedt, O. Ciccarelli, J. Cohen-Adad, A. Curt, N. Evangelou, M. G. Fehlings, M. Filippi, B. J. Kelley, S. Kollias, A. Mackay, C. a Porro, S. Smith, S. M. Strittmatter, P. Summers, and I. Tracey, “The current state-of-the-art of spinal cord imaging: methods.,” Neuroimage, vol. 84, pp. 1070–81, Jan. 2014.

2. R. Topfer, K.-M. Lo, K. Metzemaekers, D. Jette, H. P. Hetherington, P. Starewicz, and J. Cohen-Adad, “A 24-channel shim array for real-time shimming of the human spinal cord: Characterization and proof-of-concept experiment,” in Proc Intl Soc Mag Reson Med, 2015, vol. 23, p. 3083.

3. H. S. Abdul-Rahman, M. a Gdeisat, D. R. Burton, M. J. Lalor, F. Lilley, and C. J. Moore, “Fast and robust three-dimensional best path phase unwrapping algorithm.,” Appl. Opt., vol. 46, no. 26, pp. 6623–35, Sep. 2007.

4. B. De Leener, S. Kadoury, and J. Cohen-adad, “Robust, accurate and fast automatic segmentation of the spinal cord,” Neuroimage, vol. 98, pp. 528–536, 2014.

5. J. P. Stockmann, T. Witzel, B. Keil, J. R. Polimeni, A. Mareyam, C. LaPierre, K. Setsompop, and L. L. Wald, “A 32-channel combined RF and B0 shim array for 3T brain imaging,” Magn Reson Med, Early View, 2015.

6. T. Verma and J. Cohen-Adad, “Effect of respiration on the B0 field in the human spinal cord at 3T.,” Magn Reson Med, vol. 72, no. 6, pp. 1629–36, 2014.

Figures

The 24-channel shim coil, with 8-channel transceiver atop it, inserts into the patient bed-table, in close proximity to the subject’s spine. Hydraulic cooling layers lie above and below the shim coils, which are powered independently of the scanner with a dedicated amplifier (MXD-32, Resonance Research Inc., Billerica, MA, USA.)

Methods and results. Standard shimming was performed on the scanner over the green box. The spinal cord was segmented to define the shim VOI used in optimizing the spine shim. Black arrows point to regions of the initial field map (second from right) where the field homogeneity has been markedly improved after spine shim optimization (right). Subject 2 was asked to maintain an inspired apnea for the initial field map acquisition, and to do so again immediately prior to the second acquisition.



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