Probabilistic image reconstruction for radio interferometers

Abstract

We present a novel, general-purpose method for deconvolving and denoising images from gridded radio interferometric visibilities using Bayesian inference based on a Gaussian process model. The method automatically takes into account incomplete coverage of the uv-plane, signal mode coupling due to the primary beam, and noise mode coupling due to uv sampling. Our method uses Gibbs sampling to efficiently explore the full posterior distribution of the underlying signal image given the data. We use a set of widely diverse mock images with a realistic interferometer setup and level of noise to assess the method. Compared to results from a proxy for point source- based CLEAN method we find that in terms of RMS error and signal-to-noise ratio our approach performs better than traditional deconvolution techniques, regardless of the structure of the source image in our test suite. Our implementation scales as O(np log np), provides full statistical and uncertainty information of the reconstructed image, requires no supervision, and provides a robust, consistent framework for incorporating noise and parameter marginalizations and foreground removal.

12 pages, 10 figures, MNRAS accepted. arXiv admin note: text overlap with arXiv:1109.4640

References in corpus (27)

• A multi-scale multi-frequency deconvolution algorithm for synthesis imaging in radio interferometry
• Multi-Scale CLEAN deconvolution of radio synthesis images
• The Expanded Very Large Array -- a New Telescope for New Science
• The Precision Array for Probing the Epoch of Reionization: 8 Station Results
• Science with ASKAP - the Australian Square Kilometre Array Pathfinder
• The Murchison Widefield Array: Design Overview
• Compressed sensing imaging techniques for radio interferometry
• Molecular and Atomic Gas in the Local Group Galaxy M33
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• Application of Monte Carlo Algorithms to the Bayesian Analysis of the Cosmic Microwave Background
• The application of compressive sampling to radio astronomy I: Deconvolution
• Scale sensitive deconvolution of interferometric images I. Adaptive Scale Pixel (Asp) decomposition
• Radio Continuum Surveys with Square Kilometre Array Pathfinders
• Unveiling radio halos in galaxy clusters in the LOFAR era
• A Fast Gridded Method for the Estimation of the Power Spectrum of the CMB from Interferometer Data with Application to the Cosmic Background Imager
• Sparsity Averaging for Compressive Imaging
• Compressed sensing for wide-field radio interferometric imaging
• Positive deconvolution for superimposed extended source and point sources
• The ATCA HI Galactic Center Survey
• HI Power Spectra and the Turbulent ISM of Dwarf Irregular Galaxies
• Deep Radio Continuum Imaging of the Dwarf Irregular Galaxy IC10: Tracing Star Formation and Magnetic Fields
• Compressed sensing reconstruction of a string signal from interferometric observations of the cosmic microwave background
• Optimal Image Reconstruction in Radio Interferometry
• Revisiting the spread spectrum effect in radio interferometric imaging: a sparse variant of the w-projection algorithm
• Bayesian angular power spectrum analysis of interferometric data
• A multiple-beam CLEAN for imaging intra-day variable radio sources
• Radio interferometric imaging of spatial structure that varies with time and frequency

Cited by in corpus (19)

• Empirical Covariance Modeling for 21 cm Power Spectrum Estimation: A Method Demonstration and New Limits from Early Murchison Widefield Array 128-Tile Data
• LOFAR Sparse Image Reconstruction
• 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
• RESOLVE: A new algorithm for aperture synthesis imaging of extended emission in radio astronomy
• 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
• Mapmaking for Precision 21 cm Cosmology
• Comparison of classical and Bayesian imaging in radio interferometry
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• Matrix-free Large Scale Bayesian inference in cosmology
• Sparse representations and convex optimization as tools for LOFAR radio interferometric imaging
• Bayesian Inference for Radio Observations
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations
• Detecting Unspecified Structure in Low-Count Images
• High resolution VLBI polarisation imaging of AGN with the Maximum Entropy Method
• fast-resolve: Fast Bayesian Radio Interferometric Imaging
• SSSpaNG! Stellar Spectra as Sparse, data-driven, Non-Gaussian processes
• It's Always Darkest Before the Cosmic Dawn: Early Results from Novel Tools and Telescopes for 21 cm Cosmology