bonn_imaging_examples.py

File contents

###"Frequency" definition of image
taskname           = "clean"
vis                =  "ngc5921.demo.src.split.ms.contsub"
imagename          =  "ngc5921_freq"
mode               =  "frequency"
psfmode            =  "clark"
niter              =  6000
gain               =  0.1
threshold          =  8.0
mask               =  ""
nchan              =  20
start              =  "1.412849GHz"
width               =  "52kHz"
interpolation      ='cubic'
imsize             =  [256, 256]
cell               =  ['15.0arcsec', '15.0arcsec']
stokes             =  "I"
field              =  "0"
spw                =  ""
weighting          =  "briggs"
robust             =  0.5
timerange          =  ""
restfreq           =  ""
npixels            =  0
phasecenter        =  ""
async              =  False

###GALACTIC coordinates
taskname           = "clean"
vis                =  "ngc5921.demo.src.split.ms.contsub"
imagename          =  "ngc5921_gal"
mode               =  "velocity"
psfmode            =  "hogbom"
niter              =  6000
gain               =  0.1
threshold          =  8.0
nchan              =  38
start              =  "1595km/s"
width              =  "-5.5km/s"
imsize             =  [256, 256]
cell               =  ['15.0arcsec', '15.0arcsec']
stokes             =  "I"
field              =  "0"
spw                =  ""
weighting          =  "briggs"
robust             =  0.5
phasecenter        =  "GALACTIC 8deg 47.9deg"
async              =  False


#####1950

taskname           = "clean"
vis                =  "ngc5921.demo.src.split.ms.contsub"
imagename          =  "ngc5921_B1950"
mode               =  "velocity"
psfmode            =  "hogbom"
niter              =  6000
gain               =  0.1
threshold          =  8.0
nchan              =  38
start              =  "1595km/s"
width              =  "-5.5km/s"
imsize             =  [256, 256]
cell               =  ['15.0arcsec', '15.0arcsec']
stokes             =  "I"
field              =  "0"
spw                =  ""
weighting          =  "briggs"
robust             =  0.5
noise              =  "0.0Jy"
phasecenter        =  "B1950 15h19m29 5d15m9.73"
async              =  False


#### try clark mosaic

taskname           = "clean"
vis                =  "orion.ms"
imagename          =  "clark_mosaic"
mode               =  "mfs"
psfmode            =  "clark"
imsize             =  [300, 300]
cell               =  ['2arcsec', '2arcsec']
phasecenter        =  "6"
stokes             =  "I"
niter              =  7000
gain               =  0.1
threshold          =  0.0
field              =  "2~10"
spw                =  "0,1"
modelimage         =  ""
weighting          =  "natural"
mosweight          =  False
imagermode         =  "mosaic"
ftmachine          =  "mosaic"
cyclefactor        =  1.5
cyclespeedup       =  -1
scaletype          =  "SAULT"
async              =  False

#################################
##muliscale
##################################
taskname           = "clean"
vis                =  "orion.ms"
imagename          =  "mscale_mosaic"
mode               =  "mfs"
imsize             =  [300, 300]
cell               =  ['2arcsec', '2arcsec']
phasecenter        =  "6"
stokes             =  "I"
niter              =  500
gain               =  0.3
threshold          =  0.0
field              =  "2~10"
spw                =  "0,1"
weighting          =  "natural"
imagermode         = "mosaic"
mosweight          =  False
ftmachine          =  "mosaic"
scaletype          =  "SAULT"
multiscale         = [0,3,10]
async              =  False


############################################
###playing with different scales
###########################################
taskname           = "clean"
vis                =  "orion.ms"
imagename          =  "mscale_mosaic3"
mode               =  "mfs"
alg                =  "multiscale"
imsize             =  [300, 300]
cell               =  ['2arcsec', '2arcsec']
phasecenter        =  "6"
stokes             =  "I"
niter              =  2500
gain               =  0.3
threshold          =  0.0
field              =  "2~10"
spw                =  "0,1"
weighting          =  "natural"
imagermode         =  "mosaic"
ftmachine          =  "mosaic"
multiscale         =  [0, 5, 10, 20, 30]
negcomponent       =  -1






### feather multiscale with 
taskname           = "feather"
imagename          =  "testfeather.im"
highres            =  "mscale_mosaic.image"
lowres             =  "orion.gbt.im"
#feather(imagename="testfeather.im",highres="mscale_mosaic.image",lowres="orion.gbt.im")

###check out total flux of the two images

#################################
##muliscale with SD image as starting model image...then compare
### with no SD and with feathered image
##################################
taskname           = "clean"
vis                =  "orion.ms"
imagename          =  "mscale_mosaic_startsd"
mode               =  "mfs"
imsize             =  [300, 300]
cell               =  ['2arcsec', '2arcsec']
phasecenter        =  "6"
stokes             =  "I"
niter              =  500
gain               =  0.3
threshold          =  0.0
field              =  "2~10"
spw                =  "0,1"
weighting          =  "natural"
imagermode         = "mosaic"
mosweight          =  False
ftmachine          =  "mosaic"
scaletype          =  "SAULT"
multiscale         = [0,3,10]
modelimage         =  "orion.gbt.im"
async              =  False