CN115376316B - 一种基于隐式马尔可夫链的危险跟驰行为识别方法 - Google Patents
一种基于隐式马尔可夫链的危险跟驰行为识别方法 Download PDFInfo
- Publication number
- CN115376316B CN115376316B CN202210997427.8A CN202210997427A CN115376316B CN 115376316 B CN115376316 B CN 115376316B CN 202210997427 A CN202210997427 A CN 202210997427A CN 115376316 B CN115376316 B CN 115376316B
- Authority
- CN
- China
- Prior art keywords
- vehicle
- dangerous
- state
- following
- time
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000013139 quantization Methods 0.000 claims abstract description 67
- 230000006399 behavior Effects 0.000 claims abstract description 54
- 230000003993 interaction Effects 0.000 claims abstract description 8
- 230000001133 acceleration Effects 0.000 claims description 33
- 239000011159 matrix material Substances 0.000 claims description 16
- 239000013598 vector Substances 0.000 claims description 13
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 230000003321 amplification Effects 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 5
- 239000012634 fragment Substances 0.000 claims description 4
- 238000011002 quantification Methods 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000007637 random forest analysis Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- 238000007619 statistical method Methods 0.000 abstract 1
- 238000011160 research Methods 0.000 description 5
- 206010039203 Road traffic accident Diseases 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000012549 training Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000008447 perception Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 231100001267 hazard identification Toxicity 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/17—Terrestrial scenes taken from planes or by drones
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/182—Network patterns, e.g. roads or rivers
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/017—Detecting movement of traffic to be counted or controlled identifying vehicles
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/065—Traffic 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明提出了一种基于隐式马尔可夫链的危险跟驰行为识别方法,首先通过无人机以及路侧协同设备获取路段车辆信息,建立考虑车辆交互行为和宏观交通流特征的多维数据集;然后,构建时空综合维度的危险量化指标DMI,通过对各危险量化指标的局限性及相关性分析,基于熵的思想与数学统计方法,从时空综合维度对跟驰行为进危险量化;最后根据量化标定结果与隐式马尔可夫链理论,对车辆危险跟驰这类不可观测的驾驶状态进行识别,一定程度上对危险跟驰行为进行提前预警,提高行车安全性。
Description
技术领域
本发明涉及交通安全、危险识别领域,具体涉及一种基于隐式马尔可夫链的危险跟驰行为识别方法。
背景技术
随着我国城市化水平的进一步提高和道路***的不断完善,汽车已成为人们必不可少的交通工具。根据历史数据统计,我国每年有超过10万人死于交通事故,其中人为因素导致的交通事故占到总事故数的80%,其外在表现就是“危险驾驶行为”。因此,对行驶过程中潜在的危险行为进行识别,对道路交通安全和运行效率的提升具有现实意义和社会效益。
目前,针对危险驾驶行为识别的研究,主要分为三个方面:一方面,研究者们对驾驶行为的表征参数进行了研究,驾驶人生理参数获取方式为侵入式,在实时性和应用性上有所欠缺,车辆动力学参数的研究大多只针对自车信息进行提取,未能充分考虑潜在危险行为附近的微观交通流条件以及车辆与周围更多车辆的交互作用;另一方面,专家学者对危险行为的标定方法进行了讨论,研究多使用单一安全替代指标进行描述,并以固定阈值进行危险判定,忽略了车辆动态特征以及车辆在空间和时间两个维度的潜在危险;最后,机器学习算法常被用于车辆换道、转弯等驾驶行为的识别,意图识别模型比较丰富,对危险行为的识别研究较少。
发明内容
本发明的技术解决问题:解决以上现有的危险行为识别方法所存在的问题,提供一种基于隐式马尔可夫链的危险跟驰行为识别方法,大大提高了驾驶安全性,有效避免交通事故的产生。且本发明的识别方法与驾驶人实际的跟驰行为和危险感知更为相符,有着更强的理论依据和更好的普适性。
本发明的技术方案:一种基于隐式马尔可夫链的危险跟驰行为识别方法,其特征在于,包括以下步骤:
步骤1:采集路段车辆信息,提取车辆跟驰行为片段,建立多维度数据集S,多维度数据集S包括描述自车与周围车辆交互行为的数据集Scfw以及描述路段交通流特征的数据集Scfh;所述数据集Scfw共包含自车速度、前车速度、自车加速度、前车加速度、车头间距、车头时距、前后两车相对速度和加速度8个特征;数据集Scfh共包含路段交通量、车流密度、路段车辆时间平均车速、平均加速度、速度标准差以及加速度标准差6个特征;
Scfw={v,vp,a,ap,SH,TH,Δv,Δa}
其中,v为自车速度;vp为前车速度;a为自车加速度;ap为前车加速度;SH为车头间距;TH为车头时距;Δv为自车与前车相对速度;Δa为自车与前车相对加速度;Q为路段交通量;K为路段车流密度;为路段车辆时间平均车速;/>为路段车辆平均加速度;σv为路段车速标准差;σa为路段加速度标准差;
步骤2:基于步骤1中的数据集Scfw,分别计算车辆的时间危险量化指标、距离危险量化指标和减速度危险量化指标;所述时间危险量化指标包括碰撞时间TTC和改进碰撞时间MTTC;距离危险量化指标包括可接受最小停车距离MSD和紧急减速碰撞指数PICUD;所述减速度危险量化指标是指避免碰撞减速度DRAC;使用各指标熵值和模糊规则,对这五种危险量化指标进行融合,得到时空综合维度的危险量化指标DMI;
步骤3:利用步骤2的时空综合维度的危险量化指标DMI,计算得到不同道路交通服务水平下的车辆每一时刻危险量化值的累积分布曲线,当累积分布<90%时为低风险,DMI累积分布≥90%时为高风险,并根据车辆每一时刻的风险状态标定得到观测序列内车辆的跟驰状态;
步骤4:对步骤1中的多维度数据集S中的各个特征进行统计学扩增,计算各个特征的最大值、最小值、平均值和差分值,构成新的数据集S′,S′包括特征扩增后的自车与周围车辆交互行为的数据集S′cfw以及特征扩增后的路段交通流特征的数据集S′cfh;采用随机森林算法分别对S′cfw与S′cfh数据集中影响危险跟驰行为识别结果的特征进行重要度排序及筛选,得到危险跟驰行为识别模型的输入F;
步骤5:利用步骤3得到的车辆跟驰状态和步骤4的输入F,结合高斯函数与隐马尔可夫模型,构建并训练高斯混合隐马尔可夫模型GMM-HMM,由此实现危险跟驰行为的识别。
进一步,所述步骤2中,对五种危险量化指标进行融合,得到时空综合维度的危险量化指标DMI的步骤如下:
a.构建原始评价矩阵D=(dij)n×m,其中,n为车辆跟驰片段中的总时间戳个数;m为危险量化指标的个数;dij表示在车辆跟驰片段中第i个时间戳对应的第j个危险量化指标,共有TTC、MTTC、MSD、PICUD、DRAC五种危险量化指标;
b.对五种危险量化指标进行归一化处理;
c.计算第j个危险量化指标的熵值ej;
d.计算第j个危险量化指标的权重ωj;
e.基于历史数据相关性及模糊规则获得独立性权重δj,最终确定各指标的综合权重Gj,计算得到时空综合维度的危险量化指标DMI;
其中,δj为第j个危险量化指标基于模糊规则获得的独立性权重;ωj为第j个危险量化指标基于熵值法获得的客观权重;为各指标客观权重的平均值;Gj为各危险量化指标的最终权重;DMIi为车辆跟驰片段中第i个时间戳的危险量化值。
进一步,所述步骤3中,对车辆跟驰状态的标定方法具体为:首先,以DMI累积分布90%为车辆风险状态阈值,得到交通服务水平分别在一、二、三、四等级下的车辆风险状态P;其次,以1个连续时间戳作为一个车辆跟驰行为的观测序列,T为观测序列长度,该观测序列车辆跟驰状态G的标定方法如下:
其中,Pi为i时刻车辆的风险状态;DMIthreshold,k为交通服务水平为k时危险量化指标的阈值;Δt为该帧的时长;n为该观测序列的时间戳总个数即总帧数;T为观测序列的长度。
进一步,所述步骤5中,所述高斯混合隐马尔可夫模型GMM-HMM模型为:
λG={π,A,C,μ,U}
π=(πi)
A=[aij]N×N
aij=P(qt+1=Sj|qt=Si)
C=[cjg]N×G
μ=(μj)
U=(Uj)
其中,λG为模型参数集合;
π为起始状态概率向量,πi表示初始时刻处于状态Si的概率,状态变量分别为危险跟驰状态和安全跟驰状态,π是一个二维向量;
A为状态转移概率矩阵,aij表示qt时刻隐状态为Si的情况下,qt+1时刻状态为Sj的概率,N为状态变量个数,N=2,所以A是一个2×2的矩阵;
C为危险跟驰状态和安全跟驰状态下混合高斯模型的分布权重系数矩阵,cjg表示第j个隐状态下第g个单高斯模型的权重,G为隐状态的时高斯分量的数目;
μ是M个混合成分的均值向量的集合,μj表示的是j个隐状态的特征变量均值向量;
U是N×G个混合成分的协方差矩阵的集合,Uj表示的是第j个隐状态的特征变量协方差矩阵。
本发明与现有技术相比的优点在于:
本发明的方法在考虑车辆动力学和交通流参数的基础上,建立多维度驾驶人跟驰行为集;基于熵值法与数学统计方法,从时空综合维度对车辆跟驰行为进行动态风险评估和危险量化;并根据量化标定结果与隐式马尔可夫链理论,对车辆危险跟驰这类不可观测的驾驶状态进行识别,在很大程度上对危险跟驰行为进行提前预警,大大提高驾驶安全性,有效避免交通事故的产生。且本发明的方法与驾驶人实际的跟驰行为和危险感知更为相符,有着更强的理论依据和更好的普适性。
附图说明
图1为本发明的整体思路框图;
图2为危险跟驰行为识别过程图。
具体实施方式
下面结合附图和实施方案对本发明进行详细的描述。应理解该实例仅用于说明本发明而不用于限制本发明的范围。
如图1所示,本发明的一种基于隐式马尔可夫链的危险跟驰行为识别方法,具体包括以下步骤:
步骤1:采集路段车辆信息,提取车辆跟驰行为片段,建立多维度数据集S。
(1)车辆跟驰行为是指在连续的交通流中的车辆不进行换道操作,而是跟驰前车行驶的行为,后车的驾驶行为主要受到前车的影响,因此对于跟驰状态的提取多基于车头时距、车间距及持续时间等指标。为了避免选取的驾驶行为片段车辆处于自由流行驶状态或换道准备阶段,本发明对跟驰行为片段的提取标准做以下定义:车头间距SH<120m,持续时T>10s。
(2)根据采集得到的路段车辆信息,提取车辆跟驰行为片段,建立多维度数据集S,多维度数据集S包括描述自车与周围车辆交互行为的数据集Scfw以及描述路段交通流特征的数据集Scfh;所述数据集Scfw共包含自车速度、前车速度、自车加速度、前车加速度、车头间距、车头时距、前后两车相对速度和加速度8个特征;数据集Scfh共包含路段交通量、车流密度、路段车辆时间平均车速、平均加速度、速度标准差以及加速度标准差6个特征;
Scfw={v,vp,a,ap,SH,TH,Δv,Δa}
其中,v为自车速度;vp为前车速度;a为自车加速度;ap为前车加速度;SH为车头间距;TH为车头时距;Δv为自车与前车相对速度;Δa为自车与前车相对加速度;Q为路段交通量;K为路段车流密度;为路段车辆时间平均车速;/>为路段车辆平均加速度;σv为路段车速标准差;σa为路段加速度标准差;
数据集中特征v,vp,a,ap,SH可以通过无人机视频数据直接获得,其他特征的计算公式如下:
Δv=v-vp
Δa=a-ap
K=u/L
其中,u为无人机拍摄路段第i个时间戳车辆数目;vr为第i个时间戳路段上第r辆车的速度;ar为第i个时间戳路段上第r辆车的加速度;L为无人机拍摄路段长度。
步骤2:基于步骤1中的数据集Scfw,分别计算车辆的时间危险量化指标、距离危险量化指标和减速度危险量化指标;所述时间危险量化指标包括碰撞时间TTC和改进碰撞时间MTTC;距离危险量化指标包括可接受最小停车距离MSD和紧急减速碰撞指数PICUD;所述减速度危险量化指标是指避免碰撞减速度DRAC;使用各指标熵值和模糊规则,对这五种危险量化指标进行融合,得到时空综合维度的危险量化指标DMI;
(1)计算五种危险量化指标
a.构建原始评价矩阵D=(dij)n×m,其中,n为车辆跟驰片段中的总时间戳个数;m为危险量化指标的个数;dij表示在车辆跟驰片段中第i个时间戳对应的第j个危险量化指标,共有TTC、MTTC、MSD、PICUD、DRAC五种危险量化指标;
其中,amax为车辆最大减速度,一般取5m/s2;
b.对数据进行归一化处理;
Index(Ij)={TTC,MTTC,MSD,PICUD,DRAC} j=1,2,3,4,5
(2)建立时空综合维度的危险量化指标DMI
a.计算第j个危险量化指标的熵值ej;
其中,pij为第j个危险指标中第i个时间戳所占的比重;ej为第j个危险指标的熵值;
b.计算第j个危险量化指标的权重ωj;
c.基于历史数据相关性及模糊规则获得独立性权重δj,最终确定各指标的综合权重Gj,计算得到时空综合维度的危险量化指标DMI;
其中,δj为第j个危险量化指标基于模糊规则获得的独立性权重;ωj为第j个危险量化指标基于熵值法获得的客观权重;为各指标客观权重的平均值;Gj为各危险量化指标的最终权重;DMIi为车辆跟驰片段中第i个时间戳的危险量化值。
步骤3:利用步骤2的时空综合维度的危险量化指标DMI,计算得到不同道路交通服务水平下的车辆每一时刻危险量化值的累积分布曲线,当累积分布<90%时为低风险,DMI累积分布≥90%时为高风险,并根据车辆每一时刻的风险状态标定得到观测序列内车辆的跟驰状态;
根据中国高速公路服务水平划分标准,可基于宏观跟驰行为集Scfh确定自车所处的交通环境,分别获取服务水平为一、二、三、四级下车辆每一时刻危险量化值的累积分布曲线,当累积分布<90%时为低风险,DMI累积分布≥90%时为高风险。其次,以l个连续时间戳作为一个车辆跟驰行为的观测序列,T为观测序列长度,该观测序列车辆跟驰状态G的标定方法如下:
其中,Pi为i时刻车辆的风险状态;DMIthreshold,k为交通服务水平为k时危险量化指标的阈值;Δt为该帧的时长;n为该观测序列的总帧数;T为观测序列的长度。
最终可获得:一级服务水平下安全与危险跟驰状态、二级服务水平下安全与危险跟驰状态、三级服务水平下安全与危险跟驰状态、四级服务水平下安全与危险跟驰状态,共计8种跟驰状态。
步骤4:对步骤1中的多维度数据集S中的各个特征进行统计学扩增,计算各个特征的最大值、最小值、平均值和差分值,构成新的数据集S′,S′包括特征扩增后的自车与周围车辆交互行为的数据集S′cfw以及特征扩增后的路段交通流特征的数据集S′cfh;采用随机森林算法分别对S′cfw与S′cfh数据集中影响危险跟驰行为识别结果的特征进行重要度排序及筛选,得到危险跟驰行为识别模型的输入F;
步骤5:利用步骤3得到的车辆跟驰状态和步骤4的输入F,结合高斯函数与隐马尔可夫模型,构建并训练高斯混合隐马尔可夫模型GMM-HMM,由此实现危险跟驰行为的识别。
(1)模型建立,使用混合高斯函数来拟合连续性的观测值,混合高斯概率密度的形式为:
其中,bj(O)为隐藏状态j下第g个单高斯模型观察到的观测向量为O的概率;cjg表示混合权重,即隐状态j下第g个单高斯模型的权重,同时也满足G为隐状态的时高斯分量的数目;P(O,μjg,Ujg)表示的是高斯概率密度函数;μjg表示的是隐状态的第g个特征变量的均值向量;Ujg表示的是隐状态的第g个特征变量的协方差矩阵;N为隐藏状态个数。
GMM-HMM的模型参数λG可以表示为:
λG={π,A,C,μ,U}
π=(πi)
A=[aij]N×N
aij=P(qt+1=Sj|qt=Si)
C=[cjg]N×G
μ=(μj)
U=(Uj)
其中,π为起始状态概率向量,πi表示初始时刻处于状态Si的概率,状态变量分别为危险跟驰状态和安全跟驰状态,π是一个二维向量;A为状态转移概率矩阵,aij表示qt时刻隐状态为Si的情况下,qt+1时刻状态为Sj的概率,N=2,所以A是一个2×2的矩阵;C为危险跟驰状态和安全跟驰状态下混合高斯模型的分布权重系数矩阵,μ是M个特征变量的均值向量的集合;U是N×G个特征变量的协方差矩阵的集合。
(2)模型求解,实现跟驰状态识别。使用前向-后向算法降低计算复杂度并利用Baum-Welch算法进行含隐状态的参数学***下的模型训练、识别为例,具体过程如下:
a.确定初始参数:本发明中隐状态为安全跟驰状态和危险跟驰状态,所以N取值为2。由于单高斯分布已经可以很好地描述大多数情况下的概率分布情况,因此将两个隐马尔可夫模型的M值均设为1,相应地,高斯分布权重系数C也设为1。一般情况下,π和A的初值选取对模型的影响不大,采用随机数生成的方法,分别给两个模型的π和A设定初值。单高斯分布的均值μ和协方差矩阵U的初值通过抽取数据集中的部分样本估计得到;
b.取训练集中的观测序列输入GMM-HMM模型,利用前向-后向算法,计算出前向概率和后向概率;
c.通过Baum-Welch算法进行模型参数重估,更新模型参数λG;
d.再次利用前向-后向算法,计算隐马尔可夫模型输出的似然概率,若概率增加,则需要重复步骤b和c,直到参数逐渐收敛不变,保存安全跟驰的识别模型参数λG1、危险跟驰的识别模型参数λG2。
经过反复迭代,最终获得危险跟驰以及安全跟驰两个高斯混合隐马尔可夫模型的似然概率,通过对比获得识别结果,由此实现危险跟驰行为的准确识别。
以上虽然描述了本发明的具体实施方法,但是本领域的技术人员应当理解,这些仅是举例说明,在不背离本发明原理和实现的前提下,可以对这些实施方案做出多种变更或修改,因此,本发明的保护范围由所附权利要求书限定。
Claims (2)
1.一种基于隐式马尔可夫链的危险跟驰行为识别方法,其特征在于,包括以下步骤:
步骤1:采集路段车辆信息,提取车辆跟驰行为片段,建立多维度数据集S,多维度数据集S包括描述自车与周围车辆交互行为的数据集Scfw以及描述路段交通流特征的数据集Scfh;所述数据集Scfw共包含自车速度、前车速度、自车加速度、前车加速度、车头间距、车头时距、前后两车相对速度和加速度8个特征;数据集Scfh共包含路段交通量、车流密度、路段车辆时间平均车速、平均加速度、速度标准差以及加速度标准差6个特征;
Scfw={v,vp,a,ap,SH,TH,Δv,Δa}
其中,v为自车速度;vp为前车速度;a为自车加速度;ap为前车加速度;SH为车头间距;TH为车头时距;Δv为自车与前车相对速度;Δa为自车与前车相对加速度;Q为路段交通量;K为路段车流密度;为路段车辆时间平均车速;/>为路段车辆平均加速度;σv为路段车速标准差;σa为路段加速度标准差;
步骤2:基于步骤1中的数据集Scfw,分别计算车辆的时间危险量化指标、距离危险量化指标和减速度危险量化指标;所述时间危险量化指标包括碰撞时间TTC和改进碰撞时间MTTC;距离危险量化指标包括可接受最小停车距离MSD和紧急减速碰撞指数PICUD;所述减速度危险量化指标是指避免碰撞减速度DRAC;使用各指标熵值和模糊规则,对这五种危险量化指标进行融合,得到时空综合维度的危险量化指标DMI;
步骤3:利用步骤2的时空综合维度的危险量化指标DMI,计算得到不同道路交通服务水平下的车辆每一时刻危险量化值的累积分布曲线,当累积分布<90%时为低风险,DMI累积分布≥90%时为高风险,并根据车辆每一时刻的风险状态标定得到观测序列内车辆的跟驰状态;
步骤4:对步骤1中的多维度数据集S中的各个特征进行统计学扩增,计算各个特征的最大值、最小值、平均值和差分值,构成新的数据集S′,S′包括特征扩增后的自车与周围车辆交互行为的数据集S′cfw以及特征扩增后的路段交通流特征的数据集S′cfh;采用随机森林算法分别对S′cfw与S′cfh数据集中影响危险跟驰行为识别结果的特征进行重要度排序及筛选,得到危险跟驰行为识别模型的输入F;
步骤5:利用步骤3得到的车辆跟驰状态和步骤4的输入F,结合高斯函数与隐马尔可夫模型,构建并训练高斯混合隐马尔可夫模型GMM-HMM,由此实现危险跟驰行为的识别;
所述步骤2中,对五种危险量化指标进行融合,得到时空综合维度的危险量化指标DMI的步骤如下:
a.构建原始评价矩阵D=(dij)n×m,其中,n为车辆跟驰片段中的总时间戳个数;m为危险量化指标的个数;dij表示在车辆跟驰片段中第i个时间戳对应的第j个危险量化指标,共有TTC、MTTC、MSD、PICUD、DRAC五种危险量化指标;
b.对五种危险量化指标进行归一化处理;
c.计算第j个危险量化指标的熵值ej;
d.计算第j个危险量化指标的权重ωj;
e.基于历史数据相关性及模糊规则获得独立性权重δj,最终确定各指标的综合权重Gj,计算得到时空综合维度的危险量化指标DMI;
其中,δj为第j个危险量化指标基于模糊规则获得的独立性权重;ωj为第j个危险量化指标基于熵值法获得的客观权重;为各指标客观权重的平均值;Gj为各危险量化指标的最终权重;DMIi为车辆跟驰片段中第i个时间戳的危险量化值;
所述步骤5中,所述高斯混合隐马尔可夫模型GMM-HMM模型为:
λG={π,A,C,μ,U}
π=(πi)
A=[aij]N×N
aij=P(qt+1=Sj|qt=Si)
C=[cjg]N×G
μ=(μj)
U=(Uj)
其中,λG为模型参数集合;
π为起始状态概率向量,πi表示初始时刻处于状态Si的概率,状态变量分别为危险跟驰状态和安全跟驰状态,π是一个二维向量;
A为状态转移概率矩阵,aij表示qt时刻隐状态为Si的情况下,qt+1时刻状态为Sj的概率,N为状态变量个数,N=2,所以A是一个2×2的矩阵;
C为危险跟驰状态和安全跟驰状态下混合高斯模型的分布权重系数矩阵,cjg表示第j个隐状态下第g个单高斯模型的权重,G为隐状态的时高斯分量的数目;
μ是M个混合成分的均值向量的集合,μj表示的是j个隐状态的特征变量均值向量;
U是N×G个混合成分的协方差矩阵的集合,Uj表示的是第j个隐状态的特征变量协方差矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210997427.8A CN115376316B (zh) | 2022-08-19 | 2022-08-19 | 一种基于隐式马尔可夫链的危险跟驰行为识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210997427.8A CN115376316B (zh) | 2022-08-19 | 2022-08-19 | 一种基于隐式马尔可夫链的危险跟驰行为识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115376316A CN115376316A (zh) | 2022-11-22 |
CN115376316B true CN115376316B (zh) | 2023-05-30 |
Family
ID=84065274
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210997427.8A Active CN115376316B (zh) | 2022-08-19 | 2022-08-19 | 一种基于隐式马尔可夫链的危险跟驰行为识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115376316B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008127465A1 (en) * | 2007-04-11 | 2008-10-23 | Nec Laboratories America, Inc. | Real-time driving danger level prediction |
EP3035251A2 (en) * | 2014-12-15 | 2016-06-22 | Siemens Aktiengesellschaft | Markov chain generation apparatus |
CN107958269A (zh) * | 2017-11-28 | 2018-04-24 | 江苏大学 | 一种基于隐马尔可夫模型的行车危险度预测方法 |
CN113611157A (zh) * | 2021-08-10 | 2021-11-05 | 北京航空航天大学合肥创新研究院(北京航空航天大学合肥研究生院) | 一种高速公路车辆追尾碰撞风险估计方法 |
CN113761715A (zh) * | 2021-08-11 | 2021-12-07 | 江苏大学 | 一种基于高斯混合和隐马尔可夫的个性化车辆跟驰模型的建立方法 |
-
2022
- 2022-08-19 CN CN202210997427.8A patent/CN115376316B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008127465A1 (en) * | 2007-04-11 | 2008-10-23 | Nec Laboratories America, Inc. | Real-time driving danger level prediction |
EP3035251A2 (en) * | 2014-12-15 | 2016-06-22 | Siemens Aktiengesellschaft | Markov chain generation apparatus |
CN107958269A (zh) * | 2017-11-28 | 2018-04-24 | 江苏大学 | 一种基于隐马尔可夫模型的行车危险度预测方法 |
CN113611157A (zh) * | 2021-08-10 | 2021-11-05 | 北京航空航天大学合肥创新研究院(北京航空航天大学合肥研究生院) | 一种高速公路车辆追尾碰撞风险估计方法 |
CN113761715A (zh) * | 2021-08-11 | 2021-12-07 | 江苏大学 | 一种基于高斯混合和隐马尔可夫的个性化车辆跟驰模型的建立方法 |
Non-Patent Citations (1)
Title |
---|
基于隐马尔可夫模型的驾驶行为预测方法研究;肖献强;任春燕;王其东;;中国机械工程(第21期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115376316A (zh) | 2022-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107958269B (zh) | 一种基于隐马尔可夫模型的行车危险度预测方法 | |
CN110949398B (zh) | 一种车辆编队行驶中头车驾驶员异常驾驶行为检测方法 | |
CN109733390B (zh) | 一种基于驾驶人特性的自适应换道预警方法 | |
CN110569783B (zh) | 一种驾驶人换道意图识别方法及*** | |
CN109542081B (zh) | 一种基于线下车辆减速度曲线的线上行车危险估计方法 | |
CN110796859A (zh) | 基于交通流的实时交通状态辨识及事故风险预警方法 | |
CN112613786B (zh) | 基于行车安全事件的个性差异化驾驶风险评价方法 | |
CN111563555A (zh) | 一种司机驾驶行为分析方法及*** | |
CN114037015A (zh) | 基于驾驶模式原语的驾驶行为安全性分析方法 | |
CN113581188A (zh) | 一种基于车联网数据的商用车驾驶员驾驶风格识别方法 | |
CN113761715B (zh) | 一种基于高斯混合和隐马尔可夫的个性化车辆跟驰模型的建立方法 | |
CN116279484A (zh) | 一种融合演化博弈与机器学习的多风格驾驶者强制换道预测方法 | |
CN115376316B (zh) | 一种基于隐式马尔可夫链的危险跟驰行为识别方法 | |
Chandrasiri et al. | Driving skill analysis using machine learning the full curve and curve segmented cases | |
CN114613131B (zh) | 一种基于安全裕度的个性化前向碰撞预警方法 | |
CN112990563A (zh) | 一种高速公路追尾事故风险实时预测方法 | |
CN115482662B (zh) | 一种危险工况下驾驶员避撞行为预测方法及*** | |
Angkititrakul et al. | Stochastic mixture modeling of driving behavior during car following | |
CN116975671A (zh) | L3级自动驾驶车辆驾驶员信任度评估方法及*** | |
CN115171389B (zh) | 基于gmm-hmm的高速公路他车超车换道意图识别方法 | |
CN114386737B (zh) | 一种移动互联时代下网约车乘客安全监测方法 | |
Safont et al. | Multichannel signal processing for road surface identification | |
Yu et al. | Drivers' car-following correlative behavior with preceding vehicles in multilane driving | |
Wang et al. | A data-driven estimation of driving style using deep clustering | |
CN112319479B (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 |