基于兴趣点时空残差神经网络的多方向交通流量预测方法

1.本发明涉及交通流预测领域，尤其涉及一种基于兴趣点时空残差神经网络的多方向交通流量预测方法。


背景技术：

2.交通数据作为一种普遍存在的时空数据类型，同时在时间和空间维度上表现出相关性。此外，过去几十年的交通预测研究主要集中在出租车、私家车和共享单车上，而对电动自行车的研究较少。另外，又由于电动自行车行驶时间短、行驶距离短等原因，导致电动自行车行驶轨迹数据周期性不明显，多种原因使得电动自行车的预测比传统的交通预测更具挑战性。近年来，电动自行车由于价格低廉、体积小，逐渐成为人们日常生活中不可缺少的交通工具。据统计，我国电动自行车保有量近3亿辆，几乎家家户户都至少拥有一辆电动自行车。而大量的电动自行车会导致许多交通问题，如非机动车车道拥堵等。因此，电动自行车交通流预测已成为解决部分交通问题的一项有价值的课题。此外，电动自行车流量的实时预测还可以为各种交通服务提供大量的数据支持：1)提高电动自行车的导航精度；2)移动充电桩的实施；3)为乘客提供更详细的交通信息，以便乘客制定更好的出行计划。
3.综上所述，准确预测电动自行车的流量非常具有挑战性，主要受到以下4 个因素的影响：
4.1、空间相关性：由于空间的局部相关性，在相邻位置获得的观测值相互关联。
5.2、时间相关性：在连续时间间隔获得的观测值高度相关，即，如果从最近几个连续时间间隔内的观测值中检测到上升趋势，则下一时刻的观测值很有可能会增加。
6.3、地理环境：地理环境对电动自行车流量预测有显著影响，即购物中心和公司集群在早高峰时间段内具有更大的吸力，而住宅集群在早高峰时间段内具有更大的斥力。
7.4、异质性：交通数据在空间和时间两个维度上都表现出异质性，即两个地点之间的相互影响随距离和地理环境而变化。同样，时间影响的强度也因地区而异。


技术实现要素：

