WSCLEAN: an implementation of a fast, generic wide-field imager for radio astronomy

WSCLEAN: an implementation of a fast, generic wide-field imager for radio astronomy

Astronomical wide-field imaging of interferometric radio data is computationally expensive, especially for the large data volumes created by modern non-coplanar many-element arrays. We present a new wide-field interferometric imager that uses the w-stacking algorithm and can make use of the w-snapsh...

Full description

Bibliographic Details
Main Authors: Offringa, A., McKinley, B., Hurley-Walker, Natasha, Briggs, F., Wayth, Randall, Kaplan, D., Bell, M., Feng, L., Neben, A., Hughes, J., Rhee, J., Murphy, T., Bhat, Ramesh, Bernardi, G., Bowman, J., Cappallo, R., Corey, B., Deshpande, A., Emrich, David, Ewall-Wice, A., Gaensler, B., Goeke, R., Greenhill, L., Hazelton, B., Hindson, L., Johnston-Hollitt, M., Jacobs, D., Kasper, J., Kratzenberg, E., Lenc, E., Lonsdale, C., Lynch, M., McWhirter, S., Mitchell, D., Morales, M., Morgan, E., Kudryavtseva, N., Oberoi, D., Ord, S., Pindor, B., Procopio, P., Prabu, T., Riding, J., Roshi, A., Shankar, N., Srivani, K., Subrahmanyan, R., Tingay, Steven, Waterson, Mark, Webster, R., Whitney, A., Williams, A., Williams, C.
Format: Journal Article
Published: Oxford University Press 2014
Online Access:http://hdl.handle.net/20.500.11937/30694
Description
Summary:Astronomical wide-field imaging of interferometric radio data is computationally expensive, especially for the large data volumes created by modern non-coplanar many-element arrays. We present a new wide-field interferometric imager that uses the w-stacking algorithm and can make use of the w-snapshot algorithm. The performance dependences of CASA's w-projection and our new imager are analysed and analytical functions are derived that describe the required computing cost for both imagers. On data from the Murchison Widefield Array, we find our new method to be an order of magnitude faster than w-projection, as well as being capable of full-sky imaging at full resolution and with correct polarization correction. We predict the computing costs for several other arrays and estimate that our imager is a factor of 2–12 faster, depending on the array configuration. We estimate the computing cost for imaging the low-frequency Square Kilometre Array observations to be 60 PetaFLOPS with current techniques. We find that combining w-stacking with the w-snapshot algorithm does not significantly improve computing requirements over pure w-stacking. The source code of our new imager is publicly released.

Similar Items