我有一个数据立方体 fitfile,由常规 (512x512) 2d 网格的强度组成。然而,还有第三个轴;该轴上的每个值对应一个频率。在每个频率上,都会拍摄一张这样的 2d 图片。
我一直试图在 fit 文件中找到这些频率值,但到目前为止我还没有运气。我知道有 140 个不同的频道,我通过查看 fit 文件的标题知道单位等。问题在于找到这 140 个通道的值。
顺便说一句,在 Python 中编程。
它通常以与轴“n”的坐标相同的方式给出。参考值“CRVALn”、参考位置“CRPIXn”、每个像素/通道的增量“CDELTn”以及在这些关键字“CTYPEn”中给出的坐标类型。因此,您可以从中计算出立方体中每个“通道”的值
(channel_number - CRPIXn) * CRDELTn + CRVALn
既然您说您解决了它,也许您可以发布您的解决方案以造福他人?
编辑/添加:
有关如何在 FITS 文件中表达的更多信息:http: //fits.gsfc.nasa.gov/fits_wcs.html
特别是论文“FITS 中光谱坐标的表示”Greisen 等人。(2005 年)(NASA ADS 论文链接)