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

Fuente: arXiv
Gespeichert in:
Bibliographische Detailangaben
Hauptverfasser: Kultima, Jaakko, Ramlau, Ronny, Sahlström, Teemu, Tarvainen, Tanja
Format: Preprint
Veröffentlicht: 2025
Schlagworte:
Online-Zugang:
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
_version_ 1866918081212710912
author Kultima, Jaakko
Ramlau, Ronny
Sahlström, Teemu
Tarvainen, Tanja
author_facet Kultima, Jaakko
Ramlau, Ronny
Sahlström, Teemu
Tarvainen, Tanja
contents 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^*$ 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.
format Preprint
id arxiv_https___arxiv_org_abs_2507_02401
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Fast reconstruction approaches for photoacoustic tomography with smoothing Sobolev/Matérn priors
Kultima, Jaakko
Ramlau, Ronny
Sahlström, Teemu
Tarvainen, Tanja
Numerical Analysis
Mathematical Physics
35R30, 65M12, 92C55, 62P10
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^*$ 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.
title Fast reconstruction approaches for photoacoustic tomography with smoothing Sobolev/Matérn priors
topic Numerical Analysis
Mathematical Physics
35R30, 65M12, 92C55, 62P10
url https://arxiv.org/abs/2507.02401