Synopsis
Quantitative T2
mapping provides diagnostic capabilities complementing standard qualitative
imaging. However, conventional fitting algorithms to estimate T2 are
prone to bias. In this work, we propose a fitting method that remains
applicable to existing datasets while addressing many of the imperfections and
shortcomings of current methods. Our proposed method is an extension of
stimulated echo correction that highly constrains the estimated transmit field.
It was evaluated using simulated and experimental data. We found that variance
in the T2 estimate could be reduced by ~25% in certainly realistic
conditions while maintaining full accuracy relative to the current stimulated
echo corrected fit.
Transverse
relaxometry, a quantitative T2 mapping has shown superior diagnostic
capabilities compare with qualitative maps for neurological diseases. However,
the conventional fitting
Quantitative T2
mapping provides diagnostic capabilities complementing standard qualitative
imaging. However, conventional fitting algorithms to estimate T2 are
prone to bias. In this work, we propose a fitting method that remains
applicable to existing datasets while addressing many of the imperfections and
shortcomings of current methods. Our proposed method is an extension of
stimulated echo correction that highly constrains the estimated transmit field.
It was evaluated using simulated and experimental data. We found that variance
in the T2 estimate could be reduced by ~25% in certainly realistic
conditions while maintaining full accuracy relative to the current stimulated
echo corrected fit.Purpose
Transverse relaxometry, a quantitative MRI technique that
measures the signal decay time (T
2), has shown promise in detection
of subtle abnormalities in neurological diseases
1. However,
historically, this technique has produced inconsistent results, due primarily
to sub-optimal data fitting. Quantitative T
2 mapping is highly
sensitive to transmit field heterogeneities (B
1+)
2,
and in the case of slice selective imaging, to radio-frequency (RF) pulse shapes
which cause imperfections in signals generated with a multi-echo spin echo
sequence. Conventional fitting methods ignore transmit imperfections (and the
resulting echo oscillations). One strategy is to simply discard early echo
times to improve fitting
3, but this risks compromising tissue
characterization
4. A recently proposed fitting method, called
stimulated echo correction (SEC)
5, estimates major confounds
associated with fitting errors in the transmit field and returns less biased
results. SEC is a one-step least square method in which 3 parameters (T
2,
amplitude, and B
1+) are estimated. While T
2
and amplitude are not necessarily guaranteed to have strong spatial
correlations, the B
1+ field is not expected to vary
rapidly and can be highly constrained. Our aim is to develop and employ an
improved SEC fit (called iSEC) that constrains B
1+ in
order to reduce inconsistencies associated with T
2 estimation.
Methods
Our
proposed iSEC method is a two-step procedure: the first pass consists of the
standard SEC method to provide an initial estimate of T
2, amplitude,
and B
1+. The B
1+ field is then spatially
filtered by convolving with a Gaussian window with 3mm FWHM. The second pass
uses the filtered B
1+ maps as a pre-estimated input and
re-fits the other 2 parameters (T
2 and amplitude). We compared and investigated
reliability between the proposed iSEC and the standard SEC fit with simulated
and in-vivo data. Simulated data were
generated with the extended phase graph algorithm assuming nominal acquisition
and tissue parameters: 16 echoes, 10 ms echo spacing, T
1/T
2
of 3000/100 ms, and B
1+ reduced to 0.75 of its ideal
value. Gaussian noise (10<SNR<100) was added to mimic real data; 1000
noise realizations were performed at each SNR value. T2 values were
estimated using the original SEC and our proposed iSEC algorithm. Experimental
verification was performed by acquiring 10 sequentially repeated scans in one
volunteer. The sequence consisted of a multi-echo spin echo with 12 slices, 16
echoes, and 9.2 ms echo spacing on a 3T GE MR750 scanner T
2 maps generated with
both methods were evaluated based on variance across trials.
Results
Simulations indicate that both SEC and iSEC methods provide
equally accurate T
2 estimation over a wide range of SNR, Figure 1.
However, the restricted B
1+ field in the iSEC method
translates to a reduced variance at all SNR values: T
2 values estimated
with iSEC had ~25% lower standard deviation than those estimated with SEC. iSEC
was found to be particularly beneficial in low SNR regions (<35) as well as
within regions with low B
1+, Figure 2. With the nominal
parameters used in the simulation, tissues with T
2 values below
approximately 110 ms (corresponding to white and gray matter
6)
benefited the most with our proposed fit, Figure 3. Our experimental MR images
also produce similar results. Estimated T
2 maps illustrate
relatively similar average T
2 values for regions with high SNR (>
50) for both SEC and iSEC. Repeated scans indicate that standard deviations within
the regions with low B1+ fraction (about 0.75) are
lowered 24% in iSEC compared to SEC, Figure 4. In regions with near-ideal B1+
(above 0.9), the T
2 sensitivity to mis-estimates of B
1+
are minor and less than 10% reduction in variance little benefit is observed. Overall,
simulated and experimental results suggest that a B
1+
constrained fit is able to improve fits under all circumstances but is
especially beneficial when the transmit field is low.
Discussion and Conclusions
Simulated and experimental data indicate that precision in T
2
estimation is always improved with iSEC relative to SEC. The improvement
observed with iSEC is greatest in low SNR regions and in situations when system
imperfections are severe. This translates to reliable T
2 estimations
with thinner slices or with higher spatial resolution. This may also be
important if the refocusing angles are intentionally lowered to reduce SAR or
to achieve shorter echo spacing. This situation will mimic an unintentionally
low B
1+ field, where we observe benefit in a B
1+
constrained fit.
Ultimately,
iSEC is expected to translate into more reliable T
2 maps than can
currently be generated, which may provide reliable detection of pathology in
neurological disorders particularly temporal lobe epilepsy in which extensive
advantages of T
2 relaxometry-aided diagnoses have already been shown
1.
Acknowledgements
1. University of Calgary, AB Canada
2. Seaman MRI Center, Foothills Hospital, AB Canada
3. The Natural Sciences and Engineering Research Council of Canada (NSERC)
References
1. Sumar, I., Kosior, R., Frayne, R., & Federico, P.
(2011). Hippocampal T2 abnormalities in healthy adults. Epilepsy Research,
273-276.
2. Collins CM, Liu W, Schreiber W, Yang QX, Smith MB.
(2005). Central brightening due to constructive interference with, without, and
despite dielectric resonance. J Magn Reson Imaging, 21:192–196.
3. Maier CF, Tan SG, Hariharan H, Potter HG. (2003). T2
quantitation of articular cartilage at 1.5 T. J Magn Reson Imaging,
17:358–364.
4. Crawley, A. P. and
Henkelman, R. M. (1987), Errors in T2 estimation using multislice
multiple-echo imaging. Magn Reson Med, 4: 34–47.
5. Lebel, R., & Wilman, A. (2010). Transverse
relaxometry with stimulated echo compensation. Magnetic Resonance in
Medicine, 1005-1014.
6. Wansapura, J., Holland, S., Dunn, R., & Ball, W.
(1999). NMR relaxation times in the human brain at 3.0 tesla. J Magn Reson
Imaging, 531-538.