在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')
主要计算依据如下: