CN116303479A - 一种用于轨迹数据的实时压缩方法 - Google Patents

一种用于轨迹数据的实时压缩方法 Download PDF

Info

Publication number
CN116303479A
CN116303479A CN202310572745.4A CN202310572745A CN116303479A CN 116303479 A CN116303479 A CN 116303479A CN 202310572745 A CN202310572745 A CN 202310572745A CN 116303479 A CN116303479 A CN 116303479A
Authority
CN
China
Prior art keywords
straight line
time
point
data
segmentation
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
Application number
CN202310572745.4A
Other languages
English (en)
Other versions
CN116303479B (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.)
Ningbo Zhongjing Technology Development Co ltd
Zhejiang University of Science and Technology ZUST
Original Assignee
Ningbo Zhongjing Technology Development Co ltd
Zhejiang University of Science and Technology ZUST
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 Ningbo Zhongjing Technology Development Co ltd, Zhejiang University of Science and Technology ZUST filed Critical Ningbo Zhongjing Technology Development Co ltd
Priority to CN202310572745.4A priority Critical patent/CN116303479B/zh
Publication of CN116303479A publication Critical patent/CN116303479A/zh
Application granted granted Critical
Publication of CN116303479B publication Critical patent/CN116303479B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/22—Indexing; Data structures therefor; Storage structures
    • G06F16/2228—Indexing structures
    • G06F16/2272—Management thereof
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24—Querying
    • G06F16/245—Query processing
    • G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
    • G06F16/2474—Sequence data queries, e.g. querying versioned data
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/24—Querying
    • G06F16/245—Query processing
    • G06F16/2458—Special types of queries, e.g. statistical queries, fuzzy queries or distributed queries
    • G06F16/2477—Temporal data queries
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29—Geographical information databases
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M7/00—Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
    • H03M7/30—Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
    • H03M7/3059—Digital compression and data reduction techniques where the original information is represented by a subset or similar information, e.g. lossy compression
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Computational Linguistics (AREA)
  • Remote Sensing (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

本发明公开了一种用于轨迹数据的实时压缩方法，涉及轨迹数据压缩领域，其基于分割初始时刻
Figure ZY_3
与最大误差
Figure ZY_7
利用DConnPLA算法分割各维度数据序列，得到各维度数据序列对应的可行分割直线，实现了对可行直线空间的完全搜索，另外，通过获取各可行分割直线
Figure ZY_11
中最短可行分割直线对应的分割结束时刻
Figure ZY_4
，通过可行分割直线的表达式
Figure ZY_8
计算各维度数据序列在分割初始时刻
Figure ZY_12
对应的初始值
Figure ZY_15
和在分割结束时刻
Figure ZY_1
对应的结束值
Figure ZY_5
，通过各维度数据序列对应的初始值
Figure ZY_9
和结束值
Figure ZY_13
得到目标起点和目标终点，通过目标起点和目标终点分割当前的
Figure ZY_2
维实时轨迹数据
Figure ZY_6
，并以分割结束时刻
Figure ZY_10
加1作为下一段
Figure ZY_14
维实时轨迹数据的分割初始时刻，实现了三维及三维以上空间数据的最优分割。

Description

一种用于轨迹数据的实时压缩方法
技术领域
本发明涉及轨迹数据压缩领域，尤其涉及一种用于轨迹数据的实时压缩方法。
背景技术
随着传感器和物联网技术的普及发展，在地理信息系统中，大量的轨迹数据被收集并用于特定的分析任务。近些年，轨迹数据呈现爆炸式的增长，这给其传输、存储和分析带来了前所未有的挑战。作为传统的数据存储解决方案，数据压缩技术是应对这些挑战的可行手段之一。目前，针对轨迹数据的有损压缩是主流，而无损压缩在压缩率以及压缩概要（即压缩后的数据）的可用性上存在一定的局限性。在众多有损压缩算法中，线性分段拟合(PLA)技术具有简单直观的特点，因此被广泛的研究和应用。
早期针对轨迹数据的PLA算法尽管可以获得最优化（即最小）的压缩规模，但时空复杂度较高，不适合应对流数据压缩的应用需求。而目前流行的三种次优化解决方案，即批处理算法、在线算法和实时算法，均采用在经纬度空间下定义的垂直误差或同步误差进行限定，不能在线性时间复杂度下搜索所有可行直线，所以很难获得最优化的结果。
针对轨迹数据的压缩，现有的PLA压缩技术也不能在线性时间复杂度获得最优化的压缩结果，即现有技术中是通过空间欧式距离（也称同步误差）进行的度量，只能在二维空间进行划分无法对三维及三维以上的空间完成算法处理，其原因在于，这些算法所采用的垂直误差或同步误差均不能完成对可行直线空间的完全搜索，因而无法获得最优化的分割结果，对此，本发明提出了一种用于轨迹数据的实时压缩方法。
发明内容
为了对三维及三维以上空间的数据实现最优分割，本发明提出了一种用于轨迹数据的实时压缩方法，包括步骤：
S1：获取
Figure SMS_3
维实时轨迹数据/>
Figure SMS_2
并设定最大误差/>
Figure SMS_10
，其中，
Figure SMS_8
，/>
Figure SMS_14
表示数据点/>
Figure SMS_16
的时间戳索引，/>
Figure SMS_17
为常数，/>
Figure SMS_7
表示在/>
Figure SMS_9
时刻的/>
Figure SMS_5
个轨迹观测值；获取/>
Figure SMS_15
维实时轨迹数据对应的维度数据序列
Figure SMS_1
，其中，/>
Figure SMS_12
表示第/>
Figure SMS_6
维的数据序列，/>
Figure SMS_11
为/>
Figure SMS_4
维/>
Figure SMS_13
时刻的轨迹观测值；
S2：初始化分割序号
Figure SMS_18
与分割初始时刻/>
Figure SMS_19
；
S3：基于分割初始时刻
Figure SMS_20
与最大误差/>
Figure SMS_21
利用DConnPLA算法分别分割各维度数据序列，得到各维度数据序列对应的可行分割直线，所述可行分割直线的表达式为/>
Figure SMS_22
，其中
Figure SMS_23
表示时刻；
S4：获取各可行分割直线
Figure SMS_25
中最短可行分割直线对应的分割结束时刻/>
Figure SMS_29
，通过可行分割直线的表达式/>
Figure SMS_32
计算各维度数据序列在分割初始时刻/>
Figure SMS_26
对应的初始值/>
Figure SMS_28
和在分割结束时刻/>
Figure SMS_31
对应的结束值/>
Figure SMS_34
，通过各维度数据序列对应的初始值/>
Figure SMS_24
和结束值/>
Figure SMS_27
得到目标起点/>
Figure SMS_30
和目标终点/>
Figure SMS_33
；
S5：通过目标起点
Figure SMS_35
和目标终点/>
Figure SMS_36
分割当前的/>
Figure SMS_37
维实时轨迹数据/>
Figure SMS_38
，设定
Figure SMS_39
、/>
Figure SMS_40
，并返回S3步骤继续下一段/>
Figure SMS_41
维实时轨迹数据的分割。
进一步地，所述S2步骤中，分割序号
Figure SMS_42
与分割初始时刻/>
Figure SMS_43
的初始值均等于1；
所述S3步骤中，利用DConnPLA算法分割维度数据序列，具体包括：
S31：初始化时间序号
Figure SMS_44
，初始化维度数据序列对应的上凸壳、下凸壳以及第一直线与第二直线的斜率：
所述时间序号
Figure SMS_45
的初始化公式为：/>
Figure SMS_46
；
所述上凸壳的表达式为：
Figure SMS_47
；
所述下凸壳的表达式为：
Figure SMS_48
；
所述第一直线的表达式为：
Figure SMS_49
；
所述第二直线的表达式为：
Figure SMS_50
；
所述第一直线对应斜率的获取公式为：
Figure SMS_51
；
所述第二直线对应斜率的获取公式为：
Figure SMS_52
；其中：
Figure SMS_53
，/>
Figure SMS_54
；
Figure SMS_55
，/>
Figure SMS_56
；
式中，
Figure SMS_62
表示分割初始时刻/>
Figure SMS_58
对应的时间点，/>
Figure SMS_66
表示第/>
Figure SMS_61
维数据序列中/>
Figure SMS_69
时间点对应的轨迹观测值，/>
Figure SMS_70
表示最大误差，/>
Figure SMS_73
、/>
Figure SMS_63
分别表示第一直线的上容忍点与下容忍点，/>
Figure SMS_65
、/>
Figure SMS_64
分别表示第二直线的上容忍点与下容忍点，/>
Figure SMS_71
表示第一直线的斜率，/>
Figure SMS_60
表示第二直线的斜率，/>
Figure SMS_67
表示上凸壳，/>
Figure SMS_57
表示下凸壳；/>
Figure SMS_72
表示上一段第/>
Figure SMS_59
维数据序列分割时输出的时间序号/>
Figure SMS_68
；
S32：通过维度数据序列中轨迹观测值
Figure SMS_74
对应的时间点/>
Figure SMS_75
与第一直线获取第一数据值，通过轨迹观测值/>
Figure SMS_76
对应的时间点/>
Figure SMS_77
与第二直线获取第二数据值；
S33：判断轨迹观测值
Figure SMS_78
是否满足由第一数据值构建的第一预设条件，或满足由第二数据值构建的第二预设条件，若存在任意一个预设条件成立，则跳转至S36步骤，若否，则进入下一步骤；
S34：根据轨迹观测值
Figure SMS_79
相对第一直线/>
Figure SMS_80
与第二直线/>
Figure SMS_81
的位置更新第一直线与第二直线的上下容忍点、斜率以及上凸壳与下凸壳；
S35：更新时间序号：
Figure SMS_82
，并返回S32步骤；
S36：输出时间序号
Figure SMS_83
，通过更新后的第一直线与第二直线的交点与斜率生成第/>
Figure SMS_84
维数据序列对应的可行分割直线。
进一步地，所述S32步骤中：
所述第一数据值的获取公式为：
Figure SMS_85
；
所述第二数据值的获取公式为：
Figure SMS_86
；
式中，
Figure SMS_87
表示第一直线上容忍点/>
Figure SMS_88
对应的时间点，/>
Figure SMS_89
表示第二直线下容忍点
Figure SMS_90
对应的时间点，/>
Figure SMS_91
表示第一数据值，/>
Figure SMS_92
表示第二数据值。
进一步地，所述S33步骤中：
所述第一预设条件的表达式为：
Figure SMS_93
；
所述第二预设条件的表达式为：
Figure SMS_94
。
进一步地，所述S34步骤，具体包括：
判断
Figure SMS_96
是否位于第二直线/>
Figure SMS_102
的上方，若是，则从上凸壳/>
Figure SMS_106
中获取点/>
Figure SMS_98
，使得/>
Figure SMS_101
最小，更新第二直线的下容忍点/>
Figure SMS_105
、第二直线的上容忍点
Figure SMS_108
，并删除上凸壳/>
Figure SMS_95
中位于点/>
Figure SMS_99
之前的点，更新/>
Figure SMS_103
，将/>
Figure SMS_107
加至下凸壳/>
Figure SMS_97
的末端并采用三角检验法更新下凸壳/>
Figure SMS_100
；其中，/>
Figure SMS_104
。
进一步地，所述S34步骤，还包括：
判断
Figure SMS_111
是否位于第一直线/>
Figure SMS_116
的下方，若是，则从下凸壳/>
Figure SMS_120
中获取点
Figure SMS_110
，使得/>
Figure SMS_114
最大，更新第一直线的上容忍点/>
Figure SMS_118
、第一直线的下容忍点
Figure SMS_122
，并删除下凸壳/>
Figure SMS_112
中位于点/>
Figure SMS_113
之前的点，更新/>
Figure SMS_117
，将/>
Figure SMS_121
加至上凸壳/>
Figure SMS_109
的末端并采用三角检验法更新上凸壳/>
Figure SMS_115
；其中，/>
Figure SMS_119
。
进一步地，所述S36步骤中，可行分割直线的生成具体包括：
获取更新后的第一直线与第二直线的交点；
通过更新后的
Figure SMS_123
与/>
Figure SMS_124
获取目标斜率/>
Figure SMS_125
；
通过目标斜率以及更新后第一直线与第二直线的交点生成可行分割直线
Figure SMS_126
；所述目标斜率的获取公式为：/>
Figure SMS_127
。
与现有技术相比，本发明至少含有以下有益效果：
本发明基于分割初始时刻
Figure SMS_130
与最大误差/>
Figure SMS_134
利用DConnPLA算法分别分割各维度数据序列，得到各维度数据序列对应的可行分割直线，实现了对可行直线空间的完全搜索，另外，本发明通过获取各可行分割直线/>
Figure SMS_138
中最短可行分割直线对应的分割结束时刻/>
Figure SMS_131
，通过可行分割直线的表达式/>
Figure SMS_135
计算各维度数据序列在分割初始时刻/>
Figure SMS_139
对应的初始值
Figure SMS_142
和在分割结束时刻/>
Figure SMS_128
对应的结束值/>
Figure SMS_133
，通过各维度数据序列对应的初始值
Figure SMS_137
和结束值/>
Figure SMS_141
得到目标起点和目标终点，通过目标起点和目标终点分割当前的/>
Figure SMS_129
维实时轨迹数据/>
Figure SMS_132
，并以分割结束时刻/>
Figure SMS_136
加1作为下一段/>
Figure SMS_140
维实时轨迹数据的分割初始时刻，其实现了三维及三维以上空间数据的最优分割，同时实现了对轨迹数据的不间断压缩即在线压缩。
附图说明
图1为一种用于轨迹数据的实时压缩方法流程图；
图2为维度数据序列对应的直线与凸壳示意图；
图3为FSW算法极限线构造示意图；
图4为CISED算法圆锥构造示意图。
具体实施方式
以下是本发明的具体实施例并结合附图，对本发明的技术方案作进一步的描述，但本发明并不限于这些实施例。
实施例一
为了实现对可行直线空间的完全搜索，以实现对三维及三维以上空间的最优分割，如图1所示，本发明提出了一种用于轨迹数据的实时压缩方法，包括步骤：
S1：获取
Figure SMS_145
维实时轨迹数据/>
Figure SMS_148
并设定最大误差/>
Figure SMS_153
，其中，
Figure SMS_146
，/>
Figure SMS_156
表示数据点/>
Figure SMS_155
的时间戳索引，/>
Figure SMS_159
为常数，/>
Figure SMS_143
表示在/>
Figure SMS_157
时刻的/>
Figure SMS_149
个轨迹观测值；获取/>
Figure SMS_152
维实时轨迹数据对应的维度数据序列
Figure SMS_144
，其中，/>
Figure SMS_151
表示第/>
Figure SMS_150
维的数据序列，/>
Figure SMS_158
为/>
Figure SMS_147
维/>
Figure SMS_154
时刻的轨迹观测值；
需要解释的是，在地理信息系统中，轨迹数据不仅包含经度和纬度信息，还额外记录了某些应用领域的关键数值指标，譬如速度、温度或湿度等等，从而轨迹数据中包括了多个维度的数据（即多维空间数据），而现有技术中，是通过空间欧式距离（也称同步误差）进行的度量，只能在二维空间进行划分无法对三维及三维以上的空间完成算法处理，因为这些算法所采用的垂直误差或同步误差均不能完成对可行直线空间的完全搜索，因而无法获得最优化的分割结果，对此，本发明为了实现三维及三维以上空间数据的最优分割，采用分治-融合的思想，对每一个维度的数据序列单独实施一维最优化
Figure SMS_160
压缩算法即DConnPLA算法（分治），通过DConnPLA算法确保每个分割能拟合最多的数据点，然后通过各维度数据序列对应的可行分割直线获取目标起点与目标终点（融合），通过目标起点与目标终点完成分割，从而实现了对可行直线空间的完全搜索，完成了三维及三维以上空间数据的最优分割。
特别需要说明的是，DConnPLA算法在线性时间复杂度下可以在线地处理一维时序数据，并获得最少的不连续分割个数，为了最小化生成的分割数，DConnPLA算法确保每个分割能拟合最多的数据点。这是通过在添加新点的过程中，增量式地调整所有表示线（目标起点和目标终点形成的直线为分割维度数据序列的一个可行直线，简称表示线）的范围来完成的。当在最大误差范围内没有可行直线可以近似新的点时，构造过程停止，新的分割开始。在生成分割的过程中，随着新数据的引入，两条极值线（第一直线与第二直线）被迭代更新。为了有效地更新极值线，需要增量地保持两个凸壳结构，以避免大量的搜索代价。
S2：初始化分割序号
Figure SMS_161
与分割初始时刻/>
Figure SMS_162
；
所述S2步骤中，分割序号
Figure SMS_163
与分割初始时刻/>
Figure SMS_164
的初始值均等于1；
S3：基于分割初始时刻
Figure SMS_165
与最大误差/>
Figure SMS_166
利用DConnPLA算法分别分割各维度数据序列，得到各维度数据序列对应的可行分割直线，所述可行分割直线的表达式为/>
Figure SMS_167
，其中
Figure SMS_168
表示时刻；
需要说明的是，本实施例中可利用DConnPLA算法串行或并行分割各维度数据序列。本发明既保持了DConnPLA算法的线性时间复杂度和在线处理数据的能力，又实现了对可行直线空间的完全搜索。
所述S3步骤中，利用DConnPLA算法分割维度数据序列，具体包括：
S31：初始化时间序号
Figure SMS_169
，初始化维度数据序列对应的上凸壳、下凸壳以及第一直线与第二直线的斜率：
所述时间序号
Figure SMS_170
的初始化公式为：/>
Figure SMS_171
；
所述上凸壳的表达式为：
Figure SMS_172
；
所述下凸壳的表达式为：
Figure SMS_173
；
所述第一直线的表达式为：
Figure SMS_174
；
所述第二直线的表达式为：
Figure SMS_175
；
所述第一直线对应斜率的获取公式为：
Figure SMS_176
；
所述第二直线对应斜率的获取公式为：
Figure SMS_177
；其中：
Figure SMS_178
，/>
Figure SMS_179
；
Figure SMS_180
，/>
Figure SMS_181
；
式中，
Figure SMS_188
表示分割初始时刻/>
Figure SMS_185
对应的时间点，/>
Figure SMS_193
表示第/>
Figure SMS_183
维数据序列中/>
Figure SMS_190
时间点对应的轨迹观测值，/>
Figure SMS_191
表示最大误差，/>
Figure SMS_198
、/>
Figure SMS_189
分别表示第一直线的上容忍点与下容忍点，/>
Figure SMS_197
、/>
Figure SMS_182
分别表示第二直线的上容忍点与下容忍点，/>
Figure SMS_196
表示第一直线的斜率，/>
Figure SMS_184
表示第二直线的斜率，/>
Figure SMS_195
表示上凸壳，/>
Figure SMS_186
表示下凸壳；/>
Figure SMS_192
表示上一段第/>
Figure SMS_187
维数据序列分割时输出的时间序号/>
Figure SMS_194
；
需要解释的是，第一次分割时，
Figure SMS_199
，/>
Figure SMS_200
，加3的原因在于在每一次的分割过程中，会使用前两个时间点/>
Figure SMS_201
、/>
Figure SMS_202
对应的轨迹观测值初始化上凸壳、下凸壳以及第一直线与第二直线，因此在S32步骤中是从第三个时间点对应的轨迹观测值开始处理。
S32：通过维度数据序列中轨迹观测值
Figure SMS_203
对应的时间点/>
Figure SMS_204
与第一直线获取第一数据值，通过轨迹观测值/>
Figure SMS_205
对应的时间点/>
Figure SMS_206
与第二直线获取第二数据值；
所述S32步骤中：
所述第一数据值的获取公式为：
Figure SMS_207
；
所述第二数据值的获取公式为：
Figure SMS_208
；
式中，
Figure SMS_209
表示第一直线上容忍点/>
Figure SMS_210
对应的时间点，/>
Figure SMS_211
表示第二直线下容忍点
Figure SMS_212
对应的时间点，/>
Figure SMS_213
表示第一数据值，/>
Figure SMS_214
表示第二数据值。
S33：判断轨迹观测值
Figure SMS_215
是否满足由第一数据值构建的第一预设条件，或满足由第二数据值构建的第二预设条件，若存在任意一个预设条件成立，则跳转至S36步骤，若否，则进入下一步骤；
所述S33步骤中：
所述第一预设条件的表达式为：
Figure SMS_216
；
所述第二预设条件的表达式为：
Figure SMS_217
。
需要说明的是，当轨迹观测值
Figure SMS_218
满足第一预设条件或满足第二预设条件时表明轨迹观测值/>
Figure SMS_219
对应的点在第一直线与第二直线的外侧也即处于最大误差范围以外，所述最大误差范围为：/>
Figure SMS_220
。
S34：根据轨迹观测值
Figure SMS_221
相对第一直线/>
Figure SMS_222
与第二直线/>
Figure SMS_223
的位置更新第一直线与第二直线的上下容忍点、斜率以及上凸壳与下凸壳；
所述S34步骤，具体包括：
判断
Figure SMS_227
是否位于第二直线/>
Figure SMS_228
的上方，若是，则从上凸壳/>
Figure SMS_232
中获取点/>
Figure SMS_224
，使得/>
Figure SMS_230
最小，更新第二直线的下容忍点/>
Figure SMS_234
、第二直线的上容忍点
Figure SMS_236
，并删除上凸壳/>
Figure SMS_225
中位于点/>
Figure SMS_231
之前的点，更新/>
Figure SMS_235
，将/>
Figure SMS_237
加至下凸壳/>
Figure SMS_226
的末端并采用三角检验法更新下凸壳/>
Figure SMS_229
；其中，/>
Figure SMS_233
。
所述S34步骤，还包括：
判断
Figure SMS_238
是否位于第一直线/>
Figure SMS_243
的下方，若是，则从下凸壳/>
Figure SMS_247
中获取点
Figure SMS_241
，使得/>
Figure SMS_244
最大，更新第一直线的上容忍点/>
Figure SMS_248
、第一直线的下容忍点
Figure SMS_250
，并删除下凸壳/>
Figure SMS_240
中位于点/>
Figure SMS_245
之前的点，更新/>
Figure SMS_249
，将/>
Figure SMS_251
加至上凸壳/>
Figure SMS_239
的末端并采用三角检验法更新上凸壳/>
Figure SMS_242
；其中，/>
Figure SMS_246
。
直线与凸壳之间的位置关系可见图2所示。
需要解释的是，采用三角检验法更新下凸壳
Figure SMS_252
，具体为通过迭代地检查最近连续的三个点来更新下凸壳/>
Figure SMS_253
，其更新规则为：如果中间点在其他两点形成的线的下方或线上，则删除中间点，迭代这个过程，直到中间点不再被移除（中间点在线的上方则不再被移除）。
上凸壳
Figure SMS_254
的更新规则为：如果中间点在其他两点形成的线的上方或线上，则删除中间点，迭代这个过程，直到中间点不再被移除（中间点在线的下方则不再被移除）。
S35：更新时间序号：
Figure SMS_255
，并返回S32步骤；
S36：输出时间序号
Figure SMS_256
，通过更新后的第一直线与第二直线的交点与斜率生成第/>
Figure SMS_257
维数据序列对应的可行分割直线。
所述S36步骤中，可行分割直线的生成具体包括：
获取更新后的第一直线与第二直线的交点；
通过更新后的
Figure SMS_258
与/>
Figure SMS_259
获取目标斜率/>
Figure SMS_260
；
通过目标斜率以及更新后第一直线与第二直线的交点生成可行分割直线
Figure SMS_261
；所述目标斜率的获取公式为：/>
Figure SMS_262
。
S4：获取各可行分割直线
Figure SMS_264
中最短可行分割直线对应的分割结束时刻/>
Figure SMS_269
，通过可行分割直线的表达式/>
Figure SMS_272
计算各维度数据序列在分割初始时刻/>
Figure SMS_265
对应的初始值/>
Figure SMS_268
和在分割结束时刻/>
Figure SMS_271
对应的结束值/>
Figure SMS_274
，通过各维度数据序列对应的初始值/>
Figure SMS_263
和结束值/>
Figure SMS_267
得到目标起点/>
Figure SMS_270
和目标终点/>
Figure SMS_273
；所述目标起点与目标终点为在/>
Figure SMS_266
维空间中的数据点；
S5：通过目标起点
Figure SMS_275
和目标终点/>
Figure SMS_276
分割当前的/>
Figure SMS_277
维实时轨迹数据/>
Figure SMS_278
，设定
Figure SMS_279
、/>
Figure SMS_280
，并返回S3步骤继续下一段/>
Figure SMS_281
维实时轨迹数据的分割。
为了证明本发明提出的实时压缩方法可以得到最优的划分结果即本发明最终得到的
Figure SMS_282
值为最优值，本实施例推理如下：
定理：对于多维实时轨迹数据
Figure SMS_283
和给定的最大误差/>
Figure SMS_284
，本发明提出的实时压缩方法可以获得最优的分割结果。
证明：在给定最大误差
Figure SMS_285
下，记本发明获取/>
Figure SMS_286
个分割的时间索引集合为
Figure SMS_287
，显然，/>
Figure SMS_288
时，该定理是成立的；
如果
Figure SMS_293
，假设/>
Figure SMS_290
不是最优的，那么存在/>
Figure SMS_299
Figure SMS_294
个分割覆盖数据/>
Figure SMS_300
。这意味着这/>
Figure SMS_292
个分割中存在某个分割(假设该分割为第/>
Figure SMS_297
个)满足/>
Figure SMS_305
，其中，/>
Figure SMS_307
、/>
Figure SMS_289
分别表示/>
Figure SMS_302
个分割中第/>
Figure SMS_296
个分割对应的分割初始时刻与分割结束时刻，/>
Figure SMS_298
、/>
Figure SMS_304
分别表示/>
Figure SMS_306
个分割中第/>
Figure SMS_291
个分割对应的分割初始时刻与分割结束时刻；故在各个维度数据序列上映射后得到的分割同样满足该特性。由于这些映射分割中，至少存在一个由DConnPLA算法生成的完整分割（各可行分割直线/>
Figure SMS_301
中最短的可行分割直线表示为一次完整的分割），该分割被保证可以近似最多的数据点（即一维的DConnPLA算法能保证可以近似最多的数据点，使分割达到最长，分割初始时刻与分割结束时刻之间的时间差达到最大），故与
Figure SMS_295
矛盾，假设不成立，/>
Figure SMS_303
是最优的。
因而本发明提出的实时压缩方法可以得到最优的划分结果，其最终得到的
Figure SMS_308
值为最优值。
本发明基于分割初始时刻
Figure SMS_311
与最大误差/>
Figure SMS_313
利用DConnPLA算法分别分割各维度数据序列，得到各维度数据序列对应的可行分割直线，实现了对可行直线空间的完全搜索，另外，本发明通过获取各可行分割直线/>
Figure SMS_317
中最短可行分割直线对应的分割结束时刻/>
Figure SMS_312
，通过可行分割直线的表达式/>
Figure SMS_315
计算各维度数据序列在分割初始时刻/>
Figure SMS_319
对应的初始值
Figure SMS_322
和在分割结束时刻/>
Figure SMS_309
对应的结束值/>
Figure SMS_314
，通过各维度数据序列对应的初始值
Figure SMS_318
和结束值/>
Figure SMS_321
得到目标起点和目标终点，通过目标起点和目标终点分割当前的/>
Figure SMS_310
维实时轨迹数据/>
Figure SMS_316
，并以分割结束时刻/>
Figure SMS_320
加1作为下一段/>
Figure SMS_323
维实时轨迹数据的分割初始时刻，其实现了三维及三维以上空间数据的最优分割，同时实现了对轨迹数据的不间断压缩即在线压缩。
实施例二
为了说明本发明提出的实时压缩方法的技术优势，将其实例化即设定
Figure SMS_324
，并与CISED算法在存储、时间消耗和内存占用方面进行类比。本实施例将本发明提出的实时压缩方法定义为MDisPLA算法，以便后续的描述。
需要说明的是，在现存算法中，CISED算法是用于轨迹压缩领域上非常优秀的算法之一，它能在线性时间复杂度内，以常数级的内存占用，生成连续直线表示的分割。CISED算法的思想源于FSW算法，它将FSW算法处理一维时序数据的能力推广至二维。FSW算法通过维护两条极限线来构造可行直线表示的区域，当该动态更新的区域为空时，一段分割产生并开始下一段分割。如图3所示，在最大误差
Figure SMS_332
限定下，FSW算法在t-y坐标系内通过固定的起始点/>
Figure SMS_339
以及给定数据点/>
Figure SMS_347
构造极限线/>
Figure SMS_328
和/>
Figure SMS_334
，用于划定可同时近似/>
Figure SMS_341
和/>
Figure SMS_349
的可行直线表示区域；同样，在确定/>
Figure SMS_327
和/>
Figure SMS_333
的极限线/>
Figure SMS_330
和/>
Figure SMS_340
后，划定可同时近似/>
Figure SMS_326
和/>
Figure SMS_337
的可行直线表示区域，通过判定两个区域存在公共区域/>
Figure SMS_329
和/>
Figure SMS_338
可表明存在直线表示可在最大误差下近似/>
Figure SMS_352
、/>
Figure SMS_357
和/>
Figure SMS_355
三个数据点。而CISED算法将FSW算法进行推广，用于压缩二维轨迹数据，如图4所示，在CISED中，最大误差被推广为同步误差（三维空间
Figure SMS_358
中，球体的半径），数据点/>
Figure SMS_325
和/>
Figure SMS_335
的可行直线表示区域推广为圆锥/>
Figure SMS_343
，数据点/>
Figure SMS_351
和/>
Figure SMS_348
的可行直线表示区域推广为圆锥/>
Figure SMS_356
；通过判断/>
Figure SMS_342
和/>
Figure SMS_350
有交集可表明数据点/>
Figure SMS_345
、/>
Figure SMS_353
和/>
Figure SMS_346
可在同步误差/>
Figure SMS_354
下用一条直线表示近似。其中，/>
Figure SMS_336
表示时间序号，
Figure SMS_344
、/>
Figure SMS_331
为自然数，用于表示时间索引。
尽管MDisPLA算法和CISED算法都可以用于对轨迹数据的压缩，但是MDisPLA和CISED算法存在三点不同：
第一点：处理问题不同，MDisPLA算法是在限定各个维度数据序列最大误差的条件下，对多维数据进行压缩，最终生成不连续的直线表示；而CISED算法只针对二维数据，通过限定三维空间中数据的同步误差，用于生成连续的直线表示。CISED算法向多维数据压缩问题上的推广不能显然地得到；
第二点：压缩存储结果不同，MDisPLA算法可生成最小个数的不连续分割；而CISED并不是最优化的算法，不能生成最少的连续分割。一般来说，相同分割数的限定下，不连续表示要比连续表示的压缩存储大。但是，由于CISED算法进行了初始点的固定，以及同步误差的限定，因此MDisPLA算法比CISED算法生成的分割数要小的多；
第三点：时间和内存消耗不同，在时间消耗方面，尽管MDisPLA算法和CISED算法都具有线性时间复杂度，但是二者因为计算机制的不同，在时间消耗方面存在明显的差异性。MDisPLA算法的主要计算消耗在于对各个维度数据序列进行压缩过程中凸壳的更新；而CISED算法的主要计算消耗在对圆锥可行区域的交集判断以及更新方面。在内存消耗方面，MDisPLA算法在压缩过程中，仍然需要缓存当前分割的凸壳，因此和DConnPLA算法保持一致；而CISED算法的临时内存占用，依赖于在可行区域的交集确定过程中逼近圆的正多边形边的个数，它也是算法的输入参数之一。
本实施例挑选了两个具有代表性的GPS数据集—Roma数据集和Mopsi数据集来展现测试结果。Roma数据集记录了320辆出租车在罗马市区一个月的行驶数据，本实施例筛选出编号为146的出租车数据，记为Roma146；Mopsi数据集记录了51位参与者在芬兰约恩苏地区的活动轨迹数据，本实施例挑选了数据规模最大的7号子集作为测试对象之一，记为Mopsi7。测试平台为一台笔记本图形工作站，具有2.60GHz的CPU和40GB的内存。为了公平性，本实施例在Eclipse平台下采用Java语言对本发明提出的MDisPLA算法和CISED算法进行统一实现。
对于
Figure SMS_359
个不连续分割的存储问题，本实施例采用如下存储策略：将整型或者浮点型数据看作一个存储单元，那么/>
Figure SMS_360
个分割包含/>
Figure SMS_361
个/>
Figure SMS_362
维的数据点（一次分割保存一个目标起点和一个目标终点），需要/>
Figure SMS_363
个存储单元。又由于当前分割的分割结束时刻/>
Figure SMS_364
和下一个分割的分割初始时刻/>
Figure SMS_365
存在如下关系：
Figure SMS_366
因此，相邻两个不连续的分割中，前一个分割的分割结束时刻
Figure SMS_367
和后一个分割的分割初始时刻/>
Figure SMS_368
可共用时间索引，每增加一个新的分割可减少1个存储单元。故本实施例中存储/>
Figure SMS_369
个分割总共只需要/>
Figure SMS_370
个存储单元。
在算法参数设置方面，CISED算法中同步误差的单位为米，而MDisPLA中最大误差的单位为度，所以依据
Figure SMS_371
进行近似转换；CISED中正多边形的边数设定为16。
本实施例以同步误差60米和最大误差0.0006846度下MDisPLA和CISED的输出结果为例进行说明，该结果表明在分割数、存储、时间和最大内存方面，MDisPLA都要优于CISED，输出结果见下表1所示：
Figure SMS_372
。
通过表1中的实验数据得出：
分割数和存储方面：在Roma数据集和Mopsi数据集上，MDisPLA比CISED的分割数分别减少了约50%和42%，而在存储上分别减少了17%和4%左右。这些差异源于两个方面的原因：(1)CISED中同步误差的限定、固定初始点以及多边形的边数设定都会影响算法对可行直线表示的搜索能力，所以在同样分割数下，MDisPLA比CISED近似的点更多；(2)一般来说，相同分割数下，不连续直线表示的存储要比连续直线表示的存储多。上述事实以及MDisPLA的存储策略共同造成了两个算法从分割数差异到存储差异的变化。
时间消耗和最大内存方面：时间消耗上，MDisPLA比CISED在Roma146上快约8倍，而在Mopsi7上快20倍左右；在最大内存占用上，MDisPLA在两个数据集上分别为0.1403kb和0.2456kb，而CISED均为0.3125kb。上述实验结果的原因在于：(1)MDisPLA的主要计算在凸壳的更新，需要临时存储凸壳，由于轨迹数据中经纬度数据局部片段差异较小以及最大误差设定较小，造成临时存储凸壳的规模较小，故计算消耗和内存占用均偏低；(2)CISED中，需要在三维空间中判断圆锥是否存在公共区域，而在这过程中圆形采用正多边形逼近的策略，故计算消耗和内存占用依赖于正多边形边数，故计算消耗和内存占用相对MDisPLA要高。
本实施例通过与优秀的轨迹压缩算法CISED进行比较，验证了MDisPLA算法是稳健的，其整体性能要优于CISED算法。具体地，在两个数据集上，MDisPLA以更低的内存消耗降低了约46%的分割个数和10.5%的存储个数，并提高了14倍左右的处理速度。
需要说明，本发明实施例中所有方向性指示（诸如上、下、左、右、前、后……）仅用于解释在某一特定姿态下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。
另外，在本发明中如涉及“第一”、“第二”、“一”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
在本发明中，除非另有明确的规定和限定，术语“连接”、“固定”等应做广义理解，例如，“固定”可以是固定连接，也可以是可拆卸连接，或成一体；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系，除非另有明确的限定。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
另外，本发明各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。

Claims (7)

1.一种用于轨迹数据的实时压缩方法，其特征在于，包括步骤：
S1：获取
Figure QLYQS_7
维实时轨迹数据/>
Figure QLYQS_3
并设定最大误差/>
Figure QLYQS_15
，其中，
Figure QLYQS_1
，/>
Figure QLYQS_12
表示数据点/>
Figure QLYQS_4
的时间戳索引，/>
Figure QLYQS_16
为常数，/>
Figure QLYQS_5
表示在/>
Figure QLYQS_13
时刻的/>
Figure QLYQS_8
个轨迹观测值；获取/>
Figure QLYQS_11
维实时轨迹数据对应的维度数据序列
Figure QLYQS_6
，其中，/>
Figure QLYQS_9
表示第/>
Figure QLYQS_10
维的数据序列，/>
Figure QLYQS_17
为/>
Figure QLYQS_2
维/>
Figure QLYQS_14
时刻的轨迹观测值；
S2：初始化分割序号
Figure QLYQS_18
与分割初始时刻/>
Figure QLYQS_19
；
S3：基于分割初始时刻
Figure QLYQS_20
与最大误差/>
Figure QLYQS_21
利用DConnPLA算法分别分割各维度数据序列，得到各维度数据序列对应的可行分割直线，所述可行分割直线的表达式为/>
Figure QLYQS_22
，其中/>
Figure QLYQS_23
表示时刻；
S4：获取各可行分割直线
Figure QLYQS_26
中最短可行分割直线对应的分割结束时刻/>
Figure QLYQS_28
，通过可行分割直线的表达式/>
Figure QLYQS_31
计算各维度数据序列在分割初始时刻/>
Figure QLYQS_25
对应的初始值/>
Figure QLYQS_29
和在分割结束时刻/>
Figure QLYQS_32
对应的结束值/>
Figure QLYQS_34
，通过各维度数据序列对应的初始值/>
Figure QLYQS_24
和结束值/>
Figure QLYQS_27
得到目标起点/>
Figure QLYQS_30
和目标终点/>
Figure QLYQS_33
；
S5：通过目标起点
Figure QLYQS_35
和目标终点/>
Figure QLYQS_36
分割当前的/>
Figure QLYQS_37
维实时轨迹数据/>
Figure QLYQS_38
，设定/>
Figure QLYQS_39
、
Figure QLYQS_40
，并返回S3步骤继续下一段/>
Figure QLYQS_41
维实时轨迹数据的分割。
2.根据权利要求1所述的一种用于轨迹数据的实时压缩方法，其特征在于，所述S2步骤中，分割序号
Figure QLYQS_42
与分割初始时刻/>
Figure QLYQS_43
的初始值均等于1；
所述S3步骤中，利用DConnPLA算法分割维度数据序列，具体包括：
S31：初始化时间序号
Figure QLYQS_44
，初始化维度数据序列对应的上凸壳、下凸壳以及第一直线与第二直线的斜率：
所述时间序号的初始化公式为：
Figure QLYQS_45
；
所述上凸壳的表达式为：
Figure QLYQS_46
；
所述下凸壳的表达式为：
Figure QLYQS_47
；
所述第一直线的表达式为：
Figure QLYQS_48
；
所述第二直线的表达式为：
Figure QLYQS_49
；
所述第一直线对应斜率的获取公式为：
Figure QLYQS_50
；
所述第二直线对应斜率的获取公式为：
Figure QLYQS_51
；其中：
Figure QLYQS_52
，/>
Figure QLYQS_53
；
Figure QLYQS_54
，/>
Figure QLYQS_55
；
式中，
Figure QLYQS_62
表示分割初始时刻/>
Figure QLYQS_57
对应的时间点，/>
Figure QLYQS_67
表示第/>
Figure QLYQS_58
维数据序列中/>
Figure QLYQS_71
时间点对应的轨迹观测值，/>
Figure QLYQS_65
表示最大误差，/>
Figure QLYQS_72
、/>
Figure QLYQS_59
分别表示第一直线的上容忍点与下容忍点，/>
Figure QLYQS_69
、/>
Figure QLYQS_56
分别表示第二直线的上容忍点与下容忍点，/>
Figure QLYQS_68
表示第一直线的斜率，/>
Figure QLYQS_61
表示第二直线的斜率，/>
Figure QLYQS_64
表示上凸壳，/>
Figure QLYQS_63
表示下凸壳；/>
Figure QLYQS_70
表示上一段第/>
Figure QLYQS_60
维数据序列分割时输出的时间序号/>
Figure QLYQS_66
；
S32：通过维度数据序列中轨迹观测值
Figure QLYQS_73
对应的时间点/>
Figure QLYQS_74
与第一直线获取第一数据值，通过轨迹观测值/>
Figure QLYQS_75
对应的时间点/>
Figure QLYQS_76
与第二直线获取第二数据值；
S33：判断轨迹观测值
Figure QLYQS_77
是否满足由第一数据值构建的第一预设条件，或满足由第二数据值构建的第二预设条件，若存在任意一个预设条件成立，则跳转至S36步骤，若否，则进入下一步骤；
S34：根据轨迹观测值
Figure QLYQS_78
相对第一直线/>
Figure QLYQS_79
与第二直线/>
Figure QLYQS_80
的位置更新第一直线与第二直线的上下容忍点、斜率以及上凸壳与下凸壳；
S35：更新时间序号
Figure QLYQS_81
：，并返回S32步骤；
S36：输出时间序号
Figure QLYQS_82
，通过更新后的第一直线与第二直线的交点与斜率生成第/>
Figure QLYQS_83
维数据序列对应的可行分割直线。
3.根据权利要求2所述的一种用于轨迹数据的实时压缩方法，其特征在于，所述S32步骤中：
所述第一数据值的获取公式为：
Figure QLYQS_84
；
所述第二数据值的获取公式为：
Figure QLYQS_85
；
式中，
Figure QLYQS_86
表示第一直线上容忍点/>
Figure QLYQS_87
对应的时间点，/>
Figure QLYQS_88
表示第二直线下容忍点/>
Figure QLYQS_89
对应的时间点，/>
Figure QLYQS_90
表示第一数据值，/>
Figure QLYQS_91
表示第二数据值。
4.根据权利要求3所述的一种用于轨迹数据的实时压缩方法，其特征在于，所述S33步骤中：
所述第一预设条件的表达式为：
Figure QLYQS_92
；
所述第二预设条件的表达式为：
Figure QLYQS_93
。
5.根据权利要求4所述的一种用于轨迹数据的实时压缩方法，其特征在于，所述S34步骤，具体包括：
判断
Figure QLYQS_95
是否位于第二直线/>
Figure QLYQS_100
的上方，若是，则从上凸壳/>
Figure QLYQS_104
中获取点/>
Figure QLYQS_96
，使得/>
Figure QLYQS_101
最小，更新第二直线的下容忍点/>
Figure QLYQS_105
、第二直线的上容忍点
Figure QLYQS_107
，并删除上凸壳/>
Figure QLYQS_94
中位于点/>
Figure QLYQS_98
之前的点，更新/>
Figure QLYQS_102
，将/>
Figure QLYQS_106
加至下凸壳/>
Figure QLYQS_97
的末端并采用三角检验法更新下凸壳/>
Figure QLYQS_99
；其中，/>
Figure QLYQS_103
。
6.根据权利要求5所述的一种用于轨迹数据的实时压缩方法，其特征在于，所述S34步骤，还包括：
判断
Figure QLYQS_109
是否位于第一直线/>
Figure QLYQS_115
的下方，若是，则从下凸壳/>
Figure QLYQS_119
中获取点/>
Figure QLYQS_111
，使得/>
Figure QLYQS_112
最大，更新第一直线的上容忍点/>
Figure QLYQS_116
、第一直线的下容忍点/>
Figure QLYQS_120
，并删除下凸壳/>
Figure QLYQS_108
中位于点/>
Figure QLYQS_114
之前的点，更新/>
Figure QLYQS_118
，将/>
Figure QLYQS_121
加至上凸壳
Figure QLYQS_110
的末端并采用三角检验法更新上凸壳/>
Figure QLYQS_113
；其中，/>
Figure QLYQS_117
。
7.根据权利要求6所述的一种用于轨迹数据的实时压缩方法，其特征在于，所述S36步骤中，可行分割直线的生成具体包括：
获取更新后的第一直线与第二直线的交点；
通过更新后的
Figure QLYQS_122
与/>
Figure QLYQS_123
获取目标斜率/>
Figure QLYQS_124
；
通过目标斜率以及更新后第一直线与第二直线的交点生成可行分割直线
Figure QLYQS_125
；所述目标斜率的获取公式为：/>
Figure QLYQS_126
。
CN202310572745.4A 2023-05-22 2023-05-22 一种用于轨迹数据的实时压缩方法 Active CN116303479B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310572745.4A CN116303479B (zh) 2023-05-22 2023-05-22 一种用于轨迹数据的实时压缩方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310572745.4A CN116303479B (zh) 2023-05-22 2023-05-22 一种用于轨迹数据的实时压缩方法

Publications (2)

Publication Number Publication Date
CN116303479A true CN116303479A (zh) 2023-06-23
CN116303479B CN116303479B (zh) 2023-08-22

Family

ID=86789177

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310572745.4A Active CN116303479B (zh) 2023-05-22 2023-05-22 一种用于轨迹数据的实时压缩方法

Country Status (1)

Country Link
CN (1) CN116303479B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2410395A2 (en) * 2010-07-21 2012-01-25 Palo Alto Research Center Incorporated System and method for real-time mapping of an indoor environment using mobile robots with limited sensing
CN109443363A (zh) * 2018-11-09 2019-03-08 厦门大学 分治确定性路径寻优算法
CN110570501A (zh) * 2018-06-06 2019-12-13 腾讯科技（深圳）有限公司 一种线条动画绘制方法及其设备、存储介质、电子设备
CN111299815A (zh) * 2020-02-13 2020-06-19 西安交通大学 一种用于低灰度橡胶垫视觉检测与激光切割轨迹规划方法
US20230028461A1 (en) * 2021-07-15 2023-01-26 Robert Bosch Gmbh Identification of connection patterns on the basis of trajectory data

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2410395A2 (en) * 2010-07-21 2012-01-25 Palo Alto Research Center Incorporated System and method for real-time mapping of an indoor environment using mobile robots with limited sensing
CN110570501A (zh) * 2018-06-06 2019-12-13 腾讯科技（深圳）有限公司 一种线条动画绘制方法及其设备、存储介质、电子设备
CN109443363A (zh) * 2018-11-09 2019-03-08 厦门大学 分治确定性路径寻优算法
CN111299815A (zh) * 2020-02-13 2020-06-19 西安交通大学 一种用于低灰度橡胶垫视觉检测与激光切割轨迹规划方法
US20230028461A1 (en) * 2021-07-15 2023-01-26 Robert Bosch Gmbh Identification of connection patterns on the basis of trajectory data

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YASIR LATIF等: ""Go straight, turn right: Pose graph reduction through trajectory segmentation using line segments"", 《2013 EUROPEAN CONFERENCE ON MOBILE ROBOTS》, pages 1 - 4 *
张延玲等: ""移动对象子轨迹段分割与聚类算法"", 《计算机工程与应用》, pages 65 - 68 *

Also Published As

Publication number Publication date
CN116303479B (zh) 2023-08-22

Similar Documents

Publication Publication Date Title
US9256693B2 (en) Recommendation system with metric transformation
US11188738B2 (en) System and method associated with progressive spatial analysis of prodigious 3D data including complex structures
US8391610B2 (en) Rapid, spatial-data viewing and manipulating including data partition and indexing
US10593110B2 (en) Method and device for computing a path in a game scene
CN112417965B (zh) 激光点云处理方法、电子装置和存储介质
US9286312B2 (en) Data coreset compression
Lehoczky et al. Present Software of photogrammetric processing of digital images
Gong et al. An efficient point cloud management method based on a 3D R-tree
CN113190515A (zh) 基于异构并行计算的城市级海量点云坐标转换方法
US20200349759A1 (en) Pyramid generation via depth-first walk
JP2015132539A (ja) 日照量計算装置、経路提案装置および日照量計算方法
WO2021103512A1 (en) Method and apparatus for generating electronic map
CN112817545B (zh) 在线即分析影像与格网立方体数据存储管理方法及系统
US10019456B2 (en) Recovering free space in nonvolatile storage with a computer storage system supporting shared objects
CN109191484B (zh) 一种从机载激光雷达点云中快速提取平面片的方法
CN109255004A (zh) 面向地理空间数据的精度保证的在线交互式可视化方法
CN116303479B (zh) 一种用于轨迹数据的实时压缩方法
CN116797747A (zh) 水下检测数据可视化方法、装置、计算机设备和存储介质
CN111190893A (zh) 建立特征索引的方法和装置
CN109947884B (zh) 一种全球ICESat/GLAS点云的高效组织查询方法
CN109597865B (zh) 基于嵌入式平台的海量地理信息数据的存储与检索方法
CN113822768A (zh) 社区网络的处理方法、装置、设备及存储介质
CN106504325B (zh) 一种基于cuda的dem特征点提取并行化方法
CN111797880A (zh) 数据处理方法、装置、存储介质及电子设备
CN113496550B (zh) 一种dsm的计算方法、装置、计算机设备及存储介质

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