8.为了克服上述现有技术的不足，本发明提供一种基于兴趣点时空残差神经网络的电瓶车流量预测方法，可有效解决上述问题。本发明具体采用的技术方案如下：
9.一种基于兴趣点时空残差神经网络的多方向交通流量预测方法，其包括以下步骤：
10.s1、获取待预测区域中待预测时刻前的历史车流数据，所述历史车流数据中包含待预测区域中不同车辆在不同时间的位置以及车辆行驶方向；从所述历史车流数据中按预设的时间片间隔提取若干个车流数据时间片；
11.s2、对待预测区域进行栅格化处理从而划分为一系列栅格，每一个车流数据时间片中的车辆均按照其坐标映射至待预测区域的对应栅格中，并将车辆行驶方向定义为车辆在栅格中的移动状态，移动状态包括向上、向下、向左、向右四种；对每个时间片内每个栅格
中所包含的每一种移动状态的车辆总数进行统计并作为栅格值，将栅格值映射为矩阵元素，从而针对每个时间片中的不同移动状态分别构建车流二维矩阵，每一种移动状态在所有时间片中的车流二维矩阵按时间维度叠加形成车流三维矩阵；
12.s3、从待预测时刻中提取小时字段和分钟字段，拼接形成一个二元的时间信号向量；
13.s4、获取所有兴趣点(poi)空间地理位置，并将不同功能类别的兴趣点均映射到所述待预测区域的栅格中，对每个栅格中每一组功能类别的兴趣点总数进行统计并作为栅格值，将栅格值映射为矩阵元素，从而针对每一组功能类别的兴趣点分别构建一张二维矩阵形式的兴趣点切片，所有功能类别的兴趣点切片叠加形成三维张量形式的兴趣点信号；
14.s5、将所述车流三维矩阵、所述时间信号向量和所述兴趣点信号作为经过训练的多方向交通流量预测模型的输入，所述多方向交通流量预测模型包含第一全连接神经网络、第二全连接神经网络、三维残差卷积网络和重校正层，由第一全连接神经网络根据所述时间信号向量输出时间信号矩阵，由第二全连接神经网络根据所述兴趣点信号输出兴趣点信号矩阵，由三维残差卷积网络根据所述车流三维矩阵、所述兴趣点信号矩阵和所述时间信号矩阵的融合特征输出结果矩阵，最终将结果矩阵在重校正层中经过带权压缩操作后，得到多方向交通流量预测结果。
15.作为优选，在s1中，所述待预测区域为矩形区域；所述历史车流数据的时间跨度为[t,t+(m-1)*τ]，且其按τ分钟为时间片间隔提取m个车流数据时间片。
[0016]
进一步的，所述s2的具体实现步骤如下：
[0017]
s21、对待预测区域进行栅格化处理，共划分为i*j个栅格，第i行第j列的栅格为p
ij
；
[0018]
s22、每一个车流数据时间片中的车辆均按照其坐标映射至待预测区域的对应栅格中，并将车辆行驶方向定义为车辆在栅格中的移动状态，移动状态包括向上、向下、向左、向右四种；
[0019]
s23、对每个时间片内每个栅格中所包含的每一种移动状态的车辆总数进行统计，将时间片t内栅格p
ij
内移动状态为d的车流量记为将所有i*j个栅格对应的构建为时间片t内整个待预测区域中移动状态为d的车流二维矩阵车流二维矩阵中第i行第j列的元素值为
[0020]
s24、将所有m个车流数据时间片的车流二维矩阵按时间维度进行拼接，形成车流三维矩阵
[0021]
进一步的，所述s3的具体实现步骤如下：
[0022]
将待预测时刻t
pred
转化为含有小时字段t
pred_hour
和分钟字段t
pred_minute
两个元素的时间信号向量h
t
＝[t
pred_hour
,t
pred_minute
]。
[0023]
作为优选，所述s4的具体实现步骤如下：
[0024]
s41、获取所有兴趣点的空间地理位置，并将所有兴趣点按照其位置均映射到所述待预测区域的各p
ij
栅格中；
[0025]
s42、将所有兴趣点按照不同功能类别分成n组，统计每个兴趣点分组g在栅格p
ij
内
的兴趣点数量并作为栅格p
ij
的栅格值将每个兴趣点分组g对应的所有栅格的栅格值构建为兴趣点分组g对应的兴趣点切片υg，υg大小为i*j；
[0026]
s43、将所有n组兴趣点分组对应的兴趣点切片拼接，得到兴趣点信号ψ＝[υ1,υ2,
…
,υn]，其大小为n*i*j。
[0027]
作为优选，所述s5的具体实现步骤如下：
[0028]
s51：将所述车流三维矩阵所述时间信号向量h
t
和所述兴趣点信号ψ输入经过训练的多方向交通流量预测模型中，所述多方向交通流量预测模型中包含第一全连接神经网络、第二全连接神经网络、三维残差卷积网络和重校正层；
[0029]
s52、所述时间信号h
t
被输入至包含l
ts
层全连接层级联的第一全连接神经网络中，其中第1层全连接层的输入为时间信号h
t
，下一层全连接层的输入为上一层全连接层的输出，最后一层全连接层的输出为将向量逐元素映射到大小为(i*j)的一个矩阵当中，得到大小为(i,j)的时间信号矩阵h
t
；
[0030]
s53、所述兴趣点信号ψ＝[υ1,υ2,
…
,υn]输入后，先对其中的每个兴趣点切片υg获取其平均自权重zg：
[0031][0032]
得到兴趣点信号ψ的平均自权重矩阵z＝{z1,z2,
…
,zn}，其中n表示兴趣点分组的个数；
[0033]
然后将平均自权重矩阵z输入包含l
ps
层全连接层级联的第二全连接神经网络逐层计算，其中第1层全连接层的输入为平均自权重矩阵z，下一层全连接层的输入为上一层全连接层的输出，最后一层全连接层的输出为
[0034]
再后，采用门机制将输出映射为0到1之间的变量，其计算过程如下：
[0035][0036]
其中f
si
为激活函数；
[0037]
最后，通过计算获取兴趣点信号矩阵
[0038][0039]
其中，
⊙
表示矩阵点乘；
[0040]
s54、将车流三维矩阵兴趣点信号矩阵和时间信号矩阵h
t
进行特征融合，中的第k个车流数据时间片的融合特征计算公式如下：
[0041][0042]
其中，其中，和均为可训练的参数，m为中车流数据时间片的个数；
[0043]
最终得到融合后的车流矩阵
[0044]
s55、将融合后的车流矩阵x
γ
输入包含lc层三维残差卷积层级联的三维残差卷积网络中逐层计算，其中第一层三维残差卷积层的输入车流矩阵x
γ
，每一层三维残差卷积层经过三维残差卷积操作后得到的结果作为下一层三维残差卷积层的输入，三维残差卷积网络中所有三维残差卷积层的输出拼接形成最终的结果矩阵x
st
；其中对于任意第l层三维残差卷积层，其中执行的三维残差卷积操作如下：
[0045]
首先，对当前三维残差卷积层的输入进行三维卷积操作得到卷积结果：
[0046][0047]
其中，cov3d表示三维卷积操作，表示第l-1层三维残差卷积层的输出，其中1层三维残差卷积层的输出，其中和为第l层三维卷积操作中可训练的参数，fc为激活函数；
[0048]
然后，对三维卷积操作输出的卷积结果中每个元素进行批正则化操作，得到批正则结果
[0049]
最后，将批正则结果再与上一层的输出矩阵相加，得到第l层三维残差卷积层的输出矩阵公式为：
[0050][0051]
s56、在重校正层中对最终输出的结果矩阵x
st
中所有维进行带权压缩操作，得到预测结果其计算公式如下：
[0052][0053]
其中是可学习的参数矩阵。
[0054]
作为优选，所述多方向交通流量预测模型预先通过训练数据进行训练，训练过程中通过损失函数不断迭代预测结果与真实结果φ之间的损失值，当达到迭代终止条件时，输出多方向交通流量预测模型用于进行实际的预测。
[0055]
作为优选，所述的多方向交通流量预测模型的损失函数loss公式为：
[0056][0057]
其中为矩阵中各个元素的值，为矩阵φ中各个元素的值，m为训练样本数。
[0058]
作为优选，所述的m取5，时间片间隔τ取4。
[0059]
作为优选，所述的车辆在栅格中的移动状态根据车辆行驶轨迹经过该栅格时的车辆行驶方向确定。
[0060]
本发明提出的预测模型从空间和时间两个视图对车流变化信息进行了提取，同时加入了时间信号和兴趣点信号提升了预测的准确度。相比于传统的交通流量预测方法，本
法明具有如下收益：1、采用四种新的移动状态来描述电瓶车流的变化趋势；2、融入了时间信号来加强数据的时间特征；3、使用了兴趣点注意力机制，可以自适应地捕获不同的兴趣点对车流变化产生的影响，进一步加强了数据的空间上的地理特征；4、使用残差网络解决了梯度消失问题。
附图说明
[0061]
图1为基于兴趣点时空残差神经网络的多方向交通流量预测方法的流程图；
[0062]
图2为本发明使用的车流量栅格化计算方法示意图。
具体实施方式
[0063]
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
[0064]
相反，本发明涵盖任何由权利要求定义的在本发明的精髓和范围上做的替代、修改、等效方法以及方案。进一步，为了使公众对本发明有更好的了解，在下文对本发明的细节描述中，详尽描述了一些特定的细节部分。对本领域技术人员来说没有这些细节部分的描述也可以完全理解本发明。
[0065]
作为本发明的一种较佳实现形式，提供了一种基于兴趣点时空残差神经网络的多方向交通流量预测方法，包括以下步骤：
[0066]
s1、获取待预测区域中待预测时刻前的历史车流数据，所述历史车流数据中包含待预测区域中不同车辆在不同时间的位置以及车辆行驶方向；从所述历史车流数据中按预设的时间片间隔提取若干个车流数据时间片。
[0067]
在该较佳实现形式中，在s1中，所述待预测区域为矩形区域；所述历史车流数据的时间跨度为[t,t+(m-1)*τ]，且其按τ分钟为时间片间隔提取m个车流数据时间片。
[0068]
s2、对待预测区域进行栅格化处理从而划分为一系列栅格，每一个车流数据时间片中的车辆均按照其坐标映射至待预测区域的对应栅格中，并将车辆行驶方向定义为车辆在栅格中的移动状态，移动状态包括向上、向下、向左、向右四种；对每个时间片内每个栅格中所包含的每一种移动状态的车辆总数进行统计并作为栅格值，将栅格值映射为矩阵元素，从而针对每个时间片中的不同移动状态分别构建车流二维矩阵，每一种移动状态在所有时间片中的车流二维矩阵按时间维度叠加形成车流三维矩阵。其中，车辆在栅格中的移动状态可以根据车辆行驶轨迹经过该栅格时的车辆行驶方向确定。
[0069]
在该较佳实现形式中，所述s2的具体实现步骤如下：
[0070]
s21、对待预测区域进行栅格化处理，共划分为i*j个栅格，第i行第j列的栅格为p
ij
；
[0071]
s22、每一个车流数据时间片中的车辆均按照其坐标映射至待预测区域的对应栅格中，并将车辆行驶方向定义为车辆在栅格中的移动状态，移动状态包括向上、向下、向左、向右四种；
[0072]
s23、对每个时间片内每个栅格中所包含的每一种移动状态的车辆总数进行统计，将时间片t内栅格p
ij
内移动状态为d的车流量记为将所有i*j个栅格对应的构建
为时间片t内整个待预测区域中移动状态为d的车流二维矩阵车流二维矩阵中第i行第j列的元素值为
[0073]
s24、将所有m个车流数据时间片的车流二维矩阵按时间维度进行拼接，形成车流三维矩阵
[0074]
s3、从待预测时刻中提取小时字段和分钟字段，拼接形成一个二元的时间信号向量。
[0075]
在该较佳实现形式中，所述s3的具体实现步骤如下：
[0076]
将待预测时刻t
pred
转化为含有小时字段t
pred_hour
和分钟字段t
pred_minute
两个元素的时间信号向量h
t
＝[t
pred_hour
,t
pred_minute
]。
[0077]
s4、获取所有兴趣点空间地理位置，并将不同功能类别的兴趣点均映射到所述待预测区域的栅格中，对每个栅格中每一组功能类别的兴趣点总数进行统计并作为栅格值，将栅格值映射为矩阵元素，从而针对每一组功能类别的兴趣点分别构建一张二维矩阵形式的兴趣点切片，所有功能类别的兴趣点切片叠加形成三维张量形式的兴趣点信号。
[0078]
在该较佳实现形式中，所述s4的具体实现步骤如下：
[0079]
s41、获取所有兴趣点的空间地理位置，并将所有兴趣点按照其位置均映射到所述待预测区域的各p
ij
栅格中；
[0080]
s42、将所有兴趣点按照不同功能类别分成n组，统计每个兴趣点分组g在栅格p
ij
内的兴趣点数量并作为栅格p
ij
的栅格值将每个兴趣点分组g对应的所有栅格的栅格值构建为兴趣点分组g对应的兴趣点切片υg，υg大小为i*j；
[0081]
s43、将所有n组兴趣点分组对应的兴趣点切片拼接，得到兴趣点信号ψ＝[υ1,υ2,
…
,υn]，其大小为n*i*j。
[0082]
s5、将所述车流三维矩阵、所述时间信号向量和所述兴趣点信号作为经过训练的多方向交通流量预测模型的输入，所述多方向交通流量预测模型包含第一全连接神经网络、第二全连接神经网络、三维残差卷积网络和重校正层，由第一全连接神经网络根据所述时间信号向量输出时间信号矩阵，由第二全连接神经网络根据所述兴趣点信号输出兴趣点信号矩阵，由三维残差卷积网络根据所述车流三维矩阵、所述兴趣点信号矩阵和所述时间信号矩阵的融合特征输出结果矩阵，最终将结果矩阵在重校正层中经过带权压缩操作后，得到多方向交通流量预测结果。
[0083]
在该较佳实现形式中，所述s5的具体实现步骤如下：
[0084]
s51：将所述车流三维矩阵所述时间信号向量h
t
和所述兴趣点信号ψ输入经过训练的多方向交通流量预测模型中，所述多方向交通流量预测模型中包含第一全连接神经网络、第二全连接神经网络、三维残差卷积网络和重校正层；
[0085]
s52、所述时间信号h
t
被输入至包含l
ts
层全连接层级联的第一全连接神经网络中，其中第1层全连接层的输入为时间信号h
t
，下一层全连接层的输入为上一层全连接层的输出，最后一层全连接层的输出为将向量逐元素映射到大小为(i*j)的一个矩阵当中，得到大小为(i,j)的时间信号矩阵h
t
；
[0086]
s53、所述兴趣点信号ψ＝[υ1,υ2,
…
,υn]输入后，先对其中的每个兴趣点切片
υg获取其平均自权重zg：
[0087][0088]
得到兴趣点信号ψ的平均自权重矩阵z＝{z1,z2,
…
,zn}，其中n表示兴趣点分组的个数；
[0089]
然后将平均自权重矩阵z输入包含l
ps
层全连接层级联的第二全连接神经网络逐层计算，其中第1层全连接层的输入为平均自权重矩阵z，下一层全连接层的输入为上一层全连接层的输出，最后一层全连接层的输出为
[0090]
再后，采用门机制将输出映射为0到1之间的变量，其计算过程如下：
[0091][0092]
其中f
si
为激活函数；
[0093]
最后，通过计算获取兴趣点信号矩阵
[0094][0095]
其中，
⊙
表示矩阵点乘；
[0096]
s54、将车流三维矩阵兴趣点信号矩阵和时间信号矩阵h
t
进行特征融合，中的第k个车流数据时间片的融合特征计算公式如下：
[0097][0098]
其中，其中，和均为可训练的参数，m为中车流数据时间片的个数；
[0099]
最终得到融合后的车流矩阵
[0100]
s55、将融合后的车流矩阵x
γ
输入包含lc层三维残差卷积层级联的三维残差卷积网络中逐层计算，其中第一层三维残差卷积层的输入车流矩阵x
γ
，每一层三维残差卷积层经过三维残差卷积操作后得到的结果作为下一层三维残差卷积层的输入，三维残差卷积网络中所有三维残差卷积层的输出拼接形成最终的结果矩阵x
st
；其中对于任意第l层三维残差卷积层，其中执行的三维残差卷积操作如下：
[0101]
首先，对当前三维残差卷积层的输入进行三维卷积操作得到卷积结果：
[0102][0103]
其中，cov3d表示三维卷积操作，表示第l-1层三维残差卷积层的输出，其中1层三维残差卷积层的输出，其中和为第l层三维卷积层(即前述三维卷积操作 cov3d)可训练的参数，fc为激活函数；
[0104]
然后，对三维卷积层输出的卷积结果中每个元素进行批正则化操作，得到批
正则结果其公式如下：
[0105][0106]
其中，e[x]表示每维矩阵的均值，var[x]为每维矩阵的方差，∈是为防止方差为0设置的常数，υ与β为可学习参数；
[0107]
最后，将批正则结果再与上一层的输出矩阵相加，得到第l层三维残差卷积层的输出矩阵公式为：
[0108][0109]
s56、在重校正层中对最终输出的结果矩阵x
st
中所有维进行带权压缩操作，得到预测结果其计算公式如下：
[0110][0111]
其中是可学习的参数矩阵。
[0112]
需要注意的是，上述s5中，多方向交通流量预测模型预先通过训练数据进行训练，训练过程中通过损失函数不断迭代预测结果与真实结果φ之间的损失值，当达到迭代终止条件时，输出多方向交通流量预测模型用于进行实际的预测。
[0113]
作为一种多方向交通流量预测模型的损失函数loss的实现形式，其公式为：
[0114][0115]
其中为矩阵中各个元素的值，为矩阵φ中各个元素的值，m为训练样本数。
[0116]
下面将上述s1～s5所描述的基于兴趣点时空残差神经网络的多方向交通流量预测方法应用于一个具体的实例中，通过实施例进一步展示其技术效果。
[0117]
实施例
[0118]
本实施例所用的原始数据为郑州的真实电瓶车轨迹数据集。通过预处理，每条电瓶车轨迹记录保留了车辆编号、经度、纬度、方向和时间这5个字段。至此形成许多条数据σ＝《车辆编号，经度，维度，方向，时间》，其时间单位为分钟。将每条σ按照时间，映射至每个时间的地图中。将地图按照每方格 400米*400米大小进行栅格化。统计每个栅格相同行驶状态下所有车辆总和，再统计每个栅格相同分组下的兴趣点总和，分别生成车流矩阵和兴趣点信号，再将待预测时刻处理成时间信号。最后将车流矩阵、兴趣点信号和时间信号按照上述的方法输入至模型中进行训练，使用训练好的模型进行结果预测。如图 1所示，本实施例中的基于兴趣点时空残差神经网络的多方向交通流量预测方法具体包括以下步骤：
[0119]
步骤1)定义待预测区域：将需要进行预测的待预测区域定义为边长为i*j 的矩形区域；
[0120]
步骤2)定义车流移动状态：对待预测区域内的所有移动车辆区分状态，每辆车的
移动状态表示为d，d包括向上、向下、向左、向右四种移动状态，即移动状态d取leftward，rightward，upward和downward四种状态；
[0121]
步骤3)对车流进行栅格化处理：将待预测区域自坐标系原点置于二维坐标系中，按照i*j的栅格对待预测区域进行分割，栅格在x轴与y轴方向上的序号分别为i与j，将栅格记为p
ij
；
[0122]
步骤4)对每个栅格进行车流赋值，并计算区域车流二维矩阵：将时间范围划分出为多个时长ξ的时间片，对时间片t内每个栅格内包含的所有车辆进行区分移动状态求和，即时间片t内每个栅格p
ij
内移动状态为d的车流量记为时间片t内待预测区域内移动状态为d的车流量二维矩阵记为大小为(i*j)；
[0123]
步骤5)获取区域车流三维矩阵：对预测时刻前m*τ分钟的车流数据按τ分钟间隔依次执行m次步骤2到步骤6步骤，得到m个i*j的时间片，将其记作车流三维矩阵大小为(m*i*j)，如图2所示；
[0124]
步骤6)获取时间信号：将预测时刻t
pred
转化为含有两个元素的时间信号，表示为h
t
＝[t
pred_hour
,t
pred_minute
]；
[0125]
步骤7)获取兴趣点信号，具体包含以下子步骤：
[0126]
步骤71)获取所有兴趣点空间地理位置，将兴趣点映射到平面上，然后同步骤3对平面进行栅格划分，将栅格记为p
ij
；
[0127]
步骤72)将兴趣点按照功能性分成n组，对其中的兴趣点分组g，计算栅格 p
ij
内的兴趣点数量和得到栅格值依次对每个栅格执行求和操作后得到第g 组兴趣点的兴趣点切片υg，其大小为i*j；
[0128]
步骤73)对兴趣点分组内的每个分组迭代执行步骤72)操作，得到n个 i*j的兴趣点切片，将其记作兴趣点信号ψ＝[υ1,υ2,
…
,υn]，其大小为(n*i*j)；
[0129]
步骤8)训练多方向交通流量预测模型，具体包含以下子步骤：
[0130]
步骤81)获取时间信号矩阵，包含以下子步骤：
[0131]
步骤811)首先将时间信号h
t
迭代输入至包含l
ts
个全连接层的第一全连接神经网络，其中下一个全连接层的输入为上一个全连接层的输出，网络逐层计算的公式如下：
[0132][0133]
其中l
ts
是第一全连接神经网络中全连接层的数量，和分别表示第l-1层和第l层的输出，和是第l层的两个可学习参数，f
ts
表示激活函数，全连接层的第一层输入为时间信号h
t
，最后一个全连接层的输出为
[0134]
步骤812)将向量逐元素映射到大小为(i*j)的一个矩阵当中，得到时间信号矩阵h
t
，其大小为(i*j)；该逐元素映射过程实际上相当于的元素重构，将的前j个元素作为矩阵的第1行，前j+1到2j个元素作为矩阵的第2行，以此类推，从而形成i*j大小的时间信号矩阵h
t
。
[0135]
步骤82)获取兴趣点信号矩阵，包含以下子步骤：
[0136]
步骤821)输入兴趣点信号ψ＝[υ1,υ2,
…
,υn]，对于其中的每个兴趣点切片υg，获取其平均自权重zg：
[0137][0138]
对所有兴趣点切片执行上述操作，得到兴趣点信号ψ的平均自权重矩阵 z＝{z1,z2,
…
,zn}，其中n表示兴趣点分组的个数；
[0139]
步骤822)将平均自权重矩阵z输入包含l
ps
个全连接层的第二全连接神经网络，其中下一个全连接层的输入为上一个全连接层的输出，网络逐层计算的公式为：
[0140][0141]
其中l
ps
是第二全连接神经网络中全连接层的数量，和分别表示第l-1层和第l层的输出，和是第l层的两个可学习参数，f
ps
表示激活函数，第一层全连接层的输入为平均自权重矩阵z，最后一个全连接层的输出为
[0142]
步骤823)采用门机制将步骤822)中第二全连接神经网络的输出映射为0到1之间的变量，其计算过程如下：
[0143][0144]
其中f
si
为激活函数；
[0145]
步骤824)获取兴趣点信号矩阵其过程为：
[0146][0147]
其中，
⊙
表示矩阵点乘，υg表示第g组兴趣点的兴趣点切片；
[0148]
步骤83)将车流三维矩阵兴趣点信号矩阵和时间信号矩阵h
t
进行特征融合，对中的每个时间片的融合过程，其中中的第k个车流数据时间片的融合特征计算公式如下：
[0149][0150]
其中，和均为可训练的参数，m为中时间片的个数；
[0151]
最终得到融合后的车流矩阵
[0152]
步骤84)将融合后的车流矩阵x
γ
依次输入包含lc个三维残差卷积层的三维残差卷积网络，最终得到结果矩阵x
st
，其具体过程如下：
[0153]
步骤841)将x
γ
输入第l个三维残差卷积层进行三维卷积操作，其卷积结果的计算公式如下：
[0154]
(l需初始化为1)
[0155]
其中，cov3d表示三维卷积操作，表示第l-1层三维残差卷积层的输出，1层三维残差卷积层的输出，和为第l层可训练的参数，fc为激活函数sigmoid。
[0156]
步骤842)然后对中每个元素进行批正则化操作，其公式如下：
[0157][0158]
其中，e[x]表示每维矩阵的均值，var[x]为每维矩阵的方差，∈是为防止方差为0设置的常数，υ与β为可学习参数；
[0159]
步骤843)最后将经过批正则函数后的结果与上一层的输出结果相加，得到该层三维残差卷积层的输出矩阵公式为：
[0160][0161]
其中和分别表示第l个三维残差卷积层和第l-1个三维残差卷积层的输出经过步骤842)和步骤843)操作后得到的输出矩阵；
[0162]
步骤844)对于剩余l
c-1层三维残差卷积层，不断循环执行l
c-1次步骤 841)到步骤843)，其中下一个三维残差卷积层的输入为上一个三维残差卷积层的输出，最后完成所有三维残差卷积层的特征提取后，将lc层三维残差卷积层的输出矩阵拼接，最终输出的结果矩阵记为x
st
；
[0163]
步骤85)对结果矩阵x
st
进行带权压缩操作，得到预测结果其计算公式如下：
[0164][0165]
其中是可学习的参数矩阵；
[0166]
步骤86)通过损失函数不断迭代预测结果与真实结果φ之间的损失值，当达到迭代终止条件时，输出多方向交通流量预测模型
[0167]
本实施例中模型的损失函数loss公式为：
[0168][0169]
其中为矩阵中各个元素的值，为矩阵φ中各个元素的值，m为训练样本总数。
[0170]
本实施例中采用adam作为优化器，多方向交通流量预测模型的迭代终止条件可以设置为迭代轮次达到10次或者损失函数值小于0.1％。
[0171]
步骤9)预测结果：在推理阶段，与步骤5)、步骤6)和步骤7)操作一致，分别获取待预测时间前m*τ分钟的车流矩阵、时间信号向量和兴趣点信号矩阵，输入至多方向交通流量预测模型后得到车流预测结果矩阵
[0172]
在本实施例中，时长ξ取1分钟。m取5，时间片的间隔τ取4min。第一全连接网络中，全连接层数l
ts
取2，激活函数f
ts
取relu。第二全连接网络中，全连接层数l
ps
取2，激活函数f
ps
取relu。三维残差卷积网络中，三维残差卷积层的层数lc取32，激活函数fc取relu。
[0173]
另外，兴趣点poi是实现城市功能的地理实体，体现了不同出发地和目的地对交通量的变化的影响。例如，餐饮poi在午餐和晚餐时间影响周边地区的交通量，而旅游景点poi主要影响周末和节假日的交通量。在本实施例中，兴趣点分组数n取9，具体的类别形式如表1所示。
[0174]
表1.兴趣点及其类别的详细信息。
[0175][0176]
为了便于对比本发明的效果，本实施例中将本发明方法(简称mdstn) 与若干传统预测方法进行了比较，作为对照的预测方法有：(1)ha：使用对应时间最近的历史数据的平均值作为预测值；(2)arima：差分整合移动平均自回归模型，又称整合移动平均自回归模型；(3)lstm：长短是记忆神经网络，循环神经网络的一个变种。(4)gru：门控循环单元，循环神经网络的一个变种。(5)convlstm：卷积长短时记忆神经网络，一种能够捕获时空相关性的增强型lstm。(6)st-resnet：一种深层神经网络模型，融合了卷积操作和残差单元的网络。(7)st-3dnets：一种改进的深度神经网络的模型，利用三维卷积和残差单元来探索时空相关性。本实施例以均方误差(mse)与平均百分比误差(mape)作为预测准确度指标，具体定义如下：
[0177][0178][0179]
其中为真实值，为预测值，ij表示区域，t表示时间，m为样本的数量。
[0180]
表1预测结果比较
[0181][0182]
实验结果如表所示，结果表明本发明提出的方法预测精度(mse、mape) 明显优于传统数学方法与机器学习相关方法，证明了本发明方法的有效性。
[0183]
以上所述的实施例只是本发明的一种较佳的方案，然其并非用以限制本发明。有关技术领域的普通技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型。因此凡采取等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。

