- 积分
- 1254
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-1-30
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
这样写有什么问题吗??想把2.5*2.5的风场插值到2*2,画出来的图很奇怪
lat3=ncread('d:/pre/kaiti/data2/uwnd.mon.mean.nc','lat');
lon3=ncread('d:/pre/kaiti/data2/uwnd.mon.mean.nc','lon');
u3=ncread('d:/pre/kaiti/data2/uwnd.mon.mean.nc','uwnd');
uip=u3(:,:,3,:);uip=squeeze(uip);
uip=uip(:,:,25:624); %1950.1-1999.12
uip=double(uip); uip(abs(uip)>900)=NaN;
[m2,n2,tt]=size(uip);
clear u3;
u4(:,:,:,:)=reshape(uip,m2,n2,12,50);
lat3=ncread('d:/pre/kaiti/data2/vwnd.mon.mean.nc','lat');
lon3=ncread('d:/pre/kaiti/data2/vwnd.mon.mean.nc','lon');
v3=ncread('d:/pre/kaiti/data2/vwnd.mon.mean.nc','vwnd');
vip=v3(:,:,3,:);vip=squeeze(vip);
vip=vip(:,:,25:624);
vip=double(vip); vip(abs(vip)>900)=NaN;
[m2,n2,tt]=size(vip);
clear v3;
v4(:,:,:,:)=reshape(vip,m2,n2,12,50);
uo=mean(u4,4);uo1=(uo(:,:,6)+uo(:,:,7)+uo(:,:,8))/3;uo1=squeeze(uo1);
vo=mean(u4,4);vo1=(vo(:,:,6)+vo(:,:,7)+vo(:,:,8))/3;vo1=squeeze(vo1);
[xx,yy]=meshgrid(lat3,lon3);
[xi,yi]=meshgrid(-90:2:90,0:2:358);
ucha1=interp2(xx,yy,uo1,xi,yi);
vcha1=interp2(xx,yy,vo1,xi,yi);
|
-
|