Radio Galaxy Detection in the Visibility Domain

Abstract

We explore a new Bayesian method of detecting galaxies from radio interferometric data of the faint sky. Working in the Fourier domain, we fit a single, parameterised galaxy model to simulated visibility data of star-forming galaxies. The resulting multimodal posterior distribution is then sampled using a multimodal nested sampling algorithm such as MultiNest. For each galaxy, we construct parameter estimates for the position, flux, scale-length and ellipticities from the posterior samples. We first test our approach on simulated SKA1-MID visibility data of up to 100 galaxies in the field of view, considering a typical weak lensing survey regime (SNR $\ge 10$) where 98% of the input galaxies are detected with no spurious source detections. We then explore the low SNR regime, finding our approach reliable in galaxy detection and providing in particular high accuracy in positional estimates down to SNR $\sim 5$. The presented method does not require transformation of visibilities to the image domain, and requires no prior knowledge of the number of galaxies in the field of view, thus could become a useful tool for constructing accurate radio galaxy catalogs in the future.

arXiv comment11 pages, 11 figures. Accepted for publication in MNRAS

References in corpus (27)

• MultiNest: an efficient and robust Bayesian inference tool for cosmology and particle physics
• Multimodal nested sampling: an efficient and robust alternative to MCMC methods for astronomical data analysis
• X-ray spectral modelling of the AGN obscuring region in the CDFS: Bayesian model selection and catalogue
• Bayes in the sky: Bayesian inference and model selection in cosmology
• Importance Nested Sampling and the MultiNest Algorithm
• Multi-Scale CLEAN deconvolution of radio synthesis images
• PolyChord: nested sampling for cosmology
• PolyChord: next-generation nested sampling
• Bayesian Galaxy Shape Measurement for Weak Lensing Surveys - III. Application to the Canada-France-Hawaii Telescope Lensing Survey
• Shapelets: I. A Method for Image Analysis
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• UVMULTIFIT: A versatile tool for fitting astronomical radio interferometric data
• Shapelets: II. A Method for Weak Lensing Measurements
• Revisiting the radio interferometer measurement equation. II. Calibration and direction-dependent effects
• Scale sensitive deconvolution of interferometric images I. Adaptive Scale Pixel (Asp) decomposition
• A Bayesian approach to discrete object detection in astronomical datasets
• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• Weak Lensing by Large-Scale Structure with the FIRST Radio Survey
• Revisiting the radio interferometer measurement equation. III. Addressing direction-dependent effects in 21 cm WSRT observations of 3C 147
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• On statistical uncertainty in nested sampling
• Bayesian Inference for Radio Observations
• Montblanc: GPU accelerated Radio Interferometer Measurement Equations in support of Bayesian Inference for Radio Observations
• Cluster detection in weak lensing surveys
• Radio Weak Lensing Shear Measurement in the Visibility Domain - I. Methodology
• Radio Weak Lensing Shear Measurement in the Visibility Domain - II. Source Extraction
• Radio Galaxy Shape Measurement with Hamiltonian Monte Carlo in the Visibility Domain