如何从FITS文件访问图像数据？

Question

如何从FITS文件访问图像子样本？

我无法确定文件中的数据到底在哪里！文件标题显示它位于第二个标题中，尺寸为1024 * 1024！但是当我尝试访问第二个标题时，我没有获得图像数组，而是获得了＃34; nonetype＆＃34;文件！我相信在这里做错了！

ModSecurity

Answer 1

即使你在这里访问第二个HDU：

hdu=headerList[1]

您没有显示该HDU的data，但您从第一个HDU中选择了imgData：

imgData = headerList[0].data
... 
plt.imshow(imgData)

修复只是显示hdu.data：

plt.imshow(hdu.data)

Answer 2

您可能希望使用fits.open()和fits.getheader('img.fits')，而不是使用fits.getdata('img.fits')。如果你需要的只是标题或图像数据，这些方法通常更方便。

在你的情况下

from astropy.io import fits
imgData = fits.getdata('AIA20100630_0000_0211.fits', 1)

应该直接加载图像数据而不通过hdu。