Validation of Data Acquisition and Phase Estimation for Quantitative Susceptibility Mapping with a Rotating-Tube Phantom

Abstract

AbstractPurposeQuantitative Susceptibility Mapping (QSM) is an MRI tool with the potential to reveal pathological changes from magnetic susceptibility measurements. Before phase data can be used to recover susceptibility (Δχ), the QSM process begins with two steps: data acquisition and phase estimation. We assess the performance of these steps, when applied without user intervention, on several variations of a phantom imaging task.ApproachWe used a rotating-tube phantom with five tubes ranging from Δχ=0.05 ppm to Δχ=0.336 ppm. MRI data was acquired at nine angles of rotation for four different pulse sequences. The images were processed by 10 phase estimation algorithms including Laplacian, region-growing, branch-cut, temporal unwrapping and maximum-likelihood methods. We analyzed errors between measured and expected phase using the probability mass function and Cumulative Distribution Function.ResultsRepeatable acquisition and estimation methods were identified based on the probability of relative phase errors. For single-echo GRE and segmented EPI sequences, a region-growing method was most reliable with Pr(relative error<0.1)=0.95 and 0.90 respectively. For multi-echo sequences, a Maximum-Likelihood method was most reliable with Pr(relative error<0.1)=0.97. The most repeatable multi-echo methods outperformed the most repeatable single-echo methods.ConclusionsWe found a wide range of repeatability and reproducibility for off-the-shelf MRI acquisition and phase estimation approaches. The error was dominated in many cases by spatially discontinuous phase unwrapping errors. Any post-processing applied on erroneous phase estimates, such as QSM’s background field removal and dipole inversion, would suffer from error propagation. Our paradigm identifies methods that yield consistent and accurate phase estimates that would ultimately yield consistent and accurate Δχ estimates.