Open Access Open Access  Restricted Access Subscription or Fee Access

A fully Bayesian approach to shape estimation of objects from tomography data using MFS forward solutions

Robert G. Aykroyd, Daniel Lesnic, Andreas Karageorghis

Abstract


It is possible to characterise the aim of many practical inverse geometric problems as one of identifying the shape of an object within some domain of interest using non-intrusive measurements collected on the boundary of the domain. In the problem considered here the object is a rigid inclusion within a homogeneous background medium of constant conductivity, and the data are potential and current flux measurements made on the boundary of the region. The rigid inclusion is described using a geometric parametrization in terms of a star-shaped object. A Bayesian modelling approach is used to combine data likelihood and prior information, and posterior estimation is based on a Markov chain Monte Carlo algorithm which provides measures of uncertainty, as well as point estimates. This means that the inverse problem is never solved directly, but the cost is that instead the forward solution must be found many thousands of times. The forward problem is solved using the method of fundamental solutions (MFS) which is an efficient meshless alternative to the more common finite element or boundary element methods. This paper is the first to apply Bayesian modelling to a problem using the MFS, with numerical results demonstrating that for appropriate choices of prior distributions accurate results are possible. Further, it demonstrates that a fully Bayesian approach is possible where all prior smoothing parameters are estimated. It is important to note that the geometric modelling and statistical estimation approach are not limited to the inverse rigid inclusion/cavity under study and hence, the general technique can be easily applied to other inverse problems. A great benefit
of the approach is that it allows an intuitive model description and directly interpretable output. The methods are illustrated using numerical simulations.

Keywords


Bayesian inversion, Inverse conductivity problem, Markov chain Monte Carlo, Method of fundamental solutions, Posterior estimation, Statistical modelling.

Full Text:

PDF


Disclaimer/Regarding indexing issue:

We have provided the online access of all issues and papers to the indexing agencies (as given on journal web site). It’s depend on indexing agencies when, how and what manner they can index or not. Hence, we like to inform that on the basis of earlier indexing, we can’t predict the today or future indexing policy of third party (i.e. indexing agencies) as they have right to discontinue any journal at any time without prior information to the journal. So, please neither sends any question nor expects any answer from us on the behalf of third party i.e. indexing agencies.Hence, we will not issue any certificate or letter for indexing issue. Our role is just to provide the online access to them. So we do properly this and one can visit indexing agencies website to get the authentic information. Also: DOI is paid service which provided by a third party. We never mentioned that we go for this for our any journal. However, journal have no objection if author go directly for this paid DOI service.