Belief-propagation reconstruction for discrete tomography

We consider the reconstruction of a two-dimensional discrete image from a set of tomographic measurements corresponding to the Radon projection. Assuming that the image has a structure where neighbouring pixels have a larger probability to take the same value, we follow a Bayesian approach and introduce a fast message-passing reconstruction algorithm based on belief propagation. For numerical results, we specialize to the case of binary tomography. We test the algorithm on binary synthetic images with different length scales and compare our results against a more usual convex optimization approach. We investigate the reconstruction error as a function of the number of tomographic measurements, corresponding to the number of projection angles. The belief propagation algorithm turns out to be more efficient than the convex-optimization algorithm, both in terms of recovery bounds for noise-free projections, and in terms of reconstruction quality when moderate Gaussian noise is added to the projections.

[1]  E. Lieb,et al.  Mathematical Physics in One Dimension: Exactly Soluble Models of Interacting Particles , 1967 .

[2]  G. Herman,et al.  Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. , 1970, Journal of theoretical biology.

[3]  D. Haar,et al.  Statistical Physics , 1971, Nature.

[4]  G. N. Ramachandran,et al.  Three-dimensional reconstruction from radiographs and electron micrographs: application of convolutions instead of Fourier transforms. , 1971, Proceedings of the National Academy of Sciences of the United States of America.

[5]  R. Feynman.  Statistical Mechanics, A Set of Lectures , 1972 .

[6]  B. F. Logan,et al.  The Fourier reconstruction of a head section , 1974 .

[7]  B. Logan,et al.  Reconstructing Interior Head Tissue from X-Ray Transmissions , 1974 .

[8]  O. Seeberg.  Statistical Mechanics. — A Set of Lectures , 1975 .

[9]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[10]  P. McMullen,et al.  On Hammer's X-Ray Problem , 1980 .

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

[12]  Judea Pearl,et al.  Reverend Bayes on Inference Engines: A Distributed Hierarchical Approach , 1982, AAAI.

[13]  Jean Serra,et al.  Image Analysis and Mathematical Morphology , 1983 .

[14]  H. Kudo,et al.  FROM A SMALL NUMBER OF PROJECTIONS , 1988 .

[15]  R. Gardner.  Geometric Tomography: Parallel X-rays of planar convex bodies , 2006 .

[16]  S. Kak.  Information, physics, and computation , 1996 .

[17]  Attila Kuba,et al.  Binary vectors partially determined by linear equation systems , 1997, Discret. Math..

[18]  Yoram Bresler,et al.  Globally convergent edge-preserving regularized reconstruction: an application to limited-angle tomography , 1998, IEEE Trans. Image Process..

[19]  G. Herman,et al.  Discrete tomography : foundations, algorithms, and applications , 1999 .

[20]  Olga Veksler,et al.  Fast approximate energy minimization via graph cuts , 2001, Proceedings of the Seventh IEEE International Conference on Computer Vision.

[21]  Bruno M. Carvalho,et al.  Binary Tomography for Triplane Cardiography , 1999, IPMI.

[22]  M. Glas,et al.  Principles of Computerized Tomographic Imaging , 2000 .

[23]  Brendan J. Frey,et al.  Factor graphs and the sum-product algorithm , 2001, IEEE Trans. Inf. Theory.

[24]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[25]  Gabor T. Herman,et al.  Automated estimation of the parameters of Gibbs priors to be used in binary tomography , 2001, Electron. Notes Theor. Comput. Sci..

[26]  Olga Veksler,et al.  Fast Approximate Energy Minimization via Graph Cuts , 2001, IEEE Trans. Pattern Anal. Mach. Intell..

[27]  William T. Freeman,et al.  Understanding belief propagation and its generalizations , 2003 .

[28]  Christoph Schnörr,et al.  A Linear Programming Relaxation for Binary Tomography with Smoothness Priors , 2003, Electron. Notes Discret. Math..

[29]  Christoph Schnörr,et al.  Binary Tomography by Iterating Linear Programs from Noisy Projections , 2004, IWCIA.

