fastRESOLVE: fast Bayesian imaging for aperture synthesis in radio astronomy

Abstract

The standard imaging algorithm for interferometric radio data, CLEAN, is optimal for point source observations, but suboptimal for diffuse emission. Recently, RESOLVE, a new Bayesian algorithm has been developed, which is ideal for extended source imaging. Unfortunately, RESOLVE is computationally very expensive. In this paper we present fastRESOLVE, a modification of RESOLVE based on an approximation of the interferometric likelihood that allows us to avoid expensive gridding routines and consequently gain a factor of roughly 100 in computation time. Furthermore, we include a Bayesian estimation of the measurement uncertainty of the visibilities into the imaging, a procedure not applied in aperture synthesis before. The algorithm requires little to no user input compared to the standard method CLEAN while being superior for extended and faint emission. We apply the algorithm to VLA data of Abell 2199 and show that it resolves more detailed structures.

16 pages, 8 figures

Cited by in corpus (6)

• 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
• Comparison of classical and Bayesian imaging in radio interferometry
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation