从 grib 天气模型中提取数据

Question

我已经grib1从GFS下载了模型数据，我使用的是 Mac OS X 并且能够wgrib2从NOAA构建文件，但无法提取数据。我还没有找到grib1适用于 MAC的实用程序。

然后我运行到这个程序http://www.giss.nasa.gov/tools/panoply/能够读取文件并可以看到数据，但我有一些问题。

1. 这些值有其他单位（如 K 表示温度，其他与降雨量的毫米不同）
2. 我可以导出CSV但只能导出值而不是坐标

我想做的是：

1. 通过经纬度通过命令行从grib文件中提取数据
2. 提取数据移动到MySQL数据库，以便能够进行一些数据聚合（总和、最大值、最小值）等
3. 能够绘制/映射数据（我可能会使用CartoDB服务）

我已经可以grib从某个地区提取或下载，我希望这些数据能够在电子表格上看到，例如纬度、经度、温度、降雨量、风等……然后我可以从那里把它带到数据库中，或 sum /avg 等或操作数据。

示例grib数据文件：https : //dl.dropboxusercontent.com/u/104462/neavilag_rain_wind_pressure.grb

根据我的情况，处理我的需求的最佳方法是什么？你能建议我怎么做吗？

Answer 1

我建议使用cfgrib基于 python 项目xarray和eccodesECMWF 软件的基于 python 的 grib 阅读器。

这是非常容易使用：

import cfgrib
cfgrib.open_datasets('your_grib_data.grib2')

而且 xarray 带有很多有用的内置函数。玩得开心！

Answer 2

如果我进行这种操作，我通常使用 CDO，因为它可以正确处理各种本机网格。例如，如果您想计算某个区域上某个字段的总和或平均值，则简单地使用常规经纬度或简化高斯网格计算算术平均值是不正确的，因为它没有考虑变化的网格 -当你向两极移动时，细胞的大小。

因此我会选择一个区域

cdo -f nc sellonlatbox,lon1,lon2,lat1,lat2 in.grb out.nc

这里我使用“-f nc”选项同时转换为netcdf。许多软件包可以轻松读取 netcdf 来创建绘图（ncl、gdl/idl、grads、R、python、panoply 等）。

要计算文件的字段统计信息，就像简单一样

cdo fldmean in.nc out.nc
cdo fldmax in.nc out.nc
cdo fldsum in.nc out.nc 

给出一些您要求的例子。您还可以在时间维度（timmean 等）进行统计。查看 CDO 帮助页面了解详细信息。我使用 macports 在我的 MAC 上安装了 CDO。