RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy

Abstract

We present RESOLVE, a new algorithm for radio aperture synthesis imaging of extended and diffuse emission in total intensity. The algorithm is derived using Bayesian statistical inference techniques, estimating the surface brightness in the sky assuming a priori log-normal statistics. RESOLVE not only estimates the measured sky brightness in total intensity, but also its spatial correlation structure, which is used to guide the algorithm to an optimal reconstruction of extended and diffuse sources. For a radio interferometer, it succeeds in deconvolving the effects of the instrumental point spread function during this process. Additionally, RESOLVE provides a map with an uncertainty estimate of the reconstructed surface brightness. Furthermore, with RESOLVE we introduce a new, optimal visibility weighting scheme that can be viewed as an extension to robust weighting. In tests using simulated observations, the algorithm shows improved performance against two standard imaging approaches for extended sources, Multiscale-CLEAN and the Maximum Entropy Method.

22 pages, 8 figures, submitted to Astronomy & Astrophysics

References in corpus (11)

• Multi-Scale CLEAN deconvolution of radio synthesis images
• The non-coplanar baselines effect in radio interferometry: The W-Projection algorithm
• Compressed sensing imaging techniques for radio interferometry
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• Deep 1.4 GHz VLA Observations of the Radio Halo and Relic in Abell 2256
• Revisiting the radio interferometer measurement equation. II. Calibration and direction-dependent effects
• Reconstruction of Gaussian and log-normal fields with spectral smoothness
• Reconstruction of signals with unknown spectra in information field theory with parameter uncertainty
• Log-transforming the matter power spectrum
• Bayesian Inference of Polarized CMB Power Spectra from Interferometric Data

Cited by in corpus (44)

• CASA, the Common Astronomy Software Applications for Radio Astronomy
• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
• Frankenstein: Protoplanetary disc brightness profile reconstruction at sub-beam resolution with a rapid Gaussian process
• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• IMAGINE: A comprehensive view of the interstellar medium, Galactic magnetic fields and cosmic rays
• Hybrid Very Long Baseline Interferometry Imaging and Modeling with Themis
• Resolving the extended atmosphere and the inner wind of Mira ($o$ Ceti) with long ALMA baselines
• Information theory for fields
• Comparison of classical and Bayesian imaging in radio interferometry
• A novel approach to visibility-space modelling of interferometric gravitational lens observations at high angular resolution
• Unified Radio Interferometric Calibration and Imaging with Joint Uncertainty Quantification
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• A new approach to multi-frequency synthesis in radio interferometry
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations
• Wideband Super-resolution Imaging in Radio Interferometry via Low Rankness and Joint Average Sparsity Models (HyperSARA)
• Multi-scale and Multi-directional VLBI Imaging with CLEAN
• Search for quasi-periodic signals in magnetar giant flares
• A joint deconvolution algorithm to combine single dish and interferometer data for wideband multi-term and mosaic imaging
• High-resolution ALMA observations of compact discs in the wide-binary system Sz 65 and Sz 66
• Bayesian self-calibration and imaging in Very Long Baseline Interferometry
• Radio Astronomical Image Formation using Constrained Least Squares and Krylov Subspaces
• Imaging a ring-like structure and the extended jet of M87 at 86 GHz
• Tomography of the Galactic free electron density with the Square Kilometer Array
• High resolution VLBI polarisation imaging of AGN with the Maximum Entropy Method
• Denoising, deconvolving and decomposing multi-domain photon observations- The D4PO algorithm
• Noisy independent component analysis of auto-correlated components
• Multi-frequency analysis of the ALMA and VLA high resolution continuum observations of the substructured disc around CI Tau. Preference for sub-mm-sized low-porosity amorphous carbon grains
• fast-resolve: Fast Bayesian Radio Interferometric Imaging
• Gaussian Process Modelling for Improved Resolution in Faraday Depth Reconstruction
• Radio Imaging With Information Field Theory
• PolyCLEAN: Atomic Optimization for Super-Resolution Imaging and Uncertainty Estimation in Radio Interferometry
• A probabilistic approach to phase calibration: I. Effects of source structure on fringe-fitting
• How to make CLEAN variants faster? Using clustered components informed by the autocorrelation function
• D2O - a distributed data object for parallel high-performance computing in Python
• Diffuse radio sky models using large scale shapelets
• Estimating statistics of sky brightness using radio interferometric observations
• Exploring the MeV Sky with a Combined Coded Mask and Compton Telescope: The Galactic Explorer with a Coded Aperture Mask Compton Telescope (GECCO)
• Towards Bayesian Data Compression
• Enhancing Compton telescope imaging with maximum a posteriori estimation: a modified Richardson-Lucy algorithm for the Compton Spectrometer and Imager
• Convolutional Deep Denoising Autoencoders for Radio Astronomical Images
• Fast algorithms to approximate the position-dependent point spread function responses in radio interferometric wide-field imaging
• Accelerating the CLEAN algorithm of radio interferometry with convex optimization
• Prospects for radio weak lensing: studies using LOFAR observations in the ELAIS-N1 field
• Photometry on Structured Backgrounds: Local Pixelwise Infilling by Regression