A Bayesian level-set inversion method for simultaneous reconstruction of absorption and diffusion coefficients in diffuse optical tomography

Anuj Abhishek,Thilo Strauss,Taufiquar Khan
2024-04-18
Abstract:In this article, we propose a non-parametric Bayesian level-set method for simultaneous reconstruction of piecewise constant diffusion and absorption coefficients in diffuse optical tomography (DOT). We show that the Bayesian formulation of the corresponding inverse problem is well-posed and the posterior measure as a solution to the inverse problem satisfies a Lipschitz estimate with respect to the measured data in terms of Hellinger distance. We reduce the problem to a shape-reconstruction problem and use level-set priors for the parameters of interest. We demonstrate the efficacy of the proposed method using numerical simulations by performing reconstructions of the original phantom using two reconstruction methods. Posing the inverse problem in a Bayesian paradigm allows us to do statistical inference for the parameters of interest whereby we are able to quantify the uncertainty in the reconstructions for both the methods. This illustrates a key advantage of Bayesian methods over traditional algorithms.
Applications,Numerical Analysis,Statistics Theory
What problem does this paper attempt to address?
This paper proposes a reconstruction method for Diffuse Optical Tomography (DOT) based on the Bayesian Hierarchical Model, which aims to simultaneously recover the absorption and diffusion coefficients. The paper demonstrates that the inverse problem is well-posed within the Bayesian framework, and quantifies the uncertainty in the reconstruction by using hierarchical priors and parameter sampling. The study highlights the advantages of this method over traditional algorithms in quantifying uncertainty. Additionally, it presents the current research status in the related field and presents numerical simulation results to support the theoretical findings.