CN114879234B - 一种复杂gps轨迹中的重要地点挖掘方法及装置 - Google Patents

一种复杂gps轨迹中的重要地点挖掘方法及装置 Download PDF

Info

Publication number
CN114879234B
CN114879234B CN202111200912.XA CN202111200912A CN114879234B CN 114879234 B CN114879234 B CN 114879234B CN 202111200912 A CN202111200912 A CN 202111200912A CN 114879234 B CN114879234 B CN 114879234B
Authority
CN
China
Prior art keywords
point
track
points
characteristic
index
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
CN202111200912.XA
Other languages
English (en)
Other versions
CN114879234A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202111200912.XA priority Critical patent/CN114879234B/zh
Publication of CN114879234A publication Critical patent/CN114879234A/zh
Application granted granted Critical
Publication of CN114879234B publication Critical patent/CN114879234B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42—Determining position

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明涉及本发明涉及GPS技术领域，提供了一种复杂GPS轨迹中的重要地点挖掘方法及装置。其主旨在于解决目前基于密度距离测量的方法来挖掘重要地点存在的当各个簇的局部密度不统一时，用户难以用一个全局距离半径EPS值来区分出这些密度不同的簇的技术问题。主要方案为通过获取轨迹点的索引领域IN(pi)及索引领域速度NV(pi)获取特征点cpi，然后根据特征点cpi构建索引领域图，然后根据建索引领域图内的边的时空权重w(cpi,cpj)来更新通过连连接的取特征点cpi的标签，最后把标签相同的点聚集为同一个簇，然后通过用户指定的时间阈值minDuration过滤持续时间小于minDuration的簇得到最终的重要地点。

Description

