在CASA中将JY/BEAM转化为K

代码如下

// An highlighted block
inputname='IRAS16293_spw0_image' #input name of fits file without ".fits" 
imagename=inputname + '.fits'
outfile=inputname + '_K'
imhead(imagename=imagename,mode='list')

#please input rest frequency in unit of Hz in image header#
restfreq=input()       

#please input beammajor in unit of arcsec in image header#
bmaj=input()

#please input beamminor in unit of arcsec in image header#
bmin=input()

constant=1000*1.36*(3e10/restfreq)**2/bmaj/bmin

print constant

expr='IM0*' + str(constant)

os.system('rm -rf ' + outfile )

immath(imagename=imagename,mode='evalexpr',expr=expr,outfile=outfile)

imhead(imagename=outfile,mode='put',hdkey='bunit',hdvalue='K')

os.system('rm -rf ' + outfile + 'fits')

exportfits(imagename=outfile, fitsimage=outfile + '.fits')

主要计算依据如下：