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

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

Info

Publication number
CN109035784B
CN109035784B CN201811083739.8A CN201811083739A CN109035784B CN 109035784 B CN109035784 B CN 109035784B CN 201811083739 A CN201811083739 A CN 201811083739A CN 109035784 B CN109035784 B CN 109035784B
Authority
CN
China
Prior art keywords
traffic
observable
section
flow
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.)
Active
Application number
CN201811083739.8A
Other languages
English (en)
Other versions
CN109035784A (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 Zhitong Traffic Technology Co ltd
Original Assignee
Jiangsu Zhitong Traffic Technology 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 Zhitong Traffic Technology Co ltd filed Critical Jiangsu Zhitong Traffic Technology 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

Images

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、通过车辆完整的行驶轨迹信息，获得其在路网上的传播情况，即其经过路网上检测断面的实际时间，从而获得相应断面处的断面流量
Figure BDA0001801490580000021
和
Figure BDA0001801490580000022
其中
Figure BDA0001801490580000023
表示t时段内，第i个进口道的由可观测部分构成的断面流量，i∈Lc，Lc表示路网中可观测进口道集合；
Figure BDA0001801490580000024
表示t时段内，第i个路段由可观测OD需求构成的断面流量，i∈L，L表示路网中全部可观测路段集合；
S33、基于实际的断面交通流量对流量的观测值进行拆分，即将实际断面交通流量减去可观测的断面流量数值用于估计OD需求不可观测部分的断面流量
Figure BDA0001801490580000031
和
Figure BDA0001801490580000032
其中
Figure BDA0001801490580000033
表示t时段内，第i个路段由不可观测OD需求构成的断面流量，i∈L；
Figure BDA0001801490580000034
表示t时段内，第i个进口道的由不可观测部分构成的断面流量，i∈Lc，Lc表示路网中可观测进口道集合。
进一步地，步骤S4具体为，
S41、基于Logit模型建立路径选择模型，将行程时间作为效益函数，即OD出行中选择各路径的比例如下：
Figure BDA0001801490580000035
其中，
Figure BDA0001801490580000036
为在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条路径上的行程时间；
Figure BDA0001801490580000039
为有效路径集的平均行程时间；R(o，d)，t表示在t时段内由交通小区o出发且目的地为交通小区d的出行路径集合，其中o，d∈Q且o≠d；
S42、针对车流在具体路径上的传播过程得到路径流量与交叉口进口道断面流量之间的动态映射关系
Figure BDA0001801490580000037
数值；
S43、得到某一具体路段或交叉口进口道处的断面流量观测值与待估计的OD需求之间的映射关系，从而建立OD分配矩阵，其中OD需求之间的映射关系为：
Figure BDA0001801490580000038
式中，
Figure BDA0001801490580000041
为在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达路段i的比例；
Figure BDA0001801490580000042
为在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达进口道i的比例；
Figure BDA0001801490580000043
表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s并在h时段内到达路段b的需求量占路径s流量的比例；
Figure BDA0001801490580000044
表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s并在h时段内到交叉口进口道c的需求量占路径s流量的比例；
Figure BDA0001801490580000045
表示在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需求
Figure BDA0001801490580000046
其中
Figure BDA0001801490580000047
为在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Δ，则
Figure BDA0001801490580000052
若γl＜(t-1)Δ≤γf＜tΔ，则
Figure BDA0001801490580000053
若(t-1)Δ≤γl＜tΔ≤γf＜(t+1)Δ，则
Figure BDA0001801490580000054
式中，
Figure BDA0001801490580000055
表示在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估计值，具体公式如下：
Figure BDA0001801490580000051
式中，x(o，d)，t表示由交通小区o在t时段内出发，其目的地为交通小区d的实际需求量，其中o，d∈Q且o≠d；
Figure BDA0001801490580000056
表示由交通小区o在t时段内出发，其目的地为交通小区d的观测需求量，其中o，d∈Q且o≠d；
Figure BDA0001801490580000057
表示由交通小区o在t时段内出发，其目的地为交通小区d的待估计需求量，其中o，d∈Q且o≠d。
进一步地，步骤S51中，建立出状态转移方程为，
Figure BDA0001801490580000061
式中，
Figure BDA0001801490580000065
为t时段内路网不可观测OD需求的列向量，为m×1维向量组；m为待估计的OD点对的个数；Wt为m×1维向量组，即为状态转移偏差ωt的向量化；
S512、基于视频车牌识别设备在路网上的分布对部分交通小区的流入路网总需求进行观测，并根据可观测小区的需求产生量构建状态线性约束条件方程；
S513、建立待估计OD需求与实际路网中采集到的断面流量之间映射关系的观测方程。
进一步地，步骤S512中，建立状态变量的线性约束方程为，
Figure BDA0001801490580000066
其中，
Figure BDA0001801490580000067
为s×1维列向量，表示
Figure BDA0001801490580000068
的向量化数值，其中
Figure BDA0001801490580000069
s为可进行需求全样本观测的交通小区的个数；Dt为s×m维矩阵，即为交通小区的OD分配矩阵，
Figure BDA00018014905800000610
为t时段内路网不可观测OD需求的列向量，为m×1维向量组。
进一步地，步骤S513中，建立观测方程为，
Figure BDA0001801490580000062
其中，Yt为(n+θ)×1维列向量，是OD需求向量数
Figure BDA00018014905800000611
和
Figure BDA00018014905800000612
的简化合并，即
Figure BDA00018014905800000613
Ah为(n+θ)×m维矩阵，即到达比例向量数Ph和
Figure BDA00018014905800000614
的合并，
Figure BDA00018014905800000615
Ut为(n+θ)×1维列向量，即检测误差μt和
Figure BDA00018014905800000616
的合并，
Figure BDA00018014905800000617
其中
Figure BDA00018014905800000618
和
Figure BDA00018014905800000619
分别为n×1维向量组，是li，t和
Figure BDA00018014905800000620
的向量化数值，即为：
Figure BDA0001801490580000063
Figure BDA0001801490580000064
其中，Ph为n×m维矩阵，是
Figure BDA0001801490580000073
的向量化，m为待估计的OD点对数量；
Figure BDA0001801490580000074
为θ×m矩阵，是
Figure BDA0001801490580000075
的向量化；μt和
Figure BDA0001801490580000076
分别为θ×1维向量组，是
Figure BDA0001801490580000077
和
Figure BDA0001801490580000078
的向量化，
Figure BDA0001801490580000079
为x(o，d)，h的向量化；针对li，t和
Figure BDA00018014905800000710
具体为：
Figure BDA00018014905800000711
式中：li，t表示t时段内，第i个路段的观测值，i∈L；
Figure BDA00018014905800000712
表示t时段内，第i个路段由不可观测OD需求构成的断面流量，i∈L；
Figure BDA00018014905800000713
表示t时段内，第i个路段由可观测OD需求构成的断面流量，i∈L；L表示路网中全部可观测路段集合；
Figure BDA00018014905800000714
表示t时段内，第i个进口道的观测值，i∈Lc；
Figure BDA00018014905800000715
表示t时段内，第i个进口道的由不可观测部分构成的断面流量，i∈Lc；
Figure BDA00018014905800000716
表示t时段内，第i个进口道的由可观测部分构成的断面流量，i∈Lc；Lc表示路网中可观测进口道集合；
Figure BDA00018014905800000717
和
Figure BDA00018014905800000718
均为检测误差；同时，
Figure BDA0001801490580000071
Figure BDA0001801490580000072
式中，
Figure BDA00018014905800000719
表示在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达路段i的比例；
Figure BDA00018014905800000720
表示在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达进口道i的比例；x(o，d)，t表示由交通小区o在t时段内出发，其目的地为交通小区d的实际需求量，其中o，d∈Q且o≠d，
Figure BDA00018014905800000721
和
Figure BDA00018014905800000722
分别为路段流量和交叉口进口道流量的观测误差。
本发明的有益效果是：
一、该种基于多源异构数据的动态车流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.通过车辆完整的行驶轨迹信息，获得其在路网上的传播情况，即其经过路网上检测断面的实际时间，从而获得相应断面处的断面流量
Figure BDA0001801490580000091
和
Figure BDA0001801490580000092
其中
Figure BDA0001801490580000093
表示t时段内，第i个进口道的由可观测部分构成的断面流量，i∈Lc；
Figure BDA0001801490580000094
表示t时段内，第i个路段由可观测OD需求构成的断面流量，i∈L；L表示路网中全部可观测路段集合；Lc表示路网中可观测进口道集合。
S33.基于实际的断面交通流量对流量的观测值进行拆分，即将实际断面交通流量减去可观测的断面流量数值用于估计OD需求不可观测部分的断面流量
Figure BDA0001801490580000095
和
Figure BDA0001801490580000096
其中
Figure BDA0001801490580000097
表示t时段内，第i个路段由不可观测OD需求构成的断面流量，i∈L；
Figure BDA0001801490580000098
表示t时段内，第i个进口道的由不可观测部分构成的断面流量，i∈Lc。
S4.基于动态交通分配理论建立OD流量与断面流量之间的动态映射关系，建立分配矩阵。
S41.基于Logit模型建立路径选择模型，将行程时间作为效益函数，即OD出行中选择各路径的比例如下：
Figure BDA0001801490580000099
其中，
Figure BDA00018014905800000910
为在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s的比例，其中o，d∈Q且o≠d，s∈R(o，d)，t；ti为有效路径集合R(o，d)，t中，车辆在第i条路径上的行程时间；
Figure BDA00018014905800000911
为有效路径集的平均行程时间；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.针对车流在具体路径上的传播过程得到路径流量与交叉口进口道断面流量之间的动态映射关系
Figure BDA0001801490580000101
数值。
具体来说，将具体某一时段内选择某一路径的OD需求
Figure BDA0001801490580000102
抽象成若干车辆构成的车辆元胞进行处理，则在t时段内由交通小区o出发到达交通小区d，选择路径s并在h时段内到达路段c的需求量占路径s流量的比例如下：
若(t-1)Δ≤γl＜γf＜tΔ，则
Figure BDA0001801490580000103
若γl＜(t-1)Δ≤γf＜tΔ，则
Figure BDA0001801490580000104
若(t-1)Δ≤γl＜tΔ≤γf＜(t+1)Δ，则
Figure BDA0001801490580000105
式中，
Figure BDA0001801490580000106
表示在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需求之间的映射关系，即：
Figure BDA0001801490580000111
式中，
Figure BDA0001801490580000112
为在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达路段i的比例；
Figure BDA0001801490580000113
为在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达进口道i的比例。
S5.基于卡尔曼滤波对不可观测OD进行估计，进一步与可观测OD整合，对路网的动态需求进行求解。
S51.基于步骤S3拆分得到可观测的OD需求和不可观测的OD需求，在带约束的卡尔曼滤波构建动态OD估计模型基础下对不可观测的OD需求进行估计，建立出状态转移方程、状态线性约束方程和观测方程。
S511.基于随机漫步建立不可观测OD需求的状态转移方程，即：
Figure BDA0001801490580000114
式中：t表示第t个时段，t＝1，2，…，K；K为时段数；ωt是状态转移偏差，这里假设其为期望为0的白噪声；
Figure BDA0001801490580000115
由交通小区o在t时段内出发，其目的地为交通小区d的待估计需求量，其中o，d∈Q且o≠d。
进一步向量化得到新的状态转移方程，即：
Figure BDA0001801490580000116
式中：
Figure BDA0001801490580000117
为t时段内路网不可观测OD需求的列向量，为m×1维向量组；m为待估计的OD点对的个数；Wt为m×1维向量组。
S512.基于视频车牌识别设备在路网上的分布对部分交通小区的流入路网总需求进行观测，并根据可观测小区的需求产生量构建状态线性约束条件方程。
具体来说，假设对某一交通小区在某一时间段内流入路网的总需求
Figure BDA0001801490580000121
的观测是确定的，不存在检测误差，则其与各OD点对之间的交通需求满足：
Figure BDA0001801490580000122
式中
Figure BDA0001801490580000123
为t时段内，由交通小区o出发的不可观测部分的总需求量，o∈Q。
进一步向量化得到状态变量的线性约束条件，即：
Figure BDA0001801490580000124
其中
Figure BDA0001801490580000125
为s×1维列向量，s为可进行需求全样本观测的交通小区的个数；Dt为s×m维矩阵。
S513.建立待估计OD需求与实际路网中采集到的断面流量之间映射关系的观测方程。
具体来说，首先基于视频车检器得到交通小区进入路网的实际交通需求流量作为状态变量的约束条件，建立状态变量线性约束，即：
Figure BDA0001801490580000126
Figure BDA00018014905800001212
式中，
Figure BDA0001801490580000127
和
Figure BDA0001801490580000128
分别为路段流量和交叉口进口道流量的观测误差。
而对于断面流量的OD量应包含可观测OD需求构成的部分和不可观测OD需求构成的部分，即：
Figure BDA0001801490580000129
式中，
Figure BDA00018014905800001210
和
Figure BDA00018014905800001211
均为检测误差。
进一步充分考虑车流在路网上的传播行为，针对每一个可观测的断面流量数据，假设某一时段的OD需求需经历δ个时段才能从路网上完全消失—即某一具体的断面流量由之前δ个时段的OD需求所构成。其中，
Figure BDA0001801490580000134
和
Figure BDA0001801490580000135
分别为路段流量和交叉口进口道流量的观测误差，且均为白噪声；参数δ需要基于所研究路网的实际运行状况确定。通过向量化处理，得到向量化观测方程，即：
Figure BDA0001801490580000131
Figure BDA0001801490580000132
其中，
Figure BDA0001801490580000136
μt分别为n×1维向量组，n为可观测路段流量的个数；
Figure BDA0001801490580000137
分别为θ×1维向量组，θ为可观测的交叉口进口道个数；Ph为n×m维矩阵，m为待估计的OD点对数量；
Figure BDA0001801490580000138
为θ×m矩阵；
对上述两式左右两边变量进行合并，进一步化简后的观测方程如下：
Figure BDA0001801490580000133
其中Yt为(n+θ)×1维列向量，
Figure BDA0001801490580000139
Ah为(n+θ)×m维矩阵，
Figure BDA00018014905800001310
Ut为(n+θ)×1维列向量，
Figure BDA00018014905800001311
S52.将可观测OD需求和估计出的不可观测OD需求结合，构建出基于多源数据的动态OD估计模型，从而求解出路网动态OD估计值，具体公式如下：
Figure BDA00018014905800001312
式中：x(o，d)，t表示由交通小区o在t时段内出发，其目的地为交通小区d的实际需求量，其中o，d∈Q且o≠d；
Figure BDA00018014905800001313
表示由交通小区o在t时段内出发，其目的地为交通小区d的观测需求量，其中o，d∈Q且o≠d；
Figure BDA00018014905800001314
表示由交通小区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 (8)

1.一种基于多源异构数据的动态车流OD估计方法，其特征在于：将视频号牌识别数据和车辆GPS数据进行整合，将路网断面流量拆分为可观测OD和不可观测OD，同时基于卡尔曼滤波对不可观测OD进行估计，进而与可观测OD整合得到动态车流OD信息，包括以下步骤，
S1、基于用地性质、出行热点、视频号牌识别设备位置及路网拓扑结构将路网划分为若干个路网级的交通小区；
S2、基于视频号牌识别数据和车辆GPS数据提取出各交通小区的车辆轨迹信息；
S3、对缺失的轨迹数据进行补充和重构，根据车辆的完整行驶轨迹以及实际交通流量实现对路网断面流量观测值的拆分为可观测OD和不可观测OD；
S4、基于动态交通分配理论建立OD流量与断面流量之间的动态映射关系，建立OD分配矩阵；
S5、基于卡尔曼滤波对不可观测OD进行估计，进一步与可观测OD整合，对路网的动态需求进行求解；步骤S5具体为，
S51、基于步骤S3拆分得到可观测的OD需求和不可观测的OD需求，在带约束的卡尔曼滤波构建动态OD估计模型基础下对不可观测的OD需求进行估计，建立出状态转移方程、状态线性约束方程和观测方程；
S52、将可观测OD需求和估计出的不可观测OD需求结合，构建出基于多源数据的动态OD估计模型，从而求解出路网动态OD估计值，具体公式如下：
Figure FDA0002816565860000011
式中，x(o,d),t表示由交通小区o在t时段内出发，其目的地为交通小区d的实际需求量，其中o,d∈Q且o≠d；
Figure FDA0002816565860000012
表示由交通小区o在t时段内出发，其目的地为交通小区d的观测需求量，其中o,d∈Q且o≠d；
Figure FDA0002816565860000013
表示由交通小区o在t时段内出发，其目的地为交通小区d的待估计需求量，其中o,d∈Q且o≠d。
2.如权利要求1所述的基于多源异构数据的动态车流OD估计方法，其特征在于：步骤S3具体为，
S31、基于最短路算法对由视频车牌识别数据获得的车流OD缺失部分进行补全，同时基于GPS定位数据提取出车辆完整的行驶轨迹；
S32、通过车辆完整的行驶轨迹信息，获得其在路网上的传播情况，即其经过路网上检测断面的实际时间，从而获得相应断面处的断面流量
Figure FDA0002816565860000021
和
Figure FDA0002816565860000022
其中
Figure FDA0002816565860000023
表示t时段内，第i个进口道的由可观测部分构成的断面流量，i∈Lc，Lc表示路网中可观测进口道集合；
Figure FDA0002816565860000024
表示t时段内，第i个路段由可观测OD需求构成的断面流量，i∈L，L表示路网中全部可观测路段集合；
S33、基于实际的断面交通流量对流量的观测值进行拆分，即将实际断面交通流量减去可观测的断面流量数值用于估计OD需求不可观测部分的断面流量
Figure FDA0002816565860000025
和
Figure FDA0002816565860000026
其中
Figure FDA0002816565860000027
表示t时段内，第i个路段由不可观测OD需求构成的断面流量，i∈L；
Figure FDA0002816565860000028
表示t时段内，第i个进口道的由不可观测部分构成的断面流量，i∈Lc，Lc表示路网中可观测进口道集合。
3.如权利要求1所述的基于多源异构数据的动态车流OD估计方法，其特征在于：步骤S4具体为，
S41、基于Logit模型建立路径选择模型，将行程时间作为效益函数，即OD出行中选择各路径的比例如下：
Figure FDA0002816565860000029
其中，
Figure FDA00028165658600000210
为在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条路径上的行程时间；
Figure FDA00028165658600000211
为有效路径集的平均行程时间；R(o,d),t表示在t时段内由交通小区o出发且目的地为交通小区d的出行路径集合，其中o,d∈Q且o≠d；
S42、针对车流在具体路径上的传播过程得到路径流量与交叉口进口道断面流量之间的动态映射关系
Figure FDA00028165658600000212
数值；
S43、得到某一具体路段或交叉口进口道处的断面流量观测值与待估计的OD需求之间的映射关系，从而建立OD分配矩阵，其中OD需求之间的映射关系为：
Figure FDA0002816565860000031
式中，
Figure FDA0002816565860000032
为在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达路段i的比例；
Figure FDA0002816565860000033
为在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达进口道i的比例；
Figure FDA0002816565860000034
表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s并在h时段内到达路段b的需求量占路径s流量的比例；
Figure FDA0002816565860000035
表示在t时段内由交通小区o出发且目的地为交通小区d，其选择路径s并在h时段内到交叉口进口道c的需求量占路径s流量的比例；
Figure FDA0002816565860000036
表示在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需求
Figure FDA0002816565860000037
其中
Figure FDA0002816565860000038
为在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Δ，则
Figure FDA0002816565860000041
若γl＜(t-1)Δ≤γf＜tΔ，则
Figure FDA0002816565860000042
若(t-1)Δ≤γl＜tΔ≤γf＜(t+1)Δ，则
Figure FDA0002816565860000043
式中，
Figure FDA0002816565860000044
表示在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估计方法，其特征在于：步骤S51中，建立出状态转移方程为，
Figure FDA0002816565860000045
式中，
Figure FDA0002816565860000046
为t时段内路网不可观测OD需求的列向量，为m×1维向量组；m为待估计的OD点对的个数；Wt为m×1维向量组，即为状态转移偏差ωt的向量化；
S512、基于视频车牌识别设备在路网上的分布对部分交通小区的流入路网总需求进行观测，并根据可观测小区的需求产生量构建状态线性约束条件方程；
S513、建立待估计OD需求与实际路网中采集到的断面流量之间映射关系的观测方程。
7.如权利要求6所述的基于多源异构数据的动态车流OD估计方法，其特征在于：步骤S512中，建立状态变量的线性约束方程为，
Figure FDA0002816565860000047
其中，
Figure FDA0002816565860000048
为s×1维列向量，表示
Figure FDA0002816565860000049
的向量化数值，其中
Figure FDA00028165658600000410
s为可进行需求全样本观测的交通小区的个数；Dt为s×m维矩阵，即为交通小区的OD分配矩阵，
Figure FDA00028165658600000411
为t时段内路网不可观测OD需求的列向量，为m×1维向量组。
8.如权利要求6所述的基于多源异构数据的动态车流OD估计方法，其特征在于：步骤S513中，建立观测方程为，
Figure FDA0002816565860000051
其中，Yt为(n+θ)×1维列向量，是OD需求向量数
Figure FDA0002816565860000052
和
Figure FDA0002816565860000053
的简化合并，即
Figure FDA0002816565860000054
Ah为(n+θ)×m维矩阵，即到达比例向量数Ph和
Figure FDA0002816565860000055
的合并，
Figure FDA0002816565860000056
Ut为(n+θ)×1维列向量，即检测误差μt和
Figure FDA0002816565860000057
的合并，
Figure FDA0002816565860000058
其中
Figure FDA0002816565860000059
和
Figure FDA00028165658600000510
分别为n×1维向量组，是li,t和
Figure FDA00028165658600000511
的向量化数值，即为：
Figure FDA00028165658600000512
Figure FDA00028165658600000513
其中，Ph为n×m维矩阵，是
Figure FDA00028165658600000514
的向量化，m为待估计的OD点对数量；
Figure FDA00028165658600000515
为θ×m矩阵，是
Figure FDA00028165658600000516
的向量化；μt和
Figure FDA00028165658600000517
分别为θ×1维向量组，是
Figure FDA00028165658600000518
和
Figure FDA00028165658600000519
的向量化，
Figure FDA00028165658600000520
为x(o,d),h的向量化；针对li,t和
Figure FDA00028165658600000521
具体为：
Figure FDA00028165658600000522
式中：li,t表示t时段内，第i个路段的观测值，i∈L；
Figure FDA00028165658600000523
表示t时段内，第i个路段由不可观测OD需求构成的断面流量，i∈L；
Figure FDA00028165658600000524
表示t时段内，第i个路段由可观测OD需求构成的断面流量，i∈L；L表示路网中全部可观测路段集合；
Figure FDA00028165658600000525
表示t时段内，第i个进口道的观测值，i∈Lc；
Figure FDA00028165658600000526
表示t时段内，第i个进口道的由不可观测部分构成的断面流量，i∈Lc；
Figure FDA00028165658600000527
表示t时段内，第i个进口道的由可观测部分构成的断面流量，i∈Lc；Lc表示路网中可观测进口道集合；
Figure FDA00028165658600000528
和
Figure FDA00028165658600000529
均为检测误差；同时，
Figure FDA00028165658600000530
Figure FDA0002816565860000061
式中，
Figure FDA0002816565860000062
表示在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达路段i的比例；
Figure FDA0002816565860000063
表示在t时段内由交通小区o出发且目的地为交通小区d，在h时段内到达进口道i的比例；x(o,d),t表示由交通小区o在t时段内出发，其目的地为交通小区d的实际需求量，其中o,d∈Q且o≠d，
Figure FDA0002816565860000064
和
Figure FDA0002816565860000065
分别为路段流量和交叉口进口道流量的观测误差。
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 CN109035784A (zh) 2018-12-18
CN109035784B true 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)

