Iterative CT reconstruction via minimizing adaptively reweighted total variation.

BACKGROUND Iterative reconstruction via total variation (TV) minimization has demonstrated great successes in accurate CT imaging from under-sampled projections. When projections are further reduced, over-smoothing artifacts appear in the current reconstruction especially around the structure boundaries. OBJECTIVE We propose a practical algorithm to improve TV-minimization based CT reconstruction on very few projection data. METHOD Based on the theory of compressed sensing, the L-0 norm approach is more desirable to further reduce the projection views. To overcome the computational difficulty of the non-convex optimization of the L-0 norm, we implement an adaptive weighting scheme to approximate the solution via a series of TV minimizations for practical use in CT reconstruction. The weight on TV is initialized as uniform ones, and is automatically changed based on the gradient of the reconstructed image from the previous iteration. The iteration stops when a small difference between the weighted TV values is observed on two consecutive reconstructed images. RESULTS We evaluate the proposed algorithm on both a digital phantom and a physical phantom. Using 20 equiangular projections, our method reduces reconstruction errors in the conventional TV minimization by a factor of more than 5, with improved spatial resolution. CONCLUSIONS By adaptively reweighting TV in iterative CT reconstruction, we successfully further reduce the projection number for the same or better image quality.

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

[2]  Jie Tang,et al.  Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets. , 2008, Medical physics.

[3]  Guang-Hong Chen,et al.  Reduced image noise at low-dose multidetector CT of the abdomen with prior image constrained compressed sensing algorithm. , 2011, Radiology.

[4]  Mário A. T. Figueiredo,et al.  Gradient Projection for Sparse Reconstruction: Application to Compressed Sensing and Other Inverse Problems , 2007, IEEE Journal of Selected Topics in Signal Processing.

[5]  Armando Manduca,et al.  Practical nonconvex Compressive Sensing reconstruction of highly-accelerated 3D parallel MR angiograms , 2009, 2009 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[6]  Michael Elad,et al.  Image Denoising Via Sparse and Redundant Representations Over Learned Dictionaries , 2006, IEEE Transactions on Image Processing.

[7]  Lei Zhu,et al.  Shading correction for on-board cone-beam CT in radiation therapy using planning MDCT images. , 2010, Medical physics.

[8]  B. Fahimian,et al.  Enhancement of four-dimensional cone-beam computed tomography by compressed sensing with Bregman iteration. , 2013, Journal of X-ray science and technology.

[9]  Sajid Abbas,et al.  Super-sparsely view-sampled cone-beam CT by incorporating prior data. , 2013, Journal of X-ray science and technology.

[10]  Klaus Mueller,et al.  IOP PUBLISHING PHYSICS IN MEDICINE AND BIOLOGY , 2007 .

[11]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

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

[13]  Steve B. Jiang,et al.  A comprehensive study on the relationship between the image quality and imaging dose in low-dose cone beam CT , 2011, Physics in medicine and biology.

[14]  Lei Xing,et al.  Using total-variation regularization for intensity modulated radiation therapy inverse planning with field-specific numbers of segments , 2008, Physics in medicine and biology.

[15]  Lei Zhang,et al.  Low-Dose X-ray CT Reconstruction via Dictionary Learning , 2012, IEEE Transactions on Medical Imaging.

[16]  Xiaochuan Pan,et al.  Image reconstruction from few views by non-convex optimization , 2007, 2007 IEEE Nuclear Science Symposium Conference Record.

[17]  Lei Zhu,et al.  Faster STORM using compressed sensing , 2012, Nature Methods.

[18]  M. Kachelriess,et al.  Improved total variation-based CT image reconstruction applied to clinical data , 2011, Physics in medicine and biology.

[19]  Steve B. Jiang,et al.  Low-dose CT reconstruction via edge-preserving total variation regularization. , 2010, Physics in medicine and biology.

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

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

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

[23]  Lei Xing,et al.  Total-Variation Regularization Based Inverse Planning for Intensity Modulated Arc Therapy , 2012, Technology in cancer research & treatment.

[24]  Liang Li,et al.  A few-view reweighted sparsity hunting (FRESH) method for CT image reconstruction. , 2013, Journal of X-ray science and technology.

[25]  Lei Zhu,et al.  Accelerated barrier optimization compressed sensing (ABOCS) reconstruction for cone-beam CT: Phantom studies. , 2012, Medical physics.

[26]  Xiao Han,et al.  Optimization-based reconstruction of sparse images from few-view projections , 2012, Physics in medicine and biology.

[27]  Jin Sung Kim,et al.  Fast compressed sensing-based CBCT reconstruction using Barzilai-Borwein formulation for application to on-line IGRT. , 2012, Medical physics.

[28]  G C Sharp,et al.  GPU-based streaming architectures for fast cone-beam CT image reconstruction and demons deformable registration , 2007, Physics in medicine and biology.

[29]  Tao Yang,et al.  GPU based iterative cone-beam CT reconstruction using empty space skipping technique. , 2013, Journal of X-ray science and technology.

[30]  Lei Zhu,et al.  Scatter Correction Method for X-Ray CT Using Primary Modulation: Theory and Preliminary Results , 2006, IEEE Transactions on Medical Imaging.

[31]  Li Zeng,et al.  A new iterative reconstruction algorithm for 2D exterior fan-beam CT. , 2010, Journal of X-ray science and technology.

[32]  Emmanuel J. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[33]  Jun Zhao,et al.  TV-based image reconstruction of multiple objects in a fixed source-detector geometry. , 2012, Journal of X-ray science and technology.

[34]  Søren Holdt Jensen,et al.  Implementation of an optimal first-order method for strongly convex total variation regularization , 2011, ArXiv.

[35]  Jiehua Zhu,et al.  Convergence of block cyclic projection and Cimmino algorithms for compressed sensing based tomography. , 2010, Journal of X-ray science and technology.

[36]  Radhe Mohan,et al.  Scatter kernel estimation with an edge-spread function method for cone-beam computed tomography imaging , 2008, Physics in medicine and biology.

[37]  Lei Zhu,et al.  Scatter correction for full-fan volumetric CT using a stationary beam blocker in a single full scan. , 2011, Medical physics.

[38]  J Trzasko,et al.  Nonconvex prior image constrained compressed sensing (NCPICCS): theory and simulations on perfusion CT. , 2011, Medical physics.

[39]  Lei Zhu,et al.  Quantitative cone-beam CT imaging in radiation therapy using planning CT as a prior: first patient studies. , 2012, Medical physics.

[40]  Michael Elad,et al.  Double Sparsity: Learning Sparse Dictionaries for Sparse Signal Approximation , 2010, IEEE Transactions on Signal Processing.

[41]  Xin Jin,et al.  Experimental studies on few-view reconstruction for high-resolution micro-CT. , 2013, Journal of X-ray science and technology.