CN115409256B - Route recommendation method for avoiding congestion area based on travel time prediction - Google Patents

Route recommendation method for avoiding congestion area based on travel time prediction Download PDF

Info

Publication number
CN115409256B
CN115409256B CN202211021689.7A CN202211021689A CN115409256B CN 115409256 B CN115409256 B CN 115409256B CN 202211021689 A CN202211021689 A CN 202211021689A CN 115409256 B CN115409256 B CN 115409256B
Authority
CN
China
Prior art keywords
path
travel time
edge
algorithm
network
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
Application number
CN202211021689.7A
Other languages
Chinese (zh)
Other versions
CN115409256A (en
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.)
Jilin Institute of Chemical Technology
Original Assignee
Jilin Institute of Chemical Technology
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 Jilin Institute of Chemical Technology filed Critical Jilin Institute of Chemical Technology
Priority to CN202211021689.7A priority Critical patent/CN115409256B/en
Publication of CN115409256A publication Critical patent/CN115409256A/en
Application granted granted Critical
Publication of CN115409256B publication Critical patent/CN115409256B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G06Q10/047Optimisation of routes or paths, e.g. travelling salesman problem
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/08Learning methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/40Business processes related to the transportation industry
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0125Traffic data processing
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions
    • G08G1/0137Measuring and analyzing of parameters relative to traffic conditions for specific applications
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Marketing (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Biomedical Technology (AREA)
  • Data Mining & Analysis (AREA)
  • Operations Research (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Game Theory and Decision Science (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Quality & Reliability (AREA)
  • Evolutionary Computation (AREA)
  • Development Economics (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Primary Health Care (AREA)
  • Traffic Control Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了基于行程时间预测的拥堵区避让的路径推荐方法,包括基于DGCN的路段行程时间预测模型、基于行程时间预测的路网区域状态识别模型和基于Floyd算法的路径推荐模型三大步,本发明基于行程时间预测及拥堵区避让的路径推荐模型,基于历史卡口数据和路网拓扑结构,使用动态图神经网络(DGCN)对城市路网的路段行程时间进行短时预测,并通过结合路网动态实际状况和流量特征,建立基于加权GN算法的路网分区模型,根据预测的行程时间结果对路网进行区域状态识别;最后结合预测的行程时间和路网的拥堵状态,基于Floyd算法做避让的路径推荐,一定程度上,减缓拥堵路段的拥堵程度,对出行者提供最短路径推荐,提高出行体验感。

Figure 202211021689

The present invention provides a route recommendation method based on travel time prediction to avoid congestion areas, including three steps: a DGCN-based road segment travel time prediction model, a road network area state recognition model based on travel time prediction, and a route recommendation model based on the Floyd algorithm. The present invention is based on the route recommendation model of travel time prediction and congested area avoidance, based on historical checkpoint data and road network topology structure, uses dynamic graph neural network (DGCN) to carry out short-term prediction on road section travel time of urban road network, and combines Based on the dynamic actual conditions and traffic characteristics of the road network, a road network partition model based on the weighted GN algorithm is established, and the regional state of the road network is identified according to the predicted travel time results; finally, combined with the predicted travel time and the congestion status of the road network, based on the Floyd algorithm Recommendations for avoidance routes can, to a certain extent, reduce the degree of congestion on congested roads, provide travelers with shortest route recommendations, and improve travel experience.

Figure 202211021689

Description

基于行程时间预测的拥堵区避让的路径推荐方法Route recommendation method for avoiding congestion areas based on travel time prediction

技术领域Technical Field

本发明属于智能交通技术领域,涉及基于行程时间预测的拥堵区避让的路径推荐方法。The invention belongs to the technical field of intelligent transportation and relates to a path recommendation method for avoiding congestion areas based on travel time prediction.

背景技术Background Art

随着汽车数量逐年递增,部分一线城市,交通拥堵已由中心城区向边缘城区扩散,且车辆出行时间的可靠性随着交通拥堵的不确定性也在降低。路线推荐是当今使用最广泛的基于位置的服务之一,它对于良好的驾驶体验和顺畅的公共交通至关重要,可以根据给定的标准为用户提供最佳旅行体验的路线。因此,出行者和交通管理人员对于车辆出行路径问题也极为关注。As the number of cars increases year by year, traffic congestion in some first-tier cities has spread from the central urban areas to the peripheral urban areas, and the reliability of vehicle travel time is also decreasing with the uncertainty of traffic congestion. Route recommendation is one of the most widely used location-based services today. It is essential for a good driving experience and smooth public transportation. It can provide users with the best travel experience based on given standards. Therefore, travelers and traffic managers are also very concerned about the problem of vehicle travel paths.

车辆出行路径选择问题是智能交通系统研究的重点,动态交通信息的获取,也是行程时间预测和优化路径推荐的前提。杨枫等,针对应急车辆路径选择问题,结合动态路网和急救援建立最大路径可靠性第一目标和最短行程时间为第二目标的两阶段优化模型,并设计混合布谷鸟搜索算法,缩小算法全局搜索空间,提高优化模型全局寻优能力。为降低AGV的运行成本和路径成本,提出一种基于Dijkstra蚂蚁群优化算法(ACO-DA)的机场AGV路径优化模型,考虑了有障碍物的环境,通过蚁群优化进行行李提取排序,并与Dijkstra算法集成进行AGV路径规划。Wang Chaoxiong等;提出了基于A*最快路径推荐的神经网络模型,将用于路线行程时间预测的交通预测模型和用于OD行程时间估计的多任务表示学习模型,结合A*算法,无需实际路径信息即可实现高精度的行程时间的预测,同时找到实时最快的路线。Li Ning为了避免船舶碰撞,基于粒子群优化算法、遗传算法和非线性规划理论,建立了船舶防撞路径规划的优化模型。Liu Binbin;分析了物流配送路径确定方案的研究现状,并应用递归模糊神经网络算法选择电子商务物流配送路径方案。Dengkai HOU针对基于实时交通信息的多路径动态多车场冷藏车辆路径问题,基于预优化后实时调整的思想,建立以总成本最小化为目标的两阶段优化模型,设计了一种具有可变邻域搜索(HCGAVNS)的混合混沌遗传算法来生成初始路径。安冬冬等从驼背运输的实际情况出发,用交通网络图模拟路网中节点与线路的关系,考虑时间成本、费用成本和交通环境影响成本做函数分析;在算法上,将运输成本加到虚拟路线上,优化最短路路径Floyd算法,求得最优路径,为联运企业路径选择提高参考思路。于泉等提出基于时间最短、距离最短、拥挤度最低的多目标路径选择模型,从路径选择和路径规划两方面出发,结合系统中实时反馈的路况信息,实现实时的多目标动态重规划出行路径。Wang Jingyuan等将基于注意的RNN网络建模从从起点到候选位置的成本和使用价值网络来估计从候选位置到目的地的成本两部分进行集成,提出使用神经网络自动学习经典启发式算法的代价函数,将A*算法用于PRR任务中,以获得搜索候选位置的更准确成本,智路平等从微观层面对随机动态路网进行研究,建立路段行程时间可靠性计算模型,针对简单网络和复杂网络的车辆路径选择问题,提出使用行程时间可靠性作为关键控制变量的路径选择算法,提高路径选择效率。张传琪等根据路网中车辆的流量特征,将路网根据拥挤程度分为不同时段,考虑车辆服务的顺序、路径以及最优出行时间等特征,建立动态拥挤的配送路径优化模型,结合改进的遗传算法,求解最小配送成本的目标函数。为解决路径在搜索过程中很难利用有用的上下文信息问题,Ying Shen等,将改进的DBSCAN算法集成到热点提取中,以及基于密度的ε距离轨迹聚类算法来识别推荐的潜在路径,驾驶时间、速度、距离和端点吸引力等因素定义加权树,帮助巡游出租车司机找到具有最佳路线的潜在乘客。The problem of vehicle travel path selection is the focus of intelligent transportation system research. The acquisition of dynamic traffic information is also the premise of travel time prediction and optimized path recommendation. Yang Feng et al., for the problem of emergency vehicle path selection, combined with dynamic road network and emergency rescue, established a two-stage optimization model with maximum path reliability as the first goal and shortest travel time as the second goal, and designed a hybrid cuckoo search algorithm to narrow the global search space of the algorithm and improve the global optimization ability of the optimization model. In order to reduce the operating cost and path cost of AGV, an airport AGV path optimization model based on Dijkstra ant colony optimization algorithm (ACO-DA) was proposed. Considering the environment with obstacles, the baggage collection sorting was carried out through ant colony optimization, and the AGV path planning was integrated with the Dijkstra algorithm. Wang Chaoxiong et al. proposed a neural network model based on A* fastest path recommendation. The traffic prediction model for route travel time prediction and the multi-task representation learning model for OD travel time estimation were combined with the A* algorithm. High-precision travel time prediction can be achieved without actual path information, and the fastest real-time route can be found. In order to avoid ship collision, Li Ning established an optimization model for ship collision avoidance path planning based on particle swarm optimization algorithm, genetic algorithm and nonlinear programming theory. Liu Binbin; analyzed the research status of logistics distribution path determination scheme, and applied recursive fuzzy neural network algorithm to select e-commerce logistics distribution path scheme. Dengkai HOU aimed at the multi-path dynamic multi-yard refrigerated vehicle path problem based on real-time traffic information. Based on the idea of pre-optimization and real-time adjustment, a two-stage optimization model with the goal of minimizing the total cost was established, and a hybrid chaos genetic algorithm with variable neighborhood search (HCGAVNS) was designed to generate the initial path. An Dongdong et al., based on the actual situation of humpback transportation, used the traffic network diagram to simulate the relationship between nodes and routes in the road network, and considered the time cost, expense cost and traffic environment impact cost for function analysis; in terms of algorithm, the transportation cost was added to the virtual route, the shortest path Floyd algorithm was optimized, and the optimal path was obtained, which provided reference ideas for the path selection of intermodal enterprises. Yu Quan et al. proposed a multi-objective path selection model based on the shortest time, shortest distance and lowest congestion. Starting from the two aspects of path selection and path planning, combined with the real-time feedback of road conditions in the system, real-time multi-objective dynamic re-planning of travel paths is realized. Wang Jingyuan et al. integrated the cost from the starting point to the candidate location and the cost from the candidate location to the destination using the value network based on attention RNN network modeling, proposed to use neural networks to automatically learn the cost function of the classic heuristic algorithm, and used the A* algorithm in the PRR task to obtain a more accurate cost for searching candidate locations. Zhilu Pingping studied random dynamic road networks from a microscopic level, established a road section travel time reliability calculation model, and proposed a path selection algorithm using travel time reliability as a key control variable for the vehicle path selection problem of simple and complex networks to improve the efficiency of path selection. Zhang Chuanqi et al. divided the road network into different time periods according to the degree of congestion based on the traffic characteristics of vehicles in the road network, considered the characteristics of vehicle service order, path and optimal travel time, established a dynamic congested distribution path optimization model, and combined with an improved genetic algorithm to solve the objective function of the minimum distribution cost. To solve the problem that it is difficult to utilize useful contextual information during the path search process, Ying Shen et al. integrated the improved DBSCAN algorithm into hotspot extraction and the density-based ε-distance trajectory clustering algorithm to identify recommended potential paths. Factors such as driving time, speed, distance, and endpoint attractiveness defined a weighted tree to help cruising taxi drivers find potential passengers with the best route.

综上,对于城市路网车辆路径选择问题,学者多从行程时间、出行成本或者路网行程时间可靠性的角度研究,对动态路网行程时间的预测和通过最短路径算法中求得的最短路径不一定最优等问题考虑不足,本发明充分考虑了中小型城市路网交通的动态行程时间的实时变化,通过加权GN算法对路网按流量加权进行分区,并用动态图神经网络对行程时间进行预测,根据预测所得的行程时间结果结合平均行程速度,进行区域状态识别,最后根据预测的行程时间和拥堵区识别结果,采用Floyd算法进行避让拥堵区的路径推荐。In summary, for the problem of vehicle path selection in urban road networks, scholars have mostly studied it from the perspectives of travel time, travel cost or road network travel time reliability, and have not given enough consideration to the prediction of dynamic road network travel time and the fact that the shortest path obtained by the shortest path algorithm is not necessarily the optimal one. The present invention fully considers the real-time changes of the dynamic travel time of road network traffic in small and medium-sized cities, partitions the road network by traffic weights through the weighted GN algorithm, and predicts the travel time with a dynamic graph neural network. The regional status is identified based on the predicted travel time results combined with the average travel speed. Finally, based on the predicted travel time and the congested area identification results, the Floyd algorithm is used to recommend paths that avoid congested areas.

发明内容Summary of the invention

本发明的目的在于提供基于行程时间预测的拥堵区避让的路径推荐方法,以解决上述背景技术中提出的问题。The purpose of the present invention is to provide a route recommendation method for avoiding congestion areas based on travel time prediction, so as to solve the problems raised in the above background technology.

本发明的目的可通过下列技术方案来实现:基于行程时间预测的拥堵区避让的路径推荐方法,其特征在于,包括基于DGCN的路段行程时间预测模型、基于行程时间预测的路网区域状态识别模型和基于Floyd算法的路径推荐模型三大步,具体如下:The object of the present invention can be achieved by the following technical scheme: a route recommendation method for avoiding congested areas based on travel time prediction, characterized in that it includes three steps: a road section travel time prediction model based on DGCN, a road network area state recognition model based on travel time prediction, and a route recommendation model based on Floyd algorithm, which are specifically as follows:

基于DGCN的路段行程时间预测模型包括拉普拉斯矩阵潜网络模块和基于图卷积网络的交通预测模块,其中:拉普拉斯矩阵预测单元包含三部分:The road section travel time prediction model based on DGCN includes a Laplace matrix latent network module and a traffic prediction module based on a graph convolutional network, where the Laplace matrix prediction unit consists of three parts:

特征采样:对每天最近邻15min、30min、60min的数据进行采样,在减少特征的数据维度的同时对最近的交通进行采样;空间注意力机制:为建立动态的交通路网的空间关系,采用注意力机制对路网当前近邻时间的邻接矩阵进行估计;LSTM单元:采用LSTM来学习时间相关,探索邻接矩阵Ld的序列之间的内在关联;Feature sampling: Sampling the data of the nearest neighbor 15min, 30min, and 60min every day, sampling the nearest traffic while reducing the data dimension of the feature; Spatial attention mechanism: In order to establish the spatial relationship of the dynamic traffic road network, the attention mechanism is used to estimate the adjacency matrix of the current neighbor time of the road network; LSTM unit: LSTM is used to learn time correlation and explore the intrinsic correlation between the sequences of the adjacency matrix Ld;

基于GCN的交通预测模块包含两部分:The GCN-based traffic prediction module consists of two parts:

(1)时间卷积层:主要包含四个卷积核为1x3的二维卷积层,用来提取数据中高维局部时间,交通数据段X(1:k*T)=(X1........Xt......Xk*T)的时间卷积为:TC=Conv1×t(X(1:k*T));(1) Temporal convolution layer: It mainly includes four two-dimensional convolution layers with a convolution kernel of 1x3, which are used to extract high-dimensional local time in the data. The temporal convolution of the traffic data segment X(1:k*T)=(X1........Xt......Xk*T) is: TC=Conv 1×t (X(1:k*T));

(2)图时间卷积层:将GCN与TCL堆叠为时空模块,但存在计算较多的问题,通过将这两个函数集成一个GTCL可很好的解决这一问题,替换之后的GTCL可表示为:

Figure GDA0004176717260000031
(2) Graph-time convolutional layer: GCN and TCL are stacked into a spatiotemporal module, but there is a problem of more calculations. This problem can be well solved by integrating these two functions into a GTCL. The GTCL after replacement can be expressed as:
Figure GDA0004176717260000031

基于行程时间预测的路网区域状态识别模型:采用基于层次划分复杂网络的社区划分算法中GN算法,社团发现是指复杂网络中各个节点的有向连接,并按划分评判将所有节点划分到社区中,相同社团中的节点存在类似的特征,GN算法效率较高,将所有节点及边的有向关系考虑到网络划分中,实验结果的划分结果具有层次性,路网分区结果的评判标准是模块度,模块度越大,则划分的社团结构也越明显,模块度Q的定义为:Road network regional state recognition model based on travel time prediction: The GN algorithm in the community partition algorithm based on hierarchical partitioning of complex networks is adopted. Community discovery refers to the directed connection of each node in a complex network, and all nodes are divided into communities according to the partition judgment. Nodes in the same community have similar characteristics. The GN algorithm is more efficient. The directed relationship of all nodes and edges is taken into account in the network partitioning. The partitioning results of the experimental results are hierarchical. The judgment standard of the road network partitioning results is modularity. The larger the modularity, the more obvious the community structure of the partition. The modularity Q is defined as:

Figure GDA0004176717260000041
Figure GDA0004176717260000041

式中:avw为节点v与w间边的权重,kw为节点w的度;Where: a vw is the weight of the edge between nodes v and w, k w is the degree of node w;

根据复杂网络中边属性,即GN算法中边值,将路网中路段的流量属性进行加权后,在进行最大边介数移除时,特别是在边介数相等时,根据当前时间段内各边之间流量比值得到对应边的权重,从而得到边权比,可得到路网中各边的相关性,得到的社团划分结果更优,与实际路网中两条边之间的相关性越可靠。According to the edge attributes in the complex network, that is, the edge value in the GN algorithm, the flow attributes of the road sections in the road network are weighted. When the maximum edge betweenness is removed, especially when the edge betweenness is equal, the weight of the corresponding edge is obtained according to the flow ratio between the edges in the current time period, thereby obtaining the edge weight ratio, and the correlation of the edges in the road network can be obtained. The obtained community division result is better and the correlation between two edges in the actual road network is more reliable.

基于Floyd算法的路径推荐模型:采用Floyd算法,解决任意两点之间的最短路径算法,解决有向图或是无向图的最短路径问题,算法的核心是通过局部最优求解全局最优,进行动态规划。首先寻找出目标OD的最短路径长度,然后记录下该长度的路径,即可寻找到推荐的路线,将路径中处于重度拥堵状态区域中的路段进行移除,再重复以上步骤,即可寻到基于拥堵区避让的推荐路径矩阵和路由矩阵;Path recommendation model based on Floyd algorithm: Floyd algorithm is used to solve the shortest path algorithm between any two points, and to solve the shortest path problem of directed or undirected graphs. The core of the algorithm is to solve the global optimum through local optimum and perform dynamic planning. First, find the shortest path length of the target OD, and then record the path of this length to find the recommended route. Remove the sections in the path that are in a severely congested state, and repeat the above steps to find the recommended path matrix and routing matrix based on congestion avoidance;

基于拥堵区避让的Floyd算法的思想是:The idea of the Floyd algorithm based on congestion avoidance is:

(1)初始状态下:根据DGCN预测所得的行程时间数据,写出图G初始基于距离的邻接矩阵W和初始路由矩阵R0=[ri j]n×n。对于每一对顶点vi和vj,如果v到v存在边,那么长度就是该边的权,如果没边就设长度为无穷大Q;(1) Initial state: Based on the travel time data predicted by DGCN, write the initial distance-based adjacency matrix W and the initial routing matrix R 0 = [ ri j ] n×n of the graph G. For each pair of vertices v i and v j , if there is an edge from v to v, then the length is the weight of the edge. If there is no edge, the length is set to infinity Q;

Figure GDA0004176717260000042
Figure GDA0004176717260000042

Figure GDA0004176717260000043
Figure GDA0004176717260000043

(2)k=0:即对于每一对顶点vi和vj,途径顶点的下标不大于k,实际上这里只能经过vo,该路径可分为两段,即(vi,vo)和(vo,vj),这一长度就是两段路径的长度之和,比较这一新路径和之前的路径(vi,vj),就可以确定vi到vj途经下标不大于k的最短路径;(2) k = 0: that is, for each pair of vertices vi and vj , the subscript of the passing vertex is not greater than k. In fact, only vo can be passed here. The path can be divided into two segments, namely ( vi , vo ) and ( vo , vj ). The length of this path is the sum of the lengths of the two paths. By comparing this new path with the previous path ( vi , vj ), we can determine the shortest path from vi to vj with a subscript not greater than k.

(3)k=1:同理,该路径可拆成(vi,……,vk)和(vk,..,vj)两段,这两段的长度在k=0的时候就确定了,再比较新路径和前面已知的路径(vi,vj)就可以确定途经下标不大于k的最短路径,此时k=1;(3) k = 1: Similarly, the path can be split into two segments (v i , …, v k ) and (v k , …, v j ). The lengths of these two segments are determined when k = 0. By comparing the new path with the previously known path (v i , v j ), the shortest path with an index not greater than k can be determined. At this time, k = 1;

(4)重复以上步骤,直到k=n-1为止,此时已经确定了从vi和vj所有可能的最短行程时间路径和对应的路由矩阵;(4) Repeat the above steps until k = n-1, at which point all possible shortest travel time paths from vi and vj and the corresponding routing matrices have been determined;

(5)选取的部分固定起终点的路径,去除路径中包含有重度拥堵路段的边,组成新的路网图G’,重复步骤(1)-(4),得到基于拥堵区避让的路径距离矩阵和路由矩阵,得出避让后所有的最短行程时间路径。(5) Select some fixed start-end paths, remove the edges containing heavily congested sections, form a new road network graph G’, repeat steps (1)-(4), obtain the path distance matrix and routing matrix based on congestion avoidance, and obtain all the shortest travel time paths after avoidance.

在上述的基于行程时间预测的拥堵区避让的路径推荐方法中,所述加权GN算法的具体流程为:In the above-mentioned congestion area avoidance route recommendation method based on travel time prediction, the specific process of the weighted GN algorithm is as follows:

第一步:初始状态下,每一个节点被看为一个独立的社区,即该状态下社区数与节点数相同;Step 1: In the initial state, each node is regarded as an independent community, that is, the number of communities in this state is the same as the number of nodes;

第二步:忽略数据中两个节点组成边的权重,按照无权网络计算方法求得每条边的边介数;Step 2: Ignore the weights of the edges between two nodes in the data, and calculate the edge betweenness of each edge using the unweighted network calculation method;

第三步:将边介数除以对应边的权重得到边权比;Step 3: Divide the edge betweenness by the weight of the corresponding edge to get the edge weight ratio;

第四步:移除边权比最大的边,计算当前网络的模块度Q;Step 4: Remove the edge with the largest edge-weight ratio and calculate the modularity Q of the current network;

第五步:对其余边重复第二步到第四步,并计算每一步的边权比和模块度,直至网络中所有边均被移除;Step 5: Repeat steps 2 to 4 for the remaining edges, and calculate the edge weight ratio and modularity of each step until all edges in the network are removed;

第六步:运行结束,取Q最大时对应的社团划分数量和分区结果。Step 6: After the operation is completed, the number of community divisions and the partition results corresponding to the maximum Q are obtained.

与现有技术相比,本发明基于行程时间预测的拥堵区避让的路径推荐方法的优点为:基于行程时间预测及拥堵区避让的路径推荐模型。基于历史卡口数据和路网拓扑结构,使用动态图神经网络(DGCN)对城市路网的路段行程时间进行短时预测,并通过结合路网动态实际状况和流量特征,建立基于加权GN算法的路网分区模型,根据预测的行程时间结果对路网进行区域状态识别;最后结合预测的行程时间和路网的拥堵状态,基于Floyd算法做避让的路径推荐。以某城市部分区域路网为例,将该区域的卡口数据和路网GIS-T数据应用于模型中。分析结果可得,对车辆进行避让拥堵区的路径推荐,从全局的路网进行考虑,一定程度上,减缓拥堵路段的拥堵程度,对出行者提供最短路径推荐,提高出行体验感Compared with the prior art, the advantages of the route recommendation method for avoiding congested areas based on travel time prediction of the present invention are: a route recommendation model based on travel time prediction and congested area avoidance. Based on historical checkpoint data and road network topology, a dynamic graph neural network (DGCN) is used to make short-term predictions on the travel time of sections of the urban road network, and by combining the actual dynamic conditions of the road network and traffic characteristics, a road network partitioning model based on the weighted GN algorithm is established, and the regional status of the road network is identified according to the predicted travel time results; finally, based on the predicted travel time and the congestion status of the road network, avoidance routes are recommended based on the Floyd algorithm. Taking the road network of a certain area in a certain city as an example, the checkpoint data and road network GIS-T data of the area are applied to the model. The analysis results show that the route recommendation for vehicles to avoid congested areas takes into account the global road network, which, to a certain extent, reduces the congestion level of congested sections, provides travelers with the shortest route recommendation, and improves the travel experience.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1是本发明基于行程时间预测及避让拥堵区的路径推荐模型示意图。FIG1 is a schematic diagram of a route recommendation model based on travel time prediction and congestion avoidance according to the present invention.

图2是本发明预测未来12个时段行程时间热力图。FIG. 2 is a heat map of travel time predicted for the next 12 time periods according to the present invention.

图3是本发明基于行程时间的路网邻接矩结果示意图。FIG3 is a schematic diagram of the road network adjacency moment result based on travel time according to the present invention.

图4是本发明第39轮路径邻接矩阵图(a),第39轮路由矩阵图(b)。FIG. 4 is a diagram of the 39th round path adjacency matrix (a) and a diagram of the 39th round routing matrix (b) of the present invention.

图5是本发明路网实际卡口点位示意图。FIG. 5 is a schematic diagram of actual checkpoint locations of the road network of the present invention.

具体实施方式DETAILED DESCRIPTION

以下是本发明的具体实施例并结合附图,对本发明的技术方案作进一步的描述,但本发明并不限于这些实施例。The following are specific embodiments of the present invention and the accompanying drawings to further describe the technical solution of the present invention, but the present invention is not limited to these embodiments.

如图1所示,基于DGCN的路段行程时间预测模型:As shown in Figure 1, the road section travel time prediction model based on DGCN:

更为准确的行程时间的预测是进行路径推荐的基础,而采样图神经网络对行程时间进行预测,特别是交通数据这种非结构化数据,可以很好的解决两个节点之间的相关性以及路网的复杂结构和动态性能的学习问题。因此,采用动态图卷积网络(DGCN)模型对行程时间进行预测。模型主要包含:拉普拉斯矩阵潜网络(LMLN)模块和基于图卷积网络(GCN)的交通预测模块。拉普拉斯矩阵预测单元(LMLN)包括三部分:(1)特征采样:对每天最近邻15min、30min、60min的数据进行采样,在减少特征的数据维度的同时对最近的交通进行采样;(2)空间注意力机制:为建立动态的交通路网的空间关系,采用注意力机制对路网当前近邻时间的邻接矩阵进行估计;(3)LSTM单元:采用LSTM来学习时间相关,探索邻接矩阵Ld的序列之间的内在关联。基于GCN的交通预测模块包含:(1)时间卷积层:主要包含四个卷积核为(1x3)的二维卷积层,用来提取数据中高维局部时间,交通数据段X(1:k*T)=(X1........Xt......Xk*T)的时间卷积为:TC=Conv1xt(X(1:k*T))。(2)图时间卷积层:将GCN与TCL堆叠为时空模块,但存在计算较多的问题,通过将这两个函数集成一个GTCL可很好的解决这一问题,替换之后的GTCL可表示为:

Figure GDA0004176717260000071
基于行程时间预测的路网区域状态识别模型More accurate prediction of travel time is the basis for route recommendation. Sampled graph neural network can predict travel time, especially for unstructured data such as traffic data, which can well solve the correlation between two nodes and the learning problem of complex structure and dynamic performance of road network. Therefore, the dynamic graph convolutional network (DGCN) model is used to predict travel time. The model mainly includes: Laplacian matrix latent network (LMLN) module and traffic prediction module based on graph convolutional network (GCN). The Laplacian matrix prediction unit (LMLN) consists of three parts: (1) Feature sampling: sampling the data of the nearest neighbor 15min, 30min, and 60min every day, sampling the nearest traffic while reducing the data dimension of the feature; (2) Spatial attention mechanism: in order to establish the spatial relationship of the dynamic traffic network, the attention mechanism is used to estimate the adjacency matrix of the current neighbor time of the road network; (3) LSTM unit: LSTM is used to learn time correlation and explore the intrinsic correlation between the sequences of the adjacency matrix Ld. The traffic prediction module based on GCN includes: (1) Temporal convolution layer: It mainly includes four two-dimensional convolution layers with convolution kernels of (1x3), which are used to extract high-dimensional local time in the data. The temporal convolution of the traffic data segment X(1:k*T)=( X1 ........ Xt. .....Xk *T ) is: TC=Conv1xt(X(1:k*T)). (2) Graph temporal convolution layer: GCN and TCL are stacked into a spatiotemporal module, but there is a problem of more calculations. This problem can be well solved by integrating these two functions into a GTCL. The GTCL after replacement can be expressed as:
Figure GDA0004176717260000071
Road network area status identification model based on travel time prediction

交通网络具有复杂网络的结构和特征,对路网进行划分的模型包含网格划分区、依照路网分区、按照研究对象密度分区、依照参考点分区等,本研究采用基于层次划分复杂网络的社区划分算法中GN算法。社团发现是指复杂网络中各个节点的有向连接,并按划分评判将所有节点划分到社区中,相同社团中的节点存在类似的特征,GN算法效率较高,将所有节点及边的有向关系考虑到网络划分中,实验结果的划分结果具有层次性。路网分区结果的评判标准是模块度,模块度越大,则划分的社团结构也越明显。模块度的定义为:The traffic network has the structure and characteristics of a complex network. The models for dividing the road network include grid division areas, division according to the road network, division according to the density of the research object, and division according to the reference point. This study adopts the GN algorithm in the community division algorithm based on hierarchical division of complex networks. Community discovery refers to the directed connection of each node in a complex network, and all nodes are divided into communities according to the division judgment. Nodes in the same community have similar characteristics. The GN algorithm is more efficient and takes the directed relationships of all nodes and edges into account in the network division. The division results of the experimental results are hierarchical. The criterion for judging the results of road network partitioning is modularity. The larger the modularity, the more obvious the structure of the divided community. The definition of modularity is:

Figure GDA0004176717260000072
Figure GDA0004176717260000072

式中:avw为节点v与w间边的权重,kw为节点w的度;Where: a vw is the weight of the edge between nodes v and w, k w is the degree of node w;

根据复杂网络中边属性,即GN算法中边值,将路网中路段的流量属性进行加权后,在进行最大边介数移除时,特别是在边介数相等时,根据当前时间段内各边之间流量比值得到对应边的权重,从而得到边权比,可得到路网中各边的相关性,得到的社团划分结果更优,与实际路网中两条边之间的相关性越可靠。According to the edge attributes in the complex network, that is, the edge value in the GN algorithm, the flow attributes of the road sections in the road network are weighted. When the maximum edge betweenness is removed, especially when the edge betweenness is equal, the weight of the corresponding edge is obtained according to the flow ratio between the edges in the current time period, thereby obtaining the edge weight ratio, and the correlation of the edges in the road network can be obtained. The obtained community division result is better and the correlation between two edges in the actual road network is more reliable.

加权GN算法的具体流程为:①初始状态下,每一个节点被看为一个独立的社区,即该状态下社区数与节点数相同;②忽略数据中两个节点组成边的权重,按照无权网络计算方法求得每条边的边介数;③将边介数除以对应边的权重得到边权比;④移除边权比最大的边,计算当前网络的模块度Q;⑤对其余边重复步骤②到步骤④,并计算每一步的边权比和模块度,直至网络中所有边均被移除;⑥运行结束,取Q最大时对应的社团划分数量和分区结果。The specific process of the weighted GN algorithm is as follows: ① In the initial state, each node is regarded as an independent community, that is, the number of communities is the same as the number of nodes in this state; ② Ignore the weights of the edges composed of two nodes in the data, and calculate the edge betweenness of each edge according to the unweighted network calculation method; ③ Divide the edge betweenness by the weight of the corresponding edge to obtain the edge weight ratio; ④ Remove the edge with the largest edge weight ratio and calculate the modularity Q of the current network; ⑤ Repeat steps ② to ④ for the remaining edges, and calculate the edge weight ratio and modularity of each step until all edges in the network are removed; ⑥ At the end of the operation, take the number of community divisions and partition results corresponding to the maximum Q.

区域平均速度求解方法:Method for calculating regional average velocity:

(1)根据该评价方法求某路段的平均行程速度Vi(1) According to the evaluation method, the average travel speed V i of a certain road section is calculated as follows:

Figure GDA0004176717260000081
Figure GDA0004176717260000081

式中:L为路段长度单位为km;t车辆i通过区间路段的时间单位为h;n为测定车辆数。Where: L is the length of the road section in km; t is the time it takes for vehicle i to pass through the section in h; n is the number of measured vehicles.

(2)计算路网或某区域内道路平均行程速度:(2) Calculate the average travel speed of a road network or a certain area:

Figure GDA0004176717260000082
Figure GDA0004176717260000082

式中:V1........Vn表示路网道路平均行程速度,单位为km/h,n为路段数量。Where: V 1 ........V n represents the average travel speed of the road network, in km/h, and n is the number of road sections.

基于Floyd算法的路径推荐模型Route recommendation model based on Floyd algorithm

根据需求可以将最优路径问题可分为距离最短和权重最小,其中,基于距离最短方法中可分为自由路径和限制路径,本研究采用Floyd(弗洛伊德)算法,是解决任意两点之间的最短路径算法,可解决有向图或是无向图的最短路径问题,该算法的核心是通过局部最优求解全局最优,进行动态规划。首先寻找出目标OD的最短路径长度,然后记录下该长度的路径,即可寻找到推荐的路线,将路径中处于重度拥堵状态区域中的路段进行移除,再重复以上步骤,即可寻到基于拥堵区避让的推荐路径矩阵和路由矩阵。According to the needs, the optimal path problem can be divided into the shortest distance and the minimum weight. Among them, the shortest distance method can be divided into free path and restricted path. This study uses the Floyd algorithm, which is an algorithm to solve the shortest path between any two points. It can solve the shortest path problem of directed or undirected graphs. The core of this algorithm is to solve the global optimum through local optimum and perform dynamic planning. First, find the shortest path length of the target OD, and then record the path of this length to find the recommended route, remove the sections in the path that are in a severely congested state, and repeat the above steps to find the recommended path matrix and routing matrix based on congestion avoidance.

基于拥堵区避让的Floyd算法的思想是:The idea of the Floyd algorithm based on congestion avoidance is:

(1)初始状态下:根据DGCN预测所得的行程时间数据,写出图G初始基于距离的邻接矩阵W和初始路由矩阵R0=[ri j]n×n,对于每一对顶点vi和vj,如果v到v存在边,那么长度就是该边的权,如果没边就设长度为无穷大Q。(1) Initial state: Based on the travel time data predicted by DGCN, write the initial distance-based adjacency matrix W and the initial routing matrix R 0 = [ ri j ] n×n of the graph G. For each pair of vertices vi and v j , if there is an edge from v to v, then the length is the weight of the edge. If there is no edge, the length is set to infinity Q.

Figure GDA0004176717260000091
Figure GDA0004176717260000091

Figure GDA0004176717260000092
Figure GDA0004176717260000092

(2)k=0:即对于每一对顶点vi和vj,途径顶点的下标不大于k,实际上这里只能经过vo,该路径可分为两段,即(vi,vo)和(vo,vj),这一长度就是两段路径的长度之和,比较这一新路径和之前的路径(vi,vj),就可以确定vi到vj途经下标不大于k的最短路径。(2) k = 0: that is, for each pair of vertices vi and vj , the subscript of the vertex passed through is not greater than k. In fact, it can only pass through vo . The path can be divided into two segments, namely ( vi , vo ) and ( vo , vj ). The length of this path is the sum of the lengths of the two paths. By comparing this new path with the previous path ( vi , vj ), we can determine the shortest path from vi to vj with a subscript not greater than k.

(3)k=1:同理,该路径可拆成(vi,……,vk)和(vk,..,vj)两段,这两段的长度在k=0的时候就确定了,再比较新路径和前面已知的路径(vi,vj)就可以确定途经下标不大于k的最短路径,此时k=1。(3) k = 1: Similarly, the path can be divided into two segments: (v i ,…, v k ) and (v k ,…, v j ). The lengths of these two segments are determined when k = 0. By comparing the new path with the previously known path (v i , v j ), the shortest path with an index no greater than k can be determined. At this time, k = 1.

(4)重复以上步骤,直到k=n-1为止,此时已经确定了从vi和vj所有可能的最短行程时间路径和对应的路由矩阵。(4) Repeat the above steps until k = n-1, at which point all possible shortest travel time paths from vi and vj and the corresponding routing matrices have been determined.

(5)选取的部分固定起终点的路径,去除路径中包含有重度拥堵路段的边,组成新的路网图G’,重复步骤(1)-(4),得到基于拥堵区避让的路径距离矩阵和路由矩阵,得出避让后所有的最短行程时间路径。(5) Select some fixed start-end paths, remove the edges containing heavily congested sections, form a new road network graph G’, repeat steps (1)-(4), obtain the path distance matrix and routing matrix based on congestion avoidance, and obtain all the shortest travel time paths after avoidance.

实施例一Embodiment 1

基于DGCN的行程时间预测Travel time prediction based on DGCN

3.1.1数据集介绍3.1.1 Dataset Introduction

所使用的数据集来源于中心城区公开的卡口数据集,openITS平台整理发布,该区域内的总道路长度为174km,包含2017年12月3日至9日全天的卡口过车数据和卡口GIS-T布点数据,卡口过车数据中一天的过车量数据大约有83万条,删除未识别的车牌号后数据量约74万,根据求得的路段行程时间经过筛选后数据量约56万,具体字段如下表所示。日全天行程时间数据,包含一个交通特征即路段行程时间,在模型中将预测行程时间作为输出;时隙划分15分钟,每天的数据包含96个时隙,每小时所有数据有4个样本,因此设置T=4,选取数据集的60%作为训练集、20%作为验证集、20%作为测试集。同时,选择验证集上最佳模型参数,并在测试集上评估所有的模型,每个路段的数据样本用归一化The dataset used is from the public checkpoint dataset in the central urban area, which is organized and released by the openITS platform. The total road length in the area is 174km, including the checkpoint vehicle passing data and checkpoint GIS-T distribution data from December 3 to 9, 2017. There are about 830,000 vehicle passing data in one day in the checkpoint vehicle passing data. After deleting the unrecognized license plate numbers, the data volume is about 740,000. After screening according to the obtained section travel time, the data volume is about 560,000. The specific fields are shown in the following table. The daily travel time data contains a traffic feature, namely the section travel time. The predicted travel time is used as the output in the model; the time slot is divided into 15 minutes, and the data of each day contains 96 time slots. There are 4 samples of all data per hour, so T=4 is set, and 60% of the data set is selected as the training set, 20% as the validation set, and 20% as the test set. At the same time, the best model parameters on the validation set are selected, and all models are evaluated on the test set. The data samples of each section are normalized.

Figure GDA0004176717260000101
Figure GDA0004176717260000101

参数设置Parameter settings

实验运用于Linux操作系统ubuntu18.04,GPU为RTX3090*2,显存为48G,采用python3.7代码编写,模型架构采用基于pytorch1.2.0深度学习根据进行开发。The experiment was applied to the Linux operating system ubuntu18.04, the GPU was RTX3090*2, the video memory was 48G, the code was written in python3.7, and the model architecture was developed based on pytorch1.2.0 deep learning.

在图卷积和时间卷积中使用4个1×3的卷积核,预测时间步长c为12,学习率设为0.0005,每一轮的衰减率为0.92,批次大小为8,使用l2_loss为损失函数,利用Adam作为优化器进行优化,在GTCL中设置切比雪夫多项式M=3,卷积核ts=3,多头注意力机制k=4,预测时间间隔Tp=12,即未来一小时,即我们采用60%的样本进行训练和预测12个样本。Four 1×3 convolution kernels are used in graph convolution and temporal convolution, the prediction time step c is 12, the learning rate is set to 0.0005, the decay rate of each round is 0.92, the batch size is 8, l2_loss is used as the loss function, and Adam is used as the optimizer for optimization. In GTCL, the Chebyshev polynomial M=3, the convolution kernel ts=3, the multi-head attention mechanism k=4, and the prediction time interval Tp=12, that is, one hour in the future, that is, we use 60% of the samples for training and predicting 12 samples.

行程时间预测结果Travel time prediction results

为了估计不同输入的影响,我们还使用最近邻的数据建立了我们的输入,用DGCN_R表示。为了进一步评估不同拉普拉斯矩阵对GCN的效率,特别是GAT,我们将我们的方法与四种基于GCN的方法进行了比较(1)ASTGCN,其中使用了注意拉普拉斯矩阵;(2)DGCN_Mask,使用了Mask拉普拉斯矩阵;(3)GCN_Res,使用剩余拉普拉斯矩阵;(4)DGCN_GAT,一种用GAT代替模型的空间特征层GTCL的方法,因此与其他方法具有可比性。所有方法的性能都由三个指标来衡量,平均绝对百分比误差(MAPE)、平均绝对误差(MAE)和均方误差(RMSE),公式如下:To estimate the impact of different inputs, we also built our input using the nearest neighbor data, denoted by DGCN_R. To further evaluate the efficiency of different Laplacian matrices for GCN, especially GAT, we compared our method with four GCN-based methods: (1) ASTGCN, which uses the attention Laplacian matrix; (2) DGCN_Mask, which uses the Mask Laplacian matrix; (3) GCN_Res, which uses the residual Laplacian matrix; (4) DGCN_GAT, a method that replaces the spatial feature layer GTCL of the model with GAT, thus making it comparable to other methods. The performance of all methods is measured by three metrics, mean absolute percentage error (MAPE), mean absolute error (MAE), and mean square error (RMSE), as follows:

Figure GDA0004176717260000102
Figure GDA0004176717260000102

Figure GDA0004176717260000103
Figure GDA0004176717260000103

Figure GDA0004176717260000104
Figure GDA0004176717260000104

在数据集上平均一小时交通预测精度如下表所示,结果表明,DGCN_Res模型在所有指标是具有最佳性能。拉普拉斯潜在矩阵可以更好的提取道路网络的动态空间关,同时,DGCN_Res的精度优于其他模型,因为用经验拉普拉斯矩阵代替了全局优化的残差拉普拉斯矩阵。预测未来12个时段的结果,预测精度结果The average one-hour traffic prediction accuracy on the dataset is shown in the following table. The results show that the DGCN_Res model has the best performance in all indicators. The Laplace latent matrix can better extract the dynamic spatial relationship of the road network. At the same time, the accuracy of DGCN_Res is better than other models because the empirical Laplace matrix is used instead of the globally optimized residual Laplace matrix. The results of predicting the next 12 periods, the prediction accuracy results

Figure GDA0004176717260000121
Figure GDA0004176717260000121

加权GN算法的路网分区和区域状态识别结果Results of road network partitioning and regional status recognition using the weighted GN algorithm

实验的数据集来源于卡口GIS-T布点数据,其中包括39个卡口点,结合行程时间数据选取60个路段作为边,采用传统GN算法和加权GN算法对城市路网进行划分,不同分区数量和对应模块度如下表所示,传统GN算法和加权GN算法,随着分区数量的增加,模块度都在逐渐降低,说明实际路网中,各边之间的相关性也在降低,结合拥堵区避让的条件和实际路网中路段相关度,选择分区数量为15,分区结果为:[['11','1','2','12','10','26'],['35','31','3'],['7','4','13','14','9'],['5','36','34'],['39','24','8','6'],['30','16','15'],['18','17'],['21','28','29','19'],['20'],['22'],['27','23'],['25','33'],['32'],['37'],['38']],并根据分区结果,按顺序将区域划分为G1—G15,以便表示下一节拥堵区识别结果。The experimental data set comes from the checkpoint GIS-T distribution data, which includes 39 checkpoints. Combined with the travel time data, 60 road sections are selected as edges. The traditional GN algorithm and the weighted GN algorithm are used to divide the urban road network. The number of partitions and the corresponding modularity are shown in the following table. For the traditional GN algorithm and the weighted GN algorithm, as the number of partitions increases, the modularity gradually decreases, indicating that in the actual road network, the correlation between the edges is also decreasing. Combined with the conditions for avoiding congested areas and the correlation of road sections in the actual road network, the number of partitions is selected as 15, and the partition results are: [['11','1','2','12','10','26'],['35 ','31','3'],['7','4','13','14','9'],['5','36','34'],['39','24','8','6'],['30','16','15'],['18','17'],['21','28','29','19'],['20'],['22'],['27','23'],['25','33'],['32'],['37'],['38']], and according to the partition results, the area is divided into G1-G15 in sequence to represent the congestion area identification results in the next section.

Figure GDA0004176717260000131
Figure GDA0004176717260000131

根据上一节行程时间预测的结果,可得预测的60个路段的行程时间数据,结合路网分区结果和区域状态识别的步骤,得部分区域交通状态如下表所示:According to the results of travel time prediction in the previous section, the travel time data of 60 road sections can be obtained. Combined with the results of road network partitioning and the steps of regional status identification, the traffic status of some regions is shown in the following table:

Figure GDA0004176717260000132
Figure GDA0004176717260000132

基于Floyd算法的避让拥堵区的路径推荐结果Route recommendation results for avoiding congestion areas based on Floyd algorithm

初始状态下基于行程时间的路网邻接矩阵如下图3所示,两点之间的权重代表基于DGCN预测未来15分钟后该路段的行程时间,记为v(i,j)的值,将v(i,j)记录为∞,同时规定其的权值为0。依次遍历v(1,j)和v(i,1),比较v(1,j)、v(i,1)和v(i,j)的大小,如果v(1,j)+v(i,1)<v(i,j)则进行替换,遍历矩阵39次,记录每一次的矩阵结果,最后得到的最短行程时间的邻接矩阵和对应的路由矩阵如下图4(a)和(b)和图5所示。The road network adjacency matrix based on travel time in the initial state is shown in Figure 3 below. The weight between two points represents the travel time of the road section predicted by DGCN in the next 15 minutes, which is recorded as the value of v(i,j). v(i,j) is recorded as ∞, and its weight is set to 0. Traverse v(1,j) and v(i,1) in turn, compare the sizes of v(1,j), v(i,1) and v(i,j), and replace if v(1,j)+v(i,1)<v(i,j). Traverse the matrix 39 times, record the matrix results each time, and finally obtain the adjacency matrix of the shortest travel time and the corresponding routing matrix as shown in Figures 4(a), (b) and 5.

最短路径不一定最优,仅考虑路径最短进行路径推荐,对于全局路网架构,会造成被推荐的路径产生拥堵,若对部分车辆进行拥堵区避让的路径推荐,一定程度上,可减缓拥堵区的交通压力。以安徽省宣城市路网为例,将流量这一交通特征作为权重,并结合边介数组成边权比,通过加权GN算法结合实际路网情况进行分区;将历史卡口数据应用到模型中,对每天最近邻15min,30min,45min的采样,引入最近邻数据的拉普拉斯潜在矩阵,更好的学习路网的动态结构,对于行程时间的预测比在Pems04的MAE和RMSE指标上高出18.7%和30%;同时,求得各卡口点之间的最短路径和通过对拥堵区进行避让后的最短路径结果进行对比,得出避让拥堵区的路径推荐在出行所需行程时间上更少。本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。The shortest path is not necessarily the best. Only considering the shortest path for path recommendation will cause congestion in the recommended path for the global road network architecture. If the path recommendation for avoiding congested areas is performed for some vehicles, the traffic pressure in the congested areas can be alleviated to a certain extent. Taking the road network of Xuancheng City, Anhui Province as an example, the traffic feature of flow is used as the weight, and the edge weight ratio is formed by combining the edge betweenness, and the weighted GN algorithm is used to divide the area according to the actual road network situation; the historical checkpoint data is applied to the model, and the Laplace potential matrix of the nearest neighbor data is introduced for sampling the nearest neighbor 15min, 30min, and 45min every day, so as to better learn the dynamic structure of the road network, and the prediction of travel time is 18.7% and 30% higher than the MAE and RMSE indicators of Pems04; at the same time, the shortest path between each checkpoint is obtained and the shortest path result after avoiding the congested area is compared, and it is concluded that the path recommendation for avoiding the congested area is less in the travel time required for travel. The content not described in detail in this specification belongs to the prior art known to professional and technical personnel in this field. The specific embodiments described herein are merely examples of the spirit of the present invention. Those skilled in the art may make various modifications or additions to the specific embodiments described or replace them in similar ways, but they will not deviate from the spirit of the present invention or exceed the scope defined by the appended claims.

Claims (2)

1. The route recommendation method for avoiding the congestion area based on travel time prediction is characterized by comprising three steps of a road section travel time prediction model based on DGCN, a road network area state identification model based on travel time prediction and a route recommendation model based on Floyd algorithm, and specifically comprising the following steps of:
the road section travel time prediction model based on the DGCN comprises a Laplacian matrix submarine network module and a traffic prediction module based on a graph rolling network, wherein: the laplacian matrix prediction unit includes three parts:
(1) Feature sampling: sampling data of 15min,30min and 60min nearest to each other every day, and sampling the nearest traffic while reducing the data dimension of the features;
(2) Spatial attention mechanism: in order to establish a space relation of a dynamic traffic network, estimating an adjacent matrix of the current neighbor time of the road network by adopting an attention mechanism;
(3) LSTM unit: learning a time correlation by using LSTM, and exploring an internal correlation between sequences of the adjacency matrix Ld;
the GCN-based traffic prediction module includes two parts:
(1) Time convolution layer: comprising four two-dimensional convolution layers with a convolution kernel of 1X3, for extracting high-dimensional local time in the data, traffic data segment X (1:kχt) = (X1.. Xt...a. Xk X T) is:
TC=Conv 1×t (X(1:k*T))(1)
TC=Conv 1×t (X (1: kT)) is a time convolution
(2) Graph time convolution layer: stacking GCN and TCL as space-time modules, but there is a problem of large calculation amount, this problem can be solved well by integrating these two functions into one GTCL, and the GTCL after replacement can be expressed as:
Figure FDA0004176717250000011
road network area state identification model based on travel time prediction: the community discovery refers to directed connection of all nodes in the complex network by adopting a GN algorithm in a community division algorithm based on a hierarchical division complex network, all nodes in the same community are divided into communities according to division judgment, the nodes in the same community have similar characteristics, the GN algorithm is high in efficiency, the directed relationship between all nodes and edges is considered in network division, the division result of an experimental result has hierarchy, the judgment standard of the road network division result is modularity, the greater the modularity is, the more obvious the divided community structure is, and the definition of the modularity Q is as follows:
Figure FDA0004176717250000021
wherein: a, a vw Is the weight of the edge between the nodes v and w, k w The degree of node w;
according to the edge attribute in the complex network, namely the edge value in the GN algorithm, after weighting the flow attribute of the road section in the road network, when the maximum edge betweenness is removed, when the edge betweenness is equal, the weight of the corresponding edge is obtained according to the flow ratio between the edges in the current time period, so that the edge weight ratio is obtained, the relevance of each edge in the road network can be obtained, the obtained community division result is better, and the relevance between the community division result and the two edges in the actual road network is more reliable;
path recommendation model based on Floyd algorithm: the Floyd algorithm is also called an insertion point method, the shortest path algorithm between any two points can be solved, the problem of the shortest path of a directed graph or an undirected graph is solved, the core of the algorithm is that global optimization is solved through local optimization, dynamic planning is carried out, the shortest path length of a target OD is found out firstly, then the path of the length is recorded, a recommended route can be found, the road sections in the severe congestion state area in the path are removed, the algorithm is repeated, and a recommended path matrix and a route matrix based on congestion area avoidance can be found;
the idea of the Floyd algorithm based on congestion area avoidance is that:
(1) In the initial state, according to the travel time data predicted by DGCN, writing out an adjacent matrix W and an initial route matrix R based on the distance of the initial graph G 0 =[r ij ] n×n For each pair of vertices v i And v j If v i And v j If there is an edge, the edge is weighted asRoad segment length; if no edge exists, setting the length to be infinite;
Figure FDA0004176717250000022
Figure FDA0004176717250000023
(2) k=0, i.e. for each pair of vertices v i And v j The subscript of the path vertex is not greater than k, where in practice only v can pass o The path may be divided into two segments, i.e. (v) i ,v o ) Sum (v) o ,v j ) This length is the sum of the lengths of the two paths, and the new path is compared with the previous path (v i ,v j ) V can be determined i To v j A shortest path with a path index not greater than k;
(3) k=1, and the path is detachable (v) i ,……,v k ) Sum (v) k ,..,v j ) Two sections whose length is determined when k=0, and comparing the new path with the previously known path (v i ,v j ) The shortest path with a path index no greater than k can be determined, where k=1;
(4) Repeating the above steps until k=n-1, at which time the slave v has been determined i And v j All possible shortest travel time paths and corresponding routing matrices;
(5) Selecting a part of paths with fixed starting and ending points, removing edges of the paths containing the severely congested road sections to form a new road network graph G', and repeating the steps (1) - (4) to obtain a path distance matrix and a route matrix based on congestion area avoidance, so as to obtain all the shortest travel time paths after the avoidance.
2. The route recommendation method for congestion area avoidance based on travel time prediction according to claim 1, wherein the specific flow of the weighted GN algorithm is as follows:
the first step: in the initial state, each node is regarded as an independent community, namely the number of communities is the same as the number of nodes in the state;
and a second step of: neglecting the weight of the edges formed by two nodes in the data, and obtaining the edge betweenness of each edge according to an unauthorized network calculation method;
and a third step of: dividing the edge betweenness by the weight of the corresponding edge to obtain an edge weight ratio;
fourth step: removing the edge with the largest edge weight ratio, and calculating the modularity Q of the current network;
fifth step: repeating the second step to the fourth step for the rest edges, and calculating the edge weight ratio and the modularity of each step until all edges in the network are removed;
sixth step: and (5) after the operation is finished, taking the corresponding community division number and partition result when Q is maximum.
CN202211021689.7A 2022-08-24 2022-08-24 Route recommendation method for avoiding congestion area based on travel time prediction Active CN115409256B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211021689.7A CN115409256B (en) 2022-08-24 2022-08-24 Route recommendation method for avoiding congestion area based on travel time prediction

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211021689.7A CN115409256B (en) 2022-08-24 2022-08-24 Route recommendation method for avoiding congestion area based on travel time prediction

Publications (2)

Publication Number Publication Date
CN115409256A CN115409256A (en) 2022-11-29
CN115409256B true CN115409256B (en) 2023-05-12

Family

ID=84161268

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211021689.7A Active CN115409256B (en) 2022-08-24 2022-08-24 Route recommendation method for avoiding congestion area based on travel time prediction

Country Status (1)

Country Link
CN (1) CN115409256B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115952364B (en) * 2023-03-07 2023-05-23 之江实验室 Route recommendation method and device, storage medium and electronic equipment
CN116403410B (en) * 2023-06-06 2023-08-22 中南大学 Highway mixed path induction model construction method considering congestion traffic sources
CN117273590B (en) * 2023-10-19 2024-07-09 苏州大学 Neural combination optimization method and system for solving vehicle path optimization problem
CN117669982B (en) * 2023-12-18 2024-05-17 山东建筑大学 Park comprehensive energy system capacity configuration method considering user multi-comfort
CN118333525A (en) * 2024-05-10 2024-07-12 杭州安那其科技有限公司 Intelligent logistics warehouse management system
CN119129298A (en) * 2024-11-14 2024-12-13 中国电子科技集团公司第二十八研究所 A virtual-real combined unmanned autonomous algorithm verification method and system

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106918345A (en) * 2017-03-27 2017-07-04 中国农业大学 A kind of optimization method and device in scenic region guide path
CN107944605A (en) * 2017-11-10 2018-04-20 河海大学常州校区 A kind of dynamic traffic paths planning method based on data prediction
CN109215343A (en) * 2018-09-20 2019-01-15 山东交通学院 Road network congestion regions dynamic identifying method and system based on community discovery

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108562301A (en) * 2018-05-21 2018-09-21 北京石油化工学院 A kind of method and device for planning of driving path
CN112444256A (en) * 2019-08-27 2021-03-05 扬州盛世云信息科技有限公司 Method for time shortest path based on road traffic flow
CN112785085B (en) * 2021-02-08 2022-11-29 日日顺供应链科技股份有限公司 Distribution path optimization method and device

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106918345A (en) * 2017-03-27 2017-07-04 中国农业大学 A kind of optimization method and device in scenic region guide path
CN107944605A (en) * 2017-11-10 2018-04-20 河海大学常州校区 A kind of dynamic traffic paths planning method based on data prediction
CN109215343A (en) * 2018-09-20 2019-01-15 山东交通学院 Road network congestion regions dynamic identifying method and system based on community discovery

