本文整理汇总了Python中AIPSTask.AIPSTask.optype方法的典型用法代码示例。如果您正苦于以下问题：Python AIPSTask.optype方法的具体用法？Python AIPSTask.optype怎么用？Python AIPSTask.optype使用的例子？那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类AIPSTask.AIPSTask的用法示例。
在下文中一共展示了AIPSTask.optype方法的6个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: runsetjy
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import optype [as 别名]
def runsetjy(indata, sources, bif, eif, zerosp, optype):
setjy = AIPSTask('SETJY')
setjy.indata = indata
setjy.sources[1:] = sources
setjy.bif = bif
setjy.eif = eif
setjy.zerosp[1:] = zerosp, 0, 0, 0
setjy.optype = optype
setjy.go()示例2: runsnplt
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import optype [as 别名]
def runsnplt(indata, inext, invers, sources, nplots, optype, dotv):
snplt = AIPSTask('SNPLT')
snplt.indata = indata
snplt.inext= inext
snplt.invers = invers
snplt.sources[1:] = sources
snplt.nplots = nplots
snplt.optype = optype
snplt.dotv = dotv
snplt.go()示例3: round
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import optype [as 别名]
avspc.ichansel[3][1] = nAverage
avspc.avoption='SUBS'
avspc.go()
#now have fewer channels, with broader frequencies
nChan = round(nChan/nAverage)
dNu = nAverage*dNu
refChan = refChan/nAverage
# Now make an image using the last entry in the catalog
sdgrd.indisk=mydisk
sdgrd.outdisk=mydisk
sdgrd.inname=AIPSCat()[mydisk][-1].name
sdgrd.inclass=AIPSCat()[mydisk][-1].klass
sdgrd.inseq=AIPSCat()[mydisk][-1].seq
sdgrd.optype='-GLS'
sdgrd.xtype=-12
sdgrd.ytype=-12
sdgrd.reweight[1] = 0
sdgrd.reweight[2] = 0.025
# must break up RA into hours minutes seconds
sdgrd.aparm[1]=math.floor(raDeg/15.)
sdgrd.aparm[2]=math.floor(((raDeg/15.)-sdgrd.aparm[1])*60.)
sdgrd.aparm[3]=round(((((raDeg/15.)-sdgrd.aparm[1])*60.)-sdgrd.aparm[2])*60.)
#now break up degrees, but must preserve sign
decSign = 1.
if decDeg < 0.:
decSign = -1.
decDeg = -1. * decDeg
sdgrd.aparm[4]=math.floor(decDeg)示例4: AIPSTask
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import optype [as 别名]
multi = AIPSTask("MULTI")
multi.indata = uvdata
multi.outdisk = indisk
multi.outname = "POINT"
multi.outclass = "UVDATA"
multi.go()
uvdata = AIPSUVData("POINT", "UVDATA", indisk, 1)
tabed = AIPSTask("TABED")
tabed.indata = uvdata
tabed.indisk = indisk
tabed.outdisk = indisk
tabed.inext = "SU"
tabed.optype = "REPL"
tabed.aparm[1:] = 2, 0, 0, 3, 0
tabed.keystrng = "POINT"
tabed.go()
if doscal == 1:
uvdata = AIPSUVData("POINT", "UVDATA", indisk, 1)
calib = AIPSTask("CALIB")
clcal = AIPSTask("CLCAL")
snplt = AIPSTask("SNPLT")
for i in range(1, itercal + 1):
if i > 1:
nchan = uvdata.header.naxis[2]
imagr = AIPSTask("IMAGR")
imagr.nchav = nchan # use imagr to get a clean model!示例5: make_cube
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import optype [as 别名]
def make_cube(args):
print_header("Making image cube")
average_channels(args.average)
sdgrd = AIPSTask("sdgrd")
# Now make an image using the last entry in the catalog
sdgrd.indisk = DISK_ID
sdgrd.outdisk = DISK_ID
sdgrd.baddisk[1] = BADDISK
last = cat.last_entry()
sdgrd.inname = last.name
sdgrd.inclass = last.klass
sdgrd.inseq = last.seq
sdgrd.optype = "-GLS"
sdgrd.reweight[1] = 0
spectra = cat.get_uv(last)
if "center" in args and args.center:
raDeg, decDeg = args.center
else:
raDeg, decDeg = spectra.header.crval[3], spectra.header.crval[4]
# must break up RA into hours minutes seconds
hh, mm, ss = ra_deg2hms(raDeg)
sdgrd.aparm[1] = hh
sdgrd.aparm[2] = mm
sdgrd.aparm[3] = ss
# now break up degrees, but must preserve sign
decSign = 1
if decDeg < 0:
decSign = -1
decDeg = -1 * decDeg
hh, mm, ss = dec_deg2hms(decDeg)
sdgrd.aparm[4] = hh
sdgrd.aparm[5] = mm
sdgrd.aparm[6] = ss
# deal with degrees and/or minutes == 0
if decSign < 0.0:
sdgrd.aparm[4] = -1 * sdgrd.aparm[4]
if sdgrd.aparm[4] == 0:
sdgrd.aparm[5] = -1 * sdgrd.aparm[5]
if sdgrd.aparm[5] == 0:
sdgrd.aparm[6] = -1 * sdgrd.aparm[6]
print raDeg, decDeg, "->", sdgrd.aparm[1:7]
# transfer cellsize
cellsize = round(spectra.header.cdelt[4] * 3600.0)
sdgrd.cellsize[1] = sdgrd.cellsize[2] = cellsize
# sdgrd.xtype=-16 # sync/bessel convolving type
sdgrd.xtype = -12 # gaussian convolving type
# sync/bessel function parameters
if sdgrd.xtype == -16:
sdgrd.xparm[1] = 3.0 * cellsize
sdgrd.xparm[2] = 2.5 * cellsize
sdgrd.xparm[3] = 1.5 * cellsize
sdgrd.xparm[4] = 2
sdgrd.reweight[2] = 0.01
# gaussian parameters
if sdgrd.xtype == -12:
sdgrd.xparm[1] = 5.0 * cellsize
sdgrd.xparm[2] = 1.5 * cellsize # Parameter sets Gaussian FWHM
sdgrd.xparm[3] = 2
sdgrd.xparm[4] = 0
sdgrd.reweight[2] = -1.0e-6
# always make a circuluar convolving function
sdgrd.ytype = sdgrd.xtype
if "size" in args and args.size:
imxSize, imySize = args.size
else:
imxSize = (2 * round(spectra.header.crpix[3] / 1.95)) + 20
imySize = (2 * round(spectra.header.crpix[4] / 1.95)) + 20
print "Ra, Dec : {0}, {1}".format(raDeg, decDeg)
print "Image size (X,Y) : {0}, {1}".format(imxSize, imySize)
print "Cell size : {0}".format(cellsize)
sdgrd.imsize[1] = imxSize
sdgrd.imsize[2] = imySize
sdgrd.go()
seqno, restFreqHz = update_header(args)
outcube = write_image_cube(restFreqHz, args.uniqueid)
return seqno, restFreqHz, outcube示例6: AIPSTask
# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import optype [as 别名]
multi = AIPSTask('MULTI')
multi.indata = uvdata
multi.outdisk = 2
multi.outname = 'POINT'
multi.outclass = 'UVDATA'
multi.go()
uvdata = AIPSUVData('POINT','UVDATA',2,1)
tabed = AIPSTask('TABED')
tabed.indata = uvdata
tabed.indisk=2
tabed.outdisk=2
tabed.inext = 'SU'
tabed.optype = 'REPL'
tabed.aparm[1:] = 2, 0, 0, 3, 0
tabed.keystrng = 'POINT'
tabed.go()
if doscal == 1:
uvdata = AIPSUVData('POINT','UVDATA',2,1)
calib = AIPSTask('CALIB')
clcal = AIPSTask('CLCAL')
snplt = AIPSTask('SNPLT')
for i in range(1,itercal+1):
if i >1:
nchan = uvdata.header.naxis[2]
imagr = AIPSTask('IMAGR')
imagr.nchav = nchan #use imagr to get a clean model!