Windows10下处理grib数据

2019-05-06  本文已影响0人  沐辰老爹

2019年9月10日修改!
以下安装方式采用离线安装,主要考虑写博文时anaconda的清华镜像无法使用,目前清华镜像已获得anaconda授权重新开张,因此下述离线安装方式已不可取,读者可自行跳过。
通过清华镜像帮助来配置好channel,之后直接使用如下命令'conda install eccodes'即可!

介绍

grib1和grib2数据一直以来都是气象领域重要的数据格式,尤其对于数值预报的传输而言尤其重要,其他相关内容可进行百度。
今天介绍的是如何在Windows10下利用欧洲中心的eccodes库)进行grib数据的读取和可视化。

安装

首先我们这里主要使用2个Python库:

当然还包括eccodes

eccodes

使用anaconda平台的情况下可以使用conda install eccodes来进行安装,但是随着国内主要镜像关闭,下载和安装速度非常慢。所以我们去官网搜索eccodes并下载文件进行离线安装。

set "ECCODES_SAMPLES_PATH=c:\ProgramData\Anaconda3\Library\share\eccodes\samples"
set "ECCODES_DEFINITION_PATH=c:\ProgramData\Anaconda3\Library\share\eccodes\definitions"
至此我们算是完全安装可以使用eccodes了!

cfgrib

xarray

使用

首先需要下载示例数据,完成后保存到本地供后续使用。

读取

启动Python交互界面并输入下面的示例命令。

import xarray as xr
filein = 'D:\\Download\\era5-levels-members.grib'
image.png

导入相关库,下面读取数据并打印简表

data = xr.open_dataset(filein, engine='cfgrib')
data
image.png

至此我们看到已经读入到xarray的dataset类型下了。

可视化

好了,下面进行可视化。

import matplotlib.pyplot as plt
data.isel(number=1, time=1, isobaricInhPa=1).t.plot()
plt.show()
image.png

如此会弹出对话框如下:


成功

结语

错别字不该了,就酱。

上一篇 下一篇

猜你喜欢

热点阅读