CN110427373A - 一种轨迹数据处理方法及装置 - Google Patents

一种轨迹数据处理方法及装置 Download PDF

Info

Publication number
CN110427373A
CN110427373A CN201910669733.7A CN201910669733A CN110427373A CN 110427373 A CN110427373 A CN 110427373A CN 201910669733 A CN201910669733 A CN 201910669733A CN 110427373 A CN110427373 A CN 110427373A
Authority
CN
China
Prior art keywords
point
track
ith
track point
compressed
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.)
Pending
Application number
CN201910669733.7A
Other languages
English (en)
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 University of Aeronautics and Astronautics
Original Assignee
Beijing University of Aeronautics and Astronautics
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 University of Aeronautics and Astronautics filed Critical Beijing University of Aeronautics and Astronautics
Priority to CN201910669733.7A priority Critical patent/CN110427373A/zh
Publication of CN110427373A publication Critical patent/CN110427373A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/22Indexing; Data structures therefor; Storage structures
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/20Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
    • G06F16/29Geographical information databases

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Software Systems (AREA)
  • Image Analysis (AREA)

Abstract

本申请公开一种轨迹数据处理方法及装置,包括:获取起始轨迹点Ps及第i个轨迹点Pi的位置,以及起始轨迹点Ps至第i‑1个轨迹点Pi‑1之间的所有轨迹点形成的多边形交集R′i‑1;其中,i为大于等于1且小于等于N的整数,N表示所有轨迹点的数量;根据起始轨迹点Ps及第i个轨迹点Pi的位置,得到第i个轨迹点Pi与起始轨迹点Ps形成的时空锥及其内接正多边形Ri,并根据多边形交集R′i‑1及内接正多边形Ri,计算得到新的多边形交集R′i,多边形交集R′i用于判断第i个轨迹点Pi是否可作为当前压缩轨迹线段的终点;遍历i的所有取值,当i=N时,得到起始轨迹点Ps至第i个轨迹点Pi之间的所有轨迹点形成的第N条轨迹线段对应的压缩轨迹线段。

Description

