TICMR: Total Image Constrained Material Reconstruction via Nonlocal Total Variation Regularization for Spectral CT

This work develops a material reconstruction method for spectral CT, namely Total Image Constrained Material Reconstruction (TICMR), to maximize the utility of projection data in terms of both spectral information and high signal-to-noise ratio (SNR). This is motivated by the following fact: when viewed as a spectrally-integrated measurement, the projection data can be used to reconstruct a total image without spectral information, which however has a relatively high SNR; when viewed as a spectrally-resolved measurement, the projection data can be utilized to reconstruct the material composition, which however has a relatively low SNR. The material reconstruction synergizes material decomposition and image reconstruction, i.e., the direct reconstruction of material compositions instead of a two-step procedure that first reconstructs images and then decomposes images. For material reconstruction with high SNR, we propose TICMR with nonlocal total variation (NLTV) regularization. That is, first we reconstruct a total image using spectrally-integrated measurement without spectral binning, and build the NLTV weights from this image that characterize nonlocal image features; then the NLTV weights are incorporated into a NLTV-based iterative material reconstruction scheme using spectrally-binned projection data, so that these weights serve as a high-SNR reference to regularize material reconstruction. Note that the nonlocal property of NLTV is essential for material reconstruction, since material compositions may have significant local intensity variations although their structural information is often similar. In terms of solution algorithm, TICMR is formulated as an iterative reconstruction method with the NLTV regularization, in which the nonlocal divergence is utilized based on the adjoint relationship. The alternating direction method of multipliers is developed to solve this sparsity optimization problem. The proposed TICMR method was validated using both simulated and experimental data. In comparison with FBP and total-variation-based iterative method, TICMR had improved image quality, e.g., contrast-to-noise ratio and spatial resolution.

[1]  Jun Zhao,et al.  United Iterative Reconstruction for Spectral Computed Tomography , 2015, IEEE Transactions on Medical Imaging.

[2]  Jong Chul Ye,et al.  Sparse-View Spectral CT Reconstruction Using Spectral Patch-Based Low-Rank Penalty , 2015, IEEE Transactions on Medical Imaging.

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

[4]  Xavier Bresson,et al.  Bregmanized Nonlocal Regularization for Deconvolution and Sparse Reconstruction , 2010, SIAM J. Imaging Sci..

[5]  Jian-Feng Cai,et al.  Adaptive tight frame based medical image reconstruction: a proof-of-concept study for computed tomography , 2013 .

[6]  Axel Thran,et al.  Statistical Reconstruction of Material Decomposed Data in Spectral CT , 2013, IEEE Transactions on Medical Imaging.

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

[8]  Steve B. Jiang,et al.  Cine Cone Beam CT Reconstruction Using Low-Rank Matrix Factorization: Algorithm and a Proof-of-Principle Study , 2012, IEEE Transactions on Medical Imaging.

[9]  Ge Wang,et al.  Dual-dictionary learning-based iterative image reconstruction for spectral computed tomography application , 2012, Physics in medicine and biology.

[10]  Joachim Hornegger,et al.  Guided Noise Reduction for Spectral CT with Energy-Selective Photon Counting Detectors , 2013 .

[11]  Darin P Clark,et al.  Spectral diffusion: an algorithm for robust material decomposition of spectral CT data , 2014, Physics in medicine and biology.

[12]  Predrag Sukovic,et al.  Penalized weighted least-squares image reconstruction for dual energy X-ray transmission tomography , 2000, IEEE Transactions on Medical Imaging.

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

[14]  L. Xing,et al.  4D cone beam CT via spatiotemporal tensor framelet. , 2012, Medical physics.

[15]  Mark A Anastasio,et al.  Sparsity-regularized image reconstruction of decomposed K-edge data in spectral CT , 2014, Physics in Medicine and Biology.

[16]  Daniel A Low,et al.  Megavoltage CT imaging quality improvement on TomoTherapy via tensor framelet. , 2013, Medical physics.

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

[18]  Jean-Michel Morel,et al.  A Review of Image Denoising Algorithms, with a New One , 2005, Multiscale Model. Simul..

[19]  K. Taguchi,et al.  Vision 20/20: Single photon counting x-ray detectors in medical imaging. , 2013, Medical physics.

[20]  Sabee Molloi,et al.  Breast tissue characterization with photon-counting spectral CT imaging: a postmortem breast study. , 2014, Radiology.

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

[22]  C. Mistretta,et al.  Noise reduction in spectral CT: reducing dose and breaking the trade-off between image noise and energy bin selection. , 2011, Medical physics.

[23]  Taly Gilat Schmidt,et al.  Optimal "image-based" weighting for energy-resolved CT. , 2009, Medical physics.

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

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

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

[27]  Hyo-Min Cho,et al.  Characterization of energy response for photon-counting detectors using x-ray fluorescence. , 2014, Medical physics.

[28]  Bo Zhao,et al.  Tight-frame based iterative image reconstruction for spectral breast CT. , 2013, Medical physics.

[29]  Michaela C. C. Weigel,et al.  High-resolution spiral CT of the breast at very low dose: concept and feasibility considerations , 2011, European Radiology.

[30]  Yuxiang Xing,et al.  Multienergy CT acquisition and reconstruction with a stepped tube potential scan. , 2014, Medical physics.

[31]  Hongkai Zhao,et al.  Robust principal component analysis-based four-dimensional computed tomography , 2011, Physics in medicine and biology.

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

[33]  Hyo-Min Cho,et al.  A high-resolution photon-counting breast CT system with tensor-framelet based iterative image reconstruction for radiation dose reduction. , 2014, Physics in medicine and biology.

[34]  Hao Gao.  Fast parallel algorithms for the x-ray transform and its adjoint. , 2012, Medical physics.

[35]  M. Danielsson,et al.  Photon-counting spectral computed tomography using silicon strip detectors: a feasibility study , 2010, Physics in medicine and biology.

[36]  Guy Gilboa,et al.  Nonlocal Operators with Applications to Image Processing , 2008, Multiscale Model. Simul..

[37]  Zhiqiang Chen,et al.  A tensor PRISM algorithm for multi-energy CT reconstruction and comparative studies. , 2014, Journal of X-ray science and technology.