Transverse Relaxometry with B1+ Constrained Stimulated Echo Correction
Reza Basiri1, Marc Lebel2, and Paolo Federico3

1Biomedical Enginnering, University of Calgary, Calgary, AB, Canada, 2Alberta Children's Hospital Research Institute, Calgary, AB, Canada, 3Foothills Hospital, Calgary, AB, Canada

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 (T2), has shown promise in detection of subtle abnormalities in neurological diseases1. However, historically, this technique has produced inconsistent results, due primarily to sub-optimal data fitting. Quantitative T2 mapping is highly sensitive to transmit field heterogeneities (B1+) 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 fitting3, but this risks compromising tissue characterization4. 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 (T2, amplitude, and B1+) are estimated. While T2 and amplitude are not necessarily guaranteed to have strong spatial correlations, the B1+ 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 B1+ in order to reduce inconsistencies associated with T2 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 T2, amplitude, and B1+. The B1+ field is then spatially filtered by convolving with a Gaussian window with 3mm FWHM. The second pass uses the filtered B1+ maps as a pre-estimated input and re-fits the other 2 parameters (T2 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, T1/T2 of 3000/100 ms, and B1+ 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 T2 maps generated with both methods were evaluated based on variance across trials.

Results

Simulations indicate that both SEC and iSEC methods provide equally accurate T2 estimation over a wide range of SNR, Figure 1. However, the restricted B1+ field in the iSEC method translates to a reduced variance at all SNR values: T2 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 B1+, Figure 2. With the nominal parameters used in the simulation, tissues with T2 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 T2 maps illustrate relatively similar average T2 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 T2 sensitivity to mis-estimates of B1+ are minor and less than 10% reduction in variance little benefit is observed. Overall, simulated and experimental results suggest that a B1+ 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 T2 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 T2 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 B1+ field, where we observe benefit in a B1+ constrained fit. Ultimately, iSEC is expected to translate into more reliable T2 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 T2 relaxometry-aided diagnoses have already been shown1.

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.

Figures

Figure 1: Average (solid central lines) +- standard deviation (shaded gray areas) of the percent error in estimated T2 measurements at different noise realizations. Overall, iSEC is predicted to reduce fitting variance by 25% relative to the original SEC.

Figure 2: iSEC reduces variance in the estimated T2 value relative to SEC fitting. The benefit is most important in regions where the transmit field is low (below about 0.8)

Figure 3: iSEC reduces variance in the estimated T2 value relative to SEC fitting, especially in tissues with T2 values below 110 ms, which corresponds to white and gray matter at 1.5 and 3.0 T.

Figure 4: Similar T2 maps in terms of accuracies of the estimate. Improved precision ~24% in regions with low B1+ values are shown in red. PDF graph demonstrates precision improvement when standard devision are compared (σSEC - σiSEC)



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