traj_dist 笔记 源代码解析（python部分）

1distance.py

1.1 METRIC_DIC

不同实现方法对应的函数路径

1.2 sspd

功能：

• 计算轨迹 traj_1 和 traj_2 之间的对称化段路径距离。

参数：

• traj_1：一个二维 numpy 数组，代表第一个轨迹。
• traj_2：一个二维 numpy 数组，代表第二个轨迹。
• type_d：字符串，指定距离的类型（"euclidean" 或 "geographical"）。
  • "euclidean" 表示使用欧几里得距离，适用于任何维度的轨迹。
  • "geographical" 表示使用地理距离，假设轨迹是二维的，第一维是经度，第二维是纬度。
• implementation：字符串，指定实现方式（"python"、"cython" 或 "auto"）。
  • 如果轨迹是二维的，默认使用 Cython 实现；否则使用 Python 实现，除非明确指定为 "python"。

返回值：

• 返回一个浮点数，表示两个轨迹之间的对称化段路径距离。

1.3  sowd_grid

两个轨迹之间的对称化单向距离（Symmetrized One-Way Distance）

大体逻辑和1.2 相同，这边多了两个内容

• converted：布尔值，如果为 True，则表示轨迹已经转换为网格表示形式
• precision：Geohash 精度，用于地理坐标的转换。
• 'sowd_grid' 计算网格表示中轨迹之间的距离。
  • 如果坐标是地理的，可以根据 geohash 编码进行转换。如果是这样，需要 geohash 的“精度”。
  • 否则，轨迹被认为是以单元格表示编码的。

1.4 frechet

1.5 discret_frechet

1.6 hausdorff

1.7 dtw

逻辑和sspd一样（1.2节），略过

1.8 lcss

1.9 edr

逻辑和sspd一样（1.2节），只是多了一个eps（阈值），略过

1.10 erp

逻辑和sspd一样（1.2节），只是多了一个g（比较基准点），略过

1.11 pdist

用于计算轨迹列表（traj_list）中每对轨迹之间的成对距离

1.12 cdist 

用于计算两个轨迹列表（traj_list_1 和 traj_list_2）中每对轨迹之间的距离

实现逻辑和pdist几乎一模一样，唯一的区别就是一个返回的是一维向量，一个是二维数组

2 pdist/basic_euclidean.py 欧氏距离轨迹相关基本操作

2.1 eucl_dist(x,y) 两个坐标点的差的二范数

简单的两个坐标点的差的二范数

2.2 point_to_seg 点到线段的距离

2.3 point_to_trajectory  一个点与一条轨迹的最短距离

2.4  circle_line_intersection 圆和线段的交点

 2.4.1 二次方程的推导

2.4.2  推导坐标

3 basic_geographical.py  地理空间相关的操作

3.1 great_circle_distance 两个经纬度的大圆距离

3.2 initial_bearing 一个点到另一个点的初始方位角 

3.3 cross_track_distance  点到大圆路径的距离

3.4 along_track_distance  从轨迹上点到不在轨迹上点的沿轨迹距离

3.5 point_to_path  一个点  到一条由两点 和  定义的路径的点到路径距离

4 pdist/sspd.py 一条轨迹上的点到零一条轨迹的距离

4.1 e_spd 轨迹 t1 上的所有点到轨迹 t2 的   点到轨迹的距离 

4.2 e_sspd 两个轨迹之间相互的 spd-距离的平均值

4.3   g_spd 轨迹 t1 上所有点到轨迹 t2 的  平均点到路径距离（？）

这里代码的注释是轨迹 t1 上所有点到轨迹 t2 的  平均点到路径距离，我觉得应该是t2到t1？

4.4 g_sspd 两个轨迹之间相互的 spd-距离的和

注：和欧几里得的版本不同，几何空间版本是求和，而不是取平均

5 pdist/lcss.py 最长公共子序列

5.1 e_lcss

大体上和文巾解题1143. 最长公共子序列-CSDN博客是一样的，仅有一点点区别

5.2 g_lcss

基本上都是一样的，把计算欧几里得距离换成了计算大圆距离 

6 pdist/hausdorff.py 豪斯多夫距离

6.1 e_directed_hausdorff t1到t2的豪斯多夫距离

 6.2 e_hausdorff 双向豪斯多夫

t1到t2和t2到t1的单向豪斯多夫距离中较大的那个

6.3 g_directed_hausdorff

6.4 g_hausdorff

和欧氏距离的唯一区别是，point_to_trajectory改成了point_to_path

7  pdist/erp.py

地理坐标类的就是将欧氏距离改成大圆距离

8 pdist/edr.py

8 pdist/dtw.py

 

7数学笔记/scipy 笔记：豪斯多夫距离（Hausdorff ）_python 豪斯多夫距离-CSDN博客