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

【问题标题】：Extracting third axes, data cube pyfits提取第三个轴，数据立方体 pyfits
【发布时间】：2013-03-16 18:01:55
【问题描述】：

我有一个数据立方体拟合文件，由常规 (512x512) 二维网格的强度组成。然而，还有第三个轴；该轴上的每个值对应一个频率。在每个频率上，都会拍摄一张这样的 2d 图片。

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

顺便用 Python 编程。

【问题讨论】：

【解决方案1】：

通常以与轴“n”的坐标相同的方式给出。参考值“CRVALn”、参考位置“CRPIXn”、每个像素/通道的增量“CDELTn”以及在这些关键字“CTYPEn”中给出的坐标类型。因此，您可以将立方体中每个“通道”的值计算为

(channel_number - CRPIXn) * CRDELTn + CRVALn

既然你说你解决了，也许你可以发布你的解决方案以造福他人？

编辑/添加：

有关这在 FITS 文件中如何表达的更多信息： http://fits.gsfc.nasa.gov/fits_wcs.html

尤其是论文“FITS 中光谱坐标的表示”Greisen 等人。（2005 年）（link 在 NASA ADS 上发表论文）

【讨论】：