DeepAI AI Chat
Log In Sign Up

Bayesian inversion for electromyography using low-rank tensor formats

by   Anna Rörich, et al.

The reconstruction of the structure of biological tissue using electromyographic data is a non-invasive imaging method with diverse medical applications. Mathematically, this process is an inverse problem. Furthermore, electromyographic data are highly sensitive to changes in the electrical conductivity that describes the structure of the tissue. Modeling the inevitable measurement error as a stochastic quantity leads to a Bayesian approach. Solving the discretized Bayes-inverse problem means drawing samples from the posterior distribution of parameters, e.g., the conductivity, given measurement data. Using, e.g., a Metropolis-Hastings algorithm for this purpose involves solving the forward problem for different parameter combinations which requires a high computational effort. Low-rank tensor formats can reduce this effort by providing a data-sparse representation of all occurring linear systems of equations simultaneously and allow for their efficient solution. The application of Bayes' theorem proves the well-posedness of the Bayes-inverse problem. The derivation and proof of a low-rank representation of the forward problem allow for the precomputation of all solutions of this problem under certain assumptions, resulting in an efficient and theory-based sampling algorithm. Numerical experiments support the theoretical results, but also indicate that a high number of samples is needed to obtain reliable estimates for the parameters. The Metropolis-Hastings sampling algorithm, using the precomputed forward solution in a tensor format, draws this high number of samples and therefore enables solving problems which are infeasible using classical methods.


page 1

page 2

page 3

page 4


Alternating linear scheme in a Bayesian framework for low-rank tensor approximation

Multiway data often naturally occurs in a tensorial format which can be ...

A parameter-dependent smoother for the multigrid method

The solution of parameter-dependent linear systems, by classical methods...

A model reduction approach for inverse problems with operator valued data

We study the efficient numerical solution of linear inverse problems wit...

Low-rank tensor methods for Markov chains with applications to tumor progression models

Continuous-time Markov chains describing interacting processes exhibit a...

Active Sampling for Accelerated MRI with Low-Rank Tensors

Magnetic resonance imaging (MRI) is a powerful imaging modality that rev...

A Low-rank Solver for the Stochastic Unsteady Navier-Stokes Problem

We study a low-rank iterative solver for the unsteady Navier-Stokes equa...

A low-rank tensor method to reconstruct sparse initial states for PDEs with Isogeometric Analysis

When working with PDEs the reconstruction of a previous state often prov...