一种复杂GPS轨迹中的重要地点挖掘方法及装置
技术领域
本发明涉及GPS技术领域，提供了一种复杂GPS轨迹中的重要地点挖掘方法及装置。
背景技术
现有的重要地点挖掘方法主要分为基于速度特征的重要地点挖掘方法，基于方向特征的重要地点挖掘方法和基于概率的重要地点挖掘方法三类。基于速度特征的重要地点挖掘方法以速度为特征，将一定领域范围内速度较低的轨迹点聚类到一起，得到重要地点。然而，对于对速度不敏感的应用场景，该方法会产生误报。基于方向特征的重要地点挖掘方法首先计算出所有点的方向角，再根据相邻点的方向角改变量对轨迹点进行聚类得到结果。基于速度特征的重要地点挖掘方法虽然在一定程度上能够增加对速度不敏感的应用场景的适应性，但高度依赖方向变化，难以扩展到日常行进方向基本不变的应用场景中。基于概率的重要地点挖掘方法通过定义高斯核函数算出每个点的停止概率，然后根据停止概率将重要地点挖掘出来，但停止概率难以指定。目前的方法阈值参数较多，对用户指定的阈值参数高度敏感，同时存在挖掘准确度不高，可扩展性弱，挖掘效率较低的问题。
现有技术问题存在分析：
问题1：需要用户提前指定的阈值参数过多(5个以上)，没有先验知识或者经验的用户难以在短时间内获得最优的阈值参数组合；
问题2：现有的基于密度的聚类算法在一定条件下产生了有意义且充分的结果，但当存在不同密度的聚类时，其结果并不令人满意
对于C1上的每个物体p来说，物体p和它最近的邻居之间的距离大于02和C2之间的距离。因此，在存在不同密度的簇时，我们难以选择合适的距离半径EPS值。如果距离半径Eps值小于02到C2的距离，则C1中的某些对象被指定为噪声对象。如果距离半径Eps值大于02到C2的距离，则02不被指定为噪声对象。
总结一下：当各个簇的局部密度不统一时，用户难以用一个全局距离半径EPS值来区分出这些密度不同的簇。
问题3：存在性能问题，比如速度较慢，目前存在的方法的最快时间复杂度最快是0(nlogn)，而我们的方法的时间复杂度接近0(n)，几乎是线性。在描述算法复杂度时，经常用到0(1)，0(n)，0(logn)，0(nlogn)来表示对应复杂度程度。
问题4：没有充分考虑轨迹点的时间特征和空间特征，导致重要地点挖掘的精度较低。
发明内容
本发明的目的在于解决目前基于密度距离测量的方法来挖掘重要地点存在的当各个簇的局部密度不统一时，用户难以用一个全局距离半径EPS值来区分出这些密度不同的簇的技术问题。
为了实现上述目的本发明采用以下技术手段：
本发明提供了一种复杂GPS轨迹中的重要地点挖掘方法，其特征在于，
步骤1：获取GPS定位数据，得到轨迹TR＝{p1，p2，...，pi，...，pn}，其中轨迹点pi包含纬度：pi.Lat，经度：pi.Long和时间戳：pi.T，i＝1,2,...,n且p1.T＜p2.T＜...＜pn.T。
步骤2：根据领域索引半径indexR，获取每一个轨迹点pi的索引领域IN(pi)；
步骤3：根据索引领域IN(pi)获取每一个轨迹点pi索引领域速度NV(pi)；
步骤4：获取轨迹TR中索引领域速度NV(pi)小于等于给定速度阈值minVelocity的轨迹点，记为特征点cpi，从轨迹TR中提取特征点后得到特征点集合CPS＝{cp1，cp2,...，cpi，...，cpm}；
步骤5：根据索引领域半径indexR，获取每个特征点cpi特征索引领域CPIN(cpi)＝{cpa，cpa+1，...，cpi,...,cpb-1,cpb}；
步骤6：将一条轨迹上所有特征点cpi与在特征点cpi索引领域半径indexR内的其他特征点cpj连边，构建索引领域图；
步骤7：获取特征点cpi的在特征索引领域CPIN(cpi)的停留时间NST(cpi)，NST(cpi)＝cpb.T-cpa.T，其中，cpb.T和cpa.T分别是两个特征点cpb和cpa的时间戳。两个特征点cpi和特征点cpj间的互停留时间MST(cpi，cpj)表示为
Figure BDA0003303504090000031
Figure BDA0003303504090000032
步骤8：通过数据场理论计算特征点集合CPS中特征点cpi的电势，得到步骤6中连边的特征点cpi和特征点cpj间的相对电势RP(cpi，cpj)，然后根据相对电势RP(cpi，cpj)和互停留时间MST(cpi，cpj)得到特征点cpi和cpj间的边的时空权重w(cpi,cpj)；
步骤9：对特征点集合CPS中的每一个特征点cpi分配一个唯一的标签，从特征点cpi中选取未更新过的电势P(cpi)最大的特征点cpi，然后根据最大的时空权重w(cpi,cpj)确定对应的特征点cpj，将特征点cpi的标签更新为特征点cpj的标签，重复上述步骤直到所有特征点的标签不再改变时停止，最后把标签相同的点聚集为同一个簇，然后通过用户指定的时间阈值minDuration过滤持续时间小于minDuration的簇得到最终的重要地点，簇的持续时间定义为簇中点的最大时间戳与最小时间戳的差。
定义任意两个轨迹点pi和pj间的距离
Figure BDA0003303504090000033
其中1≤i＜j≤n，dist(pk，pk+1)定义了两个相邻的轨迹点pk和pk+1的欧式距离；
在上述技术方案的基础上，步骤2中索引领域
Figure BDA0003303504090000034
Figure BDA0003303504090000035
其中indexR是一个控制索引领域IN(pi)的半径大小的正数pk中的k就对应了索引下标k。
在上述技术方案的基础上，步骤3中，获取每一个轨迹点pi索引领域速度NV(pi)，轨迹点pi的索引领域IN(pi)＝{pa，pa+1，...，pi，...，pb-1，pb}，则
Figure BDA0003303504090000036
Figure BDA0003303504090000037
其中Dist(pa，ph)为索引领域NV(pi)中所有相邻轨迹点的欧式距离的和，pb.T-pa.T为最大时间差。
在上述技术方案的基础上，步骤5具体为：获取每个特征点cpi特征索引领域CPIN(cpi)，
Figure BDA0003303504090000041
其中indexR是一个控制特征索引领域CPIN(cpi)半径大小的正数，值与提取特征点中的索引领域半径indexR大小一致。
在上述技术方案的基础上，步骤8中：
特征点cpi的电势
Figure BDA0003303504090000042
相对电势RP(cpi，cpj)＝RP(cpi，cpj)＝-|P(cpi)-P(cpj)|。
时空权重
Figure BDA0003303504090000043
其中Dist(cpi，cpj)是顶点cpi到cpj中所有相邻轨迹点的欧式距离的和，
Figure BDA0003303504090000044
Figure BDA0003303504090000045
其中1≤a＜b≤m，dist(cpk，cpk+1)，为两个相邻的轨迹点cpk和cPk+1的欧式距离。
本发明还提供了一种复杂GPS轨迹中的重要地点挖掘装置，其特征在于，
GPS数据获取模块：获取GPS定位数据，得到轨迹TR＝{p1，p2，...，pi，...，pn}，其中轨迹点pi包含纬度：pi.Lat，经度：pi.Long和时间戳：pi.T，i＝1,2，...，n且p1.T＜p2.T＜...＜pn.T；
轨迹点索引领域模块：根据领域索引半径indexR，获取每一个轨迹点pi的索引领域IN(pi)；
轨迹点引领域速度模块：根据索引领域IN(pi)获取每一个轨迹点pi索引领域速度NV(pi)；
特征点模块：获取轨迹TR中索引领域速度NV(pi)小于等于给定速度阈值minVelocity的轨迹点，记为特征点cpi，从轨迹TR中提取特征点后得到特征点集合CPS＝{cp1，cp2，...，cpi，...，cpm}；
特征点索引领域模块：根据索引领域半径indexR，获取每个特征点cpi特征索引领域CPIN(cpi)＝{cpa，cpa+1，...，cPi，...，cpb-1，cpb}；
特征点连边模块：将一条轨迹上所有特征点cpi与在特征点cpi索引领域半径indexR内的其他特征点cpj连边；
停留时间模块：获取特征点cpi的在特征索引领域CPIN(cpi)的停留时间NST(cpi)，NST(cpi)＝cpb.T-cpa.T，其中，cpb.T和cpa.T分别是两个特征点cpb和cpa的时间戳，两个特征点cpi和特征点cpj间的互停留时间MST(cpi，cpj)表示为
Figure BDA0003303504090000051
Figure BDA0003303504090000052
时空权重模块：通过数据场理论计算特征点集合CPS中特征点cpi的电势，得到特征点连边模块中连边的特征点cpi和特征点cpj间的相对电势RP(cpi，cpj)，然后根据相对电势RP(cpi，cpj)和互停留时间MST(cpi，cpj)得到特征点cpi和cpj间的边的时空权重w(cpi,cpj)；
结果模块：对特征点集合CPS中的每一个特征点cpi分配一个唯一的标签，从特征点cpi中选取未更新过的电势P(cpi)最大的特征点cpi，然后根据最大的时空权重w(cpi,cpj)确定对应的特征点cpj，将特征点cpi的标签更新为特征点cpj的标签，重复上述步骤直到所有特征点的标签不再改变时停止，最后把标签相同的点聚集为同一个簇，然后通过用户指定的时间阈值minDuration过滤持续时间小于minDuration的簇得到最终的重要地点，簇的持续时间定义为簇中点的最大时间戳与最小时间戳的差。
在上述技术方案的基础上，轨迹点索引领域模块中索引领域
Figure BDA0003303504090000053
Figure BDA0003303504090000054
其中indexR是一个控制索引领域IN(pi)的半径大小的正数pk中的k就对应了索引下标k。
在上述技术方案的基础上，获取每一个轨迹点pi索引领域速度NV(pi)，轨迹点pi的索引领域IN(pi)＝{pa，pa+1，...，pi，...，pb-1,pb}，则
Figure BDA0003303504090000061
其中Dist(pa，pb)为索引领域NV(pi)中所有相邻轨迹点的欧式距离的和，Pb.T-pa.T为最大时间差。
在上述技术方案的基础上，获取每个特征点cpi特征索引领域CPIN(cpi)，
Figure BDA0003303504090000062
其中indexR是一个控制特征索引领域CPIN(cPi)半径大小的正数，值与提取特征点中的索引领域半径indexR大小一致。
在上述技术方案的基础上，特征点cpi的电势
Figure BDA0003303504090000063
σTR为整条轨迹TR所有相邻点pi，pi-1的欧式距离的方差；
相对电势RP(cpi，cpj)＝RP(cpi，cpj)＝-|P(cpi)-P(cpj)|，
时空权重
Figure BDA0003303504090000064
其中Dist(cpI，cpj)是顶点cpi到cpj中所有相邻轨迹点的欧式距离的和，
Figure BDA0003303504090000065
Figure BDA0003303504090000066
其中1≤a＜b≤m，dist(cpk，cpk+1)为两个相邻的轨迹点cpk和cPk+1的欧式距离。
因为本发明采用上述技术手段，因此具备以下有益效果：
(1)通过提取特征点构建索引领域图，当轨迹上的点离的比较远时，我们通过索引半径indexR得到的点间的相对距离半径Eps就比较大，当轨迹上的点距离比较近时，我们通过索引半径indexR得到的点间的相对距离半径Eps就比较小，这样我们就能动态的捕获不同的局部密度，能用一个全局参数indexR较好的区分不同的簇，解决了问题2。同时，我们减少了阈值参数，只保留了3个必要的阈值参数，方便没有先验知识的用户更方便地获得最优的阈值参数组合，解决了问题1。
(2)通过定义时空权重更新措施，充分考虑了特征点的时间和空间特征，扩展了标签传播算法LPA挖掘重要地点的能力，提高了重要地点挖掘结果精确度，解决了问题4。
(3)此外，其他的重要地点挖掘方法的时间复杂度在特定情况下最快可达O(nlogn)，而我们首先提取了特征点，大大减少了需要参与运算的轨迹点的数量，降低了计算复杂度，其次通过时空权重更新措施扩展后的标签传播算法ELPA几乎是线性时间复杂度O(n)，因此重要地点挖掘结果的效率得到进一步的提升，解决了问题3。
附图说明
图1为特征点提取过程示意图；
图2为特征点索引领域图CPING；
图3为不同局部密度的簇的示例图。
具体实施方式
以下将对本发明的实施例给出详细的说明。尽管本发明将结合一些具体实施方式进行阐述和说明，但需要注意的是本发明并不仅仅只局限于这些实施方式。相反，对本发明进行的修改或者等同替换，均应涵盖在本发明的权利要求范围当中。
另外，为了更好的说明本发明，在下文的具体实施方式中给出了众多的具体细节。本领域技术人员将理解，没有这些具体细节，本发明同样可以实施。
提取特征点
一条轨迹TR定义为TR＝{p1，P2，...,pi,...,pn}，其中轨迹点pi包含纬度pi.Lat，经度pi.Long和时间戳Pi.T，i＝1,2，...，n且p1.T＜p2.T＜...＜pn.T。
一个轨迹点pi的索引领域
Figure BDA0003303504090000071
其中indexR是一个控制索引领域IN(pi)的半径大小的正数，TR是Trajectory的缩写，这里就用TR来定义了什么是轨迹。pk中的k就对应了索引下标k。假设有n个轨迹点放到集合中，得到{p1，p2，p3，p4，p5，p6，.....，pn}，它们然后依次对应索引下标为1，2，3，4，5，6，....，n，比如现在要求第四个点p4的索引领域，也就是索引下标为4对应的点p4，假设索引半径为2，那么索引领域就是{p4-2，...，p4+2}＝{p2，...，p6}，也就是索引下标为2，3，4，5，6对应的点p2，p3，p4，p5，p6都是下标为4对应的点p4的索引领域。
一个轨迹点pi的索引领域速度NV(pi)定义为索引领域NV(pi)中所有点的欧式距离和与最大时间差的比值，最大时间差是指在索引领域NV(pi)中的最大时间戳pi.T的最大值减去最小值。
一个特征点cpi定义为索引领域速度NV(pi)小于等于给定速度阈值minVelocity的轨迹点。
如图1所示为特征点提取示意图，其中横轴和纵轴对应的纬度和经度转换到x，y平面直角坐标系的坐标值。在上图中p1到p12为轨迹点，轨迹点下方的数字(0.61，0.50……0.51，0.46)为索引领域速度，假设给定的速度阈值MinVelocity＝0.5，可以得到cpl，cp2，cp3，cp4，cp5这5个标记的特征点。
根据特征点构建特征点索引领域，从轨迹TR中提取特征点后得到特征点集合CPS＝{cp1，cp2,...，cpi，...,cpm}。
一个特征点CPi的特征索引领域CPIN(cpi)定义为
Figure BDA0003303504090000081
Figure BDA0003303504090000082
其中indexR是一个控制特征索引领域CPIN(cpi)半径大小的正数，值与提取特征点中的索引领域半径indexR大小一致。
本发明可以有两种实施方式，一种是直接用特征点cpi特征点索引领域图CP ING，第二种是将cpi看成顶点vi，现对第二种方法进行如下说明：
将一条轨迹上所有特征点cpi看成顶点vi，即vi＝cpi，然后根据indexR将顶点vi与其特征索引领域中的其他顶点连边，这样我们就可以得到对应的特征点索引领域图CPING。将轨迹数据建模为索引领域图可以隔离在时间和空间上彼此距离太远的点，降低计算复杂度，并保存轨迹中重要的结构信息。特征点索引领域图CPING如图2所示。
为方便本领域技术人员理解本发明的技术构思，对根据indexR将顶点vi与其特征索引领域中的其他顶点连边的具体实施细节做进一步表述，
如图2中，在一个特征点集合CPS＝{cp1，cp2,cp3,cp4,cp5}，则对应有{v1，v2，v3，v4，v5}5个顶点，此时特征索引领域CPIN(vi)领域半径indexR＝3，然后根据indexR将特征索引领域内的顶点与在其领域半径indexR＝3内的点进行连接，
如v1左边没有点，如v1右边有点v2，v3，v4，则连接v1v2，v1v3，v1v4；
如v2左边有点v1，如v2右边有点v3，v4，v5，则连接v1v2，v2v3，v2v4，v2v5；
如v3左边有点v1，v2，如v3右边有点v4，v5，则连接v1v3，v2v3，v3v4，v3v4；
如v4左边有点v1，v2，v3，如v4右边有点v5，则连接v1v4，v2v4，v3v4，v4v5；
如v5左边有点v1，v2，v3，v4，如v5右边没有点，则连接v1v5，v2v5，v3v5，v4v5；
上述涉及的连接关系，连接一次即可，本领域技术人员也应当得知，也可二连接，不影响效果。从上述例子可以看出，每个顶点对应的领域半径indexR内，仅取存在的顶点连接，如v1左边没有顶点，则不用连接。
定义时空权重更新措施
从时间角度，我们首先定义顶点vi的领域停留时间NST(vi)，其中NST(vi)为vi特征索引领域中的最大时间差，一开始轨迹点是按照时间数据排序的，假设vi的索引领域中的特征点为cp1，cp2，...，cp7，那么这cp1到cp7也是按照时间顺序排序的，即cp1对应的时间最早，cp7对应的时间最晚，那么最大时间差为cp7.T-cp1.T，其中时间戳cpi.T。接着，两个顶点vi和vj间的互停留时间
Figure BDA0003303504090000091
Figure BDA0003303504090000092
从空间角度，给定vi的特征索引领域CPIN(vi)＝{vl,vl+1，...，vj-1，vj}，我们通过数据场理论定义点vi的电势
Figure BDA0003303504090000101
其中dist(vi，vi+1)为两个顶点的欧式距离，σTR为整条轨迹TR所有相邻点特征点(cpi，cpi+1)，即对应的相邻的顶点(vi，vi+1)的欧式距离的方差。接着，顶点vi与其索引领域半径indexR内的其他vj间的相对电势RP(vi，vj)表示为RP(vi，vj)＝-|P(vi)-P(vj)|。
则有两个顶点vi和vj间的一种时空权重更新措施，即vivj连接的边的权重
Figure BDA0003303504090000102
其中Dist(vi，vj)是顶点vi到vj中所有相邻轨迹点的欧式距离的和。Dist(vi，vj)＝dist(vi,vi+1)+dist(vi+i，vi+2)+...dist(vj-1,vj)；
标签传播聚类
基于索引领域图，通过定义的时空权重更新措施扩展经典的标签传播算法LPA，使之能够更加准确高效地挖掘重要地点。首先对所有顶点分配不一样的标签，比如第一个顶点分配标签1，第二个顶点分配标签2，以此类推。然后计算所有点的电势P(vi)，并按照从小到大的顺序排序。接着选取尚未更新过的电势P(vi)最大的顶点开始更新。假设选定的顶点为vi，通过时空权重更新措施计算该点出发的所有边的权重w(vi，vj)，然后选取权重最大的边vivj，把顶点vi的标签更新为顶点vj的标签，并记录vi的标签已经更新过了。重复上述步骤直到所有点的标签不再改变时停止。最后把标签相同的点聚集为同一个簇，然后通过用户指定的时间阈值minDuration过滤持续时间小于minDuration的簇得到最终的重要地点。簇的持续时间定义为簇中点的最大时间戳与最小时间戳的差。

