爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 3914|回复: 0

[讨论] 插值问题

[复制链接]
发表于 2016-4-8 12:42:52 | 显示全部楼层 |阅读模式

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

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

x
Ktemp温度数据在nc文件中是以altitude*event的形式存储的,但是在高度方向分辨率为三四百米,我想插值到100m。试了好多次还是错,求大神赐教!
%打开单个nc文件并读取数据
Ktemp = ncread('H:\毕设\毕设资料\数据\SABER\data\SABER_chh_Dec2012_v2.0_1458740929.nc','ktemp');

%处理异常数据
[M,N]=size(Ktemp);
for m=1:M
    for n=1:N
        %id=find( LonData>1000);%选出异常数据的项
        id=find(Ktemp>1000);%令异常数据为0
        Ktemp(id)=0;
    end
end

Ktemp=double(Ktemp);%将数据改为双精度类型

%几何高度插值
[xx,yy]=meshgrid(20:10:100,1:1:N);
[xi,yi]=meshgrid(20:0.1:100,1:1:N);
Ktemp=permute(Ktemp, [2 1]);
Ktempi =interp2(xx,yy,Ktemp,xi,yi);
mesh(Ktemp);

出错如图

V8F5JMDC@J80$01)NT7EDDY.png
8]6%@8RR`QTE[`9DP1J$2O3.png
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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