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

Proper evaluation of spatially correlated noise in interferometric images

The spatial correlation problem of noise in imaging deblurring and its solution.

Influence Of Noise On Detection Accuracy Of Joint Transform Optical Correlator

Mean Square Convergence Unbiased Estimation of Thermal Light Correlated Imaging

Spatial Uncertainty of Measurement Errors in Digital Image Correlation

Seismic interferometry from correlated noise sources

Spatially Adaptive Column Fixed-Pattern Noise Correction in Infrared Imaging System Using 1D Horizontal Differential Statistics

Increasing the achievable contrast of infrared interferometry with an error correlation model

High-dimensional Statistical Analysis and Its Application to an ALMA Map of NGC 253

On the variance of radio interferometric calibration solutions: Quality-based Weighting Schemes

Sensitivity of an imaging space infrared interferometer

The autocorrelated noise filtering problem: the ISMC filter in a specific case of distance measuring

A statistical framework for recovering intensity mapping autocorrelations from crosscorrelations

Confidence Level and Sensitivity Limits in High Contrast Imaging

Bias-free model fitting of correlated data in interferometry

High Dimensional Statistical Analysis and its Application to ALMA Map of NGC 253

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

Particle image velocimetry correlation signal-to-noise ratio metrics and measurement uncertainty quantification

Revisiting a flux recovery systematic error arising from common deconvolution methods used in aperture-synthesis imaging

Forecasts and Statistical Insights for Line Intensity Mapping Cross-Correlations: A Case Study with 21cm x [CII]

Efficient modeling of correlated noise I. Statistical significance of periodogram peaks