LAUSR

Multi-material decomposition (MMD) technique decomposes the CT images into basis material images and has been promising in clinical practice for material composition quantification within the human body. MMD could be implemented using the image data acquired from spectral CT or its special case, dual-energy CT (DECT) while the spectral CT data acquisition usually requires a hardware modification. In this paper, we propose an image domain MMD method using single energy CT (SECT). The proposed objective function applies a least square data fidelity term to enforce the minimization between the linear combination of decomposed material image and the measured SECT image, and an edge-preserving (EP) regularization term to meet the piecewise constant property of the material image. We apply the optimization transfer principle to form a pixel-wise separable quadratic surrogate (PWSQS) function in each iteration to decrease the objective function. The pixelwise direct inversion method assisted by the two-material assumption (TMA) is applied to obtain a good initial value. The proposed method is evaluated using a digital phantom, a Catphan phantom and the clinical data. A low-pass filtration method is implemented for a comparison purpose. In the phantom study, the proposed TMA method achieves high volume fraction accuracy (VFA) of 79.64% and the proposed EP method further increases the VFA by 15.56% and decreases the decomposition standard deviation (STD) by 81.51% compared with the TMA method. At the comparable noise level, the proposed EP method increases spatial resolution by an overall factor of 1.01 when the modulation transfer function magnitude is decreased to 50% compared with the low-pass filtration method. In clinical data study, the virtual non-contrast image generated by the proposed method achieves the root-mean-squared-relative error of 2.93% compared with the contrast-free ground-truth image.