CN110562261A - 一种基于马尔可夫模型检测驾驶员风险等级的方法 - Google Patents

一种基于马尔可夫模型检测驾驶员风险等级的方法 Download PDF

Info

Publication number
CN110562261A
CN110562261A CN201910919365.7A CN201910919365A CN110562261A CN 110562261 A CN110562261 A CN 110562261A CN 201910919365 A CN201910919365 A CN 201910919365A CN 110562261 A CN110562261 A CN 110562261A
Authority
CN
China
Prior art keywords
driver
risk level
markov model
risk
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.)
Granted
Application number
CN201910919365.7A
Other languages
English (en)
Other versions
CN110562261B (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.)
Ziyun Intelligent Technology Shanghai Co ltd
Original Assignee
Changan University
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 Changan University filed Critical Changan University
Priority to CN201910919365.7A priority Critical patent/CN110562261B/zh
Publication of CN110562261A publication Critical patent/CN110562261A/zh
Application granted granted Critical
Publication of CN110562261B publication Critical patent/CN110562261B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/08Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to drivers or passengers
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • B60W40/08Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to drivers or passengers
    • B60W40/09Driving style or behaviour

Abstract

本发明公开了一种基于马尔可夫模型检测驾驶员风险等级的方法,本发明通过采集车辆的速度数据以及驾驶员的不同的报警类型,对驾驶行为进行判断并对驾驶员风险等级进行分类,利用训练好的马尔可夫模型对驾驶员风险等级进行识别,再利用识别结果判断驾驶员风险等级。本发明主要服务于运输企业安全管理系统,当识别到驾驶员是高风险等级的驾驶员时,可采取相应的管理培训措施,以提高其安全性。本发明的实用化可以降低因交通事故导致的人员伤亡和财产损失,提高交通系统的整体安全性。

Description

