astropy对fits文件的基础操作

astropy读取fits文件汇总
from astropy.io import fits
hdu_list = fits.open(image_file)
hdu_list.info()
image_data = hdu_list[0].data

print(type(image_data))

print(image_data.shape)
hdu_list.close() # 关闭，否则占用内存

如果你不需要查看fits头文件，可以用fits.getdata代替前面的步骤

image_data = fits.getdata(image_file)
print(type(image_data))
print(image_data.shape)
plt.imshow(image_data, cmap=‘gray’) # 不指定颜色，默认为彩色
plt.colorbar()
#打印最大值最小值平均值和标准差
print(‘Min:’, np.min(image_data))
print(‘Max:’, np.max(image_data))
print(‘Mean:’, np.mean(image_data))
print(‘Stdev:’, np.std(image_data))

绘出柱状图

NBINS = 1000
histogram = plt.hist(image_data.flatten(), NBINS)

hdu=header data unit
hdul.info()
hdul[0].header
hdul[0].header[‘DATE’]
hdul[0].header.comments[‘fitsver’] #显示关键字的注释
list(hdul[0].header) # 显示所有关键字
data = hdul[1].data
data.shape
data[0] # 第一行的数据
data.field(0) # 第一列的数据
data.field(‘utobs’)
data.field(‘freq’)
data.field(‘CHAN_BW’)
hdul[1].header
data[2][3] # 访问第二行第三列数据

data.names # 查看字段，即field
data.field(‘DATA’)

hdr[‘targname’] = (‘NGC121-a’, ‘the observation target’) 为头文件添加关键字，内容和注释

参考文献：https://www.cnblogs.com/cuiastro/p/13525114.html