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 mm
3; FOV=55x220x220
mm
3 (S1), 46x280x280 mm
3 (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 mm
2
in-plane, with 8 sagittal slices of thickness 3.0 mm, for an effective FOV=24x176x282
mm
3. 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
maps
2 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
cm
3 and 15 cm
3 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.