Also Published As

Publication number Publication date
CN115409256A (en) 2022-11-29

Similar Documents

Publication Publication Date Title
CN115409256B (en) Route recommendation method for avoiding congestion area based on travel time prediction
Jindal et al. Optimizing taxi carpool policies via reinforcement learning and spatio-temporal mining
Sadeghi-Niaraki et al. Short-term traffic flow prediction using the modified elman recurrent neural network optimized through a genetic algorithm
CN106251625B (en) Three-dimensional urban road network global state prediction technique under big data environment
CN114519932B (en) An Integrated Forecasting Method of Regional Traffic Conditions Based on Spatiotemporal Relationship Extraction
CN114202120B (en) Urban traffic journey time prediction method for multi-source heterogeneous data
CN110836675B (en) Decision tree-based automatic driving search decision method
CN112017436B (en) Method and system for predicting urban traffic travel time
CN114881356A (en) Prediction method of urban traffic carbon emission based on particle swarm optimization optimization of BP neural network
Ziat et al. Joint prediction of road-traffic and parking occupancy over a city with representation learning
CN117350519A (en) Charging station planning method and system based on new energy passenger vehicle charging demand prediction
Russo et al. A topological method to choose optimal solutions after solving the multi-criteria urban road network design problem
Menguc et al. A data driven approach to forecasting traffic speed classes using extreme gradient boosting algorithm and graph theory
Li et al. Ridesplitting demand prediction via spatiotemporal multi-graph convolutional network
CN110674990B (en) Instant distribution path selection method and system with sliding window updating mechanism
Xing et al. Recommendation of urban vehicle driving routes under traffic congestion: A traffic congestion regulation method considering road network equilibrium
CN112183871A (en) Urban Traffic Guidance System Based on Air Index
Liu et al. Learning to route via theory-guided residual network
CN111325414A (en) A Data-driven Parallel Traffic Prediction Method
Shi et al. An adaptive route guidance model considering the effect of traffic signals based on deep reinforcement learning
Ajani et al. Dynamic path planning approaches based on artificial intelligence and machine learning
Yao et al. A path planning model based on spatio-temporal state vector from vehicles trajectories
Li et al. A Vehicle Path Planning and Prediction Algorithm Based on Attention Mechanism for Complex Traffic Intersection Collaboration in Intelligent Transportation
CN113343358A (en) Electric vehicle charging load space-time distribution modeling method considering road information
Patanè et al. An Explainable Model Framework for Vehicle-to-Grid Available Aggregated Capacity Prediction

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