爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
12
返回列表 发新帖

[程序设计] matlab读取grd错误

[复制链接]
发表于 2013-12-30 11:02:24 | 显示全部楼层
楼主你好 ,我想利用matlab里的detrend来扣除我资料的线性趋势,最后输出grd,您能够帮忙看一下这段代码有没有问题吗 谢谢
>> fid=fopen('e:\lww\u850jp.grd','rb');
>> ssh=fread(fid,inf,'float32');
>> fclose(fid);
>> ssh=reshape(ssh,81,9,744);
>> for i=1,81
for j=1,9
for k=1,744
x(k)=ssh(i,j,k)
end
for k=1,744
y(i,j,k)=detrend(x)
end
end
end
密码修改失败请联系微信:mofangbao
发表于 2014-1-11 15:28:55 | 显示全部楼层
哈哈,虽然之前在网上也找到了读取grd格式的程序,但是一直都读不出来。看到楼主的程序以后,受到了启发,找到了问题所在,现在能够读出来啦。。。。嘿嘿,谢谢啦~~~~
密码修改失败请联系微信:mofangbao
发表于 2016-5-9 10:11:39 | 显示全部楼层
你们有人用micaps来读grd的文件吗,我读不出来
密码修改失败请联系微信:mofangbao
发表于 2016-7-26 17:56:04 | 显示全部楼层
你好请问你的数据读取是怎么解决的,我的读取方式和你的是一样的但是,好像数据统计后有些错误
密码修改失败请联系微信:mofangbao
发表于 2016-8-4 13:01:20 | 显示全部楼层
楼主,你好我写的读取程序和你的是一样的为甚我读出来的数据有错误?
  1. clear;clc;
  2. fclose all;
  3. fid=fopen('H-2008071418.grd','rb');
  4. data=fread(fid,inf,'float');
  5. data1=reshape(data,700,440,1,2);
  6. data2=reshape(data1(:,:,1,1),700,440);
  7. fclose all;
复制代码

这是我写的读取文件的代码
麻烦指点一下

这是我的描述文件信息

这是我的描述文件信息
4%J(SK6VMVE`JN001URB1$J.png
密码修改失败请联系微信:mofangbao
发表于 2016-8-4 19:02:12 | 显示全部楼层
感谢楼主的分享,解决了问题
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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