Estimating chromophore distributions from multiwavelength photoacoustic images.

Biomedical photoacoustic tomography (PAT) can provide qualitative images of biomedical soft tissue with high spatial resolution. However, whether it is possible to give accurate quantitative estimates of the spatially varying concentrations of the sources of photoacoustic contrast-endogenous or exogenous chromophores-remains an open question. Even if the chromophores' absorption spectra are known, the problem is nonlinear and ill-posed. We describe a framework for obtaining such quantitative estimates. When the optical scattering distribution is known, adjoint and gradient-based optimization techniques can be used to recover the concentration distributions of the individual chromophores that contribute to the overall tissue absorption. When the scattering distribution is unknown, prior knowledge of the wavelength dependence of the scattering is shown to be sufficient to overcome the absorption-scattering nonuniqueness and allow both distributions of chromophore concentrations and scattering to be recovered from multiwavelength photoacoustic images.

[1]  Zhuang Liu,et al.  Carbon nanotubes as photoacoustic molecular imaging agents in living mice. , 2008, Nature nanotechnology.

[2]  Lihong V. Wang,et al.  Monte Carlo Modeling of Light Transport in Tissues , 1995 .

[3]  William H. Press,et al.  Numerical recipes in C. The art of scientific computing , 1987 .

[4]  L Wang,et al.  MCML--Monte Carlo modeling of light transport in multi-layered tissues. , 1995, Computer methods and programs in biomedicine.

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

[6]  Lihong V. Wang,et al.  Noninvasive laser-induced photoacoustic tomography for structural and functional in vivo imaging of the brain , 2003, Nature Biotechnology.

[7]  Per Christian Hansen,et al.  Deconvolution and Regularization with Toeplitz Matrices , 2002, Numerical Algorithms.

[8]  Bo Liu,et al.  Photoacoustic spectroscopic imaging of intra-tumor heterogeneity and molecular identification , 2006, SPIE BiOS.

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

[10]  Qiang Wang,et al.  Reconstruction of optical absorption coefficient maps of heterogeneous media by photoacoustic tomography coupled with diffusion equation based regularized Newton method. , 2007, Optics express.

[11]  Paul C. Beard,et al.  Quantitative photoacoustic image reconstruction for molecular imaging , 2006, SPIE BiOS.

[12]  Lihong V. Wang,et al.  Technical considerations in quantitative blood oxygenation measurement using photoacoustic microscopy in vivo , 2006, SPIE BiOS.

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

[14]  P. C. Beard,et al.  Quantitative photoacoustic imaging: fitting a model of light transport to the initial pressure distribution , 2005, SPIE BiOS.

[15]  A. N. Bashkatov,et al.  Optical properties of human skin, subcutaneous and mucous tissues in the wavelength range from 400 to 2000 nm , 2005 .

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

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

[18]  S. Thennadil,et al.  Optical properties of human skin in the near infrared wavelength range of 1000 to 2200 nm. , 2001, Journal of biomedical optics.

[19]  B T Cox,et al.  Fast calculation of pulsed photoacoustic fields in fluids using k-space methods. , 2005, The Journal of the Acoustical Society of America.

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

[21]  M. Schweiger,et al.  Diffuse optical tomography with spectral constraints and wavelength optimization. , 2005, Applied optics.

[22]  L. O. Svaasand,et al.  Collimated light sources in the diffusion approximation. , 2000, Applied optics.

[23]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[24]  Jan Laufer,et al.  Quantitative spatially resolved measurement of tissue chromophore concentrations using photoacoustic spectroscopy: application to the measurement of blood oxygenation and haemoglobin concentration , 2007, Physics in medicine and biology.

[25]  Qiang Wang,et al.  Quantitative photoacoustic tomography: recovery of optical absorption coefficient maps of heterogeneous media , 2007, SPIE BiOS.

[26]  R. Cubeddu,et al.  Optical Tomography , 1998, Technical Digest. 1998 EQEC. European Quantum Electronics Conference (Cat. No.98TH8326).

[27]  M. Schweiger,et al.  A finite element approach for modeling photon transport in tissue. , 1993, Medical physics.

[28]  Lihong V. Wang,et al.  Effects of wavelength-dependent fluence attenuation on the noninvasive photoacoustic imaging of hemoglobin oxygen saturation in subcutaneous vasculature in vivo , 2007, SPIE BiOS.

[29]  Ashleyj . Welch,et al.  Optical-Thermal Response of Laser-Irradiated Tissue , 1995 .

[30]  Jan Laufer,et al.  Backward-mode multiwavelength photoacoustic scanner using a planar Fabry-Perot polymer film ultrasound sensor for high-resolution three-dimensional imaging of biological tissues. , 2008, Applied optics.

[31]  Geng Ku,et al.  Noninvasive imaging of hemoglobin concentration and oxygenation in the rat brain using high-resolution photoacoustic tomography. , 2006, Journal of biomedical optics.

[32]  Paul C. Beard,et al.  2D Quantitative photoacoustic image reconstruction of absorption distributions in scattering media using a simple recursive method , 2006 .

[33]  Qiang Wang,et al.  Tomographic imaging of absolute optical absorption coefficient in turbid media using combined photoacoustic and diffusing light measurements. , 2007, Optics letters.

[34]  B. Hooper.  Optical-thermal response of laser-irradiated tissue , 1996 .

[35]  Peter Kuchment,et al.  On reconstruction formulas and algorithms for the thermoacoustic and photoacoustic tomography , 2007 .

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

[37]  F. D. de Mul,et al.  Three-dimensional photoacoustic imaging of blood vessels in tissue. , 1998, Optics letters.

[38]  Anthony J. Durkin,et al.  In vivo determination of skin near-infrared optical properties using diffuse optical spectroscopy. , 2008, Journal of biomedical optics.

[39]  Paul C. Beard,et al.  Simultaneous estimation of chromophore concentration and scattering distributions from multiwavelength photoacoustic images , 2008, SPIE BiOS.

[40]  M. Schweiger,et al.  The finite element method for the propagation of light in scattering media: boundary and source conditions. , 1995, Medical physics.

[41]  L. Trefethen,et al.  Numerical linear algebra , 1997 .

[42]  Debasish Roy,et al.  Quantitative photoacoustic tomography from boundary pressure measurements: noniterative recovery of optical absorption coefficient from the reconstructed absorbed energy map. , 2008, Journal of the Optical Society of America. A, Optics, image science, and vision.

[43]  B T Cox,et al.  k-space propagation models for acoustically heterogeneous media: application to biomedical photoacoustics. , 2007, The Journal of the Acoustical Society of America.

[44]  Jan Laufer,et al.  In vitro measurements of absolute blood oxygen saturation using pulsed near-infrared photoacoustic spectroscopy: accuracy and resolution. , 2005, Physics in medicine and biology.

[45]  Lihong V. Wang,et al.  Limitations of quantitative photoacoustic measurements of blood oxygenation in small vessels , 2007, Physics in medicine and biology.