爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8900|回复: 13

[讨论] matlab 用quiverm与grads画的风场图感觉有些不同,请问这是为什么?

[复制链接]

新浪微博达人勋

发表于 2014-11-20 15:03:25 | 显示全部楼层 |阅读模式

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

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

x

matlab画的图

matlab画的图

grads画的图

grads画的图

左图为matlab画的图,右图为grads画的图,虽然很多大值区两者差不多,但是很明显有很多地方两者的结果还是不一样,比如中国大部,请问这是什么原因?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2014-11-20 16:29:05 | 显示全部楼层
楼主能不能把matlab的图大一点,看不太清…把分辨率调高一点
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-11-20 16:34:53 | 显示全部楼层
这个还可以画图,受教了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-11-20 16:42:54 | 显示全部楼层
首先,reference arrow就不一样啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-11-20 16:45:16 | 显示全部楼层
楼主,我想问下matlab的右下角图例箭头怎么加上去的呢,能不能贴下代码分享下
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-11-20 16:53:53 | 显示全部楼层
二爷名声在外 发表于 2014-11-20 16:29
楼主能不能把matlab的图大一点,看不太清…把分辨率调高一点

点击图片可以看到大图
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-11-20 16:55:34 | 显示全部楼层
一水天下 发表于 2014-11-20 16:45
楼主,我想问下matlab的右下角图例箭头怎么加上去的呢,能不能贴下代码分享下

http://bbs.06climate.com/forum.php?mod=viewthread&tid=21699 我也是参考了该帖子的方法。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-11-20 16:57:04 | 显示全部楼层
空谷幽竹 发表于 2014-11-20 16:42
首先,reference arrow就不一样啊

请忽略我的arrow大小值,为什么会在有些地方方向就不同,而大部分区域又一致呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 成长值: 32430
发表于 2014-11-20 17:02:47 | 显示全部楼层
贴上程序?
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2014-11-20 17:12:10 | 显示全部楼层

alatg = [ 90:-2.5:-90 ];along = [ 0:2.5:360-2.5];
[axllon,axllat]=meshgrid(double(along),double(alatg));
load coast
figure('Color', 'w')
axesm('MapProjection', 'Eqdcylin', 'MapLonLimit', [60,150], 'MapLatLimit', [0,60])
framem on; gridm off; plabel on; mlabel('South')
setm(gca, 'MLineLocation', 10, 'PLineLocation', 10)
setm(gca, 'MLabelLocation', 10, 'PLabelLocation', 10)
tightmap; axis off; hold on
linem(lat, long, 'Color', 'k', 'LineWidth', 1.5)
hold on;
u=double(uu);v=double(vv);
Handle=quiverm(latitude',longitude',v, u,  2.0);
m_grid('xtick',7,'tickdir','none','ytick',7,'linest','none');
labelx=[140 141]; labely=[2 3]; % location of label
labelu=[10 0]; labelv=[0 0]; % set 25m/s as a standard
scale_auto=Get_Autoscale(longitude, latitude, v, u);
scale_label=Get_Autoscale(labelx, labely,  labelu, labelv);
scale_factor=scale_auto/scale_label;
Handle=quiverm(labely, labelx, labelv, labelu, 2.0*scale_factor); % <----- size: 2.0*scale_factor
set(Handle, 'Color', 'r', 'LineWidth', 2.5)
textm(labely(1)-1.5, labelx(1), [num2str(labelu(1)) ' m/s'],...
    'FontWeight', 'bold', 'Color', 'r', 'FontName', 'Helvetica', 'FontSize', 12)
% title
title('\fontsize{14}Wind quiver with label', 'FontWeight', 'bold', 'FontName', 'Helvetica')
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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