一种轨迹数据处理方法及装置
技术领域
本申请涉及数据处理技术领域,尤其涉及一种轨迹数据处理方法及装置。
背景技术
原始轨迹数据的数据量通常比较大,有必要对原始轨迹数据进行压缩。目前,比较常用是采用折线简化压缩方法对原始轨迹数据进行压缩,折线简化压缩方法的做法是:每采集一个轨迹点的数据,就将新增的轨迹点加入之前已采集的轨迹线段中,在进行压缩计算时,对轨迹线段中的所有轨迹点进行压缩计算。
在折线简化压缩方法中,每增加一个轨迹点都要对之前轨迹线段中的所有轨迹点进行计算,即,较早采集的轨迹点会经过多次计算,这使得计算过程复杂度较高,导致数据压缩效率低下,浪费资源。
发明内容
为解决上述技术问题,本申请实施例提供了一种轨迹数据处理方法及装置。
本申请实施例提供的轨迹数据处理方法,包括:
获取起始轨迹点Ps及第i个轨迹点Pi的位置,以及所述起始轨迹点Ps至第i-1个轨迹点Pi-1之间的所有轨迹点形成的多边形交集R′i-1;其中,i为大于等于1且小于等于N的整数,N表示所有轨迹点的数量;
根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,得到所述第i个轨迹点Pi与所述起始轨迹点Ps形成的时空锥及其内接正多边形Ri,并根据所述多边形交集R′i-1及所述内接正多边形Ri,计算得到新的多边形交集R′i,所述多边形交集R′i用于判断所述第i个轨迹点Pi是否可作为当前压缩轨迹线段的终点;
遍历i的所有取值,当i=N时,得到所述起始轨迹点Ps至所述第i个轨迹点Pi之间的所有轨迹点形成的第N条轨迹线段对应的压缩轨迹线段。
在本申请的一种可选实施方式中,所述根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,得到所述第i个轨迹点Pi与所述起始轨迹点Ps形成的时空锥,包括:
根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,在三维时空中,以所述起始轨迹点Ps为坐标原点,构建直角坐标系;
在所述直角坐标系中,以所述第i个轨迹点Pi所在的时间平面t=Pi.t上,以所述第i个轨迹点Pi为圆心,以预设阈值为半径构造同步圆;
以所述同步圆为底面圆,以所述起始轨迹点Ps为顶点,构造时空锥。
在本申请的一种可选实施方式中,所述预设阈值为∈/2或∈;其中,
所述预设阈值为∈/2的情况下,所述同步圆记为O(Pi,∈/2),所述时空锥记为C(Ps,O(Pi,∈/2));
所述预设阈值为∈的情况下,所述同步圆记为O(Pi,∈),所述时空锥记为C(Ps,O(Pi,∈))。
在本申请的一种可选实施方式中,所述方法还包括:
根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,确定时空锥底面圆的投影平面为所述起始轨迹点Ps之后的第一个轨迹点Ps+1所在的时间平面t=Ps+1.t;
对于所述第i个轨迹点Pi,如果i>s+1,则将所述同步圆向所述投影平面t=Ps+1.t上进行投影,得到投影圆 分别表示投影圆的圆心和半径,满足以下公式:Pi c.x=Ps.x+c·(Pi.x-Ps.x),c·(Pi.x-Ps.x),Pi c.t=tc以及其中,
在本申请的一种可选实施方式中,所述方法还包括:通过以下方式计算所述内接正多边形Ri
根据所述投影圆计算其内接正多边形Ri(V,E),其中:
V={v1,...,vm}是多边形的顶点集合,m为多边形的边数,所述多边形Ri(V,E)的顶点在一个以所述第i个轨迹点Pi为坐标原点的极坐标系中定义:
为所述多边形的边集合,所述边集合中的各条边使用其起点进行标记。
在本申请的一种可选实施方式中,所述方法还包括:通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
若所述多边形交集R′i为空集,则所述第i个轨迹点Pi不能作为当前压缩轨迹线段的终点,选择第i-1个轨迹点Pi-1作为压缩轨迹线段的终点;
输出第i条轨迹线段经过压缩后得到的压缩轨迹线段并确定新的起始轨迹点Psn为所述第i-1个轨迹点Pi-1
在本申请的一种可选实施方式中,所述方法还包括:通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
将所述压缩轨迹线段的终点及新的起始轨迹点都确定为插值点Q,其中:
如果所述第i个轨迹点Pi位于所述投影平面t=Pi.t上的多边形交集R′i内,则选取所述第i个轨迹点Pi作为所述插值点Q;
如果且所述第i个轨迹点Pi不在所述投影平面tc=Pi.t上的多边形交集R′i内,则选取所述多边形交集R′i的中心点作为所述插值点Q。
在本申请的一种可选实施方式中,所述方法还包括:通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
若所述第i个轨迹点Pi落在之前轨迹点的时空锥的交集之外,则所述第i个轨迹点Pi不能作为当前压缩轨迹线段的终点,选择第i-1个轨迹点Pi-1作为压缩轨迹线段的终点;
输出第i条轨迹线段经过压缩后得到的压缩轨迹线段并确定新的起始轨迹点Psn为所述第i-1个轨迹点Pi-1
在本申请的一种可选实施方式中,所述方法还包括:所述多边形交集R′i的计算采用如下规则:令表示向量和向量的叉乘积,表示向量左侧的开放半平面;
第一规则:如果或者 成立,则按以下规则移动s步:
其中,g(e)表示边e的编号;
第二规则:如果或者 成立,则边直接移动到与边编号相同的那条边之后的一个位置。
本申请实施例提供的轨迹数据处理装置,包括:
获取单元,用于获取起始轨迹点Ps及第i个轨迹点Pi的位置,以及所述起始轨迹点Ps至第i-1个轨迹点Pi-1之间的所有轨迹点形成的多边形交集R′i-1;其中,i为大于等于1且小于等于N的整数,N表示所有轨迹点的数量;
处理单元,用于根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,得到所述第i个轨迹点Pi与所述起始轨迹点Ps形成的时空锥及其内接正多边形Ri,并根据所述多边形交集R′i-1及所述内接正多边形Ri,计算得到新的多边形交集R′i,所述多边形交集R′i用于判断所述第i个轨迹点Pi是否可作为当前压缩轨迹线段的终点;遍历i的所有取值,当i=N时,得到所述起始轨迹点Ps至所述第i个轨迹点Pi之间的所有轨迹点形成的第N条轨迹线段对应的压缩轨迹线段。
在本申请的一种可选实施方式中,所述处理单元,还用于根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,在三维时空中,以所述起始轨迹点Ps为坐标原点,构建直角坐标系;在所述直角坐标系中,以所述第i个轨迹点Pi所在的时间平面t=Pi.t上,以所述第i个轨迹点Pi为圆心,以预设阈值为半径构造同步圆;以所述同步圆为底面圆,以所述起始轨迹点Ps为顶点,构造时空锥。
在本申请的一种可选实施方式中,所述预设阈值为∈/2或∈;其中,
所述预设阈值为∈/2的情况下,所述同步圆记为O(Pi,∈/2),所述时空锥记为C(Ps,O(Pi,∈/2));
所述预设阈值为∈的情况下,所述同步圆记为O(Pi,∈),所述时空锥记为C(Ps,O(Pi,∈))。
在本申请的一种可选实施方式中,所述处理单元,还用于根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,确定时空锥底面圆的投影平面为所述起始轨迹点Ps之后的第一个轨迹点Ps+1所在的时间平面t=Ps+1.t;对于所述第i个轨迹点Pi,如果i>s+1,则将所述同步圆向所述投影平面t=Ps+1t上进行投影,得到投影圆 分别表示投影圆的圆心和半径,满足以下公式:Pi c.x=Ps.x+c·(Pi.x-Ps.x),Pi c.t=tc以及其中,
在本申请的一种可选实施方式中,所述处理单元,还用于通过以下方式计算所述内接正多边形Ri
根据所述投影圆计算其内接正多边形Ri(V,E),其中:
V={v1,...,vm}是多边形的顶点集合,m为多边形的边数,所述多边形Ri(V,E)的顶点在一个以所述第i个轨迹点Pi为坐标原点的极坐标系中定义:
为所述多边形的边集合,所述边集合中的各条边使用其起点进行标记。
在本申请的一种可选实施方式中,所述处理单元,还用于通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
若所述多边形交集R′i为空集,则所述第i个轨迹点Pi不能作为当前压缩轨迹线段的终点,选择第i-1个轨迹点Pi-1作为压缩轨迹线段的终点;
输出第i条轨迹线段经过压缩后得到的压缩轨迹线段并确定新的起始轨迹点Psn为所述第i-1个轨迹点Pi-1
在本申请的一种可选实施方式中,所述处理单元,还用于通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
将所述压缩轨迹线段的终点及新的起始轨迹点都确定为插值点Q,其中:
如果所述第i个轨迹点Pi位于所述投影平面t=Pi.t上的多边形交集R′i内,则选取所述第i个轨迹点Pi作为所述插值点Q;
如果且所述第i个轨迹点Pi不在所述投影平面tc=Pi.t上的多边形交集R′i内,则选取所述多边形交集R′i的中心点作为所述插值点Q。
在本申请的一种可选实施方式中,所述处理单元,还用于通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
若所述第i个轨迹点Pi落在之前轨迹点的时空锥的交集之外,则所述第i个轨迹点Pi不能作为当前压缩轨迹线段的终点,选择第i-1个轨迹点Pi-1作为压缩轨迹线段的终点;
输出第i条轨迹线段经过压缩后得到的压缩轨迹线段并确定新的起始轨迹点Psn为所述第i-1个轨迹点Pi-1
在本申请的一种可选实施方式中,其特征在于,所述多边形交集R′i的计算采用如下规则:令表示向量和向量的叉乘积,表示向量左侧的开放半平面;
第一规则:如果或者 成立,则按以下规则移动s步:
其中,g(e)表示边e的编号;
第二规则:如果或者 成立,则边直接移动到与边编号相同的那条边之后的一个位置。
采用本申请实施例提供轨迹数据处理方法,针对每个新增的轨迹点数据进行压缩时,只需要获取到前一个轨迹点与起始轨迹点之间的所有轨迹点形成的轨迹线段的内接正多边形交集R′i-1,在R′i-1的基础上,根据新增的轨迹点与初始轨迹点形成的内接正多边形,即可得到新增的轨迹点与起始轨迹点之间的所有轨迹点形成的经过压缩后的轨迹线段,压缩计算的过程只需对新增的轨迹点进行处理,轨迹片段中的所有轨迹点只需要处理一次,大大降低了计算复杂度,提升了数据压缩效率。
附图说明
图1为本申请实施例提供的轨迹数据处理方法的流程示意图;
图2为本申请实施例提供的同步圆和时空锥的示意图;
图3为本申请实施例提供的投影平面与投影圆的示意图;
图4为本申请实施例提供的投影圆的内接正多边形的示意图;
图5为本申请实施例提供的多边形交集计算的示意图;
图6为本申请实施例提供的轨迹数据处理装置的结构组成示意图。
具体实施方式
为了能够更加详尽地了解本申请实施例的特点与技术内容,下面结合附图对本申请实施例的实现进行详细阐述,所附附图仅供参考说明之用,并非用来限定本申请实施例。
图1为本申请实施例提供的轨迹数据处理方法的流程示意图,如图1所示,所述轨迹数据处理方法包括以下步骤:
步骤101:获取起始轨迹点Ps及第i个轨迹点Pi的位置,以及所述起始轨迹点Ps至第i-1个轨迹点Pi-1之间的所有轨迹点形成的多边形交集R′i-1;其中,i为大于等于1且小于等于N的整数,N表示所有轨迹点的数量。
这里,所述起始轨迹点Ps至第i-1个轨迹点Pi-1之间的所有轨迹点形成i-1个内接正多边形交集R′i-1
需要说明的是,所述R′i-1表示所述Ps至第i-1个轨迹点之间的所有轨迹点形成的内接正多边形的交集。
步骤102:根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,得到所述第i个轨迹点Pi与所述起始轨迹点Ps形成的时空锥及其内接正多边形Ri,并根据所述多边形交集R′i-1及所述内接正多边形Ri,计算得到新的多边形交集R′i,所述多边形交集R′i用于判断所述第i个轨迹点Pi是否可作为当前压缩轨迹线段的终点。
本申请实施例中,所述时空锥的形成(或者说构建)可以通过以下方式实现:根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,在三维时空中,以所述起始轨迹点Ps为坐标原点,构建直角坐标系;在所述直角坐标系中,以所述第i个轨迹点Pi所在的时间平面t=Pi.t上,以所述第i个轨迹点Pi为圆心,以预设阈值为半径构造同步圆;以所述同步圆为底面圆,以所述起始轨迹点Ps为顶点,构造时空锥。
进一步,在一种可选实施方式中,所述预设阈值为∈/2或∈;其中,(1)所述预设阈值为∈/2的情况下,所述同步圆记为O(Pi,∈/2),所述时空锥记为C(Ps,O(Pi,∈/2));(2)所述预设阈值为∈的情况下,所述同步圆记为O(Pi,∈),所述时空锥记为C(Ps,O(Pi,∈))。
在一个示例中,根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,在三维时空中,以所述起始轨迹点Ps为坐标原点,构建直角坐标系;在所述直角坐标系中,以所述第i个轨迹点Pi所在的时间平面t=Pi.t上,以所述第i个轨迹点Pi为圆心,以预设阈值∈/2为半径构造同步圆O(Pi,∈/2);以所述同步圆O(Pi,∈/2)为底面圆,以所述起始轨迹点Ps为顶点,构造时空锥C(Ps,O(Pi,∈/2)。
在另一个示例中,根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,在三维时空中,以所述起始轨迹点Ps为坐标原点,构建直角坐标系;在所述直角坐标系中,以所述第i个轨迹点Pi所在的时间平面t=Pi.t上,以所述第i个轨迹点Pi为圆心,以预设阈值∈为半径构造同步圆O(Pi,∈);以所述同步圆O(Pi,∈)为底面圆,以所述起始轨迹点Ps为顶点,构造时空锥C(Ps,O(Pi,∈))。
上述方案中,三维时空具体是x-y-t三维时空,具体地,三维时空的三个维度分别为x,y,t。x和y代表两个空间维度,t代表时间维度,图2给出了同步圆和时空锥的示意图。
本申请实施例中,还包括投影平面的确定及投影圆的计算方法,具体地,根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,确定时空锥底面圆的投影平面为所述起始轨迹点Ps之后的第一个轨迹点Ps+1所在的时间平面t=Ps+1.t;对于所述第i个轨迹点Pi,如果i>s+1,则将所述同步圆向所述投影平面t=Ps+1.t上进行投影,得到投影圆 分别表示投影圆的圆心和半径,满足以下公式:Pi c.x=Ps.x+c·(Pi.x-Ps.x), Pi c.t=tc以及其中,图3给出了投影平面与投影圆的示意图。需要说明的是,所述同步圆可以是O(Pi,∈/2)或者O(Pi,∈)。
进一步,所述方法还包括:通过以下方式计算所述内接正多边形Ri
根据所述投影圆计算其内接正多边形Ri(V,E),其中:
V={v1,...,vm}是多边形的顶点集合,m为多边形的边数,所述多边形Ri(V,E)的顶点在一个以所述第i个轨迹点Pi为坐标原点的极坐标系中定义:
为所述多边形的边集合,所述边集合中的各条边使用其起点进行标记。图4给出了投影圆的内接正多边形的示意图。
本申请实施例中,还提供了多边形交集的加速计算方法,具体地,参数图5,所述多边形交集R′i的计算采用如下规则:令表示向量和向量的叉乘积,表示向量左侧的开放半平面;
第一规则:如果或者 成立,则按以下规则移动s步:
其中,g(e)表示边e的编号;
第二规则:如果或者 成立,则边直接移动到与边编号相同的那条边之后的一个位置。
步骤103:遍历i的所有取值,当i=N时,得到所述起始轨迹点Ps至所述第i个轨迹点Pi之间的所有轨迹点形成的第N条轨迹线段对应的压缩轨迹线段。
本申请实施例中,所述方法还包括:确定压缩轨迹线段的终点及新的起始轨迹点,具体地,可以采用以下任意一种方式来确定压缩轨迹线段的终点及新的起始轨迹点:
方式一:若所述多边形交集R′i为空集,则所述第i个轨迹点Pi不能作为当前压缩轨迹线段的终点,选择第i-1个轨迹点Pi-1作为压缩轨迹线段的终点;输出第i条轨迹线段经过压缩后得到的压缩轨迹线段并确定新的起始轨迹点Psn为所述第i-1个轨迹点Pi-1
方式二:将所述压缩轨迹线段的终点及新的起始轨迹点都确定为插值点Q,其中:策略1)如果所述第i个轨迹点Pi位于所述投影平面t=Pi.t上的多边形交集R′i内,则选取所述第i个轨迹点Pi作为所述插值点Q;策略2)如果且所述第i个轨迹点Pi不在所述投影平面tc=Pi.t上的多边形交集R′i内,则选取所述多边形交集R′i的中心点作为所述插值点Q。
需要说明的是,一般情况下,如果所述投影平面t=Ps+1.t≠Pi.t,则首先把多边形R′i映射到平面t=Pi.t上,然后,再使用前述策略1)或策略2)来确定所述插值点Q。
方式三:若所述第i个轨迹点Pi落在之前轨迹点的时空锥的交集之外,则所述第i个轨迹点Pi不能作为当前压缩轨迹线段的终点,选择第i-1个轨迹点Pi-1作为压缩轨迹线段的终点;输出第i条轨迹线段经过压缩后得到的压缩轨迹线段并确定新的起始轨迹点Psn为所述第i-1个轨迹点Pi-1
采用本申请实施例提供轨迹数据处理方法,针对每个新增的轨迹点数据进行压缩时,只需要获取到前一个轨迹点与起始轨迹点之间的所有轨迹点形成的轨迹线段的内接正多边形交集R′i-1,在R′i-1的基础上,根据新增的轨迹点与初始轨迹点形成的内接正多边形,即可得到新增的轨迹点与起始轨迹点之间的所有轨迹点形成的经过压缩后的轨迹线段,压缩计算的过程只需对新增的轨迹点进行处理,轨迹片段中的所有轨迹点只需要处理一次,大大降低了计算复杂度,提升了数据压缩效率。
图6为本申请实施例提供的轨迹数据处理装置的结构组成示意图,如图6所示,所述轨迹数据处理装置包括:
获取单元601,用于获取起始轨迹点Ps及第i个轨迹点Pi的位置,以及所述起始轨迹点Ps至第i-1个轨迹点Pi-1之间的所有轨迹点形成的多边形交集R′i-1;其中,i为大于等于1且小于等于N的整数,N表示所有轨迹点的数量;
处理单元602,用于根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,得到所述第i个轨迹点Pi与所述起始轨迹点Ps形成的时空锥及其内接正多边形Ri,并根据所述多边形交集R′i-1及所述内接正多边形Ri,计算得到新的多边形交集R′i,所述多边形交集R′i用于判断所述第i个轨迹点Pi是否可作为当前压缩轨迹线段的终点;遍历i的所有取值,当i=N时,得到所述起始轨迹点Ps至所述第i个轨迹点Pi之间的所有轨迹点形成的第N条轨迹线段对应的压缩轨迹线段。
在本申请一种可选实施方式中,所述处理单元602,还用于根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,在三维时空中,以所述起始轨迹点Ps为坐标原点,构建直角坐标系;在所述直角坐标系中,以所述第i个轨迹点Pi所在的时间平面t=Pi.t上,以所述第i个轨迹点Pi为圆心,以预设阈值为半径构造同步圆;以所述同步圆为底面圆,以所述起始轨迹点Ps为顶点,构造时空锥。
进一步,可选地,所述预设阈值为∈/2或∈;其中,所述预设阈值为∈/2的情况下,所述同步圆记为O(Pi,∈/2),所述时空锥记为C(Ps,O(Pi,∈/2));所述预设阈值为∈的情况下,所述同步圆记为O(Pi,∈),所述时空锥记为C(Ps,O(Pi,∈))。
在本申请一种可选实施方式中,所述处理单元602,还用于根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,确定时空锥底面圆的投影平面为所述起始轨迹点Ps之后的第一个轨迹点Ps+1所在的时间平面t=Ps+1.t;对于所述第i个轨迹点Pi,如果i>s+1,则将所述同步圆向所述投影平面t=Ps+1.t上进行投影,得到投影圆 分别表示投影圆的圆心和半径,满足以下公式:Pi c.x=Ps.x+c·(Pi.x-Ps.x), Pi ct=tc以及其中,
在本申请一种可选实施方式中,所述处理单元602,还用于通过以下方式计算所述内接正多边形Ri
根据所述投影圆计算其内接正多边形Ri(V,E),其中:
V={v1,...,vm}是多边形的顶点集合,m为多边形的边数,所述多边形Ri(V,E)的顶点在一个以所述第i个轨迹点Pi为坐标原点的极坐标系中定义:
为所述多边形的边集合,所述边集合中的各条边使用其起点进行标记。
在本申请一种可选实施方式中,所述处理单元602,还用于通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
若所述多边形交集R′i为空集,则所述第i个轨迹点Pi不能作为当前压缩轨迹线段的终点,选择第i-1个轨迹点Pi-1作为压缩轨迹线段的终点;
输出第i条轨迹线段经过压缩后得到的压缩轨迹线段并确定新的起始轨迹点Psn为所述第i-1个轨迹点Pi-1
在本申请一种可选实施方式中,所述处理单元602,还用于通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
将所述压缩轨迹线段的终点及新的起始轨迹点都确定为插值点Q,其中:
如果所述第i个轨迹点Pi位于所述投影平面t=Pi.t上的多边形交集R′i内,则选取所述第i个轨迹点Pi作为所述插值点Q;
如果且所述第i个轨迹点Pi不在所述投影平面tc=Pi.t上的多边形交集R′i内,则选取所述多边形交集R′i的中心点作为所述插值点Q。
在本申请一种可选实施方式中,所述处理单元602,还用于通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
若所述第i个轨迹点Pi落在之前轨迹点的时空锥的交集之外,则所述第i个轨迹点Pi不能作为当前压缩轨迹线段的终点,选择第i-1个轨迹点Pi-1作为压缩轨迹线段的终点;
输出第i条轨迹线段经过压缩后得到的压缩轨迹线段并确定新的起始轨迹点Psn为所述第i-1个轨迹点Pi-1
上述方案中,所述多边形交集R′i的计算采用如下规则:令表示向量和向量的叉乘积,表示向量左侧的开放半平面;
第一规则:如果或者 成立,则按以下规则移动s步:
其中,g(e)表示边e的编号;
第二规则:如果或者 成立,则边直接移动到与边编号相同的那条边之后的一个位置。
本领域技术人员应当理解,图6所示的轨迹数据处理装置中的各单元的实现功能可参照前述轨迹数据处理方法的相关描述而理解。图6所示的轨迹数据处理装置中的各单元的功能可通过运行于处理器上的程序而实现,也可通过具体的逻辑电路而实现。具体实现时,所述轨迹数据处理装置中的各单元可由所述终端中的处理器,比如中央处理器(Central Processing Unit,CPU)、数字信号处理器(Digital Signal Processor,DSP)、微控制单元(Microcontroller Unit,MCU)或可编程门阵列(Field-Programmable GateArray,FPGA)等实现。
本申请实施例还提供了一种计算机存储介质,具体为计算机可读存储介质。其上存储有计算机指令,作为第一种实施方式,在计算机存储介质位于终端时,该计算机指令被处理器执行时实现上述轨迹数据处理方法中的任意步骤。
需要说明的是:本申请实施例所记载的技术方案之间,在不冲突的情况下,可以任意组合。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。

