1

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

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

顺便说一句,在 Python 中编程。

4

1 回答 1

1

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

(channel_number - CRPIXn) * CRDELTn + CRVALn

既然您说您解决了它,也许您可​​以发布您的解决方案以造福他人?


编辑/添加:

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

特别是论文“FITS 中光谱坐标的表示”Greisen 等人。(2005 年)(NASA ADS 论文链接)

于 2013-03-25T13:47:25.240 回答