CN113706867B - 一种高速公路相邻边缘设备节点交通异常概率评估方法 - Google Patents

一种高速公路相邻边缘设备节点交通异常概率评估方法 Download PDF

Info

Publication number
CN113706867B
CN113706867B CN202110995620.3A CN202110995620A CN113706867B CN 113706867 B CN113706867 B CN 113706867B CN 202110995620 A CN202110995620 A CN 202110995620A CN 113706867 B CN113706867 B CN 113706867B
Authority
CN
China
Prior art keywords
traffic flow
expressway
traffic
time
adjacent edge
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110995620.3A
Other languages
English (en)
Other versions
CN113706867A (zh
Inventor
郭军
蒙龙
高赫妍
闫泽宇
曹佳慧
王馨悦
张斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northeastern University China
Original Assignee
Northeastern University China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Northeastern University China filed Critical Northeastern University China
Priority to CN202110995620.3A priority Critical patent/CN113706867B/zh
Publication of CN113706867A publication Critical patent/CN113706867A/zh
Application granted granted Critical
Publication of CN113706867B publication Critical patent/CN113706867B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/065Traffic control systems for road vehicles by counting the vehicles in a section of the road or in a parking area, i.e. comparing incoming count with outgoing count

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明提供一种高速公路相邻边缘设备节点交通异常概率评估方法,涉及智慧交通及云计算技术领域。云服务器从高速公路交通流量时间序列库中提取交通流量计数结果,并基于形态相似性距离计算相邻边缘设备节点交通流量平均到达时间;生成高速公路相邻边缘设备节点的交通流量差时间序列;采用基于ADF的高速公路相邻边缘设备节点交通流量差时间序列平稳性检验方法,检验相差平均到达时间的两个相邻边缘设备节点的交通流量差时间序列是否平稳;最后根据相邻边缘设备节点交通流量差时间序列的平稳性判定高速公路相邻边缘设备节点的交通异常情况。本发明方法充分考虑了交通异常和交通流量之间相互影响、相互作用,提高了高速公路交通异常判断的准确性。

Description

一种高速公路相邻边缘设备节点交通异常概率评估方法
技术领域
本发明涉及智慧交通及云计算技术领域,尤其涉及一种高速公路相邻边缘设备节点交通异常概率评估方法。
背景技术
随着高速公路的建设发展,道路上私家车和公共汽车的数量迅速增加,交通拥堵,车辆碰撞等交通异常问题越发严重。我国高速公路交通异常基本特点主要有:
(1)高速公路交通事故重大、恶性事故多。高速公路上的车辆行驶速度为60-120千米/时,车辆行驶速度过快,发生车辆碰撞时,冲击破坏力大,容易造成车辆与部件损坏,多为恶性事故。
(2)夜间交通异常事件发生较多。高速公路上夜间光线较暗,驾驶员在驾驶车辆时周围的视野条件与白天相差较大,容易造成视线模糊,大大增加了行车危险性,交通异常事件发生概率较大。同时夜间行车时,驾驶员容易疲劳驾驶,同样容易造成交通异常。
(3)异常天气交通异常事件较多。天气状况对高速公路上车辆的行车安全有很大影响,在晴天天气下交通异常的发生概率远小于异常(雾霾,雨雪等)天气下的交通异常发生概率。雾霾天气下驾驶员视线可见范围减小,雨雪天气下道路路面摩擦阻力减小,增加了行车危险性,增大高速公路交通异常概率。
交通异常发生通常造成不可挽回的损失,因此高速公路交通安全始终是人们关注的重点。当前智能交通系统中用于高速公路交通异常监测的常用技术为视频监测技术,但是视频检测设备受外界环境(如光照不足或者雾霾天气)影响较大,并且视频连续传输使网络传输带宽压力增大。通过车辆碰撞的声音检测交通异常是最直接的检测方法,但是事故地点距离音频设备较远时,音频设备可能无法采集到碰撞声音从而判断失误。
随着物联网技术的快速发展,终端设备海量数据的不断产生,边缘计算的提出和发展,促进了智能交通系统(Intelligent Traffic System,ITS)进一步改进和完善。目前,有许多学者通过边缘计算设备对高速公路交通异常情况进行分析,但前人对于高速公路交通异常的研究忽略了交通异常和交通流量之间相互影响、相互作用,导致效果并不理想。
发明内容
本发明要解决的技术问题是针对上述现有技术的不足,提供一种高速公路相邻边缘设备节点交通异常概率评估方法,实现对高速公路交通异常情况进行评估。
为解决上述技术问题,本发明所采取的技术方案是:一种高速公路相邻边缘设备节点交通异常概率评估方法,包括以下步骤:
步骤1:云服务器从高速公路交通流量时间序列库中提取高速公路交通流量计数结果,并基于形态相似性距离计算高速公路相邻边缘设备节点交通流量平均到达时间;
步骤1.1:设定高速公路交通流量平均到达时间表示交通流量序列从第i个边缘设备节点到第i+1个边缘设备节点所用时间,车辆在高速公路上行驶速度为60~120km/h,同时设定相邻两个边缘设备节点的距离为L,则交通流量平均到达时间范围为
Figure GDA0003645077240000021
步骤1.2:计算两个高速公路交通流量序列的形态相似距离,如下公式所示:
Figure GDA0003645077240000022
Figure GDA0003645077240000023
Figure GDA0003645077240000024
其中,D(Xi+1,t,Xi,t-Δt)表示Xi+1,t和Xi,t-Δt之间的形态相似距离,Xi+1,t表示为t时刻第i+1个边缘设备节点的高速公路交通流量序列,Δt表示为
Figure GDA0003645077240000025
时间段内的滑动时间间隔,且
Figure GDA0003645077240000026
Xi,t-Δt表示为t-Δt时刻第i个边缘设备节点高速公路交通流量序列,DED表示为Xi+1,t和Xi,t-Δt之间的欧氏距离,DHD表示为Xi+1,t和Xi,t-Δt的哈夫曼距离,n为高速公路交通流量序列上限;
步骤1.3:初始化两个高速公路交通流量序列最小形态相似距离和两个高速公路交通流量序列最小形态相似距离对应的时刻,如下公式所示:
Figure GDA0003645077240000027
Figure GDA0003645077240000028
其中,dist表示为两个高速公路交通流量序列最小形态相似距离,Tdist_min表示为两个高速公路交通流量序列最小形态相似距离对应的时刻;
步骤1.4:更新高速公路交通流量序列最小形态相似距离和高速公路交通流量序列最小形态相似距离对应的时刻;
步骤1.4.1:当tt时刻在
Figure GDA0003645077240000031
之间时,读取tt时刻第i个边缘设备节点高速公路交通流量序列Xi,tt
步骤1.4.2:计算t时刻第i+1个边缘设备节点的高速公路交通流量序列与tt时刻第i个边缘设备节点高速公路交通流量序列Xi,tt的形态相似距离disti+1,i=D(Xi+1,t,Xi,tt);
步骤1.4.3:根据以下条件更新高速公路交通流量序列最小形态相似距离dist和高速公路交通流量序列最小形态相似距离对应的时刻Tdist_min
当disti+1,i<dist时,令dist=disti+1,i,Tdist_min=tt,然后执行步骤1.5;
当disti+1,i≥dist时,令tt=tt+Δt,重新执行步骤1.4.1;
步骤1.5:基于高速公路交通流量序列最小形态相似距离对应的时刻,得到高速公路交通流量平均到达时间T,如下公式所示:
T=t-Tdist_min (6)
步骤2:生成高速公路相邻边缘设备节点的交通流量差时间序列
Figure GDA0003645077240000032
具体如下公式所示:
Zt=Xi+1,t-Xi,t-T (7)
Figure GDA0003645077240000033
其中,Zt表示为t时刻的高速公路相邻节点交通流量差时间序列,
Figure GDA0003645077240000034
表示t时刻第i+1个边缘设备节点的交通流量时间序列,
Figure GDA0003645077240000035
表示t时刻第i+1个边缘设备节点记录的第n个交通流量值,
Figure GDA0003645077240000036
表示为t-T时刻第i个边缘计算设备的交通流量时间序列,
Figure GDA0003645077240000037
表示t-T时刻第i个边缘设备节点记录的第n个交通流量值;
步骤3:采用基于ADF(Augment Dickey-Fuller)的高速公路相邻边缘设备节点交通流量差时间序列平稳性检验方法,检验相差平均到达时间的两个高速公路相邻边缘设备节点的交通流量差时间序列是否平稳;
步骤3.1:设定高速公路相邻边缘设备节点交通流量差时间序列Zt是一个m阶自回归过程生成的,则Zt表示为:
Figure GDA0003645077240000041
Zt=ρZt-11ΔZt-12ΔZt-2+…+ξmΔZt-m+1t (10)
其中,
Figure GDA0003645077240000042
表示自回归系数,m表示为自回归阶数,εt表示为随机误差项,且误差序列服从均值为0,方差为σ2
Figure GDA0003645077240000043
ΔZt-1,ΔZt-2,...,ΔZt-m+1表示为Zt-1,Zt-2,...,Zt-m+1的一阶差分;
步骤3.2:采用ADF检验的三种模型对高速公路相邻边缘设备节点交通流量差时间序列Zt进行单位根检验,则t时刻高速公路相邻边缘设备节点交通流量差时间序列的三个检验模型顺序如下:
不含常数项:
Figure GDA0003645077240000044
含常数项:
Figure GDA0003645077240000045
含趋势项:
Figure GDA0003645077240000046
其中,ΔZt表示为Zt的一阶差分,δ表示为常数,βiΔZt-i表示为Zt-i的一阶差分,α表示为常数项系数,βt表示为时间趋势项,m'表示为Zt的滞后阶数;
步骤3.3:设定在平稳性检验过程中,模型假设检验中原假设为H0:δ=0,即高速公路相邻边缘设备节点交通流量差时间序列是非平稳时间序列,备择假设H1:δ<0,即高速公路相邻边缘设备节点交通流量差时间序列是平稳时间序列;
步骤3.4:基于t时刻高速公路相邻边缘设备节点交通流量差时间序列的三个检验模型判断高速公路相邻边缘设备节点交通流量差时间序列是否为平稳时间序列;
步骤3.4.1:若三个检验模型中有任意一个模型拒绝原假设时,则t时刻高速公路相邻边缘设备节点交通流量差时间序列是平稳时间序列,满足如下条件:
Figure GDA0003645077240000047
Figure GDA0003645077240000048
Figure GDA0003645077240000051
γ(t,s)=E[(Ztt)(Zss)] (17)
其中,E(Zt)表示为高速公路相邻边缘设备节点交通流量差时间序列Zt的均值,μ表示为常数,Var(Zt)表示为高速公路相邻节点交通流量差时间序列Zt的方差,σ2表示为常数,γ(t,s)表示为自协方差高速公路相邻节点交通流量差时间序列Zt的延迟k自协方差函数,Zs表示为s时刻的高速公路相邻节点交通流量差时间序列,μt和μs分别表示为t和s时刻的高速公路相邻节点交通流量差时间序列均值;
步骤3.4.2:若三个检验模型中所有的检验结果都不拒绝原假设时,则t时刻高速公路相邻边缘设备节点交通流量差时间序列是非平稳时间序列;
步骤4:根据高速公路相邻边缘设备节点交通流量差时间序列的平稳性判定高速公路相邻边缘设备节点的交通异常情况;
若高速公路相邻边缘设备节点交通流量差时间序列为平稳时间序列,则判定高速公路相邻边缘设备节点之间交通正常;
若高速公路相邻边缘设备节点交通流量差时间序列是不平稳时间序列时,则判定高速公路相邻边缘设备节点之间交通可能出现异常,计算交通异常概率,并设置交通异常阈值,通过比较交通异常概率与交通异常阈值大小判定高速公路交通异常情况;
所述交通异常概率采用高速公路相邻节点交通流量差时间序列中的连续递减子序列占原时间序列的百分比表示,如下公式所示:
Figure GDA0003645077240000052
其中,Pi+1,t表示为第i+1个高速公路边缘设备节点在t时刻发生交通异常的概率,
Figure GDA0003645077240000053
表示为高速公路第i+1个和第i个节点交通流量差时间序列中的连续递减子序列元素个数,
Figure GDA0003645077240000054
表示为高速公路第i+1个和第i个节点交通流量差时间序列中元素个数;
若交通异常概率大于设定的交通异常阈值时,则判定当前高速公路相邻边缘设备节点之间出现交通异常,向相邻的边缘设备发送交通异常信息;
若交通异常概率小于设定的交通异常阈值时,则判定当前高速公路相邻边缘设备节点之间交通正常。
采用上述技术方案所产生的有益效果在于:本发明提供的一种高速公路相邻边缘设备节点交通异常概率评估方法,由于高速公路上的任一位置都具有空间和时间可达性,车辆在高速公路上行驶一段路程所用的时间是可以估计的,路程越短行驶时间的估计会越准确,本发明通过对高速公路相邻边缘设备节点的交通流量进行时空约束的分析,得出相邻节点交通流量平均到达时间,通过ADF方法检验相差平均到达时间的相邻边缘设备节点的交通流量差时间序列是否平稳,计算交通异常的概率大小,充分考虑了交通异常和交通流量之间相互影响、相互作用,提高了高速公路交通异常判断的准确性。
附图说明
图1为本发明实施例提供的一种高速公路相邻边缘设备节点交通异常概率评估方法的流程图;
图2为本发明实施例提供的相邻边缘设备节点交通流量时间序列相似距离变化图;
图3为本发明实施例提供的高速公路相邻边缘设备节点交通流量时间序列图;
图4为本发明实施例提供的高速公路相邻边缘设备节点交通流量差时间序列图;
图5为本发明实施例提供的交通正常时高速公路交通流量差时间序列图;
图6为本发明实施例提供的交通异常时高速公路交通流量差时间序列图;
图7为本发明实施例提供的高速公路交通流量差时间序列平稳性检验结果示意图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
本实施例以某两省间的高速公路为例,采用本发明的一种高速公路相邻边缘设备节点交通异常概率评估方法对该高速公路各相邻边缘设备节点间的交通异常情况进行评估。
本实施例中,一种高速公路相邻边缘设备节点交通异常概率评估方法,如图1所示,包括以下步骤:
步骤1:云服务器从高速公路交通流量时间序列库中提取高速公路交通流量计数结果,并基于形态相似性距离计算高速公路相邻边缘设备节点交通流量平均到达时间;
步骤1.1:设定高速公路交通流量平均到达时间表示交通流量序列从第i个边缘设备节点到第i+1个边缘设备节点所用时间,车辆在高速公路上行驶速度为60~120km/h,同时设定相邻两个边缘设备节点的距离为L,则交通流量平均到达时间范围为
Figure GDA0003645077240000061
步骤1.2:计算两个高速公路交通流量序列的形态相似距离,如下公式所示:
Figure GDA0003645077240000071
Figure GDA0003645077240000072
Figure GDA0003645077240000073
其中,D(Xi+1,t,Xi,t-Δt)表示Xi+1,t和Xi,t-Δt之间的形态相似距离,Xi+1,t表示为t时刻第i+1个边缘设备节点的高速公路交通流量序列,Δt表示为
Figure GDA0003645077240000074
时间段内的滑动时间间隔,且
Figure GDA0003645077240000075
Xi,t-Δt表示为t-Δt时刻第i个边缘设备节点高速公路交通流量序列,DED表示为Xi+1,t和Xi,t-Δt之间的欧氏距离,DHD表示为Xi+1,t和Xi,t-Δt的哈夫曼距离,n为高速公路交通流量序列上限;
步骤1.3:初始化两个高速公路交通流量序列最小形态相似距离和两个高速公路交通流量序列最小形态相似距离对应的时刻,如下公式所示:
Figure GDA0003645077240000076
Figure GDA0003645077240000077
其中,dist表示为两个高速公路交通流量序列最小形态相似距离,Tdist_min表示为两个高速公路交通流量序列最小形态相似距离对应的时刻;
步骤1.4:更新高速公路交通流量序列最小形态相似距离和高速公路交通流量序列最小形态相似距离对应的时刻;
步骤1.4.1:当tt时刻在
Figure GDA0003645077240000078
之间时,读取tt时刻第i个边缘设备节点高速公路交通流量序列Xi,tt
步骤1.4.2:计算t时刻第i+1个边缘设备节点的高速公路交通流量序列与tt时刻第i个边缘设备节点高速公路交通流量序列Xi,tt的形态相似距离disti+1,i=D(Xi+1,t,Xi,tt);
步骤1.4.3:根据以下条件更新高速公路交通流量序列最小形态相似距离dist和高速公路交通流量序列最小形态相似距离对应的时刻Tdist_min
当disti+1,i<dist时,令dist=disti+1,i,Tdist_min=tt,然后执行步骤1.5;
当disti+1,i≥dist时,令tt=tt+Δt,重新执行步骤1.4.1;
步骤1.5:基于高速公路交通流量序列最小形态相似距离对应的时刻,得到高速公路交通流量平均到达时间T,如下公式所示:
T=t-Tdist_min (6)
步骤2:生成高速公路相邻边缘设备节点的交通流量差时间序列
Figure GDA0003645077240000081
具体如下公式所示:
Zt=Xi+1,t-Xi,t-T (7)
Figure GDA0003645077240000082
其中,Zt表示为t时刻的高速公路相邻节点交通流量差时间序列,
Figure GDA0003645077240000083
表示t时刻第i+1个边缘设备节点的交通流量时间序列,
Figure GDA0003645077240000084
表示t时刻第i+1个边缘设备节点记录的第n个交通流量值,
Figure GDA0003645077240000085
表示为t-T时刻第i个边缘计算设备的交通流量时间序列,
Figure GDA0003645077240000086
表示t-T时刻第i个边缘设备节点记录的第n个交通流量值;
步骤3:采用基于ADF(Augment Dickey-Fuller)的高速公路相邻边缘设备节点交通流量差时间序列平稳性检验方法,检验相差平均到达时间的两个高速公路相邻边缘设备节点的交通流量差时间序列是否平稳;
步骤3.1:设定高速公路相邻边缘设备节点交通流量差时间序列Zt是一个m阶自回归过程生成的,则Zt表示为:
Figure GDA0003645077240000087
Zt=ρZt-11ΔZt-12ΔZt-2+…+ξmΔZt-m+1t (10)
其中,
Figure GDA0003645077240000088
表示自回归系数,m表示为自回归阶数,εt表示为随机误差项,且误差序列服从均值为0,方差为σ2
Figure GDA0003645077240000089
ΔZt-1,ΔZt-2,…,ΔZt-m+1表示为Zt-1,Zt-2,…,Zt-m+1的一阶差分;
步骤3.2:采用ADF检验常用三种模型对高速公路相邻边缘设备节点交通流量差时间序列Zt进行单位根检验,则t时刻高速公路相邻边缘设备节点交通流量差时间序列的三个检验模型顺序如下:
不含常数项:
Figure GDA0003645077240000091
含常数项:
Figure GDA0003645077240000092
含趋势项:
Figure GDA0003645077240000093
其中,ΔZt表示为Zt的一阶差分,δ表示为常数,βiΔZt-i表示为Zt-i的一阶差分,α表示为常数项系数,βt表示为时间趋势项,m'表示为Zt的滞后阶数;
步骤3.3:设定在平稳性检验过程中,模型假设检验中原假设为H0:δ=0,即高速公路相邻边缘设备节点交通流量差时间序列是非平稳时间序列,备择假设H1:δ<0,即高速公路相邻边缘设备节点交通流量差时间序列是平稳时间序列;
步骤3.4:基于t时刻高速公路相邻边缘设备节点交通流量差时间序列的三个检验模型判断高速公路相邻边缘设备节点交通流量差时间序列是否为平稳时间序列;
步骤3.4.1:若三个检验模型中有任意一个模型拒绝原假设时,则t时刻高速公路相邻边缘设备节点交通流量差时间序列是平稳时间序列,满足如下条件:
Figure GDA0003645077240000094
Figure GDA0003645077240000095
Figure GDA0003645077240000096
γ(t,s)=E[(Ztt)(Zss)] (17)
其中,E(Zt)表示为高速公路相邻边缘设备节点交通流量差时间序列Zt的均值,μ表示为常数,Var(Zt)表示为高速公路相邻节点交通流量差时间序列Zt的方差,σ2表示为常数,γ(t,s)表示为自协方差高速公路相邻节点交通流量差时间序列Zt的延迟k自协方差函数,Zs表示为s时刻的高速公路相邻节点交通流量差时间序列,μt和μs分别表示为t和s时刻的高速公路相邻节点交通流量差时间序列均值;
步骤3.4.2:若三个检验模型中所有的检验结果都不拒绝原假设时,则t时刻高速公路相邻边缘设备节点交通流量差时间序列是非平稳时间序列;
步骤4:根据高速公路相邻边缘设备节点交通流量差时间序列的平稳性判定高速公路相邻边缘设备节点的交通异常情况;
若高速公路相邻边缘设备节点交通流量差时间序列为平稳时间序列,则判定高速公路相邻边缘设备节点之间交通正常;
若高速公路相邻边缘设备节点交通流量差时间序列是不平稳时间序列时,则判定高速公路相邻边缘设备节点之间交通可能出现异常,计算交通异常概率,并设置交通异常阈值,通过比较交通异常概率与交通异常阈值大小判定高速公路交通异常情况;
所述交通异常概率采用高速公路相邻节点交通流量差时间序列中的连续递减子序列占原时间序列的百分比表示,如下公式所示:
Figure GDA0003645077240000101
其中,Pi+1,t表示为第i+1个高速公路边缘设备节点在t时刻发生交通异常的概率,
Figure GDA0003645077240000102
表示为高速公路第i+1个和第i个节点交通流量差时间序列中的连续递减子序列元素个数,
Figure GDA0003645077240000103
表示为高速公路第i+1个和第i个节点交通流量差时间序列中元素个数;
若交通异常概率大于设定的交通异常阈值时,则判定当前高速公路相邻边缘设备节点之间出现交通异常,向相邻的边缘设备发送交通异常信息;
若交通异常概率小于设定的交通异常阈值时,则判定当前高速公路相邻边缘设备节点之间交通正常。
本实施例以python为编程语言,以HPZ820工作站作为云服务器,HPZ820工作站的配置参数如表1所示。
表1 HPZ820服务器配置参数表
Figure GDA0003645077240000104
Figure GDA0003645077240000111
本实施例累积8个小时的某两省间的高速公路交通流量时间序列,边缘设备采集的交通流量是5s的数据,云服务器以12个交通流量计数结果为一组,按每分钟的交通流量计数结果进行分析。
本实施例主要对高速公路交通流量时间序列进行分析,首先从高速公路相邻节点的交通流量时间序列中通过计算相似距离求得交通流量平均到达时间,然后得到高速公路相邻节点交通流量差时间序列,通过交通流量时间差序列的平稳性得到交通异常的概率,具体过程如下:
(1)计算交通流量平均到达时间
滑动窗口选择长度为20,以相邻边缘设备节点连续20分钟的交通流量时间序列,t时刻任意选取的边缘设备节点的交通流量时间序列为7,6,9,10,6,9,8,9,11,6,7,7,4,6,6,8,8,7,6,9,其前一节点的交通流量时间序列为4,5,6,7,9,10,7,8,9,10,10,5,7,8,5,6,6,7,8,7,7,8,5,4,6,10,7,6,9,…,滑动窗口每移动一个长度,计算相邻边缘设备节点交通流量时间序列相似距离,通过最小相似距离找到交通流量平均到达时间,滑动窗口移动10次的实验结果如图2所示。
由图2可得,随着滑动窗口的移动,相邻边缘设备节点高速公路交通流量时间序列的相似距离在不断变化,但是在一段时间内会有一个明显的下降、上升趋势,存在一个最小值。在移动过程中相似距离最小值时,高速公路相邻边缘设备节点的交通流量时间序列最为相似,此时对应的时间即为高速公路交通流量平均到达时间。
(2)高速公路相邻边缘设备交通流量差时间序列
上述结果得到了相邻边缘设备节点交通流量平均到达时间,之后从高速公路交通流量时间序列库中,通过该高速公路相邻边缘设备节点的交通流量对应时间序列得到高速公路交通流量差时间序列,具体如图3、图4所示。
由图3可得,高速公路交通相邻边缘设备节点流量时间序列变化趋势是相似的,在车辆行驶过程中,由于车辆的速度不同,所以边缘设备节点1的高速公路交通流量值在节点2的高速公路交通流量值基础上有微小波动。
由图4可得,高速公路相邻边缘设备节点交通流量差时间序列中的交通流量差的值是围绕0上下波动的,并且波动幅度较小。
(3)高速公路交通异常判断
通过对高速公路相邻边缘设备节点交通流量差时间序列平稳性检验,判断高速公路交通异常概率,ADF平稳性检验的假设是存在单位根,即高速公路交通流量时间序列不平稳。交通正常时高速公路相邻边缘设备节点交通流量差时间序列如图5所示,交通异常时高速公路相邻边缘设备节点交通流量差时间序列如图6所示,对两种情况下的高速公路交通流量差时间序列进行平稳性检验,结果如图7所示。
由图5、图7可得,交通正常情况下,通过数据图示法得出,高速公路相邻边缘设备节点交通流量差时间序列围绕0值上下波动,较平稳。通过ADF检验方法检验时间序列的稳定性,结果拒绝时间序列不平稳的原假设。
由图6、图7可得,交通异常情况下,通过数据图示法得出,一段高速公路交通流量差时间序列在某一时刻前后趋势变化差异大,递减趋势越来越明显。高速公路相邻边缘设备节点之间出现交通异常,导致通过相邻边缘设备节点的交通流量发生明显差异。通过ADF检验方法对该时间序列进行平稳性检验,结果接受时间序列不平稳的原假设。计算得到高速公路交通异常概率为55%,大于设定高速公路交通异常阈值50%。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。

Claims (7)

1.一种高速公路相邻边缘设备节点交通异常概率评估方法,其特征在于:包括以下步骤:
步骤1:云服务器从高速公路交通流量时间序列库中提取高速公路交通流量计数结果,并基于形态相似性距离计算高速公路相邻边缘设备节点交通流量平均到达时间;
步骤2:生成高速公路相邻边缘设备节点的交通流量差时间序列;
步骤3:采用基于ADF的高速公路相邻边缘设备节点交通流量差时间序列平稳性检验方法,检验相差平均到达时间的两个高速公路相邻边缘设备节点的交通流量差时间序列是否平稳;
步骤4:根据高速公路相邻边缘设备节点交通流量差时间序列的平稳性判定高速公路相邻边缘设备节点的交通异常情况;
若高速公路相邻边缘设备节点交通流量差时间序列为平稳时间序列,则判定高速公路相邻边缘设备节点之间交通正常;
若高速公路相邻边缘设备节点交通流量差时间序列是不平稳时间序列时,则判定高速公路相邻边缘设备节点之间交通可能出现异常,计算交通异常概率,并设置交通异常阈值,通过比较交通异常概率与交通异常阈值大小判定高速公路交通异常情况;
若交通异常概率大于设定的交通异常阈值时,则判定当前高速公路相邻边缘设备节点之间出现交通异常,向相邻的边缘设备发送交通异常信息;
若交通异常概率小于设定的交通异常阈值时,则判定当前高速公路相邻边缘设备节点之间交通正常。
2.根据权利要求1所述的一种高速公路相邻边缘设备节点交通异常概率评估方法,其特征在于:所述步骤1的具体方法为:
步骤1.1:设定高速公路交通流量平均到达时间表示交通流量序列从第i个边缘设备节点到第i+1个边缘设备节点所用时间,车辆在高速公路上行驶速度为60~120km/h,同时设定相邻两个边缘设备节点的距离为L,则交通流量平均到达时间范围为
Figure FDA0003645077230000011
步骤1.2:计算两个高速公路交通流量序列的形态相似距离,如下公式所示:
Figure FDA0003645077230000012
Figure FDA0003645077230000013
Figure FDA0003645077230000021
其中,D(Xi+1,t,Xi,t-Δt)表示Xi+1,t和Xi,t-Δt之间的形态相似距离,Xi+1,t表示为t时刻第i+1个边缘设备节点的高速公路交通流量序列,Δt表示为
Figure FDA0003645077230000022
时间段内的滑动时间间隔,且
Figure FDA0003645077230000023
Xi,t-Δt表示为t-Δt时刻第i个边缘设备节点高速公路交通流量序列,DED表示为Xi+1,t和Xi,t-Δt之间的欧氏距离,DHD表示为Xi+1,t和Xi,t-Δt的哈夫曼距离,n为高速公路交通流量序列上限;
步骤1.3:初始化两个高速公路交通流量序列最小形态相似距离和两个高速公路交通流量序列最小形态相似距离对应的时刻,如下公式所示:
Figure FDA0003645077230000024
Figure FDA0003645077230000025
其中,dist表示为两个高速公路交通流量序列最小形态相似距离,Tdist_min表示为两个高速公路交通流量序列最小形态相似距离对应的时刻;
步骤1.4:更新高速公路交通流量序列最小形态相似距离和高速公路交通流量序列最小形态相似距离对应的时刻;
步骤1.5:基于高速公路交通流量序列最小形态相似距离对应的时刻,得到高速公路交通流量平均到达时间T,如下公式所示:
T=t-Tdist_min (6)。
3.根据权利要求2所述的一种高速公路相邻边缘设备节点交通异常概率评估方法,其特征在于:所述步骤1.4的具体方法为:
步骤1.4.1:当tt时刻在
Figure FDA0003645077230000026
之间时,读取tt时刻第i个边缘设备节点高速公路交通流量序列Xi,tt
步骤1.4.2:计算t时刻第i+1个边缘设备节点的高速公路交通流量序列与tt时刻第i个边缘设备节点高速公路交通流量序列Xi,tt的形态相似距离disti+1,i=D(Xi+1,t,Xi,tt);
步骤1.4.3:根据以下条件更新高速公路交通流量序列最小形态相似距离dist和高速公路交通流量序列最小形态相似距离对应的时刻Tdist_min
当disti+1,i<dist时,令dist=disti+1,i,Tdist_min=tt,然后执行步骤1.5;
当disti+1,i≥dist时,令tt=tt+Δt,重新执行步骤1.4.1。
4.根据权利要求3所述的一种高速公路相邻边缘设备节点交通异常概率评估方法,其特征在于:所述步骤2生成的高速公路相邻边缘设备节点的交通流量差时间序列为
Figure FDA0003645077230000031
具体如下公式所示:
Zt=Xi+1,t-Xi,t-T (7)
Figure FDA0003645077230000032
其中,Zt表示为t时刻的高速公路相邻节点交通流量差时间序列,
Figure FDA0003645077230000033
表示t时刻第i+1个边缘设备节点的交通流量时间序列,
Figure FDA0003645077230000034
表示t时刻第i+1个边缘设备节点记录的第n个交通流量值,
Figure FDA0003645077230000035
表示为t-T时刻第i个边缘计算设备的交通流量时间序列,
Figure FDA0003645077230000036
表示t-T时刻第i个边缘设备节点记录的第n个交通流量值。
5.根据权利要求4所述的一种高速公路相邻边缘设备节点交通异常概率评估方法,其特征在于:所述步骤3的具体方法为:
步骤3.1:设定高速公路相邻边缘设备节点交通流量差时间序列Zt是一个m阶自回归过程生成的,则Zt表示为:
Figure FDA0003645077230000037
Zt=ρZt-11ΔZt-12ΔZt-2+…+ξmΔZt-m+1t (10)
其中,
Figure FDA0003645077230000038
表示自回归系数,m表示为自回归阶数,εt表示为随机误差项,且误差序列服从均值为0,方差为σ2
Figure FDA0003645077230000039
ΔZt-1,ΔZt-2,...,ΔZt-m+1表示为Zt-1,Zt-2,...,Zt-m+1的一阶差分;
步骤3.2:采用ADF检验的三种模型对高速公路相邻边缘设备节点交通流量差时间序列Zt进行单位根检验,则t时刻高速公路相邻边缘设备节点交通流量差时间序列的三个检验模型顺序如下:
不含常数项:
Figure FDA0003645077230000041
含常数项:
Figure FDA0003645077230000042
含趋势项:
Figure FDA0003645077230000043
其中,ΔZt表示为Zt的一阶差分,δ表示为常数,βiΔZt-i表示为Zt-i的一阶差分,α表示为常数项系数,βt表示为时间趋势项,m'表示为Zt的滞后阶数;
步骤3.3:设定在平稳性检验过程中,模型假设检验中原假设为H0:δ=0,即高速公路相邻边缘设备节点交通流量差时间序列是非平稳时间序列,备择假设H1:δ<0,即高速公路相邻边缘设备节点交通流量差时间序列是平稳时间序列;
步骤3.4:基于t时刻高速公路相邻边缘设备节点交通流量差时间序列的三个检验模型判断高速公路相邻边缘设备节点交通流量差时间序列是否为平稳时间序列。
6.根据权利要求5所述的一种高速公路相邻边缘设备节点交通异常概率评估方法,其特征在于:所述步骤3.4的具体方法为:
步骤3.4.1:若三个检验模型中有任意一个模型拒绝原假设时,则t时刻高速公路相邻边缘设备节点交通流量差时间序列是平稳时间序列,满足如下条件:
Figure FDA0003645077230000044
Figure FDA0003645077230000045
Figure FDA0003645077230000046
γ(t,s)=E[(Ztt)(Zss)] (17)
其中,E(Zt)表示为高速公路相邻边缘设备节点交通流量差时间序列Zt的均值,μ表示为常数,Var(Zt)表示为高速公路相邻节点交通流量差时间序列Zt的方差,σ2表示为常数,γ(t,s)表示为自协方差高速公路相邻节点交通流量差时间序列Zt的延迟k自协方差函数,Zs表示为s时刻的高速公路相邻节点交通流量差时间序列,μt和μs分别表示为t和s时刻的高速公路相邻节点交通流量差时间序列均值;
步骤3.4.2:若三个检验模型中所有的检验结果都不拒绝原假设时,则t时刻高速公路相邻边缘设备节点交通流量差时间序列是非平稳时间序列。
7.根据权利要求6所述的一种高速公路相邻边缘设备节点交通异常概率评估方法,其特征在于:步骤4所述交通异常概率采用高速公路相邻节点交通流量差时间序列中的连续递减子序列占原时间序列的百分比表示,如下公式所示:
Figure FDA0003645077230000051
其中,Pi+1,t表示为第i+1个高速公路边缘设备节点在t时刻发生交通异常的概率,
Figure FDA0003645077230000052
表示为高速公路第i+1个和第i个节点交通流量差时间序列中的连续递减子序列元素个数,
Figure FDA0003645077230000053
表示为高速公路第i+1个和第i个节点交通流量差时间序列中元素个数。
CN202110995620.3A 2021-08-27 2021-08-27 一种高速公路相邻边缘设备节点交通异常概率评估方法 Active CN113706867B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110995620.3A CN113706867B (zh) 2021-08-27 2021-08-27 一种高速公路相邻边缘设备节点交通异常概率评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110995620.3A CN113706867B (zh) 2021-08-27 2021-08-27 一种高速公路相邻边缘设备节点交通异常概率评估方法

Publications (2)

Publication Number Publication Date
CN113706867A CN113706867A (zh) 2021-11-26
CN113706867B true CN113706867B (zh) 2022-08-30

Family

ID=78655990

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110995620.3A Active CN113706867B (zh) 2021-08-27 2021-08-27 一种高速公路相邻边缘设备节点交通异常概率评估方法

Country Status (1)

Country Link
CN (1) CN113706867B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107798871A (zh) * 2017-10-27 2018-03-13 云南大学 一种高速公路收费站交通流量预测方法及系统
CN110765980A (zh) * 2019-11-05 2020-02-07 中国人民解放军国防科技大学 一种异常驾驶的检测方法及装置
CN110956387A (zh) * 2019-11-27 2020-04-03 西南交通大学 基于var模型的交通安全与经济损失关系计算方法
CN112150801A (zh) * 2020-08-21 2020-12-29 同济大学 一种基于准车辆轨迹的快速路实时事故风险预测方法
CN112529299A (zh) * 2020-12-11 2021-03-19 东南大学 基于arima和lstm混合神经网络的短交通流预测方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110320500A1 (en) * 2010-06-28 2011-12-29 SRM Institute of Science and Technology System and method for constructing a university model graph
IT1402100B1 (it) * 2010-09-14 2013-08-28 Le Gemme S P A Sistema per rilevare e segnalare condizioni particolari in una carreggiata stradale, tipicamente condizioni di pericolo
US9134928B2 (en) * 2011-01-28 2015-09-15 Ricoh Production Print Solutions LLC Print job processing in an automated document factory environment
CN103337161A (zh) * 2013-07-11 2013-10-02 上海济安交通工程咨询有限公司 基于实时仿真模型的交叉口动态综合评价及信控系统优化方法
US9455921B2 (en) * 2014-02-05 2016-09-27 Verizon Patent And Licensing Inc. End to end network management based on quality of service adjustments
CN103903452B (zh) * 2014-03-11 2015-10-07 东南大学 交通流短时预测方法
CN104517159A (zh) * 2014-12-18 2015-04-15 上海交通大学 一种公交短时客流的预测方法
US9558601B2 (en) * 2015-04-24 2017-01-31 The Boeing Company System and method for detecting vehicle system faults
CN105469611B (zh) * 2015-12-24 2017-12-01 大连理工大学 一种短时交通流预测模型方法
CN110991842A (zh) * 2019-11-22 2020-04-10 广东省交通集团有限公司 高速公路货物运输量和地区gdp耦合关系分析的方法
CN112330158B (zh) * 2020-11-06 2023-10-13 北京建筑大学 基于自回归差分移动平均-卷积神经网络的识别交通指数时间序列的方法
CN112669599B (zh) * 2020-12-14 2024-10-25 成都易书桥科技有限公司 一种基于时序分析和残差匹配的短时交通流预测方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107798871A (zh) * 2017-10-27 2018-03-13 云南大学 一种高速公路收费站交通流量预测方法及系统
CN110765980A (zh) * 2019-11-05 2020-02-07 中国人民解放军国防科技大学 一种异常驾驶的检测方法及装置
CN110956387A (zh) * 2019-11-27 2020-04-03 西南交通大学 基于var模型的交通安全与经济损失关系计算方法
CN112150801A (zh) * 2020-08-21 2020-12-29 同济大学 一种基于准车辆轨迹的快速路实时事故风险预测方法
CN112529299A (zh) * 2020-12-11 2021-03-19 东南大学 基于arima和lstm混合神经网络的短交通流预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
一种利用不完整数据检测交通异常的方法;王玉玲等;《计算机科学》;20160630;第43卷(第6A期);第425-429页 *

Also Published As

Publication number Publication date
CN113706867A (zh) 2021-11-26

Similar Documents

Publication Publication Date Title
WO2020244288A1 (zh) 一种基于gps轨迹数据的大货车驾驶行为评估方法及装置
CN110077398B (zh) 一种用于智能驾驶的危险处理方法
CN114783183A (zh) 一种基于交通态势算法的监控方法与系统
CN1862620A (zh) 高速公路汽车行驶智能检测预警方法及其预警系统
CN110766982B (zh) 基于车载传感器的车辆碰撞检测系统
CN103544753B (zh) 一种道闸控制方法及系统
CN104809878A (zh) 利用公交车gps数据检测城市道路交通异常状态的方法
CN111243338A (zh) 基于车辆加速度的碰撞风险评价方法
CN113127466B (zh) 一种车辆轨迹数据预处理方法及计算机存储介质
CN105160326A (zh) 一种高速公路上停车的自动检测方法与装置
CN111862598B (zh) 一种基于高清卡口数据及事故风险的可变限速控制方法
CN111524365A (zh) 一种利用多地磁传感器进行车型分类的方法
Li et al. Large car-following data based on lyft level-5 open dataset: Following autonomous vehicles vs. human-driven vehicles
CN116246229B (zh) 一种基于车联网的骑行状态监控系统及方法
CN113313357A (zh) 一种基于高斯过程回归分析的交通道路安全评价方法
Djukanović et al. A dataset for audio-video based vehicle speed estimation
CN113706867B (zh) 一种高速公路相邻边缘设备节点交通异常概率评估方法
CN110986994B (zh) 一种基于高噪音车辆轨迹数据的换道意图自动标注方法
CN117953693B (zh) 一种基于视频ai的高速公路交通智能监测系统
CN117994983B (zh) 一种基于物联网的高速路网态势感知方法及系统
CN116631187B (zh) 一种案件现场勘查信息智能采集分析系统
CN113706889B (zh) 基于目标检测分析的高速公路团雾测量系统及方法
CN112365721B (zh) 一种车辆异常行为识别、预警方法
CN117037081A (zh) 基于机器学习的交通监控方法、装置、设备及介质
CN105303833B (zh) 基于微波车辆检测器的高架桥突发事件判别方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant