Iterative spectral CT reconstruction based on low rank and average-image-incorporated BM3D.

In a photon counting detector integrated spectral CT scanner, the received photons are counted in several energy channels to generate the corresponding projections. Since the projection in each energy channel is generated using part of the received photons, the reconstructed channel image suffers from severe noise. Therefore, image reconstruction in spectral CT is considered to be a big challenge. Because the inter-channel images are all from the same object but in different energy bins, there exists a strong correlation among these images. Moreover, it is suggested that there are similarities among various patches of CT images in the spatial domain. In this work, we propose average-image-incorporated block-matching and 3D (aiiBM3D) filtering along with low rank regularization for iterative spectral CT reconstruction. The aiiBM3D method is based on filtered 3D data arrays formed by similar 2D blocks using the mapped version of the average image obtained from linear regression. The reconstruction procedure consists of two main steps. First, the alternating direction method of multipliers is employed to solve the problem with low rank regularization where the goal is to exploit the correlation in inter-channel images. Second, our proposed BM3D-based algorithm is applied to all the channel images to make use of the redundant information in the spatial domain and inter-channel. The two steps repeat until the stopping criteria are satisfied. The proposed method is validated on numerically simulated and preclinical datasets. Our results confirm its high performance in terms of signal to noise ratio and structural preservation.