Claims (10)

1.一种复杂GPS轨迹中的重要地点挖掘方法，其特征在于，
步骤1：获取GPS定位数据，得到轨迹TR＝{p1，p2，...，pi，...，pn},其中轨迹点pi包含纬度：pi.Lat，经度：pi.Long和时间戳：pi.T，i＝1,2,...,n且p1.T＜p2.T＜...＜pn.T；
步骤2：根据领域索引半径indexR，获取每一个轨迹点pi的索引领域IN(pi)；
步骤3：根据索引领域IN(pi)获取每一个轨迹点pi索引领域速度NV(pi)；
步骤4：获取轨迹TR中索引领域速度NV(pi)小于等于给定速度阈值minVelocity的轨迹点，记为特征点cpi，从轨迹TR中提取特征点后得到特征点集合CPS＝{cp1，cp2,...，cpi，...，cpm}；
步骤5：根据索引领域半径indexR，获取每个特征点cpi特征索引领域CPIN( cpi)＝{cpa，cpa+1，...，cpi,...,cpb-1,cpb},；
步骤6：将一条轨迹上所有特征点cpi与在特征点cpi索引领域半径indexR内的其他特征点cpj连边,构建索引领域图；
步骤7：获取特征点cpi的在特征索引领域CPIN(cpi)的停留时间NST(cpi) ，NST(cpi)＝cpb.T-cpa.T，其中,cpb.T和cpa.T分别是两个特征点cpb和cpa的时间戳；
两个特征点cpi和特征点cpj间的互停留时间MST(cpi，cpj)表示为
Figure 717691DEST_PATH_IMAGE028
；
步骤8：通过数据场理论计算特征点集合CPS中特征点cpi的电势，得到步骤6中连边的特征点cpi和特征点cpj间的相对电势RP(cpi，cpj)，然后根据相对电势RP(cpi，cpj)和互停留时间MST(cpi，cpj)得到特征点cpi和cpj间的边的时空权重w(cpi,cpj)；
步骤9：对特征点集合CPS中的每一个特征点cpi分配一个唯一的标签，从特征点cpi中选取未更新过的电势P(cpi)最大的特征点cpi，然后根据最大的时空权重w(cpi,cpj)确定对应的特征点cpj，将特征点cpi的标签更新为特征点cpj的标签，重复上述步骤直到所有特征点的标签不再改变时停止，最后把标签相同的点聚集为同一个簇，然后通过用户指定的时间阈值minDuration过滤持续时间小于minDuration的簇得到最终的重要地点，簇的持续时间定义为簇中点的最大时间戳与最小时间戳的差。
2.根据权利要求1所述的一种复杂GPS轨迹中的重要地点挖掘方法，其特征在于，步骤2中索引领域
Figure DEST_PATH_IMAGE035
, 其中indexR是一个控制索引领域IN(pi)的半径大小的正数pk中的k就对应了索引下标k。
3.根据权利要求1所述的一种复杂GPS轨迹中的重要地点挖掘方法，其特征在于，步骤3中，获取每一个轨迹点pi索引领域速度NV(pi)，轨迹点pi的索引领域IN(pi)＝{pa，pa+1，...，pi，...，pb-1，pb}，则
Figure DEST_PATH_IMAGE039
，其中Dist(pa，ph)为索引领域NV(pi)中所有相邻轨迹点的欧式距离的和，pb.T-pa.T为最大时间差。
4.根据权利要求1所述的一种复杂GPS轨迹中的重要地点挖掘方法，其特征在于，步骤5具体为：获取每个特征点cpi特征索引领域CPIN(cpi)，
Figure 121789DEST_PATH_IMAGE042
,其中indexR是一个控制特征索引领域CPIN(cpi)半径大小的正数，值与提取特征点中的索引领域半径indexR大小一致。
5.根据权利要求1所述的一种复杂GPS轨迹中的重要地点挖掘方法，其特征在于，步骤8中：
特征点cpi的电势
Figure DEST_PATH_IMAGE045
,σTR为整条轨迹TR所有相邻点pi，pi+1的欧式距离的方差；
相对电势RP(cpi，cpj)＝RP(cpi，cpj)＝-|P(cpi)-P(cpj)| ,
时空权重
Figure DEST_PATH_BDA0003303504090000043
，其中Dist(cpi，cpj)是顶点cpi到cpj中所有相邻轨迹点的欧式距离的和，
Figure DEST_PATH_BDA0003303504090000044
Figure DEST_PATH_BDA0003303504090000045
其中1≤a＜b≤m，dist(cpk，cpk+1)，为两个相邻的轨迹点cpk和cPk+1的欧式距离。
6.一种复杂GPS轨迹中的重要地点挖掘装置，其特征在于，
GPS数据获取模块：获取GPS定位数据，得到轨迹TR＝{p1，p2，...，pi，...，pn} ,其中轨迹点pi包含纬度：pi.Lat，经度：pi.Long和时间戳：pi.T，i＝1，2，...，n且p1.T＜p2.T＜...＜Pn·T；
轨迹点索引领域模块：根据领域索引半径indexR，获取每一个轨迹点pi的索引领域IN(pi)；
轨迹点索引领域速度模块：根据索引领域IN(pi)获取每一个轨迹点pi索引领域速度NV(pi)；
特征点模块：获取轨迹TR中索引领域速度NV(pi)小于等于给定速度阈值minVelocity的轨迹点，记为特征点cpi，从轨迹TR中提取特征点后得到特征点集合CPS＝{cp1，cp2，...，cpi，...，cpm}；
特征点索引领域模块：根据索引领域半径indexR，获取每个特征点cpi特征索引领域CPIN(cpi)＝{cpa，cpa+1，...，cPi，...，cpb-1，cpb}；
特征点连边模块：将一条轨迹上所有特征点cpi与在特征点cpi索引领域半径indexR内的其他特征点cpj连边，构建索引领域图；
停留时间模块：获取特征点cpi的在特征索引领域CPIN(cpi)的停留时间NST (cpi)，NST(cpi)＝cpb.T-cpa.T，其中, cpb.T和cpa.T分别是两个特征点cpb和cpa的时间戳，两个特征点cpi和特征点cpj间的互停留时间MST(cpi，cpj)表示为
Figure 385749DEST_PATH_IMAGE028
；
时空权重模块：通过数据场理论计算特征点集合CPS中特征点cpi的电势，得到特征点连边模块中连边的特征点cpi和特征点cpj间的相对电势RP(cpi，cpj)，然后根据相对电势RP(cpi，cpj)和互停留时间MST(cpi，cpj)得到特征点cpi和cpj间的边的时空权重w(cpi,cpj)；
结果模块：对特征点集合CPS中的每一个特征点cpi分配一个唯一的标签，从特征点cpi中选取未更新过的电势P(cpi)最大的特征点cpi，然后根据最大的时空权重w(cpi,cpj)确定对应的特征点cpj，将特征点cpi的标签更新为特征点cpj的标签，重复上述步骤直到所有特征点的标签不再改变时停止，最后把标签相同的点聚集为同一个簇，然后通过用户指定的时间阈值minDuration过滤持续时间小于minDuration的簇得到最终的重要地点，簇的持续时间定义为簇中点的最大时间戳与最小时间戳的差。
7.根据权利要求6所述的一种复杂GPS轨迹中的重要地点挖掘装置，其特征在于，轨迹点索引领域模块中索引领域
Figure 311854DEST_PATH_IMAGE035
, 其中indexR是一个控制索引领域IN(pi)的半径大小的正数pk中的k就对应了索引下标k。
8.根据权利要求6所述的一种复杂GPS轨迹中的重要地点挖掘装置，其特征在于，轨迹点索引领域速度模块中获取每一个轨迹点pi索引领域速度NV(pi)，轨迹点pi的索引领域IN(pi)＝{pa，pa+1，...，pi，...，pb-1,pb}，则
Figure 27055DEST_PATH_IMAGE039
，其中Dist(pa，pb)为索引领域NV(pi)中所有相邻轨迹点的欧式距离的和，Pb.T-pa.T为最大时间差。
9.根据权利要求6所述的一种复杂GPS轨迹中的重要地点挖掘装置，其特征在于，特征点索引领域模块实现如下：
获取每个特征点cpi特征索引领域CPIN(cpi)，
Figure 560859DEST_PATH_IMAGE042
,其中indexR是一个控制特征索引领域CPIN(cpi)半径大小的正数，值与提取特征点中的索引领域半径indexR大小一致。
10.根据权利要求6所述的一种复杂GPS轨迹中的重要地点挖掘装置，其特征在于，时空权重模块实现如下：
特征点cpi的电势
Figure 206735DEST_PATH_IMAGE045
,σTR为整条轨迹TR所有相邻点pi，pi+1的欧式距离的方差；
相对电势RP(cpi，cpj)＝RP(cpi，cpj)＝-|P(cpi)-P(cpj)|,
时空权重
Figure DEST_PATH_BDA0003303504090000064
其中Dist(cpi，cpj)是顶点cpi到cpj中所有相邻轨迹点的欧式距离的和，
Figure DEST_PATH_BDA0003303504090000065
Figure DEST_PATH_BDA0003303504090000066
其中1≤a＜b≤m，dist(cpk，cpk+1)为两个相邻的轨迹点cpk和cPk+1的欧式距离。
CN202111200912.XA 2021-10-14 2021-10-14 一种复杂gps轨迹中的重要地点挖掘方法及装置 Active CN114879234B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111200912.XA CN114879234B (zh) 2021-10-14 2021-10-14 一种复杂gps轨迹中的重要地点挖掘方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111200912.XA CN114879234B (zh) 2021-10-14 2021-10-14 一种复杂gps轨迹中的重要地点挖掘方法及装置