Claims (18)

1.一种轨迹数据处理方法,其特征在于,所述方法包括:
获取起始轨迹点Ps及第i个轨迹点Pi的位置,以及所述起始轨迹点Ps至第i-1个轨迹点Pi-1之间的所有轨迹点形成的多边形交集R′i-1;其中,i为大于等于1且小于等于N的整数,N表示所有轨迹点的数量;
根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,得到所述第i个轨迹点Pi与所述起始轨迹点Ps形成的时空锥及其内接正多边形Ri,并根据所述多边形交集R′i-1及所述内接正多边形Ri,计算得到新的多边形交集R′i,所述多边形交集R′i用于判断所述第i个轨迹点Pi是否可作为当前压缩轨迹线段的终点;
遍历i的所有取值,当i=N时,得到所述起始轨迹点Ps至所述第i个轨迹点Pi之间的所有轨迹点形成的第N条轨迹线段对应的压缩轨迹线段。
2.根据权利要求1所述的方法,其特征在于,所述根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,得到所述第i个轨迹点Pi与所述起始轨迹点Ps形成的时空锥,包括:
根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,在三维时空中,以所述起始轨迹点Ps为坐标原点,构建直角坐标系;
在所述直角坐标系中,以所述第i个轨迹点Pi所在的时间平面t=Pi.t上,以所述第i个轨迹点Pi为圆心,以预设阈值为半径构造同步圆;
以所述同步圆为底面圆,以所述起始轨迹点Ps为顶点,构造时空锥。
3.根据权利要求2所述的方法,其特征在于,所述预设阈值为∈/2或∈;其中,
所述预设阈值为∈/2的情况下,所述同步圆记为O(Pi,∈/2),所述时空锥记为C(Ps,O(Pi,∈/2));
所述预设阈值为∈的情况下,所述同步圆记为O(Pi,∈),所述时空锥记为C(Ps,O(Pi,∈))。
4.根据权利要求2或3所述的方法,其特征在于,所述方法还包括:
根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,确定时空锥底面圆的投影平面为所述起始轨迹点Ps之后的第一个轨迹点Ps+1所在的时间平面t=Ps+1.t;
对于所述第i个轨迹点Pi,如果i>s+1,则将所述同步圆向所述投影平面t=Ps+1.t上进行投影,得到投影圆 分别表示投影圆的圆心和半径,满足以下公式:Pi c.x=Ps.x+c·(Pi.x-Ps.x), Pi c.t=tc以及其中,
5.根据权利要求4所述的方法,其特征在于,所述方法还包括:通过以下方式计算所述内接正多边形Ri
根据所述投影圆计算其内接正多边形Ri(V,E),其中:
V={v1,...,vm}是多边形的顶点集合,m为多边形的边数,所述多边形Ri(V,E)的顶点在一个以所述第i个轨迹点Pi为坐标原点的极坐标系中定义:
为所述多边形的边集合,所述边集合中的各条边使用其起点进行标记。
6.根据权利要求1至5中任一项所述的方法,其特征在于,所述方法还包括:通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
若所述多边形交集R′i为空集,则所述第i个轨迹点Pi不能作为当前压缩轨迹线段的终点,选择第i-1个轨迹点Pi-1作为压缩轨迹线段的终点;
输出第i条轨迹线段经过压缩后得到的压缩轨迹线段并确定新的起始轨迹点Psn为所述第i-1个轨迹点Pi-1
7.根据权利要求1至5中任一项所述的方法,其特征在于,所述方法还包括:通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
将所述压缩轨迹线段的终点及新的起始轨迹点都确定为插值点Q,其中:
如果所述第i个轨迹点Pi位于所述投影平面t=Pi.t上的多边形交集R′i内,则选取所述第i个轨迹点Pi作为所述插值点Q;
如果且所述第i个轨迹点Pi不在所述投影平面tc=Pi.t上的多边形交集R′i内,则选取所述多边形交集R′i的中心点作为所述插值点Q。
8.根据权利要求1至5中任一项所述的方法,其特征在于,所述方法还包括:通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
若所述第i个轨迹点Pi落在之前轨迹点的时空锥的交集之外,则所述第i个轨迹点Pi不能作为当前压缩轨迹线段的终点,选择第i-1个轨迹点Pi-1作为压缩轨迹线段的终点;
输出第i条轨迹线段经过压缩后得到的压缩轨迹线段并确定新的起始轨迹点Psn为所述第i-1个轨迹点Pi-1
9.根据权利要求1所述的方法,其特征在于,所述多边形交集R′i的计算采用如下规则:令表示向量和向量的叉乘积,表示向量左侧的开放半平面;
第一规则:如果或者 成立,则按以下规则移动s步:
其中,g(e)表示边e的编号;
第二规则:如果或者 成立,则边直接移动到与边编号相同的那条边之后的一个位置。
10.一种轨迹数据处理装置,其特征在于,所述轨迹数据处理装置包括:
获取单元,用于获取起始轨迹点Ps及第i个轨迹点Pi的位置,以及所述起始轨迹点Ps至第i-1个轨迹点Pi-1之间的所有轨迹点形成的多边形交集R′i-1;其中,i为大于等于1且小于等于N的整数,N表示所有轨迹点的数量;
处理单元,用于根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,得到所述第i个轨迹点Pi与所述起始轨迹点Ps形成的时空锥及其内接正多边形Ri,并根据所述多边形交集R′i-1及所述内接正多边形Ri,计算得到新的多边形交集R′i,所述多边形交集R′i用于判断所述第i个轨迹点Pi是否可作为当前压缩轨迹线段的终点;遍历i的所有取值,当i=N时,得到所述起始轨迹点Ps至所述第i个轨迹点Pi之间的所有轨迹点形成的第N条轨迹线段对应的压缩轨迹线段。
11.根据权利要求10所述的装置,其特征在于,所述处理单元,还用于根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,在三维时空中,以所述起始轨迹点Ps为坐标原点,构建直角坐标系;在所述直角坐标系中,以所述第i个轨迹点Pi所在的时间平面t=Pi.t上,以所述第i个轨迹点Pi为圆心,以预设阈值为半径构造同步圆;以所述同步圆为底面圆,以所述起始轨迹点Ps为顶点,构造时空锥。
12.根据权利要求11所述的装置,其特征在于,所述预设阈值为∈/2或∈;其中,
所述预设阈值为∈/2的情况下,所述同步圆记为O(Pi,∈/2),所述时空锥记为C(Ps,O(Pi,∈/2));
所述预设阈值为∈的情况下,所述同步圆记为O(Pi,∈),所述时空锥记为C(Ps,O(Pi,∈))。
13.根据权利要求11或12所述的装置,其特征在于,所述处理单元,还用于根据所述起始轨迹点Ps及所述第i个轨迹点Pi的位置,确定时空锥底面圆的投影平面为所述起始轨迹点Ps之后的第一个轨迹点Ps+1所在的时间平面t=Ps+1.t;对于所述第i个轨迹点Pi,如果i>s+1,则将所述同步圆向所述投影平面t=Ps+1.t上进行投影,得到投影圆 分别表示投影圆的圆心和半径,满足以下公式:Pi c.x=Ps.x+c·(Pi.x-Ps.x), Pi c.t=tc以及其中,
14.根据权利要求13所述的装置,其特征在于,所述处理单元,还用于通过以下方式计算所述内接正多边形Ri
根据所述投影圆计算其内接正多边形Ri(V,E),其中:
V={v1,...,vm}是多边形的顶点集合,m为多边形的边数,所述多边形Ri(V,E)的顶点在一个以所述第i个轨迹点Pi为坐标原点的极坐标系中定义:
为所述多边形的边集合,所述边集合中的各条边使用其起点进行标记。
15.根据权利要求10至14中任一项所述的装置,其特征在于,所述处理单元,还用于通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
若所述多边形交集R′i为空集,则所述第i个轨迹点Pi不能作为当前压缩轨迹线段的终点,选择第i-1个轨迹点Pi-1作为压缩轨迹线段的终点;
输出第i条轨迹线段经过压缩后得到的压缩轨迹线段并确定新的起始轨迹点Psn为所述第i-1个轨迹点Pi-1
16.根据权利要求10至14中任一项所述的装置,其特征在于,所述处理单元,还用于通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
将所述压缩轨迹线段的终点及新的起始轨迹点都确定为插值点Q,其中:
如果所述第i个轨迹点Pi位于所述投影平面t=Pi.t上的多边形交集R′i内,则选取所述第i个轨迹点Pi作为所述插值点Q;
如果且所述第i个轨迹点Pi不在所述投影平面tc=Pi.t上的多边形交集R′i内,则选取所述多边形交集R′i的中心点作为所述插值点Q。
17.根据权利要求10至14中任一项所述的装置,其特征在于,所述处理单元,还用于通过以下方式确定压缩轨迹线段的终点及新的起始轨迹点:
若所述第i个轨迹点Pi落在之前轨迹点的时空锥的交集之外,则所述第i个轨迹点Pi不能作为当前压缩轨迹线段的终点,选择第i-1个轨迹点Pi-1作为压缩轨迹线段的终点;
输出第i条轨迹线段经过压缩后得到的压缩轨迹线段并确定新的起始轨迹点Psn为所述第i-1个轨迹点Pi-1
18.根据权利要求10所述的装置,其特征在于,所述多边形交集R′i的计算采用如下规则:令表示向量和向量的叉乘积,表示向量左侧的开放半平面;
第一规则:如果或者 成立,则按以下规则移动s步:
其中,g(e)表示边e的编号;
第二规则:如果或者 成立,则边直接移动到与边编号相同的那条边之后的一个位置。
CN201910669733.7A 2019-07-24 2019-07-24 一种轨迹数据处理方法及装置 Pending CN110427373A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910669733.7A CN110427373A (zh) 2019-07-24 2019-07-24 一种轨迹数据处理方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910669733.7A CN110427373A (zh) 2019-07-24 2019-07-24 一种轨迹数据处理方法及装置

