TY - JOUR
T1 - Theoretically exact photoacoustic reconstruction from spatially and temporally reduced data
AU - Do, N.
AU - Kunyansky, L.
N1 - Funding Information:
The authors would like to thank the anonymous referees for the numerous suggestions that helped to improve the manuscript. The second author is grateful for the partial support by the NSF through the awards NSF/DMS-1211521 and NSF/DMS-1418772.
Publisher Copyright:
© 2018 IOP Publishing Ltd.
PY - 2018/7/16
Y1 - 2018/7/16
N2 - We investigate the inverse source problem for the wave equation, arising in photo- and thermoacoustic tomography. There exist quite a few theoretically exact inversion formulas explicitly expressing the solution of this problem in terms of the measured data, under the assumption of the constant and known speed of sound. However, almost all of these formulas require data to be measured either on an unbounded surface, or on a closed surface completely surrounding the object. This is too restrictive for practical applications. The alternative approach we present, under certain restriction on geometry, yields a theoretically exact reconstruction of the standard Radon projections of the source from the data measured on a finite open surface. In addition, this technique reduces the time interval where the data should be known. In general, our method requires a pre-computation of densities of certain single-layer potentials. However, in the case of a truncated circular or spherical acquisition surface, these densities are easily obtained analytically, which leads to fully explicit asymptotically fast algorithms. We test these algorithms in a series of numerical simulations.
AB - We investigate the inverse source problem for the wave equation, arising in photo- and thermoacoustic tomography. There exist quite a few theoretically exact inversion formulas explicitly expressing the solution of this problem in terms of the measured data, under the assumption of the constant and known speed of sound. However, almost all of these formulas require data to be measured either on an unbounded surface, or on a closed surface completely surrounding the object. This is too restrictive for practical applications. The alternative approach we present, under certain restriction on geometry, yields a theoretically exact reconstruction of the standard Radon projections of the source from the data measured on a finite open surface. In addition, this technique reduces the time interval where the data should be known. In general, our method requires a pre-computation of densities of certain single-layer potentials. However, in the case of a truncated circular or spherical acquisition surface, these densities are easily obtained analytically, which leads to fully explicit asymptotically fast algorithms. We test these algorithms in a series of numerical simulations.
KW - explicit inversion formula
KW - photoacoustic tomography
KW - reduced data
KW - spherical means
KW - thermoacoustic tomography
KW - wave equation
UR - http://www.scopus.com/inward/record.url?scp=85051200369&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85051200369&partnerID=8YFLogxK
U2 - 10.1088/1361-6420/aacfac
DO - 10.1088/1361-6420/aacfac
M3 - Article
AN - SCOPUS:85051200369
VL - 34
JO - Inverse Problems
JF - Inverse Problems
SN - 0266-5611
IS - 9
M1 - 094004
ER -