Analysis of Iterative Methods in Photoacoustic Tomography with Variable Sound Speed

DOI10.1137/16M1104822MaRDI QIDQ5266392zbMATH OpenOpenAlexFDO

Authors

Publication date 2 June 2017

Published in SIAM Journal on Imaging Sciences (Search for Journal in Brave)

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

zbMATH Keywords

adjoint operatorCG methodimage reconstructioninvisibility conditioniterative regularizationLandweber's methodNesterov's methodphotoacoustic tomographyvariable sound speedvisibility condition

Mathematics Subject Classification ID

Wave front sets in context of PDEs (35A18) Inverse problems for PDEs (35R30) Iterative numerical methods for linear systems (65F10) Linear waves in solid mechanics (74J05) Biomedical imaging and signal processing (92C55)

Abstract: In this article, we revisit iterative methods for solving the inverse problem of photoacoustic tomography in free space. Recently, there have been interesting developments on explicit formulations of the adjoint operator, demonstrating that iterative methods are an attractive choice for photoacoustic image reconstruction. In this work, we propose several modifications of current formulations of the adjoint operator which help speed up the convergence and yield improved error estimates. We establish a stability analysis and show that, with our choices of the adjoint operator, the iterative methods can achieve a linear rate of convergence, in the $L^2$-norm as well as in the $H^1$-norm. In addition, we analyze the normal operator from the microlocal analysis point of view. This gives insight into the convergence speed of the iterative methods and choosing proper weights for the mapping spaces. Finally, we present numerical results using various iterative reconstruction methods for full as well as limited view data. Our results demonstrate that Nesterov's fast gradient and the CG methods converge faster than Landweber's and iterative time reversal methods in the visible as well as the invisible case.

Recommendations

• Variational iterative algorithms in photoacoustic tomography with variable sound speed
• Analysis for full-field photoacoustic tomography with variable sound speed
• Iterative methods for photoacoustic tomography in attenuating acoustic media
• A direct method for photoacoustic tomography with inhomogeneous sound speed
• On iterative algorithms for quantitative photoacoustic tomography in the radiative transport regime
• Parameterized joint reconstruction of the initial pressure and sound speed distributions for photoacoustic computed tomography
• A variational method for quantitative photoacoustic tomography with piecewise constant coefficients
• An efficient Neumann series-based algorithm for thermoacoustic and photoacoustic tomography with variable sound speed
• Uniqueness of reconstruction and an inversion procedure for thermoacoustic and photoacoustic tomography with variable sound speed
• Photoacoustic and thermoacoustic tomography with an uncertain wave speed

Cites work

