检测二维平面内一个点在多边形内的方法

最新推荐文章于 2024-04-30 20:30:00 发布

yessharing

最新推荐文章于 2024-04-30 20:30:00 发布

阅读量3.3k

点赞数

分类专栏：杂谈

杂谈专栏收录该内容

4 篇文章 0 订阅

订阅专栏

最近项目需要, 所以查了一下实现2D平面内判断一个点在多边形内的方法. 感觉这个方法非常巧妙且计算量不大.

在这个方法中, 多边形边界是用点来描述的. 多边形的凸凹性不受限制, 也就是说, 多边形可以是凸也可以是凹的.

对于任意两个构成多边形边界的点, 有如下这种情况:(这里假设这两个点连线斜率为正于是画成了下面的样子)

这里写图片描述

构成多边形的点是 $x_i, y_i)$ 和 $x_{i+1}, y_{i+1})$ , 要判断的点是 $x_0, y_0)$ . 红色的线是一条从 $x_0, y_0)$ 出发斜率为 $\infty$ 的射线.

首先判断 $x_0, y_0)$ 是否在 $x_i$ 和 $x_{i+1}$ 之间. 故该点需要满足
$x_i<x_0<x_{i+1}或者是\\ x_{i+1}<x_0<x_{i}$

接下来还需要判断 $x_0, y_0)$ 是不是在黑色的线的下面
$k=\frac{y_{i+1}-y_i}{x_{i+1} - x_i}\\ y_0<k(x_0-x_i)+y_i$
这样就完成了判断点 $x_0, y_0)$ 在黑色直线下的过程.

接下来对于构成多边形边界的所有的点 $x_i$ 和 $x_{i+1}$ 都应用上面的方法进行判断. 只要满足上述准则, 就设置计数器加1.
如果要判断的点 $x_0, y_0)$ 在多边形的范围内, 那一定满足计数器为奇数, 也就是对2取余为1.
例如, 下图计数器结果是1.

这里写图片描述

下图中结果是3.
这里写图片描述

代码如下

bool contains(double x0, double y0)
{
	int crossings = 0;
	for(int i=0; i<N; i++)
	{
		double slope = (y[i+1] - y[i])/(x[i+1] - x[i]);
		bool cond1 = (x[i] < x0) && (x[i+1] > x0);
		bool cond2 = (x[i+1] < x0) && (x[i] > x0);
		bool above = (y0 < slope * (x0 - x[i]) + y[i]);
		if((cond1 || cond2) && above)
		{
			crossings++;
		}
	}
	return (crossings %2 != 0)
}