爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4256|回复: 6

[程序设计] MATLAB画风场出现矢量重叠

[复制链接]

新浪微博达人勋

发表于 2019-3-24 11:50:36 | 显示全部楼层 |阅读模式

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

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

x
我在用matlab画图的时候出现了以下的情况,风矢量的箭头都重叠在一块儿了,高层的风场直接就是一条线了,请问有什么解决的办法吗?
以下是我的程序:
for i=1:4
    subplot(2,2,i);
    uu=u(:,:,i)';vv=v(:,:,i)';
    quiver(lat,level,uu,vv,1.5);
    hold on;
    xlim([25 35]);
    ylim([500 1000]);

end
放上了我的uu和vv的部分数据,哪位大神能够给出一点思路呀!!!谢谢!

风场

风场

uu

uu

vv

vv
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2019-3-25 18:12:33 | 显示全部楼层
风是矢量,应该需要加开始和结束的坐标
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-3-25 20:43:06 | 显示全部楼层
东郭北辰 发表于 2019-3-25 18:12
风是矢量,应该需要加开始和结束的坐标

不太明白,能说清楚一点吗?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-4-26 10:23:05 | 显示全部楼层
可以这样试一下,因为matlab中quiver(x,y,u,v)表征的矢量箭头,画出来的箭头跟x,y,u,v都相关的。
lev = level(1:levid,:);
u0 = ua(1:levid,:);
v0 = va(1:levid,:);
[~,m] = size(hr); [n,~] = size(lev);
scale_x = 1*(max(max(hr))-min(min(hr)))/m;
scale_y = 1*(max(max(lev))-min(min(lev)))/n;
xx = (hr+1)/scale_x;
yy = lev/scale_y;
%     [x,y] = meshgrid((hr+1)/scale_x,lev/scale_y);
[x,y] = meshgrid(xx,yy);
  quiver(x,y,u0,v0,1);
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-4-26 13:09:02 | 显示全部楼层
鲸鱼便子 发表于 2019-4-26 10:23
可以这样试一下,因为matlab中quiver(x,y,u,v)表征的矢量箭头,画出来的箭头跟x,y,u,v都相关的。
lev = l ...

但是这样做存在一个问题,把坐标给换了。。。
画图的时候得在视框里一个一个调,而且也不好叠加其它的图,比如阴影图啥的,坐标变了。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-4-26 14:25:53 | 显示全部楼层
zzw9701 发表于 2019-4-26 13:09
但是这样做存在一个问题,把坐标给换了。。。
画图的时候得在视框里一个一个调,而且也不好叠加其它的图 ...

matlab画风场的垂直剖面确实不太方便,要么就上NCL或者grads
坐标可以自己用tick和ticklabel更改,叠加阴影图可以用和上面一样的比例转换。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-4-26 15:54:47 | 显示全部楼层
鲸鱼便子 发表于 2019-4-26 14:25
matlab画风场的垂直剖面确实不太方便,要么就上NCL或者grads
坐标可以自己用tick和ticklabel更改,叠加 ...

哦哦哦,好的,非常感谢!!!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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