请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3684|回复: 3

[程序设计] matlab风场插值画图

[复制链接]

新浪微博达人勋

发表于 2019-6-22 10:33:29 | 显示全部楼层 |阅读模式

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

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

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);

QQ截图20190622103311.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2019-6-22 16:25:50 | 显示全部楼层
解决了。。。
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2019-6-23 14:57:13 | 显示全部楼层
怎么解决的楼主?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-7-3 16:45:01 | 显示全部楼层

应该是vo=mean(v4,4);打错了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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