ডিপথট:Flagging tiles
ডিপথট থেকে
import os, elementflag
o=os.system
ef = elementflag.station()
def visibilise(n):
global o, ef
# create MS
pref = '3C196.tilex.'+str(n)
ms = 'ms/'+pref+'.MS'
o('rm -rf '+ms)
print '... creating '+ms
o('cp -r ms/3C196.root.MS '+ms)
# predict
print '... making tile configuration default'
ef.ElementFlagDefault('ms/3C196.root.MS', ms)
print '... predicting in DATA of '+ms
o('calibrate-stand-alone -v -f '+ms+' pset/sim.parset sky/3C196.sky &> log/'+pref+'.sim')
print '... turning '+str(n)+' tiles of all the stations off in '+ms
ef.WriteElementFlag(ms, 'all', '', n)
print '... predicting with turned off tiles in MODEL_DATA of '+ms
o('calibrate-stand-alone -v -f '+ms+' pset/sim.MDATA.parset sky/3C196.sky &> log/'+pref+'.sim')
# calibrate
print '... calibrating and saving in CORRECTED_DATA of '+ms
o('calibrate-stand-alone -v -f '+ms+' pset/cal.parset sky/3C196.sky &> log/'+pref+'.cal')
# imaging with lwimager
comm = "lwimager ms="+ms+" image=dimg/"+pref+".DATA.img fits=fits/"+pref+".DATA.fits npix=900 cellsize=40arcsec data=DATA uvrange=1~10km stokes=IQUV niter=0 mode=mfs gain=0.1 threshold=0.0mJy operation=csclean cyclefactor=1.5 cyclespeedup=-1"
print "... executing %s"%comm
o(comm)
comm = "lwimager ms="+ms+" image=dimg/"+pref+".MODEL_DATA.img fits=fits/"+pref+".MOEL_DATA.fits npix=900 cellsize=40arcsec data=MODEL_DATA uvrange=1~10km stokes=IQUV niter=0 mode=mfs gain=0.1 threshold=0.0mJy operation=csclean cyclefactor=1.5 cyclespeedup=-1"
print "... executing %s"%comm
o(comm)
comm = "lwimager ms="+ms+" image=dimg/"+pref+".CORRECTED_DATA.img fits=fits/"+pref+".CORRECTED_DATA.fits npix=900 cellsize=40arcsec data=CORRECTED_DATA uvrange=1~10km stokes=IQUV niter=0 mode=mfs gain=0.1 threshold=0.0mJy operation=csclean cyclefactor=1.5 cyclespeedup=-1"
print "... executing %s"%comm
o(comm)
if __name__=='__main__':
ns = [1,3,8,12,16,20]
for i in ns:
print 'predicting, turning '+str(i)+' tiles off, again predicting and finally calibrating'
visibilise(i)