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

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.

Improved discrimination between benign and malignant LDCT screening-detected lung nodules with dynamic over static 18F-FDG PET as a function of injected dose

“…For low-dose dynamic PET imaging with simplified protocols, additional improvement of imaging methods requires further investigation, such as noise reduction, motion correction and direct parametric reconstruction. Smoothing can be used for low-dose data, which may decrease the noise level and help to obtain results closer to 100% counts (Yu et al , 2016a). The amount of body motion, which induces image blurring and underestimation of tracer uptake, in a short scan duration can be less as compared to a longer scan.…”

Delayed FDG PET Provides Superior Glioblastoma Conspicuity Compared to Conventional Image Timing

“…The table time of the 5-h imaging session in this study was 20 min, which provided adequate count statistics and subjective quality assessment. Future work may identify that shorter imaging times are sufficient, particularly given work on novel reconstruction algorithms for PET imaging with lower radiopharmaceutical doses ( 24 , 25 ).…”

Unsupervised PET Reconstruction from a Bayesian Perspective

“…Since several additive sources of error tend to form a Gaussian distribution, a Gaussian noise model is suitable for PET data. In fact, [51] shows that both Poisson and Gaussian models can be leveraged for PET reconstruction. In this article, to introduce the DeepRED framework, we adopt a Gaussian-based loglikelihood function as:…”