Non-Stationary Multi-layered Gaussian Priors for Bayesian Inversion

by   Muhammad Emzir, et al.

In this article, we study Bayesian inverse problems with multi-layered Gaussian priors. We first describe the conditionally Gaussian layers in terms of a system of stochastic partial differential equations. We build the computational inference method using a finite-dimensional Galerkin method. We show that the proposed approximation has a convergence-in-probability property to the solution of the original multi-layered model. We then carry out Bayesian inference using the preconditioned Crank–Nicolson algorithm which is modified to work with multi-layered Gaussian fields. We show via numerical experiments in signal deconvolution and computerized X-ray tomography problems that the proposed method can offer both smoothing and edge preservation at the same time.


Efficient algorithms for Bayesian Inverse Problems with Whittle–Matérn Priors

This paper tackles efficient methods for Bayesian inverse problems with ...

Bayesian linear inverse problems in regularity scales

We obtain rates of contraction of posterior distributions in inverse pro...

The SPDE Approach to Matérn Fields: Graph Representations

This paper investigates Gaussian Markov random field approximations to n...

Horseshoe priors for edge-preserving linear Bayesian inversion

In many large-scale inverse problems, such as computed tomography and im...

Bayesian inference for spatio-temporal spike-and-slab priors

In this work, we address the problem of solving a series of underdetermi...

Bayesian model selection for unsupervised image deconvolution with structured Gaussian priors

This paper considers the objective comparison of stochastic models to so...