[发明专利]车联网复杂网络中基于深度学习的车辆行驶影响因素的位置预测模型构造方法有效
申请号: | 201710702922.0 | 申请日: | 2017-08-16 |
公开(公告)号: | CN107609633B | 公开(公告)日: | 2020-08-14 |
发明(设计)人: | 程久军 | 申请(专利权)人: | 同济大学 |
主分类号: | G06N3/04 | 分类号: | G06N3/04;G06N3/08;G06K9/62;G08G1/00 |
代理公司: | 上海科律专利代理事务所(特殊普通合伙) 31290 | 代理人: | 叶凤 |
地址: | 200092 *** | 国省代码: | 上海;31 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 联网 复杂 网络 基于 深度 学习 车辆 行驶 影响 因素 位置 预测 模型 构造 方法 | ||
1.一种基于深度学习的车辆位置预测模型构造方法,其特征在于,具体方法包括如下步骤:
步骤1,定义特征集,包括车辆特征抽象与定义、道路特征抽象与定义、车辆行驶环境特征抽象与定义
步骤11、车辆特征抽象与定义
车辆特征包括车身长度L,车身宽度W,车辆速度,加速度,当前行驶方向,路口转向动作,当前位置,出发地,目的地,其中t时刻车辆速度、加速度分别标记为v(t)和a(t),其余特征分别定义如下:
定义1车辆行驶方向vDir(t),表示t时刻车辆的运动方向,用与正北方向之间顺时针夹角表示行驶方向,满足:
0≤vDir(t)<360° (1)
定义2车辆路口转向vMov(t),表示t时刻车辆在经过十字路口时的驾驶行为,以向量的形式表征,如公式(2),分为车辆直行,左转弯以及右转弯,
定义3车辆当前位置P(t),表示t时刻车辆在CA State Plane III in NAD83坐标系中的二维坐标向量,向量各项单位为英尺,位置信息定义如下:
P(t)=(x(t),y(t)) (3)
定义4车辆出发地srcP,表示车辆此次出行起点的二维坐标,定义如下:
srcP=(x,y) (4)
坐标系与单位同定义3;
定义5车辆目的地dstP,表示车辆此次出行终点的二维坐标,定义如下:
dstP=(x,y) (5)
坐标系与单位同定义3;
综上,t时刻车辆的特征集featurev(t)定义如下:
featurev(t)={L,W,v(t),a(t),vDir(t),vMov(t),P(t),srcP,dstP} (6)
步骤12,道路特征抽象与定义
将十字路口和丁字路口抽象为四边形,以其四个顶点坐标从西北方向的角开始按照顺时针方向依次标识,标记如下:
I=(x1,y1,x2,y2,x3,y3,x4,y4) (7)
定义6十字路口集合ISet表示研究区域内,所有十字路口组成的集合;
ISet=(I1,I2,...,Im,...) (8)
定义7道路段是指道路在两个相邻十字路口之间的一个路段,由路段两侧十字路口来标识该道路段,定义道路i的道路段集合如下:
RSegSet(i)=(I1I2,I2I3,...,IkIm,...) (9)
其中Ik表示编号为k的十字路口,1≤k≤n,n为最大十字路口编号;
每条道路段包含若干车道,定义十字路口Ik,Im之间的道路段车道集合如下:
IkIm=(lid1,lid2,...,lidn) (10)
其中lidk表示车道编号,1≤k≤n;
定义8车道方向lDir(x,y)表示车道允许的行进方向,其中x,y表示车道内的位置坐标;角度定义同vDir;同时车辆在当前车道直道的行驶方向vDir与车道方向相同,即若当前车辆处在i车道内,则有
vDir(t)=lDiri(x,y),其中(x,y)=P(t) (11)
定义9车道左侧可变道数量LAL为在某车道当前位置下,车辆可向左变动的车道数目;由当前车道i的当前位置(x,y)向左侧同向车道搜索,直至遇到实线表示车辆不可向左变道或者搜索到最左侧同向车道为止;搜索到的车道数目即为车道左侧可变道数量LALi(x,y);
定义10车道右侧可变道数量RAL为在某车道当前位置下,车辆可向右变动的车道数目;由当前车道i的当前位置(x,y)向右侧同向车道搜索,直至遇到实线表示车辆不可向右变道或者搜索到最右侧同向车道为止;搜索到的车道数目即为车道右侧可变道数量RALi(x,y);
定义11直行区车道允许驾驶行为包括直行,向左变道和向右变道;定义向量sld表示车道允许的基本驾驶行为如下:
直行区车道允许的驾驶行为SLD有
定义12路口准备区车道允许的驾驶行为包括直行,左转弯,右转弯和调头;定义向量pld表示车道允许的基本驾驶行为如下:
路口准备区车道允许的驾驶行为PLD表示为:
综上,位置(x,y)处的道路特征集featurer(x,y)定义如下:
featurer(x,y)={lDir(x,y),LAL(x,y),RAL(x,y)} (16)
∪{SLD(x,y),PLD(x,y)}
步骤13,车辆行驶环境特征抽象与定义
定义13路口距离,车辆i与前方十字路口Im距离表示以车辆i前边沿与当前车道前方十字路口停止线之间的距离;
定义14前方车辆与后方车辆;在同一车道上,车辆i与车辆j前方路口为Im,若满足则在t时刻车辆i为车辆j的前方车辆,车辆j为车辆i的后方车辆;若则认为车辆i与车辆j互为对方的前方车辆;
在同一道路段IkIm的不同车道上,假设t时刻车辆i的位置Pi(t)=(xi,t,yi,t),车辆j的位置Pj(t)=(xj,t,yj,t);路段方向向量若满足则在t时刻车辆i为车辆j的前方车辆,车辆j为车辆i的后方车辆;如果则认为车辆i与车辆j互为对方的前方车辆;
定义15周围车辆集NCarSeti(t)定义为车辆i在t时刻当前车道,左侧车道,右侧车道前方与后方车辆所组成的集合;
定义16周围车辆行驶状态以向量nstate来标识如下:
nstate(t)=(type,available,v(t),a(t),L,W,d)∪P(t) (17)
其中,type标识周围车辆与当前车辆的相对位置,available表示是否存在该车辆,存在为1,不存在为0;nstate中除type项外,其他项均标记为0;如车辆若不存在当前车道前方车辆,则该车辆当前车道前方车辆对应的nstate(t)=(1,0,0,0,0,0,0,0,0);v(t),a(t),P(t)分别表示车辆的速度、加速度和位置;L,W分别表示车辆长度和宽度;d表示邻居车辆节点与当前车辆节点之间的距离;
定义17周围车辆运动状态集定义为NCarStatei(t)={nstatej(t)|j∈NCarSeti(t)};
定义18t时刻,交通灯TLi信号允许动作使用向量sigi(t)表示;
定义19车辆t时刻的路口允许转向动作表示为IAM(t);该特征受制于车道允许的驾驶行为PLD和交通灯信号允许动作sigi(t);可表示为矩阵的Hadamard乘积,如公式(19):
IAM(t)=PLD(P(t))*sigi(t) (19)
综上,t时刻车辆行驶环境特征集featuree(t)定义如公式(20):
featuree(t)={VID(t),IAM(t)}∪NCarState(t) (20)
综上,将t时刻影响车辆驾驶行为的特征集feature(t)定义为
feature(t)=featurev(t)∪featurer(P(t))∪featuree(t) (21)
步骤2,车辆位置预测模型
步骤21,特征提取与数据预处理
直接获取的车辆特征包括车辆长度L,车辆宽度W,车辆速度v,车辆加速度a,车辆行驶方向vDir,车辆路口转向vMov,车辆当前位置P,车辆出发地点srcP,车辆目的地点dstP;
直接获取的道路特征包括车道方向1Dir,车道左侧可变道数量LAL,车道右侧可变道数量RAL,车道允许的驾驶行为SLD,路口准备区车道允许的驾驶行为PLD;
直接获取的车辆行驶环境特征包括车辆与当前车道前方车辆的距离df,车辆与当前车道后方车辆的距离da;
直接可获取十字路口集ISet,道路段集合RSegSet,交通灯信号允许的驾驶动作sig(t),每辆车所在的车道lane以及每辆车所在的道路段RSeg;根据前面步骤1的定义,需要提取的特征包括车辆与前方十字路口距离VID,周围车辆运动状态集NCarState,车辆路口允许转向动作IAM,训练样本的标签;
·步骤211,车辆与前方路口距离特征提取
在地图上根据车辆所在车道lane以及车辆所在车道方向IDir获取该车道在前方十字路口处停止线AB的A点坐标(xA,yA)和B点坐标(xB,yB),其中A和B是车道路口停止线的两个端点;从车辆前边缘向直线AB做垂线,求得垂线的长度length用length近似替代VID;
AB在二维坐标系中满足公式(22):
(yA-yB)·x+(xB-xA)·y+(yB·xA-xB·yA)=0 (22)
假设此时车辆位置P(t)=(xC,yC),则车辆到AB的距离length满足公式(23):
对数据集中每辆车各时刻点的数据做公式(23)运算,提取得到特征VID;
·步骤212,周围车辆运行状态集特征提取
从数据集的地图中获取车辆所在车道是否有左侧车道和右侧车道以及左侧车道与右侧车道的方向向量,同时,在数据集中,不在十字路口的车辆各时刻的车道都有标注,从而得到各时刻各车道上的车辆状态集;当前车道前后车辆通过数据集直接获取;左右车道中的前后车辆按照定义16在左右车道中查找,并将其行驶状态记录在集合NCarState中;
·步骤213,路口允许转向动作特征提取
从数据集中获取十字路口信号灯随时间变化的表格,即对于任意十字路口Ii获取sigi(t);从数据集中的地图车道标记中获取车道允许驾驶行为PLD,则IAM特征的提取按照定义19,依次找到每条车道对应的十字路口,然后对sig(t)和PLD做Hadamard乘积得到IAM特征;
·步骤214,标签提取
将车辆驾驶行为分为直行、向左变道、向右变道、路口直行、路口左转、路口右转、路口调头;其中直行、向左变道、向右变道是直行区和路口准备区的驾驶行为;路口直行、路口左转、路口右转是路口区的驾驶行为;
对驾驶行为按规则进行编码;
路口处直行、左转和右转的驾驶行为在数据集中直接获取;直行、向左变道和向右变道通过车辆当前时刻t所处车道lane(t)和未来时刻t+Δt所处车道lane(t+Δt)做对比;如果相同,则为直行;如果lane(t+Δt)属于lane(t)的左侧车道集,则为向左变道,否则为向右变道;根据提取到的驾驶行为以设定的规则进行编码并作为训练样本的标签;
将获取到的特征表示为定义的feature(t)向量,将获取到的训练样本标签定义为label(t),其中t表示时间;采用min-max归一化方法对feature(t)和label(t)进行归一化处理,如公式(24)所示:
其中f表示特征点取值,fmax表示该特征点取值的最大值,fmin表示该特征点取值的最小值,fN表示归一化之后的特征点取值;归一化处理后的feature(t)和label(t)分别记为featureN(t)和labelN(t);featureN(t)按照顺序拼接上labelN(t)即是网络训练样本和测试样本的结构;
步骤22,驾驶行为预测
前面预处理后的训练样本集记作train,train包括训练样本特征集和训练样本标签集,
训练样本特征集记为train_x,表示为
训练样本标签集记为train_y;测试样本集记为test,包括测试样本特征集和测试样本标签集;测试样本特征集记为test_x,测试样本标签集记为test_y;其中train_y与train_x按行一一对应,test_x与test_y按行一一对应,
车辆驾驶行为预测步骤如下:
(1)按照DBN网络结构利用Contrastive Divergence逐层训练RBM网络,学习特征的低维表达方式;
(2)利用学习所得的网络结构初始化BP神经网络,并在最后添加输出层;
(3)训练BP神经网络,利用训练标签集train_y计算训练误差,反向传播,微调每一层的网络参数;
(4)将待预测样本按照网络正向计算得到输出层节点取值;
(5)对输出层节点的取值做softmax回归,并根据在先设定的驾驶行为编码规则将回归后的结果转换为驾驶行为;
步骤23,位置预测模型
为挖掘车辆驾驶行为与未来一段时间车辆位置之间的关系,利用BP神经网络求解回归问题的能力,挖掘车辆行驶环境和驾驶行为对未来一段时间车辆位置的影响,并利用训练后的BP神经网络结合前面预测所得的驾驶行为进行车辆位置预测:
首先,定义训练样本的特征集为PFeature(t),表示为公式(27):
其中DrvB(t)表示t时刻车辆的驾驶行为经编码后的取值;P(t-1)代表车辆在前一时刻的位置;利用特征提取方法提取每个车辆节点在各时刻对应的PFeature(t),并利用min-max归一化方法对每条样本归一化,得出BP神经网络的输入样本集表示如式(28):
其中PFeatureN(ti)表示ti时刻PFeature的归一化样本;提取车辆节点在n秒后的地理位置P(t+n),1≤n≤3,n∈N,同样对其进行min-max归一化,得到BP神经网络的训练标签集表示如式(29):
然后,设置BP神经网络的输出层为线性输出,按照参数给定的网络结构和激活函数,训练BP神经网络;
最后,在需要进行位置预测时,先根据驾驶行为预测方法,得到车辆驾驶行为预测结果:预测结果是针对softmax分类器编码过的向量形式,将其按照一致的编码规则映射关系,重新编码为标量取值,作为DrvB特征值;同时提取对应的PFeature中其它特征值,做归一化处理,得到待预测样本p_x,样本结构同PFeatureN(t);将待预测样本按照训练所得的神经网络正向计算,得到的网络输出逆归一化结果即是预测的车辆位置;
在驾驶行为预测方法的基础上,结合BP神经网络正向传递的计算方法,车辆位置预测模型表示为递推式(30)的形式:
其中P表示预测模型给出的车辆位置预测结果,即未来时刻车辆节点的位置;
h(·)表示线性函数h(x)=x,f(·)代表sigmoid函数;代表本节训练所得BP神经网络第i层对第i+1层的连接权值,为第i层对第i+1层的偏置项,xp表示PFeatureN形式的样本。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于同济大学,未经同济大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201710702922.0/1.html,转载请声明来源钻瓜专利网。