使用R语言读取grib2格式文件的方法.docx

使用R读取grib2文件HRThursday, October 08, 2015这个笔记整理了如何读取grib2格式文件。下面代码中的cmd字符串也可以在命令行工具中直接输入，效果相同。首先需要下载wgrib2命令行程序，这里是链接。将其中的wgrib2.exe和其他所有dll文件下载到某个本地路径下，并且将这个路径添加到系统环境变量中（可能需要重启生效）。完成上面步骤后，可以在windows打开cmd，输入wgrib2，看看能否运行。下面的例子中构造的cmd命令都是在R中用system函数执行的，也可以在windows的cmd中直接输入。也可以使用windows命令行的升级版powershell。下面的例子可以参考一下，特别是输入到csv文件。显示文件清单(inventory)# 这个fname是我放文件的路径fname - F:/tfiles/grib2_test/z500.ensm.jan.cfsv2.data.grb2cmd -paste(wgrib2 , fname)print(cmd)## [1] wgrib2 F:/tfiles/grib2_test/z500.ensm.jan.cfsv2.data.grb2inv -system(cmd, intern =TRUE)print(inv[1:10])## [1] 1:0:d=1982010618:HGT:500 mb:0-1 month ave fcst: ## [2] 2:49817:d=1982010618:HGT:500 mb:1-2 month ave fcst: ## [3] 3:99481:d=1982010618:HGT:500 mb:2-3 month ave fcst: ## [4] 4:148921:d=1982010618:HGT:500 mb:3-4 month ave fcst: ## [5] 5:198423:d=1982010618:HGT:500 mb:4-5 month ave fcst: ## [6] 6:247781:d=1982010618:HGT:500 mb:5-6 month ave fcst: ## [7] 7:297284:d=1982010618:HGT:500 mb:6-7 month ave fcst: ## [8] 8:338341:d=1982010618:HGT:500 mb:7-8 month ave fcst: ## [9] 9:379422:d=1982010618:HGT:500 mb:8-9 month ave fcst: ## [10] 10:428705:d=1982010618:HGT:500 mb:9-10 month ave fcst:例如上面清单中最后一行：280d=2009010618:HGT:500 mb:9-10 month ave fcst:它的格式含义如下：(msg|msg.submsg):(byte location):d=(YYYYMMDDHH):(variable):(level):(ftime):(optional attr)msg = message/record number starting from 1submsg = submessage number starting from 1 (optional)byte location: byte location of the start of the message starting from 0YYYYMMDDHH = reference time of the grib message, YYYY = year, MM = month (01..12), DD = day, HH = hourlevel = level/layer attributes, ex. 500 mb, ground, 2 m above groundftime = time attributes, ex. 12 hour fcst, anl, 0-6 hr ave fcstoptional attr: optional attributes that are needed to define the field such as ensemble memeber, chemical type (for concentration), probability range上面的例子中没有最后的optional attr。使用正则表达式提取一部分数据cmd -paste(wgrib2, fname, -match :(9-10 month ave fcst:))