An ordered-subsets proximal preconditioned gradient algorithm for edge-preserving PET image reconstruction

Abolfazl Mehranian, Arman Rahmim, Mohammad Reza Ay, Fotis Kotasidis, Habib Zaidi

Research output: Contribution to journalArticle

Abstract

Purpose: In iterative positron emission tomography (PET) image reconstruction, the statistical variability of the PET data precorrected for random coincidences or acquired in sufficiently high count rates can be properly approximated by a Gaussian distribution, which can lead to a penalized weighted least-squares (PWLS) cost function. In this study, the authors propose a proximal preconditioned gradient algorithm accelerated with ordered subsets (PPG-OS) for the optimization of the PWLS cost function and develop a framework to incorporate boundary side information into edge-preserving total variation (TV) and Huber regularizations. Methods: The PPG-OS algorithm is proposed to address two issues encountered in the optimization of PWLS function with edge-preserving regularizers. First, the second derivative of this function (Hessian matrix) is shift-variant and ill-conditioned due to the weighting matrix (which includes emission data, attenuation, and normalization correction factors) and the regularizer. As a result, the paraboloidal surrogate functions (used in the optimization transfer techniques) end up with high curvatures and gradient-based algorithms take smaller step-sizes toward the solution, leading to a slow convergence. In addition, preconditioners used to improve the condition number of the problem, and thus to speed up the convergence, would poorly act on the resulting ill-conditioned Hessian matrix. Second, the PWLS function with a nondifferentiable penalty such as TV is not amenable to optimization using gradient-based algorithms. To deal with these issues and also to enhance edge-preservation of the TV and Huber regularizers by incorporating adaptively or anatomically derived boundary side information, the authors followed a proximal splitting method. Thereby, the optimization of the PWLS function is split into a gradient descent step (upgraded by preconditioning, step size optimization, and ordered subsets) and a proximal mapping associated with boundary weighted TV and Huber regularizers. The proximal mapping is then iteratively solved by dual formulation of the regularizers. Results: The convergence performance of the proposed algorithm was studied with three different diagonal preconditioners and compared with the state-of-the-art separable paraboloidal surrogates accelerated with ordered-subsets (SPS-OS) algorithm. In simulation studies using a realistic numerical phantom, it was shown that the proposed algorithm depicts a considerably improved convergence rate over the SPS-OS algorithm. Furthermore, the results of bias-variance and signal-to-noise evaluations showed that the proposed algorithm with anatomical edge information depicts an improved performance over conventional regularization. Finally, the proposed PPG-OS algorithm is used for image reconstruction of a clinical study with adaptively derived boundary edge information, demonstrating the potential of the algorithm for fast and edge-preserving PET image reconstruction. Conclusions: The proposed PPG-OS algorithm shows an improved convergence rate with the ability of incorporating additional boundary information in regularized PET image reconstruction.

Original languageEnglish (US)
Article number052503
JournalMedical Physics
Volume40
Issue number5
DOIs
StatePublished - May 2013

Fingerprint

Computer-Assisted Image Processing
Positron-Emission Tomography
Least-Squares Analysis
Costs and Cost Analysis
Normal Distribution
Noise

Keywords

  • image reconstruction
  • preconditioning
  • proximal gradient
  • total variation regularization

ASJC Scopus subject areas

  • Biophysics
  • Radiology Nuclear Medicine and imaging
  • Medicine(all)

Cite this

An ordered-subsets proximal preconditioned gradient algorithm for edge-preserving PET image reconstruction. / Mehranian, Abolfazl; Rahmim, Arman; Ay, Mohammad Reza; Kotasidis, Fotis; Zaidi, Habib.

In: Medical Physics, Vol. 40, No. 5, 052503, 05.2013.

Research output: Contribution to journalArticle

