[发明专利]一种通信辐射源相似运动轨迹比较方法

权利要求书

1.一种通信辐射源相似运动轨迹比较方法，其特征在于，包括以下步骤：

获取目标通信辐射源的查询轨迹和待比较轨迹；所述查询轨迹为已知的历史轨迹数据；

提取查询轨迹和待比较轨迹的时空信息，计算查询轨迹和待比较轨迹的时空相似度；所述轨迹的时空信息包括经纬度和采样时间；

提取查询轨迹和待比较轨迹的语义标签，对查询轨迹和待比较轨迹的语义标签进行匹配，并计算查询轨迹和待比较轨迹的语义相似度；所述轨迹的语义标签包括信号特征参数和地理热点地区；

设置两个加权值，分别乘以时空相似度和语义相似度，并求和得到综合相似度；

将综合相似度与预先设定的阈值相比较，若大于阈值，则判定所述待比较轨迹判断为查询轨迹的相似轨迹，输出待比较轨迹。

2.如权利要求1所述的一种通信辐射源相似运动轨迹比较方法，其特征在于，所述时空相似度的计算过程为：

将轨迹A映射到三维空间中，那么该轨迹表示为Traj＝{X₁,X₂,...,X_N}，其中X_i＝(x_i,y_i,t_i)，i＝1,2,...,N为轨迹序列上每一个点的向量表示，x_i,y_i,t_i分别表示经度值、纬度值和采样时间，假设X₁＝(x₁,y₁,t₁)和X₂＝(x₂,y₂,t₂)是Traj中任意两个不同向量，那么X₁和X₂之间的马氏距离为：

D_Mah(X₁,X₂)＝(X₁-X₂)WC^-1W^Τ(X₁-X₂)^T，

其中C为X₁,X₂之间的协方差矩阵，f_x,f_y,f_t分别为各元素的权重值，为常数值；

假设有查询轨迹Traj_A＝{X₁,X₂,...,X_m,...,X_M}和待比较轨迹Traj_B＝{Y₁,Y₂,...,Y_n,...,Y_N}，那么Traj_A和Traj_B之间的距离为：

其中d为轨迹数据点维度，表示方差，X_m分成d维表示为Y_n分成d维表示为

将Traj_A和Traj_B按时间顺序分割为k个子轨迹对，表示为(Traj_A1,Traj_B1),...,(Traj_Ak,Traj_Bk)，分别计算这k个待分析子轨迹对的距离得到对每一个子轨迹对构造相似子轨迹标记向量index＝[index(1),index(2),...,index(k)]，其中

其中th_dis为判定子轨迹对(Traj_Ai,Traj_Bi)相似的门限值；

计算得到Traj_A和Traj_B的时空相似度Sim_s-t为

其中，l＝[l₁,l₂,...,l_k]，其中l_i为每一个Traj_Bi的长度，即Traj_Bi中的数据点数量乘以采样间隔，i＝1,2,...,k，N为Traj_B的数据点数量。