我找遍了所有地方,但没有找到解决此问题的方法。我觉得它应该相当简单,但我们拭目以待。
我有一个 .FITS 格式的数据立方体,我需要将它折叠成一个 2D FITS 图像。数据立方体有两个空间维度和一个频谱/速度维度。
只是寻找一个简单的 python 例程来加载立方体并展平所有这些层(即沿光谱/速度轴整合它们)。感谢您的帮助。
最佳答案
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 轴展平/整合立方体,并且有大量资源可以帮助您决定正确的做法(希望基于某些分析框架)。
关于python - 在 python 中折叠/展平 FITS 数据立方体,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/17980365/