We analyze the mathematical model of multiwave tomography with a variable speed with integrating measurements on planes tangent to a sphere surrounding the source. We prove sharp uniqueness and stability estimates with full and partial data and propose a time reversal algorithm which recovers the visible singularities.
We develop a paradigm using microlocal analysis that allows one to characterize the visible and added singularities in a broad range of incomplete data tomography problems. We give precise characterizations for photo- and thermoacoustic tomography and Sonar, and provide artifact reduction strategies. In particular, our theorems show that it is better to arrange Sonar detectors so that the boundary of the set of detectors does not have corners and is smooth. To illustrate our results, we provide reconstructions from synthetic spherical mean data as well as from experimental photoacoustic data.
Quantitative image reconstruction in photoacoustic tomography requires the solution of a coupled physics inverse problem involvier light transport and acoustic wave propagation. In this paper we address this issue employing the radiative transfer equation as accurate model for light transport. As main theoretical results, we derive several stability and uniqueness results for the linearized inverse problem. We consider the case of single illumination as well as the case of multiple illuminations assuming full or partial data. The numerical solution of the linearized problem is much less costly than the solution of the non-linear problem. We present numerical simulations supporting the stability results for the linearized problem and demonstrate that the linearized problem already gives accurate quantitative results.
Filtered backprojection (FBP) is an efficient and popular class of tomographic image reconstruction methods. In photoacoustic tomography, these algorithms are based on theoretically exact analytic inversion formulas which results in accurate reconstructions. However, photoacoustic measurement data are often incomplete (limited detection view and sparse sampling), which results in artefacts in the images reconstructed with FBP. In addition to that, properties such as directivity of the acoustic detectors are not accounted for in standard FBP, which affects the reconstruction quality, too. To account for these issues, in this papers we propose to improve FBP algorithms based on machine learning techniques. In the proposed method, we include additional weight factors in the FBP, that are optimized on a set of incomplete data and the corresponding ground truth photoacoustic source. Numerical tests show that the learned FBP improves the reconstruction quality compared to the standard FBP.
We investigate the long time behaviour of the $L^2$-energy of solutions to wave equations with variable speed. The novelty of the approach is the combination of estimates for higher order derivatives of the coefficient with a stabilisation property.
We present a pressure sensor based on a Michelson interferometer, for use in photoacoustic tomography. Quadrature phase detection is employed allowing measurement at any point on the mirror surface without having to retune the interferometer, as is typically required by Fabry-Perot type detectors. This opens the door to rapid full surface detection, which is necessary for clinical applications. Theory relating acoustic pressure to detected acoustic particle displacements is used to calculate the detector sensitivity, which is validated with measurement. Proof-of-concept tomographic images of blood vessel phantoms have been taken with sub-millimeter resolution at depths of several millimeters.