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

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

关于线要素之间的拓扑

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

1. 两条线相等
2. 两条线相交
3. 两条线不相交

本文做了一个简化处理，将1和2两种情况当做一种。所以本文中判断线要素之间的拓扑关系划分如下：

1. 两条线相交
2. 两条线不相交

线要素之间的拓扑算法思路

算法如图：

分别遍历两条折线中的每一条直线是否相交，只要有一条相交，则返回true，否则返回false。

算法实现

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

int  lineAndLine(vector<malaPoint>& Line1, vector<malaPoint>& Line2);  //折线与折线的拓扑关系 相交返回2;相离返回3

• @ point 点要素
• @ Poly 面（区）要素

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

/*
* 折线与折线的拓扑关系
*/
int malaLogic::lineAndLine(vector<malaPoint>& Line1, vector<malaPoint>& Line2)
{
  int legth1 = Line1.size();
  int length2 = Line2.size();

  for (int i = 1; i < legth1; i++)
  {
    for (int j = 1; j < length2; j++)
      if (isLineIntersect(Line1[i - 1], Line1[i], Line2[j - 1], Line2[j]))
      {
        return 2;
      }
  }
  return 3;
}

其中函数 isLineIntersect 用于判断两条直线是否相交，可以参考文章《VC++开发GIS系统（167）选择线之判断直线是否与折线相交算法》

详细代码可以参考：GitHub