英语原文:https://www.ecse.rpi.edu/~wrf/Research/Short_Notes/pnpoly.html
基本思想是,以该点向任意方向(这里是x轴正向)作射线,若与多边形的交点个数为偶数,则在多边形之外,否则在多边形之内。
这是作者给出的源代码:
int pnpoly(int nvert, float *vertx, float *verty, float testx, float testy)
{
int i, j, c = 0;
for (i = 0, j = nvert-1; i < nvert; j = i++) {
if ( ((verty[i]>testy) != (verty[j]>testy)) &&
(testx < (vertx[j]-vertx[i]) * (testy-verty[i]) / (verty[j]-verty[i]) + vertx[i]) )
c = !c;
}
return c;
}
初看还比较难理解。
nvert表示多边形点的个数,*vertx表示多边形点的x轴数组,*verty表示多边形点的y轴数组,testx和testy分别表示测试点的x、y轴。
对于i=0, j = nvert-1; i<nvert; j=i++的解释
这句话的意思是遍历一遍多边形的两个相邻点j, i。j在i之前。当i为0时,j为i的前一个点,即nvert-1。
对于(verty[i]>testy) != (verty[j]>testy)的解释
这句话表示点i和j必须在测试点y值的两侧。这里处理非常巧妙,而且防止了测试点、多边形相邻点在同一直线上,防止了后半部分条件中除数为0
对于testx < (vertx[j]-vertx[i]) * (testy-verty[i]) / (verty[j]-verty[i]) + vertx[i]的解释
这句话利用了直线的两点式。