Python AIPSTask.outclass方法代码示例

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]
def runsplat(uvdata, outchan, tint, sbandl, sbandu, smootha, smoothb, smoothc, indisk):
	splat = AIPSTask('SPLAT')
	splat.indata = uvdata
	splat.outname = uvdata.name
	splat.outdisk = indisk
	splat.outclass = 'SPLAT'
	splat.outseq = uvdata.seq
	splat.solint = 0
	if outchan != -1 :
		aparm1 = 3
		splat.channel = uvdata.header['naxis'][2] / outchan
		splat.chinc = uvdata.header['naxis'][2] / outchan
	if tint != -1 :
		splat.solint = tint / 60
	splat.bif = sbandl # SBLANL -- Lower sub-band for extraction
	splat.eif = sbandu # SBANDU - Upper sub-band for extraction
	splat.douvcomp = -1
	splat.aparm[1] = 3
	splat.aparm[2] = 1
	# Smoothing options!?!?
	splat.smooth[1] = smootha # 1 = hanning default
	splat.smooth[2] = smoothb # 4 = hanning default
	splat.smooth[3] = smoothc # 1 = hanning default
	splat.flagver = 0
	splat.go()

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]
def psplit (aipsname,source,indisk,logfiledir='./'):
    uvdata = AIPSUVData (aipsname,'FITS',indisk,1)
    split = AIPSTask ('split')
    split.indata = uvdata
    split.outclass = 'SPLIT'
    split.docalib = 1
    stdout = sys.stdout; sys.stdout = open(logfiledir+aipsname+'.log','a')
    split.go()
    sys.stdout.close(); sys.stdout = stdout
    uvdata = AIPSUVData(source,'SPLIT',indisk,1)
    uvdata.rename(aipsname,'SPLIT',1)

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]
def average_channels(binsize):
    if binsize > 1:
        print "Averaging ", binsize, " Spectral Channels"
    else:
        print "Not Averaging Spectral Channels"
        return

    avspc = AIPSTask("avspc")

    last = cat.last_entry()
    spectra = cat.get_uv(last)
    nChan = round(spectra.header.naxis[2])
    print spectra.header.naxis, nChan

    # now average channels to reduce the image plane data volumn
    avspc.indisk = avspc.outdisk = DISK_ID
    avspc.outclass = ""
    avspc.inname = last.name
    avspc.inclass = last.klass
    avspc.inseq = last.seq
    avspc.channel = binsize
    avspc.ichansel[1][1] = 1
    avspc.ichansel[2][1] = nChan
    avspc.ichansel[3][1] = 1
    avspc.avoption = "SUBS"
    avspc.outcl = "avg"
    avspc.go()

    # now have fewer channels, with broader frequencies
    spectra.header.naxis[2] = round(nChan / binsize)  # write back to header
    print "naxis", spectra.header.naxis
    print "cdelt", spectra.header.cdelt
    dNu = spectra.header.cdelt[2]
    dNu = binsize * dNu  # write back to header !!!
    print "cdelt", spectra.header.cdelt, "dNu", dNu, "binsize", binsize

    refChan = spectra.header.crpix[2]
    print "refChan", refChan, "crpix", spectra.header.crpix
    refChan = refChan / binsize  # write back to header !!!
    print "refChan", refChan, "crpix", spectra.header.crpix
    return

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]
    bChan = 1
if bChan > nChan:
    bChan = 1
if eChan > nChan:
    eChan = nChan
if eChan < 1:
    eChan = nChan
#Run momnt task
momnt.indisk=mydisk
momnt.outdisk=mydisk
momnt.inname=AIPSCat()[mydisk][-1].name
momnt.inclass='IMLIN'
momnt.inseq=1
momnt.icut=-10000.
momnt.flux=-.00001
momnt.outclass='0'
momnt.cellsize[1] = 0
momnt.cellsize[2] = 0
momnt.blc[1]=bChan
momnt.blc[2]=0
momnt.blc[3]=0
momnt.trc[1]=eChan
momnt.trc[2]=0
momnt.trc[3]=0
momnt.go()

# prepare to zap after the copy
image = AIPSImage( AIPSCat()[mydisk][-1].name, \
                   AIPSCat()[mydisk][-1].klass, mydisk, \
                   AIPSCat()[mydisk][-1].seq)

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]

#.........这里部分代码省略.........
                    snVers = i[0]
                if i[1] == 'AIPS CL' and i[0] > clVers:
                    clVers = i[0]
            print("Deleting tables.")
            uvdata.zap_table('SN', -1)
            for k in range(clVers,1,-1):
                uvdata.zap_table('CL', k)
            snVers = 0
            clVers = 1
        if imageClean.exists():
            imageClean.clrstat()
            imageClean.zap()
        if imageDirty.exists():
            imageDirty.clrstat()
            imageDirty.zap()
        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()

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]
#print spectra.header

# set number of channels to average
if inNAve > 0:
    nAverage = inNAve
else:
    nAverage = 3

if nAverage > 1:
    print 'Averaging ',nAverage,' Spectral Channels'
else:
    print 'Not Averaging Spectral Channels'
# now average channels to reduce the image plane data volumn
avspc.indisk=mydisk
avspc.outdisk=mydisk
avspc.outclass=''
avspc.inname=AIPSCat()[mydisk][0].name
avspc.inclass=AIPSCat()[mydisk][0].klass
avspc.inseq=AIPSCat()[mydisk][0].seq
avspc.channel=nAverage
avspc.ichansel[1][1] = 1
avspc.ichansel[2][1] = nChan
avspc.ichansel[3][1] = 1
avspc.avoption='SUBS'
avspc.go()

