Convergence criterion for MBIR based on the local noise-power spectrum: Theory and implementation in a framework for accelerated 3D image reconstruction with a morphological pyramid

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Abstract

Model-based iterative reconstruction (MBIR) offers improved noise-resolution tradeoffs and artifact reduction in conebeam CT compared to analytical reconstruction, but carries increased computational burden. An important consideration in minimizing computation time is reliable selection of the stopping criterion to perform the minimum number of iterations required to obtain the desired image quality. Most MBIR methods rely on a fixed number of iterations or relative metrics on image or cost-function evolution, and it would be desirable to use metrics that are more representative of the underlying image properties. A second front for reduction of computation time is the use of acceleration techniques (e.g. subsets or momentum). However, most of these techniques do not strictly guarantee convergence of the resulting MBIR method. A data-dependent analytical model of noise-power spectrum (NPS) for penalized weighted least squares (PWLS) reconstruction is proposed as an absolute metric of image properties for the fully converged volume. Distance to convergence is estimated as the root mean squared error (RMSE) between the estimated NPS and an NPS measured on a uniform region of interest (ROI) in the evolving volume. Iterations are stopped when the RMSE falls below a threshold directly related with the properties of the target image. Further acceleration was achieved by combining the spectral stopping criterion with a morphological pyramid (mPyr) in which the minimization of the PWLS cost-function is divided in a cascade of stages. The algorithm parameters (voxel size in this work) change between stages to achieve faster evolution in early stages, and a final stage with the target parameters to guarantee convergence. Transition between stages is governed by the spectral stopping criterion.

Original languageEnglish (US)
Title of host publication15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine
EditorsSamuel Matej, Scott D. Metzler
PublisherSPIE
ISBN (Electronic)9781510628373
DOIs
StatePublished - Jan 1 2019
Event15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine, Fully3D 2019 - Philadelphia, United States
Duration: Jun 2 2019Jun 6 2019

Publication series

NameProceedings of SPIE - The International Society for Optical Engineering
Volume11072
ISSN (Print)0277-786X
ISSN (Electronic)1996-756X

Conference

Conference15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine, Fully3D 2019
CountryUnited States
CityPhiladelphia
Period6/2/196/6/19

Fingerprint

Convergence Criteria
3D Reconstruction
Pyramid
noise spectra
3D Image
Image Reconstruction
image reconstruction
Power spectrum
pyramids
Image reconstruction
Power Spectrum
power spectra
Stopping Criterion
Model-based
Cost functions
stopping
iteration
Penalized Least Squares
Weighted Least Squares
Iteration

ASJC Scopus subject areas

  • Electronic, Optical and Magnetic Materials
  • Condensed Matter Physics
  • Computer Science Applications
  • Applied Mathematics
  • Electrical and Electronic Engineering

Cite this

Sisniega, A., Stayman, J. W., Capostagno, S., Weiss, C. R., Ehtiati, T., & Siewerdsen, J. H. (2019). Convergence criterion for MBIR based on the local noise-power spectrum: Theory and implementation in a framework for accelerated 3D image reconstruction with a morphological pyramid. In S. Matej, & S. D. Metzler (Eds.), 15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine [1107209] (Proceedings of SPIE - The International Society for Optical Engineering; Vol. 11072). SPIE. https://doi.org/10.1117/12.2534896

Convergence criterion for MBIR based on the local noise-power spectrum : Theory and implementation in a framework for accelerated 3D image reconstruction with a morphological pyramid. / Sisniega, A.; Stayman, J. W.; Capostagno, S.; Weiss, C. R.; Ehtiati, T.; Siewerdsen, J. H.

15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine. ed. / Samuel Matej; Scott D. Metzler. SPIE, 2019. 1107209 (Proceedings of SPIE - The International Society for Optical Engineering; Vol. 11072).

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Sisniega, A, Stayman, JW, Capostagno, S, Weiss, CR, Ehtiati, T & Siewerdsen, JH 2019, Convergence criterion for MBIR based on the local noise-power spectrum: Theory and implementation in a framework for accelerated 3D image reconstruction with a morphological pyramid. in S Matej & SD Metzler (eds), 15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine., 1107209, Proceedings of SPIE - The International Society for Optical Engineering, vol. 11072, SPIE, 15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine, Fully3D 2019, Philadelphia, United States, 6/2/19. https://doi.org/10.1117/12.2534896
Sisniega A, Stayman JW, Capostagno S, Weiss CR, Ehtiati T, Siewerdsen JH. Convergence criterion for MBIR based on the local noise-power spectrum: Theory and implementation in a framework for accelerated 3D image reconstruction with a morphological pyramid. In Matej S, Metzler SD, editors, 15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine. SPIE. 2019. 1107209. (Proceedings of SPIE - The International Society for Optical Engineering). https://doi.org/10.1117/12.2534896
Sisniega, A. ; Stayman, J. W. ; Capostagno, S. ; Weiss, C. R. ; Ehtiati, T. ; Siewerdsen, J. H. / Convergence criterion for MBIR based on the local noise-power spectrum : Theory and implementation in a framework for accelerated 3D image reconstruction with a morphological pyramid. 15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine. editor / Samuel Matej ; Scott D. Metzler. SPIE, 2019. (Proceedings of SPIE - The International Society for Optical Engineering).
@inproceedings{008530c2824444d6b117720be26143c9,
title = "Convergence criterion for MBIR based on the local noise-power spectrum: Theory and implementation in a framework for accelerated 3D image reconstruction with a morphological pyramid",
abstract = "Model-based iterative reconstruction (MBIR) offers improved noise-resolution tradeoffs and artifact reduction in conebeam CT compared to analytical reconstruction, but carries increased computational burden. An important consideration in minimizing computation time is reliable selection of the stopping criterion to perform the minimum number of iterations required to obtain the desired image quality. Most MBIR methods rely on a fixed number of iterations or relative metrics on image or cost-function evolution, and it would be desirable to use metrics that are more representative of the underlying image properties. A second front for reduction of computation time is the use of acceleration techniques (e.g. subsets or momentum). However, most of these techniques do not strictly guarantee convergence of the resulting MBIR method. A data-dependent analytical model of noise-power spectrum (NPS) for penalized weighted least squares (PWLS) reconstruction is proposed as an absolute metric of image properties for the fully converged volume. Distance to convergence is estimated as the root mean squared error (RMSE) between the estimated NPS and an NPS measured on a uniform region of interest (ROI) in the evolving volume. Iterations are stopped when the RMSE falls below a threshold directly related with the properties of the target image. Further acceleration was achieved by combining the spectral stopping criterion with a morphological pyramid (mPyr) in which the minimization of the PWLS cost-function is divided in a cascade of stages. The algorithm parameters (voxel size in this work) change between stages to achieve faster evolution in early stages, and a final stage with the target parameters to guarantee convergence. Transition between stages is governed by the spectral stopping criterion.",
author = "A. Sisniega and Stayman, {J. W.} and S. Capostagno and Weiss, {C. R.} and T. Ehtiati and Siewerdsen, {J. H.}",
year = "2019",
month = "1",
day = "1",
doi = "10.1117/12.2534896",
language = "English (US)",
series = "Proceedings of SPIE - The International Society for Optical Engineering",
publisher = "SPIE",
editor = "Samuel Matej and Metzler, {Scott D.}",
booktitle = "15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine",

}

TY - GEN

T1 - Convergence criterion for MBIR based on the local noise-power spectrum

T2 - Theory and implementation in a framework for accelerated 3D image reconstruction with a morphological pyramid

AU - Sisniega, A.

AU - Stayman, J. W.

AU - Capostagno, S.

AU - Weiss, C. R.

AU - Ehtiati, T.

AU - Siewerdsen, J. H.

PY - 2019/1/1

Y1 - 2019/1/1

N2 - Model-based iterative reconstruction (MBIR) offers improved noise-resolution tradeoffs and artifact reduction in conebeam CT compared to analytical reconstruction, but carries increased computational burden. An important consideration in minimizing computation time is reliable selection of the stopping criterion to perform the minimum number of iterations required to obtain the desired image quality. Most MBIR methods rely on a fixed number of iterations or relative metrics on image or cost-function evolution, and it would be desirable to use metrics that are more representative of the underlying image properties. A second front for reduction of computation time is the use of acceleration techniques (e.g. subsets or momentum). However, most of these techniques do not strictly guarantee convergence of the resulting MBIR method. A data-dependent analytical model of noise-power spectrum (NPS) for penalized weighted least squares (PWLS) reconstruction is proposed as an absolute metric of image properties for the fully converged volume. Distance to convergence is estimated as the root mean squared error (RMSE) between the estimated NPS and an NPS measured on a uniform region of interest (ROI) in the evolving volume. Iterations are stopped when the RMSE falls below a threshold directly related with the properties of the target image. Further acceleration was achieved by combining the spectral stopping criterion with a morphological pyramid (mPyr) in which the minimization of the PWLS cost-function is divided in a cascade of stages. The algorithm parameters (voxel size in this work) change between stages to achieve faster evolution in early stages, and a final stage with the target parameters to guarantee convergence. Transition between stages is governed by the spectral stopping criterion.

AB - Model-based iterative reconstruction (MBIR) offers improved noise-resolution tradeoffs and artifact reduction in conebeam CT compared to analytical reconstruction, but carries increased computational burden. An important consideration in minimizing computation time is reliable selection of the stopping criterion to perform the minimum number of iterations required to obtain the desired image quality. Most MBIR methods rely on a fixed number of iterations or relative metrics on image or cost-function evolution, and it would be desirable to use metrics that are more representative of the underlying image properties. A second front for reduction of computation time is the use of acceleration techniques (e.g. subsets or momentum). However, most of these techniques do not strictly guarantee convergence of the resulting MBIR method. A data-dependent analytical model of noise-power spectrum (NPS) for penalized weighted least squares (PWLS) reconstruction is proposed as an absolute metric of image properties for the fully converged volume. Distance to convergence is estimated as the root mean squared error (RMSE) between the estimated NPS and an NPS measured on a uniform region of interest (ROI) in the evolving volume. Iterations are stopped when the RMSE falls below a threshold directly related with the properties of the target image. Further acceleration was achieved by combining the spectral stopping criterion with a morphological pyramid (mPyr) in which the minimization of the PWLS cost-function is divided in a cascade of stages. The algorithm parameters (voxel size in this work) change between stages to achieve faster evolution in early stages, and a final stage with the target parameters to guarantee convergence. Transition between stages is governed by the spectral stopping criterion.

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

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

U2 - 10.1117/12.2534896

DO - 10.1117/12.2534896

M3 - Conference contribution

AN - SCOPUS:85074295273

T3 - Proceedings of SPIE - The International Society for Optical Engineering

BT - 15th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine

A2 - Matej, Samuel

A2 - Metzler, Scott D.

PB - SPIE

ER -