CN113806463B

CN113806463B - 一种基于时空金字塔匹配的轨迹相似度计算方法

Info

Publication number: CN113806463B
Application number: CN202111039235.8A
Authority: CN
Inventors: 李莉; 张建宇; 戴帅夫; 刘丙双; 蒋志鹏; 周模
Original assignee: Beijing Information Science and Technology University
Current assignee: Beijing Information Science and Technology University
Priority date: 2021-09-06
Filing date: 2021-09-06
Publication date: 2023-04-28
Anticipated expiration: 2041-09-06
Also published as: CN113806463A

Abstract

本发明公开了一种基于时空金字塔匹配的轨迹相似度计算方法，包括以下步骤：S1、根据用户号码和给定时间查询，获取用户原始轨迹数据；S2、对原始轨迹数据预处理；S3、基于Google S2算法，对经纬度进行空间编码，并根据不同时间粒度，建立时间编码；S4、对每个时间粒度以及空间粒度，计算每个时空粒度下的轨迹相似度；S5、对不同的时空粒度赋予不同权重，计算总相似度。

Description

一种基于时空金字塔匹配的轨迹相似度计算方法

技术领域

本发明涉及时空数据分析及挖掘技术领域，更具体的说是涉及一种基于时空金字塔匹配的轨迹相似度计算方法。

背景技术

目前，随着物联网技术以及移动互联网的飞速发展，无处不在的定位装置如传感器、GPS定位装置、智能卡等产生的了海量的轨迹数据。轨迹挖掘已经成为计算机应用中城市计算领域的研究热点。基于移动用户的轨迹挖掘能够为社会公共安全领域、智能交通领域、人口统计部门等提供关键的技术支撑。这背后就是依靠时空大数据应用技术的支撑。

移动用户轨迹数据主要来源于被授权“定位”权限的APP的GPS轨迹数据以及运营商的移动信令数据。GPS轨迹数据是主动式数据，即只有在用户触发的情况下，才能产生数据，这些数据通过用户端的浏览器、电脑客户端软件、手机App等最终汇聚在互联网企业中。这就决定了互联网数据只能是片段式的，而且样本覆盖受到用户群体喜好偏差的影响。另外，由于企业间的竞争等因素，这些数据难以融为一体，分析利用这些数据容易形成片面倾向性结果，而移动信令数据更加客观、全面、长期持续在线，数据应用结果更加科学可信。由于GPS数据是均匀采样，而移动信令轨迹数据是非均匀采样，这就造成了无法用传统的欧式距离、巴氏距离等度量方法直接用来度量基于移动信令的轨迹相似度。DTW距离是根据动态规划的方法计使得两个序列中所有对应点的距离之和最小，虽然不受到轨迹点数是否相同的限制，但是计算复杂度高；最长公共子序列LCSS方法对噪声不敏感，但是和DTW类似，只考虑了轨迹的序列性，而没有考虑轨迹的时间特性。Frechet距离是一种基于空间路径相似度描述方法，其着重将路径空间距离考虑进去，使得其对于有一定空间时序的曲线相似度评价效率相比之下较高，对于信令轨迹数据，存在乒乓效应(短时间内轨迹点在相邻基站之间飘移)，该方法并不适用于信令数据。

因此，针对移动信令数据，本发明提出一种计算复杂度低、同时考虑轨迹的时间性和空间性的轨迹相似度度量方法。

发明内容

有鉴于此，本发明提供了一种基于时空金字塔匹配的轨迹相似度计算方法；基于Google S2编码将城市将经纬度转换成前缀匹配优化的字符串，从而实现城市地理信息网格化，并利用网格建立空间索引，其优势在于利用网格一方面能够快速基于距离对空间位置信息进行检索，另外一方面在一定程度上降低了基于基站的位置信息的噪声；本发明提出一种基于时空金字塔匹配的轨迹相似性测量方法，其核心思想是在时间和空间维度上对轨迹进行不同粒度上的划分，根据不同的时空粒度赋予不同的权重，粒度越粗权重越大，粒度越细权重越小，因此我们的方法同时考虑了轨迹的时间特性和空间特性。

为了实现上述目的，本发明采用如下技术方案：

一种基于时空金字塔匹配的轨迹相似度计算方法，所述方法包括以下步骤：

S1、根据用户号码和给定时间查询，获取用户原始轨迹数据；