• A direct method for photoacoustic tomography with inhomogeneous sound speed
• A dissipative time reversal technique for photoacoustic tomography in a cavity
• A family of inversion formulas in thermoacoustic tomography
• A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems
• A series solution and a fast algorithm for the inversion of the spherical mean Radon transform
• A uniform reconstruction formula in integral geometry
• An efficient Neumann series-based algorithm for thermoacoustic and photoacoustic tomography with variable sound speed
• Applied functional analysis. Functional analysis, Sobolev spaces and elliptic differential equations
• Artifacts in Incomplete Data Tomography with Applications to Photoacoustic Tomography and Sonar
• Circular integrating detectors in photo and thermoacoustic tomography
• Determining a Function from Its Mean Values Over a Family of Spheres
• Explicit inversion formulae for the spherical mean Radon transform
• scientific article; zbMATH DE number 3709065 (Why is no real title available?)
• scientific article; zbMATH DE number 3709066 (Why is no real title available?)
• scientific article; zbMATH DE number 556557 (Why is no real title available?)
• scientific article; zbMATH DE number 1181255 (Why is no real title available?)
• scientific article; zbMATH DE number 1552299 (Why is no real title available?)
• scientific article; zbMATH DE number 2164725 (Why is no real title available?)
• scientific article; zbMATH DE number 783550 (Why is no real title available?)
• scientific article; zbMATH DE number 936298 (Why is no real title available?)
• scientific article; zbMATH DE number 3353865 (Why is no real title available?)
• Introductory lectures on convex optimization. A basic course.
• Inversion of circular means and the wave equation on convex planar domains
• Inversion of Spherical Means and the Wave Equation in Even Dimensions
• Mathematics of thermoacoustic tomography
• Methods of conjugate gradients for solving linear systems
• Multiwave imaging in an enclosure with variable wave speed
• Multiwave tomography with reflectors: Landweber's iteration
• Non homogeneous boundary value problems for second order hyperbolic operators
• On artifacts in limited data spherical Radon transform: curved observation surface
• On artifacts in limited data spherical Radon transform: flat observation surfaces
• On singularities and instability of reconstruction in thermoacoustic tomography
• On the adjoint operator in photoacoustic tomography
• On the Convergence of the Conjugate Gradient Method for Singular Linear Operator Equations
• ON THE RATE OF CONVERGENCE OF THE CONJUGATE GRADIENT METHOD FOR LINEAR OPERATORS IN HILBERT SPACE
• Photo-acoustic inversion in convex domains
• Photoacoustic and thermoacoustic tomography: image formation principles
• Preconditioning discretizations of systems of partial differential equations.
• Reconstruction and time reversal in thermoacoustic tomography in acoustically homogeneous and inhomogeneous media
• Some Convergence Properties of the Conjugate Gradient Method in Hilbert Space
• Superlinear convergence of Krylov subspace methods for self-adjoint problems in Hilbert space
• The Conjugate Gradient Method for Linear and Nonlinear Operator Equations
• The Quasi-Reversibility Method for Thermoacoustic Tomography in a Heterogeneous Medium
• The Radon transform and medical imaging
• Thermoacoustic computed tomography with large planar receivers
• Thermoacoustic tomography arising in brain imaging
• Thermoacoustic tomography with variable sound speed
• Time reversal in thermoacoustic tomography—an error estimate
• Uniqueness of reconstruction and an inversion procedure for thermoacoustic and photoacoustic tomography with variable sound speed
• Universal inversion formulas for recovering a function from spherical means

Cited in

(27)

• On the convergence of the time reversal method for thermoacoustic tomography in elastic media
• A convergent numerical method to recover the initial condition of nonlinear parabolic equations from lateral Cauchy data
• Operator learning approach for the limited view problem in photoacoustic tomography
• Reconstruction algorithms for photoacoustic tomography in heterogeneous damping media
• A direct method for photoacoustic tomography with inhomogeneous sound speed
• A continuous adjoint for photo-acoustic tomography of the brain
• Iterative methods for photoacoustic tomography in attenuating acoustic media
• A Galerkin least squares approach for photoacoustic tomography
• Efficient regularization with wavelet sparsity constraints in photoacoustic tomography
• Deep learning for photoacoustic tomography from sparse data
• Recovering the initial condition of parabolic equations from lateral Cauchy data via the quasi-reversibility method
• Simultaneous reconstruction of optical absorption property and speed of sound in intravascular photoacoustic tomography
• Analysis for full-field photoacoustic tomography with variable sound speed
• Explicit inversion formulas for the two-dimensional wave equation from Neumann traces
• Well-posedness for photoacoustic tomography with Fabry-Perot sensors
• Regularized iterative weighted filtered back-projection for few-view data photoacoustic imaging
• Direct quantitative photoacoustic tomography for realistic acoustic media
• Photoacoustic tomography with direction dependent data: an exact series reconstruction approach
• Optoacoustic image reconstruction: the full inverse problem with variable bases
• Variational iterative algorithms in photoacoustic tomography with variable sound speed
• Parametrix for the inverse source problem of thermoacoustic tomography with reduced data
• Multiscale factorization of the wave equation with application to compressed sensing photoacoustic tomography
• Data-Driven Morozov Regularization of Inverse Problems
• A PINN-driven game-theoretic framework in limited data photoacoustic tomography
• On the optimal choice of the illumination function in photoacoustic tomography
• Recover damping and initial condition from a single measurement for a wave equation
• Determination of the sound speed and an initial source in photoacoustic tomography

Describes a project that uses

Uses Software

• k-Wave

This page was built for publication: Analysis of Iterative Methods in Photoacoustic Tomography with Variable Sound Speed

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