如何打开GRIB文件

时间:2017-09-05 05:42:58

标签: r grid arcgis

我有这个包含气象数据的GRIB文件:alaro.grb

我试图使用不同的软件打开此文件但到目前为止没有运气。这是我试过的:

使用R:

library(rgdal)
gribfile <- "alaro.grb"
grib <- readGDAL(gribfile)
  

alaro.grb是一个grib文件,但没有成功识别栅格数据集。

library(rNOMADS)
GribInfo(gribfile)

> "grib1 message ignored (use wgrib)"

所以我认为它是一个GRIB-1文件。

GribInfo(gribfile, file.type = "grib1")
  

系统错误(paste0(&#34; wgrib&#34;,grib.file),&#34; -s&#34;,intern = TRUE):      &#39; ignore.stdout&#39;必须为TRUE或FALSE

我直接在wgrib窗口中使用cmd来收集信息:

wgrib alaro.grb -V

导致

  

rec 1:0:date 2014061306 TMP kpds5 = 11 kpds6 = 105 kpds7 = 2 levels =(0,2)grid = 255 2 m以上gnd 2hr fcst:     TMP =温度。 [K]     时间范围0 P1 2 P2 0 TimeU 1 nx 508 ny 446 GDS网格0 num_in_ave 1缺失0     中心224子中心0过程1表1扫描:WE:SN风(网格)     latlon:纬度37.310000至55.110000乘0.040000 nxny 226568             long 2.000000 to 32.420000 by 0.060000,(508 x 446)scan 64 mode 136 bdsgrid 1     最小/最大数据276.842 304.709 num bits 12 BDS_Ref 276.842 DecScale 0 BinScale -7

     

rec 2:340424:date 2014061306 APCP kpds5 = 61 kpds6 = 1 kpds7 = 0 levels =(0,0)grid = 255 sfc 1-2hr acc:     APCP =总降水量[kg /平方公尺]     时间范围4 P1 1 P2 2 TimeU 1 nx 508 ny 446 GDS网格0 num_in_ave 1缺失0     中心224子中心0过程1表1扫描:WE:SN风(网格)     latlon:纬度37.310000至55.110000乘0.040000 nxny 226568             long 2.000000 to 32.420000 by 0.060000,(508 x 446)scan 64 mode 136 bdsgrid 1     最小/最大数据0 7.13696 num位16 BDS_Ref 0 DecScale 0 BinScale -13

     

rec 3:794132:date 2014061306 SNMR kpds5 = 171 kpds6 = 1 kpds7 = 0 levels =(0,0)grid = 255 sfc 2hr fcst:     SNMR =雪混合比[kg / kg]     时间范围0 P1 2 P2 0 TimeU 1 nx 508 ny 446 GDS网格0 num_in_ave 1缺失0     中心224子中心0过程1表1扫描:WE:SN风(网格)     latlon:纬度37.310000至55.110000乘0.040000 nxny 226568             long 2.000000 to 32.420000 by 0.060000,(508 x 446)scan 64 mode 136 bdsgrid 1     最小/最大数据0 1 num位8 BDS_Ref 0 DecScale 0 BinScale -7

     

rec 4:1021272:date 2014061306 CNWAT kpds5 = 223 kpds6 = 1 kpds7 = 0 levels =(0,0)grid = 255 sfc 1-2hr acc:     CNWAT =植物冠层表面水[kg / m ^ 2]     时间范围4 P1 1 P2 2 TimeU 1 nx 508 ny 446 GDS网格0 num_in_ave 1缺失0     中心224子中心0过程1表1扫描:WE:SN风(网格)     latlon:纬度37.310000至55.110000乘0.040000 nxny 226568             long 2.000000 to 32.420000 by 0.060000,(508 x 446)scan 64 mode 136 bdsgrid 1     最小/最大数据27039.5 2.87254e + 06 num位16 BDS_Ref 27039.5 DecScale 0 BinScale 6

有了这些信息,我回到了R并尝试了:

ReadGrib(gribfile, variables = "CNWAT", levels = c(0,0), file.type = "grib1")
> $meta.data
NULL

$value
NULL

$variables
NULL

$levels
NULL

$grib.type
[1] "grib1"

所以似乎没有数据wgrib可以阅读。

我还尝试使用ArcGIS打开文件,错误为:

  

无法打开指定的文件。

0 个答案:

没有答案