问题标签 [grib]

我正在使用 R 处理气象数据。我分两步进行：

1. 使用 ncar 命令语言中的命令行函数 ncl_convert2nc 将 grib 转换为 netcdf
2. 在 R 中使用包 ncdf 导入 netcdf 数据。

我还有一个问题：

2- 对于某些特定的 grib 文件，使用 ncar 工具进行的转换不起作用。是否有其他方法或技巧（除了转录为 netcdf）来读取 R 中的 grib 文件？

Dirk 回答的问题： 1- 我想在 R 中自动处理许多文件。我可以在 R 中调用 ncl_convert2nc 吗？（由下面的 Dirk Eddelbuettel 回答）

这里有从事气象或地球科学工作的人吗？我正在尝试使用 GRIB2 格式的数据文件。有可用的库，特别是 Unidata GRIB Java 解码器。而且我可以确切地看到如何在一个大的线性数组中提取数据，但我想做的是按纬度/经度访问值。找不到任何简单的例子来说明如何做到这一点。建议？

蒂亚！

我正在尝试使用 R 函数从 ECMWF ERA-40 网站读取具有波高的 GRIB 文件wavedata.grib 。到目前为止，这是我的源代码：

wavedata.grib文件在 1957 年 9 月 1 日至2002年 8 月 31日期间具有网格 SWH 值。每个波段指的是一对纬度/经度，并且在每天的 00h 具有一系列 16346 SWH 值（1 波段 = 特定纬度/经度的 16346 个值）。

myframe 的尺寸为 16438 x 91。注意 91 = 7rows x 13columns。而16438这个数字几乎等于乐队的数量。额外的 2 行/波段是 long 和 lat 值，所有其他列应该是对应于 16436 波段的波高。

问题是我想在纬度/经度 = 43.75,331.25 处提取 SWH（波高），但它们与我在同一纬度/经度上使用 Degrib 实用程序读取文件的值不匹配。

此外，即使 myframe[16438,3]=43.75 (lat) 和 myframe[16437,3]= 331.25（长）。为什么是这样？我想知道 myframe[i,j] 值实际对应于哪个纬度/经度，或者在此过程中是否存在一些数据导入错误。我假设 Degrib 没有错误。

如果我想在网格点之间提取值，是否有任何 R 例程可以轻松地在矩阵中插入值？更一般地说，我需要帮助来编写一个有效的 R 函数来提取这样的波高：

请帮忙。

有人知道用于处理GRIB数据的 .NET 库吗？

有谁知道那里是否有一个处理 GRIB 文件的 php 库？我想使用加拿大环境部提供的 GRIB（或 GRIB 2）文件在我的网站上显示全球天气。我使用php来开发，这就是我寻找这个的原因。谢谢....

我有一个需要导入美国国家气象局表面分析的应用程序，这些分析以grib2文件的形式分发。我想将它们拉入 PostGIS 2.0 栅格，进行一些计算和建模，并在 GeoServer 中显示数据和模型结果。

由于 grib2 是 GDAL 支持的格式，因此提供的 raster2pgsql 实用程序应该能够将 grib2 直接插入到与 PostGIS 兼容的 SQL 中，一旦它存在，GeoServer 应该能够处理它。但是，我遇到了没有明显解决方案的问题——至少对我来说并不明显！Raster2pgsql 运行，显然没有错误，生成 SQL，运行 SQL 创建看起来非常像光栅的东西。但是 GeoServer 无法显示它——特别是边界看起来很奇怪 (0,0 -1,-1) 并且“预览层”只是抛出了 NullPointerException。

有人已经走过这条路了吗？我遇到了一些基本问题，比如不知道数据的 SRID 应该是什么（也许是 4326？）。我不希望任何人为我调试我的问题，但如果有人已经让这个工具链或其中的一部分工作，我可以插入已知好的东西，看看我能发现什么。

TIA,

读写

更新：Per Mike，这是其中一个文件中的坐标系内容；我省略了“gdalinfo”输出中的其他 749 个波段。请注意，文件名不同——我通过在原始文件上运行“gdalinfo”发现它有问题，gdalinfo 无法读取它。新（35MB！）文件在这里。

Gdalinfo 输出：

我已经下载了一天的“高分辨率初始条件”气候预报数据，它的扩展名是 .tar.gz，所以我将它解压缩到我的本地目录中，我得到了附件图像中的文件。我认为，没有扩展名的文件是 GRIB 数据（因为其中的第一个单词是“GRIB”）。所以我想从大文件（GRIB 和 NetCDF 格式，包含气候数据，如网格中的温度和压力）中获取数据到我的数据库，但它们是二进制的。你能推荐我一些从这些文件中获取数据的简单方法吗？我无法在他们的网站上获得有关处理他们的数据集的任何信息。

将这些文件转换为 .csv 会很好，但我找不到转换 GRIB 文件的程序。

如果我有一个 grib2 文件，其中包含整个世界的信息（对于某些参数），并且我想使用 wgrib2根据用户（客户端软件到服务器）给出的纬度和经度从中提取数据。我尝试了以下命令，但我只获得了完整的 grib2 文件：

wgrib2.exe input.grb -undefine out-box 10:90 -10:10 -grib output.grb

请告诉我我哪里出错了？谢谢。

我正在开发 wpf 应用程序。我正在使用以下代码运行 degrib.exe

degrib.exe 从 grib 文件生成 csv 文件。这http://www.nws.noaa.gov/mdl/degrib/index.php是degrib的解释。上面函数fileName中是grib文件的路径。如果 grib 文件放置在 Program Files 以外的任何文件夹中，则上述函数会从 grib 文件中创建 csv 文件。如果我将 grib 文件放在 Program Files 或 Program Files 中的任何其他文件夹中，相同的函数将不会生成 csv 文件。同样的功能也不适用于 AppData 文件夹。为什么会这样？您能为我提供上述问题的任何解决方案吗？

我正在尝试读取此 grib 文件，但出现错误。有什么方法可以解决此问题