一种基于马尔可夫模型检测驾驶员风险等级的方法
技术领域
本发明属于交通安全领域,具体涉及一种基于马尔可夫模型检测驾驶员风险等级的方法。
背景技术
随着国民经济的快速发展和城市化进程的加快,我国机动车拥有量和道路交通量迅速增长,交通事故问题日趋突出。已有研究表明,驾驶员因素是导致交通事故的主要原因,驾驶员的驾驶风险等级不同,对交通事故的贡献也就不同,低风险等级的驾驶员可能引发较少的甚至避免交通事故,而风险等级越高的驾驶员引发的交通事故可能越严重。因此,针对驾驶员风险等级识别方法进行深入研究显得尤为重要。
目前,已有成果主要是针对驾驶员风险等级评估进行研究,针对驾驶员风险等级识别的研究较少,不能满足交通安全管理的要求。
发明内容
本发明的目的在于克服上述不足,提供一种基于马尔可夫模型检测驾驶员风险等级的方法,能够根据不同的驾驶行为识别驾驶员风险等级,为降低交通事故的发生提供参考。
为了达到上述目的,本发明包括以下步骤:
步骤一,采集若干组驾驶员的速度数据,根据驾驶员的速度数据,确定驾驶员的驾驶行为序列和驾驶员风险等级;
步骤二,建立马尔可夫模型,将同一风险等级的驾驶员驾驶行为作为马尔可夫模型的输入,确定马尔可夫模型中的各项参数;
步骤三,用确定参数后的马尔可夫模型检测待识别的驾驶行为序列;
步骤四,将检测的驾驶行为序列在预设时间段内统计作为该驾驶员预设置信度下最终的风险等级判定结果。
步骤一中,根据速度数据计算车辆的加速度数据,加速度数据的计算方法如下:
vt为t时刻的速度,vt-1为t-1时刻的速度,Δt为时间间隔,根据车辆的GPS数据获得,at为t时段的平均加速度。
步骤一中,确定驾驶行为时,首先根据驾驶员的速度数据对驾驶员的加速信息进行分析,再根据加速信息将驾驶员行为分为若干种,根据车辆内置的报警系统,采集驾驶人的报警次数,通过k-means聚类出驾驶人的风险等级。
步骤二中,马尔可夫模型的参数为P={C,π},C为观察值转移矩阵,π为初始状态矩阵。
步骤三中,用确定参数后的马尔可夫模型检测待识别的驾驶行为等级时,输出的结果为选择概率最大的马尔可夫模型所对应的风险等级作为该驾驶行为序列的识别结果,其中,其中,Pi为t时刻风险等级为j的概率;为风险等级j的初始概率分布,为风险等级为j的驾驶员t时刻驾驶行为为k,t+1时刻驾驶行为为l的转移概率,T指序列长度。
步骤四的具体方法如下:
第一步,根据T长度的驾驶行为序列识别结果判断短时间段内驾驶员的风险等级,判断方法如下:
其中,表示短时间内驾驶员风险等级的识别结果,为风险等级从低至高,P1x2、P3分别表示通过马尔可夫模型将待测驾驶行为序列识别的风险等级从低至高的概率。
第二步,根据短时间段内驾驶员风险等级的判断结果,计算长时间内三种风险等级的数量,计算公式如下:
n1是在长时间段内被识别为低风险等级的序列的数量,n2是在长时间段内被识别为中风险等级的序列数量,n3是在长时间段内被识别为高风险等级的序列的数量,n为总数,为短时间内驾驶员风险等级的识别结果,xj用作计数;
第三步,根据长时间内三种风险等级的数量,计算各风险等级的比例,计算公式如下:
第四步,用信息熵表示随机变量的不确定性,计算公式如下:
令C=1-H,表示驾驶员风险等级的置信度,C值越大,置信度越高;
第五步,在置信度为C的情况下,驾驶员长时间段内的风险等级判断方法如下:
其中,表示驾驶员长时间段风险等级识别结果,I1、I2、I3分别表示低风险等级、中风险等级和高风险等级。
与现有技术相比,本发明通过采集车辆的速度数据以及驾驶员的不同的报警类型,对驾驶行为进行判断并对驾驶员风险等级进行分类,利用训练好的马尔可夫模型对驾驶员风险等级进行识别,再利用识别结果判断驾驶员风险等级。本发明主要服务于运输企业安全管理系统,当识别到驾驶员是高风险等级的驾驶员时,可采取相应的管理培训措施,以提高其安全性。本发明的实用化可以降低因交通事故导致的人员伤亡和财产损失,提高交通系统的整体安全性。
附图说明
图1为本发明的马尔可夫模型图。
具体实施方式
下面结合附图对本发明做进一步说明。
本发明首先通过对各驾驶员的速度数据进行分析,设计加速度指标,并利用特征指标,识别得出驾驶员的驾驶行为序列。而后利用所获取的每位驾驶员的报警类型,采用K-Means均值聚类方法将驾驶员分为低风险等级驾驶员、中风险等级驾驶员以及高风险等级驾驶员。驾驶员风险等级分类是驾驶员风险等级识别的基础,其分类合理与否直接决定了识别算法的成败。利用部分数据对马尔可夫模型进行训练后,再对驾驶行为序列进行识别。
其中,驾驶行为包括快减速、慢减速、正常驾驶、慢加速、快加速共5种。
1、车辆加速度的计算;
车辆加速度是车辆短时间内的平均加速度,具体如式(1)所示。
式中:vt为t时刻的速度,vt-1为t-1时刻的速度,Δt为时间间隔,根据车辆的GPS数据获得,at为t时段的平均加速度。
2、驾驶行为以及驾驶员风险等级的确定;
用于区别不同驾驶行为的加速度阈值尚无明确规定,根据式(1)计算所得加速度确定判别方法如下:
at为t时段的平均加速度,1,2,3,4,5分别表示快减速、慢减速、正常驾驶、慢加速、快加速。具体的确定方法是按从大到小的加速度排序取30%位和70%位加减速度作为阈值。
确定驾驶员风险等级主要将不同驾驶报警类型利用K-Means聚类法,将驾驶员分为低风险等级驾驶员、中风险等级驾驶员和高风险等级驾驶员。
马尔可夫模型的参数包括:
状态数目N,状态集合为Q={q1,q2,...,qN},q1,q2,...,qN对应驾驶员的加速信息类别;观察状态数目M,观察状态集合为V={v1,v2,...vM},v1,v2,...vM对应驾驶行为风险等级;观察序列为o1,o2,...oT,T为序列长度;观察值转移概率矩阵为C,表示风险等级为qj的驾驶员观察值由vk转变为vl的概率,j=1,2,3;初始状态分布概率为
3、基于马尔可夫模型检测驾驶员风险等级;
首先建立马尔可夫模型,具体方法为将同一风险等级的驾驶员驾驶行为作为马尔可夫模型的输入,求解马尔可夫模型的模型参数,P={C,π},C为观察值转移矩阵,π为初始状态矩阵。
然后利用所建立的模型对驾驶行为序列进行识别,所述输出结果为:选择概率最大的马尔可夫模型所对应的风险等级作为该驾驶行为序列的识别结果。
所述将检测结果在一定时间段内统计得到该驾驶员一定置信度下最终的风险等级判定结果,分为以下几步:
①根据T长度的驾驶行为序列识别结果判断短时间段内驾驶员的风险等级,判断方法如式(2):
式中,表示短时间内驾驶员风险等级的识别结果,分别表示低风险等级、中风险等级和高风险等级,P1、P2、P3分别表示通过马尔可夫模型将待测驾驶行为序列识别为低风险等级、中风险等级和高风险等级的概率。
②根据短时间段内驾驶员风险等级的判断结果,计算长时间内三种风险等级的数量,计算公式如式(3):
式中,n1是在长时间段内被识别为低风险等级的序列的数量,n2是在长时间段内被识别为中风险等级的序列数量,n3是在长时间段内被识别为高风险等级的序列的数量,n为总数,为①中短时间内驾驶员风险等级的识别结果,xj只用作计数。
③根据长时间内三种风险等级的数量,计算各风险等级的比例,计算公式如式(4):
式中,k1、k2、k3分别表示长时间段内低风险等级、中风险等级、高风险等级的识别比例。
④用信息熵表示随机变量的不确定性,计算公式如式(5):
令C=1-H,表示驾驶员风险等级的置信度,C值越大,置信度越高。
⑤在置信度为C的情况下,驾驶员长时间段内的风险等级判断方法如式(6):
式中,表示驾驶员长时间段风险等级识别结果,I1、I2、I3分别表示低风险等级、中风险等级和高风险等级。
驾驶员的风险等级的识别结果准确与否,可能与诸多因素的选取有关,因此申请者基于驾驶员(n>30)的驾驶实测数据对驾驶行为的划分、风险等级的分类进行数据分析,研究表明当检测序列长度为8时,检测结果较为稳定,因此训练模型以及模型检测时,本发明选定驾驶行为序列长度为8。在考虑长时间段内驾驶员风险等级的识别中,本发明选定时段长为24小时。最终的识别结果则根据所给置信度C值作为参考。
本发明通过实时采集车辆的GPS数据以及驾驶员不同类型的报警,通过驾驶行为的划分以及风险等级的分类,训练马尔可夫模型并利用模型对驾驶行为序列进行检测,并克服基于短时段的检测结果的不准确性以及不稳定性,通过长时段的检测来考量驾驶员的风险等级,并对最终的识别结果给出相应的置信度。本发明主要服务于车辆安全系统和运输企业安全管理系统,不同风险等级的驾驶员给交通安全带来的隐患程度也不同,运输企业在驾驶员的管理层面可以合理的进行分批、分类、不同程度的教育,从而降低营运车辆的事故频数,为运输企业的管理带来便捷,提高运输企业的运营效益。另外此发明可以在驾驶员驾驶状态出现长时间波动时,对其给予合理的建议,从而提高驾驶员的安全性。本发明的实用化可以降低因交通事故导致的人员伤亡和财产损失,提高交通系统的整体安全性。

