- 积分
- 2652
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-8-6
- 最后登录
- 1970-1-1
![[Eagle__Z] 粉丝数:136 微博数:880 新浪微博达人勋](source/plugin/sina_login/img/light.png)
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
相信很多人都需要这个知识点,我就把今天的经验做个介绍,可能不是很完善,希望能帮到你。
1、之前我在家园帖子也发过相似的帖子,windows读取grib2,那时候我用的是OLDLee 河北干部培训学院李开元老师的方法
用wgrib2转,具体的大家可以在家园里搜索哟。
2、之前在群里有人问怎么用xarray+cfgrib读取,那时候刚好看到有人用这样读取,具体方法就没做介绍,刚好今天上完大夜班休息,两个小伙伴来问我,
我就搜了一下,也总结一下了。
3、第一步:xarray 相信这个没问题 conda环境底下conda install xarray
4、第二部:安装cfgrib,由于该模块依赖于ecmwf eccode 二进制库 必须安装在系统上并作为共享库访问。安装eccodes,
具体方法参照:https://www.cnpython.com/pypi/cfgrib 这里有说明,我在安装的过程中,实现没提前创建环境,导致numpy包冲突,
所以建议大家创建一个环境安装吧。安装好了以后再conda install cfgrib。
5、安装完成以后,import cfgrb 看看是否安装成功,我自己也是安装了两遍才成功的。
6、https://cloud.tencent.com/developer/article/1471285 学前班大队长在这里已经做了很好的总结了,大家可以参照一下。
7、在上面的帖子中,介绍了如何提取某一层的要素,那么我们在实际过程中,还会遇到画垂直剖面时,需要提取某一要素的所有层次资料,
经过摸索找到问题,直接把代码贴上来 大家自行调试吧。
ds = xr.open_dataset("E:/data/fnl_20190526_12_00.grib2",engine='cfgrib',
backend_kwargs={'filter_by_keys':{'typeOfLevel':'isobaricInhPa','cfVarName': 'w'}})
在这里 我下载的fnl再分析资料中,w代表的是垂直速度,当然还有t\gh\u\v\r等等,用cfVarName匹配对应的物理量。
大概就是这些吧 大家在上面的链接中先学习吧,有问题再欢迎交流。
祝您学习愉快,生活安康。
|
评分
-
查看全部评分
|