- 积分
- 164
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-3-25
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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);
出错如图
|
|