[发明专利]一种基于时空相关性的交通流量数据恢复方法

摘要

本发明公开了一种基于时空相关性的交通流量数据恢复方法。该方法包括前期交通数据转换、异常数据筛选方法和后续的交通流数据恢复方法。其中数据筛选方法，根据实际交通流数据的异常情况涉及了阈值筛选、零数据筛选和质量筛选法。然后根据交通流数据在时间上和空间上的相关性，结合时间序列法和多元线性回归法，即基于时间相关性和空间相关性的数据恢复方法，设计了基于时空相关性的综合交通流数据恢复方法。本发明方法简单快速，能够满足实时处理的要求，且所得结果的精确度高。

主权项

一种基于时空相关性的交通流量数据恢复方法，其特征在于该方法包括以下步骤：步骤（1）.交通流量数据格式转换：将检测点由车辆检测器检测得到的原始交通流量数据q_orig单位换算成标准小时流量数据，然后对该标准小时流量数据进行加权平均转换成固定周期T的交通流量数据集q；所述的固定周期T的交通流量数据集q={q(1),…,q(k),…,q(n)}；q(k)表示k时间段的交通流量数据，其中所述的标准小时流量数据表示单位为车辆/小时的交通流量数据；步骤（2）.交通流量数据筛选2.1阈值筛选：对步骤（1）得到的固定周期T的交通流量数据集q进行阈值筛选，即对交通流量数据集q中的交通流量数据q(k)进行逐个比较，根据公式（1）进行筛选，剔除交通流量数据集q中不符合公式（2）的数据：0＜q(k)＜Q_max    (1)；Q_max＝f_c×C    (2)；其中阈值Q_max的单位为veh/h；C表示车道的基本通行能力，单位为veh/h；f_c为常数；2.2零数据筛选：对步骤2.1筛选后得到的交通流量数据集进行零数据筛选，具体步骤如下：根据泊松分布原理得到时间间隔t内到达车辆数的概率P(x)见公式（3）：$P\left(x\right)=\frac{{\left(\mathrm{\λt}\right)}^x{e}^{-\mathrm{\λt}}}{x!}---\left(3\right);$式中P(x)为时间间隔t内到达x辆车的概率，x为自然数；λt为时间间隔t内平均到达的车辆数；若置信水平α＝0.01，采样间隔内有车辆通过，即x＞0的概率为：$P(x>0)=1-P\left(0\right)=1-\frac{{\left(\mathrm{\λt}\right)}^0{e}^{-\mathrm{\λt}}}{0!}=0.99---\left(4\right);$由此可得：λt≈4.6(veh/(tmin))    （5）；$Q_0=\mathrm{\λt}\×\frac{60}{t}(\mathrm{veh}/h)---\left(6\right);$其中Q₀表示时间间隔t内的标准小时交通流量数据；若时间间隔t等于固定周期T，当k时段的历史平均交通流量大于Q₀，且实际检测到的交通流量数据q(k)显示为零时，则从步骤2.1筛选后得到的交通流量数据集中筛选出该数据，如公式（7）所示：2.3质量筛选：对步骤2.2筛选后的交通流量数据集进行质量筛选，剔除不符合式（8）的数据；$\min \{\stackrel{\‾}{q}\left(k\right)-2{\mathrm{\σ}}_q,q(k-1)+\stackrel{\‾}{d}\left(k\right)-{2\mathrm{\σ}}_d\}\≤q\left(k\right)\≤\max \{\stackrel{\‾}{q}\left(k\right)+{2\mathrm{\σ}}_q,q(k-1)+\stackrel{\‾}{d}\left(k\right)+{2\mathrm{\σ}}_d\}---\left(8\right);$其中为k时段的历史平均交通流量，σ_q为的方差；为k时段与k‑1时段的流量平均差值，σ_d为的方差；步骤（3）.交通流量数据恢复：3.1利用基于时间相关性的交通流量数据恢复方法，得到缺失数据车道流量数据的预测值：首先采用简单移动平均法构造时间序列模型以得到步骤（2）筛选后的交通流量数据中缺失数据车道流量数据的预测值见公式（9）：${\hat{q}}_T\left(k\right)=\frac{q(k-1)+q(k-2)+\·\·\·+q(k-n)}{n}---\left(9\right);$其中n为移动平均个数；3.2利用基于空间相关性的交通流量数据恢复方法，得到缺失数据车道流量数据预测值：首先采用多元线性回归分析，以步骤（2）筛选后的交通流量数据中缺失数据车道流量数据检测点交叉口的上下游交叉口的交通流量数据为自变量，建立多元回归模型，得到当前交叉口交通流量预测值然后乘以数据缺失车道j在该时段k的车道占有率f_j(k)，得到缺失数据车道流量数据预测值见公式(10)；${\hat{q}}_S\left(k\right)=\hat{Q}\left(k\right)\×f_j\left(k\right)---\left(10\right);$若上游交叉口有a个方向车流驶向检测路段，同时有驶出路段的车辆分b个方向驶向下游交叉口,则共有i＝a+b个自变量Q₁(k)、Q₂(k)、……、Q_i(k)；$\hat{Q}\left(k\right)={\mathrm{\β}}_0+{\mathrm{\β}}_1{Q}_1\left(k\right)+{\mathrm{\β}}_2{Q}_2\left(k\right)+\·\·\·+{\mathrm{\β}}_i{Q}_i\left(k\right)---\left(11\right);$式中β₀、……、β_i为回归系数；$f_j\left(k\right)=\frac{{\mathrm{Hq}}_j\left(k\right)}{\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{Hq}}_j\left(k\right)}---\left(12\right);$式中Hq_j(k)表示数据缺失车道j在时段k内的历史交通流量，m为检测点车道数；3.3利用基于时空相关性的交通流量数据恢复方法，得到最终的缺失数据车道流量数据恢复值：根据公式（13）、（14）计算上述步骤3.1、3.2两模型预测值的百分比相对误差R_T和R_S，并根据公式（15）确定综合模型的权重系数θ；$R_T=\left(\right|{\hat{q}}_T\left(k\right)-q\left(k\right)|/q\left(k\right))\×100\%---\left(13\right);$$R_S=\left(\right|{\hat{q}}_S\left(k\right)-q\left(k\right)|/q\left(k\right))\×100\%---\left(14\right);$$\mathrm{\θ}=\frac{R_T}{R_T+R_S}---\left(15\right);$最后根据式（16）得到最终的缺失数据车道流量数据恢复值；$\hat{q}\left(k\right)=\mathrm{\θ}\widehat{q_T}\left(k\right)+(1-\mathrm{\θ})\widehat{q_S}\left(K\right)---\left(16\right);$步骤4.根据公式（17）将步骤3得到的最终的缺失数据车道流量数据恢复值补充至步骤2筛选后的交通流量数据集中；