我有一个代表 256x256 图像的 3D NumPy 数组(即 (10, 256, 256))。我想使用 astropy.io.fits 将此数组写入 FITS 文件,以便我可以使用 ds9 -mecube 打开文件并在帧中移动。我的尝试如下所示
export_array = numpy.array(images) #Create an array from a list of images
print export_array.shape ## (10, 256, 256)
hdu = fits.PrimaryHDU(export_array)
hdulist = fits.HDUList([hdu])
hdulist.writeto(out_file_name)
hdulist.close()
这会给我一个 FITS 文件,它实际上包含 3D 数组。但是,如果我用 ds9 -mecube 打开,我只能看到第一张图片。无论如何使用 astropy.io.fits 创建具有此功能的 FITS 文件?还是我缺少 ds9 的某些功能?