一种高空抛掷物判别系统的制作方法

[0001]
本发明涉及高空抛物检测技术领域，尤其涉及一种高空抛掷物判别系统。


背景技术：

[0002]
高空抛物被称为“悬在城市上空的痛”，一直以来高空抛物行为备受关注，作为城市不文明行为的同时，它所带来的社会危害也很大。由于其实施场所多为高楼层，少有目击者，抛物时间短，甚至有人在抛物时故意隐去身影，使得相关部门很难追究抛物者的法律责任，此类事件屡见不鲜、屡禁不止。因此对高空抛物行为进行识别和取证，已经成为当前城市管理和物业管理的迫切需要。
[0003]
抛物行为的识别和取证，可以通过对抛物目标轨迹检测的方式进行，当前对抛物目标运动轨迹的检测，主要是通过运动目标检测和轨迹拟合的方式实现的，这种方法首先通过经典的运动目标检测方法，例如背景建模或帧差法，实现连续帧中运动目标的提取，然后通过每帧图像中运动目标的坐标，拟合出连续帧中运动目标的轨迹，并与已知抛物轨迹或抛物线公式计算得到抛物轨迹进行比对，最终确定目标运动轨迹是否为抛物运动轨迹。在实际应用中，上述方式存在一定问题，由于高层建筑物屋面受大气活动的影响，屋面存在湍流，且抛物目标种类繁杂，气动外形千差万别，因此高空抛物目标的运动轨迹非常复杂，难以用一条理论公式来进行描述，而采取与已知高空抛物轨迹样本进行比对，则存在轨迹样本难以获取的问题，同时，进行轨迹比对时需要对整个样本库进行遍历，计算复杂度高，系统复杂；而采取实验的模拟的方式，由于实验场地限制，实验费用开销较大，且无法模拟实际场景中多变的风场条件，难以获取到不同气动外形的抛掷物在不同风场下的运动轨迹，这造成模型在训练过程中缺乏有效的训练样本数据，直接影响模型的判别性能。


技术实现要素：

