from astropy.io import fits

if __name__ == "__main__":
    # 读取FITS文件
    filename = '/workspace/input/CSST_MSC_MS_SCI_20290807111557_20290807111827_10109300542959_08_L1_V01.fits'
    hdu_list = fits.open(filename)

    # 打印所有HDU（Header/Data Unit）的信息
    hdu_list.info()

    # 访问第一个HDU的数据
    primary_hdu = hdu_list[1]
    data = primary_hdu.data

    # 打印数据类型
    print("Data type:", type(data))

    # 访问第一个HDU的头部信息
    header = primary_hdu.header
    print("Header keywords:", list(header.keys()))

    # 关闭FITS文件以释放资源
    hdu_list.close()