python - 在 python 中折叠/展平 FITS 数据立方体

Question

我找遍了所有地方，但没有找到解决此问题的方法。我觉得它应该相当简单，但我们拭目以待。

我有一个 .FITS 格式的数据立方体，我需要将它折叠成一个 2D FITS 图像。数据立方体有两个空间维度和一个频谱/速度维度。

只是寻找一个简单的 python 例程来加载立方体并展平所有这些层(即沿光谱/速度轴整合它们)。感谢您的帮助。

Answer 1

This tutorial on pyfits有点老了，但还是基本正确的。关键是用 pyfits(或 astropy.io.fits)打开 FITS 立方体的输出是您有一个 3 维 numpy 数组。

import pyfits
# if you are using astropy then for this example
# from astropy.io import fits as pyfits

data_cube, header_data_cube = pyfits.getdata("data_cube.fits", 0, header=True)
data_cube.shape
# (Z, X, Y) 

然后您必须决定如何沿 Z 轴展平/整合立方体，并且有大量资源可以帮助您决定正确的做法(希望基于某些分析框架)。