In this article we study the problem of recovering the initial data of the twodimensional wave equation from Neumann measurements on a convex domain Ω ⊂ R 2 with smooth boundary. We derive an explicit inversion formula of a so-called backprojection type and deduce exact inversion formulas for circular and elliptical domains. In addition, for circular domains, we show that the initial data can also be recovered from any linear combination of its solution and its normal derivative on the boundary. Numerical results of our implementation of the derived inversion formulas are presented demonstrating their accuracy and stability.

In photoacoustic tomography, one is interested to recover the initial pressure distribution inside a tissue from the corresponding measurements of the induced acoustic wave on the boundary of a region enclosing the tissue. In the limited view problem, the wave boundary measurements are given on the part of the boundary, whereas in the full view problem, the measurements are known on the whole boundary. For the full view problem, there exist various fast and robust reconstruction methods. These methods give severe reconstruction artifacts when they are applied directly to the limited view data. One approach for reducing such artefacts is trying to extend the limited view data to the whole region boundary, and then use existing reconstruction methods for the full view data. In this paper, we propose an operator learning approach for constructing an operator that gives an approximate extension of the limited view data. We consider the behavior of a reconstruction formula on the extended limited view data that is given by our proposed approach. Approximation errors of our approach are analyzed. We also present numerical results with the proposed extension approach supporting our theoretical analysis.

No abstract

We investigate resolution in photoacoustic tomography (PAT). Using Shannon theory, we investigate the theoretical resolution limit of sparse view PAT theoretically, and empirically demonstrate that all reconstruction methods used exceed this limit.

We study the inverse source problem in photoacoustic tomography (PAT) for mixed data, which denote a weighted linear combination of the acoustic pressure and its normal derivative on an observation surface. We consider in particular the case where the data are only available on finite time intervals, which accounts for real-world usage of PAT where data are only feasible within a certain time interval. Extending our previous work, we derive explicit formulas up to a smoothing integral on convex domains with a smooth boundary, yielding exact reconstruction for circular or elliptical domains. We also present numerical reconstructions of our new exact inversion formulas on finite time intervals and compare them with the reconstructions of our previous formulas for unlimited time wave measurements.