Iterative Methods for Photoacoustic Tomography in Attenuating Acoustic Media
arXiv:1704.07422 · doi:10.1088/1361-6420/aa8cba
Abstract
The development of efficient and accurate reconstruction methods is an important aspect of tomographic imaging. In this article, we address this issue for photoacoustic tomography. To this aim, we use models for acoustic wave propagation accounting for frequency dependent attenuation according to a wide class of attenuation laws that may include memory. We formulate the inverse problem of photoacoustic tomography in attenuating medium as an ill-posed operator equation in a Hilbert space framework that is tackled by iterative regularization methods. Our approach comes with a clear convergence analysis. For that purpose we derive explicit expressions for the adjoint problem that can efficiently be implemented. In contrast to time reversal, the employed adjoint wave equation is again damping and, thus has a stable solution. This stability property can be clearly seen in our numerical results. Moreover, the presented numerical results clearly demonstrate the Efficiency and accuracy of the derived iterative reconstruction algorithms in various situations including the limited view case.
24 pages, 7 figures
References in corpus (6)
- Regularization of systems of nonlinear ill-posed equations: I. Convergence Analysis
- On the Adjoint Operator in Photoacoustic Tomography
- A series solution and a fast algorithm for the inversion of the spherical mean Radon transform
- Regularization of systems of nonlinear ill-posed equations: II. Applications
- Reconstruction for multiwave imaging in attenuating media with large damping coefficient
- Thermoacoustic tomography for an integro-differential wave equation modeling attenuation
Cited by in corpus (4)
- A continuous adjoint for photo-acoustic tomography of the brain
- Recovery of Pressure and Wave Speed for Photoacoustic Imaging under a Condition of Relative Uncertainty
- Reconstruction formulas for Photoacoustic Imaging in Attenuating Media
- Well-posedness for Photoacoustic Tomography with Fabry-Perot Sensors