Binary tomography reconstruction from few projections with Total Variation regularization for bone microstructure studies.

Discrete tomography refers to a class of reconstruction methods adapted to discrete-valued images. Many different approaches have been investigated to address the binary case, when a two-phase object is considered. This reconstruction problem is very important in medical or material applications where it is crucial to reduce the number of projections. In this paper, we address the problem of binary image reconstruction for X-ray CT imaging from a small number of projections. We propose a TV (Total Variation) regularization approach and compare the results obtained with or without an additional box convex constraint. The schemes are applied to a simple disk image and to more complex bone cross-sections for various noise levels. The minimization of the regularization functional is performed with the state-of-the-art ADMM (Alternate Direction Minimization Method) algorithm. The methods perform equally well on a simple disk image. The additional box convex constraints improves the reconstruction results for complex structures with fine details.

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

[2]  Attila Kuba,et al.  Multi-resolution method for binary tomography , 2005, Electron. Notes Discret. Math..

[3]  Thierry Blu,et al.  The SURE-LET Approach to Image Denoising , 2007, IEEE Transactions on Image Processing.

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

[5]  Hengyong Yu,et al.  A soft-threshold filtering approach for reconstruction from a limited number of projections , 2010, Physics in medicine and biology.

[6]  Burkhard Schillinger Proposed combination of CAD data and discrete tomography for the detection of coking and lubricants in turbine blades or engines , 2005, Electron. Notes Discret. Math..

[7]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[8]  Christoph Schnörr,et al.  Discrete tomography by convex-concave regularization and D.C. programming , 2005, Discret. Appl. Math..

[9]  C. Vogel Computational Methods for Inverse Problems , 1987 .

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

[11]  José M. Bioucas-Dias,et al.  An Augmented Lagrangian Approach to the Constrained Optimization Formulation of Imaging Inverse Problems , 2009, IEEE Transactions on Image Processing.

[12]  Xavier Bresson,et al.  Fast Global Minimization of the Active Contour/Snake Model , 2007, Journal of Mathematical Imaging and Vision.

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

[14]  Florent Krzakala,et al.  Belief-propagation reconstruction for discrete tomography , 2012, ArXiv.

[15]  Emmanuel J. Candès,et al.  NESTA: A Fast and Accurate First-Order Method for Sparse Recovery , 2009, SIAM J. Imaging Sci..

[16]  G. Herman,et al.  Advances in discrete tomography and its applications , 2007 .

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

[18]  Kees Joost Batenburg,et al.  Generic iterative subset algorithms for discrete tomography , 2009, Discret. Appl. Math..

[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]  Xavier Bresson,et al.  Completely Convex Formulation of the Chan-Vese Image Segmentation Model , 2012, International Journal of Computer Vision.

[21]  Ralph Müller,et al.  Guidelines for assessment of bone microstructure in rodents using micro–computed tomography , 2010, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[22]  M. Defrise,et al.  An algorithm for total variation regularization in high-dimensional linear problems , 2011 .

[23]  Gabor T. Herman,et al.  Automated estimation of the parameters of Gibbs priors to be used in binary tomography , 2004, Discret. Appl. Math..

[24]  José M. Bioucas-Dias,et al.  Fast Image Recovery Using Variable Splitting and Constrained Optimization , 2009, IEEE Transactions on Image Processing.

[25]  E. Boller,et al.  Relevance of 2D radiographic texture analysis for the assessment of 3D bone micro-architecture. , 2006, Medical physics.

[26]  Jan Sijbers,et al.  The ASTRA Tomography Toolbox , 2013 .

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

[28]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[29]  Lin Ma,et al.  Comparison of approaches based on optimization and algebraic iteration for binary tomography , 2010, Comput. Phys. Commun..

[30]  Kees Joost Batenburg,et al.  A discrete tomography approach for superresolution micro-CT images: application to bone , 2010, 2010 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[31]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[32]  Kees Joost Batenburg,et al.  DART: A Practical Reconstruction Algorithm for Discrete Tomography , 2011, IEEE Transactions on Image Processing.

[33]  Yonina C. Eldar Generalized SURE for Exponential Families: Applications to Regularization , 2008, IEEE Transactions on Signal Processing.