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.

Sparse Signal Reconstruction for Overdispersed Low-photon Count Biomedical Imaging Using $\ell_p$ Total Variation

Alternating Direction Method of Multipliers for Negative Binomial Model with The Weighted Difference of Anisotropic and Isotropic Total Variation

Low dose PET reconstruction with total variation regularization.

Low Dose PET Image Reconstruction with Total Variation Using Alternating Direction Method.

Nonlocal Low-Rank and Total Variation Constrained PET Image Reconstruction

High-order Total Variation-Based Poissonian Image Deconvolution with Spatially Adapted Regularization Parameter

Off-the-grid regularisation for Poisson inverse problems

Nonlocal Total Variation Based Dynamic PET Image Reconstruction with Low-Rank Constraints

Sparse representation and dictionary learning penalized image reconstruction for positron emission tomography.

Hybrid High-Order Nonlocal Gradient Sparsity Regularization for Poisson Image Deconvolution.

Constrained Higher Degree Total P-Variation Minimization for MRI Reconstruction from Undersampled K-Space Data

SPOQ $\ell_p$-Over-$\ell_q$ Regularization for Sparse Signal Recovery applied to Mass Spectrometry

Sparse Recovery: The Square of $\ell_1/\ell_2$ Norms

Recovering Sparse Nonnegative Signals via Non-convex Fraction Function Penalty

Sparsity promoting regularization for effective noise suppression in SPECT image reconstruction

Constrained Tpv Minimization for Enhanced Exploitation of Gradient Sparsity: Application to Ct Image Reconstruction

Image reconstructions using sparse dictionary representations and implicit, non-negative mappings

Robust single-photon counting imaging with spatially correlated and total variation constraints

Nonconvex Regularization Based Sparse Recovery and Demixing with Application to Color Image Inpainting

Iterative image reconstruction for sparse-view CT via total variation regularization and dictionary learning.

Sparse Non-Negative Signal Reconstruction Using Fraction Function Penalty