我找遍了所有地方，都没有找到解决这个问题的办法。我觉得应该很直截了当，但我们拭目以待。
我有一个.FITS格式的数据立方体，我需要把它折叠成二维FITS图像。数据立方体有两个空间维度和一个光谱/速度维度。
只需寻找一个简单的Python例程，加载到立方体中并将所有这些层（即沿着光谱/速度轴积分）。谢谢你的帮助。
最佳答案
This tutorial on pyfits有点老了，但基本上还是正确的。关键是用pyfits（或astropy.io.fits）打开FITS立方体的输出是您有一个三维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轴平移/整合立方体，并且有足够的资源来帮助你决定正确的（希望在一些分析框架中）这样做。
关于python - 在python中折叠/展平FITS数据立方体，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/17980365/