读取FITS文件

Python 读取FITS文件

FITS 文件是天文图像使用的存储格式，一般读取使用库是：

from astropy.io import fits

读取fits：

hdu = fits.open(path)
data1 = hdu[0].data  #得到数据

一种批量读取文件夹下多个文件中的FITS文件的方式，使用os:

import os
import numpy as np
from astropy.io import fits
path = '/home/data/'  #文件所在的位置
filelist1 = os.listdir(path)
data = {}
for i in range(len(filelist1)):
    a = filelist1[i]
    path1 = path+a
    filelist2 = os.listdir(path1)
    image =[]
    for j in range(len(filelist2)):
        hdu = fits.open(path1+'/'+filelist2[j])
        data1 = hdu[0].data
        image.append(data1)
    c = np.array(image)
    data[a] = c   
np.save('data.npy')