fast-resolve: Fast Bayesian Radio Interferometric Imaging

Abstract

Context: Interferometric imaging is algorithmically and computationally challenging as there is no unique inversion from the measurement data back to the sky maps, and the datasets can be very large. Many imaging methods already exist, but most of them focus either on the accuracy or the computational aspect. Aims: This paper aims to reduce the computational complexity of the Bayesian imaging algorithm resolve, enabling the application of Bayesian imaging for larger datasets. Methods: By combining computational shortcuts of the CLEAN algorithm with the Bayesian imaging algorithm resolve we developed an accurate and fast imaging algorithm which we name fast-resolve. Results: We validate the accuracy of the presented fast-resolve algorithm by comparing it with results from resolve on VLA Cygnus A data. Furthermore, we demonstrate the computational advantages of fast-resolve on a large MeerKAT ESO 137-006 dataset which is computationally out of reach for resolve. Conclusions: The presented algorithm is significantly faster than previous Bayesian imaging algorithms, broadening the applicability of Bayesian interferometric imaging. Specifically for the single channel VLA Cygnus A datasets fast-resolve is about $144$ times faster than resolve. For the MeerKAT dataset with multiple channels the computational speedup of fast-resolve is even larger.

24 pages, 14 figures

References in corpus (37)

• Array Programming with NumPy
• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• A multi-scale multi-frequency deconvolution algorithm for synthesis imaging in radio interferometry
• The Expanded Very Large Array -- a New Telescope for New Science
• The non-coplanar baselines effect in radio interferometry: The W-Projection algorithm
• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
• Compressed sensing imaging techniques for radio interferometry
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• A parsec-scale Galactic 3D dust map out to 1.25 kpc from the Sun
• CubiCal - Fast radio interferometric calibration suite exploiting complex optimisation
• Scale sensitive deconvolution of interferometric images I. Adaptive Scale Pixel (Asp) decomposition
• Advances in Calibration and Imaging Techniques in Radio Interferometry
• NIFTY - Numerical Information Field Theory - a versatile Python library for signal inference
• Variable structures in M87* from space, time and frequency resolved interferometry
• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• Collimated synchrotron threads linking the radio lobes of ESO137-006
• RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
• Image reconstruction algorithms in radio interferometry: from handcrafted to learned regularization denoisers
• Comparison of classical and Bayesian imaging in radio interferometry
• Efficient wide-field radio interferometry response
• Unified Radio Interferometric Calibration and Imaging with Joint Uncertainty Quantification
• Deep Learning-based Imaging in Radio Interferometry
• Cygnus A super-resolved via convex optimisation from VLA data
• First AI for deep super-resolution wide-field imaging in radio astronomy: unveiling structure in ESO 137--006
• Deep Images of the Galactic Center with GRAVITY
• Probabilistic image reconstruction for radio interferometers
• Optimal Image Reconstruction in Radio Interferometry
• Re-Envisioning Numerical Information Field Theory (NIFTy.re): A Library for Gaussian Processes and Variational Inference
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations
• Bayesian radio interferometric imaging with direction-dependent calibration
• Wideband Super-resolution Imaging in Radio Interferometry via Low Rankness and Joint Average Sparsity Models (HyperSARA)
• A Wideband Polarization Study of Cygnus A with the JVLA. I: The Observations and Data
• Sparse interferometric Stokes imaging under polarization constraint (Polarized SARA)
• Polca SARA - Full polarization, direction-dependent calibration and sparse imaging for radio interferometry
• Cygnus A jointly calibrated and imaged via non-convex optimisation from VLA data
• Bayesian Reasoning with Trained Neural Networks
• JAXbind: Bind any function to JAX

Cited by in corpus (3)

• Accelerating the CLEAN algorithm of radio interferometry with convex optimization
• Latent-space Field Tension for Astrophysical Component Detection An application to X-ray imaging
• Very-Long Baseline Interferometry Imaging with Closure Invariants using Conditional Image Diffusion