我试图使用astropy模块fits从FITS文件中读取数据,然后使用标准的numpy数组处理。但是,由于某些原因,我收到了以下错误:
IndexError: too many indices这是我正在使用的代码:
from astropy.io import fits
import matplotlib.pyplot as plt
hdulist = fits.open('/Users/iMacHome/Downloads/spec-1959-53440-0605.fits')
hdu = hdulist[1]
data = hdu.data
flux = data[:, 1] 对flux = data[:, 1]行的错误跟踪。
loglam = data[:, 2]这也许是天文学家可以回答的一个问题(或者,特别是,天文学家熟悉SDSS中的.fits文件),但我欢迎numpy和python用户的输入!
发布于 2014-08-26 19:57:32
我刚从SDSS服务台得到以下答复:
取代:
flux = data[:,0]
loglam = data[:,1]使用
flux = data['flux']
loglam = data['loglam']这是访问Numpy记录数组中字段的正确方法。
https://stackoverflow.com/questions/25514130
复制相似问题