本文介绍GIS系统中,点与面之间拓扑关系的判断方法。
关于点面之间的拓扑
通常情况下可以将点与面的关系做如下划分:
点在面内
点在面的边界上
点在面外
本文做了一个简化处理,将1和2两种情况当做一种。所以本文中判断的点面之间的拓扑关系划分如下:
点在面内
点在面外
点面之间的拓扑算法思路
算法如图:
算法实现
在 _malaBase.h 中声明如下的函数:
bool isPointInPolygon(malaPoint point, vector& Poly);//点是否在多边形内部
@ point 点要素
@ Poly 面(区)要素
并在 _malaBase.cpp中实现上述函数,代码如下:
/*
* 判断点是否在区上
*/
bool malaLogic::isPointInPolygon(malaPoint point, vector& 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.yyPrev)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