CN115512543B - 一种基于深度逆向强化学习的车辆路径链重构方法 - Google Patents
一种基于深度逆向强化学习的车辆路径链重构方法 Download PDFInfo
- Publication number
- CN115512543B CN115512543B CN202211152406.2A CN202211152406A CN115512543B CN 115512543 B CN115512543 B CN 115512543B CN 202211152406 A CN202211152406 A CN 202211152406A CN 115512543 B CN115512543 B CN 115512543B
- Authority
- CN
- China
- Prior art keywords
- road section
- path
- time
- chain
- path chain
- 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 48
- 230000002787 reinforcement Effects 0.000 title claims abstract description 34
- 230000002441 reversible effect Effects 0.000 title claims abstract description 31
- 230000006870 function Effects 0.000 claims abstract description 65
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 22
- 230000000295 complement effect Effects 0.000 claims abstract description 14
- 238000007781 pre-processing Methods 0.000 claims abstract description 9
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 38
- 230000009471 action Effects 0.000 claims description 26
- 230000008569 process Effects 0.000 claims description 14
- 238000000354 decomposition reaction Methods 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 9
- 238000007689 inspection Methods 0.000 claims description 8
- 238000005065 mining Methods 0.000 claims description 5
- 230000007704 transition Effects 0.000 claims description 5
- 230000002159 abnormal effect Effects 0.000 claims description 4
- 230000005856 abnormality Effects 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 210000002569 neuron Anatomy 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 230000002459 sustained effect Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 239000000284 extract Substances 0.000 claims description 2
- 239000011521 glass Substances 0.000 claims description 2
- 238000011478 gradient descent method Methods 0.000 claims description 2
- 238000002360 preparation method Methods 0.000 claims description 2
- 230000007547 defect Effects 0.000 abstract description 3
- 230000000694 effects Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 3
- 238000007667 floating Methods 0.000 description 3
- 238000004140 cleaning Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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
- G08G1/0125—Traffic data processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/21—Design, administration or maintenance of databases
- G06F16/215—Improving data quality; Data cleansing, e.g. de-duplication, removing invalid entries or correcting typographical errors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/62—Text, e.g. of license plates, overlay texts or captions on TV images
- G06V20/625—License plates
-
- 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/0137—Measuring and analyzing of parameters relative to traffic conditions for specific applications
-
- 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/0137—Measuring and analyzing of parameters relative to traffic conditions for specific applications
- G08G1/0145—Measuring and analyzing of parameters relative to traffic conditions for specific applications for active traffic flow control
-
- 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
- G08G1/0175—Detecting movement of traffic to be counted or controlled identifying vehicles by photographing vehicles, e.g. when violating traffic rules
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/09—Arrangements for giving variable traffic instructions
- G08G1/0962—Arrangements for giving variable traffic instructions having an indicator mounted inside the vehicle, e.g. giving voice messages
- G08G1/0968—Systems involving transmission of navigation instructions to the vehicle
- G08G1/096833—Systems involving transmission of navigation instructions to the vehicle where different aspects are considered when computing the route
-
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Mathematical Physics (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Computing Systems (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Multimedia (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Quality & Reliability (AREA)
- Artificial Intelligence (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明提供一种基于深度逆向强化学习的车辆路径链重构方法,涉及智慧交通领域。首先对车牌识别数据进行预处理,并提取出路段行程时间,利用基于辅助信息的张量补全算法对路段行程时间补全,从而分离出完整路径链和缺失路径链。然后利用深度逆向强化学习,对路网中完整出行路径链进行挖掘,以非线性回报函数的形式拟合出潜藏的路径选择特性,指导智能体自主重构缺失路径链。本发明克服了现有算法决策依据主观性大的缺陷,在少量的示例数据下就能达到稳定性强、准确度高的重构效果,为交通需求结构分析和交通拥堵疏导等重点交通问题提供数据支撑。
Description
技术领域
本发明涉及智慧交通领域,特别是涉及一种基于深度逆向强化学习的车辆路径链重构方法。
背景技术
目前,交通信息采集技术主要包括线圈、微波、浮动车、GPS和车牌识别等方法,线圈和微波能记录下某一点的交通信息,但无法与车辆联系起来;浮动车、GPS数据虽然分布范围广,但其仅依靠装有GPS定位的出租、公交等车辆,车辆数据来源较为单一,而车牌识别数据拥有能够识别身份信息、高覆盖率、大样本量等特点,能更好的代表城市路网交通模式,具有广泛的应用前景。
随着城市交通监测系统的完善,大量车牌识别数据被存储下来,为研究城市交通提供了新机遇。但在信息采集过程中受制于成本和技术,采集的车牌数据在时空上存在着不连续性,难以还原城市路网的交通运行状态,无法满足目前智慧交通建设的需要。目前基于车牌识别数据的路径链重构方法很少考虑出行者的路径选择偏好,以至于难以应对复杂的路径缺失情况,无法为交通需求结构分析和交通拥堵疏导等重点交通问题提供数据支撑。
发明内容
本发明为克服现有技术的不足,提供一种基于深度逆向强化学习的车辆路径链重构方法。
本发明所采用的具体技术方案如下:
本发明提供了一种基于深度逆向强化学习的车辆路径链重构方法,具体如下:
步骤1:对采集到的车牌识别原始数据进行预处理,以清洗异常数据,得到车牌识别数据;
步骤2:对所述车牌识别数据按车辆经行时间排列为经行卡口链,计算相邻经行卡口间的时间差;划分时间窗,对各路段各时间窗内的时间差进行统计并作为样本量;若样本量满足阈值,则剔除异常后将样本均值作为该时间窗内的路段行程时间,反之,则视为路段行程时间缺失;
步骤3:对于步骤2中缺失的路段行程时间,采用基于辅助信息的张量补全算法进行路段行程时间补全;
步骤4:对于步骤2中原本完整的路段行程时间或经步骤3补全后的路段行程时间,对经行卡口链进行拓扑检查和行程时间阈值检查,分离出缺失路径链集合与完整路径链集合;
步骤5:将城市路网上的路径链重构问题建模为回报函数未知的马尔可夫决策过程,基于最大熵逆向强化学习,对所述完整路径链集合进行挖掘,以求解最佳回报函数;
步骤6:基于所述最佳回报函数,采用Q学习算法求解路径重构的最优策略,指导智能体进行所述缺失路径链集合重构,得到最终的路径重构方案。
作为优选,所述步骤1中,对于车牌识别原始数据中,由于车牌不明或设备漏检导致的无效数据,以及由于外界干扰或设备故障造成的错误数据,均采用直接剔除的方式进行车牌识别数据预处理;
对于车牌识别原始数据中,由于设备故障造成多条过车记录所有字段完全一致的重复数据,采用直接保留最后一条记录的方式进行车牌识别数据预处理;
对于车牌识别原始数据中,由于多检错拍或广角错拍造成车牌或卡口名称一致、记录时间稍有差别的重复数据,将重复记录按时序排列,计算连续两条记录之间的过车时间差;若过车时间差小于重复检测时间阈值,则被视为重复记录,采用保留后一条记录的方式进行车牌识别数据预处理,反之,两条记录均予以保留。
作为优选,所述步骤2具体如下:
根据所述车牌识别数据,将车辆一天的记录按经行时间排列为经行卡口链,当所述经行卡口链中相邻节点(j,j+1)拓扑相连,则计算其节点间的时间差;以5min为一个时间窗,统计各个路段各时间窗内的路段行程时间;若时间窗内样本量≤n,则样本不具代表性,视为路段行程时间缺失;若样本数量>n,则利用箱型图法进行异常值分析,异常值剔除后,将时间窗内的样本均值作为该时间窗内的路段行程时间。
作为优选,所述步骤3具体如下:
构建三个维度的相似矩阵作为辅助信息并加入张量补全模型,通过交替更新因子矩阵法对该优化问题进行求解,获得路段行程时间的最优补全张量;
所述三个维度的相似矩阵包括:
路段相似矩阵M1:基于路段属性矩阵,利用余弦相似度计算路段相似矩阵;所述路段属性包括道路拓扑属性和路段特征属性,道路拓扑属性包括路段入度、出度、邻接度以及路段特征属性,路段特征属性包括路段长度、路段车道数、路段等级和路段兴趣点数;
时间窗相似矩阵M2:基于车速属性,利用余弦相似度计算时间窗相似矩阵;
天数相似矩阵M3:基于车速和天气属性,利用余弦相似度计算天数相似矩阵;
所述余弦相似度的计算公式如下:
式中,cos_pn,n+1为对象n与n+1的余弦相似度,bn,j为对象n的第j个属性值,bn+1,j为对象n+1的第j个属性值;
加入辅助信息的所述张量补全模型公式如下:
式中,T为原始张量;为补全张量,基于tucker分解,/>其中G为tucker分解后的核心张量,U,V,W为分解后的因子矩阵;×n为张量与矩阵的模态积;||·||2 F为F-范数;L(·)为正则化项;M1,M2,M3为相似矩阵。
作为优选,所述步骤4具体如下:
根据路段行程时间ti,o,k,构建路段行程时间的阈值 其中,ti,o,k代表i路段第k天第o时间窗的行程时间,δ为调整系数;对于任意一条经行卡口链,按时间顺序对其每一组相邻节点(j,j+1)进行遍历;若每组节点均拓扑相连,且节点时间差满足(ti,min,ti,max),则将经行卡口链放入完整路径链集合中;若在某组中时间差不满足(ti,min,ti,max),则从(j,j+1)处分离,j节点前的经行卡口链放入完整路径链集合,从j+1节点处检验下一组相邻节点;若相邻节点拓扑不相连,以路段长度得出两个节点间的最短路径,计算节点间最短路径的时间总阈值(tsp,min,tsp,max),通过累加最短路径中每个路段的(ti,min,ti,max)即可;若节点时间差满足(tsp,min,tsp,max),则将(j,j+1)放入缺失路径链集合,继续检验j+1节点后的卡口链;若节点时间差不满足(tsp,min,tsp,max),则在(j,j+1)处分离,j节点前的卡口链放入完整路径链集合,从j+1节点继续向后检验。
作为优选,所述步骤5中马尔可夫决策过程的构建方法具体如下:
a)令环境E为路网;
b)令状态空间S为智能体可以到达的所有路段,每一个路段状态通过多个状态特征表示;
c)令动作空间A包含左转、右转、直行和掉头;
d)令状态转移T为智能体在执行动作a后从当前路段转向下一路段;
e)令专家示例D′为从完整路径链集合中抽取部分作为深度逆向强化学习的专家示例;
f)令策略π为状态到动作的映射学习,即智能体如何选择动作的决策过程;
g)回报函数R用于评价智能体策略的优劣,深度逆向强化学习的目标是使拟合出的回报函数与专家示例的回报函数相近。
进一步的,所述步骤b)中,路段状态特征包括路段长度、路段等级、路段中心经纬度、行程时间、行程时间标准差以及端点是否为信号交叉口。
进一步的,所述回报函数的求解具体如下:
利用全连接网络DNN来拟合回报函数,回报值r表示为路段状态特征的非线性函数,其公式如下所示:
r(s)=θT·fet(s)
式中,fet(s)代表路段s的状态特征;θ是回报函数的参数,也是DNN网络中所有神经元所关联的权值和偏置值;
从完整路径链集合中抽取部分路径链d′i,构造一组专家示例D′={d′1,d′2,…,d′n-1,d′n};利用最大熵理论将回报函数求解转化为最小化问题,利用梯度下降法对最小化问题求解,获得最佳回报函数,即挖掘得到路径效用函数;
利用最大熵理论将回报函数求解转化为最小化问题,公式如下:
式中,在深度逆向强化学习过程中会产生临时回报函数rθ(s),di是智能体在rθ(s)下执行出的样本路径链,f表示当前样本路径链对应所有状态的特征期望之和,表示专家示例{d′1,d′2,…,d′n-1,d′n}的特征期望。
作为优选,所述步骤6具体如下:
设定所述缺失路径链集合的起终点,在最佳回报函数已知时,在Q学习算法中设定Q(s,a)为状态-动作值函数,以计算出任意路段状态下的最佳重构动作,即生成重构策略π,以指导强化学习智能体做出决策,得到路径链重构的最优方案。
进一步的,采用Q学习算法求解路径重构的最优策略具体如下:
初始化所有状态动作对的Q(s,a);
在当前状态s下基于当前Q(s,·),选择动作a;
执行动作a,获得下一个状态s′和回报r;
更新Q(s,a),当多次迭代后,Q(s,a)收敛,则认为获得了最佳的Q(s,a);
其中,基于Bellman方程的Q(s,a)更新公式如下:
Qnew(s,a)=Q(s,a)+α[r(s,a)+γmaxQ′(s′,a′)-Q(s,a)]
式中,Qnew(s,a)为更新后的Q值;Q(s,a)为当前Q值;α是学习率;γ是折扣因子,用于权衡近期回报和远期回报;maxQ′(s′,a′)是下一步状态Q′(s′,·)的最大值。
在持续的迭代后,Q最终会收敛于最优动作-状态函数Q*,此时获得最优策略π*=argmaxQ*。
本发明相对于现有技术而言,具有以下有益效果:
本发明首先对车牌识别数据进行预处理,并提取出路段行程时间,利用基于辅助信息的张量补全算法对路段行程时间补全,从而分离出完整路径链和缺失路径链。然后利用深度逆向强化学习,对路网中完整出行路径链进行挖掘,以非线性回报函数的形式拟合出潜藏的路径选择特性,指导智能体自主重构缺失路径链。本方法克服了现有算法决策依据主观性大的缺陷,在少量的示例数据下就能达到稳定性强、准确度高的重构效果,为交通需求结构分析和交通拥堵疏导等重点交通问题提供数据支撑。
附图说明
图1为本发明提出的路径重构方法的主要步骤流程图;
图2为本发明中回报函数的全连接神经网络结构示意图;
图3为本发明实施例1采用的路网示意图;
图4为本发明实施例1采用的路网缺失情况模拟示意图;
图5为本发明实施例1中回报函数求解的路网回报值分布图;
图6为本发明实施例1中不同专家示例数量下的重构准确度示意图;
附图仅用于示例性说明,不能理解为对本专利的限制;为了更好说明本实施例,附图某些元素会有省略、放大或缩小,不能理解为对本专利的限制。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步阐述和说明。本发明中各个实施方式的技术特征在没有相互冲突的前提下,均可进行相应组合。
由于未配备识别设备或设备故障会导致出行信息的丢失,对路网中信息丢失的交叉口定义为虚拟采集卡口,正常工作的交叉口定义为实际采集卡口。将路径链重构问题定义为:当一条路径链存在两个路段拓扑不相连,那么该条路径链存在缺失,将这两条路段作为缺失起终点,找出其间经过的路段,形成完整的出行路径链。如图1所示,对于缺失的路径链,本发明提供了如下重构方法:
步骤1:车牌识别数据预处理。对于采集到的车牌识别原始数据,进行车牌识别数据预处理,清洗存在的异常数据,进而得到车牌识别数据。
在实际应用时,对于车牌识别原始数据,存在以下几种异常类型,如下表1所示:
表1
由于车牌不明、设备漏检导致的无效数据,以及由于外界干扰、设备故障造成的错误数据,均采用直接剔除的方式;由于设备故障造成多条过车记录所有字段完全一致的重复数据,采用直接保留最后一条记录的方式;由于多检错拍、广角错拍造成车牌、卡口名称等一致,记录时间稍有差别的重复数据,将重复记录按时序排列,计算连续两条记录之间的过车时间差;若过车时间差小于重复检测时间阈值(一般取5s),则被视为重复记录,保留后一条记录;反之,两条记录均予以保留。
步骤2:路段行程时间提取。对车牌识别数据按车辆经行时间排列为经行卡口链,计算相邻经行卡口间的时间差。划分时间窗,统计各路段各时间窗内的时间差,对时间窗内样本量满足阈值的,剔除异常后将样本均值作为该时间窗内的路段行程时间,反之,视为缺失。
由于在路径链提取中,路段行程时间不仅是划分单次路径链的重要指标,也是后续重构算法的重要特征。在实际应用时,将车辆一天的记录按时间排列为卡口链,当卡口链中相邻节点(j,j+1)拓扑相连,则计算其节点间的时间差。以5min为一个时间窗,统计各个路段各时间窗内的路段行程时间。若时间窗内样本量小于等于n,一般取3,则样本不具代表性,视为缺失;若样本数量大于n,利用箱型图法进行异常值分析,异常值剔除后,将时间窗内的样本均值作为该时间窗内的路段行程时间。
步骤3:路段行程时间补全。针对缺失的路段行程时间,采用基于辅助信息的张量补全算法进行路段行程时间补全。
在实际应用时,将路段行程时间构建为一个具有I×O×K个元素的三阶张量T,I、O、K分别代表路段、时间窗、天数三个维度,ti,o,k代表i路段第k天第o时间窗的行程时间。构建三个维度的相似矩阵作为辅助信息加入张量补全模型,通过交替更新因子矩阵法对该优化问题进行求解,获得路段行程时间的最优补全张量。
其中,三个维度的相似矩阵包括:
路段相似矩阵M1:基于路段属性矩阵,包含道路拓扑属性(如路段入度、出度、邻接度)以及路段特征属性(如路段长度、路段车道数、路段等级、路段兴趣点数(Point ofInterest,POI)),利用余弦相似度计算路段相似矩阵;
时间窗相似矩阵M2,基于车速属性利用余弦相似度计算时间窗相似矩阵;
天数相似矩阵M3:基于车速-天气属性利用余弦相似度计算天数相似矩阵;
余弦相似度计算公式如下:
式中,cos_pn,n+1为对象n与n+1的余弦相似度,bn,j为对象n的第j个属性值,bn+1,j为对象n+1的第j个属性值。
加入辅助信息的张量补全模型公式如下:
式中,T为原始张量;为补全张量,基于tucker分解,/>其中G为tucker分解后的核心张量,U,V,W为分解后的因子矩阵;×n为张量与矩阵的模态积(n-mode product);||·||2 F为F-范数;L(·)为正则化项;M1,M2,M3为相似矩阵。
步骤4:路径链分离。基于补全后的路段行程时间,对卡口链进行拓扑检查和行程时间阈值检查,分离出缺失路径链集合与完整路径链集合。
在实际应用时,根据路段行程时间ti,o,k,构建路段行程时间的阈值δ为调整系数,一般为1.15,以应对行程时间的区间波动。对于任意一条卡口链,按时间顺序对其每一组相邻节点(j,j+1)进行遍历。若每组节点均拓扑相连,且节点时间差满足(ti,min,ti,max),将卡口链放入完整路径链集合中;若在某组中时间差不满足(ti,min,ti,max),则从(j,j+1)处分离,j节点前的卡口链放入完整路径链集合,从j+1节点处检验下一组相邻节点;若相邻节点拓扑不相连,以路段长度得出两个节点间的最短路径,计算节点间最短路径的时间总阈值(tsp,min,tsp,max),通过累加最短路径中每个路段的(ti,min,ti,max)即可,若节点时间差满足(tsp,min,tsp,max),则将(j,j+1)放入缺失路径链集合,继续检验j+1节点后的卡口链;若节点时间差不满足(tsp,min,tsp,max),则在(j,j+1)处分离,j节点前的卡口链放入完整路径链集合,从j+1节点继续向后检验。
步骤5:回报函数求解。将城市路网上的路径链重构问题建模为回报函数未知的马尔可夫决策过程,基于最大熵逆向强化学习,对完整路径链进行挖掘,求解最佳回报函数,即每个路段的潜在效用函数。
在实际应用时,路径链重构问题建模包括:
将城市路网上的路径链重构建模为一个马尔可夫决策过程(MDP),其设定如下:a)环境E:路网;b)状态空间S:智能体可以到达的所有路段,每一个路段状态可以通过多个状态特征表示;c)动作空间A:包含左转、右转、直行和掉头;d)状态转移T:智能体在执行动作a后从当前路段转向下一路段;e)专家示例D′:从完整路径链集合中抽取部分作为深度逆向强化学习的专家示例;f)策略π:状态到动作的映射学习,即智能体如何选择动作的决策过程;g)回报函数R:用于评价智能体策略的优劣,深度逆向强化学习的目标是使拟合出的回报函数与专家示例的“回报函数”相近;
其中,路段状态特征包括路段长度、路段等级、路段中心经纬度、行程时间、行程时间标准差、端点是否为信号交叉口,如下表2所示:
表2
具体的,如图2所示,回报函数求解包括:
利用全连接网络DNN来拟合回报函数,回报值r表示为路段状态特征的非线性函数,其公式如下所示:
r(s)=θT·fet(s)
式中,fet(s)代表路段s的状态特征,θ是回报函数的参数,也是DNN网络中所有神经元所关联的权值(weight)和偏置值(bias);
从完整路径链集合中抽取部分路径链d′i,构造一组专家示例D′={d′1,d′2,…,d′n-1,d′n},从这组专家示例中学会未知的回报函数。利用最大熵理论将回报函数求解转化为最小化问题,即挖掘出了路径效用函数。
利用最大熵理论将回报函数求解转化为最小化问题,
式中,在深度逆向强化学习过程中会产生临时回报函数rθ(s),di是智能体在rθ(s)下执行出的样本路径链,f表示当前样本路径链对应所有状态的特征期望之和,表示示例路径链{d′1,d′2,…,d′n-1,d′n}的特征期望。
步骤6:路径链重构。基于最佳回报函数,采用Q学习算法求解路径重构的最优策略,指导智能体进行路径链重构,得到最终的路径重构方案。
在实际应用时,路径链重构包括:
设定缺失路径链的起终点,在回报函数已知时,利用Q学习算法求解路径重构的最优策略。Q学习算法中设定Q(s,a)为状态-动作值函数,记录着每种状态下采取每种动作的最大累积期望回报,用来评价在任意路段状态下采取某个重构动作的优劣。通过Q(s,a)可以知道特定状态下的最佳动作,由此生成策略π,指导强化学习智能体做出决策;
利用Q学习算法求解路径重构问题的具体步骤如下:
初始化所有状态动作对的Q(s,a);
在当前状态s下基于当前Q(s,·),选择动作a;
执行动作a,获得下一个状态s′和回报r;
更新Q(s,a),当多次迭代后,Q(s,a)收敛,则认为获得了最佳的Q(s,a);
其中,基于Bellman方程的Q(s,a)更新公式如下:
Qnew(s,a)=Q(s,a)+α[r(s,a)+γmaxQ′(s′,a′)-Q(s,a)]
式中,Qnew(s,a)为更新后的Q值;Q(s,a)为当前Q值;α是学习率;γ是折扣因子,用于权衡近期回报和远期回报;maxQ′(s′,a′)是下一步状态Q′(s′,·)的最大值。
在持续的迭代后,Q最终会收敛于最优动作-状态函数Q*,此时获得最优策略π*=argmaxQ*。在最优策略π*指导下,智能体可以快速重构出最优路径链。
实施例1
本实施例将本发明提出的车辆路径链重构方法用于图3所示的局部路网。该路网共包含43个交叉口,设有采集设备的实际采集卡口33个,剩余10个未设置采集设备,定义为虚拟采集卡口。实际采集卡口提供了车牌识别数据,记录了车辆牌照、过车时间、卡口编号、过车方向等字段。
以局部路网的车牌识别数据为基础,执行发明提出的基于深度逆向强化学习的车辆路径链重构方法,步骤如下:
车牌识别数据预处理。对于采集到的车牌识别原始数据,进行车牌识别数据预处理,清洗存在的异常数据;
路段行程时间提取。将车辆一天的记录按时间排列为卡口链,计算相邻卡口间的时间差。以5min划分时间窗,统计各路段各时间窗内的时间差,对时间窗内样本量满足大于3的,剔除异常后将样本均值作为该时间窗内的路段行程时间,反之,视为缺失。
路段行程时间补全。构建了三阶路段行程时间张量模型,并根据余弦相似度计算出三个维度的相似矩阵作为辅助信息,采用基于辅助信息的张量补全模型进行路段行程时间补全,通过交替更新因子矩阵法求解出最优路段行程时间张量。
路径链分离。基于补全后的路段行程时间,对卡口链进行拓扑检查和行程时间阈值检查,分离出缺失路径链集合与完整路径链集合。
缺失情况模拟。从完整路径链中随机取出若干条路段数大于3的路径链进行演示。演示中以从路网上删除某些实际采集卡口的方式来模拟缺失情况,将完整路径链中经过与该卡口相连的路段删除,作为缺失路径链;以完整路径链中的其他路径链作为专家示例。抽取了路网中六个节点实施缺失情况模拟,如图4所示。
回报函数求解。将城市路网上的路径链重构问题建模为回报函数未知的马尔可夫决策过程,基于最大熵逆向强化学习,对完整路径链进行挖掘,求解最佳回报函数,获得每个路段的潜在效用函数;图5显示了缺失六个节点情况下,求解出的回报值分布图,图中路段的粗细代表回报值的大小。
路径链重构。基于求解出的最佳回报函数,采用Q学习算法求解路径重构的最优策略,指导智能体进行路径链重构,得到最终的路径重构方案。
图6显示了缺失六个节点情况下,经本发明提出的路径链重构方法在200条、400条、600条、800条、1000条专家示例数量下的算法准确度,在专家示例数量小于1000条时,算法准确度能始终保持在90%以上,证明本方法在较小样本量下就能实现对路网出行路径决策的深入挖掘,实现出行路径链的准确重构。
其中,算法准确度公式如下:
式中,pub(di,di,true)为第i条重构路径链di与实际路径链di,true的最长连续公共子串数,即相同的连续路段数;sum(di,true)代表实际路径链的路段数;m为重构路径链总数。
以上所述的实施例只是本发明的一种较佳的方案,然其并非用以限制本发明。有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型。因此凡采取等同替换或等效变换的方式所获得的技术方案,均落在本发明的保护范围内。
Claims (10)
1.一种基于深度逆向强化学习的车辆路径链重构方法,其特征在于,具体如下:
步骤1:对采集到的车牌识别原始数据进行预处理,以清洗异常数据,得到车牌识别数据;
步骤2:对所述车牌识别数据按车辆经行时间排列为经行卡口链,计算相邻经行卡口间的时间差;划分时间窗,对各路段各时间窗内的时间差进行统计并作为样本量;若样本量满足阈值,则剔除异常后将样本均值作为该时间窗内的路段行程时间,反之,则视为路段行程时间缺失;
步骤3:对于步骤2中缺失的路段行程时间,采用基于辅助信息的张量补全算法进行路段行程时间补全;
步骤4:对于步骤2中原本完整的路段行程时间或经步骤3补全后的路段行程时间,对经行卡口链进行拓扑检查和行程时间阈值检查,分离出缺失路径链集合与完整路径链集合;
步骤5:将城市路网上的路径链重构问题建模为回报函数未知的马尔可夫决策过程,基于最大熵逆向强化学习,对所述完整路径链集合进行挖掘,以求解最佳回报函数;
步骤6:基于所述最佳回报函数,采用Q学习算法求解路径重构的最优策略,指导智能体进行所述缺失路径链集合重构,得到最终的路径重构方案。
2.根据权利要求1所述的一种基于深度逆向强化学习的车辆路径链重构方法,其特征在于,所述步骤1中,对于车牌识别原始数据中,由于车牌不明或设备漏检导致的无效数据,以及由于外界干扰或设备故障造成的错误数据,均采用直接剔除的方式进行车牌识别数据预处理;
对于车牌识别原始数据中,由于设备故障造成多条过车记录所有字段完全一致的重复数据,采用直接保留最后一条记录的方式进行车牌识别数据预处理;
对于车牌识别原始数据中,由于多检错拍或广角错拍造成车牌或卡口名称一致、记录时间稍有差别的重复数据,将重复记录按时序排列,计算连续两条记录之间的过车时间差;若过车时间差小于重复检测时间阈值,则被视为重复记录,采用保留后一条记录的方式进行车牌识别数据预处理,反之,两条记录均予以保留。
3.根据权利要求1所述的一种基于深度逆向强化学习的车辆路径链重构方法,其特征在于,所述步骤2具体如下:
根据所述车牌识别数据,将车辆一天的记录按经行时间排列为经行卡口链,当所述经行卡口链中相邻节点(j,j+1)拓扑相连,则计算其节点间的时间差;以5min为一个时间窗,统计各个路段各时间窗内的路段行程时间;若时间窗内样本量≤n,则样本不具代表性,视为路段行程时间缺失;若样本数量>n,则利用箱型图法进行异常值分析,异常值剔除后,将时间窗内的样本均值作为该时间窗内的路段行程时间。
4.根据权利要求1所述的一种基于深度逆向强化学习的车辆路径链重构方法,其特征在于,所述步骤3具体如下:
构建三个维度的相似矩阵作为辅助信息并加入张量补全模型,通过交替更新因子矩阵法对优化问题进行求解,获得路段行程时间的最优补全张量;
所述三个维度的相似矩阵包括:
路段相似矩阵M1:基于路段属性矩阵,利用余弦相似度计算路段相似矩阵;所述路段属性包括道路拓扑属性和路段特征属性,道路拓扑属性包括路段入度、出度、邻接度以及路段特征属性,路段特征属性包括路段长度、路段车道数、路段等级和路段兴趣点数;
时间窗相似矩阵M2:基于车速属性,利用余弦相似度计算时间窗相似矩阵;
天数相似矩阵M3:基于车速和天气属性,利用余弦相似度计算天数相似矩阵;
所述余弦相似度的计算公式如下:
式中,cos_pn,n+1为对象n与n+1的余弦相似度,bn,j为对象n的第j个属性值,bn+1,j为对象n+1的第j个属性值;
加入辅助信息的所述张量补全模型公式如下:
式中,T为原始张量;为补全张量,基于tucker分解,/>其中G为tucker分解后的核心张量,U,V,W为分解后的因子矩阵;×n为张量与矩阵的模态积;||·||2 F为F-范数;L(·)为正则化项;M1,M2,M3为相似矩阵。
5.根据权利要求1所述的一种基于深度逆向强化学习的车辆路径链重构方法,其特征在于,所述步骤4具体如下:
根据路段行程时间ti,o,k,构建路段行程时间的阈值 其中,ti,o,k代表i路段第k天第o时间窗的行程时间,δ为调整系数;对于任意一条经行卡口链,按时间顺序对其每一组相邻节点(j,j+1)进行遍历;若每组节点均拓扑相连,且节点时间差满足(ti,min,ti,max),则将经行卡口链放入完整路径链集合中;若在某组中时间差不满足(ti,min,ti,max),则从(j,j+1)处分离,j节点前的经行卡口链放入完整路径链集合,从j+1节点处检验下一组相邻节点;若相邻节点拓扑不相连,以路段长度得出两个节点间的最短路径,计算节点间最短路径的时间总阈值(tsp,min,tsp,max),通过累加最短路径中每个路段的(ti,min,ti,max)即可;若节点时间差满足(tsp,min,tsp,max),则将(j,j+1)放入缺失路径链集合,继续检验j+1节点后的卡口链;若节点时间差不满足(tsp,min,tsp,max),则在(j,j+1)处分离,j节点前的卡口链放入完整路径链集合,从j+1节点继续向后检验。
6.根据权利要求1所述的一种基于深度逆向强化学习的车辆路径链重构方法,其特征在于,所述步骤5中马尔可夫决策过程的构建方法具体如下:
a)令环境E为路网;
b)令状态空间S为智能体可以到达的所有路段,每一个路段状态通过多个状态特征表示;
c)令动作空间A包含左转、右转、直行和掉头;
d)令状态转移T为智能体在执行动作a后从当前路段转向下一路段;
e)令专家示例D′为从完整路径链集合中抽取部分作为深度逆向强化学习的专家示例;
f)令策略π为状态到动作的映射学习,即智能体如何选择动作的决策过程;
g)回报函数R用于评价智能体策略的优劣,深度逆向强化学习的目标是使拟合出的回报函数与专家示例的回报函数相近。
7.根据权利要求6所述的一种基于深度逆向强化学习的车辆路径链重构方法,其特征在于,所述步骤b)中,路段状态特征包括路段长度、路段等级、路段中心经纬度、行程时间、行程时间标准差以及端点是否为信号交叉口。
8.根据权利要求6所述的一种基于深度逆向强化学习的车辆路径链重构方法,其特征在于,所述回报函数的求解具体如下:
利用全连接网络DNN来拟合回报函数,回报值r表示为路段状态特征的非线性函数,其公式如下所示:
r(s)=θT·fet(s)
式中,fet(s)代表路段s的状态特征;θ是回报函数的参数,也是DNN网络中所有神经元所关联的权值和偏置值;
从完整路径链集合中抽取部分路径链d′i,构造一组专家示例D′={d′1,d′2,…,d′n-1,d′n};利用最大熵理论将回报函数求解转化为最小化问题,利用梯度下降法对最小化问题求解,获得最佳回报函数,即挖掘得到路径效用函数;
利用最大熵理论将回报函数求解转化为最小化问题,公式如下:
式中,在深度逆向强化学习过程中会产生临时回报函数rθ(s),di是智能体在rθ(s)下执行出的样本路径链,f表示当前样本路径链对应所有状态的特征期望之和,表示专家示例{d′1,d′2,…,d′n-1,d′n}的特征期望。
9.根据权利要求1所述的一种基于深度逆向强化学习的车辆路径链重构方法,其特征在于,所述步骤6具体如下:
设定所述缺失路径链集合的起终点,在最佳回报函数已知时,在Q学习算法中设定Q(s,a)为状态-动作值函数,以计算出任意路段状态下的最佳重构动作,即生成重构策略π,以指导强化学习智能体做出决策,得到路径链重构的最优方案。
10.根据权利要求9所述的一种基于深度逆向强化学习的车辆路径链重构方法,其特征在于,采用Q学习算法求解路径重构的最优策略具体如下:
初始化所有状态动作对的Q(s,a);
在当前状态s下基于当前Q(s,·),选择动作a;
执行动作a,获得下一个状态s′和回报r;
更新Q(s,a),当多次迭代后,Q(s,a)收敛,则认为获得了最佳的Q(s,a);
其中,基于Bellman方程的Q(s,a)更新公式如下:
Qnew(s,a)=Q(s,a)+α[r(s,a)+γmaxQ′(s′,a′)-Q(s,a)]
式中,Qnew(s,a)为更新后的Q值;Q(s,a)为当前Q值;α是学习率;γ是折扣因子,用于权衡近期回报和远期回报;maxQ′(s′,a′)是下一步状态Q′(s′,·)的最大值;
在持续的迭代后,Q最终会收敛于最优动作-状态函数Q*,此时获得最优策略π*=argmaxQ*。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211152406.2A CN115512543B (zh) | 2022-09-21 | 2022-09-21 | 一种基于深度逆向强化学习的车辆路径链重构方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211152406.2A CN115512543B (zh) | 2022-09-21 | 2022-09-21 | 一种基于深度逆向强化学习的车辆路径链重构方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115512543A CN115512543A (zh) | 2022-12-23 |
CN115512543B true CN115512543B (zh) | 2023-11-28 |
Family
ID=84504721
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211152406.2A Active CN115512543B (zh) | 2022-09-21 | 2022-09-21 | 一种基于深度逆向强化学习的车辆路径链重构方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115512543B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116071929B (zh) * | 2023-03-06 | 2023-08-01 | 深圳市城市交通规划设计研究中心股份有限公司 | 基于卡口车牌识别数据的实时路况监测系统及其方法 |
CN117711174A (zh) * | 2023-12-07 | 2024-03-15 | 山东高速集团有限公司 | 一种车辆通行信息的数据处理方法及系统 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106408960A (zh) * | 2016-10-29 | 2017-02-15 | 浙江大学 | 基于有序聚类的信号控制时段分割方法 |
CN110276950A (zh) * | 2019-06-24 | 2019-09-24 | 华南理工大学 | 一种基于卡口视频数据的城市交通出行链重构方法 |
CN110362557A (zh) * | 2019-05-27 | 2019-10-22 | 华南理工大学 | 一种基于机器学习和车牌识别数据的缺失路径修复方法 |
WO2020047907A1 (zh) * | 2018-09-07 | 2020-03-12 | 佳都新太科技股份有限公司 | 一种漏检车牌预测方法及处理终端 |
CN111243277A (zh) * | 2020-03-09 | 2020-06-05 | 山东大学 | 基于车牌识别数据的通勤车辆时空轨迹重构方法及系统 |
CN112365711A (zh) * | 2020-10-21 | 2021-02-12 | 东南大学 | 一种基于车牌识别数据的车辆轨迹重构方法 |
CN113140114A (zh) * | 2021-03-09 | 2021-07-20 | 中山大学 | 一种基于行程时间估计的车辆出行路径重构方法 |
CN113734197A (zh) * | 2021-09-03 | 2021-12-03 | 合肥学院 | 一种基于数据融合的无人驾驶的智能控制方案 |
WO2022083166A1 (zh) * | 2020-10-19 | 2022-04-28 | 山东交通学院 | 一种卡口数据缺失下的车辆行车轨迹重建方法及系统 |
-
2022
- 2022-09-21 CN CN202211152406.2A patent/CN115512543B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106408960A (zh) * | 2016-10-29 | 2017-02-15 | 浙江大学 | 基于有序聚类的信号控制时段分割方法 |
WO2020047907A1 (zh) * | 2018-09-07 | 2020-03-12 | 佳都新太科技股份有限公司 | 一种漏检车牌预测方法及处理终端 |
CN110362557A (zh) * | 2019-05-27 | 2019-10-22 | 华南理工大学 | 一种基于机器学习和车牌识别数据的缺失路径修复方法 |
CN110276950A (zh) * | 2019-06-24 | 2019-09-24 | 华南理工大学 | 一种基于卡口视频数据的城市交通出行链重构方法 |
CN111243277A (zh) * | 2020-03-09 | 2020-06-05 | 山东大学 | 基于车牌识别数据的通勤车辆时空轨迹重构方法及系统 |
WO2022083166A1 (zh) * | 2020-10-19 | 2022-04-28 | 山东交通学院 | 一种卡口数据缺失下的车辆行车轨迹重建方法及系统 |
CN112365711A (zh) * | 2020-10-21 | 2021-02-12 | 东南大学 | 一种基于车牌识别数据的车辆轨迹重构方法 |
CN113140114A (zh) * | 2021-03-09 | 2021-07-20 | 中山大学 | 一种基于行程时间估计的车辆出行路径重构方法 |
CN113734197A (zh) * | 2021-09-03 | 2021-12-03 | 合肥学院 | 一种基于数据融合的无人驾驶的智能控制方案 |
Non-Patent Citations (2)
Title |
---|
基于梯度提升决策树的城市车辆路径链重构;徐建闽;魏鑫;林永杰;卢凯;;华南理工大学学报(自然科学版)(第07期);59-68 * |
基于车牌识别数据的机动车出行轨迹提取算法;阮树斌;王福建;马东方;金盛;王殿海;;浙江大学学报(工学版)(第05期);23-31 * |
Also Published As
Publication number | Publication date |
---|---|
CN115512543A (zh) | 2022-12-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115512543B (zh) | 一种基于深度逆向强化学习的车辆路径链重构方法 | |
CN109754605B (zh) | 一种基于注意力时态图卷积网络的交通预测方法 | |
CN112216108B (zh) | 一种基于属性增强的时空图卷积模型的交通预测方法 | |
CN109285346B (zh) | 一种基于关键路段的城市路网交通状态预测方法 | |
CN110276950B (zh) | 一种基于卡口视频数据的城市交通出行链重构方法 | |
CN109670277A (zh) | 一种基于多模态数据融合与多模型集成的旅行时间预测方法 | |
CN112489426B (zh) | 一种基于图卷积神经网络的城市交通流量时空预测方案 | |
CN109785618B (zh) | 一种基于组合逻辑的短时交通流预测方法 | |
CN110223509B (zh) | 一种基于贝叶斯增强张量的缺失交通数据修复方法 | |
CN111243299A (zh) | 一种基于3dqn_pser算法的单交叉口信号控制方法 | |
CN110210169B (zh) | 一种基于lstm的盾构机故障预测方法 | |
CN114881356A (zh) | 基于粒子群算法优化bp神经网络的城市交通碳排预测方法 | |
CN114495507B (zh) | 融合时空注意力神经网络和交通模型的交通流预测方法 | |
CN115410372B (zh) | 基于贝叶斯lstm的高速公路交通流量可靠预测方法 | |
CN113591380A (zh) | 基于图高斯过程的交通流预测方法、介质及设备 | |
CN112633602B (zh) | 一种基于gis地图信息的交通拥堵指数预测方法及装置 | |
CN113537626B (zh) | 一种聚合信息差异的神经网络结合时间序列预测的方法 | |
CN112820104A (zh) | 一种基于时空聚类张量分解的交通数据补全方法 | |
CN111723990B (zh) | 基于双向长短期记忆神经网络的共享单车流量预测方法 | |
CN111311905A (zh) | 一种基于粒子群优化小波神经网络的高速公路行程时间预测方法 | |
CN113420488A (zh) | 一种基于轨迹重构的城市路网od估计方法 | |
CN116739376A (zh) | 一种基于数据挖掘的公路路面预防性养护决策方法 | |
CN110287995B (zh) | 全天高架交通拥堵情况分级的多特征学习网络模型方法 | |
CN116756825A (zh) | 一种中小跨径桥梁的群结构性能预测系统 | |
Amrani et al. | Train speed profiles optimization using a genetic algorithm based on a random-forest model to estimate energy consumption |
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 |