Fast W-Projection for Wide-field Imaging

Wide-field imaging has become a major challenge for modern radio astronomy, which uses high sensitivity acquisition systems that deal with huge amounts of data. In this paper we investigate a fast wide-field imaging solution based on the w-projection algorithm, which is intended for modern astronomy systems. The core idea of the proposed method is to reduce the computational complexity of the convolution kernel generation step, specifically by replacing the standard two-dimensional FFT by the one-dimensional Hankel transform. Experimental results show that the optimised w-projection proposed here produces equivalent dirty image results in a circular image region, at a significantly lower computational cost than standard $w$-projection. One of the main advantages of the proposed solution is its slow scaling with the number of w-planes, thus enabling more accurate output results at a lower computational cost.