python 判断点在随机多边形内_Python确定散点是否在多边形内,python,判断,内部

def is_inpolygun(ploygun_point, scatter_point, border=True):

# 判断散点是否在多边形内部，border=True表示在边界上的也作为内部，否则作为外部

# 输入的多边形为二维数组(x, y)，多行两列，散点为单个

# 输出True为在内部, False在外部

# 1、确保输入的值为numpy类型

pp = np.array(ploygun_point)

sp = np.array(scatter_point)

# 2、如果多边形的头尾相连，去掉重复的一个点

if (pp[0] == pp[-1]).all():

pp = np.delete(pp, -1, axis=0)

# 3、 确定交点，

x = pp[:, 0]

y = pp[:, 1]

px = sp[0]

py = sp[1]

# 边上的各个点的Y坐标与散点的Y坐标相减，如果相交，则相邻的两个交点的差值必然一个>=0，另一个<=0

dy = y - py

dy_next = dy.copy()

dy_next[0] = dy[-1]

dy_next[1:] = dy[:-1]

# 两者相乘，其积<=0 则为相交点

dy_mu = dy * dy_next

negative_index = np.where(dy_mu <=0)

# 如果有交点，判断散点左边的交点个数，奇数在内,偶数在外,

if len(negative_index) >0 :

isp = pp[np.where(dy_mu <= 0)]

if not border:

# 如果有交点与原多边形的点重合，会出现2个一样的行，删除重复的行

isp= np.unique(isp, axis=0)

left_num = len(isp[np.where(isp[:, 0] <= px)])  # 计算左边点个数,包括本身

else:

left_num = len(isp[np.where(isp[:, 0] < px)])  # 计算左边点个数

# 计算奇偶数

if left_num % 2 ==0 :

return False

else:

return True

else:

return False

# 以下为调用测试

is_a = is_inpolygun(np.array([x, y]).T, pa)  # 输出False

is_b = is_inpolygun(np.array([x, y]).T, pb)  # 输出True

is_c = is_inpolygun(np.array([x, y]).T, pc)  # 边上的算成内部，输出True

is_d = is_inpolygun(np.array([x, y]).T,   pc, False)  # 边上的算成外部，输出False