本文整理汇总了Python中AIPSTask.AIPSTask.antennas[1:]方法的典型用法代码示例。如果您正苦于以下问题：Python AIPSTask.antennas[1:]方法的具体用法？Python AIPSTask.antennas[1:]怎么用？Python AIPSTask.antennas[1:]使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类AIPSTask.AIPSTask的用法示例。
在下文中一共展示了AIPSTask.antennas[1:]方法的7个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: runvplot
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import antennas[1:] [as 别名]
def runvplot(indata, sources, stokes, timer, anten, basel, bif, eif, bchan, echan, docalib, gainuse, flagver, doband, bpver, aparm, bparm, refant, dotv, nplots):
vplot = AIPSTask('VPLOT')
vplot.indata = indata
vplot.sources[1:] = sources
vplot.stokes = stokes
vplot.timerang[1:] = timer
vplot.antennas[1:] = anten
vplot.baseline[1:] = basel
vplot.bchan = bchan
vplot.echan = echan
vplot.avgchan = 1
vplot.bif = bif
vplot.eif = eif
vplot.docalib = docalib
vplot.gainuse = gainuse
vplot.flagver = flagver
vplot.doband = doband
vplot.bpver = bpver
vplot.aparm[1:] = aparm
vplot.bparm[1:] = bparm
vplot.crowded = 0
vplot.refant = refant
vplot.dotv = dotv
vplot.xinc = 1
vplot.nplots = nplots
vplot.go()示例2: pfring
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import antennas[1:] [as 别名]
def pfring (aipsname,refant,antennas,source,indisk,delaywin=600,ratewin=20,\
solint=1,snr=2,logfiledir='./'):
uvdata = AIPSUVData (aipsname,'FITS',indisk,1)
fring = AIPSTask ('FRING')
fring.refant = refant
fring.indata = uvdata
fring.calsour[1:] = [source]
fring.antennas[1:] = antennas
fring.solint = solint
fring.aparm[1:] = [0,0,0,0,0,2,snr,0,0]
fring.dparm[1:] = [0,delaywin,ratewin,0,0,0,0,0,0]
fring.weightit = 1
stdout = sys.stdout; sys.stdout = open(logfiledir+aipsname+'.log','a')
fring.go()
sys.stdout.close(); sys.stdout = stdout示例3: runuvplt
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import antennas[1:] [as 别名]
def runuvplt(indata, sources, stokes, timer, anten, basel, bif, eif, docalib, gainuse, flagver, doband, bpver, aparm, bparm, doweight, refant, dotv):
uvplt = AIPSTask('UVPLT')
uvplt.indata = indata
uvplt.sources = sources
uvplt.stokes = stokes
uvplt.timerang = timer
uvplt.antennas[1:] = anten
uvplt.baseline[1:] = basel
uvplt.bif = bif
uvplt.eif = eif
uvplt.docalib = docalib
uvplt.gainuse = gainuse
uvplt.flagver = flagver
uvplt.doband = doband
uvplt.bpver = bpver
uvplt.aparm[1:] = aparm[1:]
uvplt.bparm[1:] = bparm[1:]
uvplt.doweight = doweight
uvplt.refant = refant
uvplt.dotv = dotv
uvplt.go()示例4: runuvhgm
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import antennas[1:] [as 别名]
def runuvhgm(indata, sources, anten, basel, SEFDbif, SEFDeif, SEFDbchan, SEFDechan):
uvhgm = AIPSTask('UVHGM')
uvhgm.indata = indata
uvhgm.sources[1:] = sources
uvhgm.antennas[1:] = anten
uvhgm.baseline[1:] = basel
uvhgm.bif = SEFDbif
uvhgm.eif = SEFDeif
uvhgm.flagver = 0
uvhgm.stokes ='HALF'
uvhgm.bchan = SEFDbchan
uvhgm.echan = SEFDechan
uvhgm.doall = 1
uvhgm.axtype ='H'
uvhgm.timerang[1:] = 0, 0, 0, 0, 0, 0, 0, 0
# uvhgm.timerang[1:] = 0, 0, 0, 0, 3, 0, 0, 0
# uvhgm.pixrange[1:] = -20, 20
uvhgm.dotv = -1
uvhgm.docal = 1
uvhgm.gainuse = 0
uvhgm.pixrange[1:] = -10, 10
uvhgm.doband = -1
uvhgm.nboxes = 1000
uvhgm.go()示例5: runpossm
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import antennas[1:] [as 别名]
def runpossm(indata, sources, timer, anten, basel, aparm, bparm, bif, eif, docalib, gainuse, flagver, stokes, doband, bpver, codetype, solint, nplots, dotv, freqid):
possm = AIPSTask('POSSM')
possm.indata = indata
possm.sources[1:] = sources
possm.timer[1:] = timer
possm.aparm[1:] = aparm[1:]
possm.aparm[9] = 1
possm.bparm[1:] = bparm[1:]
possm.antennas[1:] = anten
possm.baseline[1:] = basel
possm.bif = bif
possm.eif = eif
possm.docalib = docalib
possm.gainuse = gainuse
possm.flagver = flagver
possm.stokes = stokes
possm.doband = doband
possm.bpver = bpver
possm.codetype = codetype
possm.solint = solint
possm.nplots = nplots
possm.dotv = dotv
possm.freqid = 1
possm.go()示例6: __init__
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import antennas[1:] [as 别名]
#.........这里部分代码省略.........
if flag.exists():
flag.clrstat()
flag.zap()
fitldFL = AIPSTask("FITLD")
fitldFL.datain = args["flagPath"]
fitldFL.outname = self.args["name"]
fitldFL.outseq = 100
fitldFL.ncount = 1
fitldFL.go()
print("Copying flagging")
tacop = AIPSTask("TACOP")
tacop.indata = flag
tacop.invers = 3
tacop.ncount = 1
tacop.outname = self.args["name"]
tacop.outseq = self.args["inseq"]
tacop.outclass = 'UVDATA'
tacop.outdisk = 1
tacop.outver = 0
tacop.inext = 'FG'
tacop.go()
if self.args["doBP"]:
uvdata.zap_table('BP', -1) #create a new bp table for each reftelly
print("Running bandpass")
bpass = AIPSTask('BPASS')
bpass.indata = uvdata
bpass.calsour[1] = self.args["bandPassCal"]
bpass.timer = self.args["time"]
bpass.refant = self.args["refTelly"]
if self.args["excludeTelly"]:
bpass.antennas[1:] = self.args["excludedTellys"]
bpass.go()
"""
uvdata.zap_table('TY', -1)
uvdata.zap_table('GC', -1)
print("Running antab")
antab = AIPSTask('ANTAB')
antab.indata = uvdata
antab.calin = self.args["antPath"]
antab.go()
print("Running apcal")
apcal = AIPSTask('APCAL')
apcal.indata = uvdata
apcal.timer = self.args["time"]
apcal.tyver = 1
apcal.gcver = 1
apcal.go()
snVers = snVers + 1
print("apcal made snVers {0}".format(snVers))
#applies new SN table to CL
print("Running clcal.")
clcal = AIPSTask('clcal')
clcal.indata = uvdata
clcal.calsour[1] = self.args["cal"]
clcal.timerang = self.args["time"]
clcal.snver = snVers
clcal.inver = snVers
clcal.gainver = clInit #apply to original cl
clcal.timer = self.args["time"]
clcal.interpol = "ambg"示例7: __init__
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import antennas[1:] [as 别名]
def __init__(self, clVers, **args):
self.clVers = clVers
self.args = args
self.clVers = 3
AIPS.userno=self.args["user"]
uvdata = AIPSUVData(self.args["name"], 'UVDATA', 1, self.args["inseq"])
t = [] #files t with [day, hour, min, sec, day, hour, min, sec]
with open(self.args["timeList"], 'r') as f:
for l in f.readlines():
t.append(map(str,l.split()[1:9]))
for j in t:
realtime = ''.join([j[i] for i in range(len(j))])[0:7]
imgname = realtime[0:5]
#need to shorten the names for AIPS
imagr = AIPSTask('IMAGR')
imagr.indata = uvdata
imagr.sources[1] = self.args["source"]
imagr.docalib = 1
imagr.gainuse = self.clVers
imagr.bchan = self.args["SCbchan"]
imagr.echan = self.args["SCechan"]
imagr.nchav = (self.args["SCechan"] - self.args["SCbchan"] + 1)
#averages over all channels
if self.args["doBP"]:
imagr.doband = 1
imagr.bpver = 1
if self.args["excludeTelly"]:
imagr.antennas[1:] = self.args["excludedTellys"]
imagr.outname = imgname
imagr.outseq = 1
imagr.cellsize = AIPSList([0.0001,0.0001])
#imagr.imsize = AIPSList([1024,1024])
imagr.imsize = AIPSList([256,256])
imagr.nboxes = 1
imagr.clbox[1] = self.args["fitBox"]
imagr.niter = 1000
imagr.timerang[1:] = [int(j[i]) for i in range(len(j))]#AIPSList(j)
good = True
try:
imagr.go()
except RuntimeError as e:
print("\nSomething went wrong!\n => {0}".format(e))
with open(os.getcwd() + '/images_{0}/badtimes.txt'.format(self.args["date"]), 'a') as out:
out.write("{0}\n".format(imgname))
good = False
if good:
imageClean = AIPSImage(imgname, 'ICL001',1,1)
imageDirty = AIPSImage(imgname, 'IBM001',1,1)
imageClean.clrstat() #makes sure AIPS does not trip
imageDirty.clrstat()
jmfit = AIPSTask('JMFIT')
jmfit.indata = imageClean
jmfit.blc[1] = .80*self.args["fitBox"][1]#fraction to increase box size
jmfit.blc[2] = .80*self.args["fitBox"][2]
jmfit.trc[1] = 1.20*self.args["fitBox"][3]
jmfit.trc[2] = 1.20*self.args["fitBox"][4]
jmfit.niter = 1000
jmfit.doprint = 1 #CHANGED TO doprint!
address = os.getcwd() + '/images_{0}/'.format(self.args["date"]) + realtime
jmfit.fitout = address + '.crd'
jmfit.go()
RA = [None]*4 #holds locations
DEC = [None]*4
with open(address + '.crd', 'r') as f:
for l in f.readlines():
temp = l.split()
#for j in range(len(temp)):
if len(temp)>0:
if temp[0] == 'RA':
RA[0] = int(temp[1])
RA[1] = int(temp[2])
RA[2] = float(temp[3])
RA[3] = float(temp[5])
elif temp[0] =='DEC':
DEC[0] = int(temp[1])
DEC[1] = int(temp[2])
DEC[2] = float(temp[3])
DEC[3] = float(temp[5])
with open(os.getcwd() + '/images_{0}/locations.txt'.format(self.args["date"]), 'a') as out:
out.write("{0}\t{1}\t{2}\t{3:^10}\t{4:^6}\t{5}\t{6}\t{7:^7}\t{8}\n".
format(realtime, RA[0], RA[1], RA[2], RA[3], \
DEC[0], DEC[1], DEC[2], DEC[3]))
#makes contour plot
kntr = AIPSTask('KNTR')
kntr.indata = imageClean
kntr.levs = AIPSList([2,3,4,5,7,10,13,17])
kntr.dogrey = -1
kntr.dotv = -1
kntr.dovect = -1
#kntr.blc[1] .80*self.args["fitBox"][1]
#kntr.blc[2] .80*self.args["fitBox"][2]
#kntr.trc[1] 1.20*self.args["fitBox"][3]
#kntr.trc[2] 1.20*self.args["fitBox"][4]
#.........这里部分代码省略.........