#now have fewer channels, with broader frequencies
nChan = round(nChan/nAverage)
dNu = nAverage*dNu
refChan = refChan/nAverage

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]
AIPS.userno = 1999

print AIPSCat(1)

image = AIPSImage('MANDELBROT', 'MANDL', 1, 1)
if image.exists():
    image.zap()
    pass

try:
    mandl = AIPSTask('mandl')
    mandl.outdata = image
    mandl.imsize[1:] = [ 512, 512 ]
    mandl.go()

    mandl.outclass = 'TEMP'
    mandl.go()

    print AIPSCat(1)

    try:
        AIPSCat(1).zap(klaas='TEMP')
    except TypeError:
        pass
    else:
        raise AssertionError

    AIPSCat(1).zap(klass='TEMP')

    print AIPSCat(1)

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]
    uvsrt.outdisk = indisk
    uvsrt.go()

    uvdata.zap_table("CL", 1)
    indxr = AIPSTask("INDXR")
    indxr.indata = uvdata
    indxr.cparm[1] = 360
    indxr.cparm[2] = 360
    indxr.cparm[3] = 0.25
    indxr.go()

    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()

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]
eChan = lineChan + lineWidthChan
if bChan < 1:
    bChan = 1
if bChan > nChan:
    bChan = 1
if eChan > nChan:
    eChan = nChan
if eChan < 1:
    eChan = nChan
momnt.blc[1]=bChan
momnt.blc[2]=0
momnt.blc[3]=0
momnt.trc[1]=eChan
momnt.trc[2]=0
momnt.trc[3]=0
momnt.outclass='0'
momnt.go()

## Use the subim task to rename the output image
subim.indisk=mydisk
subim.inname=AIPSCat()[mydisk][-1].name
subim.inclass=AIPSCat()[mydisk][-1].klass
subim.inseq=AIPSCat()[mydisk][-1].seq
subim.outclass='22'
subim.outdi=mydisk
subim.go()

## and write the last thing now in the catalog to disk
fittp.indisk=mydisk
fittp.inname=AIPSCat()[mydisk][-1].name
fittp.inclass=AIPSCat()[mydisk][-1].klass

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]
	uvsrt.outdisk = 2
	uvsrt.go()

	uvdata.zap_table('CL',1)
	indxr = AIPSTask('INDXR')
	indxr.indata = uvdata
	indxr.cparm[1]=360
	indxr.cparm[2] = 360
	indxr.cparm[3] = 0.25
	indxr.go()

	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()

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]
from AIPS import AIPS
from AIPSTask import AIPSTask, AIPSList
from AIPSData import AIPSUVData, AIPSImage

time = AIPSList([1,8,18,30, 1,18,30,0])
AIPS.userno = 1010
uvdata = AIPSUVData('gr035', 'UVDATA', 1, 1)
if uvdata.exists():
    uvdata.zap()
fitld = AIPSTask('FITLD')
fitld.datain = 'PWD:gr035.sb5.idi'
fitld.outname = 'gr035'
fitld.outclass = 'UVDATA'
fitld.ncount = 6
fitld.doconcat = 1
#fitld.timer = time
fitld.go()

fring = AIPSTask('FRING')
fring.indata = uvdata
fring.calsour = AIPSList(['J1232-0224'])
fring.docalib = 1
fring.gainuse = 1
fring.timer = time
fring.go()

clcal = AIPSTask('clcal')
clcal.calsour = AIPSList(['J1232-0224'])
clcal.indata = uvdata
clcal.timer = time
clcal.go()

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]
    os.remove(outimage)
    print 'Removed existing file to make room for new one :',outimage

fittp.dataout='PWD:'+outimage
fittp.go()

#Run momnt task
#previously selected channels with the NH3 1-1 line
momnt.indisk=mydisk
momnt.outdisk=mydisk
momnt.inname=AIPSCat()[mydisk][-1].name
momnt.inclass='IMLIN'
momnt.inseq=1
momnt.icut=-300.
momnt.flux=-.005
momnt.outclass='0'
momnt.blc[1]=1
momnt.blc[2]=0
momnt.blc[3]=0
momnt.trc[1]=nChan
momnt.trc[2]=0
momnt.trc[3]=0
momnt.cellsize[1] = 0
momnt.cellsize[2] = 0
print momnt.blc,momnt.trc
momnt.go()

## and write the last thing now in the catalog to disk
fittp.indisk=mydisk
fittp.inname=AIPSCat()[mydisk][-1].name
fittp.inclass=AIPSCat()[mydisk][-1].klass

# 需要导入模块: from AIPSTask import AIPSTask [as 别名]
# 或者: from AIPSTask.AIPSTask import outclass [as 别名]
    bChan = 1
if bChan > nChan:
    bChan = 1
if eChan > nChan:
    eChan = nChan
if eChan < 1:
    eChan = nChan
#Run momnt task
momnt.indisk=mydisk
momnt.outdisk=mydisk
momnt.inname=AIPSCat()[mydisk][-1].name
momnt.inclass='IMLIN'
momnt.inseq=1
momnt.icut=-10000.
momnt.flux=fluxMin
momnt.outclass='0'
momnt.cellsize[1] = 0
momnt.cellsize[2] = 0
momnt.blc[1]=bChan
momnt.blc[2]=0
momnt.blc[3]=0
momnt.trc[1]=eChan
momnt.trc[2]=0
momnt.trc[3]=0
momnt.go()

# prepare to zap after the copy
image = AIPSImage( AIPSCat()[mydisk][-1].name, \
                   AIPSCat()[mydisk][-1].klass, mydisk, \
                   AIPSCat()[mydisk][-1].seq)