【walkscript】求共界线

walk.djg · 发表于 2012-10-31 13:53

在应用中，往往需要求几个地物共界线，使用walkscript可以快速批量完成这些任务。

如下图，需要求它们的共界线。
求共界线前.png

walkscript求得的结果如下图所示（红色的线）
walkscript求共界线后.png

walk.djg · 发表于 2012-10-31 13:54

//求共界线
void MutualBoundary(wkLayer layer, double eps, bool bSel)
{
int nFeaCount = bSel ? layer.getSelectionCount() : layer.getFeatureCount();

//求不同的边界线
wkClean cln;
int i = 0;
cln.create(layer.getBox());
for (i = 0; i < nFeaCount; ++i)
{
wkGeometry geom = bSel ? layer.getSelectionAt(i).getGeometry()
: layer.getFeatureAt(i).getGeometry();
cln.xorOverLine(geom, eps);
}
cln.cleanZero(eps);
cln.combineToArc();
array ps;
wkParts parts = cln.getFirstArc();
while (parts.handle())
{
ps.add(parts);
parts = cln.getNextArc();
}

//用所有地物边界减去不同的部分，得到共界的部分
for (i = 0; i < nFeaCount; ++i)
{
wkGeometry geom = bSel ? layer.getSelectionAt(i).getGeometry()
: layer.getFeatureAt(i).getGeometry();
cln.orOverLine(geom, eps);
}
for (i = 0; i < ps.getSize(); ++i)
{
cln.diffOverLine(ps[ i ], eps);
}
cln.cleanZero(eps);
int nCount = cln.combineToArc();
parts = cln.getFirstArc();
while (parts.handle())
{
if (parts.handle())
layer.addFeature(parts,false);
parts = cln.getNextArc();
}

cln.free();
}

walk.djg · 发表于 2012-10-31 13:54

测试代码如下：

wkGeoset geoset;
wkView view;
void main()
{
wkLayer layer = geoset.getEditableLayer();
double eps = 0.001;
MutualBoundary(layer, eps, false);
return;
}

复制代码

walkgis · 发表于 2012-10-31 14:00

多谢分享！

alert · 发表于 2012-10-31 14:28

多谢分享学习了！

lizuopei · 发表于 2012-10-31 16:31

感谢楼主啊

oxygen831 · 发表于 2012-10-31 20:46

感谢楼主分享，多多学习

walk.djg · 发表于 2012-11-1 08:02

//求共界线void MutualBoundary(wkLayer layer, double eps, bool bSel){ int nFeaCount = bSel ? layer.getSelectionCount() : layer.getFeatureCount(); //求不同的边界线 wkClean cln; int i = 0; cln.create(layer.getBox()); for (i = 0; i < nFeaCount; ++i) { wkGeometry geom = bSel ? layer.getSelectionAt(i).getGeometry() : layer.getFeatureAt(i).getGeometry(); cln.xorOverLine(geom, eps); } cln.cleanZero(eps); cln.combineToArc(); array ps; wkParts parts = cln.getFirstArc(); while (parts.handle()) { ps.add(parts); parts = cln.getNextArc(); } //用所有地物边界减去不同的部分，得到共界的部分 for (i = 0; i < nFeaCount; ++i) { wkGeometry geom = bSel ? layer.getSelectionAt(i).getGeometry() : layer.getFeatureAt(i).getGeometry(); cln.orOverLine(geom, eps); } for (i = 0; i < ps.getSize(); ++i) { cln.diffOverLine(ps[i], eps); } cln.cleanZero(eps); int nCount = cln.combineToArc(); parts = cln.getFirstArc(); while (parts.handle()) { if (parts.handle()) layer.addFeature(parts,false); parts = cln.getNextArc(); } cln.free(); }<div></div>

复制代码

walk.djg · 发表于 2012-11-1 08:02

wkGeoset geoset;
wkView view;
void main()
{
wkLayer layer = geoset.getEditableLayer();
double eps = 0.001;
MutualBoundary(layer, eps, false);
return;
}

复制代码

Lock · 发表于 2020-4-28 12:19

楼主大佬，您有walkscript的学习资料吗？

【walkscript】求共界线

爱心勋章

组织勋章

地信元老

灌水勋章

荣誉会员勋章

活跃勋章

地信专家组

VIP勋章

贡献勋章

成就学员勋章