Mehranian, Abolfazl ; Rahmim, Arman ; Ay, Mohammad Reza ; Kotasidis, Fotis ; Zaidi, Habib. / An ordered-subsets proximal preconditioned gradient algorithm for edge-preserving PET image reconstruction. In: Medical Physics. 2013 ; Vol. 40, No. 5.
@article{be9a4f7843094b12b0ff2c81cb6d2e65,
title = "An ordered-subsets proximal preconditioned gradient algorithm for edge-preserving PET image reconstruction",
abstract = "Purpose: In iterative positron emission tomography (PET) image reconstruction, the statistical variability of the PET data precorrected for random coincidences or acquired in sufficiently high count rates can be properly approximated by a Gaussian distribution, which can lead to a penalized weighted least-squares (PWLS) cost function. In this study, the authors propose a proximal preconditioned gradient algorithm accelerated with ordered subsets (PPG-OS) for the optimization of the PWLS cost function and develop a framework to incorporate boundary side information into edge-preserving total variation (TV) and Huber regularizations. Methods: The PPG-OS algorithm is proposed to address two issues encountered in the optimization of PWLS function with edge-preserving regularizers. First, the second derivative of this function (Hessian matrix) is shift-variant and ill-conditioned due to the weighting matrix (which includes emission data, attenuation, and normalization correction factors) and the regularizer. As a result, the paraboloidal surrogate functions (used in the optimization transfer techniques) end up with high curvatures and gradient-based algorithms take smaller step-sizes toward the solution, leading to a slow convergence. In addition, preconditioners used to improve the condition number of the problem, and thus to speed up the convergence, would poorly act on the resulting ill-conditioned Hessian matrix. Second, the PWLS function with a nondifferentiable penalty such as TV is not amenable to optimization using gradient-based algorithms. To deal with these issues and also to enhance edge-preservation of the TV and Huber regularizers by incorporating adaptively or anatomically derived boundary side information, the authors followed a proximal splitting method. Thereby, the optimization of the PWLS function is split into a gradient descent step (upgraded by preconditioning, step size optimization, and ordered subsets) and a proximal mapping associated with boundary weighted TV and Huber regularizers. The proximal mapping is then iteratively solved by dual formulation of the regularizers. Results: The convergence performance of the proposed algorithm was studied with three different diagonal preconditioners and compared with the state-of-the-art separable paraboloidal surrogates accelerated with ordered-subsets (SPS-OS) algorithm. In simulation studies using a realistic numerical phantom, it was shown that the proposed algorithm depicts a considerably improved convergence rate over the SPS-OS algorithm. Furthermore, the results of bias-variance and signal-to-noise evaluations showed that the proposed algorithm with anatomical edge information depicts an improved performance over conventional regularization. Finally, the proposed PPG-OS algorithm is used for image reconstruction of a clinical study with adaptively derived boundary edge information, demonstrating the potential of the algorithm for fast and edge-preserving PET image reconstruction. Conclusions: The proposed PPG-OS algorithm shows an improved convergence rate with the ability of incorporating additional boundary information in regularized PET image reconstruction.",
keywords = "image reconstruction, preconditioning, proximal gradient, total variation regularization",
author = "Abolfazl Mehranian and Arman Rahmim and Ay, {Mohammad Reza} and Fotis Kotasidis and Habib Zaidi",
year = "2013",
month = "5",
doi = "10.1118/1.4801898",
language = "English (US)",
volume = "40",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "5",

}

TY - JOUR

T1 - An ordered-subsets proximal preconditioned gradient algorithm for edge-preserving PET image reconstruction

AU - Mehranian, Abolfazl

AU - Rahmim, Arman

AU - Ay, Mohammad Reza

AU - Kotasidis, Fotis

AU - Zaidi, Habib

PY - 2013/5

Y1 - 2013/5

