Investigation of the preconditioner-parameter in the preconditioned Chambolle-Pock algorithm applied to optimization-based image reconstruction.

The optimization-based image reconstruction methods have been thoroughly investigated in the field of medical imaging. The Chambolle-Pock (CP) algorithm may be employed to solve these convex optimization image reconstruction programs. The preconditioned CP (PCP) algorithm has been shown to have much higher convergence rate than the ordinary CP (OCP) algorithm. This algorithm utilizes a preconditioner-parameter to tune the implementation of the algorithm to the specific application, which ranges from 0 and 2, but is often set to 1. In this work, we investigated the impact of the preconditioner-parameter on the convergence rate of the PCP algorithm when it is applied to the TV constrained, data-divergence minimization (TVDM) optimization based image reconstruction. We performed the investigations in the context of 2D computed tomography (CT) and 3D electron paramagnetic resonance imaging (EPRI). For 2D CT, we used the Shepp-Logan and two FORBILD phantoms. For 3D EPRI, we used a simulated 6-spheres phantom and a physical phantom. Study results showed that the optimal preconditioner-parameter depends on the specific imaging conditions. Simply setting the parameter equal to 1 cannot guarantee a fast convergence rate. Thus, this study suggests that one should adaptively tune the preconditioner-parameter to obtain the optimal convergence rate of the PCP algorithm.

[1]  Howard J Halpern,et al.  Maximally spaced projection sequencing in electron paramagnetic resonance imaging. , 2015, Concepts in magnetic resonance. Part B, Magnetic resonance engineering.

[2]  R. Siddon.  Fast calculation of the exact radiological path for a three-dimensional CT array. , 1985, Medical physics.

[3]  Xiaochuan Pan,et al.  Analysis of iterative region-of-interest image reconstruction for x-ray computed tomography , 2014, Journal of medical imaging.

[4]  Howard J Halpern,et al.  Comparison of parabolic filtration methods for 3D filtered back projection in pulsed EPR imaging. , 2014, Journal of magnetic resonance.

[5]  Antonin Chambolle,et al.  An introduction to continuous optimization for imaging , 2016, Acta Numerica.

[6]  E. Sidky,et al.  Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT , 2009, 0904.4495.

[7]  Yuhua Qian,et al.  3D pulse EPR imaging from sparse-view projections via constrained, total variation minimization. , 2015, Journal of magnetic resonance.

[8]  Lei Li,et al.  Optimization-based image reconstruction in computed tomography by alternating direction method with ordered subsets. , 2017, Journal of X-ray science and technology.

[9]  Boris Epel,et al.  Spectrometer manager: A versatile control software for pulse EPR spectrometers , 2005 .

[10]  Jianhua Ma,et al.  Nonlocal Prior Bayesian Tomographic Reconstruction , 2008, Journal of Mathematical Imaging and Vision.

[11]  Xiaochuan Pan,et al.  Optimization-based image reconstruction with artifact reduction in C-arm CBCT , 2016, Physics in medicine and biology.

[12]  Justin P. Haldar,et al.  Image Reconstruction From Highly Undersampled $( {\bf k}, {t})$-Space Data With Joint Partial Separability and Sparsity Constraints , 2012, IEEE Transactions on Medical Imaging.

[13]  E. Sidky,et al.  Convex optimization problem prototyping for image reconstruction in computed tomography with the Chambolle–Pock algorithm , 2011, Physics in medicine and biology.

[14]  M. Vannier,et al.  Why do commercial CT scanners still employ traditional, filtered back-projection for image reconstruction? , 2009, Inverse problems.

[15]  Qianjin Feng,et al.  Low-dose computed tomography image restoration using previous normal-dose scan. , 2011, Medical physics.

[16]  P. Joseph.  An Improved Algorithm for Reprojecting Rays through Pixel Images , 1982 .

[17]  Azeez Omotayo,et al.  Objective performance assessment of five computed tomography iterative reconstruction algorithms. , 2016, Journal of X-ray science and technology.

[18]  Xiaochuan Pan,et al.  Evaluation of sparse-view reconstruction from flat-panel-detector cone-beam CT , 2010, Physics in medicine and biology.

[19]  Xiao Han,et al.  Artifact reduction in short-scan CBCT by use of optimization-based reconstruction , 2016, Physics in medicine and biology.

[20]  Antonin Chambolle,et al.  A First-Order Primal-Dual Algorithm for Convex Problems with Applications to Imaging , 2011, Journal of Mathematical Imaging and Vision.

[21]  E. Sidky,et al.  Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization , 2008, Physics in medicine and biology.

[22]  Colin Mailer,et al.  A Versatile High Speed 250 MHz Pulse Imager for Biomedical Applications. , 2008, Concepts in magnetic resonance. Part B, Magnetic resonance engineering.

[23]  B. De Man,et al.  Distance-driven projection and backprojection in three dimensions. , 2004, Physics in medicine and biology.

[24]  Emil Y. Sidky,et al.  Investigation of optimization-based reconstruction with an image-total-variation constraint in PET , 2015, 2015 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC).