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