Publications (2)

Publication Number Publication Date
CN114879234A CN114879234A (zh) 2022-08-09
CN114879234B true CN114879234B (zh) 2023-01-06

Family

ID=82667139

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111200912.XA Active CN114879234B (zh) 2021-10-14 2021-10-14 一种复杂gps轨迹中的重要地点挖掘方法及装置

Country Status (1)

Country Link
CN (1) CN114879234B (zh)

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7747257B2 (en) * 2001-02-16 2010-06-29 Motorola, Inc. GPS assistance messages in cellular communications networks and methods therefor
JP5426742B1 (ja) * 2012-10-05 2014-02-26 株式会社小松製作所 掘削機械の表示システム及び掘削機械
KR101866405B1 (ko) * 2017-11-27 2018-06-12 주식회사 포스트에피 Gps 데이터의 압축 방법, 복원 방법 및 gps 단말 장치
US11010688B2 (en) * 2017-11-30 2021-05-18 Microsoft Technology Licensing, Llc Negative sampling
CN109190052B (zh) * 2018-08-03 2022-04-01 东北大学 一种分布式环境下基于社会感知的空间索引方法
CN109459045B (zh) * 2018-09-29 2020-10-09 杭州电子科技大学 一种针对低频gps轨迹的改进交互式投票匹配方法
CN111694032A (zh) * 2020-05-09 2020-09-22 杭州电子科技大学 一种基于聚类的大规模轨迹数据的快速地图匹配方法
CN112989223B (zh) * 2021-03-22 2022-03-01 杭州电子科技大学 基于轨迹的修正与简化的电动自行车轨迹地图匹配方法
CN113434499A (zh) * 2021-05-28 2021-09-24 西安电子科技大学 一种时空数据索引构建和查询方法、系统、数据库系统

