Morphology-adaptive total variation for the reconstruction of quantitative susceptibility map from the magnetic resonance imaging phase

Quantitative susceptibility mapping (QSM) is a magnetic resonance imaging technique that quantifies the magnetic susceptibility distribution within biological tissues. QSM calculates the underlying magnetic susceptibility by deconvolving the tissue magnetic field map with a unit dipole kernel. However, this deconvolution problem is ill-posed. The morphology enabled dipole inversion (MEDI) introduces total variation (TV) to regularize the susceptibility reconstruction. However, MEDI results still contain artifacts near tissue boundaries because MEDI only imposes TV constraint on voxels inside smooth regions. We introduce a Morphology-Adaptive TV (MATV) for improving TV-regularized QSM. The MATV method first classifies imaging target into smooth and nonsmooth regions by thresholding magnitude gradients. In the dipole inversion for QSM, the TV regularization weights are a monotonically decreasing function of magnitude gradients. Thus, voxels inside smooth regions are assigned with larger weights than those in nonsmooth regions. Using phantom and in vivo datasets, we compared the performance of MATV with that of MEDI. MATV results had better visual quality than MEDI results, especially near tissue boundaries. Preliminary brain imaging results illustrated that MATV has potential to improve the reconstruction of regions near tissue boundaries.