Abstract:SIAM Journal on Imaging Sciences, Volume 17, Issue 2, Page 1078-1117, June 2024. This paper presents a new accelerated proximal Markov chain Monte Carlo methodology to perform Bayesian inference in imaging inverse problems with an underlying convex geometry. The proposed strategy takes the form of a stochastic relaxed proximal-point iteration that admits two complementary interpretations. For models that are smooth or regularized by Moreau–Yosida smoothing, the algorithm is equivalent to an implicit midpoint discretization of an overdamped Langevin diffusion targeting the posterior distribution of interest. This discretization is asymptotically unbiased for Gaussian targets and shown to converge in an accelerated manner for any target that is [math]-strongly log-concave (i.e., requiring in the order of [math] iterations to converge, similar to accelerated optimization schemes), comparing favorably to Pereyra, Vargas Mieles, and Zygalakis [SIAM J. Imaging Sci., 13 (2020), pp. 905–935], which is only provably accelerated for Gaussian targets and has bias. For models that are not smooth, the algorithm is equivalent to a Leimkuhler–Matthews discretization of a Langevin diffusion targeting a Moreau–Yosida approximation of the posterior distribution of interest and hence achieves a significantly lower bias than conventional unadjusted Langevin strategies based on the Euler–Maruyama discretization. For targets that are [math]-strongly log-concave, the provided nonasymptotic convergence analysis also identifies the optimal time step, which maximizes the convergence speed. The proposed methodology is demonstrated through a range of experiments related to image deconvolution with Gaussian and Poisson noise with assumption-driven and data-driven convex priors. Source codes for the numerical experiments of this paper are available from https://github.com/MI2G/accelerated-langevin-imla.

Edge-promoting adaptive Bayesian experimental design for X-ray imaging

Source Design Optimization for Depth Image Reconstruction in X-ray Imaging

Enhancing Industrial X-ray Tomography by Data-Centric Statistical Methods

Introduction To Gaussian Process Regression In Bayesian Inverse Problems, With New ResultsOn Experimental Design For Weighted Error Measures

Alternating Minimization Algorithm with Automatic Relevance Determination for Transmission Tomography under Poisson Noise

Accelerated Bayesian Imaging by Relaxed Proximal-Point Langevin Sampling

A Bayesian approach for edge detection in medical ultrasound images

Horseshoe priors for edge-preserving linear Bayesian inversion

Non-intrusive optimal experimental design for large-scale nonlinear Bayesian inverse problems using a Bayesian approximation error approach

Cone Beam X-Ray Luminescence Computed Tomography Based on Bayesian Method

A Bayesian Approach for Edge Extraction in Ultrasound Images and Its Application to Image Segmentation

Bayesian Reconstructions with Pde Image Model for Emission Tomography

Optimizing Quantitative Photoacoustic Imaging Systems: The Bayesian Cramér-Rao Bound Approach

Tractable Optimal Experimental Design using Transport Maps

High-Fidelity Modeling of Detector Lag and Gantry Motion in CT Reconstruction

Nonlocal TV-Gaussian Prior for Bayesian Inverse Problems with Applications to Limited CT Reconstruction

Dynamic SPECT reconstruction with temporal edge correlation

An Alternative Bayesian Reconstruction of Sparse-View CT by Optimizing Deep Learning Parameters

Blind beam-hardening correction from Poisson measurements

Generalized Adaptive Gaussian Markov Random Field for X-Ray Luminescence Computed Tomography