一种基于OD数据增量迭代式估计路段交通流量的方法与流程

本发明属于交通数据挖掘和应用领域，涉及利用od数据计算路段交通量的方法。

背景技术：

道路交通流量是指单位时间内通过道路某断面的交通流量。用于支撑道路建设项目可行性研究阶段的现状评价、综合分析，也是确定道路建设项目的建设规模、技术等级、工程设施、经济效益评价及几何线形设计的主要依据。

目前道路交通量的主要获取技术有线圈、视频、激光、微波等，通过在道路路段上安装监测设备进行采集，安装成本高，长期维护难度大、无法开展个体车辆追踪等特征，导致无法大规模安装应用，在路网中覆盖有限，目前无法获得全路网的时空交通量。

随着互联网技术的飞速发展，以出行个体为载体的交通信息获取方式得到大面积应用，如手机/车载终端导航数据、网约车订单数据、高速公路收费车辆数据，覆盖用户数亿。此类数据的特征是包括用户的出行起点和终点位置、时间等信息，本专利将具备此特征的数据称为od数据。通过对od数据的挖掘与分析，能够获得od数据背后的路段交通量和行驶速度信息，为交通管理提供数据支撑与决策支持。

技术实现要素：

本发明目的在于提出一种基于od数据增量迭代式估计路段交通流量的方法，用于获取不同路段不同时段的交通流量。

一种基于od数据增量迭代式估计路段交通流量的方法，其特征在于包含以下步骤：

步骤1：od数据预处理，删除出发时间早于到达时间的数据，形成有效od数据集，；

步骤2：基于有效od数据，根据车辆出发和到达的节点，计算途经的最短路径长度pk，其中k表示第k条od数据，以最短路径pk由小到大排序，形成数据集合dk(to，td，p)，p表示pk的集合，to表示出发时间，td表示到达时间；

步骤3：将路网按照节点为间隔划分成若干路段单元，标记为ri，其中i表示第i条路段，最短路径pk由多个路段单元构成，将时间以一定的间隔均匀划分为若干时间段，标记为ij，其中j表示第j个时间间隔段，路段单元ri在某时间段ij的状态用集合c表示，标记为其中l表示路段单元长度，表示交通流量，表示平均速度；

步骤4：将以最短路径长短排序后的有效0d数据从小到大依次迭代计算更新计算路段单元的过程如下：

4.1在有效od数据中，抽取最短路径pk＜pf值的数据，pf为迭代计算的初始距离阈值，其选取相邻节点平均距离作为经验值，当pk小于pf时，认为车辆在该距离上，速度运行状态不会发生较大变化，计算该路径的平均速度vk，vk＝p/(td-to)，将车辆出发和到达时间to，td与ij进行比较，标记与ij重合的时段，标记该最短路径途经的所有路段单元在重合时段的状态其中值更新为n为路段单元在ij时间段的记录数量；

4.2在有效od数据中，继续抽取最短路径x1pf≤pk＜x2pf(x1≥1，x2≥2)值的数据。当最短路径途经的路段单元在重合的ij时间段均存在已标记的状态时，用已标记的状态作为权重，计算最短路径途经的路段单元的速度值计算过程如下：

式中：m为最短路径途经的路段单元数量，m为最短路径途经的第m个路段单元，为第k条od数据最短路径途经第m个路段单元的平均速度，为最短路径途经的第1个路段单元在重合时段i1下已标记的速度，和以此类推，标记最短路径途经的所有路段单元在重合时段的状态其中值更新为值更新为与原的平均值；

当最短路径途经的路段单元在重合的ij时间段部分存在已标记的状态时，根据已标记路段单元的最短路径长度，计算未标记路段的状态值，先计算该路径的平均速度vk，vk＝p/(td-to)，将车辆出发和到达时间to，td与ij进行比较，标记与ij重合的时段，标记该最短路径途经的所有路段单元在重合时段的状态将该状态赋给未标记路段单元，同时所有路段单元在重合的ij时间段值更新为n为路段单元在ij时间段的记录数量；

当最短路径途经的路段单元在重合的ij时间段均不存在已标记的状态时；计算该路径的平均速度vk，vk＝p/(td-to)，将车辆出发和到达时间出发和到达to，td与ij进行比较，标记与ij重合的时段，标记该最短路径途经的所有路段单元在重合时段的状态其中值更新为n为路段单元在ij时间段的记录数量；

4.3在余下的有效0d数据中，每次增加最短路径一个pf值的范围，按4.2的步骤完成所有0d数据的计算；

步骤5：汇总的q值，得到所有路段单元在不同时间段上的交通量；

附图说明

图1.基于od数据的路段交通量计算流程

图2.路段节点示意图

具体实施方式

