本发明涉及智能交通技术领域，尤其是一种基于多源融合数据的高速公路交通流参数修正方法。
背景技术：
随着信息技术的不断发展，交通运营管理中心不仅可以通过各种采集设备获取大量的交通流参数数据，尤其是高速公路上安装大量的固定检测器，可以实时检测交通流参数包括：交通量、速度和占有率，且能够获取大量的气象数据。但固定检测器由于通信、通电等原因常常发生数据缺失问题，为后续的交通数据挖掘带来较大困难。
现有技术中，交通流参数修复的方法包括时间序列及神经网络等各种智能算法。然而，这些算法不能够充分挖掘数据信息，并且多采用单一的数据源，缺乏考虑天气等因素对交通流的影响，不能够满足智能交通系统的发展要求。
技术实现要素：
本发明所要解决的技术问题在于，提供一种基于多源融合数据的高速公路交通流参数修正方法，可以对缺失数据进行修正，提高修复精度，满足实时处理的需求。
为解决上述技术问题，本发明提供一种基于多源融合数据的高速公路交通流参数修正方法，包括如下步骤：
(1)从相应的交通流检测设备和气象检测设备中提取交通流参数数据和气象数据，并对数据进行时间和空间维度的匹配；
(2)对缺失数据进行筛选；
(3)对数据进行修复建模；
(4)对数据进行修复；
(5)数据返归一化；将返归一化结果数据插入到矩阵X的相应位置，得到完整的数据矩阵。
优选的，步骤(1)中，提取交通流参数数据，根据检测设备的经纬度位置提取距离最近的气象监测设备中的气象数据，完成空间维度的匹配；选取时间间隔，将两种数据进行转换，使两种数据具有相同的时间间隔，时间间隔取两种不同数据的时间间隔的公倍数，完成时间维度的匹配。
优选的，步骤(2)中，对缺失数据进行筛选的具体步骤为：设数据矩阵其中该矩阵包含m行n列，其中m行表示时间段的个数，n列为相应的交通流参数和气象参数，表示第i列第j行的数据；从数据矩阵中筛选出缺失数据，将第i列第j行的缺失数据定义为xij。
优选的，步骤(3)中，对数据进行修复建模的具体步骤为：
(a)记录缺失数据所在行的编号集合M，将其从数据矩阵X中提出出来为数据矩阵Xmiss，其余数据组成新的数据矩阵X′，然后将矩阵进行归一化处理；针对数据矩阵X′，随机产生30％的交通流缺失数据，记录缺失数据所在行的编号集合M′，将其从数据矩阵X′中提取出来构成新的数据矩阵为X′test，剩余数据矩阵为X′train；
(b)基于多元时间序列模型和数据矩阵X′train分别建立各项交通流数据与气象数据的回归方程，包括模型一：流量与湿度、风速、降雨量、温度的回归方程；模型二：速度与流量、湿度、风速、降雨量、温度的回归方程；模型三：占有率与速度、流量、湿度、风速、降雨量、温度的回归方程；
所述的多元时间序列如下式：
yt＝δ+Φ1yt-1+…+Φpyt-p+εt-Θ1εt-1-…-Θpεt-p
式中：yt＝(y1t，...，ykt)′，t＝0，1，...表示k维时间序列向量；ykt表示第k个变量在t时刻的观测值；εt＝(ε1t，...，εkt)′为白噪声向量，当t≠s时，满足E(εt)＝0，E(εtε′s)＝0；Φ1，...，Φp，Θ1，...，Θp为模型中需求解的参数矩阵；
(c)将集合X′test中对应变量的数据依次带入到模型一、二、三中计算交通流参数的数值，得到交通流参数数据矩阵Y＝{yij}，并将数据按照其时间位置插入到数据矩阵X′train中得到数据矩阵X″train；
(d)对数据矩阵X″train，按照时间进行降序排列，得到数据矩阵X″′train和缺失数据编号集合M″；利用多元时间序列模型和数据矩阵X″′train分别建立各项交通流数据与气象数据的回归方程，包括模型四：流量与速度、占有率、湿度、风速、降雨量、温度的回归方程；模型五：速度与流量、占有率、湿度、风速、降雨量、温度的回归方程；模型六：占有率与速度、流量、湿度、风速、降雨量、温度回归方程；
(e)将集合M″中的数据依次带入到模型四、五、六中计算缺失的交通流参数数据，得到数据矩阵Y′＝{y′ij}；
(f)设数据矩阵X′test中交通流参数数据矩阵为利用最小二乘原理求解下式中的参数a、b、∈：
优选的，步骤(4)中，对数据进行修复的具体步骤为：将数据矩阵Xmiss归一化，其公式如下：
上式中X1为归一化序列，Xmin和Xmax分别为归一化序列中的最小值和最大值，Xnorm为归一化结果；
利用模型一、二、三求得数据矩阵模型四、五、六求得数据矩阵缺失数据可通过下式求得：
优选的，步骤(5)中，数据返归一化的具体步骤为：
上式中，为待返归一化序列，分别为待返归一化序列中的最大值和最小值，X2为返归一化结果；将返归一化结果数据插入到矩阵X的相应位置，得到完整的数据矩阵。
本发明的有益效果为：充分利用数据信息，包括缺失数据中数据的信息，对缺失数据进行修正；考虑天气因素对交通流的影响，可以提高修复精度；提出的算法简单明了，能够满足实时处理的要求。
附图说明
图1为本发明的方法流程示意图。
具体实施方式
如图1所示，一种基于多源融合数据的高速公路交通流参数修正方法，包括如下步骤：
(1)从相应的交通流检测设备和气象检测设备中提取交通流参数数据和气象数据，并对数据进行时间和空间维度的匹配。
首先提取交通流参数数据，根据检测设备的经纬度位置提取距离最近的气象监测设备中的气象数据，完成空间维度的匹配。然后选取时间间隔，将两种数据进行转换，使两种数据具有相同的时间间隔，完成时间维度的匹配。
所述的时间间隔应取两种不同数据的时间间隔的公倍数。
(2)缺失数据筛选
设数据矩阵其中该矩阵包含m行n列，其中m行表示时间段的个数，n列为相应的交通流参数和气象参数，表示第i列第j行的数据。从数据矩阵中筛选出缺失数据，将第i列第j行的缺失数据定义为xij。
(3)数据修复建模
3.1记录缺失数据所在行的编号集合M，将其从数据矩阵X中提出出来为数据矩阵Xmiss，其余数据组成新的数据矩阵X′，然后将矩阵进行归一化处理。针对数据矩阵X′，随机产生30％的交通流缺失数据，记录缺失数据所在行的编号集合M′，将其从数据矩阵X′中提取出来构成新的数据矩阵为X′test，剩余数据矩阵为X′train。
3.2基于多元时间序列模型和数据矩阵X′train分别建立各项交通流数据与气象数据的回归方程，包括模型一：流量与湿度、风速、降雨量、温度的回归方程；模型二：速度与流量、湿度、风速、降雨量、温度的回归方程；模型三：占有率与速度、流量、湿度、风速、降雨量、温度的回归方程。
所述的多元时间序列如下式：
yt＝δ+Φ1yt-1+…+Φpyt-p+εt-Θ1εt-1-…-Θpεt-p
式中：yt＝(y1t，...，ykt)′，t＝0，1，...表示k维时间序列向量。ykt表示第k个变量在t时刻的观测值。εt＝(ε1t，...，εkt)′为白噪声向量，当t≠s时，满足E(εt)＝0，E(εtε′s)＝0。Φ1，...，Φp，Θ1，...，Θp为模型中需求解的参数矩阵。
3.3将集合X′test中对应变量的数据依次带入到模型一、二、三中计算交通流参数的数值，得到交通流参数数据矩阵Y＝{yij}，并将数据按照其时间位置插入到数据矩阵X′train中得到数据矩阵X″train。
3.4对数据矩阵X″train，按照时间进行降序排列，得到数据矩阵X″′train和缺失数据编号集合M″。利用多元时间序列模型和数据矩阵X″′train分别建立各项交通流数据与气象数据的回归方程，包括模型四：流量与速度、占有率、湿度、风速、降雨量、温度的回归方程；模型五：速度与流量、占有率、湿度、风速、降雨量、温度的回归方程；模型六：占有率与速度、流量、湿度、风速、降雨量、温度回归方程。
3.5将集合M″中的数据依次带入到模型四、五、六中计算缺失的交通流参数数据，得到数据矩阵Y′＝{y′ij}。
3.6设数据矩阵X′test中交通流参数数据矩阵为利用最小二乘原理求解下式中的参数a、b、∈：
(4)数据修复
将数据矩阵Xmiss归一化，其公式如下：
上式中X1为归一化序列，Xmin和Xmax分别为归一化序列中的最小值和最大值，Xnorm为归一化结果。
利用模型一、二、三求得数据矩阵模型四、五、六求得数据矩阵缺失数据可通过下式求得：
(5)数据返归一化
上式中，为待返归一化序列，分别为待返归一化序列中的最大值和最小值，X2为返归一化结果。将返归一化结果数据插入到矩阵X的相应位置，得到完整的数据矩阵。
尽管本发明就优选实施方式进行了示意和描述，但本领域的技术人员应当理解，只要不超出本发明的权利要求所限定的范围，可以对本发明进行各种变化和修改。