Qgs开发“geos exception: topologyexception:”

最新推荐文章于 2024-05-31 18:36:38 发布

夏安An

最新推荐文章于 2024-05-31 18:36:38 发布

阅读量2.3k

点赞数

分类专栏： Qgs开发文章标签： Qgs二次开发

本文链接：https://blog.csdn.net/qq_37974957/article/details/89917692

版权

1 篇文章 0 订阅

订阅专栏

   最近在做Qgs开发的时候，遇到一些问题，各种查资料，百度，最终总算得到了初步解决，不知道自己理解的是否正确，但目前来说，没看到什么问题，在此写文记录，并在之后的学习中，继续完善和更正。

任务：当向QgsVectorLayer图层中添加特征（QgsFeature）时，需要判断是否和之前图层中已有的特征有相交，若有，则两个特征进行合并，并删除之前的特征，添加重新合并后的新特征。

中间涉及到的问题：

特征是QgsPolygon类型，两个多边形的关系判断（九交模型）：相交/不相交，其中相交又包括touches、overlaps、crosses，若仅判断其中一个，则不全面，容易漏掉，因此判断两个特征是否不disjoint，只要不满足这个关系，则判定两特征相交。
做两个特征合并（combine）的时候，程序异常警告“geos exception: topologyexception:”，其中包括两种类型，一是“side location conflict at…”,二是“ Input geom 1 is invalid: Self-intersection at or near point…”,查阅很多资料和尝试后，发现是输入的特征是不可用的，造成的可能原因是：多边形拓扑问题。目前解决途径如下：
1）对前期提取的轮廓，计算凸包，再转换为QgsFeature类型
/* 对每个轮廓计算其凸包*/
vector<vector >poly(g_vContours.size());
for (int i = 0; i < g_vContours.size(); i++)
{
approxPolyDP(Mat(g_vContours[i]), poly[i], 5,true);
}
2）对QgsFeature得到的QgsGeometry进行simplify
QgsGeometry* comGeo = newPolygonGeom->simplify(10);
qDebug()<<“aftersimplify validate is”<isGeosValid();

通过这两种方法处理得到的QgsGeometry验证其有效性是valid。**

关注