S2、对原始轨迹数据预处理；

S3、基于Google S2算法，对经纬度进行空间编码，并根据不同时间粒度，建立时间编码；

S4、对每个时间粒度以及空间粒度，计算每个时空粒度下的轨迹相似度；

S5、对不同的时空粒度赋予不同权重，计算总相似度。

优选的，所述步骤S1具体包括：

用户在给定时间段内的轨迹数据Tr为：

Tr＝{(p_i,t_i),1≤i≤n,t_i＜t_i+1}

其中，p_i为用户在t_i时刻的位置，p_i＝(lon_i,lat_i)，lon_i为用户经度，lat_i为用户纬度。

优选的，所述步骤S2具体包括：

按照时间排序、去除冗余数据、去噪、去除乒乓数据。

优选的，所述步骤S3具体包括：

S31、采用Google S2的地图网格化技术将经纬度p_i＝(lon_i,lat_i)转化为网格，进行空间编码,不同的层级，空间编码也不相同，例如点(116.34511,39.949147)的11级的编码为“51496a4”，15级的编码为“51496a3ec”；

S32、根据不同时间粒度，依据下列公式对时间进行编码

其中temporalScale为时间粒度，t为一天中的长整形时间，t∈[0,86400)，单位为秒。

S33、时空编码，利用空间编码和时间编码，每个位置点(p_i,t_i)编码为为该位置点在第s层的空间编码，是第l层的时间编码。经过3D 时空编码后，用户的轨迹可以表示为其中 n为位置点数，S为空间层数，L为时间层数。

优选的，所述步骤S4具体包括：

S41、记P和Q分别为两条轨迹，位置点数分别为n和m，

则表示轨迹P第l层的时间编码集合，标识轨迹Q第l层的时间编码集合；

S42、时间金字塔匹配：

表示轨迹P和Q在第l层相同的时间编码集合，

则在第l层，轨迹P和Q的相似度可以表示为

其中，分别表示轨迹P和Q在时间编码上的子轨迹，为轨迹P和Q在第l层时间、时间编码上的空间相似度；

S43、空间金字塔匹配：

表示轨迹P和Q的子轨迹在第s级S2编码中具有相同编码的位置点数。和表示P和Q子轨迹的位置点数，κ_S为空间权重，S2编码的级数越高，空间粒度越细，权重越高，

优选的，所述步骤S5具体包括：

轨迹P和Q的相似度为：

其中，ω₁表示第l层时间编码上的权重，ω_l的计算方法为：

经由上述的技术方案可知，与现有技术相比，本发明公开提供了一种基于时空金字塔匹配模型的轨迹相似度计算方法；基于Google S2编码将城市将经纬度转换成前缀匹配优化的字符串，从而实现城市地理信息网格化，并利用网格建立空间索引，其优势在于利用网格一方面能够快速基于距离对空间位置信息进行检索，另外一方面在一定程度上降低了基于基站的位置信息的噪声；本发明提出一种基于时空金字塔匹配的轨迹相似性测量方法，其核心思想是在时间和空间维度上对轨迹进行不同粒度上的划分，根据不同的时空粒度赋予不同的权重，粒度越粗权重越大，粒度越细权重越小，因此我们的方法同时考虑了轨迹的时间特性和空间特性。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。

图1附图为本发明提供的流程结构示意图。

图2附图为本发明提供的地图网格化结构示意图。

图3附图为本发明提供的准确率对比示意图。

图4附图为本发明提供的时间消耗对比示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明实施例公开了一种基于时空金字塔匹配的轨迹相似度计算方法，所述方法包括以下步骤：

S2、对原始轨迹数据预处理；

S5、对不同的时空粒度赋予不同权重，计算总相似度。

为进一步优化上述技术方案，步骤S1具体包括：

用户在给定时间段内的轨迹数据Tr为：

Tr＝{(p_i,t_i),1≤i≤n,t_i＜t_i+1}

为进一步优化上述技术方案，步骤S2具体包括：

按照时间排序、去除冗余数据、去噪、去除乒乓数据。

为进一步优化上述技术方案，步骤S3具体包括：

S31、采用Google S2的地图网格化技术将经纬度p_i＝(lon_i,lat_i)转化为网格，进行空间编码；

S32、根据不同时间粒度，依据下列公式对时间进行编码

为进一步优化上述技术方案，步骤S4具体包括：

