fits文件的读取与查看

fits文件读取与查看（写给自己，大神勿喷）

fits文件的结构

一个fits文件包括多个HDU（两个以上），开头一个是Primary HDU，储存着类似于观测日期等的信息，剩下的是Bin HDU，储存着观测数据

每个HDU由header 和data组成，data中储存着数据，header表示相应的信息，储存着和Primary HDU一样的信息，并且描述第几列的数据是什么，类型是什么

TTYPE1=‘obsid’表示第一列的数据是obsid，TFORM1和TDISP1表示第一列的数据的格式（format）

python读取fits文件

使用astropy包读写fits文件（astropy需要下载）

from astropy.io import fits 

#读取文件
lamost=fits.open('dr5_stellar.fits')
#查看HDU
lamost.info()
#查看第二个HDU的header
lamost[1].header
#查看第二个HDU的data
lamost[1].data
#按照关键字搜索第二个HDU的数据（例子中是查看lamost数据的赤经）
lamost[1].data.field('ra')

astropy的官方文档

python关闭fits文件

及时关闭文件以减小内存消耗

lamost.close()