m_quiver画出的风场不对劲，求教大家

fghjnn

初来乍到，参照网上程序想试试画海面风场，结果陆地上的也有，想问问该怎么去掉陆地上的风矢量？
lat= ncread('vwnd.mon.mean.nc','lat');
lon= ncread('vwnd.mon.mean.nc','lon');
v= ncread('vwnd.mon.mean.nc','vwnd');
u1= ncread('uwnd.mon.mean.nc','uwnd');
m_proj('miller','lat',[0,20],'lon',[100,125]);
m_coast('patch',[.6 .6 .6],'edgecolor','none');
m_grid('box','fancy','linestyle','none');
hold on
[x,y]=meshgrid(lon,lat);
v1=v(:,:,200);
u2=u1(:,:,200);
u2(u2==-9999)=nan;
u3=u2.*cosd(y');
v1(v1==-9999)=nan;
d =1 ;dd =5;
m_quiver(x(1:d:end,1:d:end),y(1:d:end,1:d:end),u3(1:d:end,1:d:end)'./dd,v1(1:d:end,1:d:end)'./dd,0);

minqing123

先m_quiver，再m_coast就把陆地上的风盖住了

minqing123

考虑到陆地的不规则可能挡住风矢量，还有一种办法，找一份海温数据插值成和风矢量数据一样的网格，然后做判定，把海温数据缺省的格点对应的风矢量全部设为缺省。再画图就没问题了。

夏夜

使用地形数据把(etopo5/etopo3)大于的网格点去掉，

夏夜

使用地形数据(etopo5/etopo3)把大于零的网格点去掉。

fghjnn

minqing123 发表于 2014-12-19 13:57
先m_quiver，再m_coast就把陆地上的风盖住了

恩，可以了，谢谢！

fghjnn

夏夜 发表于 2014-12-19 14:08
使用地形数据(etopo5/etopo3)把大于零的网格点去掉。

谢谢，等我水平再高点再试试你的方法，非常感谢！！

故乡月(环)

学习学习，最近也在用MATLAB

∮プξ☆

先画quiver再贴地形，还是盖不住啊

fghjnn

∮プξ☆ 发表于 2015-3-22 13:26
先画quiver再贴地形，还是盖不住啊

我试了一下 print出来的话确实盖不住 请问解决了吗  或者试一下上面说的差值的方法