VC++开发GIS系统（283）判断折线与面的拓扑关系

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

关于线面之间的拓扑

通常情况下可以将线面之间拓扑关系做如下划分：

1. 折线与面相交
2. 折线与面相离
3. 折线在面的内部

线面之间的拓扑算法思路

在前文《VC++开发GIS系统（282）判断直线与面的拓扑关系》中已经判断了直线与面的拓扑关系。这里可以很容易有一个拓展：

1. 判断折线中每一条直线在面元素内部，则在面内。
2. 只要有一条直线与面相交，则相交。
3. 判断折线中每一条直线都在面元素外部，则相离。

算法实现

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

int  polyLineAndPolygon(vector<malaPoint>&Line, vector<malaPoint>&Polygon); //折线与多边形拓扑关系 相离返回1，相交返回2，在区内返回3

• @ Line 折线要素
• @ Polygon 面（区）要素

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

/*
* 折线与多边形拓扑关系 
*/
int malaLogic::polyLineAndPolygon(vector<malaPoint>&Line, vector<malaPoint>&Polygon)
{
  int length = Line.size();
  int index = 0;
  int m = 0;
  for (int i = 1; i < length; i++)
  {
    m = lineAndPolygon(Line[i - 1], Line[i], Polygon);
    if (m == 1)
      index++;
    if (m == 2)
      return 2;
  }
  if (index == length - 1)
    return 1;
  return 3;
}

其中函数 lineAndPolygon 用于判断直线与面的拓扑关系，可以参考文章《VC++开发GIS系统（282）判断直线与面的拓扑关系》

详细代码可以参考：GitHub