w-snapshot algorithm

WSClean provides an algorithm that has approximately the same performance as the W-snapshot algorithm suggested by Cornwell et al. (2012), although the algorithm in WSClean is slightly different. In WSClean, it consists of phase rotating the visibilities to zenith, and then recentring the image during w-stacking. Mathematical details are explained in the WSClean paper. It seems to be worthwhile for MWA snapshots of a few minutes with images of 3072 x 3072 pixels at zenith angles > 20 degree, and provides a speed-up of about a factor of 3 at zenith angles > 45 degrees. For larger images the speed-up will be greater. It is not useful for long integrations.

To use this method, the measurement set should be phase rotated to the phase centre that you want to image (as it normally would be). This chapter assumes this is 08h20m00.0s -42d45m00s (which is Puppis A). To prepare for the w-snapshotting approach, first run the chgcentre tool with the following parameters:

chgcentre -minw copy.ms

This will calculate the optimal projection direction and perform the required phase rotations. Finally, during imaging the phases should be shifted to the original target, e.g.:

wsclean -shift 08h20m00.0s -42d45m00s ...

The chgcentre command rotates the coordinate system whereas WSClean’s -shift parameter shifts the coordinate system along the tangent plane. The result is that imaging is done in a different projection that is still centered on the target location, but that has smaller w-terms. The chgcentre page provides further information about the chgcentre tool.

The -shift parameter is new in WSClean 3.0. Before 3.0, this could be achieved by using the -shiftback option of chgcentre. That approach is no longer supported in WSClean.

The idg, w-gridding gridder and w-stacking gridders all support this approach. Support for the w-gridding gridder was added in WSClean 3.0.