python - 如何从 numpy 数组创建 .fits 文件

Question

我目前正在尝试将 .cr2 转换为 .fits 文件。我目前无法创建 .fits，因为我创建的任何文件由于错误的配合形状而无法打开或工作。

我已经尝试使用解决方案将 3d Numpy 数组写入带有 Astropy 的 FITS 文件，但仍然无法正常工作。

import numpy as np
from astropy.io import fits
im = imageio.imread(r'E:\FYP\cr.cr2', format="RAW-FI")
a = np.array(im)
print (a.shape)
hdu = fits.PrimaryHDU(a)
hdulist = fits.HDUList([hdu])
hdulist.writeto(r'E:\FYPcr.fits')

它成功创建了 .fits 文件，但是在打开“无法读取适合文件”时出现错误。每当我尝试打开从另一个软件创建的 fit 文件时，它都会毫无问题地打开。

标签： pythonnumpyastropyfitspython-imageio