VC++开发GIS系统(280)判断点与面的拓扑关系

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

关于点面之间的拓扑

通常情况下可以将点与面的关系做如下划分:

  1. 点在面内
  2. 点在面的边界上
  3. 点在面外

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

  1. 点在面内
  2. 点在面外

点面之间的拓扑算法思路

算法如图:

算法实现

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

bool isPointInPolygon(malaPoint point, vector<malaPoint>& Poly);//点是否在多边形内部

@ point 点要素
@ Poly 面(区)要素

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

/*
* 判断点是否在区上
*/
bool malaLogic::isPointInPolygon(malaPoint point, vector<malaPoint>& Poly)
{
  int LineNum = Poly.size();
  malaPoint leftP = point;
  malaPoint rightP;
  rightP.x = getMaxX(Poly) + 1;
  rightP.y = point.y;
  int count = 0, yPrev = Poly[LineNum - 2].y;
  malaPoint v1, v2;
  v1 = Poly[LineNum - 1];
  for (int i = 0; i < LineNum; i++)
  {
    v2 = Poly[i];

    if (isPointInLine(leftP, v1, v2))
      return true;

    if (v1.y != v2.y)
    {
      if (isLineIntersect(v1, v2, leftP, rightP))
      {
        if (isPointInLine(v1, leftP, rightP))
        {
          if (v1.y<v2.y) { if (v1.y>yPrev)count++; }
          else { if (v1.y < yPrev) count++; }
        }
        else if (!isPointInLine(v2, leftP, rightP))
        {
          count++;
        }
      }
    }
    yPrev = v1.y;

    v1 = v2;
  }
  return   (count % 2 == 1);
}

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

详细代码可以参考:GitHub

VC++开发GIS系统(279)判断点与线之间的拓扑关系 VC++开发GIS系统(281)判断线与线的拓扑关系

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

发表评论