Direct quantitative photoacoustic tomography for realistic acoustic media

DOI10.1088/1361-6420/AB091EzbMATH Open1427.65224arXiv1810.11685OpenAlexW3105699515MaRDI QIDQ5228008FDOQ5228008

Authors: Ashkan Javaherian, Sean F. Holman 

Publication date: 8 August 2019

Published in: Inverse Problems (Search for Journal in Brave)

Abstract: Quantitative photo-acoustic tomography (QPAT) seeks to reconstruct a distribution of optical attenuation coefficients inside a sample from a set of time series of pressure data that is measured outside the sample. The associated inverse problems involve two steps, namely acoustic and optical, which can be solved separately or as a direct composite problem. We adopt the latter approach for realistic acoustic media that possess heterogeneous and often not accurately known distributions for sound speed and ambient density, as well as an attenuation following a frequency power law that is evident in tissue media. We use a Diffusion Approximation (DA) model for the optical portion of the problem. We solve the corresponding composite inverse problem using three total variation (TV) regularised optimisation approaches. Accordingly, we develop two Krylov-subspace inexact-Newton algorithms that utilise the Jacobian matrix in a matrix-free manner in order to handle the computational cost. Additionally, we use a gradient-based algorithm that computes a search direction using the L-BFGS method, and applies a TV regularisation based on the Alternating Direction Method of Multipliers (ADMM) as a benchmark, because this method is popular for QPAT and direct QPAT. The results indicate the superiority of the developed inexact Newton algorithms over gradient-based Quasi-Newton approaches for a comparable computational complexity.

Full work available at URL: https://arxiv.org/abs/1810.11685

Recommendations

• A gradient-based method for quantitative photoacoustic tomography using the radiative transfer equation
• Quantitative photoacoustic imaging in the radiative transport regime
• A one-step reconstruction algorithm for quantitative photoacoustic imaging
• A hybrid reconstruction method for quantitative PAT
• Regularization approaches for quantitative photoacoustic tomography using the radiative transfer equation

zbMATH Keywords

total variationADMMphotoacoustic tomographyL-BFGSacoustic absorptionheterogeneous sound speed and density

Mathematics Subject Classification ID

Numerical optimization and variational techniques (65K10) Biomedical imaging and signal processing (92C55) Image processing (compression, reconstruction, etc.) in information and communication theory (94A08) Regularity of solutions in optimal control (49N60) Numerical computation of solutions to systems of equations (65H10) Numerical methods for inverse problems for boundary value problems involving PDEs (65N21) Biological applications of optics and electromagnetic theory (78A70) Newton-type methods (49M15) Numerical methods for inverse problems for initial value and initial-boundary value problems involving PDEs (65M32) Hydro- and aero-acoustics (76Q05)

Cites Work

• Explicit inversion formulae for the spherical mean Radon transform
• Inverse diffusion theory of photoacoustics
• Determining a Function from Its Mean Values Over a Family of Spheres
• Thermoacoustic tomography with variable sound speed
• A Nonlinear Primal-Dual Method for Total Variation-Based Image Restoration
• Quantitative photoacoustic tomography
• Inexact Newton Methods
• Iterative Methods for Total Variation Denoising
• Finite element model for the coupled radiative transfer equation and diffusion approximation
• Reconstruction and time reversal in thermoacoustic tomography in acoustically homogeneous and inhomogeneous media
• Multi-source quantitative photoacoustic tomography in a diffusive regime
• On multi-spectral quantitative photoacoustic tomography in diffusive regime
• A gradient-based method for quantitative photoacoustic tomography using the radiative transfer equation
• Photoacoustic tomography in absorbing acoustic media using time reversal
• Reconstructing absorption and scattering distributions in quantitative photoacoustic tomography
• Inverse transport theory of photoacoustics
• A one-step reconstruction algorithm for quantitative photoacoustic imaging
• A Bayesian approach to spectral quantitative photoacoustic tomography
• Single-stage reconstruction algorithm for quantitative photoacoustic tomography
• A direct method for photoacoustic tomography with inhomogeneous sound speed
• An efficient Neumann series-based algorithm for thermoacoustic and photoacoustic tomography with variable sound speed
• Time reversal in thermoacoustic tomography—an error estimate
• On the adjoint operator in photoacoustic tomography
• Analysis of Iterative Methods in Photoacoustic Tomography with Variable Sound Speed
• Photoacoustic and thermoacoustic tomography: image formation principles
• Joint reconstruction of the initial pressure and speed of sound distributions from combined photoacoustic and ultrasound tomography measurements
• Iterated preconditioned LSQR method for inverse problems on unstructured grids
• Limited-view multi-source quantitative photoacoustic tomography
• Efficient Inclusion of Total Variation Type Priors in Quantitative Photoacoustic Tomography
• A continuous adjoint for photo-acoustic tomography of the brain

Cited In (10)

• Determining the structure of a biological medium using acousto-optic probes
• A gradient-based method for quantitative photoacoustic tomography using the radiative transfer equation
• Photoacoustic tomography in absorbing acoustic media using time reversal
• A phase-field approach for detecting cavities via a Kohn-Vogelius type functional
• Regularization approaches for quantitative photoacoustic tomography using the radiative transfer equation
• Estimating absorption and scattering in quantitative photoacoustic tomography with an adaptive Monte Carlo method for light transport
• A linear state space model for photoacoustic imaging in an acoustic attenuating media
• Photoacoustic tomography with direction dependent data: an exact series reconstruction approach
• Quantitative photo-acoustic tomography with partial data
• Single-stage reconstruction algorithm for quantitative photoacoustic tomography

Uses Software

• k-Wave

This page was built for publication: Direct quantitative photoacoustic tomography for realistic acoustic media

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5228008)