alanloh / nenupy-tv Goto Github PK

Save figure as a fits file with the wcs.

Add a UVW method, taking in input only the Crosslet object.

Ra Dec / alt az pointings
Mini-array list as a bintablehdu
...

Diagnostic plot in png, to be computed each time a nenufartv dataset is ready (every 5min when the telescope is not observing int xst).

Gather the pointing infos and overplot the primary beam circle (deformed with spherical projection?) with respect to the mean observed frequency.
Would be particularly usefull while making movies of the TV to show sources disappear when they go out of the primary beam.

Briggs :

elif weighting == 'briggs':
            ## Robust weighting - we need to calculate it
            dens = numpy.abs(self.bm[0])
            size = dens.shape[0]
            from scipy.ndimage import uniform_filter
            dens = uniform_filter(dens, size=size*localFraction)
            dens /= dens.max()
            dens[numpy.where( dens < 1e-8 )] = 0
            f2 = (5*10**-robust)**2 / (dens**2).mean()
            dens = 1.0 / (1.0 + f2/dens)
            data = data/dens*dens.max()
            data[numpy.where(dens == 0)] = 0.0

Add a mini-array beam calculation in the response module

In /data/nenufar/nenuFAR-TV/log/nenufartv.log

Add a PSF class within the response module

Add a class from which NenuFAR or MiniArray inherit in order to generalize radio array handling.

Expand a little bit more the MiniArray class to take into account rotation, name, abd maybe real ground position.