A partially reflecting random walk on spheres algorithm for electrical impedance tomography

DOI10.1016/J.JCP.2015.10.005zbMATH Open1349.94030arXiv1502.04318OpenAlexW1828962914MaRDI QIDQ2374877FDOQ2374877

Authors: Sylvain Maire, Martin Simon 

Publication date: 5 December 2016

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: In this work, we develop a probabilistic estimator for the voltage-to-current map arising in electrical impedance tomography. This novel so-called partially reflecting random walk on spheres estimator enables Monte Carlo methods to compute the voltage-to-current map in an embarrassingly parallel manner, which is an important issue with regard to the corresponding inverse problem. Our method uses the well-known random walk on spheres algorithm inside subdomains where the diffusion coefficient is constant and employs replacement techniques motivated by finite difference discretization to deal with both mixed boundary conditions and interface transmission conditions. We analyze the global bias and the variance of the new estimator both theoretically and experimentally. In a second step, the variance is considerably reduced via a novel control variate conditional sampling technique.

Full work available at URL: https://arxiv.org/abs/1502.04318

Recommendations

• Reconstruction algorithm based on stochastic Galerkin finite element method for electrical impedance tomography
• Statistical inversion and Monte Carlo sampling methods in electrical impedance tomography
• Randomize-then-optimize for sampling and uncertainty quantification in electrical impedance tomography
• Stochastic Galerkin Finite Element Method with Local Conductivity Basis for Electrical Impedance Tomography
• Adaptive and stochastic algorithms for electrical impedance tomography and DC resistivity problems with piecewise constant solutions and many measurements

zbMATH Keywords

variance reductionMonte Carlo methodsrandom walk on sphereselectrical impedance tomographydiscontinuous diffusion coefficientreflecting Brownian motionrandom diffusion coefficient

Mathematics Subject Classification ID

Monte Carlo methods (65C05) Biomedical imaging and signal processing (92C55) Image processing (compression, reconstruction, etc.) in information and communication theory (94A08)

Cites Work

• Statistical and computational inverse problems.
• Uniqueness in Calderón's problem with Lipschitz conductivities
• Calderón's inverse conductivity problem in the plane
• Two-stage reconstruction of a circular anomaly in electrical impedance tomography
• Title not available (Why is that?)
• Title not available (Why is that?)
• Existence and Uniqueness for Electrode Models for Electric Current Computed Tomography
• Title not available (Why is that?)
• Electrical Impedance Tomography
• Stable determination of conductivity by boundary measurements
• Simulating diffusion processes in discontinuous media: a numerical scheme with constant time steps
• Title not available (Why is that?)
• New Monte Carlo schemes for simulating diffusions in discontinuous media
• Monte Carlo approximations of the Neumann problem
• Some Continuous Monte Carlo Methods for the Dirichlet Problem
• Statistical inversion and Monte Carlo sampling methods in electrical impedance tomography
• Simultaneous reconstruction of outer boundary shape and admittivity distribution in electrical impedance tomography
• Monte Carlo Methods for Calculating Some Physical Properties of Large Molecules
• Simulation of a stochastic process in a discontinuous layered medium
• Numerical Approximation for Functionals of Reflecting Diffusion Processes
• Anomaly Detection in Random Heterogeneous Media
• Probabilistic interpretation and random walk on spheres algorithms for the Poisson-Boltzmann equation in molecular dynamics
• Simulating diffusions with piecewise constant coefficients using a kinetic approximation
• Euler schemes and half-space approximation for the simulation of diffusion in a domain
• \(\varepsilon\)-shell error analysis for ``walk on spheres algorithms
• Efficient schemes for the weak approximation of reflected diffusions
• Sequential Control Variates for Functionals of Markov Processes
• Integral Formulation of the Boundary Value Problems and the Method of Random Walk on Spheres
• Bayesian anomaly detection in heterogeneous media with applications to geophysical tomography

Cited In (3)

• A parallel iterative probabilistic method for mixed problems of Laplace equations with the Feynman-Kac formula of killed Brownian motions
• A path integral Monte Carlo (PIMC) method based on Feynman-Kac formula for electrical impedance tomography
• An implementation of Milstein's method for general bounded diffusions

Uses Software

• EIDORS

This page was built for publication: A partially reflecting random walk on spheres algorithm for electrical impedance tomography

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2374877)