技术特征：
1.一种基于兴趣点时空残差神经网络的多方向交通流量预测方法，其特征在于，包括以下步骤：s1、获取待预测区域中待预测时刻前的历史车流数据，所述历史车流数据中包含待预测区域中不同车辆在不同时间的位置以及车辆行驶方向；从所述历史车流数据中按预设的时间片间隔提取若干个车流数据时间片；s2、对待预测区域进行栅格化处理从而划分为一系列栅格，每一个车流数据时间片中的车辆均按照其坐标映射至待预测区域的对应栅格中，并将车辆行驶方向定义为车辆在栅格中的移动状态，移动状态包括向上、向下、向左、向右四种；对每个时间片内每个栅格中所包含的每一种移动状态的车辆总数进行统计并作为栅格值，将栅格值映射为矩阵元素，从而针对每个时间片中的不同移动状态分别构建车流二维矩阵，每一种移动状态在所有时间片中的车流二维矩阵按时间维度叠加形成车流三维矩阵；s3、从待预测时刻中提取小时字段和分钟字段，拼接形成一个二元的时间信号向量；s4、获取所有兴趣点空间地理位置，并将不同功能类别的兴趣点均映射到所述待预测区域的栅格中，对每个栅格中每一组功能类别的兴趣点总数进行统计并作为栅格值，将栅格值映射为矩阵元素，从而针对每一组功能类别的兴趣点分别构建一张二维矩阵形式的兴趣点切片，所有功能类别的兴趣点切片叠加形成三维张量形式的兴趣点信号；s5、将所述车流三维矩阵、所述时间信号向量和所述兴趣点信号作为经过训练的多方向交通流量预测模型的输入，所述多方向交通流量预测模型包含第一全连接神经网络、第二全连接神经网络、三维残差卷积网络和重校正层，由第一全连接神经网络根据所述时间信号向量输出时间信号矩阵，由第二全连接神经网络根据所述兴趣点信号输出兴趣点信号矩阵，由三维残差卷积网络根据所述车流三维矩阵、所述兴趣点信号矩阵和所述时间信号矩阵的融合特征输出结果矩阵，最终将结果矩阵在重校正层中经过带权压缩操作后，得到多方向交通流量预测结果。2.如权利要求1所述的基于兴趣点时空残差神经网络的多方向交通流量预测方法，其特征在于，在s1中，所述待预测区域为矩形区域；所述历史车流数据的时间跨度为[t，t+(m-1)*τ]，且其按τ分钟为时间片间隔提取m个车流数据时间片。3.如权利要求2所述的基于兴趣点时空残差神经网络的多方向交通流量预测方法，其特征在于，所述s2的具体实现步骤如下：s21、对待预测区域进行栅格化处理，共划分为i*j个栅格，第i行第j列的栅格为p
ij
；s22、每一个车流数据时间片中的车辆均按照其坐标映射至待预测区域的对应栅格中，并将车辆行驶方向定义为车辆在栅格中的移动状态，移动状态包括向上、向下、向左、向右四种；s23、对每个时间片内每个栅格中所包含的每一种移动状态的车辆总数进行统计，将时间片t内栅格p
ij
内移动状态为d的车流量记为将所有i*j个栅格对应的构建为时间片t内整个待预测区域中移动状态为d的车流二维矩阵车流二维矩阵中第i行第j列的元素值为s24、将所有m个车流数据时间片的车流二维矩阵按时间维度进行拼接，形成车流三维矩阵
4.如权利要求3所述的基于兴趣点时空残差神经网络的多方向交通流量预测方法，其特征在于，所述s3的具体实现步骤如下：将待预测时刻t
pred
转化为含有小时字段t
pred_hour
和分钟字段t
pred_minute
两个元素的时间信号向量h
t
＝[t
pred_hour
，t
pred_minute
]。5.如权利要求4所述的基于兴趣点时空残差神经网络的多方向交通流量预测方法，其特征在于，所述s4的具体实现步骤如下：s41、获取所有兴趣点的空间地理位置，并将所有兴趣点按照其位置均映射到所述待预测区域的各p
ij
栅格中；s42、将所有兴趣点按照不同功能类别分成n组，统计每个兴趣点分组g在栅格p
ij
内的兴趣点数量并作为栅格p
ij
的栅格值将每个兴趣点分组g对应的所有栅格的栅格值构建为兴趣点分组g对应的兴趣点切片γ
g
，γ
g
大小为i*j；s43、将所有n组兴趣点分组对应的兴趣点切片拼接，得到兴趣点信号ψ＝[γ1，γ2，...，γ
n
]，其大小为n*i*j。6.如权利要求5所述的基于兴趣点时空残差神经网络的多方向交通流量预测方法，其特征在于，所述s5的具体实现步骤如下：s51：将所述车流三维矩阵所述时间信号向量h
t
和所述兴趣点信号ψ输入经过训练的多方向交通流量预测模型中，所述多方向交通流量预测模型中包含第一全连接神经网络、第二全连接神经网络、三维残差卷积网络和重校正层；s52、所述时间信号h
t
被输入至包含l
ts
层全连接层级联的第一全连接神经网络中，其中第1层全连接层的输入为时间信号h
t
，下一层全连接层的输入为上一层全连接层的输出，最后一层全连接层的输出为将向量逐元素映射到大小为(i*j)的一个矩阵当中，得到大小为(i，j)的时间信号矩阵h
t
；s53、所述兴趣点信号ψ＝[γ1，γ2，...，γ
n
]输入后，先对其中的每个兴趣点切片γ
g
获取其平均自权重z
g
：得到兴趣点信号ψ的平均自权重矩阵z＝{z1，z2，...，z
n
}，其中n表示兴趣点分组的个数；然后将平均自权重矩阵z输入包含l
ps
层全连接层级联的第二全连接神经网络逐层计算，其中第1层全连接层的输入为平均自权重矩阵z，下一层全连接层的输入为上一层全连接层的输出，最后一层全连接层的输出为再后，采用门机制将输出映射为0到1之间的变量，其计算过程如下：其中f
si
为激活函数；最后，通过计算获取兴趣点信号矩阵
其中，
⊙
表示矩阵点乘；s54、将车流三维矩阵兴趣点信号矩阵和时间信号矩阵h
t
进行特征融合，中的第k个车流数据时间片的融合特征计算公式如下：其中，和均为可训练的参数，m为中车流数据时间片的个数；最终得到融合后的车流矩阵s55、将融合后的车流矩阵x
γ
输入包含l
c
层三维残差卷积层级联的三维残差卷积网络中逐层计算，其中第一层三维残差卷积层的输入车流矩阵x
γ
，每一层三维残差卷积层经过三维残差卷积操作后得到的结果作为下一层三维残差卷积层的输入，三维残差卷积网络中所有三维残差卷积层的输出拼接形成最终的结果矩阵x
st
；其中对于任意第l层三维残差卷积层，其中执行的三维残差卷积操作如下：首先，对当前三维残差卷积层的输入进行三维卷积操作得到卷积结果：其中，cov3d表示三维卷积操作，表示第l-1层三维残差卷积层的输出，其中1层三维残差卷积层的输出，其中和为第l层三维卷积操作中可训练的参数，f
c
为激活函数；然后，对三维卷积操作输出的卷积结果中每个元素进行批正则化操作，得到批正则结果最后，将批正则结果再与上一层的输出矩阵相加，得到第l层三维残差卷积层的输出矩阵公式为：s56、在重校正层中对最终输出的结果矩阵x
st
中所有维进行带权压缩操作，得到预测结果其计算公式如下：其中是可学习的参数矩阵。7.如权利要求1所述的基于兴趣点时空残差神经网络的多方向交通流量预测方法，其特征在于，所述多方向交通流量预测模型预先通过训练数据进行训练，训练过程中通过损失函数不断迭代预测结果与真实结果φ之间的损失值，当达到迭代终止条件时，输出多方向交通流量预测模型用于进行实际的预测。8.如权利要求1所述的基于兴趣点时空残差神经网络的多方向交通流量预测方法，其特征在于，所述的多方向交通流量预测模型的损失函数loss公式为：
其中为矩阵中各个元素的值，为矩阵φ中各个元素的值，m为训练样本数。9.如权利要求1所述的基于兴趣点时空残差神经网络的多方向交通流量预测方法，其特征在于，所述的m取5，时间片间隔τ取4。10.如权利要求1所述的基于兴趣点时空残差神经网络的多方向交通流量预测方法，其特征在于，所述的车辆在栅格中的移动状态根据车辆行驶轨迹经过该栅格时的车辆行驶方向确定。

技术总结
本发明公开了一种基于兴趣点时空残差神经网络的多方向交通流量预测方向。该方法基于短时车辆轨迹数据，通过加入时间信号和兴趣点信号来增强时空特征，利用3D CNN提取车流变化随时间变化的时空特征，结合残差神经网络避免模型过拟合，最后对带有时空特征信息进行加权压缩，输出带有移动状态的车流分布矩阵，实现了对区域交通流量的预测。本发明方法良好地关联了车流的空间分布随时间变化的特征，并综合了时间信号与兴趣点信号关系，有精度高、实用性强等特点，可以为公共设施部署、交通导流、用地规划等提供决策支撑。地规划等提供决策支撑。地规划等提供决策支撑。


技术研发人员：何宏 王欣峰 孙笑笑 俞东进
受保护的技术使用者：杭州电子科技大学
技术研发日：2021.12.10
技术公布日：2022/3/8