区域闯入问题——如何判断点是否在闭合不规则多边形平面内（射线法）

区域闯入问题——如何判断一个点是否在闭合不规则多边形平面区域内（射线法）

---

---

一、问题概述

在目标检测与识别领域中，我们经常会遇到判断人员、车辆等物体的位置问题。例如，在城市交通中有一段道路因为发生塌陷需要施工，现来往车辆均不可以靠近这段危险道路。l利用搭建交通摄像头和语音报警器的方式，需要24小时不间断识别来往车辆的位置以及是否靠近此区域。
对于这个问题，主要分为两个步骤：

1. 摄像头区域中车辆位置的精准定位问题；（Yolov5等目标检测算法能够取得不错的结果。关于Yolov5模型的快速搭建，见我另一篇文章）
   
2. 得到①反馈的框定坐标后，我们得到可以推理得到框四个角的坐标，并判断它们的重要性（对于车来说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