Families Citing this family (12)

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

Citations (8)

* 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 连云港杰瑞电子有限公司 基于定点检测器和信号配时数据融合的城市干道车辆轨迹重构方法
US9564048B2 (en) * 2014-12-18 2017-02-07 Sap Se Origin destination estimation based on vehicle trajectory data
CN107110650A (zh) * 2014-10-29 2017-08-29 赛峰电子与防务公司 在能观测性方面受约束的导航状态的估计方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105976613A (zh) * 2016-05-30 2016-09-28 西华大学 基于多项式趋势滤波的od矩阵估计方法
CN106875314B (zh) * 2017-01-31 2020-11-10 东南大学 一种城市轨道交通线网客流od动态估计方法

Patent Citations (8)

* 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 连云港杰瑞电子有限公司 基于定点检测器和信号配时数据融合的城市干道车辆轨迹重构方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Quasi-Dynamic Estimation of OD Flows From Traffic Counts Without Prior OD Matrix;Dietmar Bauer;《IEEE TRANSACTIONS ON INTELLIGENT TRANSPORTATION SYSTEMS》;20180630;正文摘要 *
基于交通传播特性分析的城市路网动态OD在线估计;聂庆慧;《中国博士学位论文全文数据库》;20171130;摘要、图1-1、1-4节 *

