python 读取grib/grib2格式数据
一般利用NCL(NCAR Command Language)读取。利用Python语言的pygrib库也可以读取grib/grib2格式数据
grib/grib2 格式数据介绍
1 数据下载方法
数据可以通过NCAR/UCAR的官网FNL再分析资料进行下载
注册账号–找下下载的数据–下载相应的下载shell脚本–在服务器中运行该脚本进行下载
2 数据格式
fnl再分析资料的数据格式有netcdf 和grib grib2 这几种数据格式,用于WRF读取的为grib grib2格式的,所以一般下载grib格式的数据。其中,grib1 的数据从1999年7月30日-2007年12月6日;grib2格式的数据从2007年12月7日到现在
3 数据变量
可以利用Panoply load 相应的grib2格式的数据查看变量
可以利用wgrib读取变量
可以利用pygrib函数读取变量【下面进行讲述】
变量对应的TABLE不能在文件中查看,只能在官网上寻找相应的TABLE grib2格式数据变量TABLE
pygrib 函数介绍
所有的pygrib版本都在Pythonpygrib官网中有介绍,包括其版本,安装包和安装方法。
1 头文件读取
头文件相关信息读取
import pygrib as pg
grbs = pg.open('fnl_17030406.grib2')# 所有变量
for grb in grbs:
print grb #每一个变量的头文件
print grb.keys() #每一个变量的keys
print grb.values # 每一