The aim of this study is to perform a thorough comparison of quantitative susceptibility mapping (QSM) techniques and their dependence on the assumptions made. The compared methodologies were: two iterative single orientation methodologies minimizing the l2, l1TV norm of the prior knowledge of the edges of the object, one overdetermined multiple orientation method (COSMOS) and a newly proposed modulated closed-form solution (MCF). The performance of these methods was compared using a numerical phantom and in-vivo high resolution (0.65 mm isotropic) brain data acquired at 7 T using a new coil combination method. For all QSM methods, the relevant regularization and prior-knowledge parameters were systematically changed in order to evaluate the optimal reconstruction in the presence and absence of a ground truth. Additionally, the QSM contrast was compared to conventional gradient recalled echo (GRE) magnitude and R2* maps obtained from the same dataset. The QSM reconstruction results of the single orientation methods show comparable performance. The MCF method has the highest correlation (corr MCF = 0.95, r 2 MCF = 0.97) with the state of the art method (COSMOS) with additional advantage of extreme fast computation time. The L-curve method gave the visually most satisfactory balance between reduction of streaking artifacts and over-regularization with the latter being overemphasized when the using the COSMOS susceptibility maps as ground-truth. R2* and susceptibility maps, when calculated from the same datasets, although based on distinct features of the data, have a comparable ability to distinguish deep gray matter structures. © 2014 Elsevier Inc. All rights reserved.
IntroductionPhase imaging has shown over the last decade to offer a good contrast, both between and within brain tissues in respect to the conventional magnitude signal (Duyn et al., 2007;Rauscher et al., 2005) as well as veins and iron rich regions (Haacke et al., 2004). The effect observed in the phase is known to be non-local, it reflects the magnetic field induced by the tissues' magnetic susceptibility (Marques and Bowtell, 2005), which scales linearly with the increase of the fields strength (making it suitable at high field strengths).Several studies have been performed on the origin of the susceptibility contrast with the main modulators being iron and myelin. Iron contributes to tissue contrast especially in the deep gray matter (globus pallidus, putamen and caudate) which has histologically derived high iron concentration showing good correlation with phase and susceptibility contrast (Bilgic et al., 2012;Schweser et al., 2011;Wharton and Bowtell, 2010). The other proposed contributor to the phase contrast, particularly between white and gray matter, is myelin where pathological demyelination has shown a decreased phase contrast between gray and white matter (C. Liu et al., 2011;Lodygensky et al., 2012) and good correlation was found between myelination and phase contrast during development. (Lodygensky et al., 2012).In addition to the no...