Single-stage reconstruction algorithm for quantitative photoacoustic tomography

The development of efficient and accurate image reconstruction algorithms is one of the cornerstones of computed tomography. Existing algorithms for quantitative photoacoustic tomography currently operate in a two-stage procedure: First an inverse source problem for the acoustic wave propagation is solved, whereas in a second step the optical parameters are estimated from the result of the first step. Such an approach has several drawbacks. In this paper we therefore propose the use of single-stage reconstruction algorithms for quantitative photoacoustic tomography, where the optical parameters are directly reconstructed from the observed acoustical data. In that context we formulate the image reconstruction problem of quantitative photoacoustic tomography as a single nonlinear inverse problem by coupling the radiative transfer equation with the acoustic wave equation. The inverse problem is approached by Tikhonov regularization with a convex penalty in combination with the proximal gradient iteration for minimizing the Tikhonov functional. We present numerical results, where the proposed single-stage algorithm shows an improved reconstruction quality at a similar computational cost.

[1]  Paul C. Beard,et al.  Gradient-based quantitative photoacoustic image reconstruction for molecular imaging , 2007, SPIE BiOS.

[2]  Rakesh,et al.  The spherical mean value operator with centers on a sphere , 2007 .

[3]  P. Kuchment,et al.  On reconstruction formulas and algorithms for the thermoacoustic tomography , 2007, 0706.1303.

[4]  Markus Haltmeier,et al.  Universal Inversion Formulas for Recovering a Function from Spherical Means , 2012, SIAM J. Math. Anal..

[5]  Simon R. Arridge,et al.  A gradient-based method for quantitative photoacoustic tomography using the radiative transfer equation , 2013 .

[6]  M. Haltmeier,et al.  Temporal back-projection algorithms for photoacoustic tomography with integrating line detectors , 2007 .

[7]  L. Kunyansky,et al.  Explicit inversion formulae for the spherical mean Radon transform , 2006, math/0609341.

[8]  V. Ntziachristos,et al.  Acoustic Inversion in Optoacoustic Tomography: A Review , 2013, Current medical imaging reviews.

[9]  Dario Fasino,et al.  An inverse Robin problem for Laplace's equation: theoretical results and numerical methods , 1999 .

[10]  Rakesh,et al.  Determining a Function from Its Mean Values Over a Family of Spheres , 2004, SIAM J. Math. Anal..

[11]  C. DeWitt-Morette,et al.  Mathematical Analysis and Numerical Methods for Science and Technology , 1990 .

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

[13]  Anabela Da Silva,et al.  Considering sources and detectors distributions for quantitative photoacoustic tomography. , 2014, Biomedical optics express.

[14]  Émilie Chouzenoux,et al.  Variable Metric Forward–Backward Algorithm for Minimizing the Sum of a Differentiable Function and a Convex Function , 2013, Journal of Optimization Theory and Applications.

[15]  S. Arridge.  Optical tomography in medical imaging , 1999 .

[16]  H. Egger,et al.  An Lp theory for stationary radiative transfer , 2013, 1304.6504.

[17]  Stefan Kunis,et al.  Effective Discretization of Direct Reconstruction Schemes For Photoacoustic Imaging In Spherical Geometries , 2014, SIAM J. Numer. Anal..

[18]  Kui Ren,et al.  A one-step reconstruction algorithm for quantitative photoacoustic imaging , 2015, 1507.02776.

[19]  Markus Haltmeier,et al.  Inversion of circular means and the wave equation on convex planar domains , 2012, Comput. Math. Appl..

[20]  Eric Todd Quinto,et al.  Artifacts in Incomplete Data Tomography with Applications to Photoacoustic Tomography and Sonar , 2014, SIAM J. Appl. Math..

[21]  Simon R. Arridge,et al.  Reconstructing absorption and scattering distributions in quantitative photoacoustic tomography , 2012 .

[22]  Peter Kuchment,et al.  Mathematics of thermoacoustic tomography , 2007, European Journal of Applied Mathematics.

[23]  R. Kruger,et al.  Photoacoustic ultrasound (PAUS)--reconstruction tomography. , 1995, Medical physics.

[24]  Patrick L. Combettes,et al.  Signal Recovery by Proximal Forward-Backward Splitting , 2005, Multiscale Model. Simul..

[25]  Benar Fux Svaiter,et al.  Convergence of descent methods for semi-algebraic and tame problems: proximal algorithms, forward–backward splitting, and regularized Gauss–Seidel methods , 2013, Math. Program..

[26]  Guillaume Bal,et al.  Inverse transport theory of photoacoustics , 2009, 0908.4012.

[27]  Simon R Arridge,et al.  Two-dimensional quantitative photoacoustic image reconstruction of absorption distributions in scattering media by use of a simple iterative method. , 2006, Applied optics.

[28]  Herbert Egger,et al.  Stationary radiative transfer with vanishing absorption , 2014 .

