从python的.pha文件中提取数据

Question

Python 2.7.13

读取spectrum.pha头文件时出现问题

我尝试使用：

from astropy.io import fits

BBspectrum = fits.open('BBspectrum.pha')
hdul = fits.open(BBspectrum)
print hdul[0].header['MJD-OBS']

希望我的代码返回与给定.pha数据集的MJD-OBS头关联的值

谢谢