【问题标题】:Extracting data from grib file based on latitude and longitude根据经纬度从grib文件中提取数据
【发布时间】:2012-06-18 10:05:00
【问题描述】:

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

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

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

【问题讨论】:

    标签: grib


    【解决方案1】:

    我使用以下命令从 grib2 文件中提取信息。

    wgrib2.exe input_file.grib2 -lola LonSW:#lon:dlon LatSW:#lat:dlat file format

    假设我们有以下坐标供选择:

    顶部:(x0,y0) (x1,y0)

    底部:(x0,y1) (x1,y1)

    "LongSW"=x0, #lon = (x0~x1), "LatSW"=y0, #lat = (y0~y1)。 而dlon和dlan可以保持为1。'file'是输出文件名,格式可以是grib、csv、text等。

    在上面显示的命令中替换上面的值,你应该得到答案。

    【讨论】:

      【解决方案2】:

      如果您选择的经度是 360 度,纬度是 90 度:

      wgrib2.exe input_file.grb2 -lon 360 90 > output_file.txt

      【讨论】:

      • 你知道如何将它保存为 CSV 而不是 txt 吗?添加-csv 不会将文件限制为所传递的纬度/经度。
      • 输出到 CSV:wgrib2 intput.grb -s -irr_grid -121.4170309:47.45871 20 - | wgrib2 - -csv output.csv
      • 这个@ScottChamberlin 应该是一个可接受的答案。 -ijsmall_grib 而不是 if -irr_gird 如果您真的要“提取”而不是通过 wgrib2 进行插值
      【解决方案3】:

      我认为你也可以直接在 grib 上使用 cdo

      cdo sellonlatbox,lon1,lon2,lat1,lat2 in.grb out.grb
      

      如果 grib 文件位于简化的高斯网格上,您可能需要指定您想要常规的经纬度输出。我通常使用“-f nc”自己将输出格式转换为 netcdf,因为我发现在其他软件中处理起来更容易。

      【讨论】:

        【解决方案4】:

        这仍然是谷歌上的热门,所以即使它有点老,这里有一个更详细的解释。

        首先,您需要来自National Weather Service Climate Prediction Centrewgrib2。 (它的安装很简单,但解释得不太清楚。请参阅this page,或此gist 寻求帮助。)

        接下来,您需要使用lola 函数(用于经纬度网格)。

        你需要给wgrib2几个参数:

        • 包含数据的 grib 文件
        • 经度信息:
          • 边界框西南角的经度
          • 边界框的长度(以度为单位)
          • 此方向上的点间距
        • 纬度信息,同上
          • 边界框西南角的纬度
          • 边界框的长度(以度为单位)
          • 此方向上的点间距
        • 要写入的文件名
        • 要写入的文件格式(bin 用于二进制,text 用于简单文本,spread 用于电子表格格式,grib 用于 grib2)。

        例如:

        wgrib2 input.grb 220:100:1 20:50:1 output.grb2 grib
        

        将创建一个在两个方向上以 1 度间隔覆盖北美(220 E - 320 E;20 N - 70 N)的输出文件。

        罗伯

        【讨论】:

          猜你喜欢
          • 1970-01-01
          • 2016-10-21
          • 1970-01-01
          • 2016-11-26
          • 1970-01-01
          • 1970-01-01
          • 1970-01-01
          • 2011-04-01
          • 2016-09-07
          相关资源
          最近更新 更多