博文
利用matlab画风场，m_quiver
|||
用matlab画风场
重点1：
m_quiver(x(1:d:end,1:d:end),y(1:d:end,1:d:end),ustr_1_2(1:d:end,1:d:end)'./dd,vstr_2_1(1:d:end,1:d:end)'./dd,0);
m_quiver(30,-80,500./dd.*cosd(-80),0,0)
为了使得风场和标尺的大小一致，在m_quiver中最后一个值设置为“0”，即令scale=0
quiver(x,y,u,v,scale);
重点2：
d = 3;dd =10;
m_quiver(x(1:d:end,1:d:end),y(1:d:end,1:d:end),ustr_1_2(1:d:end,1:d:end)'./dd,vstr_2_1(1:d:end,1:d:end)'./dd,0);
通过d来调整箭头间距，dd来调整箭头长度（数值大小）
重点3：
ustr_1_2=ustr_1_1.*cosd(y');
h=m_quiver(30,-80,500./dd.*cosd(-80),0,0);
m_map 中用的是球坐标系，所以要在所有的纬向量上面乘以cosd（纬度） 注意：是cosd（纬度） cos（弧度）
其中纬度是每个纬向量所对应的纬度
在ustr_1_2=ustr_1_1.*cosd(y');里面，乘以cosd（y'） 要仔细检查，可以一步一步对应着去找。
%*********************************************************************
clear all;close all;clc
% 利用NCEP_Real-time_Marine_Data的风应力数据
%（http://www.esrl.noaa.gov/psd/data/gridded/data.nmc.marine.html）
% 绘制全球风场分布
%地形
m_proj('miller','lat',[-90,90],'lon',[0,360]);
m_coast('patch',[.6 .6 .6],'edgecolor','none');
m_grid('box','fancy','linestyle','none');
hold on
%NCEP mean风场数据
datapath_1='E:dataNCEPNCEP_Real-time_Marine_DataU-wind_Stressustr.mean.nc';
ncdata_1=netcdf.open(datapath_1,'NC_NOWRITE');
lat_1=netcdf.getVar(ncdata_1,0);
lon_1=netcdf.getVar(ncdata_1,1);
time_1=netcdf.getVar(ncdata_1,2);
ustr_1=netcdf.getVar(ncdata_1,3);
ustr_1_1=ustr_1(:,:,200);
datapath_2='E:dataNCEPNCEP_Real-time_Marine_DataV-wind_Stressvstr.mean.nc';
ncdata_2=netcdf.open(datapath_2,'NC_NOWRITE');
lat_2=netcdf.getVar(ncdata_2,0);
lon_2=netcdf.getVar(ncdata_2,1);
time_2=netcdf.getVar(ncdata_2,2);
vstr_2=netcdf.getVar(ncdata_2,3);
vstr_2_1=vstr_2(:,:,200);
[x,y]=meshgrid(lon_1,lat_1);
ustr_1_1(ustr_1_1<-1e+6)=nan;
ustr_1_2=ustr_1_1.*cosd(y');
vstr_2_1(vstr_2_1<-1e+6)=nan;
d = 3;dd =10;
m_quiver(x(1:d:end,1:d:end),y(1:d:end,1:d:end),ustr_1_2(1:d:end,1:d:end)'./dd,vstr_2_1(1:d:end,1:d:end)'./dd,0);
hold on
h=m_quiver(30,-80,500./dd.*cosd(-80),0,0);
set(h,'color','r','linewidth',1);
hold on
m_text(30,-75,'50m^2/s^2','fontsize',12);
print(gcf,'-djpeg','-r300','E:quiver_ncep');
%*********************************************************************
2012.11.07修改
% quiver问题：在不同的图中，令箭头大小一致。于是做了以下的实验：
%*********************************************************************
%% quiver test
%%
clear;clc;close all
figure
axes('position',[0.2,0.7,0.4,0.2]);
quiver(0,0,10,0,0); hold on
quiver(0,-0.5,20,0,0);
set(gca,'xlim',[0,25],'ylim',[-1.5 1.5])
axes('position',[0.2,0.4,0.4,0.2]);
quiver(0,0,20,0,0);
set(gca,'xlim',[0,25],'ylim',[-1.5 1.5])
axes('position',[0.2,0.1,0.6,0.2]);
quiver(0,0,20,0,0);
set(gca,'xlim',[0,25],'ylim',[-1.5 1.5])
print(gcf,'-djpeg','-r300','E:quiver_test');
%*********************************************************************
https://blog.sciencenet.cn/blog-419857-624284.html
下一篇：matlab新建文件夹及copy文件
