xv= [0 3 3 0 0]; %x坐标
yv= [0 0 3 3 0];%y坐标
x=4;
y=4;
in=inpolygon(x,y,xv,yv)
plot(xv,yv,x(in),y(in),‘.r‘,x(~in),y(~in),‘.b‘)
in=0;
===================from 气象家园
http://bbs.06climate.com/forum.php?mod=viewthread&tid=12340
[lon lat]=meshgrid([97:0.1:107],[21:0.1:30]);
% Z=griddata(gy_locat(:,2),gy_locat(:,3),EOF_used(:,1),lon,lat,‘v4‘);
yunnan=shaperead(‘yunnan.shp‘);
isin=inpolygon(lon,lat,yunnan.Lon,yunnan.Lat);
Z(~isin)=NaN;
contourf(lon,lat,Z,30);
shading flat
colorbar
hold on
plot(yunnan.X,yunnan.Y,‘-k‘,‘linewidth‘,3)
hold off
=============
http://bbs.06climate.com/forum.php?mod=viewthread&tid=13304
来论坛有段时间了,学了不少东西,今天也分享一个自己的经验。
今天看到一个帖子