CN109359690A - 基于卡口数据的车辆出行轨迹识别方法 - Google Patents
基于卡口数据的车辆出行轨迹识别方法 Download PDFInfo
- Publication number
- CN109359690A CN109359690A CN201811227691.3A CN201811227691A CN109359690A CN 109359690 A CN109359690 A CN 109359690A CN 201811227691 A CN201811227691 A CN 201811227691A CN 109359690 A CN109359690 A CN 109359690A
- Authority
- CN
- China
- Prior art keywords
- matrix
- vehicle
- track
- data
- bayonet
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Evolutionary Computation (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computing Systems (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明提供一种基于卡口数据的车辆出行轨迹识别方法,采用谱聚类方法,从卡口数据的时空特性角度实现出行链的自动拆分,从而识别出单次出行轨迹的起点、终点以及途经点位,具体为,生成车辆当日的完整出行链;建立相似矩阵S;计算邻接矩阵W与度矩阵D;构建拉普拉斯矩阵L,进一步确定特征矩阵F;采用DBSCASN聚类算法进行聚类分析,得到簇划分结果;确定轨迹点的簇划分结果,生成子轨迹;该方法兼顾车辆通行的时间与空间属性,与依据行程时间进行的轨迹划分方法相比更为全面。无需预设阈值,完全由实测数据驱动实现车辆通行轨迹识别,提高轨迹识别的准确性。该分析结果可用于车辆出行行为识别、出行特征分析、交通流OD估计等。
Description
技术领域
本发明涉及一种基于卡口数据的车辆出行轨迹识别方法。
背景技术
卡口设备记录通过检测点的所有车辆以及具体时间,由卡口过车检测数据能够生成车辆出行链,这在车辆跟踪、出行行为分析中具有重要作用。而如何从出行链中提取具体每一次出行的轨迹是研究和应用的关键。中国专利CN201410256101.5《一种车辆轨迹的生成、聚合方法和装置》采用了阈值法处理两条过车记录间的时间差,判断当前点位是否为尾卡口。中国专利CN201610272669.5《基于动态阈值的识别车辆轨迹中停靠路段的方法》则提出一种基于行程时间动态阈值的分析方法,确定轨迹中的停靠路段。
时间阈值法是当前卡口轨迹识别的一种主流方式。但以时间差作为两次出行的判别依据,很容易出现误判的情况。例如,在某些包含拥堵路段、路口的出行轨迹中,容易发生因拥堵延误过多而将完整的轨迹进行多次拆分的情况。即使采用动态阈值方法,将行程时间波动情况考虑在内,但也忽略了轨迹点的空间属性特征。
发明内容
本发明的目的是提供一种基于卡口数据的车辆出行轨迹识别方法,基于卡口数据所能反映的轨迹点时空特征,对单次出行轨迹自动识别,确定当日的车辆每次出行的起点、终点、途经点位以及通过时间,解决现有技术中存在的问题。
本发明的技术解决方案是:
一种基于卡口数据的车辆出行轨迹识别方法,采用谱聚类方法,从卡口数据的时空特性角度实现出行链的自动拆分,从而识别出单次出行轨迹的起点、终点以及途经点位;具体步骤为:
S1、读取卡口设备在一日内检测的过车数据,根据过车数据中的号牌号码,对于当日路网内所有通行经过卡口位置的车辆,均生成车辆当日的完整出行链;
S2、采用高斯核函数RBF建立相似矩阵S,矩阵中的元素其中i,j∈[1,n],pi、pj为轨迹点坐标,ti、tj为通过轨迹点坐标的时间,ρp、ρt为轨迹点坐标、通过轨迹点时刻的标准差;
S3、计算邻接矩阵W与度矩阵D,邻接矩阵W=S,即其中元素wi,j=si,j;度矩阵D为对角阵,D=diag(d1,…,di,…,dn),
S4、构建拉普拉斯矩阵L,进一步确定特征矩阵F;
S5、提取特征矩阵F的每一行组成1*k维样本矩阵f,对全部n个样本采用DBSCASN聚类算法进行聚类分析,得到簇划分结果{F1,···,Fi,···,Fm};其中m为聚类产生的簇的数量,簇Fi由若干个样本矩阵组成;
S6、根据每一个簇内的样本矩阵,确定轨迹点的簇划分结果Tr={tr1,···,tri,···,trm},其中tri是由若干轨迹点x组成的序列,每一个序列即为出行链拆分后的子轨迹;子轨迹的首末轨迹点即为对应的单次出行的起、终点。
进一步地,步骤S1中,生成车辆当日的完整出行链,具体为按时序排列的过车点集合X={x1,…,xi,…,xn};其中xi=(pi,ti),pi为轨迹点坐标,根据卡口设备位置确定;ti为通过该轨迹点时刻,由对应过车数据中的过车时间确定;n为该车当日的过车数据数量。
进一步地,步骤S1中,过车数据包括设备编号、过车时间、号牌号码。
进一步地,步骤S4具体为,
S41、拉普拉斯矩阵L=D-W;
S42、对拉普拉斯矩阵做标准化处理:D-1/2LD-1/2;
S43、将标准化后的矩阵最小的k个特征值所对应的特征向量组成n*k维特征矩阵F。
本发明的有益效果是:
一、该种基于卡口数据的车辆出行轨迹识别方法,运用谱聚类方法对车辆通过的轨迹点位置以及时间的信息进行处理,自动对整条出行链进行拆分,识别出单次出行轨迹。
二、本发明通过谱聚类实现对车辆出行链进行自动拆分,能够识别单次出行的轨迹,包括起点、终点、途经轨迹点以及在轨迹各点的通过时间。该分析结果可用于车辆出行行为识别、出行特征分析、交通流OD估计等。
三、本发明兼顾车辆通行的时间与空间属性,与依据行程时间进行的轨迹划分方法相比更为全面。
四、该种基于卡口数据的车辆出行轨迹识别方法,无需预设阈值,完全由实测数据驱动实现车辆通行轨迹识别,提高轨迹识别的准确性。
附图说明
图1是本发明实施例基于卡口数据的车辆出行轨迹识别方法的流程示意图。
具体实施方式
下面结合附图详细说明本发明的优选实施例。
实施例
本发明提出一种基于卡口数据的车辆出行链拆分以及出行轨迹分析方法,采用谱聚类方法,从卡口数据的时空特性角度实现出行链的自动拆分,从而识别出单次出行轨迹的起点、终点以及途经点位;如图1,具体步骤为:
S1、读取卡口设备在一日内检测的过车数据,其中过车数据包括设备编号、过车时间、号牌号码;根据过车数据中的号牌号码,对于当日路网内所有通行经过卡口位置的车辆,均生成车辆当日的完整出行链,即按时序排列的过车点集合X={x1,…,xi,…,xn};其中xi=(pi,ti),pi为轨迹点坐标,根据卡口设备位置确定;ti为通过该轨迹点时刻,由对应过车数据中的过车时间确定;n为该车当日的过车数据数量;
S2、采用高斯核函数RBF建立相似矩阵S,矩阵中的元素其中i,j∈[1,n],ρp、ρt为轨迹点坐标、通过轨迹点时刻的标准差;
S3、邻接矩阵W=S,即其中元素wi,j=si,j;度矩阵D为对角阵,D=diag(d1,…,di,…,dn),
S4、构建拉普拉斯矩阵L;进一步确定特征矩阵F;具体的,
S41、拉普拉斯矩阵L=D-W;
S42、对拉普拉斯矩阵做标准化处理:D-1/2LD-1/2;
S43、将标准化后的矩阵最小的k个特征值所对应的特征向量组成n*k维特征矩阵F;
S5、提取F的每一行组成1*k维样本矩阵f,对全部n个样本采用DBSCASN聚类算法即具有噪声的基于密度的聚类方法进行聚类分析,得到簇划分结果{F1,···,Fi,···,Fm};其中m为聚类产生的簇的数量,簇Fi由若干个样本矩阵组成;
S6、根据每一个簇内的样本矩阵,确定轨迹点的簇划分结果Tr={tr1,···,tri,···,trm},其中tri是由若干轨迹点x组成的序列,每一个序列即为出行链拆分后的子轨迹;子轨迹的首末轨迹点即为对应的单次出行的起、终点。
该种基于卡口数据的车辆出行轨迹识别方法,运用谱聚类方法对车辆通过的轨迹点位置以及时间的信息进行处理,自动对整条出行链进行拆分,识别出单次出行轨迹。实施例方法通过谱聚类实现对车辆出行链进行自动拆分,能够识别单次出行的轨迹,包括起点、终点、途经轨迹点以及在轨迹各点的通过时间。该分析结果可用于车辆出行行为识别、出行特征分析、交通流OD估计等。
该种基于卡口数据的车辆出行轨迹识别方法,兼顾车辆通行的时间与空间属性,与依据行程时间进行的轨迹划分方法相比更为全面。实施例方法,无需预设阈值,完全由实测数据驱动实现车辆通行轨迹识别,提高轨迹识别的准确性。
Claims (4)
1.一种基于卡口数据的车辆出行轨迹识别方法,其特征在于:采用谱聚类方法,从卡口数据的时空特性角度实现出行链的自动拆分,从而识别出单次出行轨迹的起点、终点以及途经点位;具体步骤为:
S1、读取卡口设备在一日内检测的过车数据,根据过车数据中的号牌号码,对于当日路网内所有通行经过卡口位置的车辆,均生成车辆当日的完整出行链;
S2、采用高斯核函数RBF建立相似矩阵S,矩阵中的元素其中i,j∈[1,n],pi、pj为轨迹点坐标,ti、tj为通过轨迹点坐标的时间,ρp、ρt为轨迹点坐标、通过轨迹点时刻的标准差;
S3、计算邻接矩阵W与度矩阵D,邻接矩阵W=S,即其中元素wi,j=si,j;度矩阵D为对角阵,D=diag(d1,…,di,…,dn),
S4、构建拉普拉斯矩阵L,进一步确定特征矩阵F;
S5、提取特征矩阵F的每一行组成1*k维样本矩阵f,对全部n个样本采用DBSCASN聚类算法进行聚类分析,得到簇划分结果{F1,…,Fi,…,Fm};其中m为聚类产生的簇的数量,簇Fi由若干个样本矩阵组成;
S6、根据每一个簇内的样本矩阵,确定轨迹点的簇划分结果Tr={tr1,…,tri,…,trm},其中tri是由若干轨迹点x组成的序列,每一个序列即为出行链拆分后的子轨迹;子轨迹的首末轨迹点即为对应的单次出行的起、终点。
2.如权利要求1所述的基于卡口数据的车辆出行轨迹识别方法,其特征在于:步骤S1中,生成车辆当日的完整出行链,具体为按时序排列的过车点集合X={x1,…,xi,…,xn};其中xi=(pi,ti),pi为轨迹点坐标,根据卡口设备位置确定;ti为通过该轨迹点时刻,由对应过车数据中的过车时间确定;n为该车当日的过车数据数量。
3.如权利要求1或2所述的基于卡口数据的车辆出行轨迹识别方法,其特征在于:步骤S1中,过车数据包括设备编号、过车时间、号牌号码。
4.如权利要求1或2所述的基于卡口数据的车辆出行轨迹识别方法,其特征在于:步骤S4具体为,
S41、拉普拉斯矩阵L=D-W;
S42、对拉普拉斯矩阵做标准化处理:D-1/2LD-1/2;
S43、将标准化后的矩阵最小的k个特征值所对应的特征向量组成n*k维特征矩阵F。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811227691.3A CN109359690B (zh) | 2018-10-19 | 2018-10-19 | 基于卡口数据的车辆出行轨迹识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811227691.3A CN109359690B (zh) | 2018-10-19 | 2018-10-19 | 基于卡口数据的车辆出行轨迹识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109359690A true CN109359690A (zh) | 2019-02-19 |
CN109359690B CN109359690B (zh) | 2021-10-22 |
Family
ID=65346094
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811227691.3A Active CN109359690B (zh) | 2018-10-19 | 2018-10-19 | 基于卡口数据的车辆出行轨迹识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109359690B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110264719A (zh) * | 2019-06-27 | 2019-09-20 | 武汉理工大学 | 一种基于多源数据的机动车动态od矩阵估计方法 |
CN110598999A (zh) * | 2019-08-21 | 2019-12-20 | 广东方纬科技有限公司 | 基于个体数据的交通出行分析方法、***及存储介质 |
CN111009123A (zh) * | 2019-11-20 | 2020-04-14 | 安徽百诚慧通科技有限公司 | 一种基于prefixspan算法的车辆频繁轨迹挖掘方法及*** |
CN111368943A (zh) * | 2020-05-27 | 2020-07-03 | 腾讯科技(深圳)有限公司 | 图像中对象的识别方法和装置、存储介质及电子装置 |
CN111696343A (zh) * | 2019-03-12 | 2020-09-22 | 北京嘀嘀无限科技发展有限公司 | 一种轨迹数据处理方法及装置 |
CN112309126A (zh) * | 2020-10-30 | 2021-02-02 | 杭州海康威视数字技术股份有限公司 | 车牌检测方法、装置、电子设备及计算机可读存储介质 |
CN112309129A (zh) * | 2020-11-03 | 2021-02-02 | 广州市交通规划研究院 | 一种基于卡口车牌识别的车辆出行特征分类的od矩阵形成方法 |
CN112348265A (zh) * | 2020-11-10 | 2021-02-09 | 交控科技股份有限公司 | 监控场景下的可行路径挖掘方法及装置 |
CN115273476A (zh) * | 2022-08-09 | 2022-11-01 | 公安部交通管理科学研究所 | 确定车辆通行起始位置位于目标区域的方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102855638A (zh) * | 2012-08-13 | 2013-01-02 | 苏州大学 | 基于谱聚类的车辆异常行为检测方法 |
CN103605362A (zh) * | 2013-09-11 | 2014-02-26 | 天津工业大学 | 基于车辆轨迹多特征的运动模式学习及异常检测方法 |
CN104504897A (zh) * | 2014-09-28 | 2015-04-08 | 北京工业大学 | 一种基于轨迹数据的交叉***通流特性分析及车辆运动预测方法 |
CN107240264A (zh) * | 2017-06-26 | 2017-10-10 | 安徽四创电子股份有限公司 | 一种车辆非有效行驶轨迹识别方法和城市道路设施规划方法 |
-
2018
- 2018-10-19 CN CN201811227691.3A patent/CN109359690B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102855638A (zh) * | 2012-08-13 | 2013-01-02 | 苏州大学 | 基于谱聚类的车辆异常行为检测方法 |
CN103605362A (zh) * | 2013-09-11 | 2014-02-26 | 天津工业大学 | 基于车辆轨迹多特征的运动模式学习及异常检测方法 |
CN104504897A (zh) * | 2014-09-28 | 2015-04-08 | 北京工业大学 | 一种基于轨迹数据的交叉***通流特性分析及车辆运动预测方法 |
CN107240264A (zh) * | 2017-06-26 | 2017-10-10 | 安徽四创电子股份有限公司 | 一种车辆非有效行驶轨迹识别方法和城市道路设施规划方法 |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111696343B (zh) * | 2019-03-12 | 2022-04-05 | 北京嘀嘀无限科技发展有限公司 | 一种轨迹数据处理方法及装置 |
CN111696343A (zh) * | 2019-03-12 | 2020-09-22 | 北京嘀嘀无限科技发展有限公司 | 一种轨迹数据处理方法及装置 |
CN110264719A (zh) * | 2019-06-27 | 2019-09-20 | 武汉理工大学 | 一种基于多源数据的机动车动态od矩阵估计方法 |
CN110598999A (zh) * | 2019-08-21 | 2019-12-20 | 广东方纬科技有限公司 | 基于个体数据的交通出行分析方法、***及存储介质 |
CN110598999B (zh) * | 2019-08-21 | 2022-05-10 | 广东方纬科技有限公司 | 基于个体数据的交通出行分析方法、***及存储介质 |
CN111009123A (zh) * | 2019-11-20 | 2020-04-14 | 安徽百诚慧通科技有限公司 | 一种基于prefixspan算法的车辆频繁轨迹挖掘方法及*** |
CN111368943A (zh) * | 2020-05-27 | 2020-07-03 | 腾讯科技(深圳)有限公司 | 图像中对象的识别方法和装置、存储介质及电子装置 |
CN111368943B (zh) * | 2020-05-27 | 2020-09-01 | 腾讯科技(深圳)有限公司 | 图像中对象的识别方法和装置、存储介质及电子装置 |
CN112309126A (zh) * | 2020-10-30 | 2021-02-02 | 杭州海康威视数字技术股份有限公司 | 车牌检测方法、装置、电子设备及计算机可读存储介质 |
CN112309126B (zh) * | 2020-10-30 | 2022-07-05 | 杭州海康威视数字技术股份有限公司 | 车牌检测方法、装置、电子设备及计算机可读存储介质 |
CN112309129B (zh) * | 2020-11-03 | 2022-04-01 | 广州市交通规划研究院 | 一种基于卡口车牌识别的车辆出行特征分类的od矩阵形成方法 |
CN112309129A (zh) * | 2020-11-03 | 2021-02-02 | 广州市交通规划研究院 | 一种基于卡口车牌识别的车辆出行特征分类的od矩阵形成方法 |
CN112348265A (zh) * | 2020-11-10 | 2021-02-09 | 交控科技股份有限公司 | 监控场景下的可行路径挖掘方法及装置 |
CN115273476A (zh) * | 2022-08-09 | 2022-11-01 | 公安部交通管理科学研究所 | 确定车辆通行起始位置位于目标区域的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109359690B (zh) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109359690A (zh) | 基于卡口数据的车辆出行轨迹识别方法 | |
CN110364008B (zh) | 路况确定方法、装置、计算机设备和存储介质 | |
CN109002905B (zh) | 预估到达时间的方法及*** | |
CN109657355B (zh) | 一种车辆道路虚拟场景的仿真方法及*** | |
US9564048B2 (en) | Origin destination estimation based on vehicle trajectory data | |
Autey et al. | Safety evaluation of right-turn smart channels using automated traffic conflict analysis | |
WO2019114432A1 (zh) | 路况生成方法、装置、设备和存储介质 | |
CN105608456B (zh) | 一种基于全卷积网络的多方向文本检测方法 | |
US20190189005A1 (en) | Method for judging highway abnormal event | |
CN104978570B (zh) | 基于增量学习的行车视频中交通标志的检测和识别方法 | |
CN111462488A (zh) | 一种基于深度卷积神经网络和交叉口行为特征模型的交叉口安全风险评估方法 | |
CN105336164B (zh) | 基于大数据分析的错误卡口位置信息自动识别方法 | |
CN106846806A (zh) | 基于Isolation Forest的城市道路交通异常检测方法 | |
CN108663065A (zh) | 一种道路类型的识别方法及装置 | |
CN106097717A (zh) | 基于两类浮动车数据融合的信号交叉口平均通行时间估计方法 | |
CN109242038A (zh) | 一种用于标签不足情况的机器人地形分类器训练方法 | |
CN108062857A (zh) | 用于出租车乘客出行目的的预测方法 | |
CN108734008A (zh) | 基于停车记录对车辆匿名的移动轨迹数据去匿名的方法 | |
CN103871077A (zh) | 一种道路车辆监控视频中的关键帧提取方法 | |
CN111949896A (zh) | 基于轨迹大数据的轨迹异常检测方法 | |
CN106327866A (zh) | 基于rfid的车辆出行od切分方法及其*** | |
CN108710637A (zh) | 基于时空关系的出租车异常轨迹实时检测方法 | |
CN109767615A (zh) | 路网交通流关键流向和关键路径分析方法 | |
CN110175654A (zh) | 一种轨迹路标的更新方法及*** | |
Žunić et al. | Innovative GPS data anomaly detection algorithm inspired by QRS complex detection algorithms in ECG signals |
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 | ||
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. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |