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.
Cite
@article{arxiv.1309.1469,
title = {Probabilistic image reconstruction for radio interferometers},
author = {P. M. Sutter and Benjamin D. Wandelt and Jason D. McEwen and Emory F. Bunn and Ata Karakci and Andrei Korotkov and Peter Timbie and Gregory S. Tucker and Le Zhang},
journal= {arXiv preprint arXiv:1309.1469},
year = {2015}
}
Comments
12 pages, 10 figures, MNRAS accepted. arXiv admin note: text overlap with arXiv:1109.4640