Also Published As

Publication number Publication date
CN109035784A (zh) 2018-12-18

Similar Documents

Publication Publication Date Title
CN109035784B (zh) 基于多源异构数据的动态车流od估计方法
CN110160550B (zh) 一种基于道路积水预测的城市路线引导方法
CN110570660A (zh) 一种实时在线交通仿真系统及方法
WO2018064931A1 (zh) 考虑出租车运营状态下的城市道路行程时间分布估算方法
US9240124B2 (en) Traffic-volume prediction device and method
CN107817000A (zh) 无人驾驶车辆的路径规划方法、装置及计算机设备
CN104021674B (zh) 一种快速准确预测车辆通过路段旅行时间的方法
CN104574967A (zh) 一种基于北斗的城市大面积路网交通感知方法
CN110335459A (zh) 低渗透率车辆轨迹数据的交叉口排队长度估计方法及装置
CN109686091B (zh) 一种基于多源数据融合的交通流量填补算法
CN105074793A (zh) 用于车辆路径规划和交通管理的车道级车辆导航
CN104567906A (zh) 一种基于北斗的城市路网车辆路径规划方法及装置
CN112017436B (zh) 城市市内交通旅行时间的预测方法及系统
CN111275975B (zh) 一种交叉口转向流量数据的获取方法、装置和存储介质
Rahmani et al. Floating car and camera data fusion for non-parametric route travel time estimation
CN110796876A (zh) 基于卡尔曼滤波的路段车辆总数估计方法
CN113724489A (zh) 基于多源数据的交通拥堵溯源方法
Ren et al. A hidden Markov model-based map-matching algorithm for wheelchair navigation
Chen et al. Prediction of dynamic freeway travel times based on vehicle trajectory construction
Grubwinkler et al. A modular and dynamic approach to predict the energy consumption of electric vehicles
Trogh et al. Map matching and lane detection based on Markovian behavior, GIS, and IMU data
CN114511938A (zh) 一种智慧停车管理方法和平台
Gan et al. Estimation of performance metrics at signalized intersections using loop detector data and probe travel times
JP2008181397A (ja) 交通流分岐率推定装置、交通流分岐率推定システムおよび交通流分岐率推定方法
CN114897306B (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