bmaj = image.header.bmaj #assume no smoothing in convolving function (sdgrd.xtype = -16) newBmaj = bmaj if sdgrd.xtype == -12: convolveMaj = sdgrd.xparm[2]/3600. # convolving function FWHM in degrees #Convolved image resolution adds in quadrature newBmaj = math.sqrt( (bmaj*bmaj) + (convolveMaj*convolveMaj)) print 'Gaussian Convolving function:' print bmaj*3600., convolveMaj*3600., '->',newBmaj*3600. if sdgrd.xtype == -16: #Convolved image resolution adds in quadrature newBmaj = bmaj print 'Sync Bessel Convolving function FWHM :', newBmaj image.header.bmaj = newBmaj image.header.bmin = newBmaj image.update() # This step does not seem to work! # Work-around: write file, update header, read. ## keep track of the latest cube squence for later processing outseq = AIPSCat()[mydisk][-1].seq gridType = image.header.ctype[0] #transfer coordinate back after gridding print 'Data Coordinate type: ', xType, yType xlen = len(xType) if (xlen < 2): xType = xType + '-' if (xlen < 3): xType = xType + '-' if (xlen < 4): xType = xType + '-'
sdgrd.aparm[4] = 14 sdgrd.aparm[5] = 31 sdgrd.aparm[6] = 30 # Imsize sdgrd.imsize[1] = 100 sdgrd.imsize[2] = 110 sdgrd.go() from Wizardry.AIPSData import AIPSImage as WizAIPSImage image = WizAIPSImage(AIPSCat()[mydisk][-1].name, \ AIPSCat()[mydisk][-1].klass, \ mydisk, AIPSCat()[mydisk][-1].seq) image.header.niter = 1 # Allow down stream IMSTATs to sum correctly image.header.update() image.update() ## 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 fittp.inseq=AIPSCat()[mydisk][-1].seq outimage = os.path.splitext(myfiles[0])[0]+'_cube.fits' if os.path.exists(outimage): os.remove(outimage) print 'Removed existing file to make room for new one :',outimage fittp.dataout='PWD:'+outimage fittp.go() print 'Wrote',outimage
bmaj = image.header.bmaj #assume no smoothing in convolving function (sdgrd.xtype = -16) newBmaj = bmaj if sdgrd.xtype == -12: convolveMaj = sdgrd.xparm[2] / 3600. # convolving function FWHM in degrees #Convolved image resolution adds in quadrature newBmaj = math.sqrt((bmaj * bmaj) + (convolveMaj * convolveMaj)) print 'Gaussian Convolving function:' print bmaj * 3600., convolveMaj * 3600., '->', newBmaj * 3600. if sdgrd.xtype == -16: #Convolved image resolution adds in quadrature newBmaj = bmaj print 'Sync Bessel Convolving function FWHM :', newBmaj image.header.bmaj = newBmaj image.header.bmin = newBmaj image.update() # This step does not seem to work! # Work-around: write file, update header, read. ## keep track of the latest cube squence for later processing outseq = AIPSCat()[mydisk][-1].seq gridType = image.header.ctype[0] #transfer coordinate back after gridding print 'Data Coordinate type: ', xType, yType xlen = len(xType) if (xlen < 2): xType = xType + '-' if (xlen < 3): xType = xType + '-' if (xlen < 4): xType = xType + '-'