grib1文件解析 python_python读取grib格式数据

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 # 每一个变量的值

头文件读取并写入txt脚本

import pygrib as pg

grbs = pg.open('fnl_