[29]  Otmar Scherzer,et al.  Quantitative Photoacoustic Tomography with Piecewise Constant Material Parameters , 2014, SIAM J. Imaging Sci..

[30]  Herbert Egger,et al.  Numerical methods for parameter identification in stationary radiative transfer , 2013, Comput. Optim. Appl..

[31]  Kui Ren,et al.  Quantitative photoacoustic imaging in the radiative transport regime , 2012, 1207.4664.

[32]  Otmar Scherzer,et al.  Filtered backprojection for thermoacoustic computed tomography in spherical geometry , 2005, Mathematical Methods in the Applied Sciences.

[33]  P. Beard.  Biomedical photoacoustic imaging , 2011, Interface Focus.

[34]  Vasilis Ntziachristos,et al.  Fast Semi-Analytical Model-Based Acoustic Inversion for Quantitative Optoacoustic Tomography , 2010, IEEE Transactions on Medical Imaging.

[35]  Victor Palamodov.  Remarks on the generalFunk transform and thermoacoustic tomography , 2010 .

[36]  J. Craggs.  Applied Mathematical Sciences , 1973 .

[37]  Markus Haltmeier,et al.  Inversion of Spherical Means and the Wave Equation in Even Dimensions , 2007, SIAM J. Appl. Math..

[38]  Huabei Jiang,et al.  Transport-based quantitative photoacoustic tomography: simulations and experiments , 2010, Physics in medicine and biology.

[39]  A. Aisen,et al.  Thermoacoustic CT with radio waves: a medical imaging paradigm. , 1999, Radiology.

[40]  Simon R. Arridge,et al.  Multiple Illumination Quantitative Photoacoustic Tomography using Transport and Diffusion Models , 2011 .

[41]  Victor Palamodov.  Remarks on the general Funk-Radon transform and thermoacoustic tomography , 2007 .

[42]  Marc Teboulle,et al.  Proximal alternating linearized minimization for nonconvex and nonsmooth problems , 2013, Mathematical Programming.

[43]  Markus Haltmeier,et al.  Experimental evaluation of reconstruction algorithms for limited view photoacoustic tomography with line detectors , 2007 .

[44]  Peter Kuchment,et al.  Mathematics of thermoacoustic and photoacoustic tomography , 2007 .

[45]  Habib Ammari,et al.  Reconstruction of the Optical Absorption Coefficient of a Small Absorber from the Absorbed Energy Density , 2011, SIAM J. Appl. Math..

[46]  W. D. Evans,et al.  PARTIAL DIFFERENTIAL EQUATIONS , 1941 .

[47]  Lihong V. Wang,et al.  Reconstructions in limited-view thermoacoustic tomography. , 2004, Medical physics.

[48]  A. Cezaro,et al.  Regularization approaches for quantitative photoacoustic tomography using the radiative transfer equation , 2013, 1307.3201.

[49]  Xiaoqun Zhang,et al.  Forward–backward splitting method for quantitative photoacoustic tomography , 2014 .

[50]  Guillaume Bal,et al.  Multi-source quantitative photoacoustic tomography in a diffusive regime , 2011 .

[51]  Roger J Zemp.  Quantitative photoacoustic tomography with multiple optical sources. , 2010, Applied optics.

[52]  S. Arridge,et al.  Quantitative spectroscopic photoacoustic imaging: a review. , 2012, Journal of biomedical optics.

[53]  Barbara Kaltenbacher,et al.  Regularization Methods in Banach Spaces , 2012, Radon Series on Computational and Applied Mathematics.

[54]  Jie Chen,et al.  Quantitative photo-acoustic tomography with partial data , 2012, 1204.2213.

[55]  G. Uhlmann,et al.  Thermoacoustic tomography with variable sound speed , 2009, 0902.1973.

[56]  Otmar Scherzer,et al.  Variational Methods in Imaging , 2008, Applied mathematical sciences.

[57]  Lihong V. Wang.  Multiscale photoacoustic microscopy and computed tomography. , 2009, Nature photonics.

[58]  V. Morozov.  The Regularization Method , 1984 .

[59]  Hao Gao,et al.  A Hybrid Reconstruction Method for Quantitative PAT , 2013, SIAM J. Imaging Sci..

[60]  Mustapha Mokhtar Kharroubi.  Mathematical Topics in Neutron Transport Theory: New Aspects , 1997 .

[61]  Guido Kanschat,et al.  Solution of Radiative Transfer Problems with Finite Elements , 2009 .

[62]  Xu Xiao.  Photoacoustic imaging in biomedicine , 2008 .

[63]  Richard Kowar,et al.  On Time Reversal in Photoacoustic Tomography for Tissue Similar to Water , 2013, SIAM J. Imaging Sci..

[64]  Mark A. Anastasio,et al.  Photoacoustic and Thermoacoustic Tomography: Image Formation Principles , 2015, Handbook of Mathematical Methods in Imaging.