【学习笔记】Hausdorff 距离的 C 语言实现

dadalaohua

于 2024-06-19 22:58:16 发布

阅读量875

点赞数 23

分类专栏：学习笔记文章标签：学习笔记 c语言人工智能程序员

本文链接：https://blog.csdn.net/u012028275/article/details/139814895

版权

学习笔记专栏收录该内容

40 篇文章 11 订阅

订阅专栏

Hausdorff 距离的 C 语言实现

简介

Hausdorff 距离，亦称为豪斯多夫距离，是由德国数学家费利克斯·豪斯多夫（Felix Hausdorff）提出的一种度量空间中两个非空集合之间远近程度的方法。这种距离度量广泛应用于计算机视觉、图像处理、形状分析、医学影像分析等领域，因为它能有效衡量两个复杂结构或形状的相似性，尤其是在存在噪声、尺度变化或轻微形变的情况下。

基本定义

给定两个集合 $A$ 和 $B$ ，它们可以是点集、几何形状的边界点集合等，在一个度量空间中（通常假设为欧几里得空间），Hausdorff距离 $H (A, B)$ 定义为两个集合间最远的最近点距离的最大值。具体来说，它可以通过两个单向Hausdorff 距离来描述：

从集合 $A$ 到集合 $B$ 的单向Hausdorff距离 $h (A, B)$ 定义为集合 $A$ 中每个点到集合 $B$ 中最近点的距离中的最大值：
$\max_{a \in A} \min_{b \in B} ||a - b||$
从集合 $B$ 到集合 $A$ 的单向Hausdorff距离 $h (B, A)$ 同理定义为集合 $B$ 中每个点到集合 $A$ 中最近点的距离中的最大值：
$\max_{b \in B} \min_{a \in A} ||b - a||$

3.双向Hausdorff距离 $H (A, B)$ 取这两个单向距离中的最大值：
$H(A, B) = \max\{h(A, B), h(B, A)\}$

这表明 Hausdorff 距离衡量的是一个集合中的任一点到另一个集合中最近点的最大距离。如果这个距离很小，说明两个集合在某种程度上是“接近”的；反之，如果这个距离很大，则表明它们之间差异显著。

特性

非对称性：虽然通常考虑的是双向 Hausdorff 距离以获得对称性质，但单向Hausdorff距离本身是非对称的，即 $\neq h(B, A)$ 可能成立。
鲁棒性：Hausdorff 距离对小的形状变化（如噪声）相对不敏感，但能捕捉到大的结构差异，适合于形状匹配和识别任务。
计算复杂性：计算Hausdorff距离可能相当复杂，特别是对于大规模数据集，因为需要计算所有点对之间的距离并找出最大值。
应用：在图像处理中，Hausdorff距离常用于比较边界轮廓、分割结果的评估、医学图像分析等。

Hausdorff距离因其强大的表达能力和实用性，在多个领域内成为一种重要的形状分析工具。

C语言实现

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

#include <stdio.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));
}

// 计算集合P到集合Q的Hausdorff距离
double hausdorff_distance(Point* P, int n, Point* Q, int m)
{
    double maxDistPQ = 0.0;
    double maxDistQP = 0.0;

    // 计算从P到Q的Hausdorff距离
    for (int i = 0; i < n; i++) {
        double minDist = INFINITY;
        for (int j = 0; j < m; j++) {
            double dist = euclidean_distance(P[i], Q[j]);
            if (dist < minDist) {
                minDist = dist;
            }
        }
        if (minDist > maxDistPQ) {
            maxDistPQ = minDist;
        }
    }

    // 计算从Q到P的Hausdorff距离
    for (int i = 0; i < m; i++) {
        double minDist = INFINITY;
        for (int j = 0; j < n; j++) {
            double dist = euclidean_distance(Q[i], P[j]);
            if (dist < minDist) {
                minDist = dist;
            }
        }
        if (minDist > maxDistQP) {
            maxDistQP = minDist;
        }
    }

    // 返回Hausdorff距离
    return fmax(maxDistPQ, maxDistQP);
}

int main()
{
    // 定义两个点集
    Point P[] = {{0, 0}, {1, 1}, {2, 2}};
    Point Q[] = {{0, 0}, {1, 2}, {3, 3}};

    // 计算Hausdorff距离
    double distance = hausdorff_distance(P, 3, Q, 3);
    printf("Hausdorff Distance: %f\n", distance);

    return 0;
}