Impact of pulse sequence, analysis method, and signal to noise ratio on the accuracy of intervertebral disc T2 measurement

Abstract

Noninvasive assessments of intervertebral disc health and degeneration are critical for addressing disc degeneration and low back pain. Magnetic resonance imaging (MRI) is exceptionally sensitive to tissue with high water content, and measurement of the MR transverse relaxation time, T2, has been applied as a quantitative, continuous, and objective measure of disc degeneration that is linked to the water and matrix composition of the disc. However, T2 measurement is susceptible to inaccuracies due to Rician noise, T1 contamination, and stimulated echo effects. These error generators can all be controlled for with proper data collection and fitting methods. The objective of this study was to identify sequence parameters to appropriately acquire MR data and to establish curve fitting methods to accurately calculate disc T2 in the presence of noise by correcting for Rician noise. To do so, we compared T2 calculated from the typical monoexponential (MONO) fits and noise corrected exponential (NCEXP) fits. We examined how the selected sequence parameters altered the calculated T2 in silico and in vivo. Typical MONO fits were frequently poor due to Rician noise, and NCEXP fits were more likely to provide accurate T2 calculations. NCEXP is particularly less biased and less uncertain at low SNR. This study showed that the NCEXP using sequences with data from 20 echoes out to echo times of ~300 ms is the best method for calculating T2 of discs. By acquiring signal data out to longer echo times and accounting for Rician noise, the curve fitting is more robust in calculating T2 despite the noise in the data. This is particularly important when considering degenerate discs or AF tissue because the SNR of these regions is lower.