轨迹相似性度量方法

1.欧式距离（Euclidean Distance）

我们定义如下两条轨迹，长度分别为n和m，则：

欧式距离要求两条轨迹的长度相同一一对应，其数学定义为：

欧氏距离的定义简单明了，就是两条轨迹对应点的空间距离的平均值，但是缺点也很明显，就是不能度量不同长度的轨迹相似性，且对噪音点敏感。

2.动态时间归整（Dynamic Time Warping，DTW）

动态时间归整（Dynamic Time Warping）的思想是自动扭曲两个序列，并在时间轴上进行局部的缩放对齐，以使其形态尽可能一致，从而得到最大可能的相似性。DTW将两条轨迹的点进行多对多的映射，从而较为高效地解决了数据不齐的问题，其动态规划算法如下：


动态时间归整算法灵活，对轨迹长度无限制，且效果较好，但是其并未对噪音点进行处理，离群点也会对结果造成较大影响。
假设有两个序列P=[1,3,2,4,2,5]，Q=[0,3,4,2,2]，直接计算两者的欧氏距离为5（这里直接用差值代替平方项）
采用DWT,首先计算两个序列的距离矩阵:

从左上角开始, 可以向右,向下,或者向右下前进, 对进行到这三个方向后的距离累加和进行比较, 易知向右累加距离和为3, 向下为4, 向右下为1, 因此选择向右下前进. 状态转移方程可以表示为:

最优路径如表中红色方块所示, 最终得到的最短距离为矩阵红色块之和.

3.最长公共字串（Longest Common Sub-Sequence，LCSS）

有一个经典的算法问题：求解两序列的最长公共子序列，不要求公共子序列中的两个连续相连，例如BDCABA和ABCBDAB的最大公共子序列为BCBA。在此基础上，很自然提出了基于最长公共子序列的轨迹相似性度量方法：

其值代表最多可被视为同一点的点数，也就是两条轨迹中满足最小距离阈值限制的轨迹点的对数。
最长公共子串距离对噪音点进行了处理，即因噪音点的偏离没有与其相近的轨迹点故不会被计算在最终结果内，这一步骤有效对抗噪音。
该算法的最小距离阈值Ɛ不好定义，还有可能返回并不相似的轨迹。
两条轨迹间的相似度，有归一化:


LCSS示意图，有两对点的距离小于e，则被视为同一点，LCSS值为2

4.编辑距离（Edit Distance on Real sequence，EDR）

给定两个长度分别为n和m的轨迹tr1和tr2，最小距离的匹配阈值e，则两条轨迹之间的EDR距离就是需要对tr1及逆行插入、删除或替换使其变为tr2的操作次数，其基于动态规划的算法如下：

其中，

与LCSS相比：①LCSS只计算匹配数，而EDR计算可以使轨迹一致的总编辑操作数，更精确；②假设替换、插入、删除的成本一致，则有EDR(A,B)=L(A)+L(B)-2LCSS(A,B)

与DTW相比：EDR将距离段玲转化为0、1值，对离群点更鲁棒。

5.豪斯多夫距离（Hausdorff Distance）

简单来说，豪斯多夫距离就是两条轨迹最近点距离的最大值：

其中，h(tr1,tr2)称为tr1到tr2的单向豪斯多夫距离，其定义如下：

6.弗雷歇距离（Fréchet Distance）

将图中的两条曲线想为两条坚硬的、不会变形的钢丝，每条钢丝都串上相同数目的n个珠子，再用可任意无线自由伸缩的橡皮筋把两串珠子上对应的珠子连接起来，最后再用橡皮筋连接两条钢丝对应的端点，此刻即为初始状态。
首先，测量初始状态下各条橡皮筋的长度，并记录下它们中的最大值d0；
任意挪动两条钢丝上的若干珠子使其到达新的位置，测量各条橡皮筋的长度，并记录下它们中的最大值d1；
再重复上步（N-1）次（N足够大），取其中的最小值作为弗雷歇距离的近似值。

7.单向距离（One Way Distance，OWD）

将一条轨迹tr1视为离散点，另一条轨迹 tr2按照与离散点的最小距离相应拆分为不同轨迹段，取点和对应轨迹段形成的面积之和为距离

其中，|tr1|表示轨迹tr1的长度，d(p,tr2)表示GPS点p到tr2的距离。为了对称，简单修改上述公式：

OWD距离的基本思想基于两条轨迹围成的面积，当面积大，说明轨迹之间距离较远，相似度就低；相反，若围成的面积为0，则说明两条轨迹重合，相似度最高。

该算法一般用于轨迹搜索，如果轨迹tr1为tr2的一条子轨迹，则OWD(tr1,tr2)的值极小，接近于0；反过来OWD(tr2,tr1)的值就很大。

8.多线位置距离（Locality In-between Polylines，LIP）

取两条轨迹交叉得到的面积，用区域周长对轨迹长度的占比作为权重，计算总的面积作为距离度量：

当某区域面积的周长占总长比重大时权重也自然就大；当Area均为0时，说明两条轨迹重合没有缝隙，LIP距离为0；当Area加权和大时，则说明两条轨迹之间缝隙较大，LIP距离也就大。此外，权重由区域周长占总长比重大决定，也一定程度对抗了噪音点的干扰。