区域闯入问题——如何判断一个点是否在闭合不规则多边形平面区域内(射线法)
一、问题概述
在目标检测与识别领域中,我们经常会遇到判断人员、车辆等物体的位置问题。例如,在城市交通中有一段道路因为发生塌陷需要施工,现来往车辆均不可以靠近这段危险道路。l利用搭建交通摄像头和语音报警器的方式,需要24小时不间断识别来往车辆的位置以及是否靠近此区域。
对于这个问题,主要分为两个步骤:
- 摄像头区域中车辆位置的精准定位问题;(Yolov5等目标检测算法能够取得不错的结果。关于Yolov5模型的快速搭建,见我另一篇文章)
- 得到①反馈的框定坐标后,我们得到可以推理得到框四个角的坐标,并判断它们的重要性(对于车来说4个角重要性平等,而对于行人来说,左下角和右下角坐标视为行人的足部位置,这两点才真实反馈了行人的位置),利用区域检测算法来确定是否闯入的问题。
对于识别目标,我们可以得到这些框定信息:
info[[[a1,b1,c1,d1],e1],[[a2,b2,c2,d2],e2]],a是左上角x坐标,b是左上角y坐标,c是右下角x坐标,d是右下角y坐标,e是置信度。
二、射线法
对于判断某个点的坐标是不是在区域范围(闭合不规则多边形)内的问题,射线法使用一种非常巧妙的思路予以解决。
假若有一坐标点P(x,y),要判斯它是否在多边形区域内,可从该坐标点向左/向右画一条水平扫描线(即射线)。然后计算此线段与区域边界的相交次数。如果次数为奇数,认为该点在多边形内;若为偶数,则该点在多边形外。
三、部分代码展示
#区域检测
def check(ver_1,ver_2,coordinate_1,coordinate_2):
#形成多边形的边界列表:1->2->3->4->5...n->1
length_1 = len(coordinate_1)
length_2 = len(coordinate_2)
frame_1 = set_edge(coordinate_1,length_1)
frame_2 = set_edge(coordinate_2,length_2)
level = 0
#先判断输入的两个坐标的y坐标是否符合常规
if ver_1[1] != ver_2[1]:
print('error!')
'''
射线法实现:偶数交点为区域外,奇数交点为区域内
逻辑判断上,禁止区>警示区>安全区
有一个点在禁止区就报禁止;否则有一个点在警示区就报警示;否则安全
'''
if isPoiInsidePoly(list(ver_1),frame_2) or isPoiInsidePoly(list(ver_2),frame_2):
#print('人员所在区域:禁止区!')
level = 2
return level
elif isPoiInsidePoly(list(ver_1),frame_1) or isPoiInsidePoly(list(ver_2),frame_1):
#print('人员所在区域:警示区!')
level = 1
return level
else:
#print('人员所在区域:安全区!')
level = 0
return level
if __name__ == '__main__':
#警示区ROI1(注意,区域坐标顺序应当满足顺时针/逆时针顺序)
ROI1 = [(1100,540),(900,540),(900,640),(1100,640)]
#危险区ROI2(注意,区域坐标顺序应当满足顺时针/逆时针顺序)
ROI2 = [(258, 60), (346, 247), (366, 478), (140, 184),(166,293)]
# info[[[a1,b1,c1,d1],e1],[[a2,b2,c2,d2],e2]],a是左上角x坐标,b是左上角y坐标,c是右下角x坐标,d是右下角y坐标,e是置信度
info = [[[401.0, 341.0, 1188.0, 653.0], 0.53], [[406.0, 326.0, 765.0, 637.0], 0.81]]
"""
在此处根据人员的位置信息info和检测区域ROI比对
判断异常等级 level
"""
level = 0
for i in range(len(info)):
f_left = (info[i][0][0], info[i][0][3])
f_right = (info[i][0][2], info[i][0][3])
level = check(f_left,f_right,ROI1,ROI2)
if level >= 1:
break
详细代码:https://download.csdn.net/download/qq_42972996/85327145