Conference paper
A First-Order Primal-Dual Reconstruction Algorithm for Few-View SPECT
A sparsity-exploiting algorithm intended for few-view Single Photon Emission Computed Tomography (SPECT) reconstruction is proposed and characterized. The algorithm models the object as piecewise constant subject to a blurring operation. Monte Carlo simulations were performed to provide more projection data of a phantom with varying smoothness across the field of view.
For all simulations, reconstructions were performed across a sweep of the two primary design parameters: the blurring parameter and the weighting of the total variation (TV) minimization term. Maximum-Likelihood Expectation Maximization (MLEM) reconstructions were performed to provide reference images.
Spatial resolution, accuracy, and signal-to-noise ratio was calculated and compared for all reconstructions. In general, increased values of the blurring parameter and TV weighting parameters reduced noise and streaking artifacts, while decreasing spatial resolution. The reconstructed images demonstrate that the algorithm introduces low-frequency artifacts in some cases, but eliminates streak artifacts due to angular undersampling.
Further, as the number of views was decreased from 60 to 9 the accuracy of images reconstructed using the proposed algorithm varied by less than 3%. Overall, the results demonstrate preliminary feasibility of a sparsity-exploiting reconstruction algorithm which may be beneficial for few-view SPECT.
Language: | English |
---|---|
Publisher: | IEEE |
Year: | 2012 |
Pages: | 2381-2385 |
Proceedings: | 2012 IEEE Nuclear Science Symposium and Medical Imaging Conference |
ISBN: | 1467320285 , 1467320293 , 1467320307 , 9781467320283 , 9781467320290 and 9781467320306 |
ISSN: | 25770829 and 10823654 |
Types: | Conference paper |
DOI: | 10.1109/NSSMIC.2012.6551542 |
ORCIDs: | Jørgensen, Jakob Heide |
MLEM reconstruction Monte Carlo methods Monte Carlo simulation TV weighting parameter angular undersampling blurring operation blurring parameter expectation-maximisation algorithm few-view Single Photon Emission Computed Tomography reconstruction field of view first-order primal-dual reconstruction algorithm image denoising image reconstruction image resolution low-frequency artifact maximum-likelihood expectation maximization reconstruction medical image processing minimisation noise reduction phantom phantoms piecewise constant subject projection data reference image signal-to-noise ratio single photon emission computed tomography sparsity-exploiting reconstruction algorithm spatial accuracy spatial resolution streak artifact elimination streaking artifact total variation minimization term