CN106646518B - 基于三阶贝塞尔曲线及插值的gps轨迹数据补全方法 - Google Patents

基于三阶贝塞尔曲线及插值的gps轨迹数据补全方法 Download PDF

Info

Publication number
CN106646518B
CN106646518B CN201611019110.8A CN201611019110A CN106646518B CN 106646518 B CN106646518 B CN 106646518B CN 201611019110 A CN201611019110 A CN 201611019110A CN 106646518 B CN106646518 B CN 106646518B
Authority
CN
China
Prior art keywords
point
tracing point
tracing
completion
track
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
CN201611019110.8A
Other languages
English (en)
Other versions
CN106646518A (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.)
Beijing Venture Group Credit Information Service Co Ltd
Original Assignee
Beijing Venture Group Credit Information Service 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 Beijing Venture Group Credit Information Service Co Ltd filed Critical Beijing Venture Group Credit Information Service Co Ltd
Priority to CN201611019110.8A priority Critical patent/CN106646518B/zh
Publication of CN106646518A publication Critical patent/CN106646518A/zh
Application granted granted Critical
Publication of CN106646518B publication Critical patent/CN106646518B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/03Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
    • G01S19/05Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing aiding data

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Navigation (AREA)

Abstract

本发明提供了一种基于三阶贝塞尔曲线及插值的GPS轨迹数据补全方法,属于GPS轨迹数据的处理技术领域。本方法包括:采集GPS轨迹数据,并进行预处理,剔除不符合要求的轨迹点,然后遍历轨迹点,计算相邻轨迹点的间隔时间,筛选出需要进行补全的轨迹点对,对每个待补全的轨迹点对,分别获取起始点和终止点的控制点,为该轨迹点对和两个控制点形成的四边形构建三阶贝塞尔曲线方程并求解,计算出该轨迹点对需要补全的轨迹点集。本发明不依赖高成本的地图实体数据,实现成本低,计算逻辑简单、过程便捷,依赖性少,支持对直线和弯道GPS轨迹数据补全,且修复补全准确率较高。

Description

