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