如何从 Grib1/.grb 文件中转储值答案

【问题标题】：How to dump values from a Grib1/.grb file如何从 Grib1/.grb 文件中转储值
【发布时间】：2014-09-02 07:12:46
【问题描述】：

我想知道，有没有办法从 grib1 文件中转储值？我的最终目标是在纬度和经度或至少一个网格点上找到单个消息的值。我使用的是linux系统。 Wgrib 似乎除了读取有关消息的元数据或重建消息之外什么都不做。

我懂一点python，所以我可以使用pygrib，但我不知道如何提取特定纬度和经度的值。

感谢您的回答，

【问题讨论】：

【解决方案1】：

如果您对来自 NOMADS 的数据感兴趣，我建议您通过他们的 THREDDS 数据服务器，这将允许您通过指定纬度/经度来访问数据，并且您可以将该数据作为 csv 文件取回，如果您希望。为此，请先访问 NOMADS TDS 网站：

您链接到的示例数据文件可以在这里找到：

找到一个你感兴趣的网格，比如 18Z 运行分析字段：

点击“NetcdfService”链接：

在该页面顶部附近，单击“作为点数据集”：

然后检查您感兴趣的参数，纬度/经度（将选择最接近该纬度/经度的网格点）以及输出格式类型。

这个网页界面基本上会生成一个访问 url，如果我想要温度超过 Boulder, CO，以 csv 格式返回，它是这样的：

从上面的 URL 可以看出，您可以很容易地生成这些并发出请求，而无需执行上述所有步骤。

这种访问方法（NetcdfSubsetService）可以很容易地与 Python 结合使用。例如，查看这个 ipython 笔记本：

特别是笔记本中的第一个和第二个单元格。

请注意，您可以在以下位置获取最近的 GFS 数据，其中整个模型运行包含在一个 grib 文件中：

这将允许您发出一个请求，就像上面的请求一样，但多次使用一个请求。

干杯，

肖恩

【讨论】：

【解决方案2】：

您可以使用 grib 工具，特别是 grib_ls 和 grib_get 从最接近指定纬度和经度的 1 个网格点或 4 个网格点获取值。因此，您可以使用最近邻或双线性插值或任何您喜欢的方法来获得价值。阅读此演示文稿，grib_ls 从第 31 页开始：

当您安装 grib 工具时，您将获得一些工具来帮助您使用 GRIB 文件。

【讨论】：