基于三阶贝塞尔曲线及插值的GPS轨迹数据补全方法
技术领域
本发明涉及GPS轨迹数据的处理,具体是指一种基于三阶贝塞尔曲线及插值的GPS轨迹数据补全方法。
背景技术
轨迹数据是时空环境下,通过对一个或多个移动对象运动过程的采样所获得的数据信息,包括采样点位置、采样时间、速度等,这些采样点数据信息根据采样先后顺序构成了轨迹数据。随着卫星、无线网络,以及定位设备的发展,大量移动物体的轨迹数据呈急速增长的趋势,如交通轨迹数据、动物迁徙数据、气候气流数据、人员移动数据等。对轨迹数据的研究,获取物体有关运动的未知知识,成为未来的研究热点和应用增长点。
GPS接收模块是获取轨迹数据是常用方式,例如智能手机上自带有GPS模块,可获取出行数据。当GPS模块出现使用故障或者GPS信号的接收不稳定时,将导致数据缺失。现有解决的方法有:通过对硬件和主系统的涉及改善和减少环境对GPS信号的干扰;通过某些优化方法对已经获取的数据进行处理,来补全数据。但现有存在的GPS点位修补方法,大都是基于地图匹配进行操作的,极大地依赖地图矢量支持,或者需要考虑到复杂的交通工具以及运动物体的运动规则,如乘客的出行心理,需要利用城市路网和轨道交通数据库,存在地图矢量实体依赖强、计算逻辑复杂和偶然性较大的问题,并且计算数据多,模型复杂,实现成本高。
发明内容
本发明针对移动设备在GPS信号较弱或者丢失时,采集的GPS轨迹数据发生偏移严重或者丢失,现有方法依赖地图实体数据,计算数据多且实现复杂的问题,提供了一种基于三阶贝塞尔曲线及插值的GPS轨迹数据补全方法,实现快捷的轨迹数据补全。
本发明提供的基于三阶贝塞尔曲线及插值的GPS轨迹数据补全方法,实现步骤如下:
步骤1、采集移动设备在时间T内产生的所有GPS轨迹数据,按照采集时间先后顺序将所有轨迹点进行排序,轨迹点包括五个属性:采集时间、经度、纬度、速度和运动方向。
所述的轨迹点的运动方向定义为:设定GPS方向正北为0°,与正北方向的顺时针夹角为运动方向。
步骤2、对步骤1中的GPS轨迹数据进行预处理,剔除不符合要求的轨迹点,包括剔除速度大于阈值的轨迹点、剔除速度为空或者为负数的轨迹点,以及当相邻轨迹点的相对运动速度大于阈值A时剔除相邻轨迹点中的后一个轨迹点。
步骤3、对步骤2预处理后得到的GPS轨迹数据进行遍历,计算相邻轨迹点的间隔时间,筛选出其中间隔时间大于阈值B的轨迹点对,获取需要进行补全的轨迹点对集合。
步骤4、对步骤3获取的每个待补全的轨迹点对,分别获取起始点和终止点的控制点;
对于轨迹点Pi,其控制点p′i的经度long′i和纬度lat′i计算如下:
long′i=longi+speedi·Taverage·cos(90-diri)·180/π
lat′i=lati+speedi·Taverage·sin(90-diri)·180/π
其中,longi和lati分别为第i个轨迹点的经度和纬度,speedi和diri分别为第i个轨迹点的速度和运动方向,Taverage为所有轨迹点的平均时间间隔。
步骤5、对每个待补全的轨迹点对和步骤4获取的两个控制点四个点组成的四边形构建三阶贝塞尔曲线方程;
设轨迹点对为(Pi,Pi+1),对应的控制点为p′i和p′i+1,所构建三阶贝塞尔曲线方程f(t)为:
f(t)=Pi(1-t)3+P′i3(1-t)2t+P′i+13(1-t)t2+Pi+1t3,t∈[0,1]
参数k取值为1,2,…,n,n为自然数。
步骤6、对每个待补全的轨迹点对,求解对应的三阶贝塞尔曲线方程f(t),计算出该轨迹点对需要补全的轨迹点集合C=(c1,c2,…,ci,…,cm);m为该轨迹点对需要补全的轨迹点数量。
步骤7、重复步骤4至步骤6,对所有待补全的轨迹点对进行轨迹点补全。
本发明的优点与积极效果在于:
(1)本发明区别于通用的基于地图匹配的GPS数据点位纠正补全方法,不依赖高成本的地图实体数据,实现成本低,计算逻辑简单、计算过程便捷,依赖性少,能取得较高的修复率。
(2)本发明算法时间与空间复杂度低,集成了流式计算,支持对移动对象GPS轨迹数据进行实时处理,在处理大规模数据实时处理时取得较好的效果。
(3)本发明采用了分段多项式的方法,支持对直线和弯道GPS轨迹数据补全,且修复补全准确率较高。
附图说明
图1是本发明的GPS轨迹数据补全方法的流程示意图;
图2是本发明实施例所获取的原始GPS轨迹图形;
图3是本发明实施例利用本发明方法补齐后的GPS轨迹图形。
具体实施方式
下面结合附图和实施例对本发明的技术方案进行详细说明。
本发明的基于三阶贝塞尔曲线及插值的GPS轨迹数据补全方法,整体流程如图1所示。下面结合某个Android设备产生的GPS轨迹数据,对对其进行修复补全来说明个实现步骤。
步骤1、通过Android端设备采集一段时间T内的GPS轨迹数据,并根据时间先后顺序将所有轨迹点进行排序,轨迹点包括:采集时间、经度、纬度、速度和运动方向五个属性。其中设定GPS方向正北为0°,与正北方向的顺时针夹角为运动方向。
GPS轨迹形式化描述为:TraT=(P1,P2,…,Pi,…,Pn);其中,n为轨迹点数量,Pi为第i个轨迹点,Pi表示为:Pi=(ti,longi,lati,speedi,diri);ti为采集时间,longi为经度,lati为纬度,speedi为速度,diri为运动方向。
本发明实施例通过Android端设备采集从2016-09-29 16:13:00到2016-09-2916:53:11的数据,包括242条点位数据,部分GPS轨迹点位数据如下表1所示。
表1采集的部分GPS轨迹数据
步骤2、对步骤1中获取的GPS轨迹数据进行预处理,剔除误差较大的轨迹点,以保证后续步骤不会受到这些轨迹点的影响,预处理规则如下:
1)剔除速度speedi大于阈值thresholdspeed和速度为空的轨迹点;
2)计算GPS轨迹数据中相邻轨迹点的相对运动速度Si,当Si大于阈值thresholdspeed时,剔除相邻轨迹点中的后一个轨迹点;其中相邻的第i个轨迹点和第i+1个轨迹点的相对运动速度Si的计算公式如下:
其中,r为地球的直径。ti+1、longi+1和lati+1分别为第i+1个轨迹点的采集时间、经度和纬度。当Si大于阈值thresholdspeed时,删除第i+1个轨迹点。
根据经验,将thresholdspeed设置为150km/h,通常在实际路况中运行的车辆速度都会小于该值。
本发明实施例中,对步骤1采集到的GPS轨迹数据按照预处理规则进行预处理,剔除误差较大的轨迹点,经处理需要剔除的三段数据如下。
第一段剔除的数据:
ID 纬度 经度 数据采集时间 速度(Km/h) 方向(度)
110 26.64786 106.633687 2016-09-29 16:31:10 32.93851 28.757118
111 26.64786 106.633687 2016-09-29 16:31:20 39.2611 19.671797
112 26.64786 106.633687 2016-09-29 16:31:30 29.043465 16.101349
113 26.64786 106.633687 2016-09-29 16:31:40 45.57331 10.807328
第二段剔除的数据:
ID 纬度 经度 数据采集时间 速度(Km/h) 方向(度)
190 26.64786 106.618452 2016-09-29 16:44:31 57.843052 0.786911
191 26.64786 106.618452 2016-09-29 16:44:41 40.973644 359.60925
第三段剔除的数据:
ID 纬度 经度 数据采集时间 速度(Km/h) 方向(度)
217 26.64786 106.618452 2016-09-29 16:49:01 43.838764 40.556023
218 26.64786 106.618452 2016-09-29 16:49:11 46.267464 46.9983
步骤3、筛选需要进行轨迹点补全的点位,对步骤2中预处理得到的有效GPS轨迹数据在时间维度上以长度为2的窗口进行遍历,计算相邻轨迹点的间隔时间,其中间隔时间大于阈值thresholdtime的点位之间需要进行补全;thresholdtime为两倍的所有轨迹点的平均时间间隔Taverage
获取所有需要进行点位补全的轨迹点位对集合。每个待补全的轨迹点对,包含首尾两个轨迹点,需要补全的是这两个轨迹点中间缺失的点位。
本发明实施例中,对步骤1采集到的GPS轨迹数据计算相邻轨迹点位时间间隔,并求取其平均时间间隔,得到Taverage=10s。
对步骤3剔除误差较大轨迹点后的GPS轨迹数据筛选需要做轨迹点补全的点位,对预处理后的GPS轨迹数据在时间维度上以长度为2的窗口进行遍历,计算相邻轨迹点的间隔时间,其中间隔时间大于20秒点位需要进行补全。其中第一段需要补全的轨迹数据段如下。
步骤4、对步骤3中获取的待补全点位对集合进行处理,对于点位对Pi和Pi+1,以Pi为起始点,通过speedi和diri确认其控制点p′i,控制点p′i的经度long′i和纬度lat′i计算如下:
long′i=longi+speedi·Taverage·cos(90-diri)·180/π
lat′i=lati+speedi·Taverage·sin(90-diri)·180/π
以Pi+1为终止点,通过speedi+1和diri+1确认其控制点p′i+1,控制点p′i+1的经度long′i+1和纬度lat′i+1计算如下:
long′i+1=longi+1+Speedi+1·Taverage·cos(90-diri+1)·180/π
lat′i+1=lati+1+Speedi+1·Taverage·sin(90-diri+1)·180/π
speedi+1和diri+1分别为第i+1个轨迹点的速度和运动方向。
本发明实施例对步骤3筛选出的需要做轨迹点补全的点位进行处理,对步骤3中第一段需要补全的轨迹数据进行处理,ID为109和114的点分别为起始点和终止点,计算它们对应的控制点的经度和纬度。
起始点109的控制点的经纬度为:
(26.640897 106.654088)
终止点114的控制点的经纬度为:
(26.641551 106.653832)
步骤5、对由轨迹点Pi、Pi+1和所获取的两个控制点p′i和p′i+1四个点组成的四边形构建三阶贝塞尔曲线方程f(t):
f(t)=Pi(1-t)3+P′i3(1-t)2t+P′i+13(1-t)t2+Pi+1t3,t∈[0,1]
其中:
n为正整数
本发明实施例以步骤4中得到的起始点、终止点、两个控制点组成的四边形构建三阶贝塞尔曲线方程如下:
步骤6、求解步骤5中的三届贝塞尔方程f(t),遍历参数k得到有效t值,带入方程f(t),计算出所有补齐轨迹点集合C=(c1,c2,…,ci,…,cm),m表示所有待补齐轨迹点数量。
本发明实施例根据步骤5中方程计算出所有补齐轨迹点信息如下表。
ID 纬度 经度 数据采集时间
110 26.64047 106.653315 2016-09-29 16:31:10
111 26.640816 106.653555 2016-09-29 16:31:20
112 26.641054500020003 106.65370300010001 2016-09-29 16:31:30
113 26.641516250000002 106.65390500000001 2016-09-29 16:31:40
从109到114之间补齐的点110~113如上表所示。
步骤7、重复步骤4至步骤6,对所有待补全轨迹点对进行点位插值补全。
如图2所示,为一个获取的原始GPS轨迹,利用本发明对该轨迹数据进行处理,所得结果如图3所示,从两幅图可以看出,利用本发明方法可有效去除噪声,对GPS数据点位进行纠正补全,修复准确率高。并且通过上面过程可知,本发明方法计算复杂度低,补全轨迹点的效率高。

