Full-Wave Iterative Image Reconstruction in Photoacoustic Tomography With Acoustically Inhomogeneous Media

Existing approaches to image reconstruction in photoacoustic computed tomography (PACT) with acoustically heterogeneous media are limited to weakly varying media, are computationally burdensome, and/or cannot effectively mitigate the effects of measurement data incompleteness and noise. In this work, we develop and investigate a discrete imaging model for PACT that is based on the exact photoacoustic (PA) wave equation and facilitates the circumvention of these limitations. A key contribution of the work is the establishment of a procedure to implement a matched forward and backprojection operator pair associated with the discrete imaging model, which permits application of a wide-range of modern image reconstruction algorithms that can mitigate the effects of data incompleteness and noise. The forward and backprojection operators are based on the k-space pseudospectral method for computing numerical solutions to the PA wave equation in the time domain. The developed reconstruction methodology is investigated by use of both computer-simulated and experimental PACT measurement data.

[1]  Lihong V. Wang,et al.  Photoacoustic computed tomography correcting for heterogeneity and attenuation. , 2012, Journal of biomedical optics.

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

[3]  Thomas L. Szabo,et al.  Time domain wave equations for lossy media obeying a frequency power law , 1994 .

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

[5]  Gerald R. Harris,et al.  Review of transient field theory for a baffled planar piston , 1981 .

[6]  T. D. Mast,et al.  A k-space method for coupled first-order acoustic propagation equations. , 2002, The Journal of the Acoustical Society of America.

[7]  Huabei Jiang,et al.  Three-dimensional finite-element-based photoacoustic tomography: reconstruction algorithm and simulations. , 2007, Medical physics.

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

[9]  Wiendelt Steenbergen,et al.  Passive element enriched photoacoustic computed tomography (PER PACT) for simultaneous imaging of acoustic propagation properties and light absorption. , 2011, Optics express.

[10]  M Tanter,et al.  Experimental demonstration of noninvasive transskull adaptive focusing based on prior computed tomography scans. , 2003, The Journal of the Acoustical Society of America.

[11]  J.A. Fessler,et al.  Penalized weighted least-squares image reconstruction for positron emission tomography , 1994, 5th IEEE EMBS International Summer School on Biomedical Imaging, 2002..

[12]  Thomas L. Szabo,et al.  Diagnostic Ultrasound Imaging: Inside Out , 2004 .

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

[14]  Mark A. Anastasio,et al.  An Imaging Model Incorporating Ultrasonic Transducer Properties for Three-Dimensional Optoacoustic Tomography , 2011, IEEE Transactions on Medical Imaging.

[15]  Lihong V. Wang,et al.  Photoacoustic tomography of water in phantoms and tissue. , 2010, Journal of biomedical optics.

[16]  Quing Zhu,et al.  In vivo photoacoustic tomography of mouse cerebral edema induced by cold injury. , 2011, Journal of biomedical optics.

[17]  Liming Nie,et al.  Photoacoustic tomography of monkey brain using virtual point ultrasonic transducers. , 2011, Journal of biomedical optics.

[18]  Huabei Jiang,et al.  Enhancing finite element-based photoacoustic tomography using total variation minimization , 2011 .

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

[20]  Vasilis Ntziachristos,et al.  The effects of acoustic attenuation in optoacoustic signals , 2011, Physics in medicine and biology.

[21]  Vasilis Ntziachristos,et al.  Statistical optoacoustic image reconstruction using a-priori knowledge on the location of acoustic distortions , 2011 .

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

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

[24]  C. Antonopoulos,et al.  A general form of perfectly matched layers for for three-dimensional problems of acoustic scattering in lossless and lossy fluid media , 2004, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[25]  M. Anastasio,et al.  Investigation of iterative image reconstruction in three-dimensional optoacoustic tomography , 2012, Physics in medicine and biology.

[26]  Mark A. Anastasio,et al.  Investigation of limited-view image reconstruction in optoacoustic tomography employing a priori structural information , 2010, Optical Engineering + Applications.

[27]  Lihong V. Wang,et al.  Effects of acoustic heterogeneities on transcranial brain imaging with microwave-induced thermoacoustic tomography. , 2008, Medical physics.

[28]  Patrick J. La Riviere,et al.  Photoacoustic image reconstruction in an attenuating medium using singular-value decomposition , 2008, 2008 IEEE Nuclear Science Symposium Conference Record.

[29]  T. D. Mast,et al.  A k-space method for large-scale models of wave propagation in tissue , 2001, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

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

[31]  D. T. Lee,et al.  Two algorithms for constructing a Delaunay triangulation , 1980, International Journal of Computer & Information Sciences.

[32]  B. Cox,et al.  Photoacoustic tomography in absorbing acoustic media using time reversal , 2010 .

[33]  Marc Teboulle,et al.  Fast Gradient-Based Algorithms for Constrained Total Variation Image Denoising and Deblurring Problems , 2009, IEEE Transactions on Image Processing.

[34]  A. Perretti,et al.  Diagnostic ultrasound imaging. , 1990, Rays.

[35]  Alexander A. Oraevsky,et al.  Detection of optoacoustic transients with a rectangular transducer of finite dimensions , 2002, SPIE BiOS.

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

[37]  Wiendelt Steenbergen,et al.  Speed-of-sound compensated photoacoustic tomography for accurate imaging. , 2012, Medical physics.

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

[39]  Robert Nuster,et al.  Compensation of acoustic attenuation for high-resolution photoacoustic imaging with line detectors , 2007, SPIE BiOS.

[40]  J. Barger,et al.  Acoustical properties of the human skull. , 1978, The Journal of the Acoustical Society of America.

[41]  Shuai Xu,et al.  Accelerating MatLab code using GPU: A review of tools and strategies , 2011, 2011 2nd International Conference on Artificial Intelligence, Management Science and Electronic Commerce (AIMSEC).

[42]  Patrick J La Rivière,et al.  Image reconstruction in optoacoustic tomography for dispersive acoustic media. , 2006, Optics letters.

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

[44]  E. Somersalo,et al.  Statistical inverse problems: discretization, model reduction and inverse crimes , 2007 .

[45]  G. Farhat,et al.  Diagnostic ultrasound Imaging : Inside out , 2004 .

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

[47]  Lihong V. Wang,et al.  Numerical investigation of the effects of shear waves in transcranial photoacoustic tomography with a planar geometry. , 2012, Journal of biomedical optics.

[48]  Chao Huang,et al.  Aberration correction for transcranial photoacoustic tomography of primates employing adjunct image data. , 2012, Journal of biomedical optics.

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

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