S41、记P和Q分别为两条轨迹，位置点数分别为n和m，

S42、时间金字塔匹配：

表示轨迹P和Q在第l层相同的时间编码集合，

则在第l层，轨迹P和Q的相似度可以表示为

S43、空间金字塔匹配：

为进一步优化上述技术方案，步骤S5具体包括：

轨迹P和Q的相似度为：

其中，ω₁表示第l层时间编码上的权重，ω_l的计算方法为：

图3附图为本申请的计算准确率与其他计算方式准确率的对比图，其中 STPM为本申请的计算准确率，从图中可以直观的得出，本申请的准确率要明显优于现有技术的准确率。

图4附图为本申请与现有技术的计算时间消耗对比图，其中STPM为本发明的方法，从该图可以看出，随着位置点数的增加，STPM的时间消耗比较稳定，而其他方法如DTW、LCSS、Frechet等方法随着位置点数的增加其时间消耗呈指数增长趋势，证明本发明的方法具有较低的计算复杂度。

因此，本发明具有以下优点：

根据不同的时空粒度赋予不同的权重，粒度越粗权重越大，粒度越细权重越小，因此我们的方法同时考虑了轨迹的时间特性和空间特性。

基于信令的轨迹数据受终端设备、信号强度、运营商等因素的影响具有低采样率、采样不均匀、噪声大等特征，本发明所提出的算法充分考虑了这些因素，利用不同空间粒度和时间粒度的权重组合来衡量轨迹之间的相似性，因此能够容忍低采样率、采样不均匀以及噪声带来的影响。

LCSS、DTW、Frechet、EDR、ERP等轨迹相似度计算方法都需要计算两两位置点的球面距离，其计算复杂度为O(|P|×|Q|)，|P|和|Q|分别为轨迹P和轨迹Q 的长度。时空金字塔匹配的计算复杂度为O(L×S)，其中L和S分别为时间和空间的层级数量，与轨迹长度无关，因此计算复杂度远远低于常用算法。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims

1.一种基于时空金字塔匹配的轨迹相似度计算方法，其特征在于，所述方法包括以下步骤：

所述步骤S1具体包括：

用户在给定时间段内的轨迹数据Tr为：

Tr＝{(p_i,t_i),1≤i≤n,t_i＜t_i+1}

其中，p_i为用户在t_i时刻的位置，p_i＝(lon_i,lat_i)，lon_i为用户经度，lat_i为用户纬度；

S2、对原始轨迹数据预处理；

S3、基于GoogleS2算法，对经纬度进行空间编码，并根据不同时间粒度，建立时间编码；

所述步骤S3具体包括：

S31、采用GoogleS2的地图网格化技术将经纬度p_i＝(lon_i,lat_i)转化为网格，进行空间编码；

S32、根据不同时间粒度，依据下列公式对时间进行编码

其中temporalScale为时间粒度，t为一天中的长整形时间，t∈[0,86400)，单位为秒；

S33、时空编码，利用空间编码和时间编码，每个位置点(p_i,t_i)编码为为该位置点在第s层的空间编码，是第l层的时间编码；经过3D时空编码后，用户的轨迹表示为其中n为位置点数，S为空间层数，L为时间层数；

所述步骤S4具体包括：

S41、记P和Q分别为两条轨迹，位置点数分别为n和m，

S42、时间金字塔匹配：

表示轨迹P和Q在第l层相同的时间编码集合，

则在第l层，轨迹P和Q的相似度表示为

其中，分别表示轨迹P和Q在时间编码t～上的子轨迹，为轨迹P和Q在第l层时间、时间编码t～上的空间相似度；

S43、空间金字塔匹配：

表示轨迹P和Q的子轨迹在第s级GoogleS2编码中具有相同编码的位置点数；和表示P和Q子轨迹的位置点数，κ_S为空间权重，GoogleS2编码的级数越高，空间粒度越细，权重越高，

S5、对不同的时空粒度赋予不同权重，计算总相似度；

所述步骤S5具体包括：

轨迹P和Q的相似度为：

其中，ω₁表示第l层时间编码上的权重，ω_l的计算方法为：

2.根据权利要求1所述的一种基于时空金字塔匹配的轨迹相似度计算方法，其特征在于，所述步骤S2具体包括：

按照时间排序、去除冗余数据、去噪、去除乒乓数据。