A multicomponent T2 relaxometry algorithm for myelin water imaging of the brain.

Abstract

Models based on a sum of damped exponentials occur in many applications, particularly in multicomponent T2 relaxometry. The problem of estimating the relaxation parameters and the corresponding amplitudes is known to be difficult, especially as the number of components increases. In this article, the commonly used non-negative least squares spectrum approach is compared to a recently published estimation algorithm abbreviated as Exponential Analysis via System Identification using Steiglitz-McBride. The two algorithms are evaluated via simulation, and their performance is compared to a statistical benchmark on precision given by the Cramér-Rao bound. By applying the algorithms to an in vivo brain multi-echo spin-echo dataset, containing 32 images, estimates of the myelin water fraction are computed. Exponential Analysis via System Identification using Steiglitz-McBride is shown to have superior performance when applied to simulated T2 relaxation data. For the in vivo brain, Exponential Analysis via System Identification using Steiglitz-McBride gives an myelin water fraction map with a more concentrated distribution of myelin water and less noise, compared to non-negative least squares. The Exponential Analysis via System Identification using Steiglitz-McBride algorithm provides an efficient and user-parameter-free alternative to non-negative least squares for estimating the parameters of multiple relaxation components and gives a new way of estimating the spatial variations of myelin in the brain.