Journal
SIAM JOURNAL ON IMAGING SCIENCES
Volume 10, Issue 2, Pages 751-781Publisher
SIAM PUBLICATIONS
DOI: 10.1137/16M1104822
Keywords
photoacoustic tomography; variable sound speed; iterative regularization; adjoint operator; Landweber's method; Nesterov's method; CG method; visibility condition; invisibility condition; image reconstruction
Categories
Funding
- NSF [DMS 1212125, DMS 1616904]
- Direct For Mathematical & Physical Scien
- Division Of Mathematical Sciences [1616904] Funding Source: National Science Foundation
Ask authors/readers for more resources
In this article, we revisit iterative methods for solving the inverse source 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 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.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available