麻辣GIS微信平台

更多 GIS 干货

微信关注不错过

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

相关阅读

麻辣GIS-Sailor

作者:

GIS爱好者,学GIS,更爱玩GIS。

声明

1.本文所分享的所有需要用户下载使用的内容(包括但不限于软件、数据、图片)来自于网络或者麻辣GIS粉丝自行分享,版权归该下载资源的合法拥有者所有,如有侵权请第一时间联系本站删除。

2.下载内容仅限个人学习使用,请切勿用作商用等其他用途,否则后果自负。

手机阅读
公众号关注
知识星球
手机阅读
麻辣GIS微信公众号关注
最新GIS干货
关注麻辣GIS知识星球
私享圈子

留言板(小编看到第一时间回复)