A TV-Gaussian prior for infinite-dimensional Bayesian inverse problems and its numerical implementations

Abstract

Many scientific and engineering problems require to perform Bayesian inferences in function spaces, in which the unknowns are of infinite dimension. In such problems, choosing an appropriate prior distribution is an important task. In particular we consider problems where the function to infer is subject to sharp jumps which render the commonly used Gaussian measures unsuitable. On the other hand, the so-called total variation (TV) prior can only be defined in a finite dimensional setting, and does not lead to a well-defined posterior measure in function spaces. In this work we present a TV-Gaussian (TG) prior to address such problems, where the TV term is used to detect sharp jumps of the function, and the Gaussian distribution is used as a reference measure so that it results in a well-defined posterior measure in the function space. We also present an efficient Markov Chain Monte Carlo (MCMC) algorithm to draw samples from the posterior distribution of the TG prior. With numerical examples we demonstrate the performance of the TG prior and the efficiency of the proposed MCMC algorithm.

References in corpus (5)

• Dimension-independent likelihood-informed MCMC
• Optimal scalings for local Metropolis--Hastings chains on nonproduct targets in high dimensions
• Discretization-invariant Bayesian inversion and Besov space priors
• An adaptive independence sampler MCMC algorithm for infinite dimensional Bayesian inferences
• A note on the Karhunen-Loève expansions for infinite-dimensional Bayesian inverse problems

Cited by in corpus (10)

• Prior normalization for certified likelihood-informed subspace detection of Bayesian inverse problems
• On an adaptive preconditioned Crank-Nicolson MCMC algorithm for infinite dimensional Bayesian inferences
• Spatiotemporal Besov Priors for Bayesian Inverse Problems
• A Hadamard fractioal total variation-Gaussian (HFTG) prior for Bayesian inverse problems
• NLTG Priors in Medical Image: Nonlocal TV-Gaussian (NLTG) prior for Bayesian inverse problems with applications to Limited CT Reconstruction
• Optimization-Based MCMC Methods for Nonlinear Hierarchical Statistical Inverse Problems
• A general fractional total variation-Gaussian (GFTG) prior for Bayesian inverse problems
• Bayesian inference and uncertainty quantification for image reconstruction with Poisson data
• Bayesian identification of discontinuous fields with an ensemble-based variable separation multiscale method
• A persistent-homology-based Bayesian prior for potential coefficient reconstruction in an elliptic PDE