爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6488|回复: 11

读取日本55年气候再分析资料,用matlab或者IDL读取grib,读出一维向量后如何变2维?

[复制链接]

新浪微博达人勋

发表于 2015-1-26 12:25:04 | 显示全部楼层 |阅读模式
GrADS
系统平台: win7 matlab IDL
问题截图:
问题概况: 可不可以请教一下,我有一个grib文件ctl显示x轴是320,y轴是640,空间分辨率是0.562500,fltarry是157792*1,怎么样可以将这个一维的数据显示为正确的二维呢? data=reshape(data,[640,320])对报错说总数不对,求教大神是哪里的问题?
我看过提问的智慧: 看过
自己思考时长(天): 3

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
可不可以请教一下,读取日本55年气候再分析资料grib文件,用matlab或者IDL读取,ctl显示x轴是320,y轴是640,空间分辨率是0.562500,fltarry是157792*1,怎么样可以将这个一维的数据显示为正确的二维呢? data=reshape(data,[640,320])对报错说总数不对,求教大神是哪里的问题? QQ图片20150126121513.png QQ图片20150126121537.png QQ图片20150126121613.jpg QQ图片20150126121633.jpg
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-1-26 14:18:44 | 显示全部楼层
有没有大神路过了~~~烦请指点一下呀~~~~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-26 14:19:06 | 显示全部楼层
你这应该是读二进制的时候设置的数据精度不对
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-1-26 14:24:02 | 显示全部楼层
asmon 发表于 2015-1-26 14:19
你这应该是读二进制的时候设置的数据精度不对

可不可以求大神再给详细解释下......刚开始接触这个,不是很懂.........代码是网上找的,读出来就是这个一维数组了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-26 14:32:16 | 显示全部楼层
weinigod 发表于 2015-1-26 14:24
可不可以求大神再给详细解释下......刚开始接触这个,不是很懂.........代码是网上找的,读出来 ...

matlab读取二进制数据用的是fread命令,你查看一下它的说明
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-26 16:47:04 | 显示全部楼层
飘过。。。
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2015-1-26 20:02:43 | 显示全部楼层
asmon 发表于 2015-1-26 14:32
matlab读取二进制数据用的是fread命令,你查看一下它的说明

还是不行.....我读取数据的时候用过的命令就是
;grib_struct=read_grib('d:\anl_surf125\anl_surf125.clim8110.mon12','invent')
>> fid=read_grib(['d:\anl_surf125\anl_surf125.clim8110.mon12'],3);
data=getfield(fid,'fltarray')
而且用下载IDL的代码读出来也是157792
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-1-26 21:01:59 | 显示全部楼层
本帖最后由 asmon 于 2015-1-26 22:18 编辑
weinigod 发表于 2015-1-26 20:02
还是不行.....我读取数据的时候用过的命令就是
;grib_struct=read_grib('d:\anl_surf125\anl_s ...

这个就不清楚了,你用的这个命令我没用过。
你有ctl的话可以直接用grads读.或者你matlab试试format命令
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-2-6 11:03:41 | 显示全部楼层
试试这个
grib_struct=read_grib('d:\anl_surf125\anl_surf125.clim8110.mon12','invent');
fid=read_grib('d:\anl_surf125\anl_surf125.clim8110.mon12',3,'ScreenDiag',0);
a=flipud(reshape(fid.fltarray,[640 320])');

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-2-27 16:32:19 | 显示全部楼层
顶一个
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表