Abstract:Total-variation (TV)-based Computed Tomography (CT) image reconstruction hasshown experimentally to be capable of producing accurate reconstructions fromsparse-view data. In particular TV-based reconstruction is very well suited forimages with piecewise nearly constant regions. Computationally, however,TV-based reconstruction is much more demanding, especially for 3D imaging, andthe reconstruction from clinical data sets is far from being close toreal-time. This is undesirable from a clinical perspective, and thus there isan incentive to accelerate the solution of the underlying optimization problem.The TV reconstruction can in principle be found by any optimization method, butin practice the large-scale systems arising in CT image reconstruction precludethe use of memory-demanding methods such as Newton's method. The simplegradient method has much lower memory requirements, but exhibits slowconvergence. In the present work we consider the use of two acceleratedgradient-based methods, GPBB and UPN, for reducing the number of gradientmethod iterations needed to achieve a high-accuracy TV solution in CT imagereconstruction. The former incorporates several heuristics from theoptimization literature such as Barzilai-Borwein (BB) step size selection andnonmonotone line search. The latter uses a cleverly chosen sequence ofauxiliary points to achieve a better convergence rate. The methods are memoryefficient and equipped with a stopping criterion to ensure that the TVreconstruction has indeed been found. An implementation of the methods (in Cwith interface to Matlab) is available for download fromhttp://www2.imm.dtu.dk/ pch/TVReg/. We compare the proposed methods with thestandard gradient method, applied to a 3D test problem with synthetic few-viewdata. We find experimentally that for realistic parameters the proposed methodssignificantly outperform the gradient method.

Perturbed conjugate gradient method: controlling semi-convergence, enforcing non-negativity, and accelerating iterations in tomographic imaging

New Iterative Gridding Algorithm Using Conjugate Gradient Method

A New Spectral Conjugate Subgradient Method with Application in Computed Tomography Image Reconstruction

Conjugate-gradient preconditioning methods for shift-variant PET image reconstruction

Regularized Total Least Squares Reconstruction for Optical Tomographic Imaging Using Conjugate Gradient Method

Modification of a conjugate gradient approach for convex constrained nonlinear monotone equations with applications in signal recovery and image restoration

Gradient Descent Provably Solves Nonlinear Tomographic Reconstruction

A three-term conjugate gradient descent method with some applications

Tensor Conjugate Gradient Methods with Automatically Determination of Regularization Parameters for Ill-Posed Problems with t-Product

Iterative Total Least-Squares Image Reconstruction Algorithm for Optical Tomography by the Conjugate Gradient Method

The conjugate gradient method with various viewpoints

A Penalized Linear and Nonlinear Combined Conjugate Gradient Method for the Reconstruction of Fluorescence Molecular Tomography

Regularized Step Directions in Nonlinear Conjugate Gradient Methods

Convergence Analysis of Nonlinear Kaczmarz Method for Systems of Nonlinear Equations with Component-wise Convex Mapping

Nonlinear tomographic reconstruction via nonsmooth optimization

Accelerated Gradient Methods for Total-Variation-based CT Image Reconstruction

Image Noise Reduction and Solution of Unconstrained Minimization Problems via New Conjugate Gradient Methods

An efficient inertial subspace minimization CG algorithm with convergence rate analysis for constrained nonlinear monotone equations

Barzilai and Borwein conjugate gradient method equipped with a non-monotone line search technique and its application on non-negative matrix factorization

Iterative Nonlinear Tikhonov Algorithm with Constraints for Electromagnetic Tomography