Publications (1)

Publication Number Publication Date
CN110427373A true CN110427373A (zh) 2019-11-08

Family

ID=68410542

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910669733.7A Pending CN110427373A (zh) 2019-07-24 2019-07-24 一种轨迹数据处理方法及装置

Country Status (1)

Country Link
CN (1) CN110427373A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111258558A (zh) * 2020-02-13 2020-06-09 苏宁智能终端有限公司 运动轨迹的着色处理方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102299718A (zh) * 2010-06-25 2011-12-28 汉王科技股份有限公司 手写终端轨迹的压缩和解压方法
CN105868553A (zh) * 2016-03-28 2016-08-17 北京航空航天大学 轨迹数据处理方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102299718A (zh) * 2010-06-25 2011-12-28 汉王科技股份有限公司 手写终端轨迹的压缩和解压方法
CN105868553A (zh) * 2016-03-28 2016-08-17 北京航空航天大学 轨迹数据处理方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
XUELIAN LIN等: "One-Pass Trajectory Simplification Using the Synchronous Euclidean Distance", 《ARXIV》 *
左一萌等: "路网感知的在线轨迹压缩方法", 《软件学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111258558A (zh) * 2020-02-13 2020-06-09 苏宁智能终端有限公司 运动轨迹的着色处理方法及装置
CN111258558B (zh) * 2020-02-13 2023-06-30 苏宁智能终端有限公司 运动轨迹的着色处理方法及装置

