A direct method for photoacoustic tomography with inhomogeneous sound speed

The standard approach for photoacoustic imaging with variable speed of sound is time reversal, which consists in solving a well-posed final-boundary value problem for the wave equation backwards in time. This paper investigates the iterative Landweber regularization algorithm, where convergence is guaranteed by standard regularization theory, notably also in cases of trapping sound speed or for short measurement times. We formulate and solve the direct and inverse problem on the whole Euclidean space, what is common in standard photoacoustic imaging, but not for time-reversal algorithms, where the problems are considered on a domain enclosed by the measurement devices. We formulate both the direct and adjoint photoacoustic operator as the solution of an interior and an exterior differential equation which are coupled by transmission conditions. The prior is solved numerically using a Galerkin scheme in space and finite difference discretization in time, while the latter consists in solving a boundary integral equation. We therefore use a BEM-FEM approach for numerical solution of the forward operators. We analyze this method, prove convergence, and provide numerical tests. Moreover, we compare the approach to time reversal.

[1]  Peter Kuchment,et al.  Mathematics of Photoacoustic and Thermoacoustic Tomography , 2009, Handbook of Mathematical Methods in Imaging.

[2]  P. Burgholzer,et al.  Photoacoustic tomography using a fiber based Fabry-Perot interferometer as an integrating line detector and image reconstruction by model-based time reversal method , 2007, European Conference on Biomedical Optics.

[3]  M. Fink,et al.  Time reversal of ultrasonic fields. I. Basic principles , 1992, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[4]  C. Lubich Convolution quadrature and discretized operational calculus. I , 1988 .

[5]  Giovanni Monegato,et al.  A space–time BIE method for nonhomogeneous exterior wave equation problems. The Dirichlet case , 2012 .

[6]  M. Fink,et al.  Time-reversal of ultrasonic fields. III. Theory of the closed time-reversal cavity , 1992, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[7]  T. Ha-Duong,et al.  On Retarded Potential Boundary Integral Equations and their Discretisation , 2003 .

[8]  Chen Zhang,et al.  Total variation based gradient descent algorithm for sparse-view photoacoustic image reconstruction. , 2012, Ultrasonics.

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

[10]  Leonid Kunyansky,et al.  Reconstruction of a function from its spherical (circular) means with the centers lying on the surface of certain polygons and polyhedra , 2010, 1009.0288.

[11]  M. Hanke Conjugate gradient type methods for ill-posed problems , 1995 .

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

[13]  Lihong V. Wang,et al.  Dark-Field Confocal Photoacoustic Microscopy , 2009 .

[14]  C. Bardos,et al.  Sharp sufficient conditions for the observation, control, and stabilization of waves from the boundary , 1992 .

[15]  A. Kirsch An Introduction to the Mathematical Theory of Inverse Problems , 1996, Applied Mathematical Sciences.

[16]  Minghua Xu,et al.  Time-domain reconstruction for thermoacoustic tomography in a spherical geometry , 2002, IEEE Transactions on Medical Imaging.

[17]  O. Scherzer Handbook of mathematical methods in imaging , 2011 .

[18]  D. L. Russell Review: J.-L. Lions, Controlabilité Exacte, Perturbations et Stabilisation de Systèmes Distribués , 1990 .

[19]  Linh V. Nguyen,et al.  Reconstruction and time reversal in thermoacoustic tomography in acoustically homogeneous and inhomogeneous media , 2008 .

[20]  F. G. Friedlander,et al.  The Wave Equation on a Curved Space-Time (Book Review) , 1976 .

[21]  P. Davies,et al.  The Wave Equation on a Curved Space–Time , 1976 .

[22]  Hongkai Zhao,et al.  An Efficient Neumann Series-Based Algorithm for Thermoacoustic and Photoacoustic Tomography with Variable Sound Speed , 2011, SIAM J. Imaging Sci..

[23]  Yuan Xu,et al.  Exact frequency-domain reconstruction for thermoacoustic tomography. I. Planar geometry , 2002, IEEE Transactions on Medical Imaging.

[24]  Giovanni Monegato,et al.  Exact nonreflecting boundary conditions for exterior wave equation problems , 2014 .

[25]  Minghua Xu,et al.  Time-domain reconstruction algorithms and numerical simulations for thermoacoustic tomography in various geometries , 2003, IEEE Transactions on Biomedical Engineering.

[26]  Peter Kuchment,et al.  Uniqueness of reconstruction and an inversion procedure for thermoacoustic and photoacoustic tomography with variable sound speed , 2007, 0706.0598.

[27]  Yiqiu Dong,et al.  An algorithm for total variation regularized photoacoustic imaging , 2015, Adv. Comput. Math..

[28]  Jean E. Roberts,et al.  Higher Order Triangular Finite Elements with Mass Lumping for the Wave Equation , 2000, SIAM J. Numer. Anal..

[29]  Yulia Hristova,et al.  Time reversal in thermoacoustic tomography—an error estimate , 2008, 0812.0606.

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

[31]  Dimple Modgil,et al.  Image reconstruction in photoacoustic tomography with variable speed of sound using a higher-order geometrical acoustics approximation. , 2010, Journal of biomedical optics.

[32]  H. Engl,et al.  Regularization of Inverse Problems , 1996 .

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

[34]  J. Lions,et al.  Non-homogeneous boundary value problems and applications , 1972 .

[35]  C. Lubich Convolution quadrature and discretized operational calculus. II , 1988 .

[36]  Lihong V. Wang Photoacoustic imaging and spectroscopy , 2009 .

[37]  Peter Monk,et al.  A new algorithm in electromagnetic inverse scattering theory with an application to medical imaging , 1997 .

[38]  Minghua Xu,et al.  Exact frequency-domain reconstruction for thermoacoustic tomography. II. Cylindrical geometry , 2002, IEEE Transactions on Medical Imaging.

[39]  Patrick Joly,et al.  Coupling discontinuous Galerkin methods and retarded potentials for transient wave propagation on unbounded domains , 2011, J. Comput. Phys..

[40]  Plamen Stefanov,et al.  Thermoacoustic tomography arising in brain imaging , 2010, 1009.1687.

[41]  B T Cox,et al.  k-Wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields. , 2010, Journal of biomedical optics.

[42]  Bernardo Cockburn Discontinuous Galerkin methods , 2003 .

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

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

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

[46]  C. W. Groetsch,et al.  The theory of Tikhonov regularization for Fredholm equations of the first kind , 1984 .

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

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

[49]  V. Komornik Exact Controllability and Stabilization: The Multiplier Method , 1995 .

[50]  Lihong V. Wang,et al.  Prospects of photoacoustic tomography. , 2008, Medical physics.