Qgs开发“geos exception: topologyexception:”

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

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

中间涉及到的问题:

  1. 特征是QgsPolygon类型,两个多边形的关系判断(九交模型):相交/不相交,其中相交又包括touches、overlaps、crosses,若仅判断其中一个,则不全面,容易漏掉,因此判断两个特征是否不disjoint,只要不满足这个关系,则判定两特征相交。
  2. 做两个特征合并(combine)的时候,程序异常警告“geos exception: topologyexception:”,其中包括两种类型,一是“side location conflict at…”,二是“ Input geom 1 is invalid: Self-intersection at or near point…”,查阅很多资料和尝试后,发现是输入的特征是不可用的,造成的可能原因是:多边形拓扑问题。目前解决途径如下:
    1) 对前期提取的轮廓,计算凸包,再转换为QgsFeature类型
    /* 对每个轮廓计算其凸包*/
    vectorpoly(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”

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

你可能感兴趣的:(Qgs开发)