Infimal convolution-based regularization for SPECT reconstruction Journal Article


Authors: Zhang, J.; Li, S.; Krol, A.; Schmidtlein, C. R.; Lipson, E.; Feiglin, D.; Xu, Y.
Article Title: Infimal convolution-based regularization for SPECT reconstruction
Abstract: Purpose: Total variation (TV) regularization is efficient in suppressing noise, but is known to suffer from staircase artifacts. The goal of this work was to develop a regularization method using the infimal convolution of the first- and the second-order derivatives to reduce or even prevent staircase artifacts in the reconstructed images, and to investigate if the advantage in noise suppression by this TV-type regularization can be translated into dose reduction. Methods: In the present work, we introduce the infimal convolution of the first- and the second-order total variation (ICTV) as the regularization term in penalized maximum likelihood reconstruction. The preconditioned alternating projection algorithm (PAPA), previously developed by the authors of this article, was employed to produce the reconstruction. Using Monte Carlo-simulated data, we evaluate noise properties and lesion detectability in the reconstructed images and compare the results with conventional total variation (TV) and clinical EM-based methods with Gaussian post filter (GPF-EM). We also evaluate the quality of ICTV regularized images obtained for lower photon number data, compared with clinically used photon number, to verify the feasibility of radiation-dose reduction to patients by use of the ICTV reconstruction method. Results: By comparison with GPF-EM reconstructed images, we have found that the ICTV-PAPA method can achieve a lower background variability level while maintaining the same level of contrast. Images reconstructed by the ICTV-PAPA method with 80,000 counts per view exhibit even higher channelized Hotelling observer (CHO) signal-to-noise ratio (SNR), as compared to images reconstructed by the GPF-EM method with 120,000 counts per view. Conclusions: In contrast to the TV-PAPA method, the ICTV-PAPA reconstruction method avoids substantial staircase artifacts, while producing reconstructed images with higher CHO SNR and comparable local spatial resolution. Simulation studies indicate that a 33% dose reduction is feasible by switching to the ICTV-PAPA method, compared with the GPF-EM clinical standard. © 2018 American Association of Physicists in Medicine
Keywords: technetium 99m; signal noise ratio; image quality; radioactivity; artifact; image reconstruction; noise reduction; kernel method; single photon emission computed tomography; maximum likelihood method; learning algorithm; radiation dose reduction; human; article; fixed-point proximity methods; noise suppression; spect reconstruction; infimal convolution; penalized maximum likelihood optimization total variation regularization; staircase artifact
Journal Title: Medical Physics
Volume: 45
Issue: 12
ISSN: 0094-2405
Publisher: American Association of Physicists in Medicine  
Date Published: 2018-12-01
Start Page: 5397
End Page: 5410
Language: English
DOI: 10.1002/mp.13226
PROVIDER: scopus
PUBMED: 30291718
DOI/URL:
Notes: Med. Phys. -- Export Date: 2 January 2019 -- Article -- CODEN: MPHYA C2 - 30291718 -- Source: Scopus
Altmetric Score
MSK Authors