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