- 积分
- 348
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-5-7
- 最后登录
- 1970-1-1
|
楼主 |
发表于 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') |
|