Multi-Channel Potts-Based Reconstruction for Multi-Spectral Computed Tomography

We consider reconstructing multi-channel images from measurements performed by photon-counting and energy-discriminating detectors in the setting of multi-spectral X-ray computed tomography (CT). Our aim is to exploit the strong structural correlation that is known to exist between the channels of multi-spectral CT images. To that end, we adopt the multi-channel Potts prior to jointly reconstruct all channels. This prior produces piecewise constant solutions with strongly correlated channels. In particular, edges are enforced to have the same spatial position across channels which is a benefit over TV-based methods. We consider the Potts prior in two frameworks: (a) in the context of a variational Potts model, and (b) in a Potts-superiorization approach that perturbs the iterates of a basic iterative least squares solver. We identify an alternating direction method of multipliers (ADMM) approach as well as a Potts-superiorized conjugate gradient method as particularly suitable. In numerical experiments, we compare the Potts prior based approaches to existing TV-type approaches on realistically simulated multi-spectral CT data and obtain improved reconstruction for compound solid bodies.

[1]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[2]  Junfeng Yang,et al.  Alternating Direction Algorithms for 1-Problems in Compressive Sensing , 2009, SIAM J. Sci. Comput..

[3]  Olga Veksler,et al.  Fast approximate energy minimization via graph cuts , 2001, Proceedings of the Seventh IEEE International Conference on Computer Vision.

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

[5]  Ronny Ramlau,et al.  A Mumford-Shah level-set approach for the inversion and segmentation of X-ray tomography data , 2007, J. Comput. Phys..

[6]  R. Ramlau,et al.  A Mumford-Shah level-set approach for the inversion andsegmentation of SPECT/CT data , 2011 .

[7]  Philip J. Withers,et al.  Joint image reconstruction method with correlative multi-channel prior for x-ray spectral computed tomography , 2018 .

[8]  Vladimir Kolmogorov,et al.  Computing geodesics and minimal surfaces via graph cuts , 2003, Proceedings Ninth IEEE International Conference on Computer Vision.

[9]  Chao Yang,et al.  Alternating direction methods for classical and ptychographic phase retrieval , 2012 .

[10]  Tianye Niu,et al.  Image‐domain multimaterial decomposition for dual‐energy CT based on prior information of material images , 2017, Medical physics.

[11]  R. Bellman,et al.  Curve Fitting by Segmented Straight Lines , 1969 .

[12]  Olaf Wittich,et al.  Complexity penalized least squares estimators: Analytical results , 2008 .

[13]  Ken D. Sauer,et al.  A local update strategy for iterative reconstruction from projections , 1993, IEEE Trans. Signal Process..

[14]  G. Herman,et al.  Superiorization of Preconditioned Conjugate Gradient Algorithms for Tomographic Image Reconstruction , 2018, 1807.10151.

[15]  Michael Unser,et al.  Joint image reconstruction and segmentation using the Potts model , 2014, 1405.5850.

[16]  P. Shikhaliev.  Energy-resolved computed tomography: first experimental results , 2008, Physics in medicine and biology.

[17]  Gabriele Steidl,et al.  Removing Multiplicative Noise by Douglas-Rachford Splitting Methods , 2010, Journal of Mathematical Imaging and Vision.

[18]  R. Ramlau,et al.  Regularization of ill-posed Mumford–Shah models with perimeter penalization , 2010 .

[19]  Ran Davidi,et al.  Superiorization: An optimization heuristic for medical physics , 2012, Medical physics.

[20]  Massimo Fornasier,et al.  Iterative Thresholding Meets Free-Discontinuity Problems , 2009, Found. Comput. Math..

[21]  Christoph Schnörr,et al.  Superiorization vs. Accelerated Convex Optimization: The Superiorized/Regularized Least-Squares Case , 2019, ArXiv.

[22]  Jeffrey D. Scargle,et al.  An algorithm for optimal partitioning of data on an interval , 2003, IEEE Signal Processing Letters.

[23]  R. Zabih,et al.  Efficient Graph-Based Energy Minimization Methods in Computer Vision , 1999 .

[24]  Ville Kolehmainen,et al.  Joint reconstruction in low dose multi-energy CT , 2019, Inverse Problems & Imaging.

[25]  Martin Storath,et al.  Smoothing for signals with discontinuities using higher order Mumford–Shah models , 2018, Numerische Mathematik.

[26]  Alan S. Willsky,et al.  A curve evolution-based variational approach to simultaneous image restoration and segmentation , 2002, Proceedings. International Conference on Image Processing.

[27]  S. Osher,et al.  Multi-energy CT based on a prior rank, intensity and sparsity model (PRISM) , 2011, Inverse problems.

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