Claims (6)

1.一种基于马尔可夫模型检测驾驶员风险等级的方法,其特征在于,包括以下步骤:
步骤一,采集若干组驾驶员的速度数据,根据驾驶员的速度数据,确定驾驶员的驾驶行为序列和驾驶员风险等级;
步骤二,建立马尔可夫模型,将同一风险等级的驾驶员驾驶行为作为马尔可夫模型的输入,确定马尔可夫模型中的各项参数;
步骤三,用确定参数后的马尔可夫模型检测待识别的驾驶行为序列;
步骤四,将检测的驾驶行为序列在预设时间段内统计作为该驾驶员预设置信度下最终的风险等级判定结果。
2.根据权利要求1所述的一种基于马尔可夫模型检测驾驶员风险等级的方法,其特征在于,步骤一中,根据速度数据计算车辆的加速度数据,加速度数据的计算方法如下:
vt为t时刻的速度,vt-1为t-1时刻的速度,Δt为时间间隔,根据车辆的GPS数据获得,at为t时段的平均加速度。
3.根据权利要求1所述的一种基于马尔可夫模型检测驾驶员风险等级的方法,其特征在于,步骤一中,确定驾驶行为时,首先根据驾驶员的速度数据对驾驶员的加速信息进行分析,再根据加速信息将驾驶员行为分为若干种,根据车辆内置的报警系统,采集驾驶人的报警次数,通过k-means聚类出驾驶人的风险等级。
4.根据权利要求3所述的一种基于马尔可夫模型检测驾驶员风险等级的方法,其特征在于,步骤二中,马尔可夫模型的参数为P={C,π},C为观察值转移矩阵,π为初始状态矩阵。
5.根据权利要求1所述的一种基于马尔可夫模型检测驾驶员风险等级的方法,其特征在于,步骤三中,用确定参数后的马尔可夫模型检测待识别的驾驶行为等级时,输出的结果为选择概率最大的马尔可夫模型所对应的风险等级作为该驾驶行为序列的识别结果,其中,Pj为t时刻风险等级为j的概率;为风险等级j的初始概率分布,为风险等级为j的驾驶员t时刻驾驶行为为k,t+1时刻驾驶行为为l的转移概率,T指序列长度。
6.根据权利要求1所述的一种基于马尔可夫模型检测驾驶员风险等级的方法,其特征在于,步骤四的具体方法如下:
第一步,根据T长度的驾驶行为序列识别结果判断短时间段内驾驶员的风险等级,判断方法如下:
其中,表示短时间内驾驶员风险等级的识别结果,为风险等级从低至高,P1、P2、P3分别表示通过马尔可夫模型将待测驾驶行为序列识别的风险等级从低至高的概率。
第二步,根据短时间段内驾驶员风险等级的判断结果,计算长时间内三种风险等级的数量,计算公式如下:
n1是在长时间段内被识别为低风险等级的序列的数量,n2是在长时间段内被识别为中风险等级的序列数量,n3是在长时间段内被识别为高风险等级的序列的数量,n为总数,xj用作计数;
第三步,根据长时间内三种风险等级的数量,计算各风险等级的比例,计算公式如下:
第四步,用信息熵表示随机变量的不确定性,计算公式如下:
令C=1-H,表示驾驶员风险等级的置信度,C值越大,置信度越高;
第五步,在置信度为C的情况下,驾驶员长时间段内的风险等级判断方法如下:
其中,表示驾驶员长时间段风险等级识别结果,I1、I2、I3分别表示低风险等级、中风险等级和高风险等级。
CN201910919365.7A 2019-09-26 2019-09-26 一种基于马尔可夫模型检测驾驶员风险等级的方法 Active CN110562261B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910919365.7A CN110562261B (zh) 2019-09-26 2019-09-26 一种基于马尔可夫模型检测驾驶员风险等级的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910919365.7A CN110562261B (zh) 2019-09-26 2019-09-26 一种基于马尔可夫模型检测驾驶员风险等级的方法

