一种联合多元时空轨迹的群体行为模式挖掘方法

1.本发明涉及一种联合多元时空轨迹的群体行为模式挖掘方法，属于社交网络分析技术领域。


背景技术：

2.随着互联网技术以及移动终端技术的高速发展，各类社交平台得到蓬勃发展，这也为犯罪分子提供了便利，犯罪团伙往往通过各种社交平台与通信手段进行沟通，进行有组织，有计划的犯罪行动。在实际案件侦破任务中，为了侦破暴力、恐怖等有预谋的犯罪案件，通常根据某个犯罪嫌疑人的信息，进一步挖掘他在社交网络中所属的犯罪团伙，这一过程需要消耗大量的人力物力。
3.随着gps定位系统、5g等先进前沿信息技术的飞速发展，通常能够获取每一个犯罪分子的多类型时空轨迹数据，大体可以分为三类：（1）基于gps定位的轨迹数据，如飞机、车辆、船舶、手机等。（2）基于感知设备的轨迹数据，如人脸识别、车牌识别、酒店入住、网吧上网等。（3）基于社交网络的轨迹数据，犯罪分子在社交网络中的发文、打卡、签到等行为都记录了相应的位置信息。然而，目前多数方法不能对犯罪分子的时空轨迹与社交关系网络进行有效建模，未能有效挖掘犯罪团伙分工作案的行为之间的隐性关系，导致最终构建的关系网络在节点分类、关系预测的准确性不高。因此，如何基于多元时空轨迹对犯罪群体进行有效建模，分析作案人员与犯罪团伙间的关系，挖掘并推断其作案的行为之间的隐性关系，是目前亟待解决的科学问题。


技术实现要素：

