Limited-angle CT reconstruction via the L1/L2 minimization

In this paper, we consider minimizing the L1/L2 term on the gradient for a limit-angle scanning problem in computed tomography (CT) reconstruction. We design a splitting framework for both constrained and unconstrained optimization models. In addition, we can incorporate a box constraint that is reasonable for imaging applications. Numerical schemes are based on the alternating direction method of multipliers (ADMM), and we provide the convergence analysis of all the proposed algorithms (constrained/unconstrained and with/without the box constraint). Experimental results demonstrate the efficiency of our proposed approaches, showing significant improvements over the state-of-the-art methods in the limit-angle CT reconstruction. Specifically worth noticing is an exact recovery of the Shepp-Logan phantom from noiseless projection data with 30 scanning angle.

[1]  Jeffrey A. Fessler,et al.  Statistical image reconstruction for polyenergetic X-ray computed tomography , 2002, IEEE Transactions on Medical Imaging.

[2]  J. Bushberg.  The Essential Physics of Medical Imaging , 2001 .

[3]  Adrian S. Lewis,et al.  The [barred L]ojasiewicz Inequality for Nonsmooth Subanalytic Functions with Applications to Subgradient Dynamical Systems , 2006, SIAM J. Optim..

[4]  Long Wei,et al.  A Joint Reconstruction and Segmentation Method for Limited-Angle X-Ray Tomography , 2018, IEEE Access.

[5]  Steve B. Jiang,et al.  GPU-based iterative cone-beam CT reconstruction using tight frame regularization , 2010, Physics in medicine and biology.

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

[7]  Gengsheng Lawrence Zeng,et al.  Medical Image Reconstruction: A Conceptual Tutorial , 2010 .

[8]  Jinchi Lv,et al.  A unified approach to model selection and sparse recovery using regularized least squares , 2009, 0905.3573.

[9]  B. F. Logan,et al.  The Fourier reconstruction of a head section , 1974 .

[10]  Mila Nikolova,et al.  Nonconvex Optimization for 3-Dimensional Point Source Localization Using a Rotating Point Spread Function , 2019, SIAM J. Imaging Sci..

[11]  Wotao Yin,et al.  Global Convergence of ADMM in Nonconvex Nonsmooth Optimization , 2015, Journal of Scientific Computing.

[12]  Feng Liu,et al.  Graph Regularized EEG Source Imaging with In-Class Consistency and Out-Class Discrimination , 2017, IEEE Transactions on Big Data.

[13]  Rick Chartrand,et al.  Exact Reconstruction of Sparse Signals via Nonconvex Minimization , 2007, IEEE Signal Processing Letters.

[14]  Stanley Osher,et al.  Image Recovery via Nonlocal Operators , 2010, J. Sci. Comput..

[15]  J. Hiriart-Urruty,et al.  Convex analysis and minimization algorithms , 1993 .

[16]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[17]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[18]  Ming Yan,et al.  Fast L1–L2 Minimization via a Proximal Operator , 2016, Journal of Scientific Computing.

[19]  Andreas Maier,et al.  Scale-Space Anisotropic Total Variation for Limited Angle Tomography , 2017, IEEE Transactions on Radiation and Plasma Medical Sciences.

[20]  M. Nikolova.  An Algorithm for Total Variation Minimization and Applications , 2004 .

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

[22]  Chao Wang,et al.  A Scale Invariant Approach for Sparse Signal Recovery , 2018, SIAM J. Sci. Comput..

[23]  Zongben Xu,et al.  $L_{1/2}$ Regularization: A Thresholding Representation Theory and a Fast Solver , 2012, IEEE Transactions on Neural Networks and Learning Systems.

[24]  Linh V. Nguyen.  How strong are streak artifacts in limited angle computed tomography? , 2014, 1407.3037.

[25]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[26]  Xuanqin Mou,et al.  Machine Learning for Tomographic Imaging , 2020 .

[27]  Per Christian Hansen,et al.  AIR Tools - A MATLAB package of algebraic iterative reconstruction methods , 2012, J. Comput. Appl. Math..

[28]  Xin Jin,et al.  A limited-angle CT reconstruction method based on anisotropic TV minimization , 2013, Physics in medicine and biology.

[29]  Qingping Zhou,et al.  Nonlocal TV-Gaussian prior for Bayesian inverse problems with applications to limited CT reconstruction , 2020, Inverse Problems & Imaging.

[30]  Ting-Zhu Huang,et al.  Truncated l1-2 Models for Sparse Recovery and Rank Minimization , 2017, SIAM J. Imaging Sci..

[31]  Li Zhang,et al.  Low-dose multiple-information retrieval algorithm for X-ray grating-based imaging , 2011 .

[32]  Guoyin Li,et al.  Global Convergence of Splitting Methods for Nonconvex Composite Optimization , 2014, SIAM J. Optim..

[33]  Zhi-Quan Luo,et al.  Convergence analysis of alternating direction method of multipliers for a family of nonconvex problems , 2014, 2015 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[34]  Jong-Shi Pang,et al.  Decomposition Methods for Computing Directional Stationary Solutions of a Class of Nonsmooth Nonconvex Optimization Problems , 2018, SIAM J. Optim..

[35]  Jack Xin,et al.  Point Source Super-resolution Via Non-convex $$L_1$$L1 Based Methods , 2016, J. Sci. Comput..

[36]  Günter Lauritsch,et al.  Simulation tools for two-dimensional experiments in x-ray computed tomography using the FORBILD head phantom , 2012, Physics in medicine and biology.

[37]  Deren Han,et al.  Convergence of alternating direction method for minimizing sum of two nonconvex functions with linear constraints , 2017, Int. J. Comput. Math..

