TY - JOUR
T1 - Accelerating iterative image reconstruction via adaptive surrogate functions
AU - Mitra, Ayan
AU - Politte, David G.
AU - O'Sullivan, Joseph A.
N1 - Publisher Copyright:
© 2018, Society for Imaging Science and Technology.
PY - 2018
Y1 - 2018
N2 - Three-dimensional statistical iterative reconstruction (SIR) algorithms have the potential to significantly reduce image artifacts by minimizing a cost function that models the physics and statistics of the data acquisition process in x-ray CT. SIR algorithms are important for a wide range of applications including nonstandard geometries arising from irregular sampling, limited angular range, missing data, and low-dose CT. For iterative image reconstruction algorithms to be deployed in clinical settings, the images must be quantitatively accurate and computed in clinically useful times. We describe an acceleration method that is based on adaptively varying an update factor of the additive step of the alternating minimization (AM) algorithm. Our implementation combines this method with other acceleration techniques like ordered subsets (OS) which was originally proposed for transmission tomography by Ahn, Fessler et. al [1]. Results on both an NCAT phantom and real clinical data from a Siemens Sensation 16 scanner demonstrate an improved convergence rate compared to the straightforward implementations of the alternating minimization (AM) algorithm of O'Sullivan and Benac [2] with a Huber-type edge-preserving penalty, originally proposed by Lange [3]. Our proposed acceleration method on average yields 2X acceleration of the convergence rate for both baseline and ordered subset implementations of the AM algorithm.
AB - Three-dimensional statistical iterative reconstruction (SIR) algorithms have the potential to significantly reduce image artifacts by minimizing a cost function that models the physics and statistics of the data acquisition process in x-ray CT. SIR algorithms are important for a wide range of applications including nonstandard geometries arising from irregular sampling, limited angular range, missing data, and low-dose CT. For iterative image reconstruction algorithms to be deployed in clinical settings, the images must be quantitatively accurate and computed in clinically useful times. We describe an acceleration method that is based on adaptively varying an update factor of the additive step of the alternating minimization (AM) algorithm. Our implementation combines this method with other acceleration techniques like ordered subsets (OS) which was originally proposed for transmission tomography by Ahn, Fessler et. al [1]. Results on both an NCAT phantom and real clinical data from a Siemens Sensation 16 scanner demonstrate an improved convergence rate compared to the straightforward implementations of the alternating minimization (AM) algorithm of O'Sullivan and Benac [2] with a Huber-type edge-preserving penalty, originally proposed by Lange [3]. Our proposed acceleration method on average yields 2X acceleration of the convergence rate for both baseline and ordered subset implementations of the AM algorithm.
UR - http://www.scopus.com/inward/record.url?scp=85052847782&partnerID=8YFLogxK
U2 - 10.2352/ISSN.2470-1173.2018.15.COIMG-101
DO - 10.2352/ISSN.2470-1173.2018.15.COIMG-101
M3 - Conference article
AN - SCOPUS:85052847782
SN - 2470-1173
JO - IS and T International Symposium on Electronic Imaging Science and Technology
JF - IS and T International Symposium on Electronic Imaging Science and Technology
M1 - S2
T2 - 16th Computational Imaging Conference, COMIG 2018
Y2 - 28 January 2018 through 1 February 2018
ER -