Fast reconstruction approaches for photoacoustic tomography with smoothing Sobolev/Matérn priors

Abstract

In photoacoustic tomography (PAT), the computation of the initial pressure distribution within an object from its time-dependent boundary measurements over time is considered. This problem can be approached from two well-established points of view: deterministically using regularisation methods, or stochastically using the Bayesian framework. Both approaches frequently require the solution of a variational problem. In the paper we elaborate the connection between these approaches by establishing the equivalence between a smoothing Mat{é}rn class of covariance operators and Sobolev embedding operator $E_s:H^s\hookrightarrow L^2$. We further discuss the use of a Wavelet-based implementation of the adjoint operator $E_s^{\ast }$ which also allows for efficient evaluations for certain Mat{é}rn covariance operators, leading to efficient implementations both in terms of computational effort as well as memory requirements. The proposed methods are validated with reconstructions for the photoacoustic problem.

24 pages, 7 figures