Maximum-a-posteriori estimation with Bayesian confidence regions
From MaRDI portal
Publication:5266375
Abstract: Solutions to inverse problems that are ill-conditioned or ill-posed may have significant intrinsic uncertainty. Unfortunately, analysing and quantifying this uncertainty is very challenging, particularly in high-dimensional problems. As a result, while most modern mathematical imaging methods produce impressive point estimation results, they are generally unable to quantify the uncertainty in the solutions delivered. This paper presents a new general methodology for approximating Bayesian high-posterior-density credibility regions in inverse problems that are convex and potentially very high-dimensional. The approximations are derived by using recent concentration of measure results related to information theory for log-concave random vectors. A remarkable property of the approximations is that they can be computed very efficiently, even in large-scale problems, by using standard convex optimisation techniques. In particular, they are available as a by-product in problems solved by maximum-a-posteriori estimation. The approximations also have favourable theoretical properties, namely they outer-bound the true high-posterior-density credibility regions, and they are stable with respect to model dimension. The proposed methodology is illustrated on two high-dimensional imaging inverse problems related to tomographic reconstruction and sparse deconvolution, where the approximations are used to perform Bayesian hypothesis tests and explore the uncertainty about the solutions, and where proximal Markov chain Monte Carlo algorithms are used as benchmark to compute exact credible regions and measure the approximation error.
Recommendations
- Scalable Bayesian uncertainty quantification in imaging inverse problems via convex optimization
- Maximum a posteriori estimates in linear inverse problems with log-concave priors are proper Bayes estimators
- Fast algorithms for Bayesian uncertainty quantification in large-scale linear inverse problems based on low-rank partial Hessian approximations
- MCMC-based image reconstruction with uncertainty quantification
- A proximal Markov chain Monte Carlo method for Bayesian inference in imaging inverse problems: when Langevin meets Moreau
Cites work
- scientific article; zbMATH DE number 1350773 (Why is no real title available?)
- scientific article; zbMATH DE number 2061729 (Why is no real title available?)
- scientific article; zbMATH DE number 2107836 (Why is no real title available?)
- A new semiblind deconvolution approach for Fourier-based image restoration: an application in astronomy
- A nonlocal Bayesian image denoising algorithm
- Alternating direction method for image inpainting in wavelet domains
- An Augmented Lagrangian Approach to the Constrained Optimization Formulation of Imaging Inverse Problems
- Bayesian computation: a summary of the current state, and samples backwards and forwards
- Compressed sensing
- Compressive Blind Image Deconvolution
- Concentration of the information in data with log-concave distributions
- Dictionary learning for noisy and incomplete hyperspectral images
- Empirical Bayes Confidence Intervals Based on Bootstrap Samples
- Large‐Scale Inverse Problems and Quantification of Uncertainty
- Learning Multiscale Sparse Representations for Image and Video Restoration
- Phase retrieval via matrix completion
- Photographing paintings by image fusion
- Proximal Markov chain Monte Carlo algorithms
- Proximal splitting methods in signal processing
- Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information
- Simultaneous statistical inference. 2nd ed
- Sparse and Redundant Modeling of Image Content Using an Image-Signature-Dictionary
- Statistical and computational inverse problems.
- Super-resolution of positive sources: the discrete setup
- The Bernstein-von Mises theorem and nonregular models
- The Entropy Per Coordinate of a Random Vector is Highly Constrained Under Convexity Conditions
- The Variational Gaussian Approximation Revisited
- Variational Bayesian Super Resolution
Cited in
(16)- Patch-based image restoration using expectation propagation
- Efficient Bayesian Computation for Low-Photon Imaging Problems
- A proximal Markov chain Monte Carlo method for Bayesian inference in imaging inverse problems: when Langevin meets Moreau
- Sparse Online Variational Bayesian Regression
- On maximum a posteriori estimation with Plug \& Play priors and stochastic gradient descent
- Bayesian inference and uncertainty quantification for medical image reconstruction with Poisson data
- Bayesian mode and maximum estimation and accelerated rates of contraction
- Bayesian imaging inverse problem with SA-roundtrip prior via HMC-pCN sampler
- Posterior-Variance-Based Error Quantification for Inverse Problems in Imaging
- Accelerating Proximal Markov Chain Monte Carlo by Using an Explicit Stabilized Method
- Normalizing constants of log-concave densities
- Efficient Bayesian computation by proximal Markov chain Monte Carlo: when Langevin meets Moreau
- Equivalence of weak and strong modes of measures on topological vector spaces
- The Split Gibbs Sampler Revisited: Improvements to Its Algorithmic Structure and Augmented Target Distribution
- Scalable Bayesian uncertainty quantification in imaging inverse problems via convex optimization
- Accelerated Bayesian imaging by relaxed proximal-point Langevin sampling
This page was built for publication: Maximum-a-posteriori estimation with Bayesian confidence regions
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5266375)
