The w-effect in interferometric imaging: from a fast sparse measurement operator to super-resolution

Abstract

Modern radio telescopes, such as the Square Kilometre Array (SKA), will probe the radio sky over large fields-of-view, which results in large w-modulations of the sky image. This effect complicates the relationship between the measured visibilities and the image under scrutiny. In algorithmic terms, it gives rise to massive memory and computational time requirements. Yet, it can be a blessing in terms of reconstruction quality of the sky image. In recent years, several works have shown that large w-modulations promote the spread spectrum effect. Within the compressive sensing framework, this effect increases the incoherence between the sensing basis and the sparsity basis of the signal to be recovered, leading to better estimation of the sky image. In this article, we revisit the w-projection approach using convex optimisation in realistic settings, where the measurement operator couples the w-terms in Fourier and the de-gridding kernels. We provide sparse, thus fast, models of the Fourier part of the measurement operator through adaptive sparsification procedures. Consequently, memory requirements and computational cost are significantly alleviated, at the expense of introducing errors on the radio-interferometric data model. We present a first investigation of the impact of the sparse variants of the measurement operator on the image reconstruction quality. We finally analyse the interesting super-resolution potential associated with the spread spectrum effect of the w-modulation, and showcase it through simulations. Our C++ code is available online on GitHub.

accepted for publication in MNRAS

References in corpus (11)

• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• The non-coplanar baselines effect in radio interferometry: The W-Projection algorithm
• Compressed sensing imaging techniques for radio interferometry
• Sparsity Averaging Reweighted Analysis (SARA): a novel algorithm for radio-interferometric imaging
• LOFAR Sparse Image Reconstruction
• MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators. A sparse deconvolution algorithm for radio interferometric imaging
• Compressed sensing for wide-field radio interferometric imaging
• Robust sparse image reconstruction of radio interferometric observations with purify
• Simulating full-sky interferometric observations
• An accelerated splitting algorithm for radio-interferometric imaging: when natural and uniform weighting meet
• An improved method for polarimetric image restoration in interferometry

Cited by in corpus (12)

• Uncertainty quantification for radio interferometric imaging: I. proximal MCMC methods
• Image reconstruction algorithms in radio interferometry: from handcrafted to learned regularization denoisers
• Uncertainty quantification for radio interferometric imaging: II. MAP estimation
• First AI for deep super-resolution wide-field imaging in radio astronomy: unveiling structure in ESO 137--006
• Cygnus A super-resolved via convex optimisation from VLA data
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): I. Algorithm and simulations
• Sparse interferometric Stokes imaging under polarization constraint (Polarized SARA)
• Scalable precision wide-field imaging in radio interferometry: I. uSARA validated on ASKAP data
• Cygnus A jointly calibrated and imaged via non-convex optimisation from VLA data
• Parallel faceted imaging in radio interferometry via proximal splitting (Faceted HyperSARA): II. Code and real data proof of concept
• $w$-stacking $w$-projection hybrid algorithm for wide-field interferometric imaging: implementation details and improvements
• Online radio interferometric imaging: assimilating and discarding visibilities on arrival