Abstract:Next-generation radio interferometers like the Square Kilometer Array have the potential to unlock scientific discoveries thanks to their unprecedented angular resolution and sensitivity. One key to unlocking their potential resides in handling the deluge and complexity of incoming data. This challenge requires building radio interferometric imaging methods that can cope with the massive data sizes and provide high-quality image reconstructions with uncertainty quantification (UQ). This work proposes a method coined QuantifAI to address UQ in radio-interferometric imaging with data-driven (learned) priors for high-dimensional settings. Our model, rooted in the Bayesian framework, uses a physically motivated model for the likelihood. The model exploits a data-driven convex prior, which can encode complex information learned implicitly from simulations and guarantee the log-concavity of the posterior. We leverage probability concentration phenomena of high-dimensional log-concave posteriors that let us obtain information about the posterior, avoiding MCMC sampling techniques. We rely on convex optimisation methods to compute the MAP estimation, which is known to be faster and better scale with dimension than MCMC sampling strategies. Our method allows us to compute local credible intervals, i.e., Bayesian error bars, and perform hypothesis testing of structure on the reconstructed image. In addition, we propose a novel blazing-fast method to compute pixel-wise uncertainties at different scales. We demonstrate our method by reconstructing radio-interferometric images in a simulated setting and carrying out fast and scalable UQ, which we validate with MCMC sampling. Our method shows an improved image quality and more meaningful uncertainties than the benchmark method based on a sparsity-promoting prior. QuantifAI's source code: <a class="link-external link-https" href="https://github.com/astro-informatics/QuantifAI" rel="external noopener nofollow">this https URL</a>.

PURIFYing real radio interferometric observations

PURIFY: a new approach to radio-interferometric imaging

Scalable splitting algorithms for big-data interferometric imaging in the SKA era

PolyCLEAN: Atomic Optimization for Super-Resolution Imaging and Uncertainty Estimation in Radio Interferometry

Multi-Step Reconstruction of Radio-Interferometric Images

Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations

Africanus III. pfb-imaging -- a flexible radio interferometric imaging suite

Scalable Bayesian uncertainty quantification with data-driven priors for radio interferometric imaging

Scalable precision wide-field imaging in radio interferometry: I. uSARA validated on ASKAP data

PolyCLEAN: When Högbom meets Bayes -- Fast Super-Resolution Imaging with Bayesian MAP Estimation

Deep learning-based deconvolution for interferometric radio transient reconstruction

Learned radio interferometric imaging for varying visibility coverage

Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): II. Code and real data proof of concept

LOFAR Sparse Image Reconstruction

A new approach for short-spacing correction of radio interferometric data sets

Comparison of classical and Bayesian imaging in radio interferometry

Super-resolution Full Polarimetric Imaging for Radio Interferometry with Sparse Modeling

Radio-astronomical Image Reconstruction with Conditional Denoising Diffusion Model

Multi-scale and Multi-directional VLBI Imaging with CLEAN

Towards radio astronomical imaging using an arbitrary basis

Uncertainty quantification for fast reconstruction methods using augmented equivariant bootstrap: Application to radio interferometry