[38]  Lars Ruthotto,et al.  PNKH-B: A Projected Newton-Krylov Method for Large-Scale Bound-Constrained Optimization , 2021, SIAM J. Sci. Comput..

[39]  Zongben Xu,et al.  Convergence of multi-block Bregman ADMM for nonconvex composite problems , 2015, Science China Information Sciences.

[40]  Per Christian Hansen,et al.  IR Tools: a MATLAB package of iterative regularization methods and large-scale test problems , 2017, Numerical Algorithms.

[41]  Jack Xin,et al.  Minimization of transformed $$L_1$$L1 penalty: theory, difference of convex function algorithm, and robust application in compressed sensing , 2014, Math. Program..

[42]  L. Feldkamp,et al.  Practical cone-beam algorithm , 1984 .

[43]  Michael Elad,et al.  Spatially-Adaptive Reconstruction in Computed Tomography Using Neural Networks , 2013, IEEE Transactions on Medical Imaging.

[44]  Feng Liu,et al.  Estimating Latent Brain Sources with Low-Rank Representation and Graph Regularization , 2018, BI.

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

[46]  HeQi,et al.  Computing Sparse Representation in a Highly Coherent Dictionary Based on Difference of L1 and L2 , 2015 .

[47]  Xiaoqiang Luo,et al.  An Image Reconstruction Method Based on Total Variation and Wavelet Tight Frame for Limited-Angle CT , 2018, IEEE Access.

[48]  Jack Xin,et al.  Minimization of Transformed L1 Penalty: Closed Form Representation and Iterative Thresholding Algorithms , 2014, ArXiv.

[49]  M. Glas,et al.  Principles of Computerized Tomographic Imaging , 2000 .

[50]  Jeffrey A. Fessler,et al.  Fast X-Ray CT Image Reconstruction Using a Linearized Augmented Lagrangian Method With Ordered Subsets , 2014, IEEE Transactions on Medical Imaging.

[51]  Hengyong Yu,et al.  A soft-threshold filtering approach for reconstruction from a limited number of projections , 2010, Physics in medicine and biology.

[52]  Bingsheng He,et al.  The direct extension of ADMM for multi-block convex minimization problems is not necessarily convergent , 2014, Mathematical Programming.

[53]  Tom Goldstein,et al.  The Split Bregman Method for L1-Regularized Problems , 2009, SIAM J. Imaging Sci..

[54]  J. Hsieh,et al.  A three-dimensional-weighted cone beam filtered backprojection (CB-FBP) algorithm for image reconstruction in volumetric CT—helical scanning , 2006, Physics in medicine and biology.

[55]  Jean-Baptiste Thibault,et al.  A three-dimensional statistical approach to improved image quality for multislice helical CT. , 2007, Medical physics.

[56]  Martin Benning,et al.  Directional sinogram inpainting for limited angle tomography , 2018, Inverse Problems.

[57]  Eric Todd Quinto,et al.  Exterior and limited-angle tomography in non-destructive evaluation , 1998 .

[58]  Eric Todd Quinto,et al.  Characterization and reduction of artifacts in limited angle tomography , 2013 .

[59]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[60]  Ming Jiang,et al.  Convergence of the simultaneous algebraic reconstruction technique (SART) , 2003, Conference Record of Thirty-Fifth Asilomar Conference on Signals, Systems and Computers (Cat.No.01CH37256).

[61]  Lubomir M. Hadjiiski,et al.  A comparative study of limited-angle cone-beam reconstruction methods for breast tomosynthesis. , 2006, Medical physics.

[62]  D. Donoho,et al.  Sparse MRI: The application of compressed sensing for rapid MR imaging , 2007, Magnetic resonance in medicine.

[63]  Yuxiang Xing,et al.  Motion-compensated reconstruction method based on rigid motion model with multi-object , 2010 .

[64]  Steve B. Jiang,et al.  GPU-based fast cone beam CT reconstruction from undersampled and noisy projection data via total variation , 2010 .

[65]  Stephen P. Boyd,et al.  Enhancing Sparsity by Reweighted ℓ1 Minimization , 2007, 0711.1612.

[66]  Stephen P. Boyd,et al.  Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers , 2011, Found. Trends Mach. Learn..

[67]  Esa Niemi,et al.  Tomographic X-ray data of a walnut , 2015, 1502.04064.

[68]  Jack Xin,et al.  Computing Sparse Representation in a Highly Coherent Dictionary Based on Difference of $$L_1$$L1 and $$L_2$$L2 , 2015, J. Sci. Comput..

[69]  Li Zhang,et al.  Volumetric imaging from a multisegment straight-line trajectory and a practical reconstruction algorithm , 2007 .

[70]  Jack Xin,et al.  A Weighted Difference of Anisotropic and Isotropic Total Variation Model for Image Processing , 2015, SIAM J. Imaging Sci..

[71]  Eero P. Simoncelli,et al.  Image quality assessment: from error visibility to structural similarity , 2004, IEEE Transactions on Image Processing.

[72]  Jack Xin,et al.  Ratio and difference of $l_1$ and $l_2$ norms and sparse representation with coherent dictionaries , 2014, Commun. Inf. Syst..

[73]  Tony F. Chan,et al.  Image processing and analysis - variational, PDE, wavelet, and stochastic methods , 2005 .

[74]  Raymond H. Chan,et al.  A Multiplicative Iterative Algorithm for Box-Constrained Penalized Likelihood Image Restoration , 2012, IEEE Transactions on Image Processing.

[75]  Renliang Gu,et al.  Polychromatic X-ray CT Image Reconstruction and Mass-Attenuation Spectrum Estimation , 2015, 1509.02193.