![](https://img-blog.csdnimg.cn/20201014180756928.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
astropy
文章平均质量分 56
nondefined
记录Latex及天文数据处理软件的一些问题
展开
-
备忘--myXCLASSFit.py
备忘录把拟合结果移动到当前文件夹##remove the old runing dataimport osimport shutilpath = '/home/zhang/software/XCLASS-Interface/run/myXCLASSFit/'files = os.listdir(path)if len(files)>0: shutil.rmtree(path)MolfitsFileName = "mol.molfit"experimentalDa原创 2020-12-11 16:38:15 · 115 阅读 · 0 评论 -
获取aplpy可直接读取画图的2Dhdu
由于ALMA数据是4维的,不能直接用aplpy去读。把下面的脚本放在数据所在目录下,调用一下就能用aplpy画ALMA的图了#主要用于连续谱图,不要用输出的hdu覆盖原本的fitsfile,会丢失header里的重要信息from astropy.io import fitsfrom astropy.wcs import WCSdef Get2DHDU(fitsfile): fits_file = fits.open(fitsfile) data = fits_file[0].data hd原创 2020-11-25 20:19:12 · 296 阅读 · 0 评论 -
从ALMA的data cube导出moment0
##### make moment 0 mapfrom spectral_cube import SpectralCube as scfrom astropy.io import fitsimport astropy.units as uimport astropy.constants as cimport matplotlib.pyplot as pltplt.ion()import osspw = ['spw35','spw37']def change_xaxis_in_velo.原创 2020-11-23 19:40:36 · 242 阅读 · 0 评论 -
画rotational-diagram
import numpy as npimport astropy as aimport astropy.units as uimport astropy.constants as cfrom scipy import optimizeimport matplotlib.pyplot as pltplt.ion()B = 24690.2h = 6.626e-27 #erg/sk = 1.38e-16 #erg/KTbg = 2.73 #Knu = [35原创 2020-11-15 15:32:23 · 151 阅读 · 0 评论 -
用spectral_cube抽取cube文件上某个点或者面源的光谱
from astropy.io import fitsfrom spectral_cube import SpectralCube as specimport numpy as npimport matplotlib.pyplot as pltimport osimport sysfrom astropy.coordinates import SkyCoordfrom astropy.wcs import WCSfrom astropy.table import Tableimport .原创 2020-07-17 12:58:36 · 642 阅读 · 0 评论 -
Herschel-sed-fitter代码
翻到去年写的sed点到点拟合的代码。乱七八糟还没写完,还有很多写得不够智能。建议直接用王科老师的higal_sedfitter.######! /usr/bin python3#This progress is wrote by Lixia Yuan#Version v2.0 modified by Chao Zhang at Dec. 26, 2019#This script requires dust_emissivity which can be installed by "pip in.原创 2020-06-22 10:31:33 · 816 阅读 · 4 评论