python 判断点在随机多边形内_pnpoly算法——判断点是否在多边形内

pnpoly算法原理：从一个目标点引出一条射线，统计这条射线与对变形的交点个数。如果有奇数个交点，则说明目标点在多边形内，若为偶数个交点，则在外。

C语言实现：(参考

判断点是否在多边形内部的pnpoly算法_geerniya的博客-CSDN博客_pnpoly算法​blog.csdn.net

)

/**
* n：多边形的定点数目
* pos_x：点的x坐标
* pos_y：点的y坐标
* vert：二维数组，多边形顶点的坐标
* 返回1，点在多边形范围内
* 返回0，点不在多边形范围内
*/
int pnpoly(int n, double pos_x, double pos_y, double vert[][2])
{
    int i, j, c = 0;
	
	for (i = 0, j = n - 1; i < n; j = i++)
	{
	      if (((vert[i][1] > pos_y) != (vert[j][1] > pos_y)) &&
		   (pos_x < (vert[j][0] - vert[i][0]) * (pos_y - vert[i][1]) / (vert[j][1] - vert[i][1]) + vert[i][0]))
	      {
		  c = !c;
	      }
	}	
	return c;
}

Python实现：

def is_pt_in_poly(pt, poly):
    '''判断点是否在多边形内部的(pnpoly 算法)
    '''
    nvert = len(poly)
    vertx = []
    verty = []
    testx = pt[0]
    testy = pt[1]
    for item in poly:
        vertx.append(item[0])
        verty.append(item[1])

    j = nvert - 1
    res = False
    for i in range(nvert):
        if (verty[j] - verty[i]) == 0:
            j = i
            continue
        x = (vertx[j] - vertx[i]) * (testy - verty[i]) / (verty[j] - verty[i]) + vertx[i]
        if ((verty[i] > testy) != (verty[j] > testy)) and (testx < x):
            res = not res
        j = i
    return res