Abstract:Photoacoustic Tomography (PAT) is an emerging biomedical "imaging from coupled physics" technique, in which the image contrast is due to optical absorption, but the information is carried to the surface of the tissue as ultrasound pulses. Many algorithms and formulae for PAT image reconstruction have been proposed for the case when a complete data set is available. In many practical imaging scenarios, however, it is not possible to obtain the full data, or the data may be sub-sampled for faster data acquisition. In such cases, image reconstruction algorithms that can incorporate prior knowledge to ameliorate the loss of data are required. Hence, recently there has been an increased interest in using variational image reconstruction. A crucial ingredient for the application of these techniques is the adjoint of the PAT forward operator, which is described in this article from physical, theoretical and numerical perspectives. First, a simple mathematical derivation of the adjoint of the PAT forward operator in the continuous framework is presented. Then, an efficient numerical implementation of the adjoint using a k-space time domain wave propagation model is described and illustrated in the context of variational PAT image reconstruction, on both 2D and 3D examples including inhomogeneous sound speed. The principal advantage of this analytical adjoint over an algebraic adjoint (obtained by taking the direct adjoint of the particular numerical forward scheme used) is that it can be implemented using currently available fast wave propagation solvers.

Variational Iterative Algorithms in Photoacoustic Tomography with Variable Sound Speed

Adaptive Imaging of Sound Source Based on Total Variation Prior and a Subspace Iteration Integrated Variational Bayesian Method

A Hierarchical Variational Bayesian Approximation Approach in Acoustic Imaging

On the Adjoint Operator in Photoacoustic Tomography

Efficient inclusion of total variation type priors in quantitative photoacoustic tomography

Study of Travel Time Tomography by Iterative Tikhonov Regularization Inversion

Implicit learning to determine variable sound speed and the reconstruction operator in photoacoustic tomography

Enhancement of Photoacoustic Tomography in the Tissue with Speed-of-sound Variance Using Ultrasound Computed Tomography

Inversions of two wave-forward operators with variable coefficients

Single-stage reconstruction algorithm for quantitative photoacoustic tomography

A Photoacoustic Imaging Reconstruction Method Based on Directional Total Variation with Adaptive Directivity

Total Variation-Based Limited-View Photoacoustic Reconstruction Method with Segmentation-Based Regularization

On Iterative Algorithms for Quantitative Photoacoustic Tomography in the Radiative Transport Regime

Determination of the sound speed and an initial source in photoacoustic tomography

Quantitative Photo-acoustic Tomography with Partial Data

A photoacoustic image reconstruction method using total variation and nonconvex optimization

Photoacoustic image reconstruction with an objective function using TGV and ESTGV as a regularization functional

Research on ADMM Reconstruction Algorithm of Photoacoustic Tomography With Limited Sampling Data

Asymptotic Expansions for Higher Order Elliptic Equations with an Application to Quantitative Photoacoustic Tomography

A high contrast and resolution reconstruction algorithm in quantitative photoacoustic tomography

Numerical Investigation of Nonlinear Sound Propagation of Photoacoustic Tomography Imaging.