本实施方案以道路od数据为例，结合附图对本发明详细描述具体实施流程。实施流程如附图1，路段节点如附图2。

步骤1：对od数据进行预处理，将出发时间早于到达时间的od数据删除，形成有效od数据集；

步骤2：基于有效od数据，根据车辆出发和到达的节点，计算途经的最短路径长度pk，其中k表示第k条od数据，以最短路径pk由小到大排序，形成数据集合dk(to，td，p)，to表示车辆出发时间，td表示车辆到达时间，如表2；

表2

步骤3：本例中将节点1→节点2路段标记为r1，节点2→节点3路段标记为r2，节点3→节点4路段标记为r3；本例以5分钟的时间间隔将全天划分为288个时段，标记为ij，如时段00:00:00-00:05:00为第1个时段，时段23:55:00-24:00:00为第288个时段；

步骤4：计算相邻路段的平均距离26720米，并以此作为经验值，pf选取为26720米；

4.1在有效od数据中，抽取最短路径长度满足pk≤pf的od数据，结果如表3：

表3

序号1的车辆在时间08:05:06进入r1，对应第86个时段，同理在第87个时段离开r1；序号2的车辆在第85个时段进入r1，在第87个时段离开r1；序号3的车辆在第86个时段进入r2，在第89个时段离开r2，序号4的车辆在第85个时段进入r2，在第88个时段离开r2，计算路径的平均速度vk，vk＝p/(td-to)：

(千米/小时)

(千米/小时)

(千米/小时)

(千米/小时)

将车辆出发和到达时间td，to与ij进行比较，标记与ij重合的时段，标记该最短路径途经的所有路段单元在重合时段的状态

序号1车辆和序号2的车辆在第86个时段和第87个时段都经过r1，因此第86个时段和第87个时段r1速度值为序号1和序号2的车辆速度的平均值，即：

(千米/小时)

序号3和序号4的车辆在第86个时段、第87个时段和第88个时段都经过r2，因此第86个时段、第87个时段和第88个时段r2的速度值为序号3和序号4的车辆速度的平均值，即：

(千米/小时)

故为(16740，2，108.60)，为(16740，2，108.60)，同理：为(16740，1，99.94)，为(26700，1，107.40)，为(26700，2，108.32)，为(26700，2，108.32)，为(26700，2，108.32)，为(26700，1，109.23)；

4.2提取路径长度满足pf＜pk≤2pf的od数据，满足条件的数据如表4；

表4

序号5的车辆在第85个时段进入r1，在第87个时段离开r2；序号6的车辆在第87个时段进入r2，在第89个时段离开r1；序号7的车辆在第91个时段进入r2，在第93个时段离开r3；

(1)序号5：r1和r2在第85-87时段内均已标记状态；根据发明内容步骤4.2中的说明，序号5的各时段速度如下：

第85时段：

(千米/小时)

(千米/小时)

第86时段：

(千米/小时)

(千米/小时)

第87时段：

(千米/小时)

(千米/小时)

第85个时段-第87个时段均有状态标记，因此序号为5的车辆经过后，r1和r2的状态值均有更新，根据计算结果得到序号5的车辆第85时段和第86时段在r1上，第86时段和第87时段在r2上，故：

(千米/小时)

(千米/小时)

(千米/小时)

故更新后的状态值：为(16740，2，130.87)，为(16740，3，135.20)，为(26700，3，151.90)，为(26700，3，151.90)。

(2)序号6：r1在第88、89时段无状态标记，r1和r2在第87时段有状态标记。根据发明内容步骤4.2中的说明，序号6的各时段速度如下：

第87时段：

(千米/小时)

(千米/小时)

第88时段：

根据第88时段r2的标记状态值计算车辆在r1上行驶的时长并计算时段ij的速度值。

od数据序号6的车辆在r2行驶650s，在r1上行驶189秒，故车辆第87、88时段在r2，第89时段在r1，故：

(千米/小时)

(千米/小时)

(千米/小时)

故更新后的状态值为为(26700，4，183.51)，为(26700，3，108.32)，为(16740，1，318.85)。

(3)序号7的车辆经过的时段r1和r2均无标记状态。

(千米/小时)

车辆第91时段和第92时段在r2上，第93时段在r3上.

故：为(26700，1，109.23)，为(26700，1，109.23)，为(26700，1，109.23)。

步骤5：汇总状态值，得到所有路段单元在不同时间段上的交通流量，结果见表5.

表5

最后应说明的是：以上示例仅用以说明本发明而并非限制本发明所描述的技术方案；因此，尽管本说明书参照上述的示例对本发明已进行了详细的说明，但是，本领域的普通技术人员应当理解，仍然可以对本发明进行修改或等同替换；而一切不脱离发明的精神和范围的技术方案及其改进，其均应涵盖在本发明的权利要求范围当中。