Non-convex primal-dual algorithm for image reconstruction in spectral CT

The work seeks to develop an algorithm for image reconstruction by directly inverting the non-linear data model in spectral CT. Using the non-linear data model, we formulate the image-reconstruction problem as a non-convex optimization program, and develop a non-convex primal-dual (NCPD) algorithm to solve the program. We devise multiple convergence conditions and perform verification studies numerically to demonstrate that the NCPD algorithm can solve the non-convex optimization program and under appropriate data condition, can invert the non-linear data model. Using the NCPD algorithm, we then reconstruct monochromatic images from simulated and real data of numerical and physical phantoms acquired with a standard, full-scan dual-energy configuration. The result of the reconstruction studies shows that the NCPD algorithm can correct accurately for the non-linear beam-hardening effect. Furthermore, we apply the NCPD algorithm to simulated and real data of the numerical and physical phantoms collected with non-standard, short-scan dual-energy configurations, and obtain monochromatic images comparable to those of the standard, full-scan study, thus revealing the potential of the NCPD algorithm for enabling non-standard scanning configurations in spectral CT, where the existing indirect methods are limited.

[1]  Xiaoming Yuan,et al.  Adaptive Primal-Dual Hybrid Gradient Methods for Saddle-Point Problems , 2013, 1305.0546.

[2]  Emil Y. Sidky,et al.  MOCCA: Mirrored Convex/Concave Optimization for Nonconvex Composite Functions , 2015, J. Mach. Learn. Res..

[3]  Yu Zou,et al.  Analysis of fast kV-switching in dual energy CT using a pre-reconstruction decomposition technique , 2008, SPIE Medical Imaging.

[4]  E. Roessl,et al.  K-edge imaging in x-ray computed tomography using multi-bin photon counting detectors , 2007, Physics in medicine and biology.

[5]  Katsuyuki Taguchi,et al.  Joint estimation of tissue types and linear attenuation coefficients for photon counting CT. , 2015, Medical physics.

[6]  Xiaochuan Pan,et al.  Image reconstruction and scan configurations enabled by optimization-based algorithms in multispectral CT , 2017, Physics in medicine and biology.

[7]  A. Macovski,et al.  Energy-selective reconstructions in X-ray computerised tomography , 1976, Physics in medicine and biology.

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

[9]  Peng Zhang,et al.  An Extended Algebraic Reconstruction Technique (E-ART) for Dual Spectral CT , 2015, IEEE Transactions on Medical Imaging.

[10]  Emil Y. Sidky,et al.  Optimization-based reconstruction for correcting non-linear partial volume artifacts in CT , 2019, Medical Imaging.

[11]  Mark A. Anastasio,et al.  Proximal ADMM for Multi-Channel Image Reconstruction in Spectral X-ray CT , 2014, IEEE Transactions on Medical Imaging.

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

[13]  Thomas Brox,et al.  iPiano: Inertial Proximal Algorithm for Nonconvex Optimization , 2014, SIAM J. Imaging Sci..

[14]  Xiaochuan Pan,et al.  Algorithm-enabled partial-angular-scan configurations for dual-energy CT. , 2018, Medical physics.

[15]  Andrew M Hernandez,et al.  Tungsten anode spectral model using interpolating cubic splines: unfiltered x-ray spectra from 20 kV to 640 kV. , 2014, Medical physics.

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

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

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

[19]  Simon Rit,et al.  Comparison of five one-step reconstruction algorithms for spectral CT , 2018, Physics in medicine and biology.

[20]  Ken D. Sauer,et al.  Model-Based Iterative Reconstruction for Dual-Energy X-Ray CT Using a Joint Quadratic Likelihood Model , 2014, IEEE Transactions on Medical Imaging.

[21]  Tuomo Valkonen,et al.  A primal–dual hybrid gradient method for nonlinear operators with applications to MRI , 2013, 1309.5032.

[22]  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.

[23]  Yoram Singer,et al.  Efficient projections onto the l1-ball for learning in high dimensions , 2008, ICML '08.

[24]  A. Mohammad-Djafari,et al.  A full-spectral Bayesian reconstruction approach based on the material decomposition model applied in dual-energy computed tomography. , 2013, Medical physics.

[25]  Xiaochuan Pan,et al.  An algorithm for constrained one-step inversion of spectral CT data , 2015, Physics in medicine and biology.

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

[27]  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.