N2 - Purpose: In iterative positron emission tomography (PET) image reconstruction, the statistical variability of the PET data precorrected for random coincidences or acquired in sufficiently high count rates can be properly approximated by a Gaussian distribution, which can lead to a penalized weighted least-squares (PWLS) cost function. In this study, the authors propose a proximal preconditioned gradient algorithm accelerated with ordered subsets (PPG-OS) for the optimization of the PWLS cost function and develop a framework to incorporate boundary side information into edge-preserving total variation (TV) and Huber regularizations. Methods: The PPG-OS algorithm is proposed to address two issues encountered in the optimization of PWLS function with edge-preserving regularizers. First, the second derivative of this function (Hessian matrix) is shift-variant and ill-conditioned due to the weighting matrix (which includes emission data, attenuation, and normalization correction factors) and the regularizer. As a result, the paraboloidal surrogate functions (used in the optimization transfer techniques) end up with high curvatures and gradient-based algorithms take smaller step-sizes toward the solution, leading to a slow convergence. In addition, preconditioners used to improve the condition number of the problem, and thus to speed up the convergence, would poorly act on the resulting ill-conditioned Hessian matrix. Second, the PWLS function with a nondifferentiable penalty such as TV is not amenable to optimization using gradient-based algorithms. To deal with these issues and also to enhance edge-preservation of the TV and Huber regularizers by incorporating adaptively or anatomically derived boundary side information, the authors followed a proximal splitting method. Thereby, the optimization of the PWLS function is split into a gradient descent step (upgraded by preconditioning, step size optimization, and ordered subsets) and a proximal mapping associated with boundary weighted TV and Huber regularizers. The proximal mapping is then iteratively solved by dual formulation of the regularizers. Results: The convergence performance of the proposed algorithm was studied with three different diagonal preconditioners and compared with the state-of-the-art separable paraboloidal surrogates accelerated with ordered-subsets (SPS-OS) algorithm. In simulation studies using a realistic numerical phantom, it was shown that the proposed algorithm depicts a considerably improved convergence rate over the SPS-OS algorithm. Furthermore, the results of bias-variance and signal-to-noise evaluations showed that the proposed algorithm with anatomical edge information depicts an improved performance over conventional regularization. Finally, the proposed PPG-OS algorithm is used for image reconstruction of a clinical study with adaptively derived boundary edge information, demonstrating the potential of the algorithm for fast and edge-preserving PET image reconstruction. Conclusions: The proposed PPG-OS algorithm shows an improved convergence rate with the ability of incorporating additional boundary information in regularized PET image reconstruction.

AB - Purpose: In iterative positron emission tomography (PET) image reconstruction, the statistical variability of the PET data precorrected for random coincidences or acquired in sufficiently high count rates can be properly approximated by a Gaussian distribution, which can lead to a penalized weighted least-squares (PWLS) cost function. In this study, the authors propose a proximal preconditioned gradient algorithm accelerated with ordered subsets (PPG-OS) for the optimization of the PWLS cost function and develop a framework to incorporate boundary side information into edge-preserving total variation (TV) and Huber regularizations. Methods: The PPG-OS algorithm is proposed to address two issues encountered in the optimization of PWLS function with edge-preserving regularizers. First, the second derivative of this function (Hessian matrix) is shift-variant and ill-conditioned due to the weighting matrix (which includes emission data, attenuation, and normalization correction factors) and the regularizer. As a result, the paraboloidal surrogate functions (used in the optimization transfer techniques) end up with high curvatures and gradient-based algorithms take smaller step-sizes toward the solution, leading to a slow convergence. In addition, preconditioners used to improve the condition number of the problem, and thus to speed up the convergence, would poorly act on the resulting ill-conditioned Hessian matrix. Second, the PWLS function with a nondifferentiable penalty such as TV is not amenable to optimization using gradient-based algorithms. To deal with these issues and also to enhance edge-preservation of the TV and Huber regularizers by incorporating adaptively or anatomically derived boundary side information, the authors followed a proximal splitting method. Thereby, the optimization of the PWLS function is split into a gradient descent step (upgraded by preconditioning, step size optimization, and ordered subsets) and a proximal mapping associated with boundary weighted TV and Huber regularizers. The proximal mapping is then iteratively solved by dual formulation of the regularizers. Results: The convergence performance of the proposed algorithm was studied with three different diagonal preconditioners and compared with the state-of-the-art separable paraboloidal surrogates accelerated with ordered-subsets (SPS-OS) algorithm. In simulation studies using a realistic numerical phantom, it was shown that the proposed algorithm depicts a considerably improved convergence rate over the SPS-OS algorithm. Furthermore, the results of bias-variance and signal-to-noise evaluations showed that the proposed algorithm with anatomical edge information depicts an improved performance over conventional regularization. Finally, the proposed PPG-OS algorithm is used for image reconstruction of a clinical study with adaptively derived boundary edge information, demonstrating the potential of the algorithm for fast and edge-preserving PET image reconstruction. Conclusions: The proposed PPG-OS algorithm shows an improved convergence rate with the ability of incorporating additional boundary information in regularized PET image reconstruction.

KW - image reconstruction

KW - preconditioning

KW - proximal gradient

KW - total variation regularization

UR - http://www.scopus.com/inward/record.url?scp=84877301372&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84877301372&partnerID=8YFLogxK

U2 - 10.1118/1.4801898

DO - 10.1118/1.4801898

M3 - Article

C2 - 23635293

AN - SCOPUS:84877301372

VL - 40

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

IS - 5

M1 - 052503

ER -