Spatially Adaptive Regularization in Total Field Inversion for Quantitative Susceptibility Mapping

Adaptive Total Field Inversion is described for quantitative susceptibility mapping (QSM) reconstruction from total field data through a spatially adaptive suppression of shadow artifacts through spatially adaptive regularization. The regularization for shadow suppression consists of penalizing low-frequency components of susceptibility in regions of small susceptibility contrasts as estimated by R2∗ derived signal intensity. Compared with a conventional local field method and two previously proposed regularized total field inversion methods, improvements were demonstrated in phantoms and subjects without and with hemorrhages. This algorithm, named TFIR, demonstrates the lowest error in numerical and gadolinium phantom datasets. In COSMOS data, TFIR performs well in matching ground truth in high-susceptibility regions. For patient data, TFIR comes close to meeting the quality of the reference local field method and outperforms other total field techniques in both clinical scores and shadow reduction.