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>.

BIPP: An efficient HPC implementation of the Bluebild algorithm for radio astronomy

Bifrost: a Python/C++ Framework for High-Throughput Stream Processing in Astronomy

CRBLASTER: A Parallel-Processing Computational Framework for Embarrassingly-Parallel Image-Analysis Algorithms

A fast and portable Re-Implementation of Piskunov and Valenti's Optimal-Extraction Algorithm with improved Cosmic-Ray Removal and Optimal Sky Subtraction

HEALPix Alchemy: Fast All-Sky Geometry and Image Arithmetic in a Relational Database for Multimessenger Astronomy Brokers

HVOX: Scalable Interferometric Synthesis and Analysis of Spherical Sky Maps

J-UBIK: The JAX-accelerated Universal Bayesian Imaging Kit

High-time resolution GPU imager for FRB searches at low radio frequencies

IOP4, the Interactive Optical Photo-Polarimetric Python Pipeline

PolyBin3D: A Suite of Optimal and Efficient Power Spectrum and Bispectrum Estimators for Large-Scale Structure

SBI++: Flexible, Ultra-fast Likelihood-free Inference Customized for Astronomical Applications

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

A blind deconvolution method for ground based telescopes and Fizeau interferometers

CMB-PAInT: An inpainting tool for the cosmic microwave background

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

The Blending ToolKit: A simulation framework for evaluation of galaxy detection and deblending

Spectral Imaging with QUBIC: building astrophysical components from Time-Ordered-Data using Bolometric Interferometry

BLAST: Beyond Limber Angular power Spectra Toolkit. A fast and efficient algorithm for 3x2 pt analysis

An Integrated System at the Bleien Observatory for Mapping the Galaxy