WProjectFT.h

Classes

WProjectFT -- An FTMachine for Gridded Fourier transforms (full description)

class WProjectFT : public FTMachine

Interface

Public Members
WProjectFT(MeasurementSet& ms, Int nFacets, Long cachesize, Int tilesize=16, Bool usezero=True)
WProjectFT(MeasurementSet& ms, Int nFacets, MDirection mTangent, MPosition mLocation, Long cachesize, Int tilesize=16, Bool usezero=True, Float padding=1.0)
WProjectFT(const RecordInterface& stateRec)
WProjectFT(const WProjectFT &other)
WProjectFT &operator=(const WProjectFT &other)
~WProjectFT()
void initializeToVis(ImageInterface<Complex>& image, const Complex& vb)
void initializeToVis(ImageInterface<Complex>& image, const Complex& vb, VisBuffer<Complex>& griddedVis, Array<Double>& uvscale)
void finalizeToVis()
void initializeToSky(ImageInterface<Complex>& image, Complex<Float>& weight, const Matrix& vb)
void finalizeToSky()
void get(VisBuffer& vb, Int row=-1)
void get(VisBuffer& vb, Cube<Complex>& degrid, Array<Complex>& griddedVis, Array<Double>& scale, Int row=-1)
void put(const VisBuffer& vb, Int row=-1, Bool dopsf=False, FTMachine::Type type=FTMachine::OBSERVED)
void makeImage(FTMachine::Type type, VisSet& vs, ImageInterface<Complex>& image, Complex<Float>& weight)
ImageInterface<Complex>& getImage(Complex<Float>&, Bool normalize=True)
void getWeightImage(ImageInterface<Float>&, Matrix<Float>&)
Bool toRecord(String& error, RecordInterface& outRec, Bool withImage=False)
Bool fromRecord(String& error, const RecordInterface& inRec)
Bool isFourier()
String name()
Protected Members
Int nint(Double val)
void findConvFunction(const ImageInterface<Complex>& image, const Complex& vb)
Array<Complex>* getDataPointer(const IPosition&, Bool)
void ok()
void init()
Bool recordOnGrid(const VisBuffer& vb, Int rownr) const
Int getIndex(const ROMSPointingColumns& mspc, const Double& time, const Double& interval)
Bool getXYPos(const VisBuffer& vb, Int row)
Bool checkCenterPix(const ImageInterface<Complex>& image)

Description

Prerequisite

Etymology

FTMachine is a Machine for Fourier Transforms. WProjectFT does Grid-based Fourier transforms.

Synopsis

The SkyEquation needs to be able to perform Fourier transforms on visibility data. WProjectFT allows efficient Fourier Transform processing using a VisBuffer which encapsulates a chunk of visibility (typically all baselines for one time) together with all the information needed for processing (e.g. UVW coordinates).

Gridding and degridding in WProjectFT are performed using a novel sort-less algorithm. In this approach, the gridded plane is divided into small patches, a cache of which is maintained in memory using a general-purpose LatticeCache class. As the (time-sorted) visibility data move around slowly in the Fourier plane, patches are swapped in and out as necessary. Thus, optimally, one would keep at least one patch per baseline.

A grid cache is defined on construction. If the gridded uv plane is smaller than this, it is kept entirely in memory and all gridding and degridding is done entirely in memory. Otherwise a cache of tiles is kept an paged in and out as necessary. Optimally the cache should be big enough to hold all polarizations and frequencies for all baselines. The paging rate will then be small. As the cache size is reduced below this critical value, paging increases. The algorithm will work for only one patch but it will be very slow!

This scheme works well for arrays having a moderate number of antennas since the saving in space goes as the ratio of baselines to image size. For the ATCA, VLBA and WSRT, this ratio is quite favorable. For the VLA, one requires images of greater than about 200 pixels on a side to make it worthwhile.

The FFT step is done plane by plane for images having less than 1024 * 1024 pixels on each plane, and line by line otherwise.

The gridding and degridding steps are implemented in Fortran for speed. In gridding, the visibilities are added onto the grid points in the neighborhood using a weighting function. In degridding, the value is derived by a weight summ of the same points, using the same weighting function.

Example

See the example for SkyModel.

Motivation

Define an interface to allow efficient processing of chunks of visibility data

To Do