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.

Preconditioned P-ULA for Joint Deconvolution-Segmentation of Ultrasound Images -- Extended Version

Joint Segmentation and Deconvolution of Ultrasound Images Using a Hierarchical Bayesian Model based on Generalized Gaussian Priors

Simultaneous Reconstruction and Segmentation of Dynamic Pet Via Low-Rank and Sparse Matrix Decomposition

CP-UNet: Contour-based Probabilistic Model for Medical Ultrasound Images Segmentation

A Unifying Approach to Inverse Problems of Ultrasound Beamforming and Deconvolution

Accelerated Bayesian Imaging by Relaxed Proximal-Point Langevin Sampling

Joint Blind Deconvolution and Robust Principal Component Analysis for Blood Flow Estimation in Medical Ultrasound Imaging

DeepUCT: Complex cascaded deep learning network for improved ultrasound tomography

Efficient Bayesian Computation by Proximal Markov Chain Monte Carlo: When Langevin Meets Moreau

A Wavelet Frame Method with Shape Prior for Ultrasound Video Segmentation.

Joint Calibrationless Reconstruction and Segmentation of Parallel MRI

Joint Segmentation and Image Reconstruction with Error Prediction in Photoacoustic Imaging using Deep Learning

Combining Bayesian and Deep Learning Methods for the Delineation of the Fan in Ultrasound Images

Vascular Segmentation of Functional Ultrasound Images using Deep Learning

Joint Motion Estimation and Segmentation from Undersampled Cardiac MR Image

A Partially Learned Algorithm for Joint Photoacoustic Reconstruction and Segmentation

PHOCUS: Physics-Based Deconvolution for Ultrasound Resolution Enhancement

Joint Segmentation and Landmark Localization of Fetal Femur in Ultrasound Volumes

Bayesian CNN for Segmentation Uncertainty Inference on 4D Ultrasound Images of the Femoral Cartilage for Guidance in Robotic Knee Arthroscopy

Enhancing super-resolution ultrasound localisation through multi-frame deconvolution exploiting spatiotemporal coherence