Also Published As

Publication number Publication date
CN114879234A (zh) 2022-08-09

Similar Documents

Publication Publication Date Title
CN110084272B (zh) 一种聚类地图创建方法及基于聚类地图和位置描述子匹配的重定位方法
Chen et al. City-scale map creation and updating using GPS collections
CN110738856B (zh) 一种基于移动聚类的城市交通拥堵精细识别方法
CN107133269A (zh) 基于移动目标的频繁位置轨迹生成方法及装置
CN112013862B (zh) 一种基于众包轨迹的行人路网提取及更新方法
Zhang et al. Vehicle re-identification for lane-level travel time estimations on congested urban road networks using video images
CN114612616A (zh) 建图方法、装置、电子设备和存储介质
Hashemi Automatic inference of road and pedestrian networks from spatial-temporal trajectories
CN109343134A (zh) 一种矿井瞬变电磁探测数据分析解释方法及系统
CN106971395A (zh) 一种基于点云平均背景差的三维声纳图像建模方法
CN114593739B (zh) 基于视觉检测与参考线匹配的车辆全局定位方法及装置
Tang et al. A novel method for road intersection construction from vehicle trajectory data
CN114879234B (zh) 一种复杂gps轨迹中的重要地点挖掘方法及装置
CN112668040B (zh) 一种基于语义的轨迹聚类隐私保护方法
Yoo et al. A vision-based indoor positioning systems utilizing computer aided design drawing
CN114186112A (zh) 一种基于贝叶斯优化多重信息增益探索策略的机器人导航方法
Sasaki et al. Road segment interpolation for incomplete road data
US11988522B2 (en) Method, data processing apparatus and computer program product for generating map data
Du et al. An enhanced particle filter algorithm with map information for indoor positioning system
CN115033915A (zh) 一种基于生成对抗网络的敏感标签轨迹数据差分隐私发布方法
CN114664104A (zh) 一种路网匹配方法和装置
Zou et al. Inertia mutation energy model to extract roads by crowdsourcing trajectories
CN118051740B (zh) 融合多数据源的地质构造分析方法及存储介质
Ezzat et al. Extracting road turns and intersections from crowd source GPS tracks
CN116255976B (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
GR01 Patent grant
GR01 Patent grant