Claims (4)

1.一种基于三阶贝塞尔曲线及插值的GPS轨迹数据补全方法,包括:
步骤1,采集移动设备在时间T内产生的所有GPS轨迹数据,按照采集时间先后顺序将所有轨迹点进行排序,轨迹点包括五个属性:采集时间、经度、纬度、速度和运动方向;
步骤2、对GPS轨迹数据进行预处理,剔除不符合要求的轨迹点,剔除速度大于阈值A的轨迹点、剔除速度为空或者为负数的轨迹点,以及当相邻轨迹点的相对运动速度大于阈值A时剔除相邻轨迹点中的后一个轨迹点;
其特征在于,
步骤3、对预处理后得到的GPS轨迹数据进行遍历,计算相邻轨迹点的间隔时间,筛选出其中间隔时间大于阈值B的轨迹点对,获取需要进行补全的轨迹点对集合;所述的阈值B取值为两倍的所有轨迹点的平均时间间隔;
步骤4、对步骤3获取的每个待补全的轨迹点对,分别获取起始点和终止点的控制点;
对于轨迹点Pi,其控制点P′i的经度long′i和纬度lat′i计算如下:
long′i=longi+speedi·Taverage·cos(90-diri)·180/π
lat′i=lati+speedi·Taverage·sin(90-diri)·180/π
其中,longi和lati分别为第i个轨迹点的经度和纬度,speedi和diri分别为第i个轨迹点的速度和运动方向,Taverage为所有轨迹点的平均时间间隔;
步骤5、对每个待补全的轨迹点对,为两个轨迹点和对应的两个控制点组成的四边形,构建三阶贝塞尔曲线方程;
设轨迹点对为(Pi,Pi+1),对应的控制点为P′i和P′i+1,所构建三阶贝塞尔曲线方程f(t)为:
f(t)=Pi(1-t)3+P′i3(1-t)2t+P′i+13(1-t)t2+Pi+1t3,t∈[0,1]
参数其中k取值为1,2,…,n,n为自然数;ti为第i个轨迹点的采集时间;ti+1为第i+1个轨迹点的采集时间;
步骤6、对每个待补全的轨迹点对,求解对应的三阶贝塞尔曲线方程,计算出该轨迹点对需要补全的轨迹点集合;
求解步骤5中的三阶贝塞尔方程f(t),遍历参数k得到有效t值,带入方程f(t),计算出所有补齐轨迹点集合C=(c1,c2,...,ci,...,cm),m表示所有待补齐轨迹点数量;
步骤7、重复步骤4至步骤6,对所有待补全的轨迹点对进行补全。
2.根据权利要求1所述的GPS轨迹数据补全方法,其特征在于,所述的轨迹点的运动方向定义为:设定GPS方向正北为0°,与正北方向的顺时针夹角为运动方向。
3.根据权利要求1所述的GPS轨迹数据补全方法,其特征在于,步骤2中所述的阈值A设置为150km/h。
4.根据权利要求1所述的GPS轨迹数据补全方法,其特征在于,步骤2中所述的相邻轨迹点的相对运动速度,记为Si,根据下式确定:
其中,r为地球的直径,ti、longi和lati分别为第i个轨迹点的采集时间、经度和纬度,ti+1、longi+1和lati+1分别为第i+1个轨迹点的采集时间、经度和纬度。
CN201611019110.8A 2016-11-18 2016-11-18 基于三阶贝塞尔曲线及插值的gps轨迹数据补全方法 Active CN106646518B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611019110.8A CN106646518B (zh) 2016-11-18 2016-11-18 基于三阶贝塞尔曲线及插值的gps轨迹数据补全方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611019110.8A CN106646518B (zh) 2016-11-18 2016-11-18 基于三阶贝塞尔曲线及插值的gps轨迹数据补全方法

