我在用 Iraf 生成的拟合文件中有一个光谱。波长轴在标头中编码为：
WAT0_001= 'system=multispec'
WAT1_001= 'wtype=multispec label=Wavelength units=angstroms'
WAT2_001= 'wtype=multispec spec1 = "1 1 2 1. 2.1919422441886 4200 0. 452.53 471'
WAT3_001= 'wtype=linear'
WAT2_002= '.60 1. 0. 3 3 1. 4200.00000000001 1313.88904209266 1365.65012876239 '
WAT2_003= '1422.67911152069 1479.0560707956 1535.24082980747 1584.94609332243'
有没有一种简单的方法可以将其加载到 python 中？
我一直在使用this https://github.com/kgullikson88/General/blob/master/readmultispec.py代码，是 Rick White 给我的。但是，那斯佩蒂尔斯 http://specutils.readthedocs.org/en/latest/specutils/read_fits.html包可能是更好的方法：
from specutils.io import read_fits
spectra_list = read_fits.read_fits_spectrum1d('mymultispec.fits')
print spectra_list[0]
Spectrum1D([ 338.06109619, 395.59234619, 326.0012207 , ...,
660.0098877 , 686.54498291, 689.58374023])
print spectra_list[1].dispersion
<Quantity [ 8293.44875263, 8293.40459999, 8293.36044556,...,
8166.53073537, 8166.48250242, 8166.43426803] Angstrom>
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
