Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods

Abstract

Uncertainty quantification is a critical missing component in radio interferometric imaging that will only become increasingly important as the big-data era of radio interferometry emerges. Since radio interferometric imaging requires solving a high-dimensional, ill-posed inverse problem, uncertainty quantification is difficult but also critical to the accurate scientific interpretation of radio observations. Statistical sampling approaches to perform Bayesian inference, like Markov Chain Monte Carlo (MCMC) sampling, can in principle recover the full posterior distribution of the image, from which uncertainties can then be quantified. However, traditional high-dimensional sampling methods are generally limited to smooth (e.g. Gaussian) priors and cannot be used with sparsity-promoting priors. Sparse priors, motivated by the theory of compressive sensing, have been shown to be highly effective for radio interferometric imaging. In this article proximal MCMC methods are developed for radio interferometric imaging, leveraging proximal calculus to support non-differential priors, such as sparse priors, in a Bayesian framework. Furthermore, three strategies to quantify uncertainties using the recovered posterior distribution are developed: (i) local (pixel-wise) credible intervals to provide error bars for each individual pixel; (ii) highest posterior density credible regions; and (iii) hypothesis testing of image structure. These forms of uncertainty quantification provide rich information for analysing radio interferometric observations in a statistically robust manner.

arXiv comment16 pages, 7 figures, see companion article in this arXiv listing

References in corpus (35)

• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• Multi-Scale CLEAN deconvolution of radio synthesis images
• The non-coplanar baselines effect in radio interferometry: The W-Projection algorithm
• Compressed sensing imaging techniques for radio interferometry
• Information field theory for cosmological perturbation reconstruction and non-linear signal analysis
• Correcting direction-dependent gains in the deconvolution of radio interferometric images
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• A Survey of Stochastic Simulation and Optimization Methods in Signal Processing
• PURIFY: a new approach to radio-interferometric imaging
• Proximal Splitting Methods in Signal Processing
• On Variable Density Compressive Sampling
• The application of compressive sampling to radio astronomy I: Deconvolution
• LOFAR Sparse Image Reconstruction
• Scale sensitive deconvolution of interferometric images I. Adaptive Scale Pixel (Asp) decomposition
• Advances in Calibration and Imaging Techniques in Radio Interferometry
• Scalable splitting algorithms for big-data interferometric imaging in the SKA era
• MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators. A sparse deconvolution algorithm for radio interferometric imaging
• Compressed sensing for wide-field radio interferometric imaging
• RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
• Spread spectrum for imaging techniques in radio interferometry
• Robust sparse image reconstruction of radio interferometric observations with purify
• The application of compressive sampling to radio astronomy II: Faraday rotation measure synthesis
• Compressed Sensing with Coherent and Redundant Dictionaries
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• Maximum-entropy image reconstruction using wavelets
• Simulating full-sky interferometric observations
• An accelerated splitting algorithm for radio-interferometric imaging: when natural and uniform weighting meet
• Probabilistic image reconstruction for radio interferometers
• Revisiting the spread spectrum effect in radio interferometric imaging: a sparse variant of the w-projection algorithm
• Disparity and Optical Flow Partitioning Using Extended Potts Priors
• The w-effect in interferometric imaging: from a fast sparse measurement operator to super-resolution
• A multiple-beam CLEAN for imaging intra-day variable radio sources
• Efficient Bayesian computation by proximal Markov chain Monte Carlo: when Langevin meets Moreau
• fastRESOLVE: fast Bayesian imaging for aperture synthesis in radio astronomy
• Online radio interferometric imaging: assimilating and discarding visibilities on arrival

Cited by in corpus (36)

• A Survey on Epistemic (Model) Uncertainty in Supervised Learning: Recent Advances and Applications
• Hybrid Very Long Baseline Interferometry Imaging and Modeling with Themis
• On Solving SAR Imaging Inverse Problems Using Non-Convex Regularization with a Cauchy-based Penalty
• Convergence Guarantees for Non-Convex Optimisation with Cauchy-Based Penalties
• Image reconstruction algorithms in radio interferometry: from handcrafted to learned regularization denoisers
• Comparison of classical and Bayesian imaging in radio interferometry
• Unified Radio Interferometric Calibration and Imaging with Joint Uncertainty Quantification
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• Spherical Bayesian mass-mapping with uncertainties: full sky observations on the celestial sphere
• A D-term Modeling Code (DMC) for simultaneous calibration and full-Stokes imaging of very long baseline interferometric data
• Sparse Bayesian mass-mapping with uncertainties: local credible intervals
• Proximal nested sampling for high-dimensional Bayesian model selection
• Multi-scale and Multi-directional VLBI Imaging with CLEAN
• Deep Probabilistic Imaging: Uncertainty Quantification and Multi-modal Solution Characterization for Computational Imaging
• Bayesian sparse reconstruction: a brute-force approach to astronomical imaging and machine learning
• Proximal Langevin Algorithm: Rapid Convergence Under Isoperimetry
• Sparse Bayesian mass-mapping with uncertainties: peak statistics and feature locations
• Langevin Monte Carlo without smoothness
• Detection of Ship Wakes in SAR Imagery Using Cauchy Regularisation
• Bayesian imaging using Plug & Play priors: when Langevin meets Tweedie
• Posterior sampling for inverse imaging problems on the sphere in seismology and cosmology
• Learned Interferometric Imaging for the SPIDER Instrument
• fast-resolve: Fast Bayesian Radio Interferometric Imaging
• PolyCLEAN: Atomic Optimization for Super-Resolution Imaging and Uncertainty Estimation in Radio Interferometry
• Uncertainty-Aware Blob Detection with an Application to Integrated-Light Stellar Population Recoveries
• A Distributed Block-Split Gibbs Sampler with Hypergraph Structure for High-Dimensional Inverse Problems
• Tomography for Plasma Imaging: a Unifying Framework for Bayesian Inference
• Distribution-free uncertainty quantification for inverse problems: application to weak lensing mass mapping
• Distributed and parallel sparse convex optimization for radio interferometry with PURIFY
• Radio Galaxy Detection in the Visibility Domain
• Online radio interferometric imaging: assimilating and discarding visibilities on arrival
• Offline and online reconstruction for radio interferometric imaging
• Quantifying Uncertainty in High Dimensional Inverse Problems by Convex Optimisation
• DoG-HiT: A novel VLBI Multiscale Imaging Approach
• Fast Posterior Estimation of Cardiac Electrophysiological Model Parameters via Bayesian Active Learning
• Sparse Bayesian mass-mapping with uncertainties: hypothesis testing of structure