Robust sparse image reconstruction of radio interferometric observations with purify

Abstract

Next-generation radio interferometers, such as the Square Kilometre Array (SKA), will revolutionise our understanding of the universe through their unprecedented sensitivity and resolution. However, to realise these goals significant challenges in image and data processing need to be overcome. The standard methods in radio interferometry for reconstructing images, such as CLEAN, have served the community well over the last few decades and have survived largely because they are pragmatic. However, they produce reconstructed inter\-ferometric images that are limited in quality and scalability for big data. In this work we apply and evaluate alternative interferometric reconstruction methods that make use of state-of-the-art sparse image reconstruction algorithms motivated by compressive sensing, which have been implemented in the PURIFY software package. In particular, we implement and apply the proximal alternating direction method of multipliers (P-ADMM) algorithm presented in a recent article. First, we assess the impact of the interpolation kernel used to perform gridding and degridding on sparse image reconstruction. We find that the Kaiser-Bessel interpolation kernel performs as well as prolate spheroidal wave functions, while providing a computational saving and an analytic form. Second, we apply PURIFY to real interferometric observations from the Very Large Array (VLA) and the Australia Telescope Compact Array (ATCA) and find images recovered by PURIFY are higher quality than those recovered by CLEAN. Third, we discuss how PURIFY reconstructions exhibit additional advantages over those recovered by CLEAN. The latest version of PURIFY, with developments presented in this work, is made publicly available.

22 pages, 10 figures, PURIFY code available at http://basp-group.github.io/purify

References in corpus (21)

• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• GaLactic and Extragalactic All-sky Murchison Widefield Array (GLEAM) survey I: A low-frequency extragalactic catalogue
• Multi-Scale CLEAN deconvolution of radio synthesis images
• Compressed sensing imaging techniques for radio interferometry
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• LOFAR Sparse Image Reconstruction
• The Australian Square Kilometre Array Pathfinder: System Architecture and Specifications of the Boolardy Engineering Test Array
• The Murchison Widefield Array 21 cm Power Spectrum Analysis Methodology
• Comparison of algorithms for determination of rotation measure and Faraday structure I. 1100 - 1400 MHz
• MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators. A sparse deconvolution algorithm for radio interferometric imaging
• Compressed sensing for wide-field radio interferometric imaging
• A generalised Measurement Equation and van Cittert-Zernike theorem for wide-field radio astronomical interferometry
• Limits on low frequency radio emission from southern exoplanets with the Murchison Widefield Array
• The first Murchison Widefield Array low frequency radio observations of cluster scale non-thermal emission: the case of Abell 3667
• Bent-Double Radio Sources as Probes of Intergalactic Gas
• Simulating full-sky interferometric observations
• Efficient implementation of the adaptive scale pixel decomposition algorithm
• ATCA observations of the MACS-Planck Radio Halo Cluster Project - I. New detection of a radio halo in PLCK G285.0-23.7
• Revisiting the radio interferometer measurement equation. IV. A generalized tensor formalism
• An improved method for polarimetric image restoration in interferometry
• Generalized Formalisms of the Radio Interferometer Measurement Equation

Cited by in corpus (26)

• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• Image reconstruction algorithms in radio interferometry: from handcrafted to learned regularization denoisers
• Deep Radio Interferometric Imaging with POLISH: DSA-2000 and weak lensing
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• Cygnus A super-resolved via convex optimisation from VLA data
• Spherical Bayesian mass-mapping with uncertainties: full sky observations on the celestial sphere
• LOFAR 150-MHz observations of SS 433 and W 50
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations
• SLITronomy: towards a fully wavelet-based strong lensing inversion technique
• Wideband Super-resolution Imaging in Radio Interferometry via Low Rankness and Joint Average Sparsity Models (HyperSARA)
• Sparse Lens Inversion Technique (SLIT): lens and source separability from linear inversion of the source reconstruction problem
• Polca SARA - Full polarization, direction-dependent calibration and sparse imaging for radio interferometry
• Scalable precision wide-field imaging in radio interferometry: I. uSARA validated on ASKAP data
• A fast and exact $w$-stacking and $w$-projection hybrid algorithm for wide-field interferometric imaging
• The w-effect in interferometric imaging: from a fast sparse measurement operator to super-resolution
• RadioGAN - Translations between different radio surveys with generative adversarial networks
• Posterior sampling for inverse imaging problems on the sphere in seismology and cosmology
• Wide-band Rotation Measure Synthesis
• Learned Interferometric Imaging for the SPIDER Instrument
• $w$-stacking $w$-projection hybrid algorithm for wide-field interferometric imaging: implementation details and improvements
• Distributed and parallel sparse convex optimization for radio interferometry with PURIFY
• Bayesian variational regularization on the ball
• CS-ROMER: A novel compressed sensing framework for Faraday depth reconstruction
• DoG-HiT: A novel VLBI Multiscale Imaging Approach
• Quantifying Uncertainty in High Dimensional Inverse Problems by Convex Optimisation
• Offline and online reconstruction for radio interferometric imaging