Similar Documents

Publication Publication Date Title
Schöps et al. Surfelmeshing: Online surfel-based mesh reconstruction
WO2020119684A1 (zh) 一种3d导航语义地图更新方法、装置及设备
JP6760957B2 (ja) 3次元モデリング方法及び装置
US20130127847A1 (en) System and Method for Interactive Image-based Modeling of Curved Surfaces Using Single-view and Multi-view Feature Curves
KR20170068462A (ko) 에지를 이용한 3차원 모델 생성
AU2020401367A1 (en) 3-D reconstruction using augmented reality frameworks
CN110276768B (zh) 图像分割方法、图像分割装置、图像分割设备及介质
CN107578467B (zh) 一种医疗器械三维建模方法及装置
CN107808388B (zh) 包含运动目标的图像处理方法、装置及电子设备
CN110648363A (zh) 相机姿态确定方法、装置、存储介质及电子设备
JP2022509329A (ja) 点群融合方法及び装置、電子機器、コンピュータ記憶媒体並びにプログラム
CN115439607A (zh) 一种三维重建方法、装置、电子设备及存储介质
KR101593316B1 (ko) 스테레오 카메라를 이용한 3차원 모델 재구성 방법 및 장치
KR102158390B1 (ko) 영상 처리 방법 및 장치
CN113628343B (zh) 三维网格的合并处理方法和装置、存储介质
CN113658203A (zh) 建筑物三维轮廓提取及神经网络的训练方法和装置
CN113920275A (zh) 三角网格构建方法、装置、电子设备及可读存储介质
CN117635875B (zh) 一种三维重建方法、装置及终端
US8264487B2 (en) Method for converting polygonal surfaces to levelsets
EP4246452A1 (en) Three-dimensional point cloud densification device, three-dimensional point cloud densification method, and program
CN110427373A (zh) 一种轨迹数据处理方法及装置
CN111402429B (zh) 一种尺度还原、三维重建方法、系统、存储介质及设备
WO2024002064A1 (zh) 三维模型构建方法、装置、电子设备及存储介质
CN114627206A (zh) 网格绘制方法、装置、电子设备及计算机可读存储介质
CN113379826A (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20191108

RJ01 Rejection of invention patent application after publication