4.本发明提供了一种联合多元时空轨迹的群体行为模式挖掘方法，以用于解决当前基于轨迹的群体行为挖掘方法无法有效利用地理实体信息，只能提取空间上的轨迹重合关系等现象；本发明增加驻留点地理实体信息关联度，引入路网信息对轨迹进行平滑降噪，去除冗余信息，同时将轨迹间的距离计算转化为路网关系查询，提升检测效率。
5.本发明的技术方案是：一种联合多元时空轨迹的群体行为模式挖掘方法，提取轨迹驻留点，联合地理实体数据建立空间场所信息索引，从而赋予轨迹场所及事件信息；依据路网信息，将轨迹拟合到路网不同路段，实现轨迹的降噪及路网索引，将轨迹相似度计算转为路网路段相似度查询；同时使用驻留点地理信息相似度和轨迹路网相似度联合构建轨迹相关性模型；使用无监督的图神经网络挖掘轨迹深层关联；所述方法的具体步骤如下：步骤s1：驻留点检测，将轨迹建模为：其中，x、y是空间坐标，t是当前点的时刻。
6.设置最大驻留半径r，按时序将轨迹中的点分别设置为锚点，向后遍历，寻找连续的s个与锚点距离小于最大驻留半径的点，公式如下：
其中，为锚点，是轨迹中锚点后序的点。
7.当连续s个点与锚点构成驻留点后，将该驻留点记作，其中空间坐标为驻留点序列内的所有点的空间几何中心坐标，t1为驻留起始时间，ts为驻留终止时间，并将p
s+1
作为下一个锚点继续搜索至轨迹终点。
8.步骤s2：驻留点轨迹替换，使用识别出的sp替换其中的所有驻留点，将轨迹替换为形如的新轨迹序列表示；为轨迹第一个点，为轨迹中由p2……
ps构成的驻留点，p
s+1
为驻留点后的第一个点，为轨迹中的第二个驻留点，为轨迹的最后一个点，此处p与sp的顺序按照时间顺序排列，公式中展示的仅为一种情况的举例。
9.步骤s3：驻留点和地点的关联，查询路网信息以及地理实体数据，将驻留点与地理实体数据建立空间索引，获取驻留点相关地理实体属性，作为该轨迹的重要特征点以及事件判断依据。
10.步骤s4： 驻留点地理信息相似性度量，依据上一步骤中查询到的驻留点地理实体属性，度量不同路径之间驻留点的相同地理实体的重合度，驻留点的真实地理实体相关度，以及相同场所的轨迹经过时间之间的重合率，由此衡量轨迹驻留点地理信息相似度，构建邻接矩阵。
11.步骤s5：基于路网信息的轨迹平滑降噪，根据单向距离，计算轨迹到路网的距离，选取最近路径，将轨迹拟合到最近邻路网上。
12.步骤s6： 轨迹分段处理，依据路网信息，将轨迹按照路网关键节点（例如红绿灯、环岛、高速出入口、收费站、转弯等）拆分，并提取轨迹经过路段起止点时的时刻信息，将轨迹模型转化为时间维度模型，将轨迹的不同时间段索引到固定路段和不同驻留点即
13.其中path代表轨迹中拟合到路网后，单端的平直轨迹及其时间信息，其中pathn代表第n条轨迹，spm代表第m个驻留点；此处path与sp的顺序按照时间顺序排列，公式中展示的仅为一种情况的举例。
14.步骤s7：轨迹时空相似性度量，根据路网信息查询轨迹之间每一段的路径相关性padis，由此衡量轨迹路网路径相似度邻接矩阵；步骤s8：时空轨迹关系图构建，根据步骤4与步骤7获取的邻接矩阵，以及步骤6中对轨迹的重构，对场景中所有轨迹构建图结构，以每条轨迹作为图中的一个节点，轨迹重构作为节点特征，此处path与sp的顺序按照时间顺序排列，公式中展示的仅为一种情况的举例。节点间相似性度量，其中α和β分别为两个权重，用来衡量迹驻留点地理信息相似度邻接矩阵与轨迹路网路径相似度邻接矩阵的权重，当驻留点地理信息相似度和轨迹
路网路径相似度关注度相同时，α和β都取值0.5，从而实现联合使用轨迹时间空间信息，与路网、地理实体信息的联合轨迹相似度建模。
15.步骤s9：时空轨迹聚类，将构成的图输入无监督的图神经网络，实现轨迹间深层关联信息的挖掘，从而实现轨迹间的聚类，完成群体行为分析及团伙行为识别。
16.优选的，在步骤s4中衡量轨迹驻留点地理信息相似度具体方法为：构建邻接矩阵，其中
17.其中，tr表示轨迹，tri表示第i条轨迹，trj表示第j条轨迹，spn为第n个驻留点，spm为第m个驻留点，为spn与spm之间的时间重合度权重，dis表示地点间相关性，表示驻留点spn与spm的地理信息相似性。
18.优选的，在步骤s5中，根据路网信息，拟合待处理轨迹到最近邻路网上，其过程可以表示为
19.其中，表示路网中的第i条道路轨迹；tr表示轨迹，p为轨迹中的点。
20.优选的，在步骤s7中轨迹路网路径相似度邻接矩阵，计算方法为：
21.其中，为pathi与pathj之间的时间重合度权重，为第i条路径pathi，与第j条路径pathj的路径相似度。
22.本发明的有益效果是：本发明与现有技术相比，主要解决了当前基于轨迹的群体行为挖掘方法无法有效利用地理实体信息，只能提取空间上的轨迹重合关系等现象，通过将驻留点索引到实际地图的地理实体，增加驻留点地理实体信息关联度，引入路网信息对轨迹进行平滑降噪，通过使用路网信息替换轨迹，去除冗余信息，将杂乱的轨迹转化为固定的路网路径，同时将轨迹间的距离计算转化为路网关系查询，使用查询操作代替复杂的计算。
附图说明
23.图1是本发明总流程图。
实施方式
实施例
24.如附图所示，一种联合多元时空轨迹的群体行为模式挖掘方法，所述方法的具体步骤如下：步骤s1：驻留点检测，将轨迹建模为：
25.其中，x、y是空间坐标，t是当前点的时刻；设置最大驻留半径r，按时序将轨迹中的点分别设置为锚点，向后遍历，寻找连续的s个与锚点距离小于最大驻留半径的点，公式如下：
26.其中，为锚点，是轨迹中锚点后序的点；当连续s个点与锚点构成驻留点后，将该驻留点记作，其中空间坐标为驻留点序列内的所有点的空间几何中心坐标，t1为驻留起始时间，ts为驻留终止时间，并将p
s+1
作为下一个锚点继续搜索至轨迹终点，例如以p1为第一个锚点，设置最大驻留半径r=500m，若p1与p2间的距离dis(p1, p2)《500m，则认定p2在驻留点内，继续计算p1与p3的距离dis(p1, p3)，若dis(p1, p3)《500m, 则认定p3在驻留点内，继续计算p1与p4的距离dis(p1, p4)，以此类推，直到p1与pn的距离dis(p1, pn)《500m，且p1与p
n+1
的距离dis(p1, p
n+1
)》500m，则认定p
n+1
不在驻留点内，且{p
1 ,p
2 ,p3……
pn}构成一个驻留点；如此操作可以将目标在小范围内的复杂活动轨迹简化为一个驻留点，使用一个点代替一段轨迹，并保留时间信息，大幅度简化轨迹。
27.步骤s2：驻留点轨迹替换，使用识别出的sp替换其中的所有驻留点，将轨迹替换为形如的新轨迹序列表示；为轨迹第一个点，为轨迹中由p2……
ps构成的驻留点，p
s+1
为驻留点后的第一个点，为轨迹中的第二个驻留点，为轨迹的最后一个点，此处p与sp的顺序按照时间顺序排列，公式中展示的仅为一种情况的举例。
28.步骤s3：驻留点和地点的关联，查询路网信息以及地理实体数据，将驻留点与地理实体数据建立空间索引，获取驻留点相关地理实体属性，作为该轨迹的重要特征点以及事件判断依据，目标在某处驻留，代表该目标在驻留处有过重点活动，通过查询驻留点的地理实体数据，可以判断目标的活动性质；例如轨迹驻留点处在学校范围内，则将教育相关属性列入轨迹特征中，则目标可能身份为教师、学生、家长等的概率更大。
29.步骤s4： 驻留点地理信息相似性度量，依据上一步骤中查询到的驻留点地理实体属性，度量不同路径之间驻留点的相同地理实体的重合度，驻留点的真实地理实体相关度，以及相同场所的轨迹经过时间之间的重合率，由此衡量轨迹驻留点地理信息相似度，构建邻接矩阵，其中：
30.其中，tr表示轨迹，tri表示第i条轨迹，trj表示第j条轨迹，spn为第n个驻留点，spm为第m个驻留点，为spn与spm之间的时间重合度权重，dis表示地点间相关性，表示驻留点spn与spm的地理信息相似性；例如spn与spm均在同一处场所内，则为1，代表两处地理信息完全相同，spn在场所内停留时间为9:00~18:00，spm在场所内停留时间为10:00~19:00，则两个驻留点的时间重合度为0.8。
31.步骤s5：基于路网信息的轨迹平滑降噪，根据单向距离，计算轨迹到路网的距离，选取最近路径，将轨迹拟合到最近邻路网上，其过程可以表示为
32.其中，表示路网中的第i条道路轨迹；tr表示轨迹，p为轨迹中的点，通过计算轨迹中的一系列点到路网路径上的距离，寻找差异最小的路网路径，并使用对应的路网路径替换复杂的轨迹，保留时间信息，大幅度简化轨迹。
33.步骤s6： 轨迹分段处理，依据路网信息，将轨迹按照路网转弯等关键节点拆分，并提取轨迹经过路段起止点时的时刻信息，将轨迹模型转化为时间维度模型，将轨迹的不同时间段索引到固定路段和不同驻留点即
34.其中path代表轨迹中拟合到路网后，单端的平直轨迹及其时间信息，其中pathn代表第n条轨迹，spm代表第m个驻留点；此处path与sp的顺序按照时间顺序排列，公式中展示的仅为一种情况的举例。
35.步骤s7：轨迹时空相似性度量，根据路网信息查询轨迹之间每一段的路径相关性padis，由此衡量轨迹路网路径相似度邻接矩阵，计算方法为：
36.其中，为pathi与pathj之间的时间重合度权重，为第i条路径pathi，与第j条路径pathj的路径相似度，其中的计算类似于步骤s4中的计算，由于轨迹路网路径是由固定的路网信息组成，故其轨迹间的相似度padis可由路网信息中直接查询得出，从而简化了轨迹间相似度的计算。
37.步骤s8：时空轨迹关系图构建，根据步骤4与步骤7获取的邻接矩阵，以及步骤6中对轨迹的重构，对场景中所有轨迹构建图结构，以每条轨迹作为图中的一个节点，轨迹重构作为节点特征，此处path
与sp的顺序按照时间顺序排列，公式中展示的仅为一种情况的举例。节点间相似性度量，其中α和β分别为两个权重，用来衡量迹驻留点地理信息相似度邻接矩阵与轨迹路网路径相似度邻接矩阵的权重，当驻留点地理信息相似度和轨迹路网路径相似度关注度相同时，α和β都取值0.5，从而实现联合使用轨迹时间空间信息，与路网、地理实体信息的联合轨迹相似度建模；步骤s9：时空轨迹聚类，将构成的图输入无监督的图神经网络，实现轨迹间深层关联信息的挖掘，从而实现轨迹间的聚类，完成群体行为分析及团伙行为识别。
38.上面结合附图对本发明的具体实施方式作了详细说明，但是本发明并不限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。