Uncertainty quantification for radio interferometric imaging: II. MAP estimation

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. 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, for massive data sizes, like those anticipated from the Square Kilometre Array (SKA), it will be difficult if not impossible to apply any MCMC technique due to its inherent computational cost. We formulate Bayesian inference problems with sparsity-promoting priors (motivated by compressive sensing), for which we recover maximum a posteriori (MAP) point estimators of radio interferometric images by convex optimisation. Exploiting recent developments in the theory of probability concentration, we quantify uncertainties by post-processing the recovered MAP estimate. Three strategies to quantify uncertainties are developed: (i) highest posterior density credible regions; (ii) local credible intervals (cf. error bars) for individual pixels and superpixels; 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. Our MAP-based methods are approximately $10^5$ times faster computationally than state-of-the-art MCMC methods and, in addition, support highly distributed and parallelised algorithmic structures. For the first time, our MAP-based techniques provide a means of quantifying uncertainties for radio interferometric imaging for realistic data volumes and practical use, and scale to the emerging big-data era of radio astronomy.

arXiv comment13 pages, 10 figures, see companion article in this arXiv listing

References in corpus (20)

• Compressed sensing imaging techniques for radio interferometry
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• PURIFY: a new approach to radio-interferometric imaging
• 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
• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• 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
• 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
• 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
• fastRESOLVE: fast Bayesian imaging for aperture synthesis in radio astronomy

Cited by in corpus (24)

• A Survey on Epistemic (Model) Uncertainty in Supervised Learning: Recent Advances and Applications
• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• Hybrid Very Long Baseline Interferometry Imaging and Modeling with Themis
• 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
• Free-moving Quantitative Gamma-ray Imaging
• Sparse Bayesian mass-mapping with uncertainties: local credible intervals
• Multi-scale and Multi-directional VLBI Imaging with CLEAN
• Proximal Langevin Algorithm: Rapid Convergence Under Isoperimetry
• Sparse Bayesian mass-mapping with uncertainties: peak statistics and feature locations
• Posterior sampling for inverse imaging problems on the sphere in seismology and cosmology
• Learned Interferometric Imaging for the SPIDER Instrument
• 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
• Sparse image reconstruction on the sphere: a general approach with uncertainty quantification
• Accelerating proximal Markov chain Monte Carlo by using an explicit stabilised method
• Distribution-free uncertainty quantification for inverse problems: application to weak lensing mass mapping
• Radio Galaxy Detection in the Visibility Domain
• Distributed and parallel sparse convex optimization for radio interferometry with PURIFY
• An Unsupervised Learning Method for Radio Interferometry Deconvolution
• Bayesian variational regularization on the ball
• DoG-HiT: A novel VLBI Multiscale Imaging Approach
• Quantifying Uncertainty in High Dimensional Inverse Problems by Convex Optimisation
• Offline and online reconstruction for radio interferometric imaging