Tensor decomposition and non-local means based spectral CT image denoising

BACKGROUNDAs one type of the state-of-the-art detectors, photon counting detectors are used in spectral computed tomography (CT) to classify the received photons into several energy channels and generate multichannel projections simultaneously. However, FBP reconstructed images contain severe noise due to the low photon counts in each energy channel.OBJECTIVEA spectral CT image denoising method based on tensor-decomposition and non-local means (TDNLM) is proposed.METHODSIn a CT image, it is widely accepted that there exists self-similarity over the spatial domain. In addition, because a multichannel CT image is obtained from the same object at different energies, images among different channels are highly correlated. Motivated by these two characteristics of the spectral CT images, tensor decomposition and non-local means are employed to recover fine structures in spectral CT images. Moreover, images in all energy channels are added together to form a high signal-to-noise ratio image, which is applied to encourage the signal preservation of the TDNLM. The combination of TD, NLM and the guidance of a high-quality image enhances the low-dose spectral CT, and a parameter selection strategy is designed to achieve the optimal image quality.RESULTSThe effectiveness of the developed algorithm is validated on both numerical simulations and realistic preclinical applications. The root mean square error (RMSE) and the structural similarity (SSIM) are used to quantitatively assess the image quality. The proposed method successfully restored high-quality images (average RMSE=0.0217 cm-1 and SSIM=0.987) from noisy spectral CT images (average RMSE=0.225 cm-1 and SSIM=0.633). In addition, RMSE of each decomposed material component is also remarkably reduced. Compared to the state-of-the-art iterative spectral CT reconstruction algorithms, the proposed method achieves comparable performance with dramatically reduced computational cost, resulting in a speedup of >50.CONCLUSIONSThe outstanding denoising performance, the high computational efficiency and the adaptive parameter selection strategy make the proposed method practical for spectral CT applications.