[29]  Andrew Blake,et al.  Comparison of the Efficiency of Deterministic and Stochastic Algorithms for Visual Reconstruction , 1989, IEEE Trans. Pattern Anal. Mach. Intell..

[30]  Andreas Weinmann,et al.  Fast Partitioning of Vector-Valued Images , 2014, SIAM J. Imaging Sci..

[31]  Simon Rit,et al.  Regularization of nonlinear decomposition of spectral x‐ray projection images , 2017, Medical physics.

[32]  Nahum Kiryati,et al.  Variational Pairing of Image Segmentation and Blind Restoration , 2004, ECCV.

[33]  A. Chambolle.  FINITE-DIFFERENCES DISCRETIZATIONS OF THE MUMFORD-SHAH FUNCTIONAL , 1999 .

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

[35]  G. Winkler,et al.  Complexity Penalized M-Estimation , 2008 .

[36]  Gabor T. Herman,et al.  Total variation superiorized conjugate gradient method for image reconstruction , 2017, 1709.04912.

[37]  Andreas Weinmann,et al.  Jump-Sparse and Sparse Recovery Using Potts Functionals , 2013, IEEE Transactions on Signal Processing.

[38]  Eric L. Miller,et al.  Tensor-Based Formulation and Nuclear Norm Regularization for Multienergy Computed Tomography , 2013, IEEE Transactions on Image Processing.

[39]  Zheng Xu,et al.  An Empirical Study of ADMM for Nonconvex Problems , 2016, ArXiv.

[40]  Martin Storath,et al.  An algorithmic framework for Mumford–Shah regularization of inverse problems in imaging , 2015 .

[41]  David S. Lalush,et al.  Full-Spectrum CT Reconstruction Using a Weighted Least Squares Algorithm With an Energy-Axis Penalty , 2011, IEEE Transactions on Medical Imaging.

[42]  Antonin Chambolle,et al.  Image Segmentation by Variational Methods: Mumford and Shah Functional and the Discrete Approximations , 1995, SIAM J. Appl. Math..

[43]  Michael K. Ng,et al.  Solving Constrained Total-variation Image Restoration and Reconstruction Problems via Alternating Direction Methods , 2010, SIAM J. Sci. Comput..

[44]  Jeffrey A. Fessler,et al.  Multi-Material Decomposition Using Statistical Image Reconstruction for Spectral CT , 2014, IEEE Transactions on Medical Imaging.

[45]  J H Siewerdsen,et al.  Spektr: a computational tool for x-ray spectral analysis and imaging system optimization. , 2004, Medical physics.

[46]  Gabor T. Herman,et al.  Derivative-free superiorization: principle and algorithm , 2019, Numerical Algorithms.

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

[48]  Yair Censor,et al.  Superiorization and Perturbation Resilience of Algorithms: A Continuously Updated Bibliography , 2015, 1506.04219.

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

[50]  Martin Storath,et al.  Iterative Potts and Blake–Zisserman minimization for the recovery of functions with discontinuities from indirect measurements , 2015, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[51]  G. Winkler,et al.  Smoothers for Discontinuous Signals , 2002 .

[52]  Brendt Wohlberg,et al.  A nonconvex ADMM algorithm for group sparsity with sparse groups , 2013, 2013 IEEE International Conference on Acoustics, Speech and Signal Processing.

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

[54]  Xiaojun Chen,et al.  Alternating Direction Method of Multipliers for a Class of Nonconvex and Nonsmooth Problems with Applications to Background/Foreground Extraction , 2015, SIAM J. Imaging Sci..

[55]  Jan Sijbers,et al.  The ASTRA Toolbox: A platform for advanced algorithm development in electron tomography. , 2015, Ultramicroscopy.

[56]  Esther Klann,et al.  A Mumford-Shah-Like Method for Limited Data Tomography with an Application to Electron Tomography , 2011, SIAM J. Imaging Sci..

[57]  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).

[58]  I E Auger,et al.  Algorithms for the optimal identification of segment neighborhoods. , 1989, Bulletin of mathematical biology.

[59]  Patrick J La Rivière,et al.  Joint reconstruction of multi-channel, spectral CT data via constrained total nuclear variation minimization , 2014, Physics in medicine and biology.

[60]  Ran Davidi,et al.  Perturbation-resilient block-iterative projection methods with application to image reconstruction from projections , 2009, Int. Trans. Oper. Res..

[61]  Andreas Weinmann,et al.  Iterative Potts Minimization for the Recovery of Signals with Discontinuities from Indirect Measurements: The Multivariate Case , 2018, Found. Comput. Math..

[62]  Hong-Kun Xu,et al.  Convergence of Bregman alternating direction method with multipliers for nonconvex composite problems , 2014, 1410.8625.