python reader agri satpy_Satpy基础系列教程(1)-FY4A AGRI L1数据处理

Satpy[1]目前支持的卫星数据[2]有50种(MSG, Himawari 8, GOES-R, MODIS, Sentinel- 1/2/3/5, SNPP等)。

本文以最近朝曦dawn[3]添加的风云4A(FY4A) AGRI L1数据为例。

Notebook[4]已放在GitHub上，供大家学习。

FY4A AGRI L1数据有两种类别：

1.全圆盘

FY4A-_AGRI--_N_DISK_1047E_L1-_FDI-_MULT_NOM_20190807060000_20190807061459_4000M_V0001.HDF

2.中国区

FY4A-_AGRI--_N_REGC_1047E_L1-_FDI-_MULT_NOM_20190807045334_20190807045750_1000M_V0001.HDF

全圆盘标识：DISK , 中国区标识：REGC.

风云数据链接

1.实时数据[5](30天内)及文档

2.历史数据[6] (2018-03-12 -- )

3.FY4A官方应用平台[7]

安装

推荐通过conda安装，简单快捷。

$ conda install -c conda-forge satpy

定标

FY4A AGRI L1有以下三种定标方式:

1.原始数字量化值 (所有通道)

2.反射率 (C01 - C06)

3.辐射率和亮温 (C07 - C14)

读取数据

只需几行即可读取数据：

import os, globfrom satpy.scene import Scene# 加载FY4A文件filenames = glob.glob('/xin/data/FY4