matlab中地图边界与掩膜(去掉边界外区域)的实现（基于shape文件）

tingling

学习学习~~

enyoyu

好东西必须来

enyoyu

地图有些省份是重复的，就plot不了。需要用不重复的地图

虫子要爱学习

非常感谢！

碧落kami

[lon lat]=meshgrid([73:0.125:105],[26:0.125:40]);
% Z=griddata(gy_locat(:,2),gy_locat(:,3),EOF_used(:,1),lon,lat,'v4');
tibetshp=shaperead('E:\Data\Tibetshp\qingzang.shp');
% tibetshp.Lon=[tibetshp( : ).X];
% tibetshp.Lat=[tibetshp( : ).Y];
tibet_X=[tibetshp(:).X];
tibet_Y=[tibetshp(:).Y];
isin=inpolygon(lon,lat,tibet_X,tibet_Y);
Z(~isin)=NaN;

contourf(lon,lat,Z,30);
shading flat
colorbar

hold on
plot(tibetshp.X,tibetshp.Y,'-k','linewidth',3)
hold off

大神帮忙看看 哪里有问题呢
isinpolygon
Error using contourf (line 61)
Z must be at least a 2x2 matrix.

Error in isinpolygon (line 11)
contourf(lon,lat,Z,30);

斥鷃

碧落kami 发表于 2016-4-11 10:27
[lon lat]=meshgrid([73:0.125:105],[26:0.125:40]);
% Z=griddata(gy_locat(:,2),gy_locat(:,3),EOF_used ...

报错说contour画图出问题了，你看看Z的维数情况吧

碧落kami

斥鷃 发表于 2016-4-11 12:33
报错说contour画图出问题了，你看看Z的维数情况吧

恩恩 解决了~

kyoslsl

{:eb502:}

白白melody

谢谢楼主的分享！！

xyan88

错误使用 contourf (line 61)
Z 必须至少是 2x2 矩阵。