python中对fits文件读取和保存

天文图片一般都是以.fits的格式保存的。在Windows上可以通过​​​​​​​SAOImage DS9​​​​​​  读取。​​

安装包 astropy

sudo pip3 install astropy

图片的读取：

from astropy.io import fits
import os
hud = fits.open('/home/yellow/1.fits')
img = hud[0].data

图片的保存

path = '/home/yellow/'
if os.path.exists(path + '2.fits'):
    os.remove(path+'2.fits')
grey=fits.PrimaryHDU(img)
greyHDU=fits.HDUList([grey])
greyHDU.writeto(path+'2.fits')