Python 提取第三轴，数据立方体pyfits

我有一个数据立方体fitsfile，由规则（512x512）2d网格的强度组成。然而，还有第三个轴；该轴上的每个值对应一个频率。在每个频率处，拍摄这样的2d图片

我一直在尝试在fits文件中查找这些频率值，但到目前为止我没有运气。通过查看fits文件的标题，我知道有140个不同的频道，我知道单位等。问题在于找到这140个通道的值

---

顺便说一下，用PYTHON编程。

通常以与轴“n”的坐标相同的方式给出。参考值“CRVALn”、参考位置“CRPIXn”、每个像素/通道的增量“CDELTn”以及这些关键字“CTYPEn”中给出的坐标类型。因此，您可以计算出多维数据集中每个“通道”的值，如下所示：

(channel_number - CRPIXn) * CRDELTn + CRVALn

既然你说你已经解决了这个问题，也许你可以发布你的解决方案，让其他人受益

---

编辑/添加：

有关如何在FITS文件中表达的更多信息：

---

特别是“FITS中光谱坐标的表示”一文（Greisen et al.（2005）（发表在NASA ADS上）

通常以与轴“n”坐标相同的方式给出。参考值“CRVALn”、参考位置“CRPIXn”、每个像素/通道的增量“CDELTn”以及这些关键字“CTYPEn”中给出的坐标类型。因此，您可以计算出多维数据集中每个“通道”的值，如下所示：

(channel_number - CRPIXn) * CRDELTn + CRVALn

既然你说你已经解决了这个问题，也许你可以发布你的解决方案，让其他人受益

---

编辑/添加：

有关如何在FITS文件中表达的更多信息：

---

尤其是格雷森等人（2005年）发表的论文《光谱坐标的拟合表示》（发表在美国宇航局的广告上）

我找到了！问题解决了，很抱歉给您带来不便。我找到了！问题已解决，不便之处，敬请原谅。