1926人阅读
R语言（8）
下面是用R语言读取grib2文件的代码。使用的文件是CFS模式输出的grib2格式的文件。
思路是首先将grib2文件转换成ncdf文件，然后读取ncdf文件。使用的package为ncdf。
注意，运行下面的代码，首先需要将wgrib2相关的程序放在系统path中，相关程序在下面的页面中下载：
首先加载ncdf包：
library(ncdf)
下面的命令将多个记录转化为nc文件mync.nc，其中读取了从第3个开始到第200个每隔10个读1个的记录：
shell('wgrib2 z500.ensm.apr.cfsv2.data.grb2 -for &3:200:10& -netcdf mync.nc')
程序显示：
3:98936:d=:HGT:500 mb:2-3 month ave fcst:
13:577923:d=:HGT:500 mb:2-3 month ave fcst:
23:1048267:d=:HGT:500 mb:2-3 month ave fcst:
33:1535324:d=:HGT:500 mb:2-3 month ave fcst:
43:2006379:d=:HGT:500 mb:2-3 month ave fcst:
53:2485247:d=:HGT:500 mb:2-3 month ave fcst:
63:2955864:d=:HGT:500 mb:2-3 month ave fcst:
73:3434510:d=:HGT:500 mb:2-3 month ave fcst:
83:3905217:d=:HGT:500 mb:2-3 month ave fcst:
93:4375576:d=:HGT:500 mb:2-3 month ave fcst:
103:4846312:d=:HGT:500 mb:2-3 month ave fcst:
113:5317129:d=:HGT:500 mb:2-3 month ave fcst:
123:5796269:d=:HGT:500 mb:2-3 month ave fcst:
133:6267130:d=:HGT:500 mb:2-3 month ave fcst:
143:6745573:d=:HGT:500 mb:2-3 month ave fcst:
153:7208123:d=:HGT:500 mb:2-3 month ave fcst:
163:7678480:d=:HGT:500 mb:2-3 month ave fcst:
173:8140100:d=:HGT:500 mb:2-3 month ave fcst:
183:8610887:d=:HGT:500 mb:2-3 month ave fcst:
193:9073298:d=:HGT:500 mb:2-3 month ave fcst:
然后使用ncdf包中的函数读取即可：
nc &- open.ncdf(&mync.nc&)
data &- get.var.ncdf( nc )
&&相关文章推荐
参考知识库
* 以上用户言论只代表其个人观点，不代表CSDN网站的观点或立场
访问：53286次
排名：千里之外
原创：15篇
评论：22条
(1)(2)(3)(5)(3)(1)用proto命令生成java文件
/p/protobuf/downloads/list 下载 protoc-2.5.0-win32，然后解压，配置protoc环境变量
2.cmd命令 &输入protoc --version查询版本好，如果查到版本号，说明可以。
libprotoc 2.5.0
3.在proto.exe同一级文件夹下写.proto文件，例如：msg.proto
message SearchRequest {
required string query = 1;
optional int32 page_number = 2;
optional int32 result_per_page = 3;
3.用cd命令打开文件所在目录，例如文件在E:\tools\Computer\protoc-2.5.0-win32下
E:\&cd tools
E:\tools&cd Computer
E:\tools\Computer&cd protoc-2.5.0-win32
E:\tools\Computer\protoc-2.5.0-win32&protoc.exe --java_out=./ msg.proto温馨提示！由于新浪微博认证机制调整，您的新浪微博帐号绑定已过期，请重新绑定！&&|&&
LOFTER精选
网易考拉推荐
用微信&&“扫一扫”
将文章分享到朋友圈。
用易信&&“扫一扫”
将文章分享到朋友圈。
阅读(4608)|
用微信&&“扫一扫”
将文章分享到朋友圈。
用易信&&“扫一扫”
将文章分享到朋友圈。
历史上的今天
loftPermalink:'',
id:'fks_',
blogTitle:'grib2 数据读取，转nc',
blogAbstract:'wget
{if x.moveFrom=='wap'}
{elseif x.moveFrom=='iphone'}
{elseif x.moveFrom=='android'}
{elseif x.moveFrom=='mobile'}
${a.selfIntro|escape}{if great260}${suplement}{/if}
{list a as x}
推荐过这篇日志的人：
{list a as x}
{if !!b&&b.length>0}
他们还推荐了：
{list b as y}
转载记录：
{list d as x}
{list a as x}
{list a as x}
{list a as x}
{list a as x}
{if x_index>4}{break}{/if}
${fn2(x.publishTime,'yyyy-MM-dd HH:mm:ss')}
{list a as x}
{if !!(blogDetail.preBlogPermalink)}
{if !!(blogDetail.nextBlogPermalink)}
{list a as x}
{if defined('newslist')&&newslist.length>0}
{list newslist as x}
{if x_index>7}{break}{/if}
{list a as x}
{var first_option =}
{list x.voteDetailList as voteToOption}
{if voteToOption==1}
{if first_option==false},{/if}&&“${b[voteToOption_index]}”&&
{if (x.role!="-1") },“我是${c[x.role]}”&&{/if}
&&&&&&&&${fn1(x.voteTime)}
{if x.userName==''}{/if}
网易公司版权所有&&
{list x.l as y}
{if defined('wl')}
{list wl as x}{/list}}