CN109791729B - 一种基于行程时间分布的交通异常检测方法 - Google Patents
一种基于行程时间分布的交通异常检测方法 Download PDFInfo
- Publication number
- CN109791729B CN109791729B CN201780050907.0A CN201780050907A CN109791729B CN 109791729 B CN109791729 B CN 109791729B CN 201780050907 A CN201780050907 A CN 201780050907A CN 109791729 B CN109791729 B CN 109791729B
- Authority
- CN
- China
- Prior art keywords
- time
- data
- traffic
- space
- real
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
-
- 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/0125—Traffic data processing
- G08G1/0129—Traffic data processing for creating historical data or processing based on historical data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
-
- 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
-
- 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/0108—Measuring and analyzing of parameters relative to traffic conditions based on the source of data
- G08G1/0112—Measuring and analyzing of parameters relative to traffic conditions based on the source of data from the vehicle, e.g. floating car data [FCD]
-
- 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/0125—Traffic data processing
- G08G1/0133—Traffic data processing for classifying traffic situation
-
- 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/0141—Measuring and analyzing of parameters relative to traffic conditions for specific applications for traffic information dissemination
Abstract
一种交通异常检测方法。利用浮动车的车载GNSS定位装置,可获取其不同时刻的空间位置信息,通过对海量浮动车轨迹信息的分析挖掘,可实现城市道路交通异常事件的智能化检测。该交通异常检测方法,采用行程时间的概率分布表征交通状态,采用概率分布差异衡量指标反映交通状态差异,与现有技术相比,具有原理明晰、实施简便、检测率高的特点。
Description
技术领域
本发明属于交通检测技术领域。特别地,本发明涉及一种城市道路交通异常实时检测方法。通过浮动车的车载GNSS定位装置,可获取其不同时刻的空间位置信息,经过数据预处理、地图匹配和数据融合,获得特定时空范围的行程时间概率分布;根据行程时间分布的变化情况,可有效识别城市道路交通异常事件。
背景技术
交通异常事件检测是城市交通管理的重要组成部分,也是智能交通系统的核心功能之一。交通异常事件主要包括交通事故、车辆抛锚、货车落物、道路交通设施损坏或故障以及其他造成交通流紊乱的特殊事件。该类事件容易造成交通拥堵、路段通行能力降低,严重时影响整个道路交通系统的正常运行。通过交通异常事件检测,可使交通管理者及时了解交通异常信息,并采取适当的诱导和控制措施,降低交通异常事件的不良影响。
交通异常事件检测可分为人工方式和自动方式。人工方式包括巡逻车、紧急电话上报和视频监控等,由于消耗人力物力且实时性差,无法满足交通管理的需要。自动方式依靠自动事件检测(AID,Automated Incidence Detection)算法实现,基本原理是通过检测不同位置道路交通流的变化来识别交通异常事件。目前常用的AID算法包括模式识别类算法(如California算法、莫妮卡算法)、统计预测类算法(如指数平滑法、卡尔曼滤波算法)、交通流模型算法(如McMaster算法)以及智能识别算法(如人工神经网络、模糊逻辑算法)。
但是目前的检测方法存在对设施的要求高、计算复杂度高、无法对异常状况的态势做进一步判断等缺点。本发明利用出租车、公交车车载GNSS定位装置回传的轨迹数据,建立历史交通状态数据库和实时交通状态数据库,通过分析两者反映的交通流特征差异,识别交通异常事件。该方法具有实时性好、可并行处理、识别率高以及对检测设施要求低等特点,适用于有实时浮动车定位数据的数据环境下城市道路交通异常事件的检测。
目前,针对交通异常事件监测,有以下代表性技术:
一件美国专利申请,US 20160148512,披露了一种交通异常事件检测和上报系统的组成原理和实施方法。该系统由传感器、通信模块、移动处理模块和用户交互模块组成。传感器用于采集车辆周边的相关数据;通信模块用于发送本车辆数据和接收周边车辆的数据;移动处理模块用于处理和分析相关车辆在某一区域内的数据并生成交通事件报告;用户交互模块能够像用户提供交通事件报告。该方案是一种基于车车和车路通讯网络的交通异常事件检测技术,能够利用传感器采集的各类信息,判别异常事件。然而,由于传感器、通信单元需要单独安装调试,实施难度较大;移动处理单元处理能力受限;同时需要移动和固定的讯息接收端,系统本身存在故障概率,可靠性不佳。
一件中国专利申请,CN 104809878 A,披露了一种利用公交车GPS数据检测城市道路交通异常状态的方法。该方案根据GPS历史数据获得路段延误时间指数,根据GPS当前数据获得瞬时速度、周期平均速度、加权滑动平均速度和多车平均速度,利用规范变量分析算法检测异常。这一方案不需要新增检测设施,实施便利。但是对于交通态势的表征过于简化,无法分析交通异常状况的特点和成因;对交通场景的划分缺乏依据,未能考虑天气等因素对交通态势变化的影响。
发明内容
为了更清晰地阐述本发明的内容,首先将涉及到的专业术语解释如下:
浮动车:也称探测车。指安装了车载定位装置并行驶在城市道路上的公交汽车和出租车。
GNSS:全球导航卫星系统(Global Navigation Satellite System)。包括GPS、GLONASS、GALILEO以及北斗卫星导航系统等。
时空子区:按照时间和空间两个维度划分的片区,反映在一段时间内,一定的空间范围内的情况。将一天划分为若干时间片段,例如0:00-0:10,0:10-0:20……,每个时间片段称之为一个时间子区;将城市道路交通异常检测的实施区域划分为若干空间片段,例如经度121.58°E-121.59°E,纬度31.16°N-31.17°N之间的区域,每个空间片段称为一个空间子区;任意一个时间子区和任意一个空间子区的交集形成的时空片段,称为时空子区,例如经度121.58°E-121.59°E,纬度31.16°N-31.17°N之间的区域在0:00-0:10的时空片段。
历史轨迹数据:历史轨迹数据是长时间积累并存储在数据库中的轨迹数据。历史轨迹数据是动态变化的数据,需要及时进行更新,并定期做重新处理和分析,以保证历史交通特征提取的准确性。每个时空子区的数据可以并行处理以提高效率。本发明中可简称为历史数据。
实时轨迹数据:实时轨迹数据是距离当前时刻最近的一个时间区段内的轨迹数据集合。本发明中可简称为实时数据。
交通态势:一定时间、一定空间内交通运行的综合情况的总称。
交通异常:交通事故、车辆抛锚、货车落物、道路交通设施损坏或故障等事件引发的交通流紊乱的情况。
交通异常严重性:即交通流紊乱的严重性,是正常状态下交通流与交通异常发生后交通流特征的差异。
交通异常指数:交通异常严重性的量度。范围为0~10,数值越大,交通异常越严重。
交通环境:作用于道路交通参与者的所有外界影响与力量的总和。包括道路状况、交通设施、地物地貌、气象条件,以及其他交通参与者的交通活动。
地图匹配:将地理坐标与城市路网关联的过程。
高峰小时流量:某城市道路断面一日内小时交通流量的最大值。
有限混合模型:一种用简单密度模拟复杂密度的数学方法。变量集合为y、成分数量为K的有限混合模型可表示为:
响应变量:根据自变量发生改变的变量,也称因变量。
贝叶斯信息准则:是在不完全情报下,对部分未知的状态用主观概率估计,然后用贝叶斯公式对发生概率进行修正的结果可靠性的评价指标。其计算方法为:
BIC=-2ln L+k·ln n
式中,L为似然函数的最大值,k为未知参数的个数,n为样本量。
似然函数:似然函数是一种关于统计模型参数的函数。给定输出x时,关于参数θ的似然函数L(θ|x)(在数值上)等于给定参数θ后变量X的概率:L(θ|x)=P(X=x|θ)。
参数估计:根据从总体中抽取的样本估计总体分布中包含的未知参数的方法。
EM算法:最大期望算法(Expectation Maximization Algorithm),是一种迭代算法,用于含有隐变量的概率参数模型的最大似然估计或极大后验概率估计。
Kullback-Leibler散度:两个概率分布P和Q差异的一种量度。
Jensen-Shannon散度:是Kullback-Leibler散度的一种对称化形式。
K-Medoids算法:一种聚类算法,每次迭代都从当前类别中选取这样一个点——它到其他所有(当前类别中的)点的距离之和最小——作为新的中心点。
本发明的目的是建立一套基于浮动车轨迹记录系统,利用历史GNSS定位数据和实时GNSS定位数据,结合交通环境信息识别道路交通异常事件的方案。为了达到上述目的,本发明提供了如下技术方案:
本发明的实施前提是:搭载GNSS轨迹记录仪的浮动车(出租车、公交车等);具有大规模存储、计算、实时任务处理能力的数据中心。
本发明的适用范围是:有上述浮动车经过的城市道路(包括地面道路和高架道路)。
本发明的实施步骤包括:
1)确定检测的时空范围和建立时空子区。
基于实际的应用需求,确定需要进行交通异常事件检测的时间范围和空间范围。时间范围可以设定为全天,即0:00-24:00;也可以设定为某一特定的时段,例如要检测17:00-20:00这个时段的交通异常时间,则将检测时间范围设定为17:00-20:00,这里只是列举一个特殊实例,还有很多其他情况,此处不再一一说明。空间范围可以按照行政区划设置为某个市域,例如北京市、上海市、黄浦区等;也可以按照城市空间结构设置为某个城市功能区,例如某市中央商务区、工业区等。
时空子区的建立是指,将检测的时间范围划分为若干个更小的时间片段,将检测的空间范围,即城市道路交通异常检测的实施区域,划分为若干个更小的空间片段。时空子区的建立,可以采用多种经验划分方法,包括等距时空划分法和非等距时空划分法。
2)数据预处理。
将GNSS定位数据进行数据清洗、数据集成、数据转换、数据归约,提高数据的结构化程度。GNSS,即全球导航卫星系统定位系统,是能在地球表面或近地空间的任何地点为提供全天候三维坐标和速度以及时间信息的空基无线电导航定位系统。它主要包括美国的GPS(Global Positioning System),俄罗斯的GLONASS(Global Navigation SatelliteSystem)、欧盟的GALILEO和中国的北斗卫星导航系统四大全球性导航定位系统,同时还包括日本的QZSS、印度的IRNSS等区域导航定位系统以及美国的WASS、日本的MSAS等卫星定位增强系统。为了在不同的导航定位系统设备中建立统一的数据分发标准,美国国家海洋电子协会制定了统一的NEMA(National Marine Electronics Association)通讯协议,以规范GNSS的数据广播。因此,GNSS中的各个成员系统,例如GPS、GLONASS等,虽然分别由不同国家和机构建立和维护,但是拥有一致的数据分发格式,因此不需要对数据格式进行变换。
选定的空间范围内,有许多安装GNSS定位设备的车辆,常见的有出租车、公交车、货运汽车、私家车等。基于当前城市交通数据应用现状,在实际应用当中,通常选用城市出租车为浮动车作为交通异常检测系统的数据来源。
采集的GNSS定位信息中包含一些不合理的信息,为了保证交通异常状态检测判别结果的准确性,首先需要进行甄别以提出异常的数据,保证数据的可靠性。这些异常数据包括:落在检测时空范围之外的数据、明显超出合理范围的空间位置跳跃。所谓“明显超出合理范围的空间位置跳跃”,下面举例说明之。若某日10:30:00时刻某辆浮动车定位设备上传的定位点记为A,当日10:30:30时刻该浮动车定位设备上传的定位点记为B,位置A与位置B的距离为1500米,那么据此计算得到该浮动车的行驶速度至少为180km/h,超出了一般常识,因此是一种异常的空间位置跳跃,数据处理中应当予以剔除。
3)快速地图匹配。
经过预处理后的GNSS定位数据,需要结合城市路网数据,通过地图匹配算法,将GNSS定位点投影到城市地图,建立定位点与路段的匹配关系,并修正定位漂移带来的误差。
目前各个地理区域的电子地图都已较为详实,这种电子地图可以来源于城市的地理信息系统,当然也可以来源自其他方式和途经。这些电子地图对城市道路信息进行了详细刻画,通过划分可以得到若干路段。通过借助距离、角度等信息,将定位点匹配到路段上,这样就实现了将定位信息匹配到实际的地理环境中。
4)浮动车路径的表示和不同车辆路径的匹配。
在给定一组起终点的前提下,车辆的路径可能不是唯一的。复杂的城市交通路网包含了若干路段,将这些不同的路段进行编号,例如,将路段表示为L1,L2等。道路可能有两个不同的行驶方向,在这种情况下,应该将两个不同的行驶方向表示为两个不同的路段,给予不同的路段编号。
给定的起点和终点,通常可采用城市路网中路段的交点。已知某浮动车行驶的路径,现需要从其他浮动车已经发送的路径信息中,选择与该浮动车路径相同的路径,从而获得起点和终点间的同路径组。
5)数据抽样。
浮动车的定位数据中,包含位置坐标、瞬时车速、记录时间等信息。在本专利提出的基于行程时间分布的交通异常检测方法中,数据抽样是指从全部的浮动车数据中筛选出部分数据进行后续的分析处理,这种筛选是基于数据中心的计算能力以及预先提出的精度要求而进行的。基于不同的计算能力和精度要求,可采用不同的数据抽样方法。例如,当数据中心的计算能力较强,且对检测的精度要求较高时,可以将全部的浮动车定位数据作为处理对象,进行全面的处理分析;而当数据中心的计算能力有限时,假定当前的数据中心能够在1分钟内,对每个空间子区处理500条数据,而实际情况是在1分钟每个空间子区能产生了2000条浮动车定位数据,那么可以从2000条数据中随机抽取500条数据进行分析,从而在数据中心的计算能力范围内,获得精度受限的处理结果。
根据对浮动车数据利用方式的不同,可以针对浮动车数据的不同属性进行采样,例如行程车速和行程时间等。本专利中提出的基于行程时间分布的交通异常检测方法,采用行程时间作为基础进行城市道路交通异常检测。因此,数据抽样是指对行程时间进行抽样。
6)历史轨迹数据分析和特征提取。
所谓历史轨迹数据,是指在长期的城市道路交通运行中积累下来的浮动车轨迹数据。利用历史浮动车轨迹数据,可以建立城市道路交通特征模型,用来反映城市交通运行的一般特性。这里所说的城市道路交通特征模型,可以指某些特定的指标,例如平均速度、加权平均速度等;也可以指各种某种统计模型,例如行程速度的概率分布。以往的很多模型,采用单一的指标来表示某个路段或区域的交通特征(如历史平均行程时间),这种方式虽然应用简便,但是精度不高,敏感性差,往往不能在交通异常状态检测中发挥良好的效果。因此,本专利提出对于每个时空子区,用交通特征变量的概率分布描述交通特征,建立交通特征模型并进行参数估计。
可采集的交通特征变量,包括行程车速和行程时间等,本专利中所述的交通特征变量的概率分布,指行程时间的概率分布。
7)实时轨迹数据分析和特征提取。
所谓实时轨迹数据,是指距离当前时刻不远的一段时间内的交通运行中浮动车的轨迹数据。利用实时浮动车轨迹数据,可以掌握交通特征的变化动态,用来反映当前交通运行的即时特性。本专利采用当前时空子区特定起止点间的行程时间描述当前交通特征。
8)异常检测。
系统状态异常检测的思想最早由Denning提出,即通过监视系统审计记录上系统使用的异常情况,可以检测出违反安全、可能引发系统异常的事件。Denning建立的这种模型独立于任何特定的系统、应用环境、系统弱点、故障类型,因而是一种普遍意义上的异常检测模型。该模型包括主体、客体、审计记录、轮廓、异常记录和活动规则5个部分。轮廓是用度量和统计模型来表示的主体相对于客体的正常行为。Denning的模型定义了3种度量,即事件计数器、间隔定时器、资源测量器,并提出了5种统计模型,即可操作模型、均值和标准差模型、多变量模型、马尔可夫过程模型和时间序列模型。Denning提出的模型通过对系统审计数据的分析,建立起系统主体的基于统计的正常行为特征轮廓,检测时,系统中的审计数据与已建立的主体的正常行为特征轮廓相比较,若相异部分超过某个阈值,就认为是一个异常事件。该模型奠定了异常检测的基础,以后发展的许多异常检测方法和系统都是以它为基础而发展起来的。
近几年在异常检测技术的发展过程中,引入了更多人工智能的方法,以提高异常检测的性能。这些人工智能的方法主要包括数据挖掘、人工神经网络、模糊证据理论等。数据挖掘的方法用来确定在大量的数据集合中什么特征是最重要的。该技术用于异常检测中主要是寻求一种正常模式更简洁的定义,而不是像传统的异常检测方法那样简单列举出所有的正常模式。数据挖掘方法的引入使得检测系统能仅通过识别正常模式中的主要特征,就能够概括性地包括训练数据中所未包括的正常模式。人工神经网络异常检测问题可被看作是一个一般的数据分类问题.在前面谈到的统计异常检测中,用户行为数据按照某种统计准则被分为两类:即异常行为和正常行为。由于基于统计的方法在提取、抽象审计实例时存在一定困难,可能造成较大误差,必须依赖于一些概率分布假设,一般需要凭经验和感觉来刻画用户行为的度量,所以引入了人工神经网络的聚类方法。人工神经网络具有自学习自适应能力,用代表正常用户行为的样本点来训练神经网络,通过反复多次学习,神经网络能从数据中提取正常的用户或系统活动的模式,并编码到网络结构中,检测时,将审计数据通过学习好的神经网络,即可判定系统是否正常。由于异常的评判标准具有一定的模糊性,所以模糊证据理论被引入到异常中,如建立一种基于模糊专家系统的入侵检测框架模型,能较好地降低漏警率和虚警率。
本专利提出一种基于统计特征的异常检测方案,基本思想是通过Jensen-Shannon散度衡量历史交通特征与实时交通特征的差异,以实现异常交通状况的检测。该方案具有可解释性好,计算负担不大的优点,既克服了采用单一统计量检测不准确、不及时的弱点,又避免了人工神经网络等方法计算负担重,硬件要求高的缺陷。
9)异常严重性量化表征及异常信息发布。
交通异常状况的严重性应该通过简洁明了的方式向公众发布,以避开可能的拥堵区域,提高城市交通的运行效率。异常状况的严重程度用交通异常指数表征,范围为0-10,其中0表示无异常,10表示高度异常。
异常的发生位置投影到电子地图上,并通过智能移动设备APP等形式公开发布。
10)系统性能评价。
系统性能的评价是指评价交通异常状态检测的准确性,其评价指标包括误报率和漏报率。误报率和漏报率越低表明系统的性能越好。
所述步骤1)中,时空子区的划分具体可以采用以下方法:
11)等距时空划分法。确定时间维度的片段尺度,时间片段跨度为固定值,通常取30min作为一个时间片段;确定空间维度的片段尺度,空间片段跨度为固定值,通常取200m×200m的空间网格作为一个空间片段;
12)基于路网密度的非等距时空划分法:基于路网密度作为判断指标,当路网密度大于等于2km/km2时,取30min的时间片段和200m×200m的空间片段;当路网密度小于2km/km2时,取30min的时间片段和400m×400m的空间片段;
13)基于高峰小时流量的非等距时空划分法:基于高峰小时流量作为判断指标,当高峰小时流量大于等于1000辆/小时时,取30min的时间片段和200m×200m的空间片段;当高峰小时流量小于1000辆/小时时,取30min的时间片段和400m×400m的空间片段。
所述步骤3)具体包含以下步骤:
31)将所需处理的空间区域划分为一定大小的格网,每个格网区域的范围可表示为As={(xs,ys)|xs∈[xr,xr+1),ys∈[yr,yr+1)},每个格网区域包含若干个路段,把这些路段的集合表示为Rs,所述路段的集合Rs中的每条路段表示为ij,并为每个路段赋予编号;
32)判定定位点所在的格网区域,并利用距离和方位角,在路段的集合Rs中搜索某定位点A所在的路段ij,匹配方案包括:
321)单点匹配方案:
搜索距离点A最近的路段,当满足点A的行驶方向角与路段ij的方向角的差值小于阈值时,即满足|θA-θij|<δθ,完成匹配,所述的阈值δθ可取2.5°,5°,10°等;若不满足|θA-θij|<δθ,在搜索空间中删除路段ij,并继续搜索其他路段,直至满足条件。匹配方法如图3所示。
322)点序列匹配方案:
本方案适用于高频浮动车数据。将浮动车GNSS数据采集频率表示为f0=1/t0,将时间上与A相邻的点P(tA-t0),P(tA+t0)定义为A的1-邻近点,P(tA-2t0),P(tA+2t0)定义为A的2-邻近点,以此类推,则P(tA-kt0),P(tA+kt0)定义为A的k-邻近点。在f0<1Hz时,取k=1或2。取距离A及A的k-邻近点距离最小的路段ij,并计算A及A的k-邻近点行驶方向角的均值若满足完成匹配;否则,搜索其他路段,直至满足
33)利用路段的直线方程(若为曲线路段则近似拆分为直线),计算GNSS定位点在路段上的投影坐标,减小因GNSS定位漂移带来的误差。具体方法采用GNSS定位点直线投影法为:
确定路段ij的直线方程(若路段为曲线,则划分为若干直线路段):
y-yi=k(x-xi)
在地图匹配过程后,结合定位点坐标的时间戳数据,将定位点匹配到时空子区。
所述步骤5)具体可以采用以下方法之一:
51)行程时间信息的全样本方案。由给定的OD对之间相同路径的所有行程时间数据,构成总体。实施方法是计算时空子区ξ内每辆车的行程时间:tξ=t1,2+t2,3+...+tn-1,n,其中t1,2...tn-1,n为时空子区ξ内的第1个和第2个GNSS定位点间的行程时间,......,第n-1个与第n个GNSS定位点间的行程时间;将每个时空子区内的数据不做筛选,构成一个集合Tξ,用于后续处理。
52)行程时间信息的时间平滑抽样方案。指定时间片段长度,设置同一时间片段数据条数上限;搜索一个时空子区内时间各时间片段内的行程时间数据,若时间片段内行程时间数据条数超过上限,随机取上限条数的数据加入待处理数据样本。实施方法是计算时空子区ξ内每辆车的行程时间:tξ=t1,2+t2,3+...+tn-1,n,其中t1,2...tn-1,n为时空子区ξ内的第1个和第2个GNSS定位点间的行程时间,......,第n-1个与第n个GNSS定位点间的行程时间;指定时间片段长度tp,同一时间片段数据条数上限pmax;搜索一个时空子区内时间第i各时间片段内的行程时间数据,若时间片段内行程时间数据条数超过上限pmax,随机取pmax条数据加入Tξ并用于后续处理。
所述步骤6)具体可以采用以下方法:
61)简单历史轨迹数据融合法。将无交通异常状况下的历史数据,作为一个整体,进行交通特征模型建立和参数估计。该方法利用有限混合模型,建立交通特征模型,并进行参数估计。具体可采用以下三种方案之一:
611)固定成分的混合高斯模型
本方案采用固定成分数量K的混合高斯模型描述行程时间的概率分布。成分数量根据行程时间在典型情况下的分布模式人工指定。为了保证概率分布的可靠性,成分数量K不能过小。一般可取K=4~6。
612)成分数量可变的混合高斯模型
本方案采用一种基于模型评价的方法来选择合适的成分数量,方法如下:
确定可能的最大成分数量K,并分别对n=1,2,…,K个成分的混合高斯模型进行参数估计;对于K个模型,通过贝叶斯信息准则(BIC)确定最佳模型。最大成分的数量一般按精度需求选取,但必须注意成分数量越多,期望最大化算法收敛越慢。这里选择的最大成分数量为K=5,即需要计算:
共5种混合模型。同时,计算5种模型的BIC,其定义为:
BIC=-2ln L+k·ln n
式中,L为最大似然函数值,k为模型中参数的个数,n为数据总量。
之后,选取BIC最小的的混合模型,记录其参数向量η、μ、σ,其中,η是所述历史交通特征模型中各个子成分所占的比例向量,μ是所述历史交通特征模型中各个子成分的均值向量,σ是所述历史交通特征模型中各个子成分的标准差向量,作为本时空子区的特征记录。混合模型的密度曲线形态在图6中示出。
613)成分数量、分布类型均可变的有限混合模型
本方案采用与612)相同的基于模型评价的方法,但子成分的分布形态和成分的数量均可变,方法如下:
选取M种概率分布模型作为子成分的分布类型,包括但不限于:正态分布、伽马分布、威布尔分布。使用正态分布时,子分布函数采用:
使用伽马分布时,子分布函数采用:
使用威布尔分布时,子分布函数采用:
假定混合模型所有子成分的分布类型相同,确定可能的最大成分数量K。对于M种子成分分布类型、K种成分数量的选择,共形成M·K种组合,分别计算BIC值,并取BIC最小的模型为最佳模型。
62)分情境的历史轨迹数据分类法。依据气温、降水量、能见度和交通管制措施,将无交通异常状况下的历史数据划分成不同的类别,分别建立模型和进行参数估计。实施方法如下:
根据气温、降水量、能见度和交通管制措施的不同,将交通环境分为5~8个类别,由历史数据对应的交通环境的不同,将历史数据归入以上各个类别中。对每个类别,分别进行如同51)所述的处理,从而建立了映射关系R(E→T),E为交通环境,T为交通态势。
63)历史数据聚类法。对于历史数据,通过时空子区两两之间的比较,获得不同时空区域的差异量化表征,并利用量化后的差异进行聚类。将气温、降水量、能见度和交通管制措施作为特征因子,进行多项Logit回归,建立交通环境与类别的映射关系。实施流程参见附图4。实施步骤如下:
631)根据51)所述的方法,建立交通特征模型,并进行参数估计。
632)根据之前的有限混合模型参数估计结果,写出时空子区在不同日期对应的行程时间分布的概率密度函数pi(x),其参数以混合高斯模型为例:
式中,K表示行程时间分布的子成分数量,η表示行程时间分布中某个子成分所占的比例,μ表示行程时间分布中某个子成分的均值,σ表示行程时间分布中某个子成分的标准差。
633)计算各分布两两之间的Jensen-Shannon散度dij:
在采用有限混合模型的情况下,其值无法显式表示,但可采用蒙特卡罗抽样方法近似计算,其计算方法是:
式中,DMC表示采用蒙特卡罗抽样方法近似计算得到的Kullback-Leibler散度,f和g表示任意两个分布函数。
634)将分布两两间的散度表示成距离矩阵:
该矩阵满足dij=dji,dij=0(i=j)。
635)将距离矩阵作为K-Medoids算法的输入,得到聚类结果,并对类别建立索引。
636)以类别索引为响应变量,将交通环境数据(包括气温、降水量、能见度等)作为自变量,进行多项Logit回归,获取交通环境E与交通态势类别T的映射关系R(E→T)。
637)将相同类别的数据进行聚合,并利用聚合后新的数据集重新建立混合模型,并进行参数估计,得到最终的历史交通特征数据集。
所述步骤7)具体可以采用以下方法:
71)简单实时数据处理法。该方法与61)同时实施。将实时交通数据进行模型建立和参数估计,获取当前交通状况的特征函数。该方法的实施步骤,与61)完全相同,只是采用的数据是实时交通数据。
72)分类处理法。该方法与62)或63)同时实施。获取交通状况的特征函数,同时获取当前的气温、降水量、能见度、交通管制措施等信息,并判断当前交通状况的类别。实施流程参见附图5。实施步骤如下:
721)计算时空子区内特定起终点间的行程时间,构成实时行程时间总体Tξ,rt;
723)将当前交通环境数据(包括气温、降水量、能见度等)作为输入参数,利用映射关系R(E→T)获得当前交通态势的所述类别T。
所述步骤8)具体包含以下步骤:
81)当采用步骤72)时,根据当前交通态势所属类别T,定位该类别下历史交通特征数据,否则不进行处理;
82)根据当前交通特征的描述参数ηrt、μrt、σrt和历史交通特征的描述参数η、μ、σ计算两个时间分布间的差异:diff[(ηrt,μrt,σrt),(η,μ,σ)]=JSD(Prt||P)。其中,ηrt是实时交通特征模型中各个子成分所占的比例向量,μrt是实时交通特征模型中各个子成分的均值向量,σrt是实时交通特征模型中各个子成分的标准差向量;η是历史交通特征模型中各个子成分所占的比例向量,μ是历史交通特征模型中各个子成分的均值向量,σ是历史交通特征模型中各个子成分的标准差向量。当历史交通特征与实时交通特征(即历史行程时间分布与实时行程时间分布)相近时,将得到较小的Jensen-Shannon散度值,即两者之间的差异较小;当历史交通特征与实时交通特征差别较大时,将得到较大的Jensen-Shannon散度值,即两者之间的差异较大,亦即存在异常的概率较大,参见附图7。
所述步骤9)具体包含以下步骤:
91)将各个时空子区的行程时间分布差异标准化为0~1的规范化数值aξ:
92)计算各个时空子区的交通异常指数Aξ=aξ×10;
93)将异常指数高于5的区域位置投影到电子地图上,并智能移动设备APP等形式向社会公开发布,以使司机避开潜在拥堵点,提高城市道路交通的通行效率。
所述步骤10)具体包含以下步骤:
101)计算交通异常状态的漏报率:
102)计算交通异常状态的误报率:
以上两式中,nl为单位时间内漏报事件总数,ne为单位时间内误报事件总数,na为单位时间内实际发生的异常时间总数。
本发明相较于同一领域的相似技术,具有以下优点:
(1)充分利用现有的浮动车运营数据(GNSS轨迹数据),通过历史交通特征提取和实时交通态势分析,检测交通状态发生的变化,可以实现城市道路交通异常事件实时性、低成本、智能化检测;
(2)将交通特征参数的概率分布作为交通特征的描述,反映的特征更加全面,避免了利用单一指数表征交通特征的片面性、不稳定性,检测的可靠性更高;
(3)针对交通特征受到交通环境(如天气状况)影响的特点,引入了聚类-多项Logit回归联合算法,建立了交通环境特征与交通态势类别的映射关系;
(4)经实际数据的检验,本发明提出的基于浮动车数据的城市道路交通异常检测技术,能够实现准确度较高的异常事件检测,检测率超过90%,漏报率低于15%,误报率低于20%,取得了良好的检测效果,可以应用于城市交通智能化管理、服务。
附图说明
本发明的具体内容及优点结合以下附图将变得明晰和易于理解,其中:
图1示出了本发明的组成要素和基本原理示意图;
图2示出了本发明在实施过程中的总体流程示意图;
图3示出了本发明快速地图匹配算法实施方式示意图;
图4示出了本发明实施历史交通特征提取方案的流程示意图;
图5示出了本发明实施实时交通特征提取方案的流程示意图;
图6示出了高斯混合模型概率分布的形态示意图;
图7示出了历史交通特征与实时交通特征比较过程中差异的衡量示意图。
具体实施方案
为了更加清晰明确地表述本发明的目的、技术方案和优势,下面对本发明的具体实施方案进行详细描述。
如附图1所示,本发明的整体系统构架包括:浮动车搭载的车载GNSS轨迹记录仪、数据中心、GNSS卫星以及通信系统。此处的GNSS包括GPS、GLONASS、GALILEO、北斗、IRNSS、QZSS等任何类似的导航卫星定位系统。出租车、公交车等浮动车搭载的GNSS轨迹记录仪,以一定的采样频率f(一般要求f>0.1Hz)记录车辆在行驶中各时点的位置信息,并通过GPRS移动通信网络(亦可采用WCDMA、TD-LTE等无线网络通信技术,但成本将相应提高)将位置信息实时发送至数据中心。数据中心通过数据预处理、数据融合,并通过特定算法建立历史道路交通特征数据库;对于最近接收的实时数据,建立实时交通特征数据库;通过历史数据库和实时数据库的映射关系,判别当前交通特征是否异常,并通过处理终端进行可视化展示并生成交通异常事件报告。
方案的总体流程参见图2,包括采集和存储GNSS轨迹数据,建立时空子区,历史交通特征提取,实时交通特征提取,异常识别等步骤。采集和存储GNSS轨迹数据,是整个方案的数据基础,由于数据量级巨大,应采用分布式存储方案,对于分布式存储目前已有成熟的技术,不是本发明的内容。建立时空子区,其基本假设是在某一特定区域、特定时段内,有着相同的交通特征,这一假设,经过长期观测,是普遍适用的。历史交通特征提取,其原理是利用GNSS轨迹数据,计算得到行程时间,利用同一时空子区大量的行程时间数据,建立行程时间的概率分布模型,并进行参数估计,用少量参数表征交通特征。实时交通特征提取,其原理是将当前时间段内的行程时间数据进行处理分析,同样建立当前的行程时间概率分布模型。异常识别是采用差异衡量指标,判断实时特征相较于历史特征的变化程度,根据其是否达到阈值,确定是否出现交通异常事件。
根据发明内容所述实施方法的组合,给出实施例如下。
实施例一
步骤11、采用等距时空划分法,确定时间维度的片段尺度,时间片段跨度为固定值,通常取30min作为一个时间片段;确定空间维度的片段尺度,空间片段跨度为固定值,通常取200m×200m的空间网格作为一个空间片段。
步骤12、进行数据预处理,将GNSS定位数据进行数据清洗、数据集成、数据转换、数据归约,提高数据的结构化程度。
步骤13、将所需处理的空间区域划分为一定大小的格网,每个格网区域的范围可表示为As={(xs,ys)|xs∈[xr,xr+1),ys∈[yr,yr+1)};判定定位点所在的格网区域,并利用距离和方位角,搜索定位点所在的路段;搜索距离点A最近的路段,取阈值δθ=2.5°,当满足点A的行驶方向角与路段ij的方向角的差值小于阈值δθ时,即满足|θA-θij|<δθ,完成匹配;若不满足|θA-θij|<δθ,在搜索空间中删除路段ij,并继续搜索其他路段,直至满足条件;利用路段的直线方程(若为曲线路段则近似拆分为直线),计算GNSS定位点在路段上的投影坐标,减小因GNSS定位漂移带来的误差,具体方法为:
确定路段ij的直线方程(若路段为曲线,则划分为若干直线路段):
y-yi=k(x-xi)
在地图匹配过程后,结合定位点坐标的时间戳数据,将定位点匹配到时空子区。
步骤14、由给定的OD对之间相同路径的所有行程时间数据,构成总体。实施方法是计算时空子区ξ内每辆车的行程时间:tξ=t1,2+t2,3+...+tn-1,n,其中t1,2...tn-1,n为时空子区ξ内的第1个和第2个GNSS定位点间的行程时间,......,第n-1个与第n个GNSS定位点间的行程时间;将每个时空子区内的数据不做筛选,构成一个集合Tξ,用于后续处理。
步骤15、将无交通异常状况下的历史数据,作为一个整体,进行交通特征模型建立和参数估计。该方法利用有限混合模型,建立交通特征模型,并进行参数估计。取最大成分数量K=5,并分别对n=1,2,…,K个成分的混合高斯模型进行参数估计;对于K个模型,通过贝叶斯信息准则(BIC)确定最佳模型。计算:
共5种混合模型。同时,计算5种模型的BIC:
BIC=-2ln L+k·ln n
式中,L为最大似然函数值,k为模型中参数的个数,n为数据总量。
之后,选取BIC最小的的混合模型,记录其参数向量η、μ、σ,作为本时空子区的特征记录。
步骤16、将实时交通数据进行模型建立和参数估计,获取当前交通状况的特征函数,方法同步骤15,记录参数向量ηrt、μrt、σrt。
步骤17、根据当前交通特征的描述参数ηrt、μrt、σrt和历史交通特征的描述参数η、μ、σ计算两个行程时间分布间的差异:diff[(ηrt,μrt,σrt),(η,μ,σ)]=JSD(Prt||P)。
步骤18、将各个时空子区的行程时间分布差异标准化为0~1的规范化数值aξi:
计算各个时空子区的交通异常指数Aξ=aξ×10。
实施例二
步骤21、采用等距时空划分法,确定时间维度的片段尺度,时间片段跨度为固定值,通常取30min作为一个时间片段;确定空间维度的片段尺度,空间片段跨度为固定值,通常取200m×200m的空间网格作为一个空间片段。
步骤22、进行数据预处理,将GNSS定位数据进行数据清洗、数据集成、数据转换、数据归约,提高数据的结构化程度。
步骤23、将所需处理的空间区域划分为一定大小的格网,每个格网区域的范围可表示为As={(xs,ys)|xs∈[xr,xr+1),ys∈[yr,yr+1)};判定定位点所在的格网区域,并利用距离和方位角,搜索定位点所在的路段;搜索距离点A最近的路段,取阈值δθ=2.5°,当满足点A的行驶方向角与路段ij的方向角的差值小于阈值δθ时,即满足|θA-θij|<δθ,完成匹配;若不满足|θA-θij|<δθ,在搜索空间中删除路段ij,并继续搜索其他路段,直至满足条件;利用路段的直线方程(若为曲线路段则近似拆分为直线),计算GNSS定位点在路段上的投影坐标,减小因GNSS定位漂移带来的误差,具体方法为:
确定路段ij的直线方程(若路段为曲线,则划分为若干直线路段):
y-yi=k(x-xi)
在地图匹配过程后,结合定位点坐标的时间戳数据,将定位点匹配到时空子区。
步骤24、计算时空子区ξ内每辆车的行程时间:tξ=t1,2+t2,3+...+tn-1,n,其中t1,2...tn-1,n为时空子区ξ内的第1个和第2个GNSS定位点间的行程时间,......,第n-1个与第n个GNSS定位点间的行程时间;指定时间片段长度tp,同一时间片段数据条数上限pmax;搜索一个时空子区内时间第i各时间片段内的行程时间数据,若时间片段内行程时间数据条数超过上限pmax,随机取pmax条数据加入Tξ。
步骤25、将无交通异常状况下的历史数据,作为一个整体,进行交通特征模型建立和参数估计。该方法利用有限混合模型,建立交通特征模型,并进行参数估计。取最大成分数量K=5,并分别对n=1,2,…,K个成分的混合高斯模型进行参数估计;对于K个模型,通过贝叶斯信息准则(BIC)确定最佳模型。计算:
共5种混合模型。同时,计算5种模型的BIC:
BIC=-2ln L+k·ln n
式中,L为最大似然函数值,k为模型中参数的个数,n为数据总量。
之后,选取BIC最小的的混合模型,记录其参数向量η、μ、σ,作为本时空子区的特征记录。
根据参数估计结果,写出时空子区在不同日期对应的行程时间分布的概率密度函数pi(x):
计算各分布两两之间的Jensen-Shannon散度dij:
在采用有限混合模型的情况下,采用蒙特卡罗抽样方法近似计算,其计算方法是:
将分布两两间的散度表示成距离矩阵:
该矩阵满足dij=dji,dij=0(i=j)。
将距离矩阵作为K-Medoids算法的输入,得到聚类结果,并对类别建立索引。
以类别索引为响应变量,将交通环境数据(包括气温、降水量、能见度等)作为自变量,进行多项Logit回归,获取交通环境E与交通态势类别T的映射关系R(E→T)。
将相同类别的数据进行聚合,并利用聚合后新的数据集重新建立混合模型,并进行参数估计,得到最终的历史交通特征数据集。
步骤26、获取交通状况的特征函数,同时获取当前的气温、降水量、能见度、交通管制措施等信息,并判断当前交通状况的类别。
计算时空子区内特定起终点间的行程时间,构成实时行程时间总体Tξ,rt;
将当前交通环境数据(包括气温、降水量、能见度等)作为输入参数,利用映射关系R(E→T)获得当前交通态势的所述类别T。
步骤27、根据当前交通态势所属类别T,定位该类别下历史交通特征数据;根据当前交通特征的描述参数ηrt、μrt、σrt和历史交通特征的描述参数η、μ、σ计算两个行程时间分布间的差异:diff|(ηrt,μrt,σrt),(η,μ,σ)|=JSD(Prt||P)。
步骤28、将各个时空子区的行程时间分布差异标准化为0~1的规范化数值aξ:
计算各个时空子区的交通异常指数Aξ=aξ×10。
实施例三
步骤31、采用非等距时空划分法,对于路网密度大于2km/km2或高峰小时流量大于1000辆/小时的城市中心区,取30min的时间片段和200m×200m的空间片段,对于路网密度小于2km/km2或高峰小时流量小于1000辆/小时的城市郊区,取30min的时间片段和400m×400m的空间片段。
步骤32、进行数据预处理,将GNSS定位数据进行数据清洗、数据集成、数据转换、数据归约,提高数据的结构化程度。
步骤33、将所需处理的空间区域划分为一定大小的格网,每个格网区域的范围可表示为As={(xs,ys)|xs∈[xr,xr+1),ys∈[yr,yr+1)};
将浮动车GNSS数据采集频率表示为f0=1/t0,将时间上与A相邻的点P(tA-t0),P(tA+t0)定义为A的1-邻近点,P(tA-2t0),P(tA+2t0)定义为A的2-邻近点,以此类推,则P(tA-kt0),P(tA+kt0)定义为A的k-邻近点。在f0<1Hz时,取k=1或2。取距离A及A的k-邻近点距离最小的路段ij,并计算A及A的k-邻近点行驶方向角的均值取阈值δθ=5°,若满足完成匹配;否则,搜索其他路段,直至满足条件。
利用路段的直线方程(若为曲线路段则近似拆分为直线),计算GNSS定位点在路段上的投影坐标,减小因GNSS定位漂移带来的误差。具体方法为:
确定路段ij的直线方程(若路段为曲线,则划分为若干直线路段):y-yi=k(x-xi)
在地图匹配过程后,结合定位点坐标的时间戳数据,将定位点匹配到时空子区。
步骤34、计算时空子区ξ内每辆车的行程时间:tξ=t1,2+t2,3+...+tn-1,n,其中t1,2...tn-1,n为时空子区ξ内的第1个和第2个GNSS定位点间的行程时间,......,第n-1个与第n个GNSS定位点间的行程时间;指定时间片段长度tp,同一时间片段数据条数上限pmax;搜索一个时空子区内时间第i各时间片段内的行程时间数据,若时间片段内行程时间数据条数超过上限pmax,随机取pmax条数据加入Tξ。
步骤35、将无交通异常状况下的历史数据,作为一个整体,进行交通特征模型建立和参数估计。该方法利用有限混合模型,建立交通特征模型,并进行参数估计。取最大成分数量K=5,并分别对n=1,2,…,K个成分的混合高斯模型进行参数估计;对于K个模型,通过贝叶斯信息准则(BIC)确定最佳模型。计算:
共5种混合模型。同时,计算5种模型的BIC:
BIC=-2ln L+k·ln n
式中,L为最大似然函数值,k为模型中参数的个数,n为数据总量。
之后,选取BIC最小的的混合模型,记录其参数向量η、μ、σ,作为本时空子区的特征记录。
根据参数估计结果,写出时空子区在不同日期对应的行程时间分布的概率密度函数pi(x):
计算各分布两两之间的Jensen-Shannon散度dij:
在采用有限混合模型的情况下,采用蒙特卡罗抽样方法近似计算,其计算方法是:
将分布两两间的散度表示成距离矩阵:
该矩阵满足dij=dji,dij=0(i=j)。
将距离矩阵作为K-Medoids算法的输入,得到聚类结果,并对类别建立索引。
以类别索引为响应变量,将交通环境数据(包括气温、降水量、能见度等)作为自变量,进行多项Logit回归,获取交通环境E与交通态势类别T的映射关系R(E→T)。
将相同类别的数据进行聚合,并利用聚合后新的数据集重新建立混合模型,并进行参数估计,得到最终的历史交通特征数据集。
步骤36、获取交通状况的特征函数,同时获取当前的气温、降水量、能见度、交通管制措施等信息,并判断当前交通状况的类别。
计算时空子区内的行程时间,构成实时行程时间总体Tξ,rt;
将当前交通环境数据(包括气温、降水量、能见度等)作为输入参数,利用映射关系R(E→T)获得当前交通态势的所述类别T。
步骤37、根据当前交通态势所属类别T,定位该类别下历史交通特征数据;根据当前交通特征的描述参数ηrt、μrt、σrt和历史交通特征的描述参数η、μ、σ计算两个行程时间分布间的差异:diff|(ηrt,μrt,σrt),(η,μ,σ)|=JSD(Prt||P)。
步骤38、将各个时空子区的行程时间分布差异标准化为0~1的规范化数值aξ:
计算各个时空子区的交通异常指数Aξ=aξ×10。
Claims (6)
1.一种基于行程时间分布的交通异常检测方法,包括如下步骤:
1)建立时空子区:将一天划分为若干时间片段,每个时间片段称为一个时间子区;将城市道路交通异常检测的实施区域划分为若干空间片段,每个空间片段称为一个空间子区;任意一个时间子区和任意一个空间子区的交集称为时空子区;
2)历史轨迹数据的预处理:将浮动车GNSS定位历史数据处理为历史轨迹的抽样行程时间数据;
实时轨迹数据的预处理:将浮动车GNSS定位实时数据处理为实时轨迹的抽样行程时间数据;
3)历史轨迹数据分析和特征提取:利用所述历史轨迹的抽样行程时间数据,建立历史行程时间概率分布,得到历史交通特征模型Ph,实施方法是:将无交通异常状况下的历史数据,作为一个整体,该方法利用有限混合模型,建立交通特征模型,并进行参数估计;
实时轨迹数据分析和特征提取:利用所述实时轨迹的抽样行程时间数据;建立实时行程时间概率分布,得到实时交通特征模型Prt;
4)异常检测:通过Jensen–Shannon散度衡量所述历史交通特征模型与所述实时交通特征模型的差异,根据实时交通特征模型的描述参数ηrt、μrt、σrt和历史交通特征模型描述参数η、μ、σ计算两个时间分布间Jensen-Shannon散度:diff[(ηrt,μrt,σrt),(η,μ,σ)]=JSD(Prt||P);其中,ηrt是实时交通特征模型中各个子成分所占的比例向量,μrt是实时交通特征模型中各个子成分的均值向量,σrt是实时交通特征模型中各个子成分的标准差向量;η是历史交通特征模型中各个子成分所占的比例向量,μ是历史交通特征模型中各个子成分的均值向量,σ是历史交通特征模型中各个子成分的标准差向量;通过利用历史交通特征模型与实时交通特征模型计算得到的Jensen-Shannon散度计算得到历史与实时交通特征差异值;
5)异常严重性量化表征:利用所述历史与实时交通特征差异值,计算交通状况异常指数;
6)系统性能评价:评价交通异常状态检测的准确性,衡量系统运行的稳定程度;步骤3)采用下述方法之一:
3a)固定成分混合高斯模型法:采用固定成分数量K的混合高斯模型描述行程时间的概率分布,成分数量根据行程时间在典型情况下的分布模式人工指定,成分数量K取4~6;将实时交通数据进行模型建立和参数估计,获取当前交通状况的特征函数;
3b)可变混合高斯模型法:采用可变的成分数量,或者可变的成分数量和可变的分布类型;采用分类处理法,获取交通状况的特征函数,同时获取当前的气温、降水量、能见度、交通管制措施信息,并判断当前交通状况的类别;
步骤3b)所述的可变混合高斯模型法包括下述方法之一:
3b1)成分数量可变的混合高斯模型法:采用基于模型评价的方法来选择合适的成分数量,方法如下:确定可能的最大成分数量K,并分别对n=1,2,…,K个成分的混合高斯模型进行参数估计;对于K个模型,通过贝叶斯信息准则确定最佳模型;
3b2)成分数量、分布类型均可变的有限混合模型法:子成分的分布形态和成分的数量均可变;
步骤3b)所述的分类处理法包括:
3b3)计算时空子区内的OD对间的行程时间tξ,rt,构成实时行程时间总体Tξ,rt;
3b4)建立实时交通特征模型式中,K表示所述实时交通特征的子成分数量,η表示所述实时交通特征的某个子成分所占的比例,μ表示所述实时交通特征模型的某个子成分的均值,σ表示所述实时交通特征模型的某个子成分的标准差,并进行参数估计,得到当前实时交通特征的描述参数ηrt、μrt、σrt,其中,ηrt是所述实时交通特征模型中各个子成分所占的比例向量,μrt是所述实时交通特征模型中各个子成分的均值向量,σrt是所述实时交通特征模型中各个子成分的标准差向量;
3b5)将当前交通环境数据,包括气温、降水量、能见度作为输入参数,利用映射关系R(E→T)获得当前交通态势的所述类别T,其中E表示所述交通环境数据。
2.如权利要求1所述的交通异常检测方法,其特征在于,步骤1)采用下述方法之一:
1a)等距时空划分法:确定时间维度的片段尺度,时间片段跨度为固定值,取30min作为一个时间片段;确定空间维度的片段尺度,空间片段跨度为固定值,取200m×200m的空间网格作为一个空间片段;
1b)基于路网密度的非等距时空划分法:基于路网密度作为判断指标,当路网密度大于等于2km/km2时,取30min的时间片段和200m×200m的空间片段;当路网密度小于2km/km2时,取30min的时间片段和400m×400m的空间片段;
1c)基于高峰小时流量的非等距时空划分法:基于高峰小时流量作为判断指标,当高峰小时流量大于等于1000辆/小时时,取30min的时间片段和200m×200m的空间片段;当高峰小时流量小于1000辆/小时时,取30min的时间片段和400m×400m的空间片段。
3.如权利要求1所述的交通异常检测方法,其特征在于,步骤2)所述的历史轨迹数据的预处理包括:
2a)数据结构化:将浮动车GNSS定位历史数据进行数据清洗、数据集成、数据转换、数据归约,得到结构化GNSS定位历史数据;
2b)快速地图匹配:结合城市路网数据,通过地图匹配算法,将结构化的GNSS定位历史数据投影到城市路网,建立所述结构化GNSS定位历史数据中的定位点与路段的匹配关系,得到所述结构化GNSS定位历史数据中的定位点与所述路段的匹配关系表,并修正定位漂移带来的误差;
2c)历史轨迹的行程时间计算和抽样:在每个时空子区内选取若干个OD对,根据所述结构化GNSS定位历史数据计算这些OD对间的行程时间,得到历史轨迹的行程时间数据,并对所述历史轨迹的行程时间数据进行数据抽样,得到历史轨迹的抽样行程时间数据。
4.如权利要求3所述的交通异常检测方法,其特征在于,步骤2b)所述的快速地图匹配包括:
2b1)将所需处理的空间区域划分为一定大小的格网,每个格网区域的范围可表示为As={(xs,ys)|xs∈[xr,xr+1),ys∈[yr,yr+1)},每个格网区域包含若干个路段,把这些路段的集合表示为Rs,所述路段的集合Rs中的每条路段表示为ij,并为每个路段赋予编号;
2b2)判定定位点所在的格网区域,并利用距离和方位角,在路段的集合Rs中搜索某定位点A所在的路段ij;
2b3)利用GNSS定位点直线投影法,计算GNSS定位点在路段上的投影坐标。
5.如权利要求4所述的交通异常检测方法,其特征在于,步骤2b2)采用下述方法之一:
2b21)单点匹配方法:搜索距离某定位点A最近的路段,实施方法是:对于路段的集合Rs中的某一路段ij,当满足点A的行驶方向角θA与路段ij的方向角θij的差值小于阈值δθ时,即满足|θA-θij|<δθ时,完成匹配;若不满足|θA-θij|<δθ,继续搜索路段集合Rs中的其他路段,直至满足
2b22)点序列匹配方法:本方案适用于高频浮动车数据;将每两条相邻时间的浮动车GNSS数据时间间隔表示为t0,将浮动车GNSS数据采集频率表示为f0=1/t0,将某定位点A的时间记录表示为tA,将时间上与所述定位点A相邻的点P(tA-t0),P(tA+t0)定义为A的1-邻近点,P(tA-2t0),P(tA+2t0)定义为某定位点A的2-邻近点,以此类推,则P(tA-kt0),P(tA+kt0)定义为某定位点A的k-邻近点;在f0<1Hz时,取k=1或2;取距离某定位点A及该定位点的k-邻近点距离最小的一个路段ij,并计算该定位点及其k-邻近点行驶方向角的均值若满足完成匹配;否则,搜索其他路段,直至满足
6.如权利要求3所述的交通异常检测方法,其特征在于,步骤2c)所述的历史轨迹的行程时间计算和抽样采用下述方法之一:
2c1)行程时间信息的全样本抽样方法:由所述OD对之间相同路径的所有行程时间数据,构成总体;实施方法是计算时空子区ξ内每辆车的行程时间:tξ=t1,2+t2,3+...+tn-1,n,其中t1,2…tn-1,n为时空子区ξ内的第1个和第2个GNSS定位点间的行程时间,……,第n-1个与第n个GNSS定位点间的行程时间;将每个时空子区内的行程时间数据不做筛选,构成一个集合Tξ,用于后续处理;
2c2)行程时间信息的时间平滑抽样方法:指定时间片段长度,设置同一时间片段数据条数上限;搜索某时空子区ξ各时间片段内的行程时间数据,若时间片段内行程时间数据条数超过上限,随机取上限条数的数据用于后续处理;实施方法是计算时空子区ξ内每辆车的行程时间:tξ=t1,2+t2,3+...+tn-1,n,其中t1,2…tn-1,n为时空子区ξ内的第1个和第2个GNSS定位点间的行程时间,……,第n-1个与第n个GNSS定位点间的行程时间;指定时间片段长度tp,同一时间片段数据条数上限pmax;搜索一个时空子区内时间第i各时间片段内的行程时间数据,若时间片段内行程时间数据条数超过上限pmax,随机取pmax条数据构成集合Tξ并用于后续处理。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/IB2016/058105 WO2018122585A1 (zh) | 2016-12-30 | 2016-12-30 | 一种基于浮动车数据的城市道路交通异常检测方法 |
IBPCT/IB2016/058105 | 2016-12-30 | ||
PCT/IB2017/058535 WO2018122805A1 (zh) | 2016-12-30 | 2017-12-30 | 一种基于行程时间分布的交通异常检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109791729A CN109791729A (zh) | 2019-05-21 |
CN109791729B true CN109791729B (zh) | 2021-10-15 |
Family
ID=59713568
Family Applications (7)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201680088595.8A Active CN109923595B (zh) | 2016-12-30 | 2016-12-30 | 一种基于浮动车数据的城市道路交通异常检测方法 |
CN201780050755.4A Active CN109643485B (zh) | 2016-12-30 | 2017-12-30 | 一种城市道路交通异常检测方法 |
CN201780050906.6A Active CN109716414B (zh) | 2016-12-30 | 2017-12-30 | 一种多模态道路交通异常检测方法 |
CN201780050719.8A Pending CN110168520A (zh) | 2016-12-30 | 2017-12-30 | 一种智能化道路交通异常检测方法 |
CN201780050907.0A Active CN109791729B (zh) | 2016-12-30 | 2017-12-30 | 一种基于行程时间分布的交通异常检测方法 |
CN201780050765.8A Active CN109844832B (zh) | 2016-12-30 | 2017-12-30 | 一种基于行程时间分布的多模态交通异常检测方法 |
CN201780050754.XA Pending CN109997179A (zh) | 2016-12-30 | 2017-12-30 | 一种非等距时空划分的道路交通异常检测方法 |
Family Applications Before (4)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201680088595.8A Active CN109923595B (zh) | 2016-12-30 | 2016-12-30 | 一种基于浮动车数据的城市道路交通异常检测方法 |
CN201780050755.4A Active CN109643485B (zh) | 2016-12-30 | 2017-12-30 | 一种城市道路交通异常检测方法 |
CN201780050906.6A Active CN109716414B (zh) | 2016-12-30 | 2017-12-30 | 一种多模态道路交通异常检测方法 |
CN201780050719.8A Pending CN110168520A (zh) | 2016-12-30 | 2017-12-30 | 一种智能化道路交通异常检测方法 |
Family Applications After (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201780050765.8A Active CN109844832B (zh) | 2016-12-30 | 2017-12-30 | 一种基于行程时间分布的多模态交通异常检测方法 |
CN201780050754.XA Pending CN109997179A (zh) | 2016-12-30 | 2017-12-30 | 一种非等距时空划分的道路交通异常检测方法 |
Country Status (3)
Country | Link |
---|---|
CN (7) | CN109923595B (zh) |
GB (11) | GB201711408D0 (zh) |
WO (7) | WO2018122585A1 (zh) |
Families Citing this family (108)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108388540B (zh) * | 2018-03-01 | 2022-09-23 | 兰州交通大学 | 基于模糊综合评价的道路网选取方法 |
CN110766936A (zh) * | 2018-07-25 | 2020-02-07 | 高德软件有限公司 | 基于多源数据融合的交通运行状态感知方法和系统 |
EP3867720A1 (en) * | 2018-10-15 | 2021-08-25 | Starship Technologies OÜ | Method and system for operating a mobile robot |
US11100793B2 (en) * | 2019-01-15 | 2021-08-24 | Waycare Technologies Ltd. | System and method for detection and quantification of irregular traffic congestion |
CN109871876B (zh) * | 2019-01-22 | 2023-08-08 | 东南大学 | 一种基于浮动车数据的高速公路路况识别与预测方法 |
CN111613047B (zh) * | 2019-02-26 | 2022-07-08 | 阿里巴巴集团控股有限公司 | 一种信息处理方法及装置 |
CN109948701B (zh) * | 2019-03-19 | 2022-08-16 | 太原科技大学 | 一种基于轨迹间时空关联性的数据聚类方法 |
CN109918368A (zh) * | 2019-03-27 | 2019-06-21 | 成都市公安科学技术研究所 | 一种通过轨迹关联度识别车辆驾驶者的系统和方法 |
CN110197293A (zh) * | 2019-04-15 | 2019-09-03 | 同济大学 | 基于浮动车数据的分时租赁汽车出行需求预测方法和系统 |
CN110322054B (zh) * | 2019-06-14 | 2023-04-28 | 中交第一公路勘察设计研究院有限公司 | 一种公路路段交通监测器的优化布设方法 |
CN110298500B (zh) * | 2019-06-19 | 2022-11-08 | 大连理工大学 | 一种基于出租车数据和城市路网的城市交通轨迹数据集生成方法 |
CN112135239B (zh) * | 2019-06-25 | 2023-03-17 | 杭州萤石软件有限公司 | 位置监测方法及装置 |
US11587448B2 (en) | 2019-07-26 | 2023-02-21 | General Electric Company | Systems and methods for manifolds learning of airline network data |
CN110634288B (zh) * | 2019-08-30 | 2022-06-21 | 上海电科智能系统股份有限公司 | 基于三元高斯混合模型的多维城市交通异常事件识别方法 |
CN110533249B (zh) * | 2019-09-02 | 2021-09-14 | 合肥工业大学 | 一种基于集成长短期记忆网络的冶金企业能耗预测方法 |
CN110716925B (zh) * | 2019-09-05 | 2023-08-04 | 中电科大数据研究院有限公司 | 一种基于轨迹分析的跨境行为识别方法 |
CN110674858B (zh) * | 2019-09-16 | 2023-06-09 | 长沙理工大学 | 一种基于时空关联与大数据挖掘的交通舆情检测方法 |
CN110795467B (zh) * | 2019-09-26 | 2024-02-27 | 腾讯大地通途(北京)科技有限公司 | 交通规则数据处理方法、装置、存储介质和计算机设备 |
CN110895598B (zh) * | 2019-10-23 | 2021-09-14 | 山东九州信泰信息科技股份有限公司 | 基于多源预测的实时异常检测并行化方法 |
JP7226248B2 (ja) * | 2019-10-31 | 2023-02-21 | トヨタ自動車株式会社 | 通信装置および異常判定装置 |
US11587433B2 (en) | 2019-10-31 | 2023-02-21 | Here Global B.V. | Method, apparatus, and system for probe anomaly detection |
CN110956803A (zh) * | 2019-11-14 | 2020-04-03 | 深圳尚桥信息技术有限公司 | 一种基于多模态的车辆检测方法及系统 |
CN110941278B (zh) * | 2019-12-20 | 2023-05-23 | 交控科技股份有限公司 | 一种站内动态安全分析方法 |
CN111192454B (zh) * | 2020-01-07 | 2021-06-01 | 中山大学 | 基于行程时间演化的交通异常识别方法、系统及存储介质 |
CN111192456A (zh) * | 2020-01-14 | 2020-05-22 | 泉州市益典信息科技有限公司 | 一种道路交通运行态势多时间尺度预测方法 |
CN111311910B (zh) * | 2020-02-19 | 2020-12-22 | 中南大学 | 多层次道路级浮动车异常轨迹探测方法 |
CN111310833B (zh) * | 2020-02-19 | 2022-11-15 | 清华大学 | 一种基于贝叶斯神经网络的出行方式识别方法 |
CN112164223B (zh) * | 2020-02-27 | 2022-04-29 | 浙江恒隆智慧科技集团有限公司 | 基于云平台的智能交通信息处理方法及装置 |
CN111540194B (zh) * | 2020-03-23 | 2021-08-10 | 深圳奇迹智慧网络有限公司 | 车辆监测数据处理方法、装置、计算机设备和存储介质 |
CN111353828B (zh) * | 2020-03-30 | 2023-09-12 | 中国工商银行股份有限公司 | 网点到店人数预测方法及装置 |
CN111475774B (zh) * | 2020-03-31 | 2022-03-18 | 清华大学 | 一种光伏电站设备异常状态检测方法及装置 |
CN111583642B (zh) * | 2020-05-06 | 2021-03-23 | 北京航空航天大学 | 交通轨迹流式大数据实时处理方法 |
CN111580500B (zh) * | 2020-05-11 | 2022-04-12 | 吉林大学 | 一种针对自动驾驶汽车安全性的评价方法 |
CN111722252B (zh) * | 2020-05-14 | 2023-04-11 | 江苏北斗卫星应用产业研究院有限公司 | 一种高精准的基于轨迹数据的作业面识别方法 |
CN111739289B (zh) * | 2020-06-02 | 2024-02-20 | 腾讯科技(深圳)有限公司 | 车辆预警信息的处理方法及装置 |
CN111667015B (zh) * | 2020-06-11 | 2023-06-27 | 深圳市兴海物联科技有限公司 | 物联网设备状态检测方法、装置及检测设备 |
CN111709378B (zh) * | 2020-06-18 | 2022-07-12 | 湖南大学 | 一种基于js散度和模糊证据理论的路况状态评估新方法 |
CN111723997A (zh) * | 2020-06-23 | 2020-09-29 | 上海电科智能系统股份有限公司 | 一种基于gan的城市重大交通事故数据样本自动生成方法 |
CN111932873B (zh) * | 2020-07-21 | 2022-10-04 | 重庆交通大学 | 一种山地城市热点区域实时交通预警管控方法及系统 |
CN112085949B (zh) * | 2020-08-13 | 2023-06-27 | 浙江工业大学 | 一种基于交通运行状况异常的路网脆弱性识别、分析与应对方法 |
CN112000653B (zh) * | 2020-08-21 | 2024-03-29 | 睿驰达新能源汽车科技(北京)有限公司 | 基于空间和时间的区域网格化驾驶行为数据预处理方法 |
CN112053561B (zh) * | 2020-09-11 | 2021-11-23 | 深兰人工智能芯片研究院(江苏)有限公司 | 非监控路段判断与定位交通事故的方法、系统及装置 |
CN111986487B (zh) * | 2020-09-11 | 2022-02-25 | 腾讯科技(深圳)有限公司 | 一种路况信息的管理方法以及相关装置 |
CN112100243B (zh) * | 2020-09-15 | 2024-02-20 | 山东理工大学 | 一种基于海量时空数据分析的异常聚集检测方法 |
CN112215261A (zh) * | 2020-09-18 | 2021-01-12 | 武汉理工大学 | 基于元学习的车辆od点聚类方法、系统、装置及存储介质 |
CN112308116B (zh) * | 2020-09-28 | 2023-04-07 | 济南大学 | 一种助老陪护机器人的自调优多通道融合方法和系统 |
CN116368544A (zh) * | 2020-10-16 | 2023-06-30 | 格步计程车控股私人有限公司 | 用于检测超速的方法、电子装置及系统 |
CN114882696B (zh) * | 2020-10-28 | 2023-11-03 | 华为技术有限公司 | 道路容量的确定方法、装置及存储介质 |
CN112419722B (zh) * | 2020-11-18 | 2022-08-30 | 百度(中国)有限公司 | 交通异常事件检测方法、交通管控方法、设备和介质 |
CN114945959B (zh) * | 2020-11-23 | 2023-06-20 | 深圳元戎启行科技有限公司 | 行驶轨迹确定方法、装置、计算机设备和存储介质 |
CN112652170B (zh) * | 2020-12-24 | 2022-04-08 | 航天科工智能运筹与信息安全研究院(武汉)有限公司 | 交通拥堵热点的定位方法和装置 |
CN112927497B (zh) * | 2021-01-14 | 2023-01-17 | 阿里巴巴集团控股有限公司 | 一种浮动车识别方法、相关方法和装置 |
CN112863185A (zh) * | 2021-01-15 | 2021-05-28 | 林安齐 | 一种道路交通设备智能管理系统及方法 |
CN112907949B (zh) * | 2021-01-20 | 2022-11-22 | 北京百度网讯科技有限公司 | 交通异常的检测方法、模型的训练方法及装置 |
CN112784963B (zh) * | 2021-01-22 | 2022-07-01 | 重庆邮电大学 | 基于模拟退火优化bp神经网络的室内外无缝定位方法 |
CN113221677B (zh) * | 2021-04-26 | 2024-04-16 | 阿波罗智联(北京)科技有限公司 | 一种轨迹异常检测方法、装置、路侧设备和云控平台 |
GB202106070D0 (en) * | 2021-04-28 | 2021-06-09 | Tomtom Navigation Bv | Methods and systems for determining estimated travel times through a navigable network |
CN113240265B (zh) * | 2021-05-11 | 2023-10-27 | 西北工业大学 | 一种基于多模式交通数据的城市空间划分方法 |
CN113470347B (zh) * | 2021-05-20 | 2022-07-26 | 上海天壤智能科技有限公司 | 结合卡口过车记录和浮动车gps数据的拥堵识别方法及系统 |
CN113312722B (zh) * | 2021-05-28 | 2023-05-05 | 广西大学 | 一种城轨列车关键部件的可靠度预测优化方法 |
CN113380028B (zh) * | 2021-06-01 | 2022-09-06 | 公安部交通管理科学研究所 | 一种智慧出行交通数据融合方法及装置 |
CN113313317B (zh) * | 2021-06-11 | 2024-04-12 | 哈尔滨工业大学 | 一种基于inla算法的共享单车使用需求预测方法及预测系统 |
CN113378301B (zh) * | 2021-06-22 | 2022-05-24 | 北京航空航天大学 | 一种基于重要度抽样的无人车超车场景关键测试案例生成方法 |
CN113535510B (zh) * | 2021-06-24 | 2024-01-26 | 北京理工大学 | 一种大规模数据中心数据采集的自适应抽样模型优化方法 |
CN113436433B (zh) * | 2021-06-24 | 2022-06-21 | 福建师范大学 | 一种高效的城市交通离群值检测方法 |
CN113343905B (zh) * | 2021-06-28 | 2022-06-14 | 山东理工大学 | 道路异常智能识别模型训练、道路异常识别的方法及系统 |
CN113554869B (zh) * | 2021-07-01 | 2022-04-05 | 华东师范大学 | 一种基于多特征融合的道路封闭检测方法 |
CN113569759B (zh) * | 2021-07-29 | 2022-06-10 | 沭阳新辰公路仪器有限公司 | 一种基于人工智能的道路掉落物识别方法及系统 |
US11828860B2 (en) | 2021-08-27 | 2023-11-28 | International Business Machines Corporation | Low-sampling rate GPS trajectory learning |
CN113779169B (zh) * | 2021-08-31 | 2023-09-05 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | 时空数据流模型自增强方法 |
CN113920728B (zh) * | 2021-10-11 | 2022-08-12 | 南京微达电子科技有限公司 | 高速公路抛洒障碍物检测与预警方法及系统 |
CN113920739B (zh) * | 2021-10-29 | 2023-05-23 | 复旦大学 | 基于信息物理融合系统的交通数据驱动框架及构建方法 |
CN114093164B (zh) * | 2021-11-15 | 2022-08-19 | 上海市城乡建设和交通发展研究院 | 基于车辆轨迹的动态交通流识别校正方法、装置及设备 |
CN114239929A (zh) * | 2021-11-30 | 2022-03-25 | 东南大学 | 一种基于随机森林的出租车交通需求特征预测方法 |
CN113888877B (zh) * | 2021-12-08 | 2022-03-08 | 南方科技大学 | 交通状态检测方法、装置、设备和存储介质 |
CN114495498B (zh) * | 2022-01-20 | 2023-01-10 | 青岛海信网络科技股份有限公司 | 一种交通数据分布有效性判别方法及装置 |
CN114626682A (zh) * | 2022-02-17 | 2022-06-14 | 华录智达科技股份有限公司 | 一种考虑聚集区域的城市公交线网规划方法 |
CN114217333A (zh) * | 2022-02-21 | 2022-03-22 | 北京交研智慧科技有限公司 | 路网拓扑异常位置定位方法、装置及相关设备 |
CN114822066B (zh) * | 2022-04-14 | 2023-06-13 | 北京百度网讯科技有限公司 | 车辆定位方法、装置、电子设备和存储介质 |
CN114492544B (zh) * | 2022-04-15 | 2022-07-26 | 腾讯科技(深圳)有限公司 | 模型训练方法及装置、交通事件发生概率评估方法及装置 |
CN115086910A (zh) * | 2022-04-27 | 2022-09-20 | 同济大学 | 一种v2x环境下基于动静态特征的实时路段划分方法 |
CN115019507B (zh) * | 2022-06-06 | 2023-12-01 | 上海旷途科技有限公司 | 城市路网行程时间可靠性实时估计方法 |
CN114996251A (zh) * | 2022-06-07 | 2022-09-02 | 武汉众智数字技术有限公司 | 一种基于电台gps轨迹数据的安全勤务管理方法及系统 |
CN115171372B (zh) * | 2022-06-20 | 2023-10-24 | 青岛海信网络科技股份有限公司 | 一种交通异常检测方法、设备及装置 |
CN115311846B (zh) * | 2022-06-24 | 2023-08-11 | 华东师范大学 | 一种结合货车任务状态的厂区道路拥堵预测方法及预测系统 |
CN114863685B (zh) * | 2022-07-06 | 2022-09-27 | 北京理工大学 | 一种基于风险接受程度的交通参与者轨迹预测方法及系统 |
CN115185780B (zh) * | 2022-07-21 | 2023-10-24 | 北京国联视讯信息技术股份有限公司 | 基于工业互联网的数据采集方法及系统 |
CN115440044B (zh) * | 2022-07-29 | 2023-10-13 | 深圳高速公路集团股份有限公司 | 一种公路多源事件数据融合方法、装置、存储介质及终端 |
CN115691111B (zh) * | 2022-09-22 | 2024-01-23 | 连云港杰瑞电子有限公司 | 适用于交通流数据采集的网联车最小渗透率计算方法 |
CN115661672B (zh) * | 2022-10-24 | 2023-03-14 | 中国人民解放军海军工程大学 | 基于GMM的PolSAR图像CFAR检测方法及系统 |
CN115798198B (zh) * | 2022-11-03 | 2024-04-05 | 公安部交通管理科学研究所 | 一种基于数据融合的城市路网行程时间分布估计方法 |
CN115985088B (zh) * | 2022-11-30 | 2024-01-26 | 东南大学 | 基于车辆碰撞时间反馈的交通流稳定性提升方法 |
CN116029736B (zh) * | 2023-01-05 | 2023-09-29 | 浙江警察学院 | 一种网约车异常轨迹实时检测和安全预警方法、系统 |
CN116165274B (zh) * | 2023-02-17 | 2023-11-14 | 哈尔滨工业大学 | 基于贝叶斯全局稀疏概率主成分分析的城市轨道损伤识别方法 |
CN115984077B (zh) * | 2023-02-24 | 2023-06-13 | 南方科技大学 | 一种交通异常流量因果检测方法及设备 |
CN116245362B (zh) * | 2023-03-07 | 2023-12-12 | 北京磁浮有限公司 | 城轨接触网风险评估方法及相关装置 |
CN117115759B (zh) * | 2023-04-12 | 2024-04-09 | 盐城工学院 | 一种基于类别引导的路侧交通目标检测系统及方法 |
CN116703004B (zh) * | 2023-07-19 | 2023-09-29 | 共享数据(福建)科技有限公司 | 一种基于预训练模型的水系流域智慧巡护方法和装置 |
CN116777120B (zh) * | 2023-08-16 | 2023-10-27 | 合肥综合性国家科学中心人工智能研究院(安徽省人工智能实验室) | 一种基于路网od对的城市道路移动源碳排放计算方法 |
CN116820014B (zh) * | 2023-08-24 | 2023-11-14 | 山西交通科学研究院集团有限公司 | 一种交通机电设备智能化监控预警方法及系统 |
CN116860840B (zh) * | 2023-09-05 | 2023-11-07 | 陕西交通电子工程科技有限公司 | 用于高速公路路面信息快速检索方法 |
CN117419732A (zh) * | 2023-10-10 | 2024-01-19 | 中国船舶集团有限公司第七〇九研究所 | 一种基于路网的感知目标定位纠偏方法、设备及存储介质 |
CN117517596B (zh) * | 2024-01-08 | 2024-03-15 | 辽宁中消安全设备有限公司 | 基于物联网的可燃及有毒有害气体实时监测方法及系统 |
CN117572470B (zh) * | 2024-01-15 | 2024-04-19 | 广东邦盛北斗科技股份公司 | 应用于人工智能的北斗系统定位更新方法及系统 |
CN117575546B (zh) * | 2024-01-17 | 2024-04-05 | 北京白龙马云行科技有限公司 | 一种网约车平台用后台管理系统 |
CN117706478B (zh) * | 2024-02-02 | 2024-05-03 | 腾讯科技(深圳)有限公司 | 定位漂移的识别方法、装置、设备及存储介质 |
CN117688505B (zh) * | 2024-02-04 | 2024-04-19 | 河海大学 | 一种植被大范围区域化负异常的预测方法及系统 |
CN117894181B (zh) * | 2024-03-14 | 2024-05-07 | 北京动视元科技有限公司 | 一种全域通行异常状况集成监测方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102637357A (zh) * | 2012-03-27 | 2012-08-15 | 山东大学 | 一种区域交通状态评价方法 |
CN103050005A (zh) * | 2012-11-16 | 2013-04-17 | 北京交通大学 | 城市道路交通状态时空分析方法与系统 |
CN103065466A (zh) * | 2012-11-19 | 2013-04-24 | 北京世纪高通科技有限公司 | 一种交通异常状况的检测方法和装置 |
CN104537833A (zh) * | 2014-12-19 | 2015-04-22 | 深圳大学 | 一种交通异常检测方法及系统 |
CN104809878A (zh) * | 2015-05-14 | 2015-07-29 | 重庆大学 | 利用公交车gps数据检测城市道路交通异常状态的方法 |
CN105608895A (zh) * | 2016-03-04 | 2016-05-25 | 大连理工大学 | 一种基于局部异常因子的城市交通拥堵路段检测方法 |
CN106023592A (zh) * | 2016-07-11 | 2016-10-12 | 南京邮电大学 | 一种基于gps数据的交通拥堵检测方法 |
Family Cites Families (64)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05250598A (ja) * | 1992-03-06 | 1993-09-28 | Matsushita Electric Ind Co Ltd | 交通情報提供装置 |
JPH07220193A (ja) * | 1994-01-28 | 1995-08-18 | Nagoya Denki Kogyo Kk | 交通状況判別装置 |
JP4200747B2 (ja) * | 2001-11-27 | 2008-12-24 | 富士ゼロックス株式会社 | 移動情報分類装置、移動情報分類方法、及び移動情報分類プログラム |
US7187800B2 (en) * | 2002-08-02 | 2007-03-06 | Computerized Medical Systems, Inc. | Method and apparatus for image segmentation using Jensen-Shannon divergence and Jensen-Renyi divergence |
DE102005002760B4 (de) * | 2004-01-20 | 2018-05-24 | Volkswagen Ag | Vorrichtung und Verfahren zur Unfallvermeidung bei Kraftfahrzeugen |
US7522995B2 (en) * | 2004-02-05 | 2009-04-21 | Nortrup Edward H | Method and system for providing travel time information |
CN102254434B (zh) * | 2006-03-03 | 2013-05-01 | 因瑞克斯有限公司 | 使用来自移动数据源的数据估算道路交通状况 |
US7460948B2 (en) * | 2006-03-10 | 2008-12-02 | Gm Global Technology Operations, Inc. | Traffic notification system for reporting traffic anomalies based on historical probe vehicle data |
JP4591395B2 (ja) * | 2006-03-31 | 2010-12-01 | アイシン・エィ・ダブリュ株式会社 | ナビゲーションシステム |
CN100492434C (zh) * | 2006-11-30 | 2009-05-27 | 上海交通大学 | 交通流状态分析所需的探测车采样量的获取方法 |
JP4539666B2 (ja) * | 2007-03-19 | 2010-09-08 | アイシン・エィ・ダブリュ株式会社 | 渋滞状況演算システム |
CN101373559B (zh) * | 2007-08-24 | 2010-08-18 | 同济大学 | 基于浮动车数据评估城市路网交通状态的方法 |
DE102007050154A1 (de) * | 2007-10-19 | 2009-04-23 | Siemens Ag | Prognosesystem zum Vorhersagen von Fahrzeiten, fahrzeuggestütztes Routenplanungssystem, Verkehrsinformationssystem und -verfahren |
CN101286269A (zh) * | 2008-05-26 | 2008-10-15 | 北京捷讯畅达科技发展有限公司 | 兼有动态实时交通数据的交通流量预测系统 |
CN101620781B (zh) * | 2008-06-30 | 2012-08-29 | 株式会社查纳位资讯情报 | 预测乘客信息的系统和搜索乘客信息的系统及其方法 |
CN201262784Y (zh) * | 2008-09-28 | 2009-06-24 | 华南理工大学 | 基于数据特征的城市信号控制路口交通状态检测和评价系统 |
CN101727747A (zh) * | 2009-12-16 | 2010-06-09 | 南京信息工程大学 | 基于流量检测的道路非正常拥堵报警方法 |
CN101794510A (zh) * | 2009-12-30 | 2010-08-04 | 北京世纪高通科技有限公司 | 一种浮动车数据处理的方法和装置 |
CN101950477B (zh) * | 2010-08-23 | 2012-05-23 | 北京世纪高通科技有限公司 | 一种交通信息处理方法及装置 |
CN101976505A (zh) * | 2010-10-25 | 2011-02-16 | 中国科学院深圳先进技术研究院 | 交通评价方法及系统 |
CN201946111U (zh) * | 2011-01-13 | 2011-08-24 | 西安邮电学院 | 高速公路交通异常事件预警系统用数据采集及发布节点 |
CN201927174U (zh) * | 2011-01-13 | 2011-08-10 | 西安邮电学院 | 一种高速公路交通异常事件预警系统 |
CN103477376B (zh) * | 2011-04-20 | 2016-05-25 | 日本电气株式会社 | 交通状况监视系统、方法和存储介质 |
CN202075864U (zh) * | 2011-04-28 | 2011-12-14 | 北京市劳动保护科学研究所 | 异常交通状态自动检测系统 |
CN102231235B (zh) * | 2011-04-29 | 2016-02-24 | 陈伟 | 一种交通流异常点检测定位方法 |
US8775059B2 (en) * | 2011-10-26 | 2014-07-08 | Right There Ware LLC | Method and system for fleet navigation, dispatching and multi-vehicle, multi-destination routing |
US20130166188A1 (en) * | 2011-12-21 | 2013-06-27 | Microsoft Corporation | Determine Spatiotemporal Causal Interactions In Data |
US9141874B2 (en) * | 2012-07-19 | 2015-09-22 | Qualcomm Incorporated | Feature extraction and use with a probability density function (PDF) divergence metric |
CN102855638B (zh) * | 2012-08-13 | 2015-02-11 | 苏州大学 | 基于谱聚类的车辆异常行为检测方法 |
CN103065468A (zh) * | 2012-12-14 | 2013-04-24 | 中国航天系统工程有限公司 | 交通信息的评估方法和装置 |
CN103903433B (zh) * | 2012-12-27 | 2016-09-14 | 南京中兴新软件有限责任公司 | 一种道路交通状态的实时动态判别方法及装置 |
CN103903436A (zh) * | 2012-12-28 | 2014-07-02 | 上海优途信息科技有限公司 | 一种基于浮动车的高速公路交通拥堵检测方法和系统 |
CN103258427B (zh) * | 2013-04-24 | 2015-03-11 | 北京工业大学 | 基于信息物理网络的城市快速路交通实时监控方法 |
CN103247177B (zh) * | 2013-05-21 | 2016-01-20 | 清华大学 | 大规模路网交通流实时动态预测系统 |
CN103309964A (zh) * | 2013-06-03 | 2013-09-18 | 广州市香港科大霍英东研究院 | 一种针对大规模交通数据的高效可视监测分析系统 |
CN103354030B (zh) * | 2013-07-29 | 2015-06-24 | 吉林大学 | 利用浮动公交车can总线信息判别道路交通状况的方法 |
CN103514743B (zh) * | 2013-09-28 | 2016-01-06 | 上海电科智能系统股份有限公司 | 一种实时指数匹配记忆区间的异常交通状态特征识别方法 |
US9582999B2 (en) * | 2013-10-31 | 2017-02-28 | Here Global B.V. | Traffic volume estimation |
CN103632546B (zh) * | 2013-11-27 | 2016-01-20 | 中国航天系统工程有限公司 | 一种基于浮动车数据的城市道路交通事故影响预测方法 |
US9240123B2 (en) * | 2013-12-13 | 2016-01-19 | Here Global B.V. | Systems and methods for detecting road congestion and incidents in real time |
CN103971521B (zh) * | 2014-05-19 | 2016-06-29 | 清华大学 | 道路交通异常事件实时检测方法及装置 |
KR101598343B1 (ko) * | 2014-09-23 | 2016-02-29 | 목원대학교 산학협력단 | 정체 시공간 패턴 자동인식 시스템 및 그 방법 |
CN104282151A (zh) * | 2014-09-30 | 2015-01-14 | 北京交通发展研究中心 | 基于高频卫星定位数据的实时浮动车路径匹配方法 |
CN104408958B (zh) * | 2014-11-11 | 2016-09-28 | 河海大学 | 一种城市动态路径行程时间预测方法 |
US9349285B1 (en) * | 2014-12-01 | 2016-05-24 | Here Global B.V. | Traffic classification based on spatial neighbor model |
CN104408924B (zh) * | 2014-12-04 | 2016-06-01 | 深圳北航新兴产业技术研究院 | 一种基于耦合隐马尔可夫模型的城市道路异常交通流检测方法 |
CN104504901B (zh) * | 2014-12-29 | 2016-06-08 | 浙江银江研究院有限公司 | 一种基于多维数据的交通异常点检测方法 |
CN104778834B (zh) * | 2015-01-23 | 2017-02-22 | 哈尔滨工业大学 | 一种基于车辆gps数据的城市道路交通拥堵判别方法 |
CN104657746B (zh) * | 2015-01-29 | 2017-09-12 | 电子科技大学 | 一种基于车辆轨迹相似性的异常检测方法 |
CN104573116B (zh) * | 2015-02-05 | 2017-11-03 | 哈尔滨工业大学 | 基于出租车gps数据挖掘的交通异常识别方法 |
KR101728219B1 (ko) * | 2015-02-23 | 2017-04-19 | 전북대학교산학협력단 | 양방향 통신을 이용한 시공간 교통량 분산 제어 방법 및 시스템 |
US11482100B2 (en) * | 2015-03-28 | 2022-10-25 | Intel Corporation | Technologies for detection of anomalies in vehicle traffic patterns |
CN104778837B (zh) * | 2015-04-14 | 2017-12-05 | 吉林大学 | 一种道路交通运行态势多时间尺度预测方法 |
CN104809787B (zh) * | 2015-04-23 | 2017-11-17 | 中电科安(北京)科技股份有限公司 | 一种基于摄像头的智能客流量统计装置 |
US9576481B2 (en) * | 2015-04-30 | 2017-02-21 | Here Global B.V. | Method and system for intelligent traffic jam detection |
CN105005760B (zh) * | 2015-06-11 | 2018-04-24 | 华中科技大学 | 一种基于有限混合模型的行人再识别方法 |
CN105261212B (zh) * | 2015-09-06 | 2018-06-19 | 中山大学 | 一种基于出租车gps数据地图匹配的出行时空分析方法 |
CN105404890B (zh) * | 2015-10-13 | 2018-10-16 | 广西师范学院 | 一种顾及轨迹时空语义的犯罪团伙判别方法 |
CN105513350A (zh) * | 2015-11-30 | 2016-04-20 | 华南理工大学 | 基于时空特性的分时段多参数短时交通流预测方法 |
CN105489008B (zh) * | 2015-12-28 | 2018-10-19 | 北京握奇智能科技有限公司 | 基于浮动车卫星定位数据的城市道路拥堵计算方法及系统 |
CN105761488B (zh) * | 2016-03-30 | 2018-11-23 | 湖南大学 | 基于融合的实时极限学习机短时交通流预测方法 |
CN106067248B (zh) * | 2016-05-30 | 2018-08-24 | 重庆大学 | 一种考虑速度离散特性的高速公路交通状态估计方法 |
CN106228808B (zh) * | 2016-08-05 | 2019-04-30 | 北京航空航天大学 | 基于浮动车时空网格数据的城市快速路旅行时间预测方法 |
CN106781468B (zh) * | 2016-12-09 | 2018-06-15 | 大连理工大学 | 基于建成环境和低频浮动车数据的路段行程时间估计方法 |
-
2016
- 2016-12-30 WO PCT/IB2016/058105 patent/WO2018122585A1/zh active Application Filing
- 2016-12-30 GB GBGB1711408.3A patent/GB201711408D0/en not_active Ceased
- 2016-12-30 CN CN201680088595.8A patent/CN109923595B/zh active Active
-
2017
- 2017-12-30 CN CN201780050755.4A patent/CN109643485B/zh active Active
- 2017-12-30 GB GBGB1909409.3A patent/GB201909409D0/en not_active Ceased
- 2017-12-30 GB GB1905907.0A patent/GB2569924B/en active Active
- 2017-12-30 GB GB2009834.9A patent/GB2582532B/en active Active
- 2017-12-30 GB GB2100341.3A patent/GB2587588B/en active Active
- 2017-12-30 GB GBGB1909407.7A patent/GB201909407D0/en not_active Ceased
- 2017-12-30 CN CN201780050906.6A patent/CN109716414B/zh active Active
- 2017-12-30 GB GBGB1909406.9A patent/GB201909406D0/en not_active Ceased
- 2017-12-30 WO PCT/IB2017/058532 patent/WO2018122802A1/zh active Application Filing
- 2017-12-30 WO PCT/IB2017/058533 patent/WO2018122803A1/zh active Application Filing
- 2017-12-30 CN CN201780050719.8A patent/CN110168520A/zh active Pending
- 2017-12-30 CN CN201780050907.0A patent/CN109791729B/zh active Active
- 2017-12-30 WO PCT/IB2017/058535 patent/WO2018122805A1/zh active Application Filing
- 2017-12-30 WO PCT/IB2017/058536 patent/WO2018122806A1/zh active Application Filing
- 2017-12-30 GB GB1909405.1A patent/GB2572717B/en active Active
- 2017-12-30 CN CN201780050765.8A patent/CN109844832B/zh active Active
- 2017-12-30 WO PCT/IB2017/058534 patent/WO2018122804A1/zh active Application Filing
- 2017-12-30 CN CN201780050754.XA patent/CN109997179A/zh active Pending
- 2017-12-30 GB GB2100340.5A patent/GB2588556B/en active Active
- 2017-12-30 GB GBGB1909408.5A patent/GB201909408D0/en not_active Ceased
- 2017-12-30 GB GB2009833.1A patent/GB2582531B/en active Active
- 2017-12-30 WO PCT/IB2017/058531 patent/WO2018122801A1/zh active Application Filing
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102637357A (zh) * | 2012-03-27 | 2012-08-15 | 山东大学 | 一种区域交通状态评价方法 |
CN103050005A (zh) * | 2012-11-16 | 2013-04-17 | 北京交通大学 | 城市道路交通状态时空分析方法与系统 |
CN103065466A (zh) * | 2012-11-19 | 2013-04-24 | 北京世纪高通科技有限公司 | 一种交通异常状况的检测方法和装置 |
CN104537833A (zh) * | 2014-12-19 | 2015-04-22 | 深圳大学 | 一种交通异常检测方法及系统 |
CN104809878A (zh) * | 2015-05-14 | 2015-07-29 | 重庆大学 | 利用公交车gps数据检测城市道路交通异常状态的方法 |
CN105608895A (zh) * | 2016-03-04 | 2016-05-25 | 大连理工大学 | 一种基于局部异常因子的城市交通拥堵路段检测方法 |
CN106023592A (zh) * | 2016-07-11 | 2016-10-12 | 南京邮电大学 | 一种基于gps数据的交通拥堵检测方法 |
Also Published As
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109791729B (zh) | 一种基于行程时间分布的交通异常检测方法 | |
Jindal et al. | A unified neural network approach for estimating travel time and distance for a taxi trip | |
US20180233035A1 (en) | Method and filter for floating car data sources | |
Liu et al. | Intersection delay estimation from floating car data via principal curves: a case study on Beijing’s road network | |
Carli et al. | Automated evaluation of urban traffic congestion using bus as a probe | |
Kianfar et al. | Placement of roadside equipment in connected vehicle environment for travel time estimation | |
Cerqueira et al. | On evaluating floating car data quality for knowledge discovery | |
Kurkcu et al. | Data visualization tool for monitoring transit operation and performance | |
Robinson | The development and application of an urban link travel time model using data derived from inductive loop detectors | |
Bing et al. | Integrating semantic zoning information with the prediction of road link speed based on taxi GPS data | |
CN115222936A (zh) | 过期兴趣点的确定方法、装置、电子设备及存储介质 | |
Chidlovskii et al. | Traffic delay detection by mining ticket validation transactions | |
Ratneswaran et al. | Extracting potential Travel time information from raw GPS data and Evaluating the Performance of Public transit-a case study in Kandy, Sri Lanka | |
CN115798212B (zh) | 一种基于出租车轨迹的交通拥堵检测方法 | |
Xiang et al. | Network-wide performance assessment of urban traffic based on probe vehicle data | |
Komar et al. | Intelligent Method for Counting Cars from Satellite Images | |
Yuan et al. | Research on long term parking node screening technology based on the threshold of freight vehicle parking speed | |
Monaco et al. | Exploring waste-collection fleet data: challenges in a real-world use case from multiple data providers | |
Ng et al. | Urban Travel Pattern Recognition Based on Clustering Techniques Using License Plate Sensing Data | |
Pu | Bus probe based urban travel time prediction | |
Pack et al. | Transportation Energy Analytics Dashboard (TEAD) | |
Shi | Bus Arrival Time Reliability Analyses and Dynamic Prediction Model Based on Multi-source Data | |
Li | Improving Estimation Accuracy of GPS-based Arterial Travel Time using k-Nearest Neighbors Algorithm by Zheng Li | |
PUGLIESE et al. | USING GPS DATA FROM A SAMPLE OF PRIVATE CARS FOR MODELLING THE URBAN TRAFFIC | |
Wang | Link Travel Speed Estimation Using Transit GPS Data |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |