我有这个包含气象数据的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
打开文件,错误为:
无法打开指定的文件。