使用python读取GRIB文件

ztftom

到现在为止，我所知道的使用python读取grib文件的方法有两个，都需要安装额外的package：
1. pygrib
http://code.google.com/p/pygrib/
基于ECMWF GRIB API。 这个package的安装我始终没有成功，大家可以继续尝试一下。

2. pyNIO
http://www.pyngl.ucar.edu/Nio.shtml
我现在在用的一个软件包，下载解压缩之后放在EPD python的library里面就可以调用了。
=========下面是一个小样例======
import numpy as np      # pyNio is based on Numpy
import Nio                    # import pyNio
fname = '~/data/fnl'       # replace with real filename
f     = Nio.open_file(fname, mode='r',\
         options=None, history='', format='grib')   # open the file
var   = f.variables[var_name]   # build up the pointer to the variable
var_value = var.get_value()     # get the value of the variable
f.close()                                 # close file
=======  END ========

GRIB文件需要独立的索引表才能读取。pyNio对主流的reanalysis data有内建的索引表（基于wgrib）。具体细节可参考上面给出的链接。下面是从reanalysis的GRIB中读到的SST画出的图：

何永利

nc文件应该也差不多吧

何永利

貌似没找到windows下的pynio

ztftom

何永利 发表于 2012-3-14 13:46

                               
登录/注册后可看大图

貌似没找到windows下的pynio

恐怕要用Cygwin

(``)nnnnnn..

                                            

海蓝蓝

有没有用fortran读取grib的？

kongfeng0824

最近处理grib，不过用的是matlab

suzhongyue

何永利 发表于 2012-3-14 13:46

                               
登录/注册后可看大图

貌似没找到windows下的pynio

好像我也要用python读nc文件了，结果就看到你的ID

zhs1019569837

{:eb502:}不错啊

天涯王丹

多谢楼主分享啊。