A fast and exact $w$-stacking and $w$-projection hybrid algorithm for wide-field interferometric imaging

Abstract

The standard wide-field imaging technique, the $w$-projection, allows correction for wide-fields of view for non-coplanar radio interferometric arrays. However, calculating exact corrections for each measurement has not been possible due to the amount of computation required at high resolution and with the large number of visibilities from current interferometers. The required accuracy and computational cost of these corrections is one of the largest unsolved challenges facing next generation radio interferometers such as the Square Kilometre Array. We show that the same calculation can be performed with a radially symmetric $w$-projection kernel, where we use one dimensional adaptive quadrature to calculate the resulting Hankel transform, decreasing the computation required for kernel generation by several orders of magnitude, whilst preserving the accuracy. We confirm that the radial $w$-projection kernel is accurate to approximately 1% by imaging the zero-spacing with an added $w$-term. We demonstrate the potential of our radially symmetric $w$-projection kernel via sparse image reconstruction, using the software package PURIFY. We develop a distributed $w$-stacking and $w$-projection hybrid algorithm. We apply this algorithm to individually correct for non-coplanar effects in 17.5 million visibilities over a $25$ by $25$ degree field of view MWA observation for image reconstruction. Such a level of accuracy and scalability is not possible with standard $w$-projection kernel generation methods. This demonstrates that we can scale to a large number of measurements with large image sizes whilst still maintaining both speed and accuracy.

9 Figures, 19 Pages. Accepted to ApJ

References in corpus (11)

• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• Hydrogen Epoch of Reionization Array (HERA)
• The non-coplanar baselines effect in radio interferometry: The W-Projection algorithm
• Correcting direction-dependent gains in the deconvolution of radio interferometric images
• Image Domain Gridding: a fast method for convolutional resampling of visibilities
• Robust sparse image reconstruction of radio interferometric observations with purify
• Simulating full-sky interferometric observations
• Cygnus A super-resolved via convex optimisation from VLA data
• 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 (8)

• High accuracy wide field imaging method in radio interferometry
• Bayesian radio interferometric imaging with direction-dependent calibration
• Cygnus A jointly calibrated and imaged via non-convex optimisation from VLA data
• Parallel implementation of w-projection wide-field imaging
• Wide-band Rotation Measure Synthesis
• Accuracy vs. Complexity: Calibrating radio interferometer arrays with non-homogeneous element patterns
• $w$-stacking $w$-projection hybrid algorithm for wide-field interferometric imaging: implementation details and improvements
• Multi-Step Reconstruction of Radio-Interferometric Images