CN110288044B - 一种基于轨迹划分与优先队列的轨迹简化方法 - Google Patents
一种基于轨迹划分与优先队列的轨迹简化方法 Download PDFInfo
- Publication number
- CN110288044B CN110288044B CN201910586296.2A CN201910586296A CN110288044B CN 110288044 B CN110288044 B CN 110288044B CN 201910586296 A CN201910586296 A CN 201910586296A CN 110288044 B CN110288044 B CN 110288044B
- Authority
- CN
- China
- Prior art keywords
- track
- point
- original
- priority
- trajectory
- 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
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000001360 synchronised effect Effects 0.000 claims abstract description 23
- 238000012545 processing Methods 0.000 claims abstract description 5
- 230000009467 reduction Effects 0.000 claims description 13
- 230000008569 process Effects 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000003825 pressing Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 claims description 2
- 238000000638 solvent extraction Methods 0.000 claims description 2
- 239000002243 precursor Substances 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 6
- 238000009825 accumulation Methods 0.000 description 4
- 230000006835 compression Effects 0.000 description 4
- 238000007906 compression Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003064 k means clustering Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 241000238366 Cephalopoda Species 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- KDYFGRWQOYBRFD-UHFFFAOYSA-N succinic acid Chemical compound OC(=O)CCC(O)=O KDYFGRWQOYBRFD-UHFFFAOYSA-N 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
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
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Train Traffic Observation, Control, And Security (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种基于轨迹划分与优先队列的轨迹简化方法,包含:S1、对原始轨迹的轨迹段的平均速度进行计算,得到整个轨迹的平均速度;S2、对整个轨迹的平均速度进行聚类划分,生成类名与速度范围表;S3、对原始轨迹中的轨迹段进行标定,形成由标签组成的轨迹段集合;S4、将标定后的轨迹段集合进行子轨迹的划分,直至结果符合预先设定的阈值;S5、基于划分好的子轨迹,生成简化率及容差;S6、将参数与子轨迹段一同传入基于优先队列的轨迹简化算法中进行简化处理;S7、将每个子轨迹的简化结果进行合并。本发明避免了原始算法的在较大的轨迹集上性能衰减的问题,同时算法在空间误差、同步欧氏距离误差方面也有较好的表现,具有很高的实用价值。
Description
技术领域
本发明涉及移动计算技术领域,特别涉及一种基于轨迹划分与优先队列的轨迹简化方法。
背景技术
2019年巴塞罗那大会上***联合酋长国的阿联酋电信(Etisalat),豪掷千金直接签下了国家级5G部署大单。日本软银(Softbank))Future Stride实验室实现了三车组编组单人驾驶,充分展示了5G通讯的优势。5G通信具有高速度低延时的特性。这些特性在无人驾驶等方面有重要的作用。但是,无论是5G还是传统卫星通信,它们的成本都很高。而通用采样的方法是在移动对象上安装传感器,并且利用现有的卫星或地面基站配合电子地图技术等,在一定时间内发送目标所在的位置信息与时间信息。由于技术条件限制,所获得的位置和时间信息是离散的。而精确的轨迹,采样率往往较高(每5秒一个采样点),造成轨迹数据集非常庞大,给轨迹数据的传输、存储、分析挖掘有效信息等方面带来极大的困难。因此,有必要在给定的误差情况下对原始轨迹数据进行简化。
移动对象轨迹数据的存储与挖掘分析是移动计算领域中的重要课题。在自动驾驶及智能交通中有广泛应用。其中,移动对象轨迹简化是该技术领域中的重要组成部分。尽管人们已提出了许许多多的轨迹简化算法,但这些算法都存在一定的局限性。其中一个经典轨迹简化算法是SQUISH算法。算法SQUISH如下所示:
输入:原始轨迹Traw,优先队列大小size;
输出:简化后的轨迹Tsim;
1、初始化一个长度为size的优先队列Buffer,同时记录点与权值;
2、forPiin Traw://遍历整个原始轨迹;
3、将Pi压入Buffer最后一个点并且置π(Pi)=∞//π(Pj)轨迹点的冗余度,权值。
4、π(Pi-1)=sed(Pi-2Pi-1Pi)//将缓冲区中的倒数第二个点赋初值为该点与其前驱后继
//点形成的同步欧几里得距离;
5、if len(Buffer)>size://缓冲区已满,len()为序列的长度;
6、找出π(Pj)0<j<size最小的点;
7、del(Pj);
8、将π(Pj)0<j<size的值直接加到π(pred(Pj))andπ(succ(Pj));
9、returnTsim=Buffer;
SQUISH算法是一个在线算法,将轨迹中的点逐个实时输入。另外需要输入的是缓冲区(优先队列大小)size,这个缓冲区被称为优先队列。最初,所有采集的轨迹点都存入缓冲区中,直到缓冲区充满。一旦数据量超出缓冲区,向缓冲区中存入任何一个点都需要删除缓冲区内已有的一个点。选出的被删除的点是包含“信息”很少或不包含“信息”的点。即信息高度冗余的点。原始的基于优先队列的轨迹压缩算法,当一个点从优先队列中移除时,增加它的邻近点的优先度。这种方法是一种很好的近似。统计上来讲运动轨迹可视为一个马尔科夫链,所以一个点只受其前驱影响,只会影响它的后继。但是实际情况复杂度远比假设情况严重,随着越来越多的点从缓冲区中移除,会不可避免的出现误差传播。所以该算法缺陷也是十分明显的。过于粗糙的冗余度分配方式导致***性的误差持续积累。实际使用中在较高的简化比时该算法的性能可能相当差。更加突出的一个问题是一般的简化算法都会让用户自己定义容差范围,来给最终的结果一个指导性的误差范围而本算法则没有。不符合通常的使用习惯,使得用户自由程度相对较低。
为了解决原始基于优先队列的轨迹压缩算法中,可能出现的误差积累,在较高的简化率或面对较大的轨迹时,可能出现的性能急剧衰减,以及缺少具有指导意义的容差指定,用户的自由度较低,研发一种基于轨迹划分与优先队列的轨迹简化算法实为必要。
发明内容
本发明的目的在提供一种基于轨迹划分与优先队列的轨迹简化方法,它是针对SQUISH算法存在的不足加以改良,首先对于输入的轨迹按照其每个轨迹段的平均速度进行聚类,将原始轨迹中的速度特征,集中地反应出来,对原始轨迹的运动状态进行大致分类,得到速度范围表。再依据这张表对原始轨迹中的轨迹段进行标定,形成一条已经标定好的轨迹,然后对这个标定好的轨迹依据,依据分类回归树进行子轨迹划分。对于划分好的子轨迹进行压缩,再经合并后成为压缩结果。本发明采用了k-means聚类算法,分类回归树算法,自适应参数生成,得到更有效的基于优先队列的轨迹压缩算法。
为了达到上述目的,本发明通过以下技术方案实现:
一种基于轨迹划分与优先队列的轨迹简化方法,该方法包含以下步骤:
S1、输入的原始轨迹经过预处理后,对原始轨迹中的每个轨迹段进行求取平均速度的操作;
S2、对步骤S1的整个轨迹的平均速度进行聚类划分,生成类名与速度范围表;
S3、基于所述类名与速度范围表,对原始轨迹中的所有轨迹段的平均速度进行标定,形成由标签组成的轨迹段集合;
S4、将标定后的轨迹段集合进行子轨迹的划分,直至划分的每个子段的基尼指数符合预先设定的阈值;
S5、基于所述步骤S4中划分好的子轨迹,生成简化率及相应的容差;
S6、将所述步骤S5中自适应生成的简化率和容差参数与子轨迹段一同传入基于优先队列的轨迹简化算法中进行简化处理;
S7、将每个子轨迹的简化结果进行合并,得到最终的结果。
优选地,所述步骤S1中,预处理的过程包含:
将原始轨迹的数据采用Pi=(xi,yi,ti)的统一形式,其中i表示原始轨迹的第i个轨迹点,Pi表示第i个轨迹点的空间位置坐标,xi,yi分别表示Pi经过投影到二维平面上的横、纵坐标,ti表示该坐标点的采集时刻即时间戳;
输入的原始轨迹表示为:
优选地,所述步骤S2中,每段轨迹的平均速度为:
整个轨迹的所有轨迹段的平均速度velocity_tab记为:
优选地,所述步骤S2中,将整个轨迹的所有轨迹段的平均速度
velocity_tab采用K-means算法进行聚类划分,包含:
通过遍历一定阈值范围内的数值,当作聚类系数k值,分别计算轮廓系数,将其中轮廓系数最小的k值作为实际划分时使用的k值,在随着k值的修正,K-means算法的结果将不断修正至最优解。
优选地,所述步骤S2中,依据所述类名与速度范围表,对整个轨迹的所有轨迹段的平均速度velocity_tab进行标记生成lable_tab,形成由标签组成的轨迹段集合。
优选地,所述步骤S4中,将标定后的所述轨迹段集合采用分类回归树的方法进行子轨迹的划分,包含:通过遍历整个轨迹,基于轨迹段集合依次在每一处进行二分划分,并计算在划分前后的基尼指数的变化,选择基尼指数减小程度最大的划分方式为最佳的划分方式;
其中,基尼指数计算公式如下:
式中,Tsample表示示例轨迹;Split(i)表示在i点处划分;T1-i表示轨迹中第1个点到第i个点组成的轨迹,Ti-n表示第i个点到第n个点组成的轨迹;Gini(T1-i)以及Gini(Ti-n)为对应轨迹的基尼指数;
公式(4)是一个递归式,递归地进行子轨迹的划分,直至每个子轨迹段的基尼指数小于所述预先设定的阈值时,结束递归过程。
优选地,所述步骤S5中,基于所述步骤S4中划分好的子轨迹,依据最小描述长度原则,生成适配的参数简化率sizeε与容差errorε;其中,sizeε表示简化轨迹的轨迹点数,errorε表示最大同步欧式距离。
优选地,所述步骤S5中,进一步包含:
S51、构建简化轨迹,其中,初始的简化轨迹为轨迹的起始点至终点组成的直线段;
S52、寻找输入的原始轨迹中与简化轨迹存在最大同步欧式距离的寻找点,该最大同步欧式距离记为errorε,其中,初始值的最大同步欧式距离记法为errormax;
S53、将步骤S52中的所述寻找点加入简化轨迹中,并计算平衡参数fix:
式中,sizeε为简化轨迹的轨迹点数,sizemax为输入的原始轨迹的轨迹点数;
S54、遍历整个原始轨迹,记录得到平衡参数fix的最小值min(fix),此时对应的sizeε与errorε为生成的轨迹最佳参数。
优选地,所述步骤S6中,所述基于优先队列的轨迹简化算法包含以下过程:
S61、初始化,构造长度为输入参数sizeε的优先队列,同时存储轨迹点与优先度,优先队列起始点与终点的优先度始终为无穷大;
S62、将轨迹点不断地压入优先队列中,每加入一个点,当前的点成为优先队列的新的终点,原先终点通过计算其与相邻的前后两个点组成的轨迹段的同步欧氏距离作为原先终点优先度;
S63、判断步骤S62中计算的原先终点的优先度是否小于输入的参数errorε,若是,则直接在优先队列中删掉该原先终点,否则,保留该原先终点;
S64、在添加新的点之后,若此时优先队列的长度大于size,则在该优先队列中寻找最小优先度的点并将其删除,同时将该点的优先度依照与其前驱点与后继点的距离的反比分配给其前驱点与后继点的优先度;
S65、重复上述步骤S63和S64添加点的过程,直至处理到输入的原始轨迹的结尾;
S66、对得到的新的优先队列再次计算每个点的同步欧氏距离并与最大同步欧式距离errorε比较,删除小于errorε的点,剩余点组成的轨迹为最终的简化轨迹。
所述步骤S64中,进一步包含:
对于轨迹点Pi、Pj、Pk,如果delPj且i<j<k:
与现有技术相比,本发明的有益效果为:原始的简化算法的初始参数与数学信息紧密相连是优先队列的大小,需要使用者在一定的经验积累之后才能较为准确的提出;而本发明新的算法所需的初始参数是基尼指数阈值这是一个对轨迹划分零散程度的感性描述;本发明只要对轨迹有初步的了解便可以直接给出;此外在压缩结果的方面,本发明的ATNS算法一定程度上避免了原始算法的在较大的轨迹集上性能衰减的问题,同时算法在空间误差、同步欧氏距离误差方面也有较好的表现;总体而言,本发明所提出的基于轨迹划分与优先队列的简化算法使用效果较好,具有很高的实用价值。
附图说明
图1为本发明的基于轨迹划分与优先队列的轨迹简化方法ATNS流程图;
图2为本发明的改进算法(ATNS)与原始算法(SQUISH)的运行时间的比较示意图;
图3为本发明的改进算法(ATNS)与原始算法(SQUISH)的同步欧式距离误(SED)差的比较示意图;
图4为本发明的改进算法(ATNS)与原始算法(SQUISH)的空间误差(PED)的比较示意图;
图5为本发明例举的一种类名与速度范围表。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。
本发明公开了一种基于轨迹划分与优先队列的轨迹简化方法ATNS,该方法包含以下步骤:
S1、输入原始轨迹以及基尼指数阈值,并对经过预处理的原始轨迹进行轨迹段的平均速度的计算;
S2、对于上述计算结果进行基于K-means算法的聚类划分处理,生成类名与速度范围表;
S3、基于步骤S2中产生的类名与速度范围表,对原始轨迹中的轨迹段的平均速度进行标定,形成一组由标签组成的轨迹段集合;
S4、对于标定后的轨迹段集合,采用分类回归树的方法进行子轨迹的划分,直至划分的每个子段的基尼指数符合事先规定好的阈值;
S5、对于步骤S4中划分好的子轨迹,依据最小描述长度原则,生成合适的参数sizeε与errorε;其中,sizeε为简化轨迹的轨迹点数(也称简化率),errorε为最大同步欧式距离(也称容差)。
S6、将步骤S5自适应生成的参数与子轨迹段一同传入基于优先队列的轨迹简化算法中,进行简化处理;
S7、将每个子轨迹的简化结果进行合并,形成结果。
所述步骤S1中,数据预处理包含以下:
对于原始轨迹的数据往往没有一种统一的形式,但是为了方便计算,本算法推荐采用形如Pi=(xi,yi,ti)的形式表示,i表示原始轨迹的第i个关键点,Pi表示第i个关键点的空间位置坐标,xi,yi分别表示Pi经过投影到二维平面上的横、纵坐标;位置信息xi,yi可以依照地理知识或规定的电子地图直接得到;最后一项ti表示该坐标点的采集时刻即时间戳;由于数据采集往往源于GPS传感器本算法推荐以世界标准时间的1980年1月6日零时作为零点,之后的时间戳均是在该零点上累加而成,自此采用统一的格式才能作为算法的输入。
所述步骤S2中,所述K-means算法(k均值聚类算法)具体包含以下:通过遍历一定阈值范围(即原始轨迹中可能包含的运动方式如:步行、跑步、乘车等的数目)内的数值,当作聚类系数k值,分别计算轮廓系数,将其中轮廓系数最小的k值作为实际划分时使用的k值(即最佳聚类系数);而在随着k值的修正,K-means算法的结果也将不断修正至最优解。
对于输入的原始轨迹Traw可以表述为:
则对于整个轨迹的平均速度有:
本实施例中,对于整个轨迹的平均速度velocity_tab使用聚类方法,无论是K-means算法还是轮廓系数的计算、选优,均可简单直接地使用TensorFlow与Python的包来实现,经过这个步骤处理之后生成类似速度范围-类名的表格,如图5所示的速度值与误差阈值的表格;依照这个表格对velocity_tab进行标记生成lable_tab,形如(000111333...),即各个轨迹段标签的组合,其中的数字均为类名。
所述S4的划分子轨迹的算法具体为:
通过遍历整个轨迹,基于lable_tab(即步骤S3中通过标定之后的由标签组成的轨迹段集合),依次在每一处进行二分划分,并计算在划分前后的基尼指数的变化,选择基尼指数减小程度最大的划分方式为最佳的划分方式。
基尼指数计算公式如下:
式中,Tsample表示示例轨迹;Split(i)为在i点处划分;T1-i表示轨迹中第1个点到第i个点组成的轨迹,Ti-n同理,表示第i个点到第n个点组成的轨迹;Gini(T1-i)以及Gini(Ti-n)为对应轨迹的基尼指数;
上述公式(4)是一个递归式,则递归地进行子轨迹的划分,直至每个子轨迹段的基尼指数小于规定的阈值(即算法开始时输入的基尼指数阈值,该阈值通过对原始轨迹想要的切分程度选取0-1范围内的数字)时,结束递归过程。
所述步骤S5中,自适应参数的生成步骤如下:
S51、构建简化轨迹;其中,初始的简化轨迹为轨迹的起始点至终点组成的直线段;
S52、寻找输入的原始轨迹中与简化轨迹有着最大同步欧式距离的点,这个最大同步欧式距离记为errorε(也称容差),其中,初始值的最大同步欧式距离记法为errormax;
S65、将步骤S52中寻找的这个点加入简化轨迹中,计算平衡参数fix:
其中,sizeε为简化轨迹的轨迹点数,sizemax为输入的原始轨迹的轨迹点数。
S54、遍历整个原始轨迹,记录得到平衡参数fix的最小值min(fix),此时对应的sizeε与errorε为生成的轨迹最佳参数。
所述步骤S6中,所述基于优先队列的轨迹简化算法如下:
S61、初始化,构造长度为输入参数sizeε的优先队列,可以同时存储轨迹点与优先度,该优先队列起始点与终点的优先度始终为无穷大;
S62、将轨迹点不断地压入优先队列中,每加入一个点,该当前的点成为优先队列的新的终点,原先终点通过计算它与紧邻的前后两个点组成的轨迹段的同步欧氏距离作为原先终点优先度;
S63、判断步骤S62中计算的原先终点的优先度是否小于输入的参数errorε,若是,则直接在优先队列中删掉原先的终点,否则,保留这个点。
S64、若在添加新的点之后优先队列的长度大于参数sizeε,则在优先队列中寻找最小优先度的点将它删除,同时将该点的优先度依照与其前驱点与后继点的距离的反比分配给其前驱点与后继点的优先度。
所述步骤S64中,对于轨迹点Pi、Pj、Pk,如果delPj且i<j<k:
S65、重复上述步骤S63和S64添加点的过程,直至处理到输入的原始轨迹的结尾;
S66、对得到的新的优先队列再次计算每个点的同步欧氏距离与最大同步欧式距离errorε比较,删除小于errorε的点,剩余点组成的轨迹为最终的简化轨迹。
如图2所示为本发明的改进算法(ATNS)与原始算法(SQUISH)的运行时间的比较示意图,横坐标为轨迹,纵坐标为运行时间。由图可见,在运行时间的数量级上两者差异并不大而改进算法相对有优势。
如图3所示为本发明的改进算法(ATNS)与原始算法(SQUISH)的同步欧式距离误差(SED)的比较示意图,横坐标为轨迹,纵坐标为同步欧式距离误差。由图可见,两个算法随数据集的增大,误差积累也在增大但是改进算法表现更优。
如图4所示为本发明的改进算法(ATNS)与原始算法(SQUISH)的空间误差(PED)的比较示意图,横坐标为轨迹,纵坐标为空间误差。由图可见,两个算法随数据集的增大,误差积累也在增大但是改进算法表现更优。
综上所述,本发明的ATNS算法一定程度上避免了原始算法的在较大的轨迹集上性能衰减的问题,同时算法在空间误差、同步欧氏距离误差方面也有较好的表现;总体而言,本发明所提出的基于轨迹划分与优先队列的简化算法使用效果较好,具有很高的实用价值。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (2)
1.一种基于轨迹划分与优先队列的轨迹简化方法,其特征在于,该方法包含以下步骤:
S1、输入的原始轨迹经过预处理后,对原始轨迹中的每个轨迹段进行求取平均速度的操作;
S2、对步骤S1的整个轨迹的平均速度进行聚类划分,生成类名与速度范围表;
S3、基于所述类名与速度范围表,对原始轨迹中的所有轨迹段的平均速度进行标定,形成由标签组成的轨迹段集合;
S4、将标定后的轨迹段集合进行子轨迹的划分,直至划分的每个子段的基尼指数符合预先设定的阈值;
S5、基于所述步骤S4中划分好的子轨迹,生成简化率及相应的容差;
S6、将所述步骤S5中自适应生成的简化率和容差参数与子轨迹段一同传入基于优先队列的轨迹简化算法中进行简化处理;
S7、将每个子轨迹的简化结果进行合并,得到最终的结果;
所述步骤S1中,预处理的过程包含:
将原始轨迹的数据采用Pi=(xi,yi,ti)的统一形式,其中i表示原始轨迹的第i个轨迹点,Pi表示第i个轨迹点的空间位置坐标,xi,yi分别表示Pi经过投影到二维平面上的横、纵坐标,ti表示该坐标点的采集时刻即时间戳;
输入的原始轨迹表示为:
所述步骤S2中,每段轨迹的平均速度为:
整个轨迹的所有轨迹段的平均速度velocity_tab记为:
所述步骤S2中,将整个轨迹的所有轨迹段的平均速度velocity_tab采用K-means算法进行聚类划分,包含:
通过遍历一定阈值范围内的数值,当作聚类系数k值,分别计算轮廓系数,将其中轮廓系数最小的k值作为实际划分时使用的k值,在随着k值的修正,K-means算法的结果将不断修正至最优解;
所述步骤S2中,依据所述类名与速度范围表,对整个轨迹的所有轨迹段的平均速度velocity_tab进行标记生成lable_tab,形成由标签组成的轨迹段集合;
所述步骤S4中,将标定后的所述轨迹段集合采用分类回归树的方法进行子轨迹的划分,包含:通过遍历整个轨迹,基于轨迹段集合依次在每一处进行二分划分,并计算在划分前后的基尼指数的变化,选择基尼指数减小程度最大的划分方式为最佳的划分方式;
其中,基尼指数计算公式如下:
式中,Tsample表示示例轨迹;Split(i)表示在i点处划分;T1-i表示轨迹中第1个点到第i个点组成的轨迹,Ti-n表示第i个点到第n个点组成的轨迹;Gini(T1-i)以及Gini(Ti-n)为对应轨迹的基尼指数;
公式(4)是一个递归式,递归地进行子轨迹的划分,直至每个子轨迹段的基尼指数小于所述预先设定的阈值时,结束递归过程;
所述步骤S5中,基于所述步骤S4中划分好的子轨迹,依据最小描述长度原则,生成适配的参数简化率sizeε与容差errorε;其中,sizeε表示简化轨迹的轨迹点数,errorε表示最大同步欧式距离;
所述步骤S5中,进一步包含:
S51、构建简化轨迹,其中,初始的简化轨迹为轨迹的起始点至终点组成的直线段;
S52、寻找输入的原始轨迹中与简化轨迹存在最大同步欧式距离的寻找点,该最大同步欧式距离记为errorε,其中,初始值的最大同步欧式距离记法为errormax;
S53、将步骤S52中的所述寻找点加入简化轨迹中,并计算平衡参数fix:
式中,sizeε为简化轨迹的轨迹点数,sizemax为输入的原始轨迹的轨迹点数;
S54、遍历整个原始轨迹,记录得到平衡参数fix的最小值min(fix),此时对应的sizeε与errorε为生成的轨迹最佳参数;
所述步骤S6中,所述基于优先队列的轨迹简化算法包含以下过程:
S61、初始化,构造长度为输入参数sizeε的优先队列,同时存储轨迹点与优先度,优先队列起始点与终点的优先度始终为无穷大;
S62、将轨迹点不断地压入优先队列中,每加入一个点,当前的点成为优先队列的新的终点,原先终点通过计算其与相邻的前后两个点组成的轨迹段的同步欧氏距离作为原先终点优先度;
S63、判断步骤S62中计算的原先终点的优先度是否小于输入的参数errorε,若是,则直接在优先队列中删掉该原先终点,否则,保留该原先终点;
S64、在添加新的点之后,若此时优先队列的长度大于sizeε,则在该优先队列中寻找最小优先度的点并将其删除,同时将该点的优先度依照与其前驱点与后继点的距离的反比分配给其前驱点与后继点的优先度;
S65、重复上述步骤S63和S64添加点的过程,直至处理到输入的原始轨迹的结尾;
S66、对得到的新的优先队列再次计算每个点的同步欧氏距离并与最大同步欧式距离errorε比较,删除小于errorε的点,剩余点组成的轨迹为最终的简化轨迹。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910586296.2A CN110288044B (zh) | 2019-07-01 | 2019-07-01 | 一种基于轨迹划分与优先队列的轨迹简化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910586296.2A CN110288044B (zh) | 2019-07-01 | 2019-07-01 | 一种基于轨迹划分与优先队列的轨迹简化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110288044A CN110288044A (zh) | 2019-09-27 |
CN110288044B true CN110288044B (zh) | 2021-06-08 |
Family
ID=68021537
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910586296.2A Active CN110288044B (zh) | 2019-07-01 | 2019-07-01 | 一种基于轨迹划分与优先队列的轨迹简化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110288044B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112668239B (zh) * | 2020-12-30 | 2022-11-15 | 山东交通学院 | 一种基于对抗学习的混合动力卡车队列经验传授方法 |
CN113052265B (zh) * | 2021-04-26 | 2024-04-05 | 上海海事大学 | 基于特征选择的移动对象轨迹简化方法 |
CN115775285B (zh) * | 2022-12-09 | 2023-09-15 | 杭州滨电信息技术有限公司 | 基于优先队列的客户端自适应地图矢量线简化方法及*** |
CN117194591B (zh) * | 2023-06-20 | 2024-05-31 | 北京和德宇航技术有限公司 | 一种球体轨迹到平面地图的显示方法、装置、设备和介质 |
CN117874155B (zh) * | 2024-01-15 | 2024-06-11 | 中国***四川省委员会党校、四川行政学院(四川长征干部学院) | 基于速度未成熟随机值的移动对象轨迹简化方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104794295A (zh) * | 2015-04-28 | 2015-07-22 | 上海海事大学 | 一种船舶轨迹实时简化方法 |
CN104794494A (zh) * | 2015-04-30 | 2015-07-22 | 中国地质大学(武汉) | 一种可扩展快速的轨迹聚类方法 |
CN104951655A (zh) * | 2015-06-08 | 2015-09-30 | 浙江大学 | 一种基于频繁模式树的车辆运动模式挖掘方法 |
CN108931795A (zh) * | 2018-05-21 | 2018-12-04 | 千寻位置网络有限公司 | 定位设备轨迹优化和边界提取方法及装置 |
CN109238271A (zh) * | 2018-08-23 | 2019-01-18 | 中国科学院电子学研究所苏州研究院 | 一种基于时间的线路拟合方法 |
CN109739585A (zh) * | 2018-12-29 | 2019-05-10 | 广西交通科学研究院有限公司 | 基于spark集群并行化计算的交通拥堵点发现方法 |
-
2019
- 2019-07-01 CN CN201910586296.2A patent/CN110288044B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104794295A (zh) * | 2015-04-28 | 2015-07-22 | 上海海事大学 | 一种船舶轨迹实时简化方法 |
CN104794494A (zh) * | 2015-04-30 | 2015-07-22 | 中国地质大学(武汉) | 一种可扩展快速的轨迹聚类方法 |
CN104951655A (zh) * | 2015-06-08 | 2015-09-30 | 浙江大学 | 一种基于频繁模式树的车辆运动模式挖掘方法 |
CN108931795A (zh) * | 2018-05-21 | 2018-12-04 | 千寻位置网络有限公司 | 定位设备轨迹优化和边界提取方法及装置 |
CN109238271A (zh) * | 2018-08-23 | 2019-01-18 | 中国科学院电子学研究所苏州研究院 | 一种基于时间的线路拟合方法 |
CN109739585A (zh) * | 2018-12-29 | 2019-05-10 | 广西交通科学研究院有限公司 | 基于spark集群并行化计算的交通拥堵点发现方法 |
Non-Patent Citations (6)
Title |
---|
Compression of Trajectory Data: A Comprehensive Evaluation and New Approach;Jonathan Muckell等;《Geoinformatica》;20140731;第1-25页 * |
Spatiotemporal Compression Techniques for Moving Point Objects;Nirvana Meratnia等;《Advances in Database Technology-EDBT 2004》;20041231;第765-782页 * |
SQUISH: An Online Approach for GPS Trajectory Compression;Jonathan Muckell等;《Proceedings of the 2nd International Conference on Computing for Geospatial Research & Applications》;20110531;第1-8页 * |
一种基于速度的移动对象轨迹简化算法;杨彪等;《计算机应用研究》;20180930;第35卷(第9期);第2666-2670页 * |
基于分段的移动对象轨迹简化算法;张甜等;《计算机应用研究》;20180424;第1-8页 * |
基于速度约束的轨迹划分方法;王敏芬;《科技创新与应用》;20131231(第26期);第3-4页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110288044A (zh) | 2019-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110288044B (zh) | 一种基于轨迹划分与优先队列的轨迹简化方法 | |
CN107103754B (zh) | 一种道路交通状况预测方法及*** | |
US8082263B2 (en) | Method, apparatus and system for multimedia model retrieval | |
US7545986B2 (en) | Adaptive resampling classifier method and apparatus | |
CN111242180B (zh) | 一种基于轻量化卷积神经网络的图像识别方法及*** | |
CN104794171B (zh) | 标记图片地理位置信息的方法及装置 | |
US20110015860A1 (en) | Systems and methods for improved generation of textual directions based on positional information | |
CN109255042B (zh) | 基于大规模街景数据的实时图片匹配地理定位方法和*** | |
CN111667693B (zh) | 用于确定预计到达时间的方法、装置、设备及介质 | |
CN112232341B (zh) | 文本检测方法、电子设备及计算机可读介质 | |
CN112050821A (zh) | 一种车道线聚合方法 | |
CN112488304A (zh) | 一种卷积神经网络中的启发式滤波器剪枝方法和*** | |
CN115358412A (zh) | 一种基于边缘辅助车辆网络的异步联邦优化方法 | |
CN114973673A (zh) | 车路协同***中结合noma和内容缓存的任务卸载方法 | |
CN105229625A (zh) | 获取有效hmm的混合散列方案 | |
CN109284409B (zh) | 基于大规模街景数据的图片组地理定位方法 | |
CN116989816B (zh) | 一种偏航识别方法、装置及电子设备 | |
CN115292962B (zh) | 基于轨迹抽稀的路径相似度匹配方法、设备及存储介质 | |
CN110411450A (zh) | 一种用于压缩轨迹的地图匹配方法 | |
CN111595353B (zh) | 一种基于GPU和Spark混合并行计算架构的实时地图匹配方法 | |
CN115334167A (zh) | 一种角度阈值自适应的轨迹压缩方法 | |
CN115329116A (zh) | 一种基于多层特征融合的图像检索方法 | |
CN115147432A (zh) | 一种基于深度残差语义分割网络的初至拾取方法 | |
CN110475198B (zh) | 一种城市道路用户轨迹纠偏处理方法及装置 | |
CN113205118A (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 |