Pyxis script for imaging prac

pyxis-imaging.py

## imaging Tut. 
## update Pyxis if you haven't done so since 05 Feb 2015.
## Also ensure that you have lwimager installed
# run as $ pyxis msname image[options] 
# options can be: cellsize,npix,gain,wprojplanes,column/data, etc 
# to make a per N channel image add im.IMAGE_CHANNELIZE=N
# Find the MSs at jake:~makhathini/ratt/imaging

import Pyxis
import ms
import im
import im.lwimager
from Pyxis.ModSupport import *

define("MSDIR","MS","Store MSs here")
MS = II("${MSDIR>/}vla_10m30s.MS")
MS_List = [MS] + [II("${MSDIR>/}")+"vla_%dh30s"%d for d in range(2,14,2)]

DESTDIR_Template = "${OURDIR>/}plots${-<MS:BASE}"
OUTFILE_Template = "${DESTDIR>/}${MS:BASE}${-<LABEL}"

im.stokes = 'I'
im.cellsize = '2arcsec'
im.npix = 2048
im.niter = 1000
im.weight = "briggs"
im.robust = 0
im.CLEAN_ALGORITHM = 'csclean'

def image(**kw):
    im.lwimager.make_image(restore_lsm=False,**kw)

def runall(**kw):
    per_ms( lambda : image(**kw) )