[30]  Gabor T. Herman,et al.  A method for reconstructing label images from a few projections, as motivated by electron microscopy , 2004, 2004 2nd IEEE International Symposium on Biomedical Imaging: Nano to Macro (IEEE Cat No. 04EX821).

[31]  E. Candès,et al.  Stable signal recovery from incomplete and inaccurate measurements , 2005, math/0503066.

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

[33]  Gabor T. Herman,et al.  A coordinate ascent approach to tomographic reconstruction of label images from a few projections , 2005, Discret. Appl. Math..

[34]  P. Cloetens,et al.  Advances in synchrotron radiation microtomography , 2006 .

[35]  C. Shaw,et al.  Noise simulation in cone beam CT imaging with parallel computing. , 2006, Physics in medicine and biology.

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

[37]  Attila Kuba,et al.  A Benchmark Evaluation of Large-Scale Optimization Approaches to Binary Tomography , 2006, DGCI.

[38]  Attila Kuba,et al.  Binary Tomography with Deblurring , 2006, IWCIA.

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

[40]  Kees Joost Batenburg.  A Network Flow Algorithm for Reconstructing Binary Images from Discrete X-rays , 2006, Journal of Mathematical Imaging and Vision.

[41]  Jiayu Song,et al.  Sparseness prior based iterative image reconstruction for retrospectively gated cardiac micro-CT. , 2007, Medical physics.

[42]  Florent Krzakala,et al.  Phase Transitions in the Coloring of Random Graphs , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[43]  G T Herman,et al.  Image reconstruction from a small number of projections , 2008, Inverse problems.

[44]  T. Gureyev,et al.  The binary dissector: phase contrast tomography of two- and three-material objects from few projections. , 2008, Optics express.

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

[46]  Gabor T. Herman,et al.  Fundamentals of Computerized Tomography: Image Reconstruction from Projections , 2009, Advances in Pattern Recognition.

[47]  David L. Donoho,et al.  Observed universality of phase transitions in high-dimensional geometry, with implications for modern data analysis and signal processing , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[48]  A. Mohammad-Djafari.  Gauss-Markov-Potts Priors for Images in Computer Tomography Resulting to Joint Optimal Reconstruction and segmentation , 2009 .

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

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

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

[52]  E. Maire,et al.  In Situ Experiments with X ray Tomography: an Attractive Tool for Experimental Mechanics , 2010 .

[53]  David L. Donoho,et al.  Precise Undersampling Theorems , 2010, Proceedings of the IEEE.

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

[55]  J. Pesquet,et al.  Parallel algorithm and hybrid regularization for dynamic PET reconstruction , 2010, IEEE Nuclear Science Symposuim & Medical Imaging Conference.

[56]  Kees Joost Batenburg,et al.  Bounds on the Difference between Reconstructions in Binary Tomography , 2011, DGCI.

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

[58]  Patrick L. Combettes,et al.  Proximal Splitting Methods in Signal Processing , 2009, Fixed-Point Algorithms for Inverse Problems in Science and Engineering.

[59]  Yannick Boursier,et al.  On the efficiency of proximal methods in CBCT and PET , 2011, 2011 18th IEEE International Conference on Image Processing.

[60]  Florent Krzakala,et al.  Statistical physics-based reconstruction in compressed sensing , 2011, ArXiv.

[61]  Nelly Pustelnik,et al.  Parallel Proximal Algorithm for Image Restoration Using Hybrid Regularization , 2009, IEEE Transactions on Image Processing.

[62]  Kees Joost Batenburg,et al.  A semi-automatic algorithm for grey level estimation in tomography , 2011, Pattern Recognit. Lett..

[63]  Florent Krzakala,et al.  Probabilistic reconstruction in compressed sensing: algorithms, phase diagrams, and threshold achieving matrices , 2012, ArXiv.

[64]  Lotfi Chaâri,et al.  Multidimensional Wavelet-based Regularized Reconstruction for Parallel Acquisition in Neuroimaging , 2012, ArXiv.

[65]  S. Roux,et al.  Tomographic reconstruction of binary fields , 2012 .

[66]  Mohamed-Jalal Fadili,et al.  A Generalized Forward-Backward Splitting , 2011, SIAM J. Imaging Sci..