Publications (2)

Publication Number Publication Date
CN106646518A CN106646518A (zh) 2017-05-10
CN106646518B true CN106646518B (zh) 2019-06-11

Family

ID=58808134

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611019110.8A Active CN106646518B (zh) 2016-11-18 2016-11-18 基于三阶贝塞尔曲线及插值的gps轨迹数据补全方法

Country Status (1)

Country Link
CN (1) CN106646518B (zh)

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109831736B (zh) * 2017-11-23 2022-01-18 腾讯科技(深圳)有限公司 一种数据处理方法、装置、服务器及客户端
CN108376415B (zh) * 2018-02-13 2022-01-21 中国联合网络通信集团有限公司 一种轨迹填充的方法及装置
CN108509576B (zh) * 2018-03-28 2021-04-20 浙江知水信息技术有限公司 一种基于电子地图的巡查轨迹错误点位的过滤优化方法
CN109031371B (zh) * 2018-06-13 2020-11-03 山东众云教育科技有限公司 一种基于持续定位信息的路线平滑方法及系统
CN109508471B (zh) * 2018-09-21 2024-04-12 长安通信科技有限责任公司 运动轨迹补全方法及装置、可读存储介质
CN109405839B (zh) * 2018-10-23 2022-04-12 南京林业大学 一种基于多路径的交通网络离线地图匹配算法
CN109726355B (zh) * 2019-01-04 2022-11-25 重庆邮电大学 一种基于向量插值的船舶轨迹修复方法
CN110727757A (zh) * 2019-10-22 2020-01-24 北京卡路里信息技术有限公司 轨迹数据的处理方法及装置、电子设备
CN113093716B (zh) * 2019-12-19 2024-04-30 广州极飞科技股份有限公司 一种运动轨迹规划方法、装置、设备及存储介质
CN111199075B (zh) * 2019-12-30 2023-04-11 四川函钛科技有限公司 一种基于时序qar参数的飞行轨迹自适应平滑的方法
CN111432347B (zh) * 2020-06-11 2020-09-08 腾讯科技(深圳)有限公司 信息处理方法、装置、存储介质及电子设备
CN113891234A (zh) * 2020-06-16 2022-01-04 大众问问(北京)信息科技有限公司 一种轨迹数据补全方法、装置、计算机设备和存储介质
CN111767361A (zh) * 2020-07-06 2020-10-13 中电万维信息技术有限责任公司 一种基于贝塞尔曲线算法的gis应用
CN112435573B (zh) * 2020-11-24 2022-03-08 北京赛目科技有限公司 一种用于自动驾驶仿真测试的地图绘制方法及装置
CN112597189B (zh) * 2020-12-28 2024-08-27 广州羊城通有限公司 基于乘车记录的出行轨迹调整方法及装置
CN112783161B (zh) * 2020-12-29 2023-04-25 广东嘉腾机器人自动化有限公司 一种基于贝塞尔曲线的agv避障方法
CN118152677B (zh) * 2024-05-10 2024-08-02 浙江大华技术股份有限公司 轨迹补全方法、终端及计算机可读存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0816643A (ja) * 1994-06-29 1996-01-19 Ricoh Co Ltd 面積計算方法
CN102750722A (zh) * 2011-08-11 2012-10-24 新奥特(北京)视频技术有限公司 一种运动模型轨迹的生成方法及系统
CN102998684A (zh) * 2012-11-21 2013-03-27 厦门雅迅网络股份有限公司 一种基于贝塞尔曲线的终端定位轨迹拟合方法
CN104335010A (zh) * 2012-09-24 2015-02-04 日本先锋公司 图像处理装置和图像处理方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0816643A (ja) * 1994-06-29 1996-01-19 Ricoh Co Ltd 面積計算方法
CN102750722A (zh) * 2011-08-11 2012-10-24 新奥特(北京)视频技术有限公司 一种运动模型轨迹的生成方法及系统
CN104335010A (zh) * 2012-09-24 2015-02-04 日本先锋公司 图像处理装置和图像处理方法
CN102998684A (zh) * 2012-11-21 2013-03-27 厦门雅迅网络股份有限公司 一种基于贝塞尔曲线的终端定位轨迹拟合方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
一种基于贝塞尔曲线的行车数据恢复策略;何振峰;《中国电子学会第十四届信息论学术年会》;20071218;第458-461页