[0004]
本发明的目的是提供一种高空抛掷物判别系统，将高空抛物目标运动检测与识别转化为一个序列分类问题，提升检测性能的同时，降低检测系统的复杂度。
[0005]
为实现上述目的，本发明提供了一种高空抛掷物判别系统，包括：
[0006]
视频图像采集模块，用于拍摄建筑物楼面；
[0007]
网络传输模块，用于传输拍摄建筑物楼面的画面；
[0008]
运动目标检测模块，用于接收来自网络传输模块的建筑物楼面的画面，以及检测并提取拍摄画面中的运动目标；所述运动目标检测模块与视频图像采集模块连接；
[0009]
图像处理模块，用于提取运动目标的运动轨迹坐标序列并输入高空抛物判别模块；所述图像处理模块与所述运动目标检测模块连接；
[0010]
高空抛物判别模块，用于判别输入的运动目标是否属于高空抛物，所述高空抛物判别模块与图像处理模块连接，所述高空抛物判别模块中建立有用于判别输入的运动目标是否属于高空抛物的神经网络判别模型。
[0011]
进一步，所述高空抛物判别模块还包括模型迭代模块，用于将属于高空抛物的运
动目标的运动特征加入到神经网络判别模型的训练集中，进行模型迭代。
[0012]
进一步，所述神经网络判别模型通过以下步骤获得：
[0013]
(s41)通过计算机仿真，获得不同气动外形的抛掷物在不同风场条件下的运动轨迹坐标序列样本；
[0014]
(s42)将抛掷物运动轨迹坐标序列样本集合作为训练集，对判别神经网络进行训练，以得到训练好的神经网络判别模型。
[0015]
进一步，所述通过计算机仿真，获得不同气动外形的抛掷物在不同模拟风场条件下的运动轨迹样本，具体包括以下步骤：
[0016]
(s411)通过计算机仿真，随机生成1个抛掷物，抛掷物形态在块、棒、板这3种形态中随机确定；
[0017]
(s412)初始化抛掷物运动状态及设定计算抛掷物几何中心点坐标的迭代时间步长δt及总迭代步数，计算并记录每个迭代时间步上的几何中心点坐标，将抛掷物在初始状态下的几何中心点坐标和每个迭代时间步上的几何中心点坐标组成1个坐标序列p
tag
保存至数据库；
[0018]
(s413)通过设置不同气动外形的抛掷物在不同风场条件下运动，获得运动轨迹样本p
tag1
、p
tag2
···
p
tagm
，生成1个轨迹样本集p＝{p
tag1
,p
tag2
,...,p
tagm
}；可选的，所述初始化抛掷物运动状态，具体包括以下内容：随机设定抛掷物初始速度矢量其角度任意选择，择，的单位为m/s；抛掷物旋转速度模拟时间与风场模拟时间同步；抛掷物运动起始点为模拟风场上部边缘正中；其中，板状抛掷物迎风角随机取值，
[0019]
进一步，所述计算并记录每个迭代时间步上的几何中心点坐标，将抛掷物在初始状态下的几何中心点坐标和每个迭代时间步上的几何中心点坐标组成1个坐标序列p
tag
保存至数据库；具体执行以下步骤：
[0020]
(q1)从模拟风场中获取在t
n-1
时刻抛掷物运动位置处风速矢量结合抛掷物形态特征参数，计算t
n-1
时刻抛掷物的表面风力和其中，和分别为t
n-1
时刻抛掷物在空间3个维度x、y和z方向上的表面风力，为t
n-1
时刻抛掷物所受翻转作用力；其中，t
n
＝t0+nδt，n大于等于1，小于等于总迭代步数；当n＝1时，t
n-1
＝t0，t0表示抛掷物初始状态下的初始时刻，t
n
表示当前经过n次迭代时间步后的当前时刻；
[0021]
(q2)根据抛掷物的表面风力和分别计算抛掷物加速度分别计算抛掷物加速度和其中，和分别为抛掷物在空间3个维度x、y和z上受力所产生的加速度，为抛掷物翻转加速度；
[0022]
(q3)通过t
n-1
时刻抛掷物加速度和t
n-1
时刻抛掷物在空间3个维度x、y和z上的速度和分别计算t
n
时刻抛掷物在空间
(x
i
,t
n-1
)在v、u和w方向上的风速脉冲分量；ξ为1组符合标准正态分布的随机数；σ则根据边界层大气状态计算得出，σ
v'
、σ
u'
和σ
w'
分别为σ在v、u和w方向上的取值；t
l
为t
n
时刻风速v(x
i
,t
n
)的脉动分量v'(x
i
,t
n
)的拉格朗日时间尺度，t
lv'
、t
lu'
和t
lw'
分别为t
l
在v、u和w方向上拉格朗日时间尺度；δt为计算时间步长；r(δt)为指数形式的相关系数。
[0039]
进一步，所述边界层条件包括不稳定边界层、中性边界层和稳定边界层，分别通过步骤(d1)至(d3)获得对应不稳定边界层、中性边界层和稳定边界层的模拟风场的数值模型v
x,t
；可选的，对于不稳定边界层，则σ
u'
、σ
v'
、σ
w'
、t
lu'
、t
lv'
和t
lw'
的公式分别为：
[0040]
σ
u
′
＝σ
v
′
＝u
*
(12+0.5z
i
/|l|)
1/3
；
[0041][0042]
t
lu
′
＝t
lv
′
＝0.15z
i
/σ
u
′
；
[0043][0044]
对于中性边界层，则σ
u'
、σ
v'
、σ
w'
、t
lu'
、t
lv'
和t
lw'
的公式分别为：
[0045]
σ
u
′
＝2u
*
exp(-3fz/u
*
)；
[0046]
σ
v
′
＝σ
w
′
＝1.3u
*
exp(-2fz/u
*
)；
[0047][0048]
对于稳定边界层，则σ
u'
、σ
v'
、σ
w'
、t
lu'
、t
lv'
和t
lw'
的公式分别为：
[0049]
σ
u
′
＝2u
*
(1-z/z
i
)；
[0050]
σ
v
′
＝σ
w
′
＝1.3u
*
(1-z/z
i
)；
[0051][0052][0053][0054]
其中，z为计算点高度，u
*
为摩擦速度，w
*
为对流特征速度。z
i
为混合层高度，l为莫宁-奥布霍夫长度，根据边界条件取值。f为科氏力长度，根据我国的纬度分布，取7.29
×
10-5
。
[0055]
进一步，训练集的构建步骤包括：
[0056]
(s421)在模拟环境中根据实际使用环境中摄像机布设的位置模拟摄像机机位，模
拟摄像机光轴与模拟楼面的夹角等于实际环境中摄像机光轴与楼面的夹角；
[0057]
(s422)模拟抛掷物运动轨迹，获得各种抛掷物形态及各种风场条件下，各种抛掷物在模拟摄像机中的轨迹图像；
[0058]
(s423)模拟非抛掷物运动轨迹，随机生成横向飞行目标、纵向上升目标、横向往返运动目标、纵向往复运动目标、随机闪烁目标，并记录其运动轨迹，获得各种类型的非抛掷物及各种风场条件下，各类非抛掷物运动目标的运动轨迹图像；
[0059]
(s424)将抛掷物运动轨迹图像样本与非抛掷物运动轨迹图像样本按1:1的比例混合；
[0060]
(s425)将混合后的样本集按训练样本数量：测试样本数量＝7:3的比例进行采样，构建训练数据集和测试数据集。
[0061]
进一步，抛掷物和非抛掷物均通过以下步骤获得各种类型抛掷物或非抛掷物及各种风场条件下的运动轨迹：
[0062]
(f1)设摄像机的坐标为(x
c
,y
c
,z
c
)，非抛掷物或抛掷物的运动轨迹通过仿真得出，设t时刻非抛掷物或抛掷物的坐标为(x
tag
,y
tag
,z
tag
)，非抛掷物或抛掷物轨迹点在摄像机成像画面的坐标(x
′
tag
,y
′
tag
，z
′
tag
)，x
′
tag
和y
′
tag
的计算公式分别为：
[0063][0064][0065]
其中，f为摄像机镜头焦距；
[0066]
(f2)重复执行步骤(f1)，直至1条轨迹上所有的非抛掷物或抛掷物位置点均在模拟摄像机中成像并形成轨迹图像，然后执行步骤(f3)；
[0067]
(f3)重复步骤(f1)和步骤(f2)，获得各种抛掷物形态各种风场条件下，各种非抛掷物或抛掷物在模拟摄像机中的轨迹图像；可选的，所述神经网络判别模型通过卷积神经网络训练得到；可选的，所述神经网络判别模型分为输入层、特征抽取层和判别输出层；输入层接收输入的运动轨迹坐标序列，并向特征抽取层馈送；特征抽取层由3个卷积池化层堆叠而成，即每个卷积层之后接1个1维池化层；判别输出层由1个全连接层和1个softmax层构成，softmax进行二分类判别，即输入轨迹是否为高空抛掷物运动轨迹；可选的，输入的运动轨迹坐标序列长度为125，超过125的序列保留从第1位到125位数据，不足125的序列在序列末尾补0；可选的，每个卷积池化层中，卷积层卷积核大小为1
×
5，卷积核步长为1；池化大小为1
×
5，步长为2。
[0068]
本发明与现有技术相比较具有以下优点：
[0069]
本发明的高空抛掷物判别系统，将高空抛物目标运动检测与识别转化为一个序列分类问题，提升检测性能的同时，降低检测系统的复杂度；还解决了在轨迹样本难以获取的问题；系统简单，大大降低了计算，降低了成本，且能够模拟实际场景中多变的风场条件，能够获取到不同气动外形的抛掷物在不同风场下的运动轨迹，使得模型在训练过程中获取到有效的训练样本数据，提升了模型的判别性能。
附图说明
[0070]
图1为本发明高空抛掷物判别方法的结构示意图；
[0071]
图2为本发明高空抛掷物判别系统的结构示意图。
[0072]
图中：
[0073]
1-视频图像采集模块；2-运动目标检测模块；3-图像处理模块；4-高空抛物判别模块；41-模型迭代模块；5-网络传输模块。
具体实施方式
[0074]
下面结合附图对本发明的具体实施方式作进一步说明。
[0075]
参见图1所示，本实施例公开一种高空抛掷物判别方法，包括以下步骤：
[0076]
(s1)摄像机对建筑物楼面进行拍摄；
[0077]
(s2)使用帧差法检测并提取摄像机画面中的运动目标；
[0078]
(s3)提取运动目标的运动轨迹坐标序列；
[0079]
(s4)将运动目标的运动轨迹坐标序列输入到预先训练好的神经网络判别模型中，判别该运动目标是否属于高空抛物。
[0080]
在本实施例中，若运动目标属于高空抛物，则将该运动目标的运动特征加入到神经网络判别模型的训练集中，进行模型迭代。
[0081]
在本实施例中，所述预先训练好的神经网络判别模型通过以下步骤获得：
[0082]
(s41)通过计算机仿真，获得不同气动外形的抛掷物在不同风场条件下的运动轨迹坐标序列样本；
[0083]
(s42)将抛掷物运动轨迹坐标序列样本集合作为训练集，对判别神经网络进行训练，以得到训练好的神经网络判别模型。
[0084]
在本实施例中，所述通过计算机仿真，获得不同气动外形的抛掷物在不同模拟风场条件下的运动轨迹样本，具体包括以下步骤：
[0085]
(s411)通过计算机仿真，随机生成1个抛掷物，抛掷物形态在块、棒、板这3种形态中随机确定；
[0086]
(s412)初始化抛掷物运动状态及设定计算抛掷物几何中心点坐标的迭代时间步长δt及总迭代步数，计算并记录每个迭代时间步上的几何中心点坐标，将抛掷物在初始状态下的几何中心点坐标和每个迭代时间步上的几何中心点坐标组成1个坐标序列p
tag
保存至数据库；
[0087]
(s413)通过设置不同气动外形的抛掷物在不同风场条件下运动，获得运动轨迹样本p
tag1
、p
tag2
···
p
tagm
，生成1个轨迹样本集p＝{p
tag1
,p
tag2
,...,p
tagm
}。
[0088]
在本实施例中，所述初始化抛掷物运动状态，具体包括以下内容：随机设定抛掷物初始速度矢量其角度任意选择，其角度任意选择，的单位为m/s；抛掷物旋转速度
[0089]
模拟时间与风场模拟时间同步；
[0090]
抛掷物运动起始点为模拟风场上部边缘正中。其中，板状抛掷物迎风角随机取值，板状抛掷物迎风角是用于输入fluent计算用。
[0091]
在本实施例中，所述计算并记录每个迭代时间步上的几何中心点坐标，将抛掷物在初始状态下的几何中心点坐标和每个迭代时间步上的几何中心点坐标组成1个坐标序列p
tag
保存至数据库；具体执行以下步骤：
[0092]
(q1)从模拟风场中获取在t
n-1
时刻抛掷物运动位置处风速矢量结合抛掷物形态特征参数，计算t
n-1
时刻抛掷物的表面风力和其中，和分别为t
n-1
时刻抛掷物在空间3个维度x、y和z方向上的表面风力，为t
n-1
时刻抛掷物所受翻转作用力；其中，t
n
＝t0+nδt，n大于等于1，小于等于总迭代步数；当n＝1时，t
n-1
＝t0，t0表示抛掷物初始状态下的初始时刻，t
n
表示当前经过n次迭代时间步后的当前时刻；
[0093]
(q2)根据抛掷物的表面风力和分别计算抛掷物加速度分别计算抛掷物加速度和其中，和分别为抛掷物在空间3个维度x、y和z上受力所产生的加速度，为抛掷物翻转加速度；
[0094]
(q3)通过t
n-1
时刻抛掷物加速度和t
n-1
时刻抛掷物在空间3个维度x、y和z上的速度和分别计算t
n
时刻抛掷物在空间3个维度x、y和z上的速度和和位移和以及抛掷物翻转速度和位移进而得到t
n
时刻下当前迭代时间步上的几何中心点坐标；
[0095]
(q4)迭代开始时，按照设定的时间步长δt依次序迭代，每次迭代均需重复执行步骤(q1)至步骤(q3)，直至达到总迭代步数停止，记录下抛掷物在每个迭代时间步上的几何中心点坐标，将抛掷物在初始状态下的几何中心点坐标和每个迭代时间步上的几何中心点坐标组成1个坐标序列p
tag
保存至数据库；其中，达到预设迭代步数即此时的抛掷物落地或飞出模拟风场边界；抛掷物落地即v
tag,y,t
＝0。
[0096]
在本实施例中，对于块状抛掷物，随机生成球体、球状多面体和立方体抛掷物，球体和球状多面体直径在[10cm,50cm]范围内随机选取，立方体边长在[10cm,30cm]范围内随机选取；
[0097]
对于棒状抛掷物，随机生成圆柱体、多棱柱抛掷物，长度在[30cm,300cm]内随机选择，截面直径在[5cm,20cm]范围内随机选取；
[0098]
对于板状抛掷物，随机生成矩形、圆形和多边形板状抛掷物，面积在[0.1m2,1m2]内随机选择；
[0099]
其中，抛掷物质量随机选取，m∈[50,5000]，质量单位为g。
[0100]
在本实施例中，所述表面风力和的计算公式分别为：
[0101][0102]
[0103][0104][0105]
其中，ρ
a
为空气密度，a为抛掷物最大参考面积，为t
n-1
时刻抛掷物运动速度矢量，为抛掷物旋转速度矢量；为t
n-1
时刻下抛掷物运动位置处风速矢量；和分别为t
n-1
时刻抛掷物在空间3个维度x、y和z方向上的表面风力系数，为t
n-1
时刻抛掷物所受翻转作用力对应表面风力系数。在本实施例中，每次迭代时，以下数据需进行更新：抛掷物的位置坐标；抛掷物的速度矢量抛掷物当前运动位置处风速矢量抛掷物在t
n-1
时刻抛掷物表面风力系数和和和是根据t
n-1
时刻下抛掷物运动位置处风速矢量并结合抛射物形态特征参数，采用rans均时应力模型、cfd计算软件fluent计算得到。t
n
时刻抛掷物的表面风力和采用rans均时应力模型及cfd计算软件fluent进行计算。
[0106]
在本实施例中，计算抛掷物加速度和的公式分别为：
[0107][0108][0109][0110][0111]
计算抛掷物在空间3个维度x、y和z上的速度和和位移和位移和以及抛掷物翻转速度和位移的公式分别为：
[0112][0113][0114][0115][0116][0117]
[0118][0119][0120]
在本实施例中，所述模拟风场生成步骤如下：
[0121]
(d1)设置一个用于模拟风场的计算空间，并设定模拟风场的边界层条件；
[0122]
(d2)根据设定模拟风场的边界层条件，计算各计算点在t
n
时刻、坐标x
i
时的风速在x,y,z方向上的风速分量v(x
i
,t
n
)、u(x
i
,t
n
)和w(x
i
,t
n
)，最终得到在t
n
时刻各计算点的合成风速v(x
i
,t
n
)，其中，t
n
＝t0+nδt，n大于等于1，小于等于总迭代步数；当n＝1时，t
n-1
＝t0，t0表示抛掷物初始状态下的初始时刻，t
n
表示当前经过n次迭代时间步后的当前时刻；
[0123]
(d3)根据设定的时间步长δt依次序迭代重复执行步骤(d1)至步骤(d2)，直到达到总迭代步数后停止计算，得到对应该模拟风场边界层条件的模拟风场的数值模型v
x,t
。
[0124]
在本实施例中，所述计算空间的大小为高120m，宽300m，纵深300m；在某些实施例中，所述计算空间还可以为其他大小，在此不作限定。各计算点划分：根据计算能力，确定风场计算点网格尺寸，模拟网格大小在1mx1m至5mx5m内可选。不稳定边界层的风场模拟选择小尺度、较密集的网格，中性和稳定边界层的模拟选择大尺度、稀疏的网格。确认风场模拟时间与抛掷物运动状态下总迭代步数相同。
[0125]
在本实施例中，在[0m/s,5m/s]之间进行随机设定。
[0126]
在本实施例中，通过公式(1)至公式(4)可以推导出x,y,z方向上的风速分量v(x
i
,t
n
)、u(x
i
,t
n
)和w(x
i
,t
n
)的计算公式，具体推算过程如下：
[0127]
模拟风场中每个计算点在t时刻的风速：
[0128][0129]
式(1)可以拆分为3个速度分量来表示：
[0130][0131]
上述3个速度分量生成合成风速的公式为：
[0132][0133]
其中，将气流涡旋过程作为1个连续过程，在遵从markov假定条件下，经过1个时间步长，计算点风速为：
[0134]
v
′
(x
i
,t
n
)＝r(δt)v
′
(x
i
,t
n-1
)+v
″
(x
i
,t
n-1
)
ꢀꢀꢀꢀꢀꢀꢀ
(2)；
[0135]
式(2)可以拆分为3个速度分量来表示：
[0136][0137]
r(δt)＝exp(-δt/t
l
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)；
[0138]
v
″
(x
i
,t
n
)＝σ[1-r2(δt)]
1/2
ξ
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)；
[0139]
综上所述，得到风速计算公式：在x,y,z方向上的风速分量v(x
i
,t
n
)、u(x
i
,t
n
)和w(x
i
,t
n
)，拆分后的v(x
i
,t
n
)、u(x
i
,t
n
)和w(x
i
,t
n
)的公式为：
[0140][0141][0142][0143]
对于v方向上，r(δt)＝exp(-δt/t
lv'
)，v”(x
i
,t
n-1
)＝σ
v
′
[1-r2(δt)]
1/2
ξ；
[0144]
对于u方向上，r(δt)＝exp(-δt/t
lu'
)，u”(x
i
,t
n-1
)＝σ
u
′
[1-r2(δt)]
1/2
ξ；
[0145]
对于w方向上，r(δt)＝exp(-δt/t
lv'
)，v”(x
i
,t
n-1
)＝σ
v
′
[1-r2(δt)]
1/2
ξ；
[0146]
其中，和为t
n
时刻的风速平均分量v(x
i
,t
n
)在v、u和w方向上的风速平均分量；v”(x
i
,t
n-1
)、u”(x
i
,t
n-1
)和w”(x
i
,t
n-1
)为t
n-1
时刻的风速脉冲分量v”(x
i
,t
n-1
)在v、u和w方向上的风速脉冲分量；ξ为1组符合标准正态分布的随机数；σ则根据边界层大气状态计算得出，σ
v'
、σ
u'
和σ
w'
分别为σ在v、u和w方向上的取值；t
l
为t
n
时刻风速v(x
i
,t
n
)的脉动分量v'(x
i
,t
n
)的拉格朗日时间尺度，t
lv'
、t
lu'
和t
lw'
分别为t
l
在v、u和w方向上拉格朗日时间尺度；δt为计算时间步长；r(δt)为指数形式的相关系数。此处平均分量原则上由1个诊断风场软件给出，但是由于抛物时间较短，风速的平均分量变化不大，因此在本方案中，每次仿真计算中，风速的平均分量均人为设定且保持恒定不变，因此和的取值在[1,10]区间内，的取值在[0,3]区间内。时间步长根据实际需要进行取值，在此不作限定，在本实施例中，时间步长取值为1s。在计算上，上面3个式子的中右边第二项需要对v'(x
i
,t
n
)、u'(x
i
,t
n
)、w'(x
i
,t
n
)进行迭代。采用上面三个式子进行描述时。在实际计算中v’的初始值为0或者1个很小的值。
[0147]
在本实施例中，所述边界层条件包括不稳定边界层、中性边界层和稳定边界层，分别通过步骤(d1)至(d3)获得对应不稳定边界层、中性边界层和稳定边界层的模拟风场的数值模型v
x,t
。
[0148]
在本实施例中，对于不稳定边界层，则σ
u'
、σ
v'
、σ
w'
、t
lu'
、t
lv'
和t
lw'
的公式分别为：
[0149]
σ
u
′
＝σ
v
′
＝u
*
(12+0.5z
i
/|l|)
1/3
；
[0150]
[0151]
t
lu
′
＝t
lv
′
＝0.15z
i
/σ
u
′
；
[0152][0153]
对于中性边界层，则σ
u'
、σ
v'
、σ
w'
、t
lu'
、t
lv'
和t
lw'
的公式分别为：
[0154]
σ
u
′
＝2u
*
exp(-3fz/u
*
)；
[0155]
σ
v
′
＝σ
w
′
＝1.3u
*
exp(-2fz/u
*
)；
[0156][0157]
对于稳定边界层，则σ
u'
、σ
v'
、σ
w'
、t
lu'
、t
lv'
和t
lw'
的公式分别为：
[0158]
σ
u
′
＝2u
*
(1-z/z
i
)；
[0159]
σ
v
′
＝σ
w
′
＝1.3u
*
(1-z/z
i
)；
[0160][0161][0162][0163]
其中，z为计算点高度，u
*
为摩擦速度，w
*
为对流特征速度。z
i
为混合层高度，l为莫宁-奥布霍夫长度(monin
–
obukhov length)，根据边界条件取值。f为科氏力长度，根据我国的纬度分布，取7.29
×
10-5
。
[0164]
在本实施例中，训练集的构建步骤包括：
[0165]
(s421)在模拟环境中根据实际使用环境中摄像机布设的位置模拟摄像机机位，模拟摄像机光轴与模拟楼面的夹角等于实际环境中摄像机光轴与楼面的夹角；
[0166]
(s422)模拟抛掷物运动轨迹，获得各种抛掷物形态及各种风场条件下，各种抛掷物在模拟摄像机中的轨迹图像；
[0167]
(s423)模拟非抛掷物运动轨迹，随机生成横向飞行目标、纵向上升目标、横向往返运动目标、纵向往复运动目标、随机闪烁目标，并记录其运动轨迹，获得各种类型的非抛掷物及各种风场条件下，各类非抛掷物运动目标的运动轨迹图像；
[0168]
(s424)将抛掷物运动轨迹图像样本与非抛掷物运动轨迹图像样本按1:1的比例混合；
[0169]
(s425)将混合后的样本集按训练样本数量：测试样本数量＝7:3的比例进行采样，构建训练数据集和测试数据集。
[0170]
在本实施例中，抛掷物和非抛掷物均通过以下步骤获得各种类型抛掷物或非抛掷物及各种风场条件下的运动轨迹：
[0171]
(f1)设摄像机的坐标为(x
c
,y
c
,z
c
)，非抛掷物或抛掷物的运动轨迹通过仿真得出，
设t时刻非抛掷物或抛掷物的坐标为(x
tag
,y
tag
,z
tag
)，非抛掷物或抛掷物轨迹点在摄像机成像画面的坐标(x
′
tag
,y
′
tag
，z
′
tag
)，x
′
tag
和y
′
tag
的计算公式分别为：
[0172][0173][0174]
其中，f为摄像机镜头焦距；采用针孔成像模型，在针孔成像模型中不对z轴上的数据进行计算。
[0175]
(f2)重复执行步骤(f1)，直至1条轨迹上所有的非抛掷物或抛掷物位置点均在模拟摄像机中成像并形成轨迹图像，然后执行步骤(f3)；
[0176]
(f3)重复步骤(f1)和步骤(f2)，获得各种抛掷物形态各种风场条件下，不同非抛掷物或抛掷物在模拟摄像机中的轨迹图像。
[0177]
在本实施例中，所述神经网络判别模型通过卷积神经网络训练得到。所述神经网络判别模型分为输入层、特征抽取层和判别输出层；输入层接收输入的坐标序列轨迹坐标序列，并向特征抽取层馈送；特征抽取层由3个卷积池化层堆叠而成，即每个卷积层之后接1个1维池化层；判别输出层由1个全连接层和1个softmax层构成，softmax进行二分类判别，即输入轨迹是否为高空抛掷物运动轨迹。
[0178]
在本实施例中，输入的运动轨迹坐标序列长度为125，超过125的序列保留从第1位到125位数据，不足125的序列在序列末尾补0；每个卷积池化层中，卷积层卷积核大小为1
×
5，卷积核步长为1；池化大小为1
×
5，步长为2。
[0179]
在本实施例中，神经网络使用binary cross entropy作为损失函数：
[0180][0181]
式中：
[0182]
y为轨迹样本的类别；
[0183]
p(yi)为神经网络给出的样本类别分类概率；
[0184]
n为样本个数。抛掷物在风场中的运动轨迹，通过摄像机拍摄后，离散化为1个轨迹坐标序列。将该序列输入1个判别网络，实现轨迹的分类判别。
[0185]
参见图2所示，本实施例还公开了一种高空抛掷物判别系统，执行上述高空抛掷物判别方法，包括：
[0186]
视频图像采集模块1，用于拍摄建筑物楼面；
[0187]
网络传输模块5，用于传输拍摄建筑物楼面的画面；
[0188]
运动目标检测模块2，用于接收来自网络传输模块5的建筑物楼面的画面，以及检测并提取拍摄画面中的运动目标；所述运动目标检测模块2与视频图像采集模块1连接；
[0189]
图像处理模块3，用于提取运动目标的运动轨迹坐标序列并输入高空抛物判别模块4；所述图像处理模块3与所述运动目标检测模块2连接；
[0190]
高空抛物判别模块4，用于判别输入的运动目标是否属于高空抛物，所述高空抛物判别模块4与图像处理模块3连接，所述高空抛物判别模块4中建立有用于判别输入的运动
目标是否属于高空抛物的神经网络判别模型。
[0191]
在本实施例中，所述高空抛物判别模块4包括模型迭代模块41，用于将属于高空抛物的运动目标的运动特征加入到神经网络判别模型的训练集中，进行模型迭代。
[0192]
本发明的高空抛掷物判别系统，将高空抛物目标运动检测与识别转化为一个序列分类问题，提升检测性能的同时，降低检测系统的复杂度；还解决了在轨迹样本难以获取的问题；系统简单，大大降低了计算，降低了成本，且能够模拟实际场景中多变的风场条件，能够获取到不同气动外形的抛掷物在不同风场下的运动轨迹，使得模型在训练过程中获取到有效的训练样本数据，提升了模型的判别性能。
[0193]
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。