Publications (2)

Publication Number Publication Date
CN110562261A true CN110562261A (zh) 2019-12-13
CN110562261B CN110562261B (zh) 2021-07-23

Family

ID=68782747

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910919365.7A Active CN110562261B (zh) 2019-09-26 2019-09-26 一种基于马尔可夫模型检测驾驶员风险等级的方法

Country Status (1)

Country Link
CN (1) CN110562261B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112158199A (zh) * 2020-09-25 2021-01-01 北京百度网讯科技有限公司 巡航控制方法、装置、设备、车辆及介质
US20210101618A1 (en) * 2019-10-02 2021-04-08 Upstream Security, Ltd. System and method for connected vehicle risk detection
CN112651666A (zh) * 2021-01-16 2021-04-13 北京工业大学 一种基于驾驶模式转移特征的驾驶员风险评估方法
CN113011712A (zh) * 2021-02-25 2021-06-22 长安大学 基于欧氏距离和信息熵的驾驶员驾驶稳定性评估方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011014037A (ja) * 2009-07-03 2011-01-20 Fuji Heavy Ind Ltd リスク予測システム
CN105528593A (zh) * 2016-01-22 2016-04-27 江苏大学 一种前向车辆驾驶人驾驶行为预测系统及预测方法
WO2017217044A1 (ja) * 2016-06-17 2017-12-21 アイシン精機株式会社 視認方向推定装置
CN107958269A (zh) * 2017-11-28 2018-04-24 江苏大学 一种基于隐马尔可夫模型的行车危险度预测方法
CN109177982A (zh) * 2018-10-31 2019-01-11 吉林大学 考虑驾驶风格的车辆行驶危险度评价方法
CN110222596A (zh) * 2019-05-20 2019-09-10 浙江零跑科技有限公司 一种基于视觉的驾驶员行为分析防作弊方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011014037A (ja) * 2009-07-03 2011-01-20 Fuji Heavy Ind Ltd リスク予測システム
CN105528593A (zh) * 2016-01-22 2016-04-27 江苏大学 一种前向车辆驾驶人驾驶行为预测系统及预测方法
WO2017217044A1 (ja) * 2016-06-17 2017-12-21 アイシン精機株式会社 視認方向推定装置
CN107958269A (zh) * 2017-11-28 2018-04-24 江苏大学 一种基于隐马尔可夫模型的行车危险度预测方法
CN109177982A (zh) * 2018-10-31 2019-01-11 吉林大学 考虑驾驶风格的车辆行驶危险度评价方法
CN110222596A (zh) * 2019-05-20 2019-09-10 浙江零跑科技有限公司 一种基于视觉的驾驶员行为分析防作弊方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20210101618A1 (en) * 2019-10-02 2021-04-08 Upstream Security, Ltd. System and method for connected vehicle risk detection
CN112158199A (zh) * 2020-09-25 2021-01-01 北京百度网讯科技有限公司 巡航控制方法、装置、设备、车辆及介质
CN112651666A (zh) * 2021-01-16 2021-04-13 北京工业大学 一种基于驾驶模式转移特征的驾驶员风险评估方法
CN113011712A (zh) * 2021-02-25 2021-06-22 长安大学 基于欧氏距离和信息熵的驾驶员驾驶稳定性评估方法
CN113011712B (zh) * 2021-02-25 2024-01-19 长安大学 基于欧氏距离和信息熵的驾驶员驾驶稳定性评估方法

Also Published As

Publication number Publication date
CN110562261B (zh) 2021-07-23

Similar Documents

Publication Publication Date Title
CN110588658B (zh) 一种基于综合模型检测驾驶员风险等级的方法
CN110562261B (zh) 一种基于马尔可夫模型检测驾驶员风险等级的方法
WO2020244288A1 (zh) 一种基于gps轨迹数据的大货车驾驶行为评估方法及装置
CN112508392B (zh) 一种山区双车道公路隐患路段交通冲突风险动态评估方法
CN110866677B (zh) 一种基于基准化分析的驾驶人相对风险评价方法
CN103996287B (zh) 一种基于决策树模型的车辆强制换道决策方法
CN110544373B (zh) 一种基于北斗车联网的货车预警信息提取与风险识别方法
CN103150900A (zh) 一种基于视频的交通拥堵事件自动检测方法
AU2020103488A4 (en) Method and device for evaluating driving behaviour of truck based on gps trajectory data
CN105809193B (zh) 一种基于kmeans算法的非法运营车辆的识别方法
CN108769104B (zh) 一种基于车载诊断系统数据的路况分析预警方法
CN111401414B (zh) 一种基于自然驾驶数据的危险场景提取及分类方法
CN114783183A (zh) 一种基于交通态势算法的监控方法与系统
US20220383738A1 (en) Method for short-term traffic risk prediction of road sections using roadside observation data
CN110781873A (zh) 一种双模态特征融合的驾驶员疲劳等级识别方法
CN113870564B (zh) 一种封闭路段交通拥堵分类方法、系统、电子设备及存储介质
CN111563555A (zh) 一种司机驾驶行为分析方法及系统
CN110705628B (zh) 一种基于隐马尔可夫模型检测驾驶员风险等级的方法
CN113643552B (zh) 一种智能限速提醒系统的整车测试方法
CN109598931A (zh) 基于交通安全风险的群体划分与差异性分析方法及系统
CN109191828B (zh) 基于集成学习的交通参与者事故风险预测方法
CN110634288B (zh) 基于三元高斯混合模型的多维城市交通异常事件识别方法
CN113581188A (zh) 一种基于车联网数据的商用车驾驶员驾驶风格识别方法
CN110473085B (zh) 一种车辆风险判别方法和装置
CN114299742B (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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230801

Address after: Room 601-19, No. 6, Lane 288, Qianfan Road, Xinqiao Town, Songjiang District, Shanghai, 201600

Patentee after: Ziyun Intelligent Technology (Shanghai) Co.,Ltd.

Address before: 710064 middle section, south two ring road, Shaanxi, Xi'an

Patentee before: CHANG'AN University