CN116680519A - 一种基于gm-hmm模型的车辆换道意图识别方法 - Google Patents
一种基于gm-hmm模型的车辆换道意图识别方法 Download PDFInfo
- Publication number
- CN116680519A CN116680519A CN202310413716.3A CN202310413716A CN116680519A CN 116680519 A CN116680519 A CN 116680519A CN 202310413716 A CN202310413716 A CN 202310413716A CN 116680519 A CN116680519 A CN 116680519A
- Authority
- CN
- China
- Prior art keywords
- vehicle
- lane
- lane change
- hmm
- intention
- 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.)
- Pending
Links
- 230000008859 change Effects 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000009499 grossing Methods 0.000 claims abstract description 24
- 239000000203 mixture Substances 0.000 claims abstract description 20
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 17
- 238000012549 training Methods 0.000 claims abstract description 15
- 230000008569 process Effects 0.000 claims abstract description 11
- 238000012360 testing method Methods 0.000 claims abstract description 8
- 230000006399 behavior Effects 0.000 claims description 35
- 238000009826 distribution Methods 0.000 claims description 9
- 238000011160 research Methods 0.000 claims description 9
- 238000010586 diagram Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 238000012512 characterization method Methods 0.000 claims description 4
- 238000004458 analytical method Methods 0.000 claims description 3
- 238000004140 cleaning Methods 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 8
- 239000011159 matrix material Substances 0.000 description 5
- 238000001914 filtration Methods 0.000 description 4
- 230000003993 interaction Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000012706 support-vector machine Methods 0.000 description 2
- 102100034112 Alkyldihydroxyacetonephosphate synthase, peroxisomal Human genes 0.000 description 1
- 101000799143 Homo sapiens Alkyldihydroxyacetonephosphate synthase, peroxisomal Proteins 0.000 description 1
- 206010039203 Road traffic accident Diseases 0.000 description 1
- 238000000848 angular dependent Auger electron spectroscopy Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000002474 experimental method Methods 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
- 238000010606 normalization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/10—Pre-processing; Data cleansing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/217—Validation; Performance evaluation; Active pattern learning techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/29—Graphical models, e.g. Bayesian networks
- G06F18/295—Markov models or related models, e.g. semi-Markov models; Markov random fields; Networks embedding Markov models
-
- 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
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明涉及一种基于GM‑HMM模型的车辆换道意图识别方法，包括以下步骤：1)基于NGSIM数据集，提取样本数据集；2)对样本数据集进行平滑处理；3)以换道过程中斜率的四分位数为界限，截取2s的目标车辆数据，并基于目标车辆轨迹特征差异显著性选取模型的输入特征，分为训练集和测试集；4)构建高斯混合隐马尔科夫模型GM‑HMM；5)将训练集输入GM‑HMM中，分别对向左换道、向右换道和车道保持三种意图识别进行训练；6)将测试集输入GM‑HMM中，进行换道意图识别，并采用Viterbi算法对换道意图识别的准确率进行验证分析。与现有技术相比，本发明具有准确识别车辆换道意图、提高行驶安全性等优点。
Description
技术领域
本发明涉及驾驶意图识别及机器学习领域，尤其是涉及一种基于GM-HMM模型的车辆换道意图识别方法。
背景技术
据统计，大多数交通事故都是由驾驶人的操作失误引起，而其中相当一部分则是发生于车辆在进行换道(超车/切入)或转向动作过程。报告显示，如果驾驶员能多0.5秒的警告时间则60％的追尾事故可以避免，如果再多1秒的警告时间则90％的追尾事故可以避免。由此来看，能精确高效识别车辆换道意图、提前预知车辆行驶状态对提升驾驶安全性、优化驾驶行为等方面尤为重要。为此，AV或新一代高级驾驶辅助系统(Advanced DrivingAssistance System,ADAS)感知决策模块需要具备识别周围环境车辆驾驶意图的能力，从而为规划安全路径、保持车辆平稳跟随状态提供基础。
近年来，国内外学者针对驾驶行为进行了大量的研究，验证了多种不同模型的实施效果，例如支持向量机模型(SVM)、神经网络、模糊推理、贝叶斯网络、隐马尔可夫模型(HMM)及其改进方法、深度学习和其他算法。综合来看，现有的研究虽然取得了良好的成果，但仍存在一些不足。首先，现有驾驶意图识别方法大多针对主车的驾驶行为，利用车载数据获取主车的详细驾驶操作来推断其驾驶意图，这限制了其在预测周围车辆驾驶意图方面的应用；其次，现有研究多采用轨迹离散数据识别换道意图，较少考虑车辆历史行驶信息中存在时序特征和空间特征，从而影响了换道意图识别的准确性和科学性。因此亟需结合车辆历史行驶信息和周边车辆驾驶行为，形成一套高效、准确的车辆换道意图识别模型，以便于实时检测周边车辆的换道意图。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于GM-HMM模型的车辆换道意图识别方法。
本发明的目的可以通过以下技术方案来实现：
一种基于GM-HMM模型的车辆换道意图识别方法，所述的方法包括以下步骤：
1)基于NGSIM数据集，提取由换道行为和非换道行为共同组成的样本数据集；
2)采用高次多项式平滑处理算法对样本数据集进行平滑处理；
3)以车辆换道过程中斜率的四分位数为界限，截取2s的目标车辆数据用以完整描述反向波动现象，并基于目标车辆轨迹特征差异显著性选取模型的输入特征，其中2/3作为训练集，剩余1/3作为测试集；
4)结合车辆运动的数据多维和连续的特点，通过高斯混合模型与隐马尔科夫模型相结合，形成用于车辆换道意图识别的高斯混合隐马尔科夫模型GM-HMM；
5)将训练集输入GM-HMM模型中，分别对向左换道、向右换道和车道保持三种意图识别进行训练，得到训练完成的GM-HMM模型；
6)将测试集输入训练完成的GM-HMM模型中，进行车辆换道意图识别，并采用Viterbi算法对车辆进行换道意图识别的准确率进行验证分析。
进一步地，所述的1)中的样本数据集包括543个换道行为和870个非换道行为。
进一步地，所述的1)具体为：
101)选取NGSIM数据集中US 101号公路上的车辆作为研究对象，对车辆行驶数据进行数据清洗和去重，并统一时间粒度为0.1s，删除缺失数据；
102)根据101)中处理后的数据，将目标车辆和周边车辆布置在同一坐标系中，其中周边车辆包括本车道前车、本车道后车、目标车道前车、目标车道后车以及目标车道并行车辆，并利用车辆行驶横纵向坐标及时间戳，确定车辆周围车辆，得到目标车辆和周边车辆的行驶轨迹数据；
103)将所有目标车辆及其周边车辆的行驶轨迹数据构成包括543个换道行为和870个非换道行为的样本数据集。
进一步地，所述的行驶轨迹数据包括车辆的长度、宽度、位置、速度和加速度。
进一步地，所述的2)具体为：综合对比不同数据平滑算法，选取高次多项式平滑处理算法对样本数据集进行平滑处理，处理后的样本数据集的时间间隔为0.1s。
进一步地，所述的3)具体为：
301)换道意图识别点确定：通过目标车辆的斜率确定换道意图识别点，向左换道的车辆以车辆斜率的上四分位数作为意图识别的观测点，向右换道的车辆则以车辆斜率的下四分位数作为意图识别的观测点，计算公式如下：
其中，kit为第i辆换道车辆在时刻t的斜率；xit为第i辆车在t时刻的横向位置；yit为第i辆车在t时刻的纵向位置，且向左换道斜率为正值，向右滑换道斜率为负值；
302)基于目标车辆的横向位置的标准差确定时间窗长度，截取2s的目标车辆数据用以完整描述反向波动现象；
303)车辆换道意图识别特征优选；将筛选出的车队轨迹数据按照目标车辆的运动轨迹划分为向左换道车队、向右换道车队以及车道保持车队运动轨迹数据，并基于箱型图对比车辆不同运动特征分布的差异性，得到优选的换道意图识别特征；
304)将换道意图识别特征归一化处理，作为模型的输入特征。
进一步地，所述的4)具体为：
401)构建隐马尔科夫模型HMM；
402)引入新的连续化表征方法高斯混合模型GMM；
403)将GMM和HMM相结合，形成用于车辆换道意图识别的高斯混合隐马尔科夫模型GM-HMM。
进一步地，所述的采用Viterbi算法对车辆进行换道意图识别的准确率进行验证分析具体为：计算车辆换道意图识别的预测值和真实值一致的车队编号数量与总车队数量的比值，比值越高，则换道意图识别的准确率越高；比值越低，则换道意图识别的准确率越低。
一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，所述处理器执行所述程序时实现如上所述的方法。
一种计算机可读存储介质，其上存储有计算机程序，所述程序被处理器执行时实现如上所述的方法。
与现有技术相比，本发明具有以下有益效果：
一、本发明首次提出基于轻微反向波动现象构建GM-HMM模型，并在原有仅考虑与目标车辆物理状态相关的特征，如位置、速度、加速度等基础上，考虑目标车辆与其邻域车辆间的交互作用，相较于利用车载数据获取主车的详细驾驶操作进行换道意图预测，更具通用性，拓展了其在预测周围车辆驾驶意图方面的应用，且换道意图识别的准确率更高。
二、本发明充分考虑车辆历史行驶信息中存在的时序特征，相较于传统采用轨迹离散数据识别换道意图，且选用多组不同时间窗的目标车辆运动轨迹数据求解GM-HMM模型，发现以目标车辆横向位置波动的标准差为指标，且以换道过程中斜率的四分位数为换道行为观测的起始点并以2s为观测长度，可以较好的描述车辆反向波动行为，更贴近实际情形。
三、本发明在应用到实际中时，可以降低自动驾驶车辆规划决策难度和失效概率、提高车辆行驶的安全性和平顺性，相比于现有算法在识别准确率上有明显提升，模型准确率99％左右，且模型能够准确识别提前1～3s的车辆换道行为。
附图说明
图1为本发明的换道行为示意图；
图2为本发明的车辆换道意图识别方法流程示意图；
图3为本发明的最小二乘法平滑的基本方法示意图；
图4为本发明的车辆轨迹平滑处理对比图；
图5为本发明的车辆轨迹速度平滑处理对比图；
图6为本发明的目标车辆以及周围车辆轨迹示意图；
图7为本发明的换道观测点示意图；
图8为本发明的换道时间窗长度示意图；
图9为本发明的目标车辆的反向波动特征示意图；
图10为本发明的向左目标车辆换道前横向位置波动标准差对比图；
图11为本发明的车辆轨迹数据对比箱型图，其中a表示本车速度分布，b表示前车速度分布，c表示前车与本车横向相对位置，d表示目标车道前车速度分布，e表示目标车道后车速度分布，f表示目标车道重叠车辆速度，g表示目标车道重叠车辆纵向位置；
图12为本发明的目标车辆换道点示意图；
图13为本发明的换道意图识别模型检测换道行为提前时间分布图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
如图1所示，当目标车辆需要进行换道时，需要本车道和目标车道的周边车辆行驶状态，本发明可利用微观层面目标车辆及其周边车辆轨迹数据，分析目标车辆在跟驰与换道行驶过程中自身行驶轨迹与周边车辆交互行为特征，结合目标车辆轨迹的反向波动现象，通过训练好的GM-HMM模型实现车辆换道意图识别。
本发明所依托的车辆轨迹数据集为NGSIM数据集，主要涉及车辆的轨迹、长度、类型、宽度、速度、加速度、车头间距和车头时距等数据，具体字段名和其含义如表1所示，且在NGSIM数据集中，可以根据数据处理得到真实驾驶环境中的交互信息。
表1NGSIM数据集字段
如图2所示，一种基于GM-HMM模型的车辆换道意图识别方法，所述的方法包括以下步骤：
1)基于NGSIM数据集，提取由543换道行为和870非换道行为共同组成的样本数据集。
所述的1)具体为：
101)选取NGSIM数据集中US 101号公路上的车辆作为研究对象，对车辆行驶数据进行数据清洗和去重，并统一时间粒度为0.1s，删除缺失数据；
102)根据101)中处理后的数据，根据NGSIM数据集中的车辆位置、速度等信息将目标车辆和周边车辆布置在同一坐标系中，其中周边车辆包括本车道前车、本车道后车、目标车道前车、目标车道后车以及目标车道并行车辆，并利用车辆行驶横纵向坐标及时间戳，确定车辆周围车辆，得到目标车辆和周边车辆的行驶轨迹数据，包括车辆的长度、宽度、位置、速度和加速度等。
103)将所有目标车辆及其周边车辆的行驶轨迹数据构成包括543个换道行为和870个非换道行为的样本数据集。
2)采用高次多项式平滑处理算法对样本数据集进行平滑处理。
由于NGSIM数据集存在一定的噪声，特别是车辆的速度抖动较大，因此需要对样本数据集进行平滑处理，本发明对比常见的卷积拟合(Savitzky-Golay滤波拟合法)和多项式拟合法，最终根据二者平滑结果选用高次多项式平滑处理算法对样本数据集平滑处理，具体为：
201)如图3所示，基于Savitzky-Golay滤波拟合法对样本数据集进行平滑处理；首先根据NDVI时间序列曲线的平均趋势，确定合适的滤波参数，用多项式实现滑动窗内的最小二乘拟合，其基本原理如下：
假设一列数据为x[n]，将一组以n＝0为中心的2M+1个数据，使用如下的多项式进行拟合，
其中，p(n)为拟合值，N为给定需要拟合的点数，n为给定的整数值，k为多项式的阶数，ak为第k项的系数；
最小二乘拟合的残差为：
其中，εN为拟合数据点和实际数据点的差值，x[n]为待拟合的数据值，M为单边数据点数；
在图3中，实线表示的是N＝2，M＝2时得到的拟合多项式，滤波的结果为：
y[0]＝p(0)＝a0
因此只需要获得拟合多项式的常数项；计算参数a0相当于对原始数据进行一次FIR滤波(非递归型滤波)，即可以利用卷积运算实现上述计算：
式中h[m]和x[n-m]为卷积的两个变量；由上式可知即对输入数据进行了加权平均，x即为加权系数，对于a0的获取，可以由微积分的基本原理推到得到，若要ε最小，则ε对各个参数的偏导数均应为0，即如下式所示：
化简后为：
此时若引入一个(2M+1)行(N+1)列的辅助矩阵A＝{an,j}。
an,j＝ni,-M≤n≤M,0≤i≤N
再设辅助矩阵B，使得：B＝ATA，经计算可知：
定义：
则可以得到：
Ba＝ATAa＝ATx
a＝(ATA)-1ATx＝Hx
H第一行行向量即为所求卷积系数，由上述推导可知，H与x无关，由N和M来决定；
202)基于多项式拟合法对样本数据集进行平滑处理；首先假设给定数据由高次多项式函数生成，选择最有可能产生这些数据的高次多项式函数，即在高次多项式函数中选择一个对已知数据以及未知数据都有很好预测能力的函数；拟合前的样本数据集中的原始数据和拟合后的轨迹及速度数据如图4和图5所示。
203)综合对比不同数据平滑方法后，本发明选取高次多项式平滑处理算法对样本数据集进行平滑处理，处理后的数据时间间隔为0.1s，平滑后的车辆轨迹以及数据格式分别如图6和表2所示。
表2平滑处理后样本数据集的数据格式
3)以车辆换道过程中斜率的四分位数为界限，截取2s的目标车辆数据用以完整描述反向波动现象，并基于目标车辆轨迹特征差异显著性选取模型的输入特征，其中2/3作为训练集，剩余1/3作为测试集，具体为：
301)换道意图识别点确定：由于HMM在语音识别方面的效果较优，即在拟合具有声波特征的数据方面具有较大的优势，因此尝试通过提取车辆在换道前具有波动特征的数据作为模型的输入。基于该数据集提取换道车队的运动状态，首次发现车辆换道过程中车辆会向换道方向的反方向波动的现象，如图7和图8所示；因此本发明拟通过目标车辆的斜率确定换道意图识别点；计算每辆车的换道过程中斜率的公式如下：
式中kit为第i辆目标车辆在时刻t的斜率；xit为第i辆车在t时刻的横向位置；yit为第i辆车在t时刻的纵向位置，且向左换道斜率为正值，向右滑换道斜率为负值；
选取三个时间单位间隔(0.3s)的点计算车辆的斜率，可以减少短时间内车辆微调造成的换到斜率较大的波动，最终向左换道的车辆以车辆斜率的上四分位数作为意图识别的观测点，向右换道的车辆则以车辆斜率的下四分位数作为意图识别的观测点；
302)基于目标车辆的横向位置的标准差确定时间窗长度，截取出时间窗长度为1s、1.2s、1.5s、1.7s、2.0s、2.2s、2.5s、2.7s和3.0s的目标车辆数据用以完整描述反向波动现象，分别计算不同时间窗长度车辆横向位置波动标准差(fluctuation_std)，如图9所示，得到取时间窗长度为2s可以使数据较为完整的包含目标车辆换道前的波动特征，如图10所示；
303)车辆换道意图识别特征优选；将筛选出的车队轨迹数据按照目标车辆的运动轨迹划分为向左换道车队、向右换道车队以及车道保持车队运动轨迹数据，并基于箱型图对比车辆不同运动特征分布的差异性，得到优选的换道意图识别特征，如图11和表3所示；
表3车辆换道意图识别优选特征
304)为提高模型的预测精度和模型的收敛速度，将换道意图识别特征归一化处理，作为模型的输入特征；由于轨迹数据的特殊性使用普通的归一化方法会缩小车队换道过程中车辆之间的交互特征，因此采取将周围车辆的运动特征与同一时刻目标车辆的对应特征进行求差之后的数据。
4)结合车辆运动的数据多维和连续的特点，通过高斯混合模型与隐马尔科夫模型相结合，形成用于车辆换道意图识别的高斯混合隐马尔科夫模型GM-HMM，具体为：
401)构建隐马尔科夫模型HMM；
HMM模型的参数可以由λ描述，λ＝{A,B,π,M,N}，M代表观测状态数，N代表隐状态数，其他参数含义见表4，参数解释如图11(a)～图11(g)所示；
表4HMM参数含义
HMM主要包含3大问题，分别为评估问题、解码问题和学习问题，在本发明中换道行为被分为3类，即向左换道、向右换道与车道保持；
首先模型参数和观测序列已知，求解该观测序列在给定模型下的条件概率称为评估问题，结合换道意图识别，即在收集时间窗长度为T的车辆数据后计算当前分属向左换道、向右换道、车道保持3种行为的概率。
其中，αt9i)为前向概率，定义为给定隐马尔科夫模型参数λ，到时刻t的观测序列为o1,o2,…,ot，且隐状态为it＝qi的概率：
αt(i)＝P(o1,o2,o3,…,ot,it＝qi|λ)
βt(i)为后向概率，给定隐马尔科夫模型参数λ时，时刻t隐状态为it＝qi下定义部分观测序列为ot+1,ot+2,…,oT的概率为后向概率：
βt(i)＝P(ot+1,ot+2,ot+3,…,oT∣it＝qi,λ)
上述提到需要计算当前状态分属3类行为的概率，意味着本发明需要3组不同的HMM参数，左换道λ1、车道保持λ2以及右换道λ3；参数学习问题在换道行为识别中表现为在给定3类模型参数初值后通过对应行为数据训练集对每组参数迭代更新，即每组确定一个描述已知观测序列的最优参数组合使得训练数据在该模型下输出的概率值最大；定义两个重要参数γi(t),ξt(i,j)，其中，γi(t)代表已知模型参数和观测序列，求解t时刻隐状态为qi的概率：
ξt(i,j)代表已知模型参数和观测序列，在t时刻隐状态为qi，t+1时刻隐状态为qj的概率：
402)由于车辆运动是一个时间连续行为，其运动观测状态不是简单地有限划分，若采用离散化的表征方式将大大减少原始数据的特性，降低最终精度，为最大化保留原始数据特征，提高模型识别精度，引入新的连续化表征方法高斯混合模型GMM；
403)将GMM和HMM相结合，形成用于车辆换道意图识别的高斯混合隐马尔科夫模型GM-HMM；
该模型参数仍可记为但其中的输出观测B不再是一个混淆矩阵,而是一组观察概率密度函数，每一个隐状态量对应的观察量由多个多维高斯函数生成，其对应的观察值概率密度函数如下：
其中，o为观测量，M为隐状j包含的混合高斯元个数，cjm为隐状态j下第m个混合高斯元的权重，N代表多维高斯概率密度函数，μjm为隐状态j下第m个混合高斯元的均值数组，Ujm为隐状态为j下第m个混合高斯元的协方差矩阵；这样高斯混合隐马尔科夫模型的参数变为λ(π,A,Cjm,μjm,Ujm，其中N表达如下：
其中，G为观测变量维数,初始概率分布π和状态转移概率矩阵A与离散隐马尔科夫模型的训练公式相同，其它高斯混合隐马尔科夫模型参数的重估公式为：
其中，γt(j,m)表示t时刻处于隐状态j，观测状态o来自于第m个高斯分量的概率：
虽然现有计算机性能相较以往得到了极大提升，但由于行为识别的时效性要求较高，因此在考虑精度的同时也要兼顾计算量；从GM-HMM的模型构成及参数计算方法可以发现，模型的计算量主要与隐状态数量、观测数据维数、观测数据时间序列长度和高斯混合数相关。这些参数越小则计算量越小，但是考虑到模型的识别效果，这些参数并不是越小越好。
本发明基于NGSIM数据集使用python中的hmmlearn包对HMM模型的参数进行拟合。模型实现的代码以及各参数的含义分别如表5和表6所示：
表5模型实现的代码
表6hmmlearn参数含义
其中隐藏状态的个数是不可自适应的参数且没有实际的物理意义，因此对该参数进行遍历寻优得到最优的准确率的隐藏状态取值(目标车辆n_components＝2、非目标车辆n_components＝1)。
5)将训练集输入GM-HMM模型中，并将车辆三种换道状态作为GM-HMM模型的输出，分别对向左换道、向右换道和车道保持三种意图识别进行训练，得到训练完成的GM-HMM模型。
6)将测试集输入训练完成的GM-HMM模型中，进行车辆换道意图识别，并采用Viterbi算法对车辆进行换道意图识别的准确率进行验证分析。
如图12所示，将向左换道、车道保持、向右换道行为记为-1、0、1，将剩余的1/3车队数据作为验证集，分别带入三种换道意图识别模型中计算该观测序列属于该模型的概率，最后将概率得分最高的模型所代表的车辆行为作为最终的预测值，如表7所示；
表7模型预测结果示例
/>
如图13所示，计算车辆换道意图识别的预测值和真实值一致的车队编号数量与总车队数量的比值，作为模型预测准确率的评价指标，比值越高，则换道意图识别的准确率越高；比值越低，则换道意图识别的准确率越低；此外，还可以计算诸如查准率、召回率、F1-score等指标；最后，为验证GM-HMM模型对车辆换道意图识别的提前预测能力，对模型预测的提前时间进行计算，如表8所示。
表8不同预测时间下模型识别性能对比
相较于已有车辆换道意图识别研究，本发明提出基于轻微反向波动现象构建的GM-HMM模型在提前2-3s预测向左换道的准确率达到98％以上，相比于其他研究在延长了预测时间的前提下，将模型预测准确率提高了2％左右；与其他研究最优的车道保持车辆预测准确率相比，模型将预测准确率提高了5％左右。
本发明提供一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，所述处理器执行所述程序时实现如上所述的方法。
本发明还提供一种计算机可读存储介质，其上存储有计算机程序，所述程序被处理器执行时实现如上所述的方法。
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。
Claims (10)
1.一种基于GM-HMM模型的车辆换道意图识别方法，其特征在于，所述的方法包括以下步骤：
1)基于NGSIM数据集，提取由换道行为和非换道行为共同组成的样本数据集；
2)采用高次多项式平滑处理算法对样本数据集进行平滑处理；
3)以车辆换道过程中斜率的四分位数为界限，截取2s的目标车辆数据用以完整描述反向波动现象，并基于目标车辆轨迹特征差异显著性选取模型的输入特征，其中2/3作为训练集，剩余1/3作为测试集；
4)结合车辆运动的数据多维和连续的特点，通过高斯混合模型与隐马尔科夫模型相结合，形成用于车辆换道意图识别的高斯混合隐马尔科夫模型GM-HMM；
5)将训练集输入GM-HMM模型中，分别对向左换道、向右换道和车道保持三种意图识别进行训练，得到训练完成的GM-HMM模型；
6)将测试集输入训练完成的GM-HMM模型中，进行车辆换道意图识别，并采用Viterbi算法对车辆进行换道意图识别的准确率进行验证分析。
2.根据权利要求1所述的一种基于GM-HMM模型的车辆换道意图识别方法，其特征在于，所述的1)中的样本数据集包括543个换道行为和870个非换道行为。
3.根据权利要求2所述的一种基于GM-HMM模型的车辆换道意图识别方法，其特征在于，所述的1)具体为：
101)选取NGSIM数据集中US 101号公路上的车辆作为研究对象，对车辆行驶数据进行数据清洗和去重，并统一时间粒度为0.1s，删除缺失数据；
102)根据101)中处理后的数据，将目标车辆和周边车辆布置在同一坐标系中，其中周边车辆包括本车道前车、本车道后车、目标车道前车、目标车道后车以及目标车道并行车辆，并利用车辆行驶横纵向坐标及时间戳，确定车辆周围车辆，得到目标车辆和周边车辆的行驶轨迹数据；
103)将所有目标车辆及其周边车辆的行驶轨迹数据构成包括543个换道行为和870个非换道行为的样本数据集。
4.根据权利要求3所述的一种基于GM-HMM模型的车辆换道意图识别方法，其特征在于，所述的行驶轨迹数据包括车辆的长度、宽度、位置、速度和加速度。
5.根据权利要求1所述的一种基于GM-HMM模型的车辆换道意图识别方法，其特征在于，所述的2)具体为：综合对比不同数据平滑算法，选取高次多项式平滑处理算法对样本数据集进行平滑处理，处理后的样本数据集的时间间隔为0.1s。
6.根据权利要求1所述的一种基于GM-HMM模型的车辆换道意图识别方法，其特征在于，所述的3)具体为：
301)换道意图识别点确定：通过目标车辆的斜率确定换道意图识别点，向左换道的车辆以车辆斜率的上四分位数作为意图识别的观测点，向右换道的车辆则以车辆斜率的下四分位数作为意图识别的观测点，计算公式如下：
其中，kit为第i辆换道车辆在时刻t的斜率；xit为第i辆车在t时刻的横向位置；yit为第i辆车在t时刻的纵向位置，且向左换道斜率为正值，向右滑换道斜率为负值；
302)基于目标车辆的横向位置的标准差确定时间窗长度，截取2s的目标车辆数据用以完整描述反向波动现象；
303)车辆换道意图识别特征优选；将筛选出的车队轨迹数据按照目标车辆的运动轨迹划分为向左换道车队、向右换道车队以及车道保持车队运动轨迹数据，并基于箱型图对比车辆不同运动特征分布的差异性，得到优选的换道意图识别特征；
304)将换道意图识别特征归一化处理，作为模型的输入特征。
7.根据权利要求1所述的一种基于GM-HMM模型的车辆换道意图识别方法，其特征在于，所述的4)具体为：
401)构建隐马尔科夫模型HMM；
402)引入新的连续化表征方法高斯混合模型GMM；
403)将GMM和HMM相结合，形成用于车辆换道意图识别的高斯混合隐马尔科夫模型GM-HMM。
8.根据权利要求1所述的一种基于GM-HMM模型的车辆换道意图识别方法，其特征在于，所述的采用Viterbi算法对车辆进行换道意图识别的准确率进行验证分析具体为：计算车辆换道意图识别的预测值和真实值一致的车队编号数量与总车队数量的比值，比值越高，则换道意图识别的准确率越高；比值越低，则换道意图识别的准确率越低。
9.一种电子设备，包括存储器和处理器，所述存储器上存储有计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1～8中任一项所述的方法。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述程序被处理器执行时实现如权利要求1～8中任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310413716.3A CN116680519A (zh) | 2023-04-18 | 2023-04-18 | 一种基于gm-hmm模型的车辆换道意图识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310413716.3A CN116680519A (zh) | 2023-04-18 | 2023-04-18 | 一种基于gm-hmm模型的车辆换道意图识别方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116680519A true CN116680519A (zh) | 2023-09-01 |
Family
ID=87784328
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310413716.3A Pending CN116680519A (zh) | 2023-04-18 | 2023-04-18 | 一种基于gm-hmm模型的车辆换道意图识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116680519A (zh) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107273805A (zh) * | 2017-05-18 | 2017-10-20 | 江苏大学 | 一种基于视觉特性的gm‑hmm预测驾驶行为方法 |
| CN109460023A (zh) * | 2018-11-09 | 2019-03-12 | 上海理工大学 | 基于隐马尔科夫模型的驾驶人换道意图识别方法 |
| CN109471436A (zh) * | 2018-11-09 | 2019-03-15 | 上海理工大学 | 基于混合高斯-隐马尔科夫模型的换道模型参数优化方法 |
| CN113859266A (zh) * | 2021-10-12 | 2021-12-31 | 北京理工大学 | 一种非结构化道路目标车辆换道轨迹预测方法及系统 |
-
2023
- 2023-04-18 CN CN202310413716.3A patent/CN116680519A/zh active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107273805A (zh) * | 2017-05-18 | 2017-10-20 | 江苏大学 | 一种基于视觉特性的gm‑hmm预测驾驶行为方法 |
| CN109460023A (zh) * | 2018-11-09 | 2019-03-12 | 上海理工大学 | 基于隐马尔科夫模型的驾驶人换道意图识别方法 |
| CN109471436A (zh) * | 2018-11-09 | 2019-03-15 | 上海理工大学 | 基于混合高斯-隐马尔科夫模型的换道模型参数优化方法 |
| CN113859266A (zh) * | 2021-10-12 | 2021-12-31 | 北京理工大学 | 一种非结构化道路目标车辆换道轨迹预测方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 白玉等: "基于NGSIM数据的车辆换道前跟驰模型研究", 《交通与运输》, pages 25 - 29 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111079590B (zh) | 一种无人驾驶车辆的周边车辆行为预判方法 | |
| CN112347567B (zh) | 一种车辆意图和轨迹预测的方法 | |
| CN111104969B (zh) | 一种无人驾驶车辆与周边车辆的碰撞可能性预判方法 | |
| CN106971194B (zh) | 一种基于改进hmm和svm双层算法的驾驶意图识别方法 | |
| CN113486822B (zh) | 一种基于驾驶意图的周围车辆轨迹预测方法及系统 | |
| CN112562328B (zh) | 一种车辆行为预测方法及装置 | |
| CN113859266B (zh) | 一种非结构化道路目标车辆换道轨迹预测方法及系统 | |
| Wirthmüller et al. | Predicting the time until a vehicle changes the lane using LSTM-based recurrent neural networks | |
| CN112793576B (zh) | 一种基于规则与机器学习融合的换道决策方法及系统 | |
| CN110928918B (zh) | 时间序列数据组成模式的提取方法、装置及终端设备 | |
| CN112116100B (zh) | 一种考虑驾驶员类型的博弈论决策方法 | |
| CN111695737A (zh) | 一种基于lstm神经网络的群目标行进趋势预测方法 | |
| CN111159642B (zh) | 一种基于粒子滤波的在线轨迹预测方法 | |
| CN113033899B (zh) | 无人驾驶邻近车辆轨迹预测方法 | |
| Martinsson et al. | Clustering vehicle maneuver trajectories using mixtures of hidden markov models | |
| CN114926802A (zh) | 基于双层lstm网络的车辆轨迹预测方法 | |
| CN111907523A (zh) | 一种基于模糊推理的车辆跟驰寻优控制方法 | |
| CN112632706B (zh) | 一种基于gmm和cart回归树的路感模拟方法 | |
| CN116680519A (zh) | 一种基于gm-hmm模型的车辆换道意图识别方法 | |
| Zhang et al. | A lane-changing prediction method based on temporal convolution network | |
| CN113761715B (zh) | 一种基于高斯混合和隐马尔可夫的个性化车辆跟驰模型的建立方法 | |
| Safont et al. | Comparison of dimensionality reduction methods for road surface identification system | |
| CN113722982A (zh) | 一种基于多模型融合的汽车传感器攻击检测与防御方法 | |
| Ma et al. | Lane change analysis and prediction using mean impact value method and logistic regression model | |
| CN115731261B (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 |