Bayesian Inference for Radio Observations

Abstract

New telescopes like the Square Kilometre Array (SKA) will push into a new sensitivity regime and expose systematics, such as direction-dependent effects, that could previously be ignored. Current methods for handling such systematics rely on alternating best estimates of instrumental calibration and models of the underlying sky, which can lead to inadequate uncertainty estimates and biased results because any correlations between parameters are ignored. These deconvolution algorithms produce a single image that is assumed to be a true representation of the sky, when in fact it is just one realization of an infinite ensemble of images compatible with the noise in the data. In contrast, here we report a Bayesian formalism that simultaneously infers both systematics and science. Our technique, Bayesian Inference for Radio Observations (BIRO), determines all parameters directly from the raw data, bypassing image-making entirely, by sampling from the joint posterior probability distribution. This enables it to derive both correlations and accurate uncertainties, making use of the flexible software MEQTREES to model the sky and telescope simultaneously. We demonstrate BIRO with two simulated sets of Westerbork Synthesis Radio Telescope data sets. In the first, we perform joint estimates of 103 scientific (flux densities of sources) and instrumental (pointing errors, beamwidth and noise) parameters. In the second example, we perform source separation with BIRO. Using the Bayesian evidence, we can accurately select between a single point source, two point sources and an extended Gaussian source, allowing for 'super-resolution' on scales much smaller than the synthesized beam.

Published in MNRAS. See https://vimeo.com/117391380 for a video of MultiNest converging to the correct source model

References in corpus (30)

• emcee: The MCMC Hammer
• MCMC using Hamiltonian dynamics
• MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
• LOFAR: The LOw-Frequency ARray
• Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
• Bayes in the sky: Bayesian inference and model selection in cosmology
• Science with the Square Kilometer Array: Motivation, Key Science Projects, Standards and Assumptions
• Compressed sensing imaging techniques for radio interferometry
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• Correcting direction-dependent gains in the deconvolution of radio interferometric images
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• Radio Interferometric Calibration Using The SAGE Algorithm
• Revisiting the radio interferometer measurement equation. II. Calibration and direction-dependent effects
• The application of compressive sampling to radio astronomy I: Deconvolution
• The MeqTrees software system and its use for third-generation calibration of radio interferometers
• Over-resolution of compact sources in interferometric observations
• Bayesian inference from photometric redshift surveys
• Revisiting the radio interferometer measurement equation. III. Addressing direction-dependent effects in 21 cm WSRT observations of 3C 147
• Bayesian modelling of clusters of galaxies from multi-frequency pointed Sunyaev--Zel'dovich observations
• Robust Radio Interferometric Calibration Using the t-Distribution
• Very Small Array observations of the Sunyaev-Zel'dovich effect in nearby galaxy clusters
• Probabilistic image reconstruction for radio interferometers
• Optimal Image Reconstruction in Radio Interferometry
• A new approach to multi-frequency synthesis in radio interferometry
• Sunyaev-Zel'dovich observations of galaxy clusters out to the virial radius with the Arcminute Microkelvin Imager
• Bayesian angular power spectrum analysis of interferometric data
• A blind detection of a large, complex, Sunyaev--Zel'dovich structure
• Improving self-calibration
• Signal inference with unknown response: Calibration-uncertainty renormalized estimator
• MeqTrees and direction-dependent effects

Cited by in corpus (15)

• 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
• DeepSource: Point Source Detection using Deep Learning
• Separating the EoR Signal with a Convolutional Denoising Autoencoder: A Deep-learning-based Method
• Multi-GPU maximum entropy image synthesis for radio astronomy
• Montblanc: GPU accelerated Radio Interferometer Measurement Equations in support of Bayesian Inference for Radio Observations
• Trajectory Based RFI Subtraction and Calibration for Radio Interferometry
• Resolving the blazar CGRaBS J0809+5341 in the presence of telescope systematics
• Far beyond stacking: Fully bayesian constraints on sub-microJy radio source populations over the XMM-LSS-VIDEO field
• Radio Astronomical Image Formation using Constrained Least Squares and Krylov Subspaces
• Radio Galaxy Shape Measurement with Hamiltonian Monte Carlo in the Visibility Domain
• Statistical estimation of full-sky radio maps from 21cm array visibility data using Gaussian constrained realizations
• A probabilistic approach to phase calibration: I. Effects of source structure on fringe-fitting
• Bayesian Source Discrimination in Radio Interferometry
• Radio Galaxy Detection in the Visibility Domain