CN113313947B - 短期交通预测图卷积网络的路况评估方法 - Google Patents
短期交通预测图卷积网络的路况评估方法 Download PDFInfo
- Publication number
- CN113313947B CN113313947B CN202110603374.2A CN202110603374A CN113313947B CN 113313947 B CN113313947 B CN 113313947B CN 202110603374 A CN202110603374 A CN 202110603374A CN 113313947 B CN113313947 B CN 113313947B
- Authority
- CN
- China
- Prior art keywords
- traffic
- data
- time
- convolution
- value
- 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
Images
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
- 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
- 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
- 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/052—Detecting movement of traffic to be counted or controlled with provision for determining speed or overspeed
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/065—Traffic control systems for road vehicles by counting the vehicles in a section of the road or in a parking area, i.e. comparing incoming count with outgoing count
Abstract
本发明涉及一种短期交通预测图卷积网络的路况评估方法,属于智能交通领域。通过将表示交通网络中节点的距离二元关系的邻接矩阵抽象为模糊关系矩阵,并通过计算出其传递闭包得到表达所有潜在节点之间联系的转移关系矩阵,并结合注意力机制的自学习矩阵得到基于当前输入的转移关系矩阵,然后将此矩阵用于图卷积过程。有益效果在于,通过权重学习表示出交通网络中数据的转移关系,提高了图卷积网络对交通数据流中动态空间依赖的捕捉能力。同时模型通过构建三类数据集考虑了交通流量中长时间跨度依赖,并通过因果空洞卷积考虑短时间跨度,提高了交通数据预测的准确性。
Description
技术领域
本发明涉及一种短期交通预测图卷积网络的路况评估方法,属于智能交通和人工智能技术领域。
背景技术
在汽车工业水平大幅进步和居民经济水平日益提高的背景下,我国居民汽车保有量持续增加,由此带来交通拥堵、环境污染等问题也日益严重。发展绿色高效的智能交通系统是国家战略也是民心所向。交通流量是智能交通系统分析的重要数据,因为区域的交通流量情况可以很好地反映出该区域当前的交通状况。使用交通流历史数据分析出交通流的潜在规律从而对未来交通流量做出预测可以帮助交通管理者有效地提前干预交通、帮助参与者选择更优的行车路线从而缓解交通拥堵。但历史交通数据本身具备的高维度,高复杂性等特点,对其进行有效预测分析是一个有挑战的课题。
现有的道路交通流预测方法大致可以分为两种,一种是基于统计模型的方法,包括历史平均回归、卡尔曼滤波法、自向量回归等,这类方法受时间序列静态假设的限制,没有考虑时空相关性。另一种是近年来兴起的深度学习模型,这类方法往往能更好的捕捉交通流数据中复杂的非线性时空关系,从而达到更好的效果。但时空关系可能是动态变化的,不同路段在不同的时间节点彼此之间的影响可能是不同的,这种动态特征通常没有被有效捕捉。对交通路况的评估往往需要对交通网络的先验知识,本发明基于历史数据估计了两种对交通流状况的评测指标,可应用于不同的数据集,提高了深度模型的可应用性。
发明内容
针对交通流数据中动态时空特征以及提高交通流数据预测的准确性,本发明提出了一种短期交通预测图卷积网络的路况评估方法,从而克服现有技术中的不足。
本发明的原理在于,通过将表示交通网络中节点的距离二元关系的邻接矩阵抽象为模糊关系矩阵,并通过计算出其传递闭包得到表达所有潜在节点之间联系的转移关系矩阵,并结合注意力机制的自学习矩阵得到基于当前输入的转移关系矩阵,然后将此矩阵用于图卷积过程。
本发明通过以下技术方案实现,包括以下步骤:
步骤一:获取交通数据并进行预处理:获取交通路段的空间距离信息,用于构建交通网络的图结构;获取的交通流数据建立交通流数据集,提取每个节点的自由行驶速度和最大自由通行流量,然后将历史交通数据按交通节点数N、特征数量F、时间序列T三个维度整理成N×F×T的多维数据集,对缺失数据进行缺失补齐,并进行归一化处理最后划分训练集和测试集;
步骤1.2:提取每个节点的自由行驶速度Vfree和最大自由通行流量MaxFfree;自由行驶速度对历史数据均匀取样后按各节点凌晨时段的历史平均速度数据去噪后平均估计得到,最大自由通行流量按各节点在自由行驶速度向上向下取整范围内的速度对应的流量数据集合取部分最大值平均估计得到;
步骤1.3:建立交通流数据集,数据集结构表示如下:
缺失数据补齐规则:在连续缺失值少于等于3个时采用线性插值,大于3个时采用上个数值、昨天数值与上周数值加权平均。数据补齐后按照采用0均值归一化方法对数据集进行归一化处理,归一化公式如下:
x=(x-mean)/std (3)
其中mean表示该交通特征数据的平均值,std表示该交通特征数据的标准差;
步骤二:对多维数据集划分为小时趋势数据、天趋势数据、周趋势数据:三类数据具有相同的目标值,小时趋势的输入值为目标值时段前几个小时时间段的数据,天趋势的输入值为前几天相同目标时间段与小时趋势同长度的前几个小时的数据,周趋势的输入值为近几周前相同目标时间段的与小时趋势同长度的前几个小时的数据;
步骤2.1周趋势数据、天趋势数据、小时趋势数据从步骤一中的数据集产生;
假设预测时间长度为p,则目标值为[Xt,Xt+1,...,Xt+p],生成规则如下:
小时趋势数据:
Xh=[Xt-n,Xt-n+1,...,Xt-1]
其中n表示取与目标值相邻的前n个时间步;
天趋势数据:
Xd
=[X(t-n)-dt*dn,X(t-n+1)-dt*dn,...,X(t-1)-dt*dn;X(t-n)-dt*(dn-1),X(t-n+1)-dt*(dn-1),...,X(t-1)-dt*(dn-1);......;X(t-n)-dt,X(t-n+1)-dt,...,X(t-1)-dt]
其中dt为一天的时间长度,dn表示取与目标值相邻的前dn天;
周趋势数据:
Xw
=[X(t-n)-wt*wn,X(t-n+1)-wt*wn,...,X(t-1)-wt*wn;X(t-n)-wt*(wn-1),X(t-n+1)-wt*(wn-1),...,X(t-1)-wt*(wn-1);......;X(t-n)-wt,X(t-n+1)-wt,...,X(t-1)-wt]
其中wt为一周的时间长度,wn表示取与目标值相邻的前wn个周;
步骤三:构建基于模糊关系的短期交通预测图卷积网络:构造时空卷积单元,每个单元包括时间注意力模块、空间注意力模块、模糊关系图卷积模块、门控空洞因果时间卷积模块;通过时空卷积单元的堆叠,完成对交通流张量的时空特征提取,最后构建全连接层通过relu非线性回归得到预测结果;
步骤3.1:时间注意力模块用于捕获输入数据之间的动态时间相关性,即找到当前输入X中相关性较高的时间步:
P=Vp·tanh((XTU1)U2(U3X)+bp) (4)
采用tanh激活函数;
P中各元素的值表示输入时间序列中不同时间步长之间的动态联系程度;然后对P使用softmax函数进行归一化处理,然后用P乘以输入X得到根据时间相关性动态调整后的输入X;
步骤3.2:空间注意力模块:
S=Vs·tanh((X W1)W2(W3XT)+bs) (5)
其中为空间注意力矩阵,N为交通节点数;Vs,W1,W2,W3,bS为可学习的参数,采用tanh激活函数;S中各元素的值表示当前输入中不同路段之间的动态联系程度;然后对S使用softmax函数进行归一化处理;
采用K阶近似切比雪夫多项式图卷积来获取交通流数据中各路段之间空间依赖:
其中⊙表示哈达玛积;
步骤3.4时间卷积模块:
采用K层门控空洞因果卷积来获取交通流数据中路段的时间依赖;
其中*表示空洞因果卷积操作,Φ表示卷积核参数,tanh是输出的激活函数,sigmiod函数决定了当前输入传递给下一层的比率;K层空洞因果卷积结构通过跳过部分序列扩大卷积的感受野,同时使序列变得更加光滑,增强模型鲁棒性;
步骤四:构建三个并行模型:三个并行模型具有相同的网络结构,每个网络结构由时空卷积单元堆叠,最后一层添加全连接层构成;三个并行模型分别处理周趋势数据、天趋势数据和时刻趋势数据,最后将得到的三个预测结果通过参数加权得到最终预测结果;融合过程如下:
步骤五:训练深度模型优化参数实现短期交通流的实时预测:将训练集作为深度模型的输入,定义模型的损失函数,根据损失函数数值利用反向传播算法不断更新模型参数;待参数模型收敛后利用训练好的模型对测试集进行实时交通流预测;
定义模型的损失函数,训练模型:
步骤六:实时交通流的路况评估:将实时交通流数据输入到上述训练好的模型,得到预测车速和流量,再根据自由行驶速度和最大自由通行流量计算行驶时间指标TTI和道路流量指标RFI,按照两个指标对交通状态的隶属度函数以及专家经验进行交通状况进行模糊综合评价;步骤6.1:TTI,RFI指标两者计算公式如下:
其中Vfree自由通行速度,Vcurrent为当前速度;
其中MaxFfree表示最大自由通行流量,Fcurrent表示当前流量;
TTI指标的理论取值范围为[1,+∞],由于自由通行速度为平均估计数据以及考虑到超速情况,对于小于1的值按向上取整处理,RFI指标的理论取值范围为[1,+∞],将两者归一化按下式归一化处理:
步骤6.2:根据模糊理论,得到因素论域U={TTI,RFI},本发明将交通等级论域:
V={完全畅通,基本畅通,轻微拥堵,中度拥堵,严重拥堵},采用TTI和RFI的隶属度函数,高斯隶属的函数的公式为:
在TTI中的三个高斯隶属度函数的k取值依次为[0.3,0.5,0.7],r的取值依次为[0.07,0.08,0.08];
在RFI中的三个高斯隶属度函数的k取值依次为[0.25,0.4,0.55],r的取值依次为[0.07,0.08,0.09];根据隶属度函数计算得到的TTI和RFI指标得到模糊关系矩阵R:
rij表示模糊关系矩阵R中第i个指标对第j个交通等级的隶属度。随后根据专家经验确定评判因素权向量A={a1,a2},ai是U中第i个因素对交通等级评判的重要性分配权重,设置为0.5,0.5;选择合适的评论合成算子,将A与R合成得到综合评判B=(b1,b2,b3,b4,b5),本发明采用算子,其中运算符表示求和取与1比较的较小值:
其中B为模糊综合评价论域,bi为评价隶属值。
所述步骤一中,F表示收集的交通流数据特征总数,至少包括车流量、平均速度。
本发明的有益效果在于,通过权重学习表示出交通网络中数据的转移关系,提高了图卷积网络对交通数据流中动态空间依赖的捕捉能力。同时模型通过构建三类数据集考虑了交通流量中长时间跨度依赖,并通过因果空洞卷积考虑短时间跨度,提高了交通数据预测的准确性。
附图说明
图1为一种基于深度模型的路况评估方法流程图。
图2为交通流预测深度模型结构图。
图3为时空卷积单元组成图。
图4为TTI隶属度函数。
图5为RFI隶属度函数。
具体实施方式
下面结合附图1至5对本发明的优选实施例作进一步说明,本发明方法流程如图1所示,主要包括:
步骤一:获取交通数据并进行预处理:获取交通路段的空间距离信息,用于构建交通网络的图结构;获取的交通流数据建立交通流数据集,提取每个节点的自由行驶速度和最大自由通行流量,然后将历史交通数据按节点数量N、特征数量F、时间序列T三个维度整理成N×F×T的多维数据集,数据特征包括路段车辆平均速度、车流量等,对缺失数据进行缺失补齐,并进行归一化处理最后划分训练集和测试集;
步骤1.1:首先获取交通各路段之间的空间距离信息,用于构建交通网络有向加权图结构,图表示如下:
其中表示图的顶点集合,即交通网络中的所有路段,ε表示图的边集合,即交通网络中路段之间的距离信息。用来表示的自循环有向加权邻接矩阵,的对角线元素为1,如果是ε中的一条边,则等于和之间的归一化距离,本发明使用最大最小归一化方法;
步骤1.2:提取每个节点的自由行驶速度Vfree和最大自由通行流量MaxFfree;自由行驶速度对历史数据均匀取样后按各节点凌晨0点到4点的历史平均速度数据去噪后平均估计得到,最大自由通行流量按各节点在自由行驶速度向上向下取整范围内的速度对应的流量数据集合取前的最大值平均估计得到;
步骤1.3:建立交通流数据集x,数据集结构表示如下:
其中N表示交通节点数,即的大小;F表示收集的交通流数据特征总数,例如车流量、平均速度、占有率等;T表示时间序列;用表示路段i在t时刻的所有交通特征数据;缺失数据补齐规则:在连续缺失值少于等于3个时采用线性插值,大于3个时采用上个数值、昨天数值与上周数值加权平均;数据补齐后按照采用0均值归一化方法对数据集进行归一化处理,归一化公式如下:
x=(x-mean)/std (3)
其中mean表示该交通特征数据的平均值,std表示该交通特征数据的标准差。
步骤二:对多维数据集划分为小时趋势数据、天趋势数据、周趋势数据:
三类数据具有相同的目标值,小时趋势的输入值为目标值时段前几个小时时间段的数据,天趋势的输入值为前几天相同目标时间段与小时趋势同长度的前几个小时的数据,周趋势的输入值为近几周前相同目标时间段的与小时趋势同长度的前几个小时的数据;
步骤2.1:周趋势数据、天趋势数据、小时趋势数据从步骤一中的数据集产生;
假设预测时间长度为p,则目标值为[Xt,Xt+1,...,Xt+p],生成规则如下:
小时趋势数据:
Xh=[Xt-n,Xt-n+1,...,Xt-1]
其中n表示取与目标值相邻的前n个时间步;
天趋势数据:
Xd
=[X(t-n)-dt*dn,X(t-n+1)-dt*dn,...,X(t-1)-dt*dn;X(t-n)-dt*(dn-1),X(t-n+1)-dt*(dn-1),...,X(t-1)-dt*(dn-1);......;X(t-n)-dt,X(t-n+1)-dt,...,X(t-1)-dt]
其中dt为一天的时间长度,dn表示取与目标值相邻的前dn天;
周趋势数据:
Xw
=[X(t-n)-wt*wn,X(t-n+1)-wt*wn,...,X(t-1)-wt*wn;X(t-n)-wt*(wn-1),X(t-n+1)-wt*(wn-1),...,X(t-1)-wt*(wn-1);......;X(t-n)-wt,X(t-n+1)-wt,...,X(t-1)-wt]
其中wt为一周的时间长度,wn表示取与目标值相邻的前wn个周。
步骤三:构建基于模糊关系的短期交通预测图卷积网络:构造时空卷积单元,其组成如图3所示,每个单元包括时间注意力模块、空间注意力模块、模糊关系图卷积模块、门控空洞因果时间卷积模块;通过时空卷积单元的堆叠,完成对交通流张量的时空特征提取,最后构建全连接层通过relu非线性回归得到预测结果;
步骤3.1:时间注意力模块用于捕获输入数据之间的动态时间相关性,即找到当前输入X中相关性较高的时间步:
P=Vp•tanh((XT U1)U2(U3X)+bp) (4)
其中为时间注意力矩阵,T为输入的时间步长。Vp,U1,U2,U3,bp为可学习的参数,采用tanh激活函数;P中各元素的值表示输入时间序列中不同时间步长之间的动态联系程度;对P使用softmax函数进行归一化处理,然后用P乘以输入X得到根据时间相关性动态调整后的输入X;
步骤3.2:空间注意力模块:
S=Vs•tanh((X W1)W2(W3XT)+bs) (5)
其中为空间注意力矩阵,N为交通节点数;Vs,W1,W2,W3,bS为可学习的参数,采用tanh激活函数;S中各元素的值表示当前输入中不同路段之间的动态联系程度;然后对S使用softmax函数进行归一化处理;
步骤3.3:模糊关系图卷积网络模块:
采用K阶近似切比雪夫多项式图卷积来获取交通流数据中各路段之间空间依赖:
其中⊙表示哈达玛积;
步骤3.4:时间卷积模块,采用K层门控空洞因果卷积来获取交通流数据中路段的时间依赖:
其中*表示空洞因果卷积操作,Φ表示卷积核参数,tanh是输出的激活函数,sigmiod函数决定了当前输入传递给下一层的比率。K层空洞因果卷积结构如图3所示的,其通过跳过部分序列扩大卷积的感受野,同时使序列变得更加光滑,增强模型鲁棒性。
步骤四:构建三个并行模型:三个并行模型具有相同的网络结构,每个网络结构由时空卷积单元堆叠,最后一层添加全连接层构成;三个并行模型分别处理周趋势数据、天趋势数据和时刻趋势数据,最后将得到的三个预测结果通过参数加权得到最终预测结果;
融合过程如下:
步骤五:训练深度模型优化参数实现短期交通流的实时预测:将训练集作为深度模型的输入,定义模型的损失函数,根据损失函数数值利用反向传播算法不断更新模型参数;待参数模型收敛后利用训练好的模型对测试集进行实时交通流预测;
定义模型的损失函数,训练模型:
步骤六:实时交通流的路况评估:将实时交通流数据输入到上述训练好的模型,得到预测车速和流量,再根据自由行驶速度和最大自由通行流量计算行驶时间指标TTI和道路流量指标RFI,按照两个指标对交通状态的隶属度函数以及专家经验进行交通状况进行模糊综合评价;
步骤6.1:)TTI,RFI指标两者计算公式如下:
其中Vfree自由通行速度,Vcurrent为当前速度;
其中MaxFfree表示最大自由通行流量,Fcurrent表示当前流量;
TTI指标的理论取值范围为[1,+∞],由于自由通行速度为平均估计数据以及考虑到超速情况,对于小于1的值按向上取整处理,RFI指标的理论取值范围为[1,+∞],将两者归一化按下式归一化处理:
步骤6.2:根据模糊理论,得到因素论域U={TTI,RFI},本发明将交通等级论域V={完全畅通,基本畅通,轻微拥堵,中度拥堵,严重拥堵},采用TTI和RFI的隶属度函数如图3图4所示,高斯隶属的函数的公式为:
在TTI中的三个高斯隶属度函数的k取值依次为[0.3,0.5,0.7],
r的取值依次为[0.07,0.08,0.08];
在RFI中的三个高斯隶属度函数的k取值依次为[0.25,0.4,0.55],
r的取值依次为[0.07,0.08,0.09];
根据隶属度函数计算得到的TTI和RFI指标得到模糊关系矩阵R:
rij表示模糊关系矩阵R中第i个指标对第j个交通等级的隶属度,随后根据专家经验确定评判因素权向量A={a1,a2},ai是U中第i个因素对交通等级评判的重要性分配权重,设置为a1=0.5,a2=0.5;选择合适的评论合成算子,将A与R合成得到综合评判B=(b1,b2,b3,b4,b5),本发明采用算子,其中运算符表示求和取与1比较的较小值:
其中B为模糊综合评价论域,bi为评价隶属值。
一种基于模糊关系的短时交通预测图卷积网络在深圳北站的实施例,包括以下步骤:
步骤1:选取实验数据实验选用深圳北站周边的12条路段,同名公路的不同方向视为不同的公路。首先根据12条路段的地理距离信息构建有向加权图,抽象出道路的空间拓扑关系。然后构建交通流数据集,本实验需要预测的数据是道路的TTI指标,道路的TTI指标数据由深圳市政府网站提供,在该数据集上无需估算每条道路TTI。其他交通数据特征包括平均车速、车流量。根据每条道路的TTI指标和车速可以反推出每条道路的自由通行速度Vfree,根据Vfree估算每条道路的最大自由通行流量,得到结果如下表所示:
路段之间最大自由通行流量存在比较大的差异是因为部分路为快速路,部分路为深圳北站周边人流量多的单行道。数据的采样间隔为10分钟,采样长度为24624。数据按步骤一中的说明进行预处理后按3:1:1的比例划分训练集、验证集和测试集。
步骤2:实验环境及超参数设定
本实验使用MXNet深度学习框架完成模型搭建,并在配置Nvidia Tesla V100GPUs的云服务器上完成模型的训练。小时趋势数据采用目标区间前两个小时数据,目标预测未来一个小时数据。周趋势数据采用前一周目标区间同一时间段的前两个小时数据,天趋势数据采用前一天与目标区间同一时间段的前两个小时数据,。模型学习率设置为0.001。在图卷积网络模块,使用3阶切比雪夫多项式64卷积核的图卷积网络。在时间卷积模块,采用两层门控空洞因果卷积,空洞数分别为0和1,也就是说模型考虑的短时间跨度为20分钟。
本实验选取绝对值均方差(MAE)、均方根误差(RMSE)作为评价指标,其计算公式分别如下所示:
模型预测结果如下表所示:
假设现有对路段1十分钟后的交通预测值为:TTI:2.26和车流量162:由预测值、路段1的最大自由通行车流量和公式(12)~(13),可得到归一化后的TTI和RFI为0.56,0.52带入隶属度函数得到关系矩阵:
根据公式(17)、(18)计算得到B={0,0,0.535,0.578,0};由此可以判断路段1未来十分钟的路况为中度拥堵。
综上所述,针对预测交通流数据转化为有效认知的问题,本发明在缺少交通网络先验知识的情况下,基于历史数据估计了两种对交通流状况的评测指标,并对不同评测指标进行了模糊综合评价给模型的预测结果提供了直观的认识。
Claims (2)
1.短期交通预测图卷积网络的路况评估方法,其特征在于包括以下步骤:
步骤一:获取交通数据并进行预处理:获取交通节点的空间距离信息,用于构建交通网络的图结构;获取的交通流数据建立交通流数据集,提取每个节点的自由行驶速度和最大自由通行流量,然后将历史交通数据按交通节点数N、特征数量F、时间序列T三个维度整理成N×F×T的多维数据集,对缺失数据进行缺失补齐,并进行归一化处理最后划分训练集和测试集;
步骤1.2:提取每个节点的自由行驶速度Vfree和最大自由通行流量MaxFfree;自由行驶速度对历史数据均匀取样后按各节点凌晨时段的历史平均速度数据去噪后平均估计得到,最大自由通行流量按各节点在自由行驶速度向上向下取整范围内的速度对应的流量数据集合取部分最大值平均估计得到;
缺失数据补齐规则:在连续缺失值少于等于3个时采用线性插值,大于3个时采用上个数值、昨天数值与上周数值加权平均,数据补齐后按照采用0均值归一化方法对数据集进行归一化处理,归一化公式如下:
x=(x-mean)/std (3)
其中mean表示该交通特征数据的平均值,std表示该交通特征数据的标准差;
步骤二:对多维数据集划分为小时趋势数据、天趋势数据、周趋势数据:三类数据具有相同的目标值,小时趋势的输入值为目标值时段前几个小时时间段的数据,天趋势的输入值为前几天相同目标时间段与小时趋势同长度的前几个小时的数据,周趋势的输入值为近几周前相同目标时间段的与小时趋势同长度的前几个小时的数据;
步骤2.1周趋势数据、天趋势数据、小时趋势数据从步骤一中的数据集产生;
假设预测时间长度为p,则目标值为[Xt,Xt+1,...,Xt+p],生成规则如下:
小时趋势数据:
Xh=[Xt-n,Xt-n+1,...,Xt-1]
其中n表示取与目标值相邻的前n个时间步;
天趋势数据:
Xd=[X(t-n)-dt*dn,X(t-n+1)-dt*dn,...,X(t-1)-dt*dn;X(t-n)-dt*(dn-1),X(t-n+1)-dt*(dn-1),...,X(t-1)-dt*(dn-1);......;X(t-n)-dt,X(t-n+1)-dt,...,X(t-1)-dt]
其中dt为一天的时间长度,dn表示取与目标值相邻的前dn天;
周趋势数据:
Xw=[X(t-n)-wt*wn,X(t-n+1)-wt*wn,...,X(t-1)-wt*wn;X(t-n)-wt*(wn-1),X(t-n+1)-wt*(wn-1),...,X(t-1)-wt*(wn-1);......;X(t-n)-wt,X(t-n+1)-wt,...,X(t-1)-wt]
其中wt为一周的时间长度,wn表示取与目标值相邻的前wn个周;
步骤三:构建基于模糊关系的短期交通预测图卷积网络:构造时空卷积单元,每个单元包括时间注意力模块、空间注意力模块、模糊关系图卷积模块、门控空洞因果时间卷积模块;通过时空卷积单元的堆叠,完成对交通流张量的时空特征提取,最后构建全连接层通过relu非线性回归得到预测结果;
步骤3.1:时间注意力模块用于捕获输入数据之间的动态时间相关性,即找到当前输入X中相关性较高的时间步:
P=Vp·tanh((XTU1)U2(U3X)+bp) (4)
P中各元素的值表示输入时间序列中不同时间步长之间的动态联系程度;然后对P使用softmax函数进行归一化处理,然后用P乘以输入X得到根据时间相关性动态调整后的输入X;
步骤3.2:空间注意力模块:
S=Vs•tanh((X W1)W2(W3XT)+bs) (5)
其中为空间注意力矩阵,N为交通节点数;Vs,W1,W2,W3,bS为可学习的参数,采用tanh激活函数;S中各元素的值表示当前输入中不同交通节点之间的动态联系程度;然后对S使用softmax函数进行归一化处理;
采用K阶近似切比雪夫多项式图卷积来获取交通流数据中各节点之间空间依赖:
其中⊙表示哈达玛积;
步骤3.4时间卷积模块:
采用K层门控空洞因果卷积来获取交通流数据中交通节点的时间依赖;
其中*表示空洞因果卷积操作,Φ表示卷积核参数,tanh是输出的激活函数,sigmiod函数决定了当前输入传递给下一层的比率;K层空洞因果卷积结构通过跳过部分序列扩大卷积的感受野,同时使序列变得更加光滑,增强模型鲁棒性;
步骤四:构建三个并行模型:三个并行模型具有相同的网络结构,每个网络结构由时空卷积单元堆叠,最后一层添加全连接层构成;三个并行模型分别处理周趋势数据、天趋势数据和时刻趋势数据,最后将得到的三个预测结果通过参数加权得到最终预测结果;融合过程如下:
步骤五:训练深度模型优化参数实现短期交通流的实时预测:将训练集作为深度模型的输入,定义模型的损失函数,根据损失函数数值利用反向传播算法不断更新模型参数;待参数模型收敛后利用训练好的模型对测试集进行实时交通流预测;
定义模型的损失函数,训练模型:
步骤六:实时交通流的路况评估:将实时交通流数据输入到上述训练好的模型,得到预测车速和流量,再根据自由行驶速度和最大自由通行流量计算行驶时间指标TTI和道路流量指标RFI,按照两个指标对交通状态的隶属度函数以及专家经验进行交通状况进行模糊综合评价;
步骤6.1:TTI,RFI指标两者计算公式如下:
其中Vfree自由通行速度,Vcurrent为当前速度;
其中MaxFfree表示最大自由通行流量,Fcurrent表示当前流量;
TTI指标的理论取值范围为[1,+∞],由于自由通行速度为平均估计数据以及考虑到超速情况,对于小于1的值按向上取整处理,RFI指标的理论取值范围为[1,+∞],将两者归一化按下式归一化处理:
步骤6.2:根据模糊理论,得到因素论域U={TTI,RFI},将交通等级论域:
V={完全畅通,基本畅通,轻微拥堵,中度拥堵,严重拥堵},采用TTI和RFI的隶属度函数,高斯隶属的函数的公式为:
在TTI中的三个高斯隶属度函数的k取值依次为[0.3,0.5,0.7],r的取值依次为[0.07,0.08,0.08];
在RFI中的三个高斯隶属度函数的k取值依次为[0.25,0.4,0.55],r的取值依次为[0.07,0.08,0.09];根据隶属度函数计算得到的TTI和RFI指标得到模糊关系矩阵R:
rij表示模糊关系矩阵R中第i个指标对第j个交通等级的隶属度,随后根据专家经验确定评判因素权向量A={a1,a2},ai是U中第i个因素对交通等级评判的重要性分配权重,设置为a1=0.5,a2=0.5;选择合适的评论合成算子,将A与R合成得到综合评判B=(b1,b2,b3,b4,b5),采用算子,其中运算符表示求和取与1比较的较小值:
其中B为模糊综合评价论域,bi为评价隶属值。
2.根据权利要求1所述的短期交通预测图卷积网络的路况评估方法,其特征在于:所述步骤一中,F表示收集的交通流数据特征总数,至少包括车流量、平均速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110603374.2A CN113313947B (zh) | 2021-05-31 | 2021-05-31 | 短期交通预测图卷积网络的路况评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110603374.2A CN113313947B (zh) | 2021-05-31 | 2021-05-31 | 短期交通预测图卷积网络的路况评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113313947A CN113313947A (zh) | 2021-08-27 |
CN113313947B true CN113313947B (zh) | 2022-04-19 |
Family
ID=77376624
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110603374.2A Active CN113313947B (zh) | 2021-05-31 | 2021-05-31 | 短期交通预测图卷积网络的路况评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113313947B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114267170A (zh) * | 2021-09-28 | 2022-04-01 | 浙江工业大学 | 基于考虑人类流动性的图时空Transformer模型的交通流预测方法 |
CN114120630B (zh) * | 2021-10-15 | 2023-04-07 | 国网浙江省电力有限公司杭州供电公司 | 一种考虑道路交通动态特征的应急电源配置方法 |
CN114038200B (zh) * | 2021-11-29 | 2022-09-20 | 东北大学 | 基于注意力机制的时空同步图卷积网络交通流量预测方法 |
CN114299728B (zh) * | 2021-12-29 | 2023-04-21 | 杭州电子科技大学 | 结合注意力机制与动态时空卷积模型的车辆流量预测方法 |
CN114743374B (zh) * | 2022-04-01 | 2023-04-07 | 清华大学深圳国际研究生院 | 基于图网络的多元交通流预测方法 |
CN114973653B (zh) * | 2022-04-27 | 2023-12-19 | 中国计量大学 | 一种基于时空图卷积网络的交通流预测方法 |
CN115482656B (zh) * | 2022-05-23 | 2023-09-26 | 汕头大学 | 一种使用空间动态图卷积网络进行交通流预测方法 |
CN114913689B (zh) * | 2022-06-01 | 2023-06-16 | 南通大学 | 一种基于先验知识时变图卷积网络的交通流预测方法 |
CN115512545B (zh) * | 2022-09-30 | 2023-06-27 | 浙江财经大学 | 一种基于时空动态图卷积网络的交通速度预测方法 |
CN117198037B (zh) * | 2023-05-07 | 2024-02-09 | 河海大学 | 一种基于图卷积回声状态网络的交通流预测方法 |
CN116413395B (zh) * | 2023-06-08 | 2023-08-25 | 山东瑞程数据科技有限公司 | 一种环境异常智能检测方法 |
CN116957166B (zh) * | 2023-08-25 | 2024-04-26 | 南京纳尼亚科技有限公司 | 一种基于鸿蒙系统的隧道交通情况预测方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109285346A (zh) * | 2018-09-07 | 2019-01-29 | 北京航空航天大学 | 一种基于关键路段的城市路网交通状态预测方法 |
CN109754605A (zh) * | 2019-02-27 | 2019-05-14 | 中南大学 | 一种基于注意力时态图卷积网络的交通预测方法 |
CN110264709A (zh) * | 2019-05-06 | 2019-09-20 | 北京交通大学 | 基于图卷积网络的道路的交通流量的预测方法 |
CN110969854A (zh) * | 2019-12-13 | 2020-04-07 | 深圳先进技术研究院 | 一种交通流量的预测方法、系统及终端设备 |
-
2021
- 2021-05-31 CN CN202110603374.2A patent/CN113313947B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109285346A (zh) * | 2018-09-07 | 2019-01-29 | 北京航空航天大学 | 一种基于关键路段的城市路网交通状态预测方法 |
CN109754605A (zh) * | 2019-02-27 | 2019-05-14 | 中南大学 | 一种基于注意力时态图卷积网络的交通预测方法 |
CN110264709A (zh) * | 2019-05-06 | 2019-09-20 | 北京交通大学 | 基于图卷积网络的道路的交通流量的预测方法 |
CN110969854A (zh) * | 2019-12-13 | 2020-04-07 | 深圳先进技术研究院 | 一种交通流量的预测方法、系统及终端设备 |
Non-Patent Citations (2)
Title |
---|
基于云模型的短时交通流预测方法研究;刘芹等;《计算机工程与设计》;20120516(第05期);第1953-1957页 * |
基于时空多图卷积网络的交通站点流量预测;荣斌等;《计算机工程》;20200531(第05期);第26-33页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113313947A (zh) | 2021-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113313947B (zh) | 短期交通预测图卷积网络的路况评估方法 | |
CN111223301B (zh) | 一种基于图注意力卷积网络的交通流量预测方法 | |
CN111612206B (zh) | 一种基于时空图卷积神经网络的街区人流预测方法及系统 | |
CN114299723B (zh) | 一种交通流量预测方法 | |
CN110570035B (zh) | 同时建模时空依赖性和每日流量相关性的人流量预测系统 | |
CN112071062B (zh) | 一种基于图卷积网络和图注意力网络的行车时间估计方法 | |
CN114802296A (zh) | 一种基于动态交互图卷积的车辆轨迹预测方法 | |
CN114202122A (zh) | 基于马尔可夫聚类图注意力网络的城市交通流量预测方法 | |
CN113762595B (zh) | 通行时间预测模型训练方法、通行时间预测方法及设备 | |
CN114495507B (zh) | 融合时空注意力神经网络和交通模型的交通流预测方法 | |
CN113537580B (zh) | 一种基于自适应图学习的公共交通客流预测方法及系统 | |
CN115376317B (zh) | 一种基于动态图卷积和时序卷积网络的交通流预测方法 | |
CN116307152A (zh) | 时空交互式动态图注意力网络的交通预测方法 | |
CN112766603A (zh) | 一种交通流量预测方法、系统、计算机设备及存储介质 | |
CN115206092A (zh) | 一种基于注意力机制的BiLSTM和LightGBM模型的交通预测方法 | |
CN114566047A (zh) | 一种基于短时路况预测算法的预警方法及系统 | |
CN112562312B (zh) | 一种基于融合特征的GraphSAGE交通路网数据预测的方法 | |
CN117456730A (zh) | 一种基于自适应部分注意力扩散图卷积神经网络的交通流量预测方法 | |
CN116578858A (zh) | 基于图神经网络的空压机故障预测与健康度评价方法及系统 | |
CN113240904B (zh) | 基于特征融合的交通流预测方法 | |
CN114254214A (zh) | 一种基于时空层次化网络的交通预测方法及系统 | |
CN113516163B (zh) | 基于网络剪枝的车辆分类模型压缩方法、装置及存储介质 | |
CN114267170A (zh) | 基于考虑人类流动性的图时空Transformer模型的交通流预测方法 | |
CN114004421B (zh) | 一种基于时空集成学习的交通数据缺失值插补方法 | |
CN115905434B (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 |