X-ray CT metal artifact reduction using wavelet domain L0 sparse regularization

X-ray computed tomography (CT) imaging of patients with metallic implants usually suffers from streaking metal artifacts. In this paper, we propose a new projection completion metal artifact reduction (MAR) algorithm by formulating the completion of missing projections as a regularized inverse problem in the wavelet domain. The Douglas-Rachford splitting (DRS) algorithm was used to iteratively solve the problem. Two types of prior information were exploited in the algorithm: 1) the sparsity of the wavelet coefficients of CT sinograms in a dictionary of translation-invariant wavelets and 2) the detail wavelet coefficients of a prior sinogram obtained from the forward projection of a segmented CT image. A pseudo- L0 synthesis prior was utilized to exploit and promote the sparsity of wavelet coefficients. The proposed L0-DRS MAR algorithm was compared with standard linear interpolation and the normalized metal artifact reduction (NMAR) approach proposed by Meyer using both simulated and clinical studies including hip prostheses, dental fillings, spine fixation and electroencephalogram electrodes in brain imaging. The qualitative and quantitative evaluations showed that our algorithm substantially suppresses streaking artifacts and can outperform both linear interpolation and NMAR algorithms.