提取第三轴，数据立方体 pyfits
我有一个数据立方体的fits文件，这个文件里包含了一个规则的（512x512）二维网格上的强度值。不过，这里还有第三个轴；这个轴上的每个值对应一个频率。在每个频率下，都会拍摄这样一幅二维图像。
我一直在尝试在fits文件中找到这些频率值，但到目前为止都没有成功。我知道有140个不同的频率通道，通过查看fits文件的头部信息，我也知道单位等信息。问题是如何找到这140个通道的值。
顺便提一下，我是在用PYTHON编程。
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的论文)