Also Published As

Publication number Publication date
CN106646518A (zh) 2017-05-10

Similar Documents

Publication Publication Date Title
CN106646518B (zh) 基于三阶贝塞尔曲线及插值的gps轨迹数据补全方法
CN104331422B (zh) 一种路段类型推测方法
CN106969764B (zh) 一种道路匹配方法、装置及车载地图采集系统
Ahmed et al. Constructing street networks from GPS trajectories
CN106595656B (zh) 一种基于车载pos的公路线形特征识别与参数计算方法
CN102102992B (zh) 基于多级网络划分的匹配道路初筛方法及地图匹配系统
CN103499828B (zh) 浮动车轨迹地图匹配的方法及装置
CN110008872B (zh) 一种结合车辆轨迹和遥感图像的路网提取方法
CN105203116B (zh) 一种基于条件随机场和低采样频率浮动车数据的地图匹配方法
CN111210612B (zh) 基于公交gps数据与站点信息提取公交线路轨迹的方法
CN108759833A (zh) 一种基于先验地图的智能车辆定位方法
Gikas et al. A novel geodetic engineering method for accurate and automated road/railway centerline geometry extraction based on the bearing diagram and fractal behavior
CN104318781B (zh) 基于rfid技术的行程速度获取方法
CN104700617A (zh) 基于低精度gps轨迹数据的高精度车道信息提取方法
CN106469505B (zh) 一种浮动车轨迹纠偏方法和装置
CN107328423B (zh) 基于地图数据的弯道识别方法及其系统
CN108253974B (zh) 浮动车定位数据自适应缓冲路径匹配系统及其方法
CN103617328A (zh) 一种飞机三维姿态解算方法
CN103592665A (zh) 基于b样条曲线的车载终端定位轨迹拟合方法
CN106441316A (zh) 一种基于历史数据的单点路网匹配方法
CN103839274A (zh) 一种基于几何比例关系的扩展目标跟踪方法
CN104504245A (zh) 一种应用gps出行调查数据识别出行和活动的方法
CN106289281A (zh) 一种基于三证据ds理论的双模式地图匹配方法
CN108447262A (zh) 一种由行车轨迹构造“数字化道路”的技术与方法
Yu et al. Improving integrity and reliability of map matching techniques

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