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

sxtylh

还没有亲自尝试，不过边界问题确实是个大问题，谢谢分享

Lighting

sxtylh 发表于 2016-9-18 18:04
还没有亲自尝试，不过边界问题确实是个大问题，谢谢分享

已经的到了完美解决，可以看一下这个帖子MATLAB对地图进行白化

sxtylh

Lighting 发表于 2016-9-18 20:53
已经的到了完美解决，可以看一下这个帖子MATLAB对地图进行白化

谢谢支持，字数限制，感谢的话说两遍，谢谢支持

/kfEspresso

这帖子真的很棒，学习了！

pbchao

感谢楼主分享！我要好好学习！

lingaochang

C:\Users\Administrator\Desktop\DQ%5DSATJB`PKR@KDH(%_7JCI.png您好，能帮我看看我的代码为什么边界没效果


%加载数据文件
load data_fj.txt;
%读书数据赋值变量
[a,b] = size(data_fj);
for i = 1:a
n(i) = data_fj(i,1);
x(i) = data_fj(i,2);
y(i) = data_fj(i,3);
rain(i) = data_fj(i,4);
end
%产生重构函数用的自变量
XI = linspace(min(x),max(x),838);
YI = linspace(min(y),max(y),838);

%产生重构函数用的自变量
[lon,lat] =meshgrid(XI,YI);
%数据格网化，对离散数据进行插值；对于一维、二维、三维规则数据点使用函数 interpl()/intep2 () /interpe3 ()进行插值，

Z = griddata(x,y,rain,lon,lat,'v4');
%
%[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('Fujian_province.shp');
R=makerefmat('RasterSize',size(Z),'Lonlim',[115 121],'Latlim',[23 30]);
MASK=vec2mtx(yunnan.Y,yunnan.X,Z,R);
MASK(find(MASK>1))=nan;
MASK(find(MASK==1))=0;

contourf(lon,lat,Z+MASK,20);
shading flat
colorbar
hold on
plot(yunnan.X,yunnan.Y,'-k','linewidth',3)
hold off

贫道敬孔

多谢楼主分享

斥鷃

lingaochang 发表于 2016-10-27 11:31
您好，能帮我看看我的代码为什么边界没效果

图没传上来，如果是全白的话查查经纬度位置

祁闯

多谢楼主分享，急需这样的

greenfirst

好东东，顶顶顶，赞