Optimization-based reconstruction of sparse images from few-view projections

In this work, we investigate optimization-based image reconstruction from few-view (i.e. less than ten views) projections of sparse objects such as coronary-artery specimens. Using optimization programs as a guide, we formulate constraint programs as reconstruction programs and develop algorithms to reconstruct images through solving the reconstruction programs. Characterization studies are carried out for elucidating the algorithm properties of 'convergence' (relative to designed solutions) and 'utility' (relative to desired solutions) by using simulated few-view data calculated from a discrete FORBILD coronary-artery phantom, and real few-view data acquired from a human coronary-artery specimen. Study results suggest that carefully designed reconstruction programs and algorithms can yield accurate reconstructions of sparse images from few-view projections.

[1]  Timothy L. Kline,et al.  Accuracy of Microvascular Measurements Obtained From Micro-CT Images , 2010, Annals of Biomedical Engineering.

[2]  B. De Man,et al.  Distance-driven projection and backprojection in three dimensions. , 2004, Physics in medicine and biology.

[3]  R. Kress,et al.  Inverse Acoustic and Electromagnetic Scattering Theory , 1992 .

[4]  M. Sezan,et al.  Image Restoration by the Method of Convex Projections: Part 2-Applications and Numerical Results , 1982, IEEE Transactions on Medical Imaging.

[5]  K W Leszczynski,et al.  On few-view tomographic reconstruction with megavoltage photon beams. , 2001, Medical physics.

[6]  Meihua Li,et al.  Improved iterative algorithm for sparse object reconstruction and its performance evaluation with micro-CT data , 2004, IEEE Transactions on Nuclear Science.

[7]  Xiaochuan Pan,et al.  A constrained, total-variation minimization algorithm for low-intensity x-ray CT. , 2011, Medical physics.

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

[9]  Xiaochuan Pan,et al.  Image reconstruction exploiting object sparsity in boundary-enhanced X-ray phase-contrast tomography. , 2010, Optics express.

[10]  R. Lewitt Alternatives to voxels for image representation in iterative reconstruction algorithms , 1992, Physics in medicine and biology.

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

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

[13]  Xiaochuan Pan,et al.  Evaluation of sparse-view reconstruction from flat-panel-detector cone-beam CT , 2010, Physics in medicine and biology.

[14]  Hiroyuki Kudo,et al.  An accurate iterative reconstruction algorithm for sparse objects: application to 3D blood vessel reconstruction from a limited number of projections. , 2002, Physics in medicine and biology.

[15]  Armand Wirgin,et al.  The inverse crime , 2004, math-ph/0401050.

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

[17]  Marc Teboulle,et al.  Fast Gradient-Based Algorithms for Constrained Total Variation Image Denoising and Deblurring Problems , 2009, IEEE Transactions on Image Processing.

[18]  T. M. Peters Algorithms for Fast Back- and Re-Projection in Computed Tomography , 1981 .

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

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

[21]  Han Wang,et al.  Image representation by blob and its application in CT reconstruction from few projections , 2011, ArXiv.

[22]  R. Chartrand,et al.  Restricted isometry properties and nonconvex compressive sensing , 2007 .

[23]  T. K. Narayan,et al.  Evaluation of task-oriented performance of several fully 3D PET reconstruction algorithms. , 1994, Physics in medicine and biology.

[24]  Robert M. Lewitt,et al.  Practical considerations for 3-D image reconstruction using spherically symmetric volume elements , 1996, IEEE Trans. Medical Imaging.

[25]  Michael Grass,et al.  High-quality 3-D coronary artery imaging on an interventional C-arm x-ray system. , 2010, Medical physics.

[26]  Emil Y. Sidky,et al.  Quantification of admissible undersampling for sparsity-exploiting iterative image reconstruction in X-ray CT , 2011 .

[27]  P. Joseph An Improved Algorithm for Reprojecting Rays through Pixel Images , 1982 .

[28]  Xiaochuan Pan,et al.  Enhanced imaging of microcalcifications in digital breast tomosynthesis through improved image-reconstruction algorithms. , 2009, Medical physics.

[29]  Emmanuel J. Candès,et al.  Templates for convex cone problems with applications to sparse signal recovery , 2010, Math. Program. Comput..

[30]  P. L. Combettes The foundations of set theoretic estimation , 1993 .

[31]  Michael Grass,et al.  Evaluation of Iterative Sparse Object Reconstruction From Few Projections for 3-D Rotational Coronary Angiography , 2008, IEEE Transactions on Medical Imaging.

[32]  Emil Y. Sidky,et al.  Algorithm-Enabled Low-Dose Micro-CT Imaging , 2011, IEEE Transactions on Medical Imaging.

[33]  Mike E. Davies,et al.  Iterative Hard Thresholding for Compressed Sensing , 2008, ArXiv.

[34]  S M Jorgensen,et al.  Three-dimensional imaging of vasculature and parenchyma in intact rodent organs with X-ray micro-CT. , 1998, The American journal of physiology.

[35]  E. Somersalo,et al.  Statistical inverse problems: discretization, model reduction and inverse crimes , 2007 .

[36]  P. Joseph An Improved Algorithm for Reprojecting Rays through Pixel Images , 1983, IEEE Transactions on Medical Imaging.

[37]  C. M. Chen,et al.  Incremental algorithm-a new fast backprojection scheme for parallel beam geometries. , 1990, IEEE transactions on medical imaging.

[38]  D. Youla,et al.  Image Restoration by the Method of Convex Projections: Part 1ߞTheory , 1982, IEEE Transactions on Medical Imaging.

[39]  Jie Tang,et al.  Performance comparison between total variation (TV)-based compressed sensing and statistical iterative reconstruction algorithms , 2009, Physics in medicine and biology.

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

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

[42]  Emil Y. Sidky,et al.  Quantifying Admissible Undersampling for Sparsity-Exploiting Iterative Image Reconstruction in X-Ray CT , 2011, IEEE Transactions on Medical Imaging.