CN114169647B - 持续记忆的自适应异质时空图卷积的交通预测方法及系统 - Google Patents
持续记忆的自适应异质时空图卷积的交通预测方法及系统 Download PDFInfo
- Publication number
- CN114169647B CN114169647B CN202210014979.2A CN202210014979A CN114169647B CN 114169647 B CN114169647 B CN 114169647B CN 202210014979 A CN202210014979 A CN 202210014979A CN 114169647 B CN114169647 B CN 114169647B
- Authority
- CN
- China
- Prior art keywords
- space
- time
- heterogeneous
- diagram
- layer
- 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
- 238000010586 diagram Methods 0.000 title claims abstract description 171
- 230000015654 memory Effects 0.000 title claims abstract description 81
- 238000000034 method Methods 0.000 title claims abstract description 37
- 239000011159 matrix material Substances 0.000 claims description 87
- 238000010276 construction Methods 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 15
- 230000008859 change Effects 0.000 claims description 14
- 238000011156 evaluation Methods 0.000 claims description 14
- 230000003044 adaptive effect Effects 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 7
- 230000004927 fusion Effects 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 6
- 238000012549 training Methods 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims description 3
- 238000010008 shearing Methods 0.000 claims description 3
- 238000012360 testing method Methods 0.000 claims description 3
- 238000010200 validation analysis Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 10
- 230000007774 longterm Effects 0.000 abstract description 5
- 238000013528 artificial neural network Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 101001095088 Homo sapiens Melanoma antigen preferentially expressed in tumors Proteins 0.000 description 1
- 102100037020 Melanoma antigen preferentially expressed in tumors Human genes 0.000 description 1
- 206010039203 Road traffic accident Diseases 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013527 convolutional neural network Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000007787 long-term memory Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- General Business, Economics & Management (AREA)
- Software Systems (AREA)
- Marketing (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Development Economics (AREA)
- Data Mining & Analysis (AREA)
- Primary Health Care (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Educational Administration (AREA)
- Operations Research (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于交通预测技术领域,具体公开了一种持续记忆的自适应异质时空图卷积的交通预测方法及系统,该方法将交通流量数据和流量数据的历史记忆,并输入记忆输入层,记忆输入层输出时间序列,并将该时间序列作为异质时空图卷积层的第一个子层的输入,异质时空图卷积层设有多个子层,前一个子层的输出是下一个子层的输入,构建不同时空异质图,并利用时空异质图完成图卷积操作,异质时空图卷积层的各层输出时间序列至时空信息融合层,得到交通流量预测数据和新的历史记忆。采用本技术方案,捕获交通流量数据的异质性,通过历史信息获取交通流量的长期时间依赖,提升预测效果。
Description
技术领域
本发明属于交通预测技术领域,涉及一种持续记忆的自适应异质时空图卷积的交通预测方法及系统。
背景技术
随着智能交通系统(I TS)的发展,交通流量数据的预测成为了I TS不可或缺的一部分,准确及时的交通预测可以帮助进行有效的交通管制。交通流量是道路状况的基础指标,如果可以有效预测交通流量数据,I TS能够更有效、更合理地规划车辆路线和交通信号灯。
交通流量数据有以下三个特点:(1)时间相关性(时间依赖),一个区域的流量与该区域之前一段时间的流量相关;(2)空间相关性(空间依赖),车辆从一个区域驶向另一个区域,导致同一时刻下,一个区域的流量会受到另一个区域流量的影响;(3)异质性,不同时间下,一个区域的流量受到其他区域流量的影响是不一样的。所以,如果可以建模交通流量数据的特点,将达到更好的预测效果。
为了检测道路的交通状况,人们将大量的传感器应用到交通路网中,并借此积累了大量关于交通流量的历史数据。许多研究人员借助历史数据挖掘交通流量数据的时空相关性,来达到预测交通流量、检测交通事故发生等研究目的。预测交通流量的研究方法在十几年间发生了明显的变化,从早期基于统计学的方法,到将机器学习的方法以及深度学习的方法应用到交通流量预测中,模型的预测性能有了显著的改善。
但现有技术中采用的时空模型在捕获时间依赖时,往往只建立同一种图结构,忽视不同时间下,一个区域的流量受到其他区域流量的影响是不一样的现象,即交通流量数据的异质性。且以往的模型采用LSTM或是GRU等模块,来获取交通流量数据的时间依赖,这样获取的时间依赖的时间跨度局限在当前输入数据的时间长度,不能记忆之前输入数据的历史信息。
发明内容
本发明的目的在于提供一种持续记忆的自适应异质时空图卷积的交通预测方法及系统,捕获交通流量数据的异质性,通过历史信息获取交通流量的长期时间依赖,提升预测效果。
为了达到上述目的,本发明的基础方案为:一种持续记忆的自适应异质时空图卷积的交通预测方法,包括如下步骤:
获取交通流量数据和流量数据的历史记忆,并输入记忆输入层;
记忆输入层输出时间序列,并将该时间序列作为异质时空图卷积层的第一个子层的输入;
异质时空图卷积层设有多个子层,前一个子层的输出是下一个子层的输入,构建不同时空异质图,并利用时空异质图完成图卷积操作;
异质时空图卷积层的各层输出时间序列至时空信息融合层,得到交通流量预测数据和新的历史记忆。
本基础方案的工作原理和有益效果在于:本方案构建不同时刻的时空图的图结构,以此捕获交通流量数据的异质性。将当前输入序列的信息保存到历史信息中,到达持续记忆的效果,可以通过历史信息获取交通流量的长期时间依赖,提升模型的预测效果。
进一步,所述记忆输入层包括记忆模块和第一全连接层,输入的交通流量数据X0和流量数据的历史记忆依次经过记忆模块和第一全连接层,T为时间长度,N为节点数量,C1为节点历史信息特征的长度,R表示实数张量矩阵。
结构简单,便于使用。
进一步,所述记忆模块包括第二全连接层和GRU模块;
交通流量数据X0先经过第二全连接层,第二全连接层的输出和历史记忆/>作为GRU模块的输入;
GRU模块更新的历史记忆信息并保留当前交通流量序列的信息,得到新的历史记忆
将新的历史记忆作为新的特征拼接到当前节点的特征上,即将M′与X1进行拼接得到新的节点特征,其计算式为:
X2再经过第一全连接层来扩展节点的特征维度,以提高网络的表达能力,得到该层的输出H(0)∈RT×N×C,H(0)作为异质时空图卷积层的第一个子层的输入,C1为节点拼接历史信息特征后的长度。
GRU模块的目的是更新当前输入序列的信息到历史记忆上,将历史输入序列的历史特征(即历史记忆M)拼接到当前的输入序列上,以便后续使用。
进一步,所述异质时空图卷积层的子层包括异质时空图构建模块和时空图卷积模块,异质图构建模块进行图信号选择和异质图结构构建,用来构建不同时空异质图;
异质时空图构建模块输出T-2个异质时空图的图信号和对应的邻接矩阵,作为时空图卷积模块的输入,时空图卷积模块中包含多个图卷积块,图卷积块用于同步捕获异质时空图的时空依赖。
操作简单,利于使用。
进一步,异质时空图构建模块构建不同时空异质图的方法如下:
把相邻的3个图信号构建一个异质时空图,对于输入的时间序列H∈RT×N×C,即T个时刻的图信号,可以构建出T-2个异质时空图;
对于第i个异质时空图,其图信号为Hi=Reshape(H[i:i+2])∈R3N×C;
构建每个异质时空图的邻接矩阵,其过程如下:
根据输入的空间邻接矩阵,得到如下初始矩阵A:
其中,I∈RN×N是单位矩阵,这样构建的原因是在保持原来空间上的路网结构的同时,将中间时刻的节点与前后时刻的同一节点连接来构建时间上的联系;A代表异质时空图的基础连接关系,As∈RN×N为根据节点的空间关系生成的空间邻接矩阵,As不会随着时间发生变化;
邻接矩阵A与一个随机初始化可学习的参数矩阵Mask∈R3N×3M做哈达玛乘积来对已有的边进行边权调整,得到一个邻域关系矩阵Aneibor,其数学公式表示为:
Aneibor=A*Mask
随机初始化构建一个可学习的参数矩阵AEmb∈RT×N×E,该参数矩阵可以自适应地学习不同时刻节点的流量变化趋势,代表时间为i时,N个节点的流量变化趋势的嵌入表示;当构建第i个异质时空图的邻接矩阵时,选取i到i+2时刻的节点嵌入表示,即reshape是一个矩阵形状变换函数,/>表示第i到i+2时刻的节点的流量变化趋势的嵌入表示;E为超参数,为流量变化趋势嵌入表示的维度大小;通过将/>与/>的转置矩阵相乘来计算时空图上任意两点之间的流量变化趋势相似性,得到全局关系矩阵/>从而获取该异质时空图的节点之间的全局潜在关系;
将邻域关系矩阵和全局关系矩阵相加,并进行ReLu操作得到最终的邻接矩阵,其计算方式表示为:
其中,A′i为第i个异质时空图的邻接矩阵,每个异质时空图由于的不同,从而学习到不同的节点间时空关系来捕获交通流量数据的异质性。
对邻接矩阵A上的非0值做边权调整,用参数矩阵AEmb来得到不同时刻的图结构,来获取交通流量数据的异质性。
进一步,时空图卷积模块捕获异质时空图的时空依赖的方法如下:
图卷积模块中包含B个图卷积块,邻接矩阵使用异质时空图构建模块构建的异质时空图邻接矩阵,对于第i个异质时空图的图信号hi和邻接矩阵A′i,第b个图卷积块的计算公式为:
其中,αb和βb是两个超参数,I是单位矩阵,Wb∈RC×C是可学习的参数矩阵,σ为激活函数;为第b个图卷积块的输出,作为第b+1个图卷积块的输出,/>为输入的图信号hi;
对最后的图卷积块的输出进行剪切操作,仅保留中间时刻的图信号,得到图卷积模块的输出:
异质时空图卷积层拼接每个异质时空图通过时空图卷积模块得到的输出,即得到该异质时空图卷积层的输出:
H(l)=[h(0),h(1),...,h(T-2)]∈R(T-2)×N×C
H(l)是第l层的输出以及下一层异质时空图卷积层的输入。
图卷积模块中包含的图卷积块,图卷积块用于同步捕获异质时空图的时空依赖,运算简单,利于使用。
进一步,利用时空信息融合层,得到交通流量预测数据和新的历史记忆的方法如下:
拼接记忆输入层和异质时空图卷积层的结果,进行Reshape操作,即:
Xt=Reshape([H0,H1,...,HL])∈RN×KC
K为各层输出结果的时间维度的特征长度之和;
拼接结果Xt作为T′个预测模块的输入,每个预测模块包括两层全连接层,第i个预测模块用于预测未来第i个时刻的交通流量数据,其计算式为:
其中,是时刻i的预测结果,/>是可学习的参数,C′是第一层全连接层输出的维度,通过这样的方式,直接进行多步预测,避免单步预测中带来的误差传递;
拼接所有时间步的预测结果,得到未来T′时刻的交通流量预测数据:
交通流量预测数据和新的历史记忆M′作为预测输出。
将拼接结果用T个预测模块进行,因为考虑到不同范围的时空信息对未来不同时刻的区域的流量的影响是不一样的,所以针对不同的时刻使用不同的预测模块。做拼接和reshape操作也是将节点在中间层学习到的不同范围下的时空信息作为节点的特征,方便作为全连接层的输入。
本发明还提供一种持续记忆的自适应异质时空图卷积的交通预测系统,包括数据获取模块和预测模块,所述数据获取模块的输出端与预测模块的输入端连接,数据获取模块用于获取交通流量数据和流量数据的历史记忆,所述预测模块执行本发明所述方法,进行交通预测。
本方案利用预测模块进行交通预测,捕获交通流量数据的异质性,通过历史信息获取交通流量的长期时间依赖,提升预测模块的预测效果。
进一步,还包括评估模块,所述评估模块提取数据集,使用过去一定时间内的流量数据预测未来一定时间内的流量数据,划分数据集为训练集、验证集和测试集;
使用平均绝对误差,平均绝对百分比误差和均方根误差作为评价系统预测性能的评价指标,其计算方式如下所示:
为预测的交通流量结果,Y为真实的交通流量,yi为某一时刻某个节点的流量的真实值,/>为yi对应时刻对应节点的流量的预测值,n为所有真实值的数量。
利用评估模块评估预测模块的预测性能,判断预测模块的可靠性,以便进行使用或改进预测模块。
附图说明
图1是本发明持续记忆的自适应异质时空图卷积的交通预测方法的结构示意图;
图2是本发明持续记忆的自适应异质时空图卷积的交通预测系统的评估模块的评估结果示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
在本发明的描述中,需要理解的是,术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
在本发明的描述中,除非另有规定和限定,需要说明的是,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是机械连接或电连接,也可以是两个元件内部的连通,可以是直接相连,也可以通过中间媒介间接相连,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。
现有技术中,预测模型旨在学习时空数据的空间依赖和时间关系。ST-ResNet,ST-3DNet等方法采用CNN等方式获取交通流量的空间依赖,但是这些模型忽视了道路的拓扑网络结构。交通路网可以更加自然地表示为图结构,其主要构建方法是将道路上的传感器作为图的节点,并根据节点间的距离或是路网结构构建图的边集,从而模拟现实的交通路网结构。
近年来随着图神经网络(GNN)的发展,人们引入图神经网络和其变体到预测交通流量的模型中来捕获空间依赖,同时使用长短期记忆网络(LSTM)或是门控递归单元(GRU)等组件捕获时间依赖,这些模型很好地完成预测任务并取得了不错的效果。最新的模型尝试同时捕获时空关系,STG2Seq使用门控残差网络将相邻时间步上的节点特征进行拼接,尝试建模数据的时空依赖。STSGCN则是将相邻时刻的图信号构建成时空图,只使用图神经网络模型同时捕获时空图的时空依赖。
但现有的时空模型在捕获时间依赖的时候,往往只建立同一种图结构,忽视不同时间下,一个区域的流量受到其他区域流量的影响是不一样的现象,即交通流量数据的异质性。另外,以往的模型采用LSTM或是GRU等模块来获取交通流量数据的时间依赖,但其获取的时间依赖的时间跨度局限在当前输入数据的时间长度,不能记忆之前输入数据的历史信息。
本发明公开了一种持续记忆的自适应异质时空图卷积的交通预测方法,通过自适应矩阵,自适应构建不同时刻的时空图的图结构,捕获交通流量数据的异质性。并使用一个记忆模块将当前输入序列的信息保存到历史信息中,将生成的历史信息用于交通流量预测,历史信息会用到下一次流量预测中并进行更新,到达持续记忆的效果。
如图1所示,本发明的交通预测方法包括如下步骤:
获取交通流量数据和流量数据的历史记忆,并输入记忆输入层。记忆输入层输出时间序列,并将该时间序列作为异质时空图卷积层的第一个子层的输入。异质时空图卷积层设有多个(L个)子层,前一个子层的输出是下一个子层的输入,构建不同时空异质图,并利用时空异质图完成图卷积操作。异质时空图卷积层的各层输出时间序列至时空信息融合层,得到交通流量预测数据和新的历史记忆。
本发明的一种优选方案中,记忆输入层包括记忆模块和第一全连接层,输入的交通流量数据X0和流量数据的历史记忆依次经过记忆模块和第一全连接层,在某个时刻i的交通流量信息用/>表示,T为时间长度,N为节点数量,节点的特征为当前时刻当前节点的交通流量数据;M初始化为全0矩阵,表示不含历史信息,模型完成一次预测时会输出新的历史记忆,用于下一时刻的交通流量预测;C1为节点历史信息特征的长度,R表示实数张量矩阵。
记忆输入层中的记忆模块包括第二全连接层和GRU模块,交通流量数据X0先经过第二全连接层,第二全连接层的输出和历史记忆/>作为GRU模块的输入。GRU模块更新的历史记忆信息并保留当前交通流量序列的信息,得到新的历史记忆/>将新的历史记忆作为新的特征拼接到当前节点的特征上,即将M′与X1进行拼接得到新的节点特征,其计算式为:
X2再经过第一全连接层来扩展节点的特征维度,以提高网络的表达能力,得到该层的输出H(0)∈RT×N×C,H(0)作为异质时空图卷积层的第一个子层的输入,C1为节点拼接历史信息特征后的长度。GRU模块其目的是更新当前输入序列的信息到历史记忆上,将历史输入序列的历史特征(即历史记忆M)拼接到当前的输入序列上,而不是提取当前输入序列的时间依赖。
本发明的一种优选方案中,异质时空图卷积层的子层包括异质时空图构建模块和时空图卷积模块,异质图构建模块进行图信号选择和异质图结构构建,用来构建不同时空异质图。异质时空图构建模块输出T-2个异质时空图的图信号和对应的邻接矩阵,作为时空图卷积模块的输入,时空图卷积模块中包含多个图卷积块,图卷积块用于同步捕获异质时空图的时空依赖。
异质时空图卷积层中的异质时空图构建模块构建不同时空异质图的方法如下:
把相邻的3个图信号构建一个异质时空图,对于输入的时间序列H∈RT×N×C,即T个时刻的图信号,可以构建出T-2个异质时空图,过去T时间的历史数据可用序列表示,模型输入该序列预测未来T′时间的数据序列
对于第i个异质时空图,其图信号为Hi=Reshape(H[i:i+2])∈R3N×C;
构建每个异质时空图的邻接矩阵,其过程如下:
根据输入的空间邻接矩阵,得到如下初始矩阵A:
其中,J∈RN×N是单位矩阵,这样构建的原因是在保持原来空间上的路网结构的同时,将中间时刻的节点与前后时刻的同一节点连接来构建时间上的联系;A代表异质时空图的基础连接关系,As∈RN×N为根据节点的空间关系生成的空间邻接矩阵,As不会随着时间发生变化;As上元素的取值过程如下所示:当前节点i与节点j之间的距离小于人为设定的阈值时,该元素的值取为1,否则为0,1表示节点i和节点j相连,0表示不相连;
邻接矩阵A与一个随机初始化可学习的参数矩阵Mask∈R3N×3N做哈达玛乘积来对已有的边进行边权调整,做哈达玛乘积两个矩阵的大小一致,A的大小是3N x 3N,所以参数矩阵Mask的大小也设定为3N x 3N,Mask中的值一开始是随机初始化的,在模型的训练中通过反向传播机制进行更新矩阵上的值。得到一个邻域关系矩阵Aneibor,其数学公式表示为:
Aneibor=A*Mask
对邻接矩阵A上的非0值做边权调整,但是之前的方法往往使用这种方式构建得到邻接矩阵A后就不再变化了,使用不变的邻接矩阵做图卷积操作难以建模流量数据的异质性,所以后面的部分是对原来的方法进行改进,用参数矩阵AEmb来得到不同时刻的图结构,来获取交通流量数据的异质性;
随机初始化构建一个可学习的参数矩阵AEmb∈RT×N×E,创建时矩阵中的值一开始是随机初始化的,在模型的训练中通过反向传播机制进行更新。其形状根据输入该子层的时间序列的形状确定。该参数矩阵可以自适应地学习不同时刻节点的流量变化趋势,代表时间为i时,N个节点的流量变化趋势的嵌入表示;当构建第i个异质时空图的邻接矩阵时,选取i到i+2时刻的节点嵌入表示,即/>reshape是一个矩阵形状变换函数,不会改变矩阵上的值;/>表示第i到i+2时刻的节点的流量变化趋势的嵌入表示,虽然形状经过reshape后发生了改变,但含义没变化;E为超参数,为流量变化趋势嵌入表示的维度大小;通过将/>与/>的转置矩阵相乘来计算时空图上任意两点之间的流量变化趋势相似性,得到全局关系矩阵/>从而获取该异质时空图的节点之间的全局潜在关系;
将邻域关系矩阵和全局关系矩阵相加,并进行ReLu操作得到最终的邻接矩阵,其计算方式表示为:
其中,A′i为第i个异质时空图的邻接矩阵,每个异质时空图由于的不同,从而学习到不同的节点间时空关系来捕获交通流量数据的异质性。
本发明的一种优选方案中,时空图卷积模块捕获异质时空图的时空依赖的方法如下:
异质时空图构建模块输出T-2个异质时空图的图信号和对应的邻接矩阵,作为时空图卷积模块的输入。图卷积模块中包含B个图卷积块,图卷积块用于同步捕获异质时空图的时空依赖。图卷积块根据现有GCNII改进而来,邻接矩阵不再使用归一化的图拉普拉斯矩阵,而是使用异质时空图构建模块构建的异质时空图邻接矩阵,对于第i个异质时空图的图信号hi和邻接矩阵A′i,第b个图卷积块的计算公式为:
其中,αb和βb是两个超参数,I是单位矩阵,Wb∈RC×C是可学习的参数矩阵,σ为激活函数;为第b个图卷积块的输出,作为第b+1个图卷积块的输出,/>为输入的图信号hi;
对最后的图卷积块(即第B个图卷积块)的输出进行剪切操作,仅保留中间时刻的图信号,得到图卷积模块的输出:
异质时空图卷积层拼接每个异质时空图通过时空图卷积模块得到的输出,即得到该异质时空图卷积层的输出:
H(l)=[h(0),h(1),...,h(T-2)]∈R(T-2)×N×C
H(l)是第l层的输出以及下一层异质时空图卷积层的输入。
本发明的一种优选方案中,利用时空信息融合层(即图1中的时空融合输出层),得到交通流量预测数据和新的历史记忆的方法如下:
拼接记忆输入层和异质时空图卷积层的结果,进行Reshape操作,即:
Xt=Reshape([H0,H1,...,HL])∈RN×KC
K为各层输出结果的时间维度的特征长度之和;
拼接结果Xt作为T′个预测模块的输入,每个预测模块包括两层全连接层,第i个预测模块用于预测未来第i个时刻的交通流量数据,其计算式为:
其中,是时刻i的预测结果,/>是可学习的参数;C′是第一层全连接层输出的维度,通过这样的方式,直接进行多步预测,避免单步预测中带来的误差传递;将拼接结果用T个预测模块进行,因为考虑到不同范围的时空信息对未来不同时刻的区域的流量的影响是不一样的,所以针对不同的时刻使用不同的预测模块,同时这也是为什么会在前面保留中间层的结果(即不同范围的时空信息)的原因;做拼接和reshape操作也是将节点在中间层学习到的不同范围下的时空信息作为节点的特征,方便作为全连接层的输入;
拼接所有时间步的预测结果,得到未来T′时刻的交通流量预测数据:
交通流量预测数据和新的历史记忆M′作为预测输出。
本发明还提供一种持续记忆的自适应异质时空图卷积的交通预测系统,包括数据获取模块和预测模块,数据获取模块的输出端与预测模块的输入端电性连接,数据获取模块用于获取交通流量数据和流量数据的历史记忆,预测模块执行本发明所述方法,进行交通预测。
本发明的一种优选方案中,交通预测系统还包括评估模块,评估模块提取数据集,数据集可采用公开的数据集PEMS03和PEMS08,两个数据集为美国加州高速公路不同区域的交通流量数据。PEMS03包含了350个节点,PEMS08包含了170个节点。使用过去一定时间内的流量数据预测未来一定时间内的流量数据,以6:2:2的比例划分数据集为训练集、验证集和测试集;
使用平均绝对误差(MAE),平均绝对百分比误差(MAPE)和均方根误差(RMSE)作为评价系统预测性能的评价指标,其计算方式如下所示:
为预测的交通流量结果,Y为真实的交通流量,yi为某一时刻某个节点的流量的真实值,/>为yi对应时刻对应节点的流量的预测值,n为所有真实值的数量;评估结果如图2所示。评价指标为mae,mape和rmse,计算如上面的公式所示,三个指标的值越小,说明预测结果越准确。对比模型为STGCN、DCRNN、STSGCN,是目前最先进的模型。本方案预测模块的模型为MASTGCN,在三个评价指标上均小于上述模型,说明了本方案的预测效果均优于上面的三个模型。
本方案:一是通过可学习的自适应矩阵AEmb为不同时刻的时空图构建不同的图结构以提取流量数据的异质性,而现有方法是只使用一种固定的图结构用于图卷积,这样难以捕获异质性。
二是提出一种持续记忆的方法,即记忆输入层中的记忆模块,其中用到的GRU是用于将当前输入序列的特征更新到历史记忆中,然后记忆模块还会将得到历史记忆拼接到当前输入序列的特征上用于当前的预测,到达持续记忆的效果,并将历史记忆作为模型的输出用于下次输入。通过历史信息获取交通流量的长期时间依赖,提升模型的预测效果。而现有方法都是利用GRU等方式提取当前输入序列中的数据的时间相关性,与之前的输入序列无关系,本方法用到的GRU并不是用于提取当前输入序列中的数据的时间相关性。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同物限定。
Claims (6)
1.一种持续记忆的自适应异质时空图卷积的交通预测方法,其特征在于,包括如下步骤:
获取交通流量数据和流量数据的历史记忆,并输入记忆输入层;
记忆输入层输出时间序列,并将该时间序列作为异质时空图卷积层的第一个子层的输入;
异质时空图卷积层设有多个子层,前一个子层的输出是下一个子层的输入,构建不同时空异质图,并利用时空异质图完成图卷积操作;
异质时空图卷积层的各层输出时间序列至时空信息融合层,得到交通流量预测数据和新的历史记忆;
所述异质时空图卷积层的子层包括异质时空图构建模块和时空图卷积模块,异质图构建模块进行图信号选择和异质图结构构建,用来构建不同时空异质图;
异质时空图构建模块输出T-2个异质时空图的图信号和对应的邻接矩阵,作为时空图卷积模块的输入,时空图卷积模块中包含多个图卷积块,图卷积块用于同步捕获异质时空图的时空依赖;
异质时空图构建模块构建不同时空异质图的方法如下:
把相邻的3个图信号构建一个异质时空图,对于输入的时间序列H∈RT×N×C,即T个时刻的图信号,可以构建出T-2个异质时空图;
对于第i个异质时空图,其图信号为Hi=Reshape(H[i:i+2])∈R3N×C;
构建每个异质时空图的邻接矩阵,其过程如下:
根据输入的空间邻接矩阵,得到如下初始矩阵A:
其中,I∈RN×N是单位矩阵,这样构建的原因是在保持原来空间上的路网结构的同时,将中间时刻的节点与前后时刻的同一节点连接来构建时间上的联系;A代表异质时空图的基础连接关系,As∈RN×N为根据节点的空间关系生成的空间邻接矩阵,As不会随着时间发生变化;
邻接矩阵A与一个随机初始化可学习的参数矩阵Mask∈R3N×3N做哈达玛乘积来对已有的边进行边权调整,得到一个邻域关系矩阵Aneibor,其数学公式表示为:
Aneibor=A*Mask
随机初始化构建一个可学习的参数矩阵AEmb∈RT×N×E,该参数矩阵可以自适应地学习不同时刻节点的流量变化趋势,代表时间为i时,N个节点的流量变化趋势的嵌入表示;当构建第i个异质时空图的邻接矩阵时,选取i到i+2时刻的节点嵌入表示,即reshape是一个矩阵形状变换函数,/>表示第i到i+2时刻的节点的流量变化趋势的嵌入表示;E为超参数,为流量变化趋势嵌入表示的维度大小;通过将/>与/>的转置矩阵相乘来计算时空图上任意两点之间的流量变化趋势相似性,得到全局关系矩阵/>从而获取该异质时空图的节点之间的全局潜在关系;
将邻域关系矩阵和全局关系矩阵相加,并进行ReLu操作得到最终的邻接矩阵,其计算方式表示为:
其中,A′i为第i个异质时空图的邻接矩阵,每个异质时空图由于的不同,从而学习到不同的节点间时空关系来捕获交通流量数据的异质性;
时空图卷积模块捕获异质时空图的时空依赖的方法如下:
图卷积模块中包含B个图卷积块,邻接矩阵使用异质时空图构建模块构建的异质时空图邻接矩阵,对于第i个异质时空图的图信号hi和邻接矩阵A′i,第b个图卷积块的计算公式为:
其中,αb和βb是两个超参数,I是单位矩阵,Wb∈RC×C是可学习的参数矩阵,σ为激活函数;为第b个图卷积块的输出,作为第b+1个图卷积块的输出,/>为输入的图信号hi;
对最后的图卷积块的输出进行剪切操作,仅保留中间时刻的图信号,得到图卷积模块的输出:
异质时空图卷积层拼接每个异质时空图通过时空图卷积模块得到的输出,即得到该异质时空图卷积层的输出:
H(l)=[h(0),h(1),...,h(T-2)]∈R(T-2)×N×C
H(l)是第l层的输出以及下一层异质时空图卷积层的输入。
2.如权利要求1所述的持续记忆的自适应异质时空图卷积的交通预测方法,其特征在于,所述记忆输入层包括记忆模块和第一全连接层,输入的交通流量数据X0和流量数据的历史记忆依次经过记忆模块和第一全连接层,T为时间长度,N为节点数量,C1为节点历史信息特征的长度,R表示实数张量矩阵。
3.如权利要求2所述的持续记忆的自适应异质时空图卷积的交通预测方法,其特征在于,所述记忆模块包括第二全连接层和GRU模块;
交通流量数据X0先经过第二全连接层,第二全连接层的输出和历史记忆作为GRU模块的输入;
GRU模块更新的历史记忆信息并保留当前交通流量序列的信息,得到新的历史记忆
将新的历史记忆作为新的特征拼接到当前节点的特征上,即将M′与X1进行拼接得到新的节点特征,其计算式为:
X2再经过第一全连接层来扩展节点的特征维度,以提高网络的表达能力,得到该层的输出H(0)∈RT×N×C,H(0)作为异质时空图卷积层的第一个子层的输入,C为节点拼接历史信息特征后的长度。
4.如权利要求1所述的持续记忆的自适应异质时空图卷积的交通预测方法,其特征在于,利用时空信息融合层,得到交通流量预测数据和新的历史记忆的方法如下:
拼接记忆输入层和异质时空图卷积层的结果,进行Reshape操作,即:
Xt=Reshape([H0,H1,...,HL])∈RN×KC
其中,K为各层输出结果的时间维度的特征长度之和,HL是第L层异质时空图卷积层的输出;
拼接结果Xt作为T′个预测模块的输入,每个预测模块包括两层全连接层,第i个预测模块用于预测未来第i个时刻的交通流量数据,其计算式为:
其中,是时刻i的预测结果,/> 是可学习的参数,C′是第一层全连接层输出的维度,通过这样的方式,直接进行多步预测,避免单步预测中带来的误差传递;
拼接所有时间步的预测结果,得到未来T′时刻的交通流量预测数据:
交通流量预测数据和新的历史记忆M′作为预测输出。
5.一种持续记忆的自适应异质时空图卷积的交通预测系统,其特征在于,包括数据获取模块和预测模块,所述数据获取模块的输出端与预测模块的输入端连接,数据获取模块用于获取交通流量数据和流量数据的历史记忆,所述预测模块执行权利要求1-4之一所述方法,进行交通预测。
6.如权利要求5所述的持续记忆的自适应异质时空图卷积的交通预测系统,其特征在于,还包括评估模块,所述评估模块提取数据集,使用过去一定时间内的流量数据预测未来一定时间内的流量数据,划分数据集为训练集、验证集和测试集;
使用平均绝对误差,平均绝对百分比误差和均方根误差作为评价系统预测性能的评价指标,其计算方式如下所示:
为预测的交通流量结果,Y为真实的交通流量,yi为某一时刻某个节点的流量的真实值,/>为yi对应时刻对应节点的流量的预测值,n为所有真实值的数量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210014979.2A CN114169647B (zh) | 2022-01-07 | 2022-01-07 | 持续记忆的自适应异质时空图卷积的交通预测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210014979.2A CN114169647B (zh) | 2022-01-07 | 2022-01-07 | 持续记忆的自适应异质时空图卷积的交通预测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114169647A CN114169647A (zh) | 2022-03-11 |
CN114169647B true CN114169647B (zh) | 2024-05-07 |
Family
ID=80489085
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210014979.2A Active CN114169647B (zh) | 2022-01-07 | 2022-01-07 | 持续记忆的自适应异质时空图卷积的交通预测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114169647B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115115094B (zh) * | 2022-05-26 | 2024-02-06 | 浙江省机电设计研究院有限公司 | 结合序列局部信息与多序列关联关系的交通流量预测方法 |
CN115240425B (zh) * | 2022-07-26 | 2023-07-04 | 西北师范大学 | 一种基于多尺度时空融合图网络的交通预测方法 |
CN115565370B (zh) * | 2022-09-20 | 2024-04-23 | 北京工商大学 | 一种局部时空图卷积交通流量预测方法及系统 |
CN116153089B (zh) * | 2023-04-24 | 2023-06-27 | 云南大学 | 基于时空卷积与动态图的交通流量预测系统及方法 |
CN117151285B (zh) * | 2023-08-29 | 2024-07-02 | 淮阴工学院 | 一种基于多要素注意力时空图卷积网络的径流预报方法 |
CN118247968A (zh) * | 2024-04-15 | 2024-06-25 | 青岛大学 | 基于异质图预测未来交通流的方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112749832A (zh) * | 2020-12-03 | 2021-05-04 | 南京邮电大学 | 基于图卷积嵌入时空长短时记忆网络的交通流量预测方法 |
CN112784121A (zh) * | 2021-01-28 | 2021-05-11 | 浙江工业大学 | 一种基于时空图表示学习的交通事故预测方法 |
CN113053115A (zh) * | 2021-03-17 | 2021-06-29 | 北京工商大学 | 一种基于多尺度图卷积网络模型的交通预测方法 |
CN113450568A (zh) * | 2021-06-30 | 2021-09-28 | 兰州理工大学 | 一种基于时空注意力机制的卷积网络交通流预测方法 |
CN113780662A (zh) * | 2021-09-14 | 2021-12-10 | 北京航空航天大学 | 一种流量预测方法、装置、设备和介质 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112101682B (zh) * | 2020-09-25 | 2024-04-09 | 北京百度网讯科技有限公司 | 流量模式预测方法、装置、服务器以及可读介质 |
-
2022
- 2022-01-07 CN CN202210014979.2A patent/CN114169647B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112749832A (zh) * | 2020-12-03 | 2021-05-04 | 南京邮电大学 | 基于图卷积嵌入时空长短时记忆网络的交通流量预测方法 |
CN112784121A (zh) * | 2021-01-28 | 2021-05-11 | 浙江工业大学 | 一种基于时空图表示学习的交通事故预测方法 |
CN113053115A (zh) * | 2021-03-17 | 2021-06-29 | 北京工商大学 | 一种基于多尺度图卷积网络模型的交通预测方法 |
CN113450568A (zh) * | 2021-06-30 | 2021-09-28 | 兰州理工大学 | 一种基于时空注意力机制的卷积网络交通流预测方法 |
CN113780662A (zh) * | 2021-09-14 | 2021-12-10 | 北京航空航天大学 | 一种流量预测方法、装置、设备和介质 |
Non-Patent Citations (3)
Title |
---|
Adaptive Spatial-Temporal Fusion Graph Convolutional Networks for Traffic Flow Forecasting;Senwen Li 等;《2022 International Joint Conference on Neural Networks》;20220930;第1-8页 * |
基于注意力机制的交通流时空预测模型研究;李林杰;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20210215(第2期);第C034-1076页 * |
面向交通流量预测的跨邻域时空图卷积算法研究;黎森文;《万方学位论文》;20231101;第1-59页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114169647A (zh) | 2022-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114169647B (zh) | 持续记忆的自适应异质时空图卷积的交通预测方法及系统 | |
CN109697852B (zh) | 基于时序交通事件的城市道路拥堵程度预测方法 | |
CN112215337B (zh) | 一种基于环境注意力神经网络模型的车辆轨迹预测方法 | |
CN110164128B (zh) | 一种城市级智能交通仿真系统 | |
CN109887282B (zh) | 一种基于层级时序图卷积网络的路网交通流预测方法 | |
CN111639787B (zh) | 一种基于图卷积网络的时空数据预测方法 | |
CN113313947B (zh) | 短期交通预测图卷积网络的路况评估方法 | |
CN114299723B (zh) | 一种交通流量预测方法 | |
CN114802296A (zh) | 一种基于动态交互图卷积的车辆轨迹预测方法 | |
CN113643528B (zh) | 信号灯控制方法、模型训练方法、系统、装置及存储介质 | |
CN111710162B (zh) | 一种城市路网交通运行状况监测方法及系统 | |
CN110737968A (zh) | 基于深层次卷积长短记忆网络的人群轨迹预测方法及系统 | |
CN110570035B (zh) | 同时建模时空依赖性和每日流量相关性的人流量预测系统 | |
CN114925836B (zh) | 基于动态多视图图神经网络的城市交通车流量推理方法 | |
CN116935649A (zh) | 多视角融合的时空动态图卷积网络城市交通流量预测方法 | |
CN115578851A (zh) | 一种基于mgcn的交通预测方法 | |
CN115376317B (zh) | 一种基于动态图卷积和时序卷积网络的交通流预测方法 | |
CN115457081A (zh) | 一种基于图神经网络的分层融合式预测方法 | |
CN115641720A (zh) | 基于时空融合图神经网络的交通预测方法及系统 | |
CN114582131B (zh) | 一种基于匝道智能控流算法的监控方法及系统 | |
CN117709602B (zh) | 一种基于社会价值取向的城市智能车辆拟人化决策方法 | |
CN117291780A (zh) | 一种周期自适应图卷积循环网络交通流预测方法 | |
CN117217381A (zh) | 一种属性增强的感知时序图卷积网络交通速度预测方法 | |
CN117456730A (zh) | 一种基于自适应部分注意力扩散图卷积神经网络的交通流量预测方法 | |
CN113312442B (zh) | 一种智慧城市电子地图生成方法和系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |