VC++开发GIS系统(284)判断面元素之间的拓扑关系

本文介绍GIS系统中面元素之间拓扑关系的判断方法。

关于面元素之间的拓扑

通常情况下可以将面元素之间拓扑关系做如下划分:

  1. 两个面完全重合
  2. 两个面相交
  3. 两个面相离
  4. 一个面内含另外一个面

其中为了简化,可以将1与4看作同一种情况。

面元素之间的拓扑算法思路

  1. 将一个面拆分成一条折线,并判断其是否与另一面相交,有一条相交则两个面相交。
  2. 如果一个面所有的点在另外一个面内,则包含
  3. 非上述情况则相离

算法实现

在 _malaBase.h 中声明如下的函数:

int  polygonAndPolygon(vector<malaPoint>&Polygon1, vector<malaPoint>&Polygon2);//面与面的拓扑关系 相交返回2,内含返回3,相离返回4
  • @ Polygon1 面要素1
  • @ Polygon2 面要素2

并在 _malaBase.cpp中实现上述函数,代码如下:

/*
* 面与面的拓扑关系  相等返回1,相交返回2,内含返回3,相离返回4,覆盖返回5
*/
int malaLogic::polygonAndPolygon(vector<malaPoint>&Polygon1, vector<malaPoint>&Polygon2)
{
  if (polyLineAndPolygon(Polygon1, Polygon2) == 2)
    return 2;
  int length = Polygon1.size();

  int index = 0;
  int i;
  for (i = 0; i < length; i++)
  {
    if (isPointInPolygon(Polygon1[i], Polygon2))
      index++;
  }
  if (index == length)
    return 3;
  int length2 = Polygon2.size();
  index = 0;

  for (i = 0; i < length2; i++)
  {
    if (isPointInPolygon(Polygon2[i], Polygon1))
      index++;
  }
  if (index == length2)
    return 5;
  return 4;
}

其中函数 polyLineAndPolygon 用于判断折线与面的拓扑关系,可以参考文章《VC++开发GIS系统(283)判断折线与面的拓扑关系isPointInPolygon 用于判断点与面的拓扑关系,可以参考文章《VC++开发GIS系统(280)判断点与面的拓扑关系

详细代码可以参考:GitHub

VC++开发GIS系统(283)判断折线与面的拓扑关系 VC++开发GIS系统(285)判断拓扑关系类的设计及实现

作者:,GIS爱好者。
分享本文,请您带上本文链接
分享到:

发表评论