本发明涉及公路交通事件检测领域，特别是涉及基于趋势性指标及波动性指标的交通事件自动检测方法。
背景技术：
名词解释：
占有率：
峰度：又称峰态系数，表征概率密度分布曲线在平均值处峰值高低的特征数；
偏度：是统计数据分布偏斜方向和程度的度量，是统计数据分布非对称程度的数字特征；
上下游：本申请中指交通线路的某一地点的上游路段和下游路段；
脏数据：是指源系统中的数据不在给定的范围内或对于实际业务毫无意义，或是数据格式非法，以及在源系统中存在不规范的编码和含糊的业务逻辑。
交通拥堵问题已经成为了阻碍城市发展的一大问题，不仅造成环境污染，还影响人们的出行与健康。除了机动车数量本身超过道路通行能力造成交通压力外，道路上的各种交通事件也是造成道路拥堵的主要原因。当交通事件发生时，交通状态发生突变，道路的服务水平陡转急下，不仅导致行车延误，还常伴有直接的经济财产损失和人员伤亡。快速有效的交通事件检测能够及时发现交通事故发生的时间、地点，可以让交通管理者在尽快的赶到现场，并做出适当的调整，减少对出行者的影响，同时也为事故伤亡人员赢得宝贵的救治时间，将人员伤亡和经济财产损失降到最低。
目前技术中，较为典型的交通事件检测算法是加州算法，于20世纪60年代由美国加州运输部开发，并得到广泛承认和应用，一般作为评价新算法的参考。该算法属双截面算法，它基于事件发生时上游检测截面占有率将增加，下游检测截面占有率将减少这一事实，通过计算平均占有率，来进行交通事件检测。但是，这种传统的时间检测方法主要利用了空间尺度上占有率的变化，包括上下游占有率的相对差值和绝对差值，在时间尺度上只考虑了占有率的相对差值。由于占有率不能直接体现出路段的是否交通事件，终端报警存在很多误报和漏报的情况，故需要进行数据探索和预处理，总结发生交通事件的占有率变化规律，再从数据中提炼出描述交通事件的特征指标。当道路中出现检测器故障时，会导致其上下游都无法利用该方法进行检测。同时，交通事件的开始时间和结束时间是表征事件的关键时间节点，在这些时间节点上，占有率和终端报警等数据也会有一定的特征变化，传统的加州算法并未考虑关键时间节点前后一定范围的数据波动性特征指标。总的来说，目前的检测方法无法有效地对交通事件进行检测，从而无法及时地获知交通状态的变化。
技术实现要素：
为了解决上述的技术问题，本发明的目的是提供基于趋势性指标及波动性指标的交通事件自动检测方法。
本发明解决其技术问题所采用的技术方案是：
基于趋势性指标及波动性指标的交通事件自动检测方法，包括步骤：
S1、通过传感器采集实时交通数据；
S2、对实时交通数据进行预处理；
S3、基于预处理后的实时交通数据，计算如下的实时特征向量：趋势性指标、波动性指标和上下游变化指标；
S4、将计算获得的实时特征向量作为训练模型的输入序列，采用训练模型计算获得对应的输出结果作为交通事件的检测结果。
进一步，还包括以下步骤：
S0、从交通数据库系统中获取历史交通数据及对应的告警数据后，构建训练模型并根据获取的数据进行模型训练和评价。
进一步，所述步骤S0包括以下步骤：
S01、从交通数据库系统中获取历史交通数据及对应的告警数据构成历史样本数据；
S02、对历史样本数据进行预处理；
S03、基于预处理后的历史样本数据，计算如下的历史特征向量：趋势性指标、波动性指标和上下游变化指标；
S04、构建训练模型，将计算获得的历史特征向量作为训练模型的输入序列，将对应的告警数据作为训练模型的输出序列，对训练模型进行训练及评价。
进一步，所述步骤S04，其具体为：
构建训练模型，将计算获得的历史特征向量作为训练模型的输入序列，将对应的告警数据作为训练模型的输出序列，采用其中80％的训练数据对训练模型进行训练，并用剩下20％的训练数据对训练模型进行评价，从而进行模型优化与重构，直到满足预设的训练参数。
进一步，所述步骤S04中所述训练模型采用LM神经网络模型或CART决策树模型。
进一步，所述交通数据包括速度、占有率和/或交通流量。
进一步，所述趋势性指标是对连续H个时刻获取到的交通数据序列进行以下计算所获得的三个参数：
一、采用最小二乘法对连续H个时刻获取到的交通数据序列拟合得到的曲线的斜率k；
二、计算交通数据序列中的递减趋势时刻数：
上式中，T表示趋势指标，D(i)表示前后两个时刻的交通数据的递减趋势，满足下式：
其中，ki表示当前时刻的交通数据，ki-1表示前一时刻的交通数据；
三、计算前H/2个时刻和后H/2个时刻的平均交通数据的差值：
上式中，davg_H表示计算获得的差值。
进一步，所述波动性指标包括传感器连续H个时刻获取到的交通数据序列的标准差、峰度和偏度。
进一步，所述上下游变化指标包括上下游交通数据序列的绝对差值和相对差值。
进一步，所述步骤S2，其具体为：
对实时交通数据剔除脏数据和误报警数据后，进行数据清洗、缺失值处理。
本发明的有益效果是：本发明的基于趋势性指标及波动性指标的交通事件自动检测方法，包括步骤：S1、通过传感器采集实时交通数据；S2、对实时交通数据进行预处理；S3、基于预处理后的实时交通数据，计算如下的实时特征向量：趋势性指标、波动性指标和上下游变化指标；S4、将计算获得的实时特征向量作为训练模型的输入序列，采用训练模型计算获得对应的输出结果作为交通事件的检测结果。本方法通过传感器采集实时交通数据，可以对交通事件进行检测判断，及时获取交通状态的变化，及时、准确发现交通事件发生的时间、地点。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是本发明的基于趋势性指标及波动性指标的交通事件自动检测方法的原理框图。
具体实施方式
实施例一
参照图1，本发明提供了一种基于趋势性指标及波动性指标的交通事件自动检测方法，包括步骤：
S1、通过传感器采集实时交通数据；
S2、对实时交通数据进行预处理；
S3、基于预处理后的实时交通数据，计算如下的实时特征向量：趋势性指标、波动性指标和上下游变化指标；
S4、将计算获得的实时特征向量作为训练模型的输入序列，采用训练模型计算获得对应的输出结果作为交通事件的检测结果。
进一步作为优选的实施方式，还包括以下步骤：
S5、根据交通事件的检测结果进行及时告警。
进一步作为优选的实施方式，还包括以下步骤：
S0、从交通数据库系统中获取历史交通数据及对应的告警数据后，构建训练模型并根据获取的数据进行模型训练和评价。
进一步作为优选的实施方式，所述步骤S0包括以下步骤：
S01、从交通数据库系统中获取历史交通数据及对应的告警数据构成历史样本数据；
S02、对历史样本数据进行预处理，具体为：对历史样本数据剔除脏数据和误报警数据后，进行数据清洗、缺失值处理；
S03、基于预处理后的历史样本数据，计算如下的历史特征向量：趋势性指标、波动性指标和上下游变化指标；
S04、构建训练模型，将计算获得的历史特征向量作为训练模型的输入序列，将对应的告警数据作为训练模型的输出序列，对训练模型进行训练及评价。对应的告警数据是指告警状态标识，例如，当某个历史特征向量对应的告警状态为出现事故，则该告警数据为1，反之，若告警状态为未出现故障，则告警数据为0。
步骤S02和S03中实现的处理过程与步骤S2和S3相同，区别仅在于处理数据的不同，步骤S02和S03是对获取的历史交通数据的处理过程，采用与实际检测过程中一样的处理方式。
进一步作为优选的实施方式，所述步骤S04，其具体为：
构建训练模型，将计算获得的历史特征向量作为训练模型的输入序列，将对应的告警数据作为训练模型的输出序列，采用其中80％的训练数据对训练模型进行训练，并用剩下20％的训练数据对训练模型进行评价，从而进行模型优化与重构，直到满足预设的训练参数。预设的训练参数一般是设定训练模型的预测误差在某一预设范围，例如5％-10％之间的任一参数。具体的，当预设的训练参数为他预测误差在6％的范围内，则当对训练模型评价，将训练数据中的历史特征向量输入训练模型，计算得到输出序列与实际告警数据的训练误差小于6％，不再模型优化与重构。
进一步作为优选的实施方式，所述步骤S04中所述训练模型采用LM神经网络模型或CART决策树模型。
进一步作为优选的实施方式，所述交通数据包括速度、占有率和/或交通流量。
进一步作为优选的实施方式，所述趋势性指标是对连续H个时刻获取到的交通数据序列进行以下计算所获得的三个参数：
一、采用最小二乘法对连续H个时刻获取到的交通数据序列拟合得到的曲线的斜率k；
二、计算交通数据序列中的递减趋势时刻数：
上式中，T表示趋势指标，D(i)表示前后两个时刻的交通数据的递减趋势，满足下式：
其中，ki表示当前时刻的交通数据，ki-1表示前一时刻的交通数据；
三、计算前H/2个时刻和后H/2个时刻的平均交通数据的差值：
上式中，davg_H表示计算获得的差值。
进一步作为优选的实施方式，所述波动性指标包括传感器连续H个时刻获取到的交通数据序列的标准差、峰度和偏度。
进一步作为优选的实施方式，所述上下游变化指标包括上下游交通数据序列的绝对差值和相对差值。
进一步作为优选的实施方式，所述步骤S2，其具体为：
对实时交通数据剔除脏数据和误报警数据后，进行数据清洗、缺失值处理。
本实施例可以通过获取交通系统的传感器采集的实时交通数据，对交通事件进行检测判断，及时、准确发现事件发生的时间、地点，从而可以用于排除事故、发布信息等交通应用领域中，帮助故障路段尽快恢复正常运行。
实施例二
本实施例是实施例一的详细实例。本实施例以采集的交通数据为占有率为例，仔细说明训练模型的预测过程。其它计算过程等，与训练过程原理相似，可参照实施例一的描述。
通过连续H个时刻获取到的占有率序列，可以计算获得如下的历史特征向量：
趋势性指标：
一、连续H个时刻获取到的占有率序列，采用最小二乘法拟合曲线，计算得到的斜率k；
二、计算占有率序列中的递减趋势时刻数，即计算连续H个时刻获取到的占有率序列内，当前时刻比前一个时刻趋势为递减的时刻数：
上式中，T表示趋势指标，D(i)表示前后两个时刻的占有率的递减趋势，满足下式：
其中，ki表示当前时刻的占有率，ki-1表示前一时刻的占有率；
三、计算前H/2个时刻和后H/2个时刻的平均占有率的差值：
上式中，davg_H表示计算获得的差值。
波动性指标：
(1)传感器连续H个时刻获取到的占有率数据序列的标准差；
(2)传感器连续H个时刻获取到的占有率数据序列的峰度；
(3)传感器连续H个时刻获取到的占有率数据序列的偏度。
上下游变化指标：
(1)上下游占有率的绝对差值：
OCCDF＝OCC(i,t)-OCC(i+1,t)；
(2)上下游占有率的相对差值：
其中，OCCDF表示上下游占有率的绝对差值，OCCRDF表示上下游占有率的相对差值，OCC(i，t)表示第i个检测站t时刻所测得的占有率。
本实施例提取了时间尺度上占有率的波动性指标，这些指标只需要利用单个传感器就可以计算得到，不需要考虑空间尺度上前后检测的占有率差值，即使在单个传感器出现故障时，也不影响路段上下游的检测。同时，这些波动性指标充分提取了事件下的占有率变化特征，对交通事件检测算法准确率的提高有很大的帮助。
以上是对本发明的较佳实施进行了具体说明，但本发明创造并不限于所述实施例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做出种种的等同变形或替换，这些等同的变型或替换均包含在本申请权利要求所限定的范围内。