////////////////////////////////////////////////////////////////
//计算点是否在多边形内
int msIntersectPointPolygon(pointObj *point, shapeObj *poly)
{//主函数1-在 0-不在
int i;
int status=MS_FALSE;
for(i=0; i<poly->numlines; i++)
{
if(msPointInPolygon(point, &poly->line[i]) == MS_TRUE) /* ok, the point is in a line */
status = !status;
}
return(status);
}
int msPointInPolygon(pointObj *p, lineObj *c)
{//算法非常奇怪!!!
int i, j, status = MS_FALSE;
for (i = 0, j = c->numpoints-1; i < c->numpoints; j = i++)
{ if ((((c->point[i].y<=p->y) && (p->y<c->point[j].y)) || ((c->point[j].y<=p->y) && (p->y<c->point[i].y))) && (p->x < (c->point[j].x - c->point[i].x) * (p->y - c->point[i].y) /(c->point[j].y - c->point[i].y) + c->point[i].x))
status = !status;
}
return status;
}
本文介绍了一种用于判断点是否位于多边形内部的算法实现。该算法通过扫描线方法,针对多边形的每条边检查点的位置,并依据交叉次数确定点是否在多边形内部。适用于地图信息系统(MGIS)等场景。
43

被折叠的 条评论
为什么被折叠?



