High resolution myelin water imaging incorporating local tissue susceptibility analysis

Quantitative myelin water imaging (MWI) from signal T2* decay acquired with multiple Gradient-Recalled Echo (mGRE) sequence has been widely used since its first report. A recent study showed that with low resolution data (2 mm isotropic voxels), direct application of complex fitting to a three-pool WM model with frequency shift terms could produce more stable parameter estimation for myelin water fraction mapping. MWI maps of higher spatial resolution resulting in more detailed tissue structures and reduced partial volume effects around white matter/gray matter (WM/GM) interface, however, is more desirable. Furthermore, as signal-to-noise ratio (SNR) of original images decreases due to reduced voxel size, the direct complex fitting procedure of myelin water imaging becomes more prone to systematic errors which severely compromised stability and reliability of the result. Instead of using the original part of T2* decay, this work presents a new method based on the WM-induced phase from tissue susceptibility calculated with the same mGRE dataset, in a three-pool WM model (water of myelin, axonal and extracellular water), to improve high resolution MWI. Compared with direct complex fitting for the higher spatial resolution case, the proposed method is shown to provide a more stable and accurate estimation of MWI parameters, and finer details near WM/GM boundaries with greatly reduced partial volume effects.