Multi-source quantitative PAT in diffusive regime

Photoacoustic tomography (PAT) is a novel hybrid medical imaging technique that aims to combine the large contrast of optical coefficients with the high resolution capabilities of ultrasound. We assume that the first step of PAT, namely the reconstruction of a map of absorbed radiation from ultrasound boundary measurement, has been done. We focus on quantitative photoacoustic tomography (QPAT), which aims at quantitatively reconstructing the optical coefficients from knowledge of the absorbed radiation map. We present a non-iterative procedure to reconstruct such optical coefficients, namely the diffusion and absorption coefficients, and the Grüneisen coefficient when the propagation of radiation is modeled by a second-order elliptic equation. We show that PAT measurements allow us to uniquely reconstruct only two out of the above three coefficients, even when data are collected using an arbitrary number of radiation illuminations. We present uniqueness and stability results for the reconstructions of such two parameters and demonstrate the accuracy of the reconstruction algorithm with numerical reconstructions from two-dimensional synthetic data.

[1]  G. Richter.  An Inverse Problem for the Steady State Diffusion Equation , 1981 .

[2]  Nicolas Lerner,et al.  Uniqueness of continuous solutions for BV vector fields , 2002 .

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

[4]  Wotao Yin,et al.  An Iterative Regularization Method for Total Variation-Based Image Restoration , 2005, Multiscale Model. Simul..

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

[6]  L. Ambrosio.  Transport equation and Cauchy problem for BV vector fields , 2004 .

[7]  S. Arridge,et al.  Optical tomography: forward and inverse problems , 2009, 0907.2586.

[8]  Jérôme Fehrenbach,et al.  Imaging by Modification: Numerical Reconstruction of Local Conductivities from Corresponding Power Density Measurements , 2009, SIAM J. Imaging Sci..

[9]  S. Arridge,et al.  Estimating chromophore distributions from multiwavelength photoacoustic images. , 2009, Journal of the Optical Society of America. A, Optics, image science, and vision.

[10]  Lihong V. Wang.  Ultrasound-Mediated Biophotonic Imaging: A Review of Acousto-Optical Tomography and Photo-Acoustic Tomography , 2004, Disease markers.

[11]  Guillaume Bal,et al.  On multi-spectral quantitative photoacoustic tomography , 2011 .

[12]  Alexandru Tamasan,et al.  Conductivity imaging with a single measurement of boundary and interior data , 2007 .

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

[14]  Otmar Scherzer,et al.  Photoacoustic Imaging Taking into Account Attenuation , 2010, 1009.4350.

[15]  Joyce R. McLaughlin,et al.  Unique identifiability of elastic parameters from time-dependent interior displacement measurement , 2004 .

[16]  Lihong V. Wang,et al.  Photoacoustic imaging in biomedicine , 2006 .

[17]  M. Hauray.  On two-dimensional hamiltonian transport equations with Llocp coefficients , 2003, 1310.0974.

[18]  Alexandru Tamasan,et al.  Recovering the conductivity from a single measurement of interior data , 2009 .

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

[20]  Vasilis Ntziachristos,et al.  Quantitative point source photoacoustic inversion formulas for scattering and absorbing media. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  Guillaume Bal,et al.  Inverse transport theory and applications , 2009 .

[22]  P. Bassanini,et al.  Elliptic Partial Differential Equations of Second Order , 1997 .

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

[24]  Gianluca Crippa,et al.  Uniqueness, Renormalization, and Smooth Approximations for Linear Transport Equations , 2006, SIAM J. Math. Anal..

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

[26]  Faouzi Triki,et al.  Uniqueness and stability for the inverse medium problem with internal data , 2010 .

[27]  Guillaume Bal,et al.  Inverse diffusion theory of photoacoustics , 2009, 0910.2503.

[28]  Habib Ammari,et al.  Mathematical Modeling in Photoacoustic Imaging of Small Absorbers , 2010, SIAM Rev..

[29]  F. Colombini,et al.  A Note on Two-Dimensional Transport with Bounded Divergence , 2006 .

[30]  Stanley Osher,et al.  Bregman methods in quantitative photoacoustic tomography , 2010 .

[31]  Andrew R. Fisher,et al.  Photoacoustic effect for multiply scattered light. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[32]  Eric Bonnetier,et al.  Electrical Impedance Tomography by Elastic Deformation , 2008, SIAM J. Appl. Math..

[33]  B. T. Cox,et al.  The challenges for quantitative photoacoustic imaging , 2009, BiOS.

[34]  G. Alessandrini.  An identification problem for an elliptic equation in two variables , 1986 .

[35]  G. Bal,et al.  Inverse scattering and acousto-optic imaging. , 2009, Physical review letters.

[36]  Guillaume Bal,et al.  Quantitative thermo-acoustics and related problems , 2011 .

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

[38]  B. Cox,et al.  Quantitative Photoacoustic Image Reconstruction using Fluence Dependent Chromophores , 2010, Biomedical optics express.

[39]  P. Lions,et al.  On the Cauchy problem for Boltzmann equations: global existence and weak stability , 1989 .