A New Numerical Algorithm for Thermoacoustic and Photoacoustic Tomography with Variable Sound Speed

We present a new algorithm for reconstructing an unknown source in Thermoacoustic and Photoacoustic Tomography based on the recent advances in understanding the theoretical nature of the problem. We work with variable sound speeds that might be also discontinuous across some surface. The latter problem arises in brain imaging. The new algorithm is based on an explicit formula in the form of a Neumann series. We present numerical examples with nontrapping, trapping and piecewise smooth speeds, as well as examples with data on a part of the boundary. These numerical examples demonstrate the robust performance of the new algorithm.

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

[2]  Lihong V. Wang,et al.  Biomedical Optics: Principles and Imaging , 2007 .

[3]  Robert A Kruger,et al.  Thermoacoustic computed tomography using a conventional linear transducer array. , 2003, Medical physics.

[4]  Yuan Xu,et al.  Effects of acoustic heterogeneity in breast thermoacoustic tomography , 2003, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[5]  Lihong V. Wang,et al.  Thermoacoustic tomography with correction for acoustic speed variations , 2006, Physics in medicine and biology.

[6]  D. Medková.  On the convergence of Neumann series for noncompact operators , 1991 .

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

[8]  Minghua Xu,et al.  Thermoacoustic and Photoacoustic Tomography of Thick Biological Tissues toward Breast Imaging , 2005, Technology in cancer research & treatment.

[9]  S. Osher,et al.  Lax-Friedrichs sweeping scheme for static Hamilton-Jacobi equations , 2004 .

[10]  L. Hörmander.  Fourier integral operators. I , 1995 .

[11]  Lihong V. Wang,et al.  Universal back-projection algorithm for photoacoustic computed tomography. , 2005 .

[12]  Xinmai Yang,et al.  Monkey brain cortex imaging by photoacoustic tomography. , 2008, Journal of biomedical optics.

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

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

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

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

[17]  Qing Huo Liu,et al.  The perfectly matched layer for acoustic waves in absorptive media , 1997 .

[18]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[19]  Maciej Zworski,et al.  Complex scaling and the distribution of scattering poles , 1991 .

[20]  Hongkai Zhao,et al.  Fast Sweeping Methods for Eikonal Equations on Triangular Meshes , 2007, SIAM J. Numer. Anal..

[21]  Otmar Scherzer,et al.  Thermoacoustic computed tomography with large planar receivers , 2004 .

[22]  Lihong V. Wang,et al.  Time reversal and its application to tomography with diffracting sources. , 2004, Physical review letters.

[23]  Wei Wang,et al.  Simultaneous Molecular and Hypoxia Imaging of Brain Tumors In Vivo Using Spectroscopic Photoacoustic Tomography , 2008, Proceedings of the IEEE.

[24]  Hongkai Zhao,et al.  A fast sweeping method for Eikonal equations , 2004, Math. Comput..

[25]  S. Patch,et al.  Thermoacoustic tomography--consistency conditions and the partial scan problem. , 2004, Physics in medicine and biology.

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

[27]  R A Kruger,et al.  Thermoacoustic computed tomography--technical considerations. , 1999, Medical physics.

[28]  Daniel Tataru,et al.  Unique continuation for operatorswith partially analytic coefficients , 1999 .

[29]  L. Hörmander,et al.  Fourier integral operators. II , 1972 .

[30]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .

[31]  Yuan Xu,et al.  Rhesus monkey brain imaging through intact skull with thermoacoustic tomography , 2006, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

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

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

[34]  Plamen Stefanov,et al.  Linearizing non-linear inverse problems and an application to inverse backscattering , 2008, 0809.0270.

[35]  D. Brandt,et al.  Multi-level adaptive solutions to boundary-value problems math comptr , 1977 .