Comparison of classical and Bayesian imaging in radio interferometry

Abstract

CLEAN, the commonly employed imaging algorithm in radio interferometry, suffers from a number of shortcomings: in its basic version it does not have the concept of diffuse flux, and the common practice of convolving the CLEAN components with the CLEAN beam erases the potential for super-resolution; it does not output uncertainty information; it produces images with unphysical negative flux regions; and its results are highly dependent on the so-called weighting scheme as well as on any human choice of CLEAN masks to guiding the imaging. Here, we present the Bayesian imaging algorithm resolve which solves the above problems and naturally leads to super-resolution. We take a VLA observation of Cygnus~A at four different frequencies and image it with single-scale CLEAN, multi-scale CLEAN and resolve. Alongside the sky brightness distribution resolve estimates a baseline-dependent correction function for the noise budget, the Bayesian equivalent of weighting schemes. We report noise correction factors between 0.4 and 429. The enhancements achieved by resolve come at the cost of higher computational effort.

arXiv comment23 pages, 16 figures, 4 tables, data published at https://doi.org/10.5281/zenodo.4267057

References in corpus (20)

• 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
• Multi-Scale CLEAN deconvolution of radio synthesis images
• The non-coplanar baselines effect in radio interferometry: The W-Projection algorithm
• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
• Information field theory for cosmological perturbation reconstruction and non-linear signal analysis
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• PURIFY: a new approach to radio-interferometric imaging
• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• The Denoised, Deconvolved, and Decomposed Fermi $\gamma$-ray sky - An application of the D${}^3$PO algorithm
• Reconstruction of signals with unknown spectra in information field theory with parameter uncertainty
• RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
• Information theory for fields
• Unified Radio Interferometric Calibration and Imaging with Joint Uncertainty Quantification
• Cygnus A super-resolved via convex optimisation from VLA data
• Probabilistic image reconstruction for radio interferometers
• A new approach to multi-frequency synthesis in radio interferometry
• Wideband Super-resolution Imaging in Radio Interferometry via Low Rankness and Joint Average Sparsity Models (HyperSARA)
• Discovery of a Luminous Radio Transient 460 pc from the Central Supermassive Black Hole in Cygnus A
• fastRESOLVE: fast Bayesian imaging for aperture synthesis in radio astronomy

Cited by in corpus (26)

• The mass distribution in the Galactic Centre from interferometric astrometry of multiple stellar orbits
• Efficient wide-field radio interferometry response
• Geometric variational inference
• Deep Images of the Galactic Center with GRAVITY
• Data Combination: Interferometry and Single-dish Imaging in Radio Astronomy
• Bayesian radio interferometric imaging with direction-dependent calibration
• Using multiobjective optimization to reconstruct interferometric data (I)
• Multi-scale and Multi-directional VLBI Imaging with CLEAN
• Trajectory Based RFI Subtraction and Calibration for Radio Interferometry
• Dynamic and polarimetric VLBI imaging with a multiscalar approach
• Bayesian self-calibration and imaging in Very Long Baseline Interferometry
• Strong-lensing source reconstruction with variationally optimised Gaussian processes
• Imaging a ring-like structure and the extended jet of M87 at 86 GHz
• Swarm intelligence for full Stokes dynamic imaging reconstruction of interferometric data
• Information Field Theory and Artificial Intelligence
• Multi-Component Imaging of the Fermi Gamma-ray Sky in the Spatio-spectral Domain
• fast-resolve: Fast Bayesian Radio Interferometric Imaging
• Gaussian Process Modelling for Improved Resolution in Faraday Depth Reconstruction
• PolyCLEAN: Atomic Optimization for Super-Resolution Imaging and Uncertainty Estimation in Radio Interferometry
• How to make CLEAN variants faster? Using clustered components informed by the autocorrelation function
• Fast-Cadence High-Contrast Imaging with Information Field Theory
• Mid-infrared extinction toward the Galactic center
• Probing circular polarization and magnetic field structure in AGN
• Toward a Robust R2D2 Paradigm for Radio-interferometric Imaging: Revisiting Deep Neural Network Training and Architecture
• Latent-space Field Tension for Astrophysical Component Detection An application to X-ray imaging
• Accelerating the CLEAN algorithm of radio interferometry with convex optimization