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

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

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

发表评论