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

Question

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

我尝试过使用解决方案Write 3d Numpy array to FITS file with Astropy，但是，它仍然不起作用。在

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文件，但在打开时显示一个错误“Cannot read fits file”。每当我试图打开从另一个软件创建的fits文件时，它会毫无问题地打开。在