【学习笔记】Fréchet距离的 C 语言实现

Fréchet 距离

Fréchet 距离，又称为弗雷歇距离，是一种衡量两条曲线（或两个路径）之间相似性的度量方法。这个概念最初在度量空间理论中被定义，后来被广泛应用于计算机科学、地理信息系统、图像处理、生物信息学等多个领域，特别是在比较和匹配形状、路径或序列时非常有用。

数学背景

首先，我们需要一个基础的数学框架。假设有两个参数化的曲线 $A (t)$ 和 $B (t)$ ，它们定义在一个闭区间 $[0, 1]$ 上，并且映射到一个度量空间 $(X, d)$ 中，其中 $d$ 是这个空间内的两点间距离函数。度量空间可以简单理解为一个集合 $X$ 加上一个定义在这个集合上任意两点间距离的方式 $d$ 。

Fréchet距离的定义

Fréchet 距离试图量化这样的场景：假设有两个人分别沿着曲线 $A$ 和 $B$ 行走，他们可以从各自的起点出发，以任意的速度前进，但要求两人始终保持同步，即每个人在自己路径上的位置与另一人在其路径上的位置之间有一个对应关系。Fréchet 距离是保持彼此距离尽可能小的情况下所需的最大距离。

Fréchet距离的数学定义可以用以下公式表达：

设 $\rightarrow X$ 和 $\rightarrow X$ 是定义在度量空间 $(X, d)$ 中的两条连续曲线，其中 $d$ 是 $X$ 中两点间的距离函数。Fréchet距离 $F (A, B)$ 定义为：

$\inf_{\alpha, \beta} \max_{t \in [0, 1]} d(A(\alpha(t)), B(\beta(t)))$

这里， $\alpha$ 和 $\beta$ 是从区间 $[0, 1]$ 到自身的连续单调递增函数，代表了沿曲线 $A$ 和 $B$ 的“行走”速度。 $\alpha(t)$ 和 $\beta(t)$ 分别表示在时间 $t$ 时，沿曲线 $A$ 和 $B$ 的位置。 $\inf$ 表示 infimum（下确界），即所有可能的 $\alpha$ 和 $\beta$ 对中最大的 $d(A(\alpha(t)), B(\beta(t)))$ 的最小值。

简单来说，这个公式表达了找到两个函数（代表两条路径）上的点对，这些点对在它们各自的路径上以某种方式“同步行走”，使得在任何给定时间点 $t$ 上，这两点之间的最大距离尽可能小，最终得到的就是Fréchet距离。

在离散情况下，如果我们将曲线 $A$ 和 $B$ 近似为点集 ${A_1, A_2, ..., A_n\}$ 和 ${B_1, B_2, ..., B_m\}$ ，则可以通过动态规划等方法来近似计算Fréchet距离，此时的公式和计算方法会有所不同，但基本思想仍然遵循上述原理。

C语言实现

以下是计算 Fréchet 距离的C语言实现示例：

#include <stdio.h>
#include <stdlib.h>
#include <math.h>// 定义点结构体
typedef struct {double x;double y;
} Point;// 计算两点之间的欧几里得距离
double euclidean_distance(Point p1, Point p2) {return sqrt(pow(p1.x - p2.x, 2) + pow(p1.y - p2.y, 2));
}// 递归计算Fréchet距离
double recursive_frechet_distance(Point* P, int i, Point* Q, int j, double** cache) {if (cache[i][j] != -1) {return cache[i][j];}double dist = euclidean_distance(P[i], Q[j]);if (i == 0 && j == 0) {cache[i][j] = dist;} else if (i == 0) {cache[i][j] = fmax(recursive_frechet_distance(P, i, Q, j - 1, cache), dist);} else if (j == 0) {cache[i][j] = fmax(recursive_frechet_distance(P, i - 1, Q, j, cache), dist);} else {double minPrevDist = fmin(fmin(recursive_frechet_distance(P, i - 1, Q, j, cache),recursive_frechet_distance(P, i - 1, Q, j - 1, cache)),recursive_frechet_distance(P, i, Q, j - 1, cache));cache[i][j] = fmax(minPrevDist, dist);}return cache[i][j];
}// 计算Fréchet距离
double frechet_distance(Point* P, int n, Point* Q, int m) {// 创建缓存数组并初始化double** cache = (double**)malloc(n * sizeof(double*));for (int i = 0; i < n; i++) {cache[i] = (double*)malloc(m * sizeof(double));for (int j = 0; j < m; j++) {cache[i][j] = -1;}}double result = recursive_frechet_distance(P, n - 1, Q, m - 1, cache);// 释放缓存数组for (int i = 0; i < n; i++) {free(cache[i]);}free(cache);return result;
}int main() {// 定义两条曲线Point P[] = {{0, 0}, {1, 1}, {2, 2}};Point Q[] = {{0, 0}, {1, 2}, {3, 3}};// 计算Fréchet距离double distance = frechet_distance(P, 3, Q, 3);printf("Fréchet Distance: %f\n", distance);return 0;
}