Matlab处理边界文件缓冲区的代码示例

在地理信息系统中(GIS)中,我们一般使用ArcGIS对点、线、面要素进行缓冲区分析。这里我提供一个可以对封闭曲线进行缓冲区操作的matlab代码,具体的函数内容如下:


clear;clc;
%% Load original boundary data
border = load('Tibet_boundary.txt');
% sta_info=importdata('sites.info');
%% Generate a buffer with any radiuses for study area
d=2.5; % 2.5 degree buffer
polyout = polybuffer([border(:,1) border(:,2)],'lines',d);
out = inpolygon(polyout.Vertices(:,1),polyout.Vertices(:,2),border(:,1),border(:,2));
edge_points=[polyout.Vertices(~out,1) polyout.Vertices(~out,2)];
d=0.25; % 0.25 degree buffer
polyout = polybuffer([border(:,1) border(:,2)],'lines',d);
out = inpolygon(polyout.Vertices(:,1),polyout.Vertices(:,2),border(:,1),border(:,2));
edge_points1=[polyout.Vertices(~out,1) polyout.Vertices(~out,2)];
%% Plot original and extended boundary
figure('color',[1 1 1])
plot(border(:,1),border(:,2));
hold on
plot(edge_points(:,1),edge_points(:,2));
hold on
plot(edge_points1(:,1),edge_points1(:,2));

得到了青藏高原边界2.5°和0.25°的缓冲区






免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删

QR Code
微信扫一扫,欢迎咨询~

联系我们
武汉格发信息技术有限公司
湖北省武汉市经开区科技园西路6号103孵化器
电话:155-2731-8020 座机:027-59821821
邮件:tanzw@gofarlic.com
Copyright © 2023 Gofarsoft Co.,Ltd. 保留所有权利
遇到许可问题?该如何解决!?
评估许可证实际采购量? 
不清楚软件许可证使用数据? 
收到软件厂商律师函!?  
想要少购买点许可证,节省费用? 
收到软件厂商侵权通告!?  
有正版license,但许可证不够用,需要新购? 
联系方式 155-2731-8020
预留信息,一起解决您的问题
* 姓名:
* 手机:

* 公司名称:

姓名不为空

手机不正确

公司不为空