Abstract:The purpose of this research is to develop an advanced reconstruction method for low-count, hence high-noise, single-photon emission computed tomography (SPECT) image reconstruction. It consists of a novel reconstruction model to suppress noise while conducting reconstruction and an efficient algorithm to solve the model. A novel regularizer is introduced as the nonconvex denoising term based on the approximate sparsity of the image under a geometric tight frame transform domain. The deblurring term is based on the negative log-likelihood of the SPECT data model. To solve the resulting nonconvex optimization problem a preconditioned fixed-point proximity algorithm (PFPA) is introduced. We prove that under appropriate assumptions, PFPA converges to a local solution of the optimization problem at a global O ( 1 / k ) convergence rate. Substantial numerical results for simulation data are presented to demonstrate the superiority of the proposed method in denoising, suppressing artifacts and reconstruction accuracy. We simulate noisy 2D SPECT data from two phantoms: hot Gaussian spheres on random lumpy warm background, and the anthropomorphic brain phantom, at high- and low-noise levels (64k and 90k counts, respectively), and reconstruct them using PFPA. We also perform limited comparative studies with selected competing state-of-the-art total variation (TV) and higher-order TV (HOTV) transform-based methods, and widely used post-filtered maximum-likelihood expectation-maximization. We investigate imaging performance of these methods using: contrast-to-noise ratio (CNR), ensemble variance images (EVI), background ensemble noise (BEN), normalized mean-square error (NMSE), and channelized hotelling observer (CHO) detectability. Each of the competing methods is independently optimized for each metric. We establish that the proposed method outperforms the other approaches in all image quality metrics except NMSE where it is matched by HOTV. The superiority of the proposed method is especially evident in the CHO detectability tests results. We also perform qualitative image evaluation for presence and severity of image artifacts where it also performs better in terms of suppressing 'staircase' artifacts, as compared to TV methods. However, edge artifacts on high-contrast regions persist. We conclude that the proposed method may offer a powerful tool for detection tasks in high-noise SPECT imaging.

Combined Transformation of Ordering SPECT Sinograms for Signal Extraction from Measurements with Poisson Noise.

A Noise Reduction Method for Non-Stationary Noise Model of SPECT Sinogram Based on Kalman Filter

3D Distance-Weighted Wiener Filter for Poisson Noise Reduction in Sinogram Space for SPECT Imaging

Nonlinear Sinogram Smoothing for Low-Dose X-ray CT

Statistical Sinogram Restoration For Single Photon Emission Computed Tomography

[Analytical Algorithms Improving the Reconstruction of SPECT Image].

A Theoretically Based Pre-Reconstructing Filter for Spatio-Temporal Noise Reduction in Gated Cardiac SPECT

Noise Reduction of Low-Dose CT Sinograms Based on EM Algorithm

Analytical Noise Treatment for Low-Dose Ct Projection Data by Penalized Weighted Least-Square Smoothing in the K-L Domain

Sparsity promoting regularization for effective noise suppression in SPECT image reconstruction

Noise reduction for cone-beam SPECT by penalized reweighted least-squares projection restoration

Penalized-Likelihood Sinogram Restoration for Computed Tomography

Nonlinear Dual Reconstruction of SPECT Activity and Attenuation Images

Potential Equivalence of Sinogram and Image-Domain Penalized Likelihood Methods

Analytical reconstruction of 4-D dynamic cardiac SPECT with noise-reduction and non-uniform attenuation compensation

The De-noising Methods for SPECT Reconstruction with Non-uniform Attenuation Based on Wavelet Transform

Simultaneous reconstruction of tissue attenuation and radioactivity maps in SPECT.

Wavelet based de-noising methods for local SPECT reconstruction with non-uniform attenuation

A Unified Analysis of Exact Methods of Inverting the 2-D Exponential Radon Transform, with Implications for Noise Control in SPECT

Robust Low-Dose CT Sinogram Preprocessing Via Exploiting Noise-Generating Mechanism.

Non-iterative Methods and Their Noise Characteristics in 2D SPECT Image Reconstruction