Gradient-flow-based semi-implicit finite-element method and its convergence analysis for image reconstruction

In this paper, we present a novel and effective L 2 -gradient-flow-based semiimplicit finite-element method for solving a variational problem of image reconstruction. The method is applicable to several data scenarios, especially for the contaminated data detected from uniformly sparse or randomly distributed projection directions. We also give a complete and rigorous proof for the convergence of the semi-implicit finite-element method, in which the convergence does not rely on the choices of the regularization parameter and the temporal step size. The experimental results show that our method has more desirable performance comparing with other reconstruction methods in solving a number of challenging reconstruction problems. (Some figures may appear in colour only in the online journal)

[1]  Constantin Popa,et al.  Constrained Kaczmarz extended algorithm for image reconstruction , 2008 .

[2]  Ming Li,et al.  Computational Inversion of Electron Tomography Images Using L2-Gradient Flows. , 2011, Journal of computational mathematics : an international journal on numerical methods, analysis and applications.

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

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

[5]  Chandrajit L. Bajaj,et al.  INVERSION OF ELECTRON TOMOGRAPHY IMAGES USING L 2 -GRADIENT FLOWS —- COMPUTATIONAL METHODS * , 2011 .

[6]  Jacques Simeon,et al.  Compact Sets in the Space L~(O, , 2005 .

[7]  A. C. Riddle,et al.  Inversion of Fan-Beam Scans in Radio Astronomy , 1967 .

[8]  Constantin Popa.  A hybrid Kaczmarz-Conjugate Gradient algorithm for image reconstruction , 2010, Math. Comput. Simul..

[9]  Jan Timmer,et al.  The gridding method for image reconstruction by Fourier transformation , 1995, IEEE Trans. Medical Imaging.

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

[11]  D. J. De Rosier,et al.  Reconstruction of Three Dimensional Structures from Electron Micrographs , 1968, Nature.

[12]  P. Gilbert.  Iterative methods for the three-dimensional reconstruction of an object from projections. , 1972, Journal of theoretical biology.

[13]  Sigurdur Helgason,et al.  The Radon Transform on ℝn , 1999 .

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

[15]  Yoram Bresler,et al.  A fast and accurate Fourier algorithm for iterative parallel-beam tomography , 1996, IEEE Trans. Image Process..

[16]  L. Vese.  A Study in the BV Space of a Denoising—Deblurring Variational Problem , 2001 .

[17]  L. Evans.  Measure theory and fine properties of functions , 1992 .

[18]  Michel Barlaud,et al.  Deterministic edge-preserving regularization in computed imaging , 1997, IEEE Trans. Image Process..

[19]  Guoliang Xu,et al.  Single-particle reconstruction using L(2)-gradient flow. , 2011, Journal of structural biology.

[20]  Constantin Popa.  Extended and constrained diagonal weighting algorithm with application to inverse problems in image reconstruction , 2010 .

[21]  Frank Natterer,et al.  Fourier reconstruction in tomography , 1985 .

[22]  Gjlles Aubert,et al.  Mathematical problems in image processing , 2001 .

[23]  Guoliang Xu,et al.  A new reconstruction algorithm in tomography with geometric feature-preserving regularization , 2010, 2010 3rd International Conference on Biomedical Engineering and Informatics.

[24]  A. Kak,et al.  Simultaneous Algebraic Reconstruction Technique (SART): A Superior Implementation of the Art Algorithm , 1984, Ultrasonic imaging.

[25]  Pierre Kornprobst,et al.  Mathematical problems in image processing - partial differential equations and the calculus of variations , 2010, Applied mathematical sciences.

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

[27]  D. DeRosier,et al.  The reconstruction of a three-dimensional structure from projections and its application to electron microscopy , 1970, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[28]  Adel Faridani,et al.  Introduction to the Mathematics of Computed Tomography , 2003 .

[29]  S. Helgason.  The Radon Transform , 1980 .

[30]  Xiaobing Feng,et al.  Analysis of total variation flow and its finite element approximations , 2003 .

[31]  Gabriele Steidl,et al.  A new linogram algorithm for computerized tomography , 2001 .

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

[33]  F. Natterer.  The Mathematics of Computerized Tomography , 1986 .

[34]  Gabor T. Herman,et al.  Image reconstruction from projections : the fundamentals of computerized tomography , 1980 .

[35]  J. Simon.  Compact sets in the spaceLp(O,T; B) , 1986 .

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

[37]  R. Temam,et al.  Pseudosolutions of the time-dependent minimal surface problem , 1978 .

[38]  R. Bracewell.  Strip Integration in Radio Astronomy , 1956 .

[39]  G. Herman,et al.  Image reconstruction from linograms: implementation and evaluation. , 1988, IEEE transactions on medical imaging.

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

[41]  J. Frank.  Three-Dimensional Electron Microscopy of Macromolecular Assemblies , 2006 .