CN109035784A - 基于多源异构数据的动态车流od估计方法 - Google Patents

基于多源异构数据的动态车流od估计方法 Download PDF

Info

Publication number
CN109035784A
CN109035784A CN201811083739.8A CN201811083739A CN109035784A CN 109035784 A CN109035784 A CN 109035784A CN 201811083739 A CN201811083739 A CN 201811083739A CN 109035784 A CN109035784 A CN 109035784A
Authority
CN
China
Prior art keywords
traffic
flow
observable
section
road
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.)
Granted
Application number
CN201811083739.8A
Other languages
English (en)
Other versions
CN109035784B (zh
Inventor
周东
李璐
吕伟韬
潘阳阳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
JIANGSU INTELLIGENT TRANSPORTATION SYSTEMS Co Ltd
Original Assignee
JIANGSU INTELLIGENT TRANSPORTATION SYSTEMS Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by JIANGSU INTELLIGENT TRANSPORTATION SYSTEMS Co Ltd filed Critical JIANGSU INTELLIGENT TRANSPORTATION SYSTEMS Co Ltd
Priority to CN201811083739.8A priority Critical patent/CN109035784B/zh
Publication of CN109035784A publication Critical patent/CN109035784A/zh
Application granted granted Critical
Publication of CN109035784B publication Critical patent/CN109035784B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/0137—Measuring and analyzing of parameters relative to traffic conditions for specific applications
    • G08G1/0145—Measuring and analyzing of parameters relative to traffic conditions for specific applications for active traffic flow control
    • 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/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • 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
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/40—Business processes related to the transportation industry

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • Development Economics (AREA)
  • Educational Administration (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明提供一种基于多源异构数据的动态车流OD估计方法，将视频号牌识别数据和车辆GPS数据进行整合，将路网断面流量拆分为可观测OD和不可观测OD，同时基于卡尔曼滤波对不可观测OD进行估计，进而与可观测OD整合得到动态车流OD信息；该种基于多源异构数据的动态车流OD估计方法，将号牌识别设备采集的过车数据与GPS数据相结合，提取出有效路径作为可观测OD需求，并划分为可观测部分和不可观测部分来实现路网OD估计，对比传统仅基于断面数据进行估计，该方法更能捕捉到OD变化趋势，更具准确性。

Description

基于多源异构数据的动态车流OD估计方法
技术领域
本发明涉及一种基于多源异构数据的动态车流OD估计方法。
背景技术
车辆起-讫点(Origin-Destination，简称OD)流量作为道路交通管控的重要数据依据以及各种中微观交通模型及仿真平台的重要输入参数，OD流量的估计问题已逐渐发展为交通领域中的一个重要研究内容。随着智能交通运输系统的快速发展，现代交通运营与管控更加强调精细化、实时性，因此，动态车流OD估计已经成为起讫点流量估计问题中的研究重点。
目前常规断面流量的动态OD估计主要有最小二乘模型、双层规划模型、极大熵模型、贝叶斯模型和状态空间模型、卡尔曼滤波模型等。尤其是卡尔曼滤波已经成为近年来动态OD的主要估计方法之一，如专利201510374265.2提出一种“基于AFC数据的公交客流OD实时估计方法”，通过挖掘公交客流OD与公交站点客流之间的映射关系，构建基于卡尔曼滤波的公交客流OD实时估计模型，实现公交客流OD的实时估计。而随着视频号牌识别技术的发展，与传统“断面型”交通检测技术相比，视频号牌识别设备可根据处理得到车辆的出行信息，从而获得某一具体点对OD时间的真是流量，但是由于路网内视频号牌识别设备未能全部覆盖因此无法获取完整的轨迹信息，只能针对高速公路、城市快速路等“封闭型”路网，对城市路网的估计准确性有待加强。
因此针对路网动态车流OD的估计问题，应该充分利用城市路网中各种交通检测设备的观测值，并融合多源交通流数据，实现路网动态车流OD的估计。
发明内容
本发明的目的是提供一种基于多源异构数据的动态车流OD估计方法解决现有技术中存在的由于路网内视频号牌识别设备未能全部覆盖因此无法获取完整的轨迹信息，只能针对高速公路、城市快速路等“封闭型”路网，对城市路网的估计准确性有待加强的问题。
该种基于多源异构数据的动态车流OD估计方法，基于车辆GPS定位数据和视频号牌数据将断面客流划分为可观测OD需求数据和不可观测OD需求数据，并基于动态交通分配理论建立OD流量与断面流量之间的动态映射关系，从而构建了带线性状态约束条件的卡尔曼滤波估计模型对不可观测OD需求进行估计求解，进而整合可观测OD需求获得路网动态OD需求数值，对比传统的动态OD估计方法，大大提高了估计值的准确度，将其更为贴近路网实际通行情况，为交通规划与管理提供了更为精准的参数，提升了规划管理的合理性。
本发明的技术解决方案是：
一种基于多源异构数据的动态车流OD估计方法，将视频号牌识别数据和车辆GPS数据进行整合，将路网断面流量拆分为可观测OD和不可观测OD，同时基于卡尔曼滤波对不可观测OD进行估计，进而与可观测OD整合得到动态车流OD信息，包括以下步骤，
S1、基于用地性质、出行热点、视频号牌识别设备位置及路网拓扑结构将路网划分为若干个路网级的交通小区；
S2、基于视频号牌识别数据和车辆GPS数据提取出各交通小区的车辆轨迹信息；
S3、对缺失的轨迹数据进行补充和重构，根据车辆的完整行驶轨迹以及实际交通流量实现对路网断面流量观测值的拆分为可观测OD和不可观测OD；
S4、基于动态交通分配理论建立OD流量与断面流量之间的动态映射关系，建立OD分配矩阵；
S5、基于卡尔曼滤波对不可观测OD进行估计，进一步与可观测OD整合，对路网的动态需求进行求解。
进一步地，步骤S3具体为，
S31、基于最短路算法对由视频车牌识别数据获得的车流OD缺失部分进行补全，同时基于GPS定位数据提取出车辆完整的行驶轨迹；
S32、通过车辆完整的行驶轨迹信息，获得其在路网上的传播情况，即其经过路网上检测断面的实际时间，从而获得相应断面处的断面流量和其中表示t时段内，第i个进口道的由可观测部分构成的断面流量，i∈Lc，Lc表示路网中可观测进口道集合；表示t时段内，第i个路段由可观测OD需求构成的断面流量，i∈L，L表示路网中全部可观测路段集合；
S33、基于实际的断面交通流量对流量的观测值进行拆分，即将实际断面交通流量减去可观测的断面流量数值用于估计OD需求不可观测部分的断面流量和其中表示t时段内，第i个路段由不可观测OD需求构成的断面流量，i∈L；表示t时段内，第i个进口道的由不可观测部分构成的断面流量，i∈Lc，Lc表示路网中可观测进口道集合。
进一步地，步骤S4具体为，
S41、基于Logit模型建立路径选择模型，将行程时间作为效益函数，即OD出行中选择各路径的比例如下：
其中，为在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s的比例，其中o，d∈Q，Q为交通小区集合；且o≠d，s∈R(o，d)，t；ts表示有效路径集合R(o，d)，t中，车辆行驶的行程时间；ti为有效路径集合R(o，d)，t中，车辆在第i条路径上的行程时间；为有效路径集的平均行程时间；R(o，d)，t表示在t时段内由交通小区o出发且目的地为交通小区d的出行路径集合，其中o，d∈Q且o≠d；
S42、针对车流在具体路径上的传播过程得到路径流量与交叉口进口道断面流量之间的动态映射关系数值；
S43、得到某一具体路段或交叉口进口道处的断面流量观测值与待估计的OD需求之间的映射关系，从而建立OD分配矩阵，其中OD需求之间的映射关系为：
式中，为在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达路段i的比例；为在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达进口道i的比例；表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s并在h时段内到达路段b的需求量占路径s流量的比例；表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s并在h时段内到交叉口进口道c的需求量占路径s流量的比例；表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s的比例，其中o，d∈Q且o≠d，b∈s，s∈R(o，d)，t，h＝t，t+1，…t+δ-1，δ为完成出行的时段数。
进一步地，步骤S41中，对于OD点对之间的有效路径集R(o，d)，t从已获得的部分车辆路径信息确定，具体采用k最短路算法计算，或在最短路行程时间TTmin的基础上，适当乘以一个权重系数ρ，ρ≥1，从而得到OD点对之间的行程时间范围[TTmin,ρTTmin]，再得到有效路径集，其中权责系数ρ根据所提取的可观测OD需求的实际行程时间和最短算法得到的路径行程时间确定。
进一步地，步骤S42具体为，将具体某一时段内选择某一路径的OD需求其中为在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s的比例，x(o，d)，t表示由交通小区o在t时段内出发，其目的地为交通小区d的实际需求量，o，d∈Q，Q为交通小区集合；
抽象成若干车辆构成的车辆元胞进行处理，则在t时段内由交通小区o出发到达交通小区d，选择路径s并在h时段内到达路段c的需求量占路径s流量的比例如下：
若(t-1)Δ≤γl＜γf＜tΔ，则
若γl＜(t-1)Δ≤γf＜tΔ，则
若(t-1)Δ≤γl＜tΔ≤γf＜(t+1)Δ，则
式中，表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s并在h时段内到交叉口进口道c的需求量占路径s流量的比例，其中o，d∈Q且o≠d，b∈s，s∈R(o，d)，t，h＝t，t+1，…t+δ-1，δ为完成出行的时段数；γf表示车辆元胞中第一辆车和最后一辆车到达某一路段或某一交叉口进口道的时刻；γl表示车辆元胞到达相应路段直至其完全消失的动态过程，且γf＞γl；Δ为每个时段的固定时长。
进一步地，步骤S5具体为，
S51、基于步骤S3拆分得到可观测的OD需求和不可观测的OD需求，在带约束的卡尔曼滤波构建动态OD估计模型基础下对不可观测的OD需求进行估计，建立出状态转移方程、状态线性约束方程和观测方程；
S52、将可观测OD需求和估计出的不可观测OD需求结合，构建出基于多源数据的动态OD估计模型，从而求解出路网动态OD估计值，具体公式如下：
式中，x(o，d)，t表示由交通小区o在t时段内出发，其目的地为交通小区d的实际需求量，其中o，d∈Q且o≠d；表示由交通小区o在t时段内出发，其目的地为交通小区d的观测需求量，其中o，d∈Q且o≠d；表示由交通小区o在t时段内出发，其目的地为交通小区d的待估计需求量，其中o，d∈Q且o≠d。
进一步地，步骤S51中，建立出状态转移方程为，
式中，为t时段内路网不可观测OD需求的列向量，为m×1维向量组；m为待估计的OD点对的个数；Wt为m×1维向量组，即为状态转移偏差ωt的向量化；
S512、基于视频车牌识别设备在路网上的分布对部分交通小区的流入路网总需求进行观测，并根据可观测小区的需求产生量构建状态线性约束条件方程；
S513、建立待估计OD需求与实际路网中采集到的断面流量之间映射关系的观测方程。
进一步地，步骤S512中，建立状态变量的线性约束方程为，
其中，为s×1维列向量，表示的向量化数值，其中s为可进行需求全样本观测的交通小区的个数；Dt为s×m维矩阵，即为交通小区的OD分配矩阵，为t时段内路网不可观测OD需求的列向量，为m×1维向量组。
进一步地，步骤S513中，建立观测方程为，
其中，Yt为(n+θ)×1维列向量，是OD需求向量数和的简化合并，即Ah为(n+θ)×m维矩阵，即到达比例向量数Ph和的合并，Ut为(n+θ)×1维列向量，即检测误差μt和的合并，
其中和分别为n×1维向量组，是li，t和的向量化数值，即为：
其中，Ph为n×m维矩阵，是的向量化，m为待估计的OD点对数量；为θ×m矩阵，是的向量化；μt和分别为θ×1维向量组，是和的向量化，为x(o，d)，h的向量化；针对li，t和具体为：
式中：li，t表示t时段内，第i个路段的观测值，i∈L；表示t时段内，第i个路段由不可观测OD需求构成的断面流量，i∈L；表示t时段内，第i个路段由可观测OD需求构成的断面流量，i∈L；L表示路网中全部可观测路段集合；表示t时段内，第i个进口道的观测值，i∈Lc；表示t时段内，第i个进口道的由不可观测部分构成的断面流量，i∈Lc；表示t时段内，第i个进口道的由可观测部分构成的断面流量，i∈Lc；Lc表示路网中可观测进口道集合；和均为检测误差；同时，
式中，表示在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达路段i的比例；表示在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达进口道i的比例；x(o，d)，t表示由交通小区o在t时段内出发，其目的地为交通小区d的实际需求量，其中o，d∈Q且o≠d，和分别为路段流量和交叉口进口道流量的观测误差。
本发明的有益效果是：
一、该种基于多源异构数据的动态车流OD估计方法，将号牌识别设备采集的过车数据与GPS数据相结合，提取出有效路径作为可观测OD需求，并划分为可观测部分和不可观测部分来实现路网OD估计，对比传统仅基于断面数据进行估计，该方法更能捕捉到OD变化趋势，更具准确性。
二、对比传统普通的卡尔曼滤波模型算法，本发明增加了状态约束条件，基于视频车牌是被设备在路网上的分布对部分交通小区流入路网的总需求进行观测，也有效提高了估计值的精度，能够较为准确的给出满足道路断面流量分布的OD需求动态估计值，结果更趋于实际情况。
附图说明
图1是本发明实施例基于多源异构数据的动态车流OD估计方法的流程示意图。
图2是实施例中交通小区划分情况以及可观测到的路段情况的示意图。
图3是采用标准卡尔曼滤波模型(KF模型)和实施例采用带约束的卡尔曼滤波模型(CKF模型)对不可观测OD需求估计得到的结果进行评估的对比示意图。
图4是将近基于断面数据得到的估计(单源数据)与实施例方法(多源数据)得到的估计进行对比的示意图。
具体实施方式
下面结合附图详细说明本发明的优选实施例。
实施例
一种基于多源异构数据的动态车流OD估计方法，将视频号牌识别数据和车辆GPS数据进行整合，将路网断面流量拆分为可观测OD和不可观测OD，同时基于卡尔曼滤波对不可观测OD进行估计，进而与可观测OD整合得到动态车流OD信息，如图1，具体步骤如下：
S1.基于用地性质、出行热点、视频号牌识别设备位置及路网拓扑结构将路网划分为若干个路网级的交通小区。
S2.基于视频号牌识别数据和车辆GPS数据提取出车辆轨迹信息。
S3.对缺失的轨迹数据进行补充和重构，根据车辆的完整行驶轨迹以及实际交通流量实现对路网断面流量观测值的拆分。
S31.基于最短路算法对由视频车牌识别数据获得的车流OD缺失部分进行补全，同时基于GPS定位数据提取出车辆完整的行驶轨迹。
S32.通过车辆完整的行驶轨迹信息，获得其在路网上的传播情况，即其经过路网上检测断面的实际时间，从而获得相应断面处的断面流量和其中表示t时段内，第i个进口道的由可观测部分构成的断面流量，i∈Lc；表示t时段内，第i个路段由可观测OD需求构成的断面流量，i∈L；L表示路网中全部可观测路段集合；Lc表示路网中可观测进口道集合。
S33.基于实际的断面交通流量对流量的观测值进行拆分，即将实际断面交通流量减去可观测的断面流量数值用于估计OD需求不可观测部分的断面流量和其中表示t时段内，第i个路段由不可观测OD需求构成的断面流量，i∈L；表示t时段内，第i个进口道的由不可观测部分构成的断面流量，i∈Lc。
S4.基于动态交通分配理论建立OD流量与断面流量之间的动态映射关系，建立分配矩阵。
S41.基于Logit模型建立路径选择模型，将行程时间作为效益函数，即OD出行中选择各路径的比例如下：
其中，为在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s的比例，其中o，d∈Q且o≠d，s∈R(o，d)，t；ti为有效路径集合R(o，d)，t中，车辆在第i条路径上的行程时间；为有效路径集的平均行程时间；R(o，d)，t表示在t时段内由交通小区o出发且目的地为交通小区d的出行路径集合，其中o，d∈Q且o≠d。
一般情况下，对于OD点对之间的有效路径集R(o，d)，t可以从已获得的部分车辆路径信息确定，具体方法如采用k最短路算法计算，或在最短路行程时间TTmin的基础上，适当乘以一个权重系数ρ(ρ≥1)，从而得到OD点对之间的行程时间范围[TTmin,ρTTmin]，再得到有效路径集，其中权责系数ρ根据所提取的可观测OD需求的实际行程时间和最短算法得到的路径行程时间确定。
S42.针对车流在具体路径上的传播过程得到路径流量与交叉口进口道断面流量之间的动态映射关系数值。
具体来说，将具体某一时段内选择某一路径的OD需求抽象成若干车辆构成的车辆元胞进行处理，则在t时段内由交通小区o出发到达交通小区d，选择路径s并在h时段内到达路段c的需求量占路径s流量的比例如下：
若(t-1)Δ≤γl＜γf＜tΔ，则
若γl＜(t-1)Δ≤γf＜tΔ，则
若(t-1)Δ≤γl＜tΔ≤γf＜(t+1)Δ，则
式中，表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s并在h时段内到交叉口进口道c的需求量占路径s流量的比例，其中o，d∈Q且o≠d，b∈s，s∈R(o，d)，t，h＝t，t+1，…t+δ-1，δ为完成出行的时段数；γf表示车辆元胞中第一辆车和最后一辆车到达某一路段或某一交叉口进口道的时刻；γl表示车辆元胞到达相应路段直至其完全消失的动态过程，且γf＞γl；Δ为每个时段的固定时长。
S43.得到某一具体路段或交叉口进口道处的断面流量观测值与待估计的OD需求之间的映射关系，即：
式中，为在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达路段i的比例；为在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达进口道i的比例。
S5.基于卡尔曼滤波对不可观测OD进行估计，进一步与可观测OD整合，对路网的动态需求进行求解。
S51.基于步骤S3拆分得到可观测的OD需求和不可观测的OD需求，在带约束的卡尔曼滤波构建动态OD估计模型基础下对不可观测的OD需求进行估计，建立出状态转移方程、状态线性约束方程和观测方程。
S511.基于随机漫步建立不可观测OD需求的状态转移方程，即：
式中：t表示第t个时段，t＝1，2，…，K；K为时段数；ωt是状态转移偏差，这里假设其为期望为0的白噪声；由交通小区o在t时段内出发，其目的地为交通小区d的待估计需求量，其中o，d∈Q且o≠d。
进一步向量化得到新的状态转移方程，即：
式中：为t时段内路网不可观测OD需求的列向量，为m×1维向量组；m为待估计的OD点对的个数；Wt为m×1维向量组。
S512.基于视频车牌识别设备在路网上的分布对部分交通小区的流入路网总需求进行观测，并根据可观测小区的需求产生量构建状态线性约束条件方程。
具体来说，假设对某一交通小区在某一时间段内流入路网的总需求的观测是确定的，不存在检测误差，则其与各OD点对之间的交通需求满足：
式中为t时段内，由交通小区o出发的不可观测部分的总需求量，o∈Q。
进一步向量化得到状态变量的线性约束条件，即：
其中为s×1维列向量，s为可进行需求全样本观测的交通小区的个数；Dt为s×m维矩阵。
S513.建立待估计OD需求与实际路网中采集到的断面流量之间映射关系的观测方程。
具体来说，首先基于视频车检器得到交通小区进入路网的实际交通需求流量作为状态变量的约束条件，建立状态变量线性约束，即：
式中，和分别为路段流量和交叉口进口道流量的观测误差。
而对于断面流量的OD量应包含可观测OD需求构成的部分和不可观测OD需求构成的部分，即：
式中，和均为检测误差。
进一步充分考虑车流在路网上的传播行为，针对每一个可观测的断面流量数据，假设某一时段的OD需求需经历δ个时段才能从路网上完全消失—即某一具体的断面流量由之前δ个时段的OD需求所构成。其中，和分别为路段流量和交叉口进口道流量的观测误差，且均为白噪声；参数δ需要基于所研究路网的实际运行状况确定。通过向量化处理，得到向量化观测方程，即：
其中，μt分别为n×1维向量组，n为可观测路段流量的个数；分别为θ×1维向量组，θ为可观测的交叉口进口道个数；Ph为n×m维矩阵，m为待估计的OD点对数量；为θ×m矩阵；
对上述两式左右两边变量进行合并，进一步化简后的观测方程如下：
其中Yt为(n+θ)×1维列向量，Ah为(n+θ)×m维矩阵，Ut为(n+θ)×1维列向量，
S52.将可观测OD需求和估计出的不可观测OD需求结合，构建出基于多源数据的动态OD估计模型，从而求解出路网动态OD估计值，具体公式如下：
式中：x(o，d)，t表示由交通小区o在t时段内出发，其目的地为交通小区d的实际需求量，其中o，d∈Q且o≠d；表示由交通小区o在t时段内出发，其目的地为交通小区d的观测需求量，其中o，d∈Q且o≠d；表示由交通小区o在t时段内出发，其目的地为交通小区d的待估计需求量，其中o，d∈Q且o≠d。
该种基于多源异构数据构建了动态车流OD估计方法，通过将OD需求拆分为可观测的确定部分和不可观测的待估计部分，其中针对不可观测部分基于状态空间模型构建了带线性状态约束条件的卡尔曼滤波估计模型并求解，并与可观测OD需求整合得到路网的动态OD估计值，为城市开放式路网提供了OD需求估计方法，可以得到更加准确，贴合实际情况的OD需求，从而为中微观交通模型及仿真平台提供合理化准确性的数据参数，提高了交通规划和管理效率。
实施例的一个具体示例如下：
选取某市某区域为研究对象，将其划分为35个交通小区，同时提取出路网号牌识别设备采集的过车信息和区域内出租车车辆GPS定位数据，提取出出租车车辆的运行轨迹及其经过每个号牌识别设备的信息。
根据实施例所述的步骤S3，基于GPS定位数据和号牌识别记录将出行轨迹补全，并将出租车车辆经过视频号牌设备的出行轨迹作为可观测OD，将剩下的交通流数据作为不可观测(待估计)OD。
交通小区划分情况以及可观测到的路段情况如图2所示：
进一步基于实施例所述的步骤S4，建立出OD需求与各断面流量之间的映射关系，并对于所述的步骤S5，基于建立的带约束卡尔曼滤波估计模型，对不可观测的OD需求进行估计，汇总得到路网OD估计值。
为评判步骤S5模型的准确性，本实施例对比采用标准卡尔曼滤波模型(KF模型)和S5中带约束的卡尔曼滤波模型(CKF模型)对不可观测OD需求估计得到的结果进行评估，具体对比结果如图3所示(选择某一小区，对比预测结果和实际结果)，发现带约束的卡尔曼滤波模型更贴近真实值。
同时将近基于断面数据得到的估计与实施例所述方法得到的估计进行对比，对比情况如图4所示(选择某一小区，对比预测结果和实际结果)，发现实施例所述模型更贴近实际数值。

Claims (9)

1.一种基于多源异构数据的动态车流OD估计方法，其特征在于：将视频号牌识别数据和车辆GPS数据进行整合，将路网断面流量拆分为可观测OD和不可观测OD，同时基于卡尔曼滤波对不可观测OD进行估计，进而与可观测OD整合得到动态车流OD信息，包括以下步骤，
S1、基于用地性质、出行热点、视频号牌识别设备位置及路网拓扑结构将路网划分为若干个路网级的交通小区；
S2、基于视频号牌识别数据和车辆GPS数据提取出各交通小区的车辆轨迹信息；
S3、对缺失的轨迹数据进行补充和重构，根据车辆的完整行驶轨迹以及实际交通流量实现对路网断面流量观测值的拆分为可观测OD和不可观测OD；
S4、基于动态交通分配理论建立OD流量与断面流量之间的动态映射关系，建立OD分配矩阵；
S5、基于卡尔曼滤波对不可观测OD进行估计，进一步与可观测OD整合，对路网的动态需求进行求解。
2.如权利要求1所述的基于多源异构数据的动态车流OD估计方法，其特征在于：步骤S3具体为，
S31、基于最短路算法对由视频车牌识别数据获得的车流OD缺失部分进行补全，同时基于GPS定位数据提取出车辆完整的行驶轨迹；
S32、通过车辆完整的行驶轨迹信息，获得其在路网上的传播情况，即其经过路网上检测断面的实际时间，从而获得相应断面处的断面流量和其中表示t时段内，第i个进口道的由可观测部分构成的断面流量，i∈Lc，Lc表示路网中可观测进口道集合；表示t时段内，第i个路段由可观测OD需求构成的断面流量，i∈L，L表示路网中全部可观测路段集合；
S33、基于实际的断面交通流量对流量的观测值进行拆分，即将实际断面交通流量减去可观测的断面流量数值用于估计OD需求不可观测部分的断面流量和其中表示t时段内，第i个路段由不可观测OD需求构成的断面流量，i∈L；表示t时段内，第i个进口道的由不可观测部分构成的断面流量，i∈Lc，Lc表示路网中可观测进口道集合。
3.如权利要求1所述的基于多源异构数据的动态车流OD估计方法，其特征在于：步骤S4具体为，
S41、基于Logit模型建立路径选择模型，将行程时间作为效益函数，即OD出行中选择各路径的比例如下：
其中，为在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s的比例，其中o，d∈Q，Q为交通小区集合；且o≠d，s∈R(o，d)，t；ts表示有效路径集合R(o，d)，t中，车辆行驶的行程时间；ti为有效路径集合R(o，d)，t中，车辆在第i条路径上的行程时间；为有效路径集的平均行程时间；R(o，d)，t表示在t时段内由交通小区o出发且目的地为交通小区d的出行路径集合，其中o，d∈Q且o≠d；
S42、针对车流在具体路径上的传播过程得到路径流量与交叉口进口道断面流量之间的动态映射关系数值；
S43、得到某一具体路段或交叉口进口道处的断面流量观测值与待估计的OD需求之间的映射关系，从而建立OD分配矩阵，其中OD需求之间的映射关系为：
式中，为在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达路段i的比例；为在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达进口道i的比例；表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s并在h时段内到达路段b的需求量占路径s流量的比例；表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s并在h时段内到交叉口进口道c的需求量占路径s流量的比例；表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s的比例，其中o，d∈Q且o≠d，b∈s，s∈R(o，d)，t，h＝t，t+1，…t+δ-1，δ为完成出行的时段数。
4.如权利要求3所述的基于多源异构数据的动态车流OD估计方法，其特征在于：步骤S41中，对于OD点对之间的有效路径集R(o，d)，t从已获得的部分车辆路径信息确定，具体采用k最短路算法计算，或在最短路行程时间TTmin的基础上，适当乘以一个权重系数ρ，ρ≥1，从而得到OD点对之间的行程时间范围[TTmin，ρTTmin]，再得到有效路径集，其中权责系数ρ根据所提取的可观测OD需求的实际行程时间和最短算法得到的路径行程时间确定。
5.如权利要求3所述的基于多源异构数据的动态车流OD估计方法，其特征在于：步骤S42具体为，将具体某一时段内选择某一路径的OD需求其中为在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s的比例，x(o，d)，t表示由交通小区o在t时段内出发，其目的地为交通小区d的实际需求量，o，d∈Q，Q为交通小区集合；
抽象成若干车辆构成的车辆元胞进行处理，则在t时段内由交通小区o出发到达交通小区d，选择路径s并在h时段内到达路段c的需求量占路径s流量的比例如下：
若(t-1)Δ≤γl＜γf＜tΔ，则
若γl＜(t-1)Δ≤γf＜tΔ，则
若(t-1)Δ≤γl＜tΔ≤γf＜(t+1)Δ，则
式中，表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s并在h时段内到交叉口进口道c的需求量占路径s流量的比例，其中o，d∈Q且o≠d，b∈s，s∈R(o，d)，t，h＝t，t+1，…t+δ-1，δ为完成出行的时段数；γf表示车辆元胞中第一辆车和最后一辆车到达某一路段或某一交叉口进口道的时刻；γl表示车辆元胞到达相应路段直至其完全消失的动态过程，且γf＞γl；Δ为每个时段的固定时长。
6.如权利要求1-5任一项所述的基于多源异构数据的动态车流OD估计方法，其特征在于：步骤S5具体为，
S51、基于步骤S3拆分得到可观测的OD需求和不可观测的OD需求，在带约束的卡尔曼滤波构建动态OD估计模型基础下对不可观测的OD需求进行估计，建立出状态转移方程、状态线性约束方程和观测方程；
S52、将可观测OD需求和估计出的不可观测OD需求结合，构建出基于多源数据的动态OD估计模型，从而求解出路网动态OD估计值，具体公式如下：
式中，x(o，d)，t表示由交通小区o在t时段内出发，其目的地为交通小区d的实际需求量，其中o，d∈Q且o≠d；表示由交通小区o在t时段内出发，其目的地为交通小区d的观测需求量，其中o，d∈Q且o≠d；表示由交通小区o在t时段内出发，其目的地为交通小区d的待估计需求量，其中o，d∈Q且o≠d。
7.如权利要求1-5任一项所述的基于多源异构数据的动态车流OD估计方法，其特征在于：步骤S51中，建立出状态转移方程为，
式中，为t时段内路网不可观测OD需求的列向量，为m×1维向量组；m为待估计的OD点对的个数；Wt为m×1维向量组，即为状态转移偏差ωt的向量化；
S512、基于视频车牌识别设备在路网上的分布对部分交通小区的流入路网总需求进行观测，并根据可观测小区的需求产生量构建状态线性约束条件方程；
S513、建立待估计OD需求与实际路网中采集到的断面流量之间映射关系的观测方程。
8.如权利要求7所述的基于多源异构数据的动态车流OD估计方法，其特征在于：步骤S512中，建立状态变量的线性约束方程为，
其中，为s×1维列向量，表示的向量化数值，其中s为可进行需求全样本观测的交通小区的个数；Dt为s×m维矩阵，即为交通小区的OD分配矩阵，为t时段内路网不可观测OD需求的列向量，为m×1维向量组。
9.如权利要求7所述的基于多源异构数据的动态车流OD估计方法，其特征在于：步骤S513中，建立观测方程为，
其中，Yt为(n+θ)×1维列向量，是OD需求向量数和的简化合并，即Ah为(n+θ)×m维矩阵，即到达比例向量数Ph和的合并，Ut为(n+θ)×1维列向量，即检测误差μt和的合并，
其中和分别为n×1维向量组，是li，t和的向量化数值，即为：
其中，Ph为n×m维矩阵，是的向量化，m为待估计的OD点对数量；为θ×m矩阵，是的向量化；μt和分别为θ×1维向量组，是和的向量化，为x(o，d)，h的向量化；针对li，t和具体为：
式中：li，t表示t时段内，第i个路段的观测值，i∈L；表示t时段内，第i个路段由不可观测OD需求构成的断面流量，i∈L；表示t时段内，第i个路段由可观测OD需求构成的断面流量，i∈L；L表示路网中全部可观测路段集合；表示t时段内，第i个进口道的观测值，i∈Lc；表示t时段内，第i个进口道的由不可观测部分构成的断面流量，i∈Lc；表示t时段内，第i个进口道的由可观测部分构成的断面流量，i∈Lc；Lc表示路网中可观测进口道集合；和均为检测误差；同时，
式中，表示在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达路段i的比例；表示在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达进口道i的比例；x(o，d)，t表示由交通小区o在t时段内出发，其目的地为交通小区d的实际需求量，其中o，d∈Q且o≠d，和分别为路段流量和交叉口进口道流量的观测误差。
CN201811083739.8A 2018-09-17 2018-09-17 基于多源异构数据的动态车流od估计方法 Active CN109035784B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811083739.8A CN109035784B (zh) 2018-09-17 2018-09-17 基于多源异构数据的动态车流od估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811083739.8A CN109035784B (zh) 2018-09-17 2018-09-17 基于多源异构数据的动态车流od估计方法

Publications (2)

Publication Number Publication Date
CN109035784A true CN109035784A (zh) 2018-12-18
CN109035784B CN109035784B (zh) 2021-01-26

Family

ID=64622531

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811083739.8A Active CN109035784B (zh) 2018-09-17 2018-09-17 基于多源异构数据的动态车流od估计方法

Country Status (1)

Country Link
CN (1) CN109035784B (zh)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110532515A (zh) * 2019-08-05 2019-12-03 北京交通大学 基于afc与视频数据的城市轨道交通乘客行程反演系统
CN111105613A (zh) * 2019-12-02 2020-05-05 北京建筑大学 一种基于多源数据的交通分配方法及系统
CN111145548A (zh) * 2019-12-27 2020-05-12 银江股份有限公司 一种基于数据场和节点压缩的重点路口识别及子区划分方法
CN111724589A (zh) * 2020-06-03 2020-09-29 重庆大学 一种基于多源数据的高速公路断面流量估计方法
WO2020207291A1 (zh) * 2019-04-10 2020-10-15 阿里巴巴集团控股有限公司 一种轨迹的处理方法、装置、及电子设备
CN111815944A (zh) * 2019-04-11 2020-10-23 阿里巴巴集团控股有限公司 数据有效性检测方法、装置、电子设备及计算机存储介质
CN111915903A (zh) * 2020-06-04 2020-11-10 东南大学 一种基于多源数据的快速路od矩阵获取方法
CN111915887A (zh) * 2020-07-10 2020-11-10 广州运星科技有限公司 一种基于多源异构交通数据的集成与处理系统及其方法
CN113034903A (zh) * 2021-03-05 2021-06-25 交通运输部公路科学研究所 基于多源信息融合的交通状态估计方法和装置
CN113450572A (zh) * 2021-09-02 2021-09-28 西南交通大学 一种交通状态估计方法、装置、设备及可读存储介质
CN113947922A (zh) * 2021-09-23 2022-01-18 重庆理工大学 基于网络切分的路网精细化动态od流量估计方法
CN114999162A (zh) * 2022-08-02 2022-09-02 北京交研智慧科技有限公司 一种道路交通流量获取方法和装置
CN116030617A (zh) * 2021-10-22 2023-04-28 交通运输部规划研究院 一种基于公路od数据进行车流量预测的方法及装置

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101373559A (zh) * 2007-08-24 2009-02-25 同济大学 基于浮动车数据评估城市路网交通状态的方法
CN102289932A (zh) * 2011-06-17 2011-12-21 同济大学 基于车辆自动识别设备的动态od矩阵估计方法
CN102592448A (zh) * 2012-01-18 2012-07-18 河海大学 一种利用当量交通冲突测评无信号交叉口安全性能的方法
CN103854483A (zh) * 2014-03-21 2014-06-11 广东新快易通智能信息发展有限公司 基于gps卫星定位及ecu的车辆动态信息系统
CN104537883A (zh) * 2014-12-18 2015-04-22 招商局重庆交通科研设计院有限公司 基于云服务平台的错位时间停车方法
CN105788252A (zh) * 2016-03-22 2016-07-20 连云港杰瑞电子有限公司 基于定点检测器和信号配时数据融合的城市干道车辆轨迹重构方法
CN105976613A (zh) * 2016-05-30 2016-09-28 西华大学 基于多项式趋势滤波的od矩阵估计方法
US9564048B2 (en) * 2014-12-18 2017-02-07 Sap Se Origin destination estimation based on vehicle trajectory data
CN106875314A (zh) * 2017-01-31 2017-06-20 东南大学 一种城市轨道交通线网客流od动态估计方法
CN107110650A (zh) * 2014-10-29 2017-08-29 赛峰电子与防务公司 在能观测性方面受约束的导航状态的估计方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101373559A (zh) * 2007-08-24 2009-02-25 同济大学 基于浮动车数据评估城市路网交通状态的方法
CN102289932A (zh) * 2011-06-17 2011-12-21 同济大学 基于车辆自动识别设备的动态od矩阵估计方法
CN102592448A (zh) * 2012-01-18 2012-07-18 河海大学 一种利用当量交通冲突测评无信号交叉口安全性能的方法
CN103854483A (zh) * 2014-03-21 2014-06-11 广东新快易通智能信息发展有限公司 基于gps卫星定位及ecu的车辆动态信息系统
CN107110650A (zh) * 2014-10-29 2017-08-29 赛峰电子与防务公司 在能观测性方面受约束的导航状态的估计方法
CN104537883A (zh) * 2014-12-18 2015-04-22 招商局重庆交通科研设计院有限公司 基于云服务平台的错位时间停车方法
US9564048B2 (en) * 2014-12-18 2017-02-07 Sap Se Origin destination estimation based on vehicle trajectory data
CN105788252A (zh) * 2016-03-22 2016-07-20 连云港杰瑞电子有限公司 基于定点检测器和信号配时数据融合的城市干道车辆轨迹重构方法
CN105976613A (zh) * 2016-05-30 2016-09-28 西华大学 基于多项式趋势滤波的od矩阵估计方法
CN106875314A (zh) * 2017-01-31 2017-06-20 东南大学 一种城市轨道交通线网客流od动态估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DIETMAR BAUER: "Quasi-Dynamic Estimation of OD Flows From Traffic Counts Without Prior OD Matrix", 《IEEE TRANSACTIONS ON INTELLIGENT TRANSPORTATION SYSTEMS》 *
聂庆慧: "基于交通传播特性分析的城市路网动态OD在线估计", 《中国博士学位论文全文数据库》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020207291A1 (zh) * 2019-04-10 2020-10-15 阿里巴巴集团控股有限公司 一种轨迹的处理方法、装置、及电子设备
CN111815944A (zh) * 2019-04-11 2020-10-23 阿里巴巴集团控股有限公司 数据有效性检测方法、装置、电子设备及计算机存储介质
CN111815944B (zh) * 2019-04-11 2022-08-02 阿里巴巴集团控股有限公司 数据有效性检测方法、装置、电子设备及计算机存储介质
CN110532515B (zh) * 2019-08-05 2021-08-03 北京交通大学 基于afc与视频数据的城市轨道交通乘客行程反演系统
CN110532515A (zh) * 2019-08-05 2019-12-03 北京交通大学 基于afc与视频数据的城市轨道交通乘客行程反演系统
CN111105613B (zh) * 2019-12-02 2021-01-26 北京建筑大学 一种基于多源数据的交通分配方法及系统
CN111105613A (zh) * 2019-12-02 2020-05-05 北京建筑大学 一种基于多源数据的交通分配方法及系统
CN111145548B (zh) * 2019-12-27 2021-06-01 银江股份有限公司 一种基于数据场和节点压缩的重点路口识别及子区划分方法
CN111145548A (zh) * 2019-12-27 2020-05-12 银江股份有限公司 一种基于数据场和节点压缩的重点路口识别及子区划分方法
CN111724589A (zh) * 2020-06-03 2020-09-29 重庆大学 一种基于多源数据的高速公路断面流量估计方法
CN111915903A (zh) * 2020-06-04 2020-11-10 东南大学 一种基于多源数据的快速路od矩阵获取方法
CN111915903B (zh) * 2020-06-04 2022-08-09 东南大学 一种基于多源数据的快速路od矩阵获取方法
CN111915887A (zh) * 2020-07-10 2020-11-10 广州运星科技有限公司 一种基于多源异构交通数据的集成与处理系统及其方法
CN113034903A (zh) * 2021-03-05 2021-06-25 交通运输部公路科学研究所 基于多源信息融合的交通状态估计方法和装置
CN113034903B (zh) * 2021-03-05 2021-11-16 交通运输部公路科学研究所 基于多源信息融合的交通状态估计方法和装置
CN113450572A (zh) * 2021-09-02 2021-09-28 西南交通大学 一种交通状态估计方法、装置、设备及可读存储介质
CN113947922A (zh) * 2021-09-23 2022-01-18 重庆理工大学 基于网络切分的路网精细化动态od流量估计方法
CN116030617A (zh) * 2021-10-22 2023-04-28 交通运输部规划研究院 一种基于公路od数据进行车流量预测的方法及装置
CN114999162A (zh) * 2022-08-02 2022-09-02 北京交研智慧科技有限公司 一种道路交通流量获取方法和装置

Also Published As

Publication number Publication date
CN109035784B (zh) 2021-01-26

Similar Documents

Publication Publication Date Title
CN109035784B (zh) 基于多源异构数据的动态车流od估计方法
CN110160550B (zh) 一种基于道路积水预测的城市路线引导方法
CN107817000A (zh) 无人驾驶车辆的路径规划方法、装置及计算机设备
CN110570660A (zh) 一种实时在线交通仿真系统及方法
US9240124B2 (en) Traffic-volume prediction device and method
CN104021674B (zh) 一种快速准确预测车辆通过路段旅行时间的方法
CN109686091B (zh) 一种基于多源数据融合的交通流量填补算法
CN113724489B (zh) 基于多源数据的交通拥堵溯源方法
CN101964941A (zh) 基于动态信息的智能导航与位置服务系统及方法
CN104615897B (zh) 基于低频gps数据的路段行程时间估算方法
CN105513359A (zh) 一种基于智能手机移动检测的城市快速路交通状态估计方法
CN112017436B (zh) 城市市内交通旅行时间的预测方法及系统
CN114861514A (zh) 一种车辆行驶方案的规划方法、装置和存储介质
Rahmani et al. Floating car and camera data fusion for non-parametric route travel time estimation
CN111275975B (zh) 一种交叉口转向流量数据的获取方法、装置和存储介质
CN104750963A (zh) 交叉口延误时长估计方法及装置
CN105279963B (zh) 一种获取路况信息的方法及装置
Ren et al. A hidden Markov model-based map-matching algorithm for wheelchair navigation
Trogh et al. Map matching and lane detection based on Markovian behavior, GIS, and IMU data
Grubwinkler et al. A modular and dynamic approach to predict the energy consumption of electric vehicles
JP2008181397A (ja) 交通流分岐率推定装置、交通流分岐率推定システムおよび交通流分岐率推定方法
CN112183871B (zh) 基于空气指数的城市交通诱导系统
CN111862657B (zh) 一种确定路况信息的方法及装置
CN107886192B (zh) 基于固定式与移动式车辆检测数据的数据与信息融合方法
CN113538902B (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
CB02 Change of applicant information

Address after: 211106 19 Su Yuan Avenue, Jiangning economic and Technological Development Zone, Nanjing, Jiangsu

Applicant after: JIANGSU ZHITONG TRAFFIC TECHNOLOGY Co.,Ltd.

Address before: 210006, Qinhuai District, Jiangsu, Nanjing should be 388 days street, Chenguang 1865 Technology Creative Industry Park E10 building on the third floor

Applicant before: JIANGSU ZHITONG TRAFFIC TECHNOLOGY Co.,Ltd.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant