Abstract:In this paper, a total variation (TV) minimization strategy is proposed to overcome the problem of sparse spatial resolution and large amounts of noise in low dose positron emission tomography (PET) imaging reconstruction. Two types of objective function were established based on two statistical models of measured PET data, least-square (LS) TV for the Gaussian distribution and Poisson-TV for the Poisson distribution. To efficiently obtain high quality reconstructed images, the alternating direction method (ADM) is used to solve these objective functions. As compared with the iterative shrinkage/thresholding (IST) based algorithms, the proposed ADM can make full use of the TV constraint and its convergence rate is faster. The performance of the proposed approach is validated through comparisons with the expectation-maximization (EM) method using synthetic and experimental biological data. In the comparisons, the results of both LS-TV and Poisson-TV are taken into consideration to find which models are more suitable for PET imaging, in particular low-dose PET. To evaluate the results quantitatively, we computed bias, variance, and the contrast recovery coefficient (CRC) and drew profiles of the reconstructed images produced by the different methods. The results show that both Poisson-TV and LS-TV can provide a high visual quality at a low dose level. The bias and variance of the proposed LS-TV and Poisson-TV methods are 20% to 74% less at all counting levels than those of the EM method. Poisson-TV gives the best performance in terms of high-accuracy reconstruction with the lowest bias and variance as compared to the ground truth (14.3% less bias and 21.9% less variance). In contrast, LS-TV gives the best performance in terms of the high contrast of the reconstruction with the highest CRC.

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

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

Alternating direction method of multipliers for nonconvex log total variation image restoration

A stochastic ADMM algorithm for large-scale ptychography with weighted difference of anisotropic and isotropic total variation

Weighted Anisotropic-Isotropic Total Variation for Poisson Denoising

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

An Inexact Majorized Proximal Alternating Direction Method of Multipliers for Diffusion Tensors

Alternating Minimization Algorithm with Automatic Relevance Determination for Transmission Tomography under Poisson Noise

An Inexact Majorized Proximal Alternating Direction Method of Multipliers for Diffusion Tensors\ast

Hybrid Priors Based on Weighted Hyper-Laplacian with Overlapping Group Sparsity for Poisson Noise Removal

Dual Alternating Direction Method of Multipliers for Inverse Imaging

Efficient Multiplicative Noise Removal Method Using Isotropic Second Order Total Variation.

An alternating direction method of multipliers for inverse lithography problem

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

Anisotropic total generalized variation model for Poisson noise removal

Alternating Direction Method of Multipliers for Nonlinear Image Restoration Problems

Alternating Direction Method of Multipliers for a Class of Nonconvex Bilinear Optimization: Convergence Analysis and Applications.

High-order Total Variation-Based Multiplicative Noise Removal with Spatially Adapted Parameter Selection.

A Convex Total Generalized Variation Regularized Model for Multiplicative Noise and Blur Removal

A HYBRID ALTERNATING MINIMIZATION ALGORITHM FOR STRUCTURED CONVEX OPTIMIZATION PROBLEMS WITH APPLICATION IN POISSONIAN IMAGE PROCESSING