Harris角点检测及代码分析（续）

最新推荐文章于 2024-11-05 09:30:01 发布

mqfcu8

最新推荐文章于 2024-11-05 09:30:01 发布

阅读量742

点赞数 1

分类专栏： opencv 图像处理

本文链接：https://blog.csdn.net/mqfcu8/article/details/45971661

版权

opencv 同时被 2 个专栏收录

2 篇文章 0 订阅

订阅专栏

图像处理

2 篇文章 0 订阅

订阅专栏

本文详细介绍了OpenCV中的cvGoodFeaturesToTrack()函数，该函数用于角点检测，并对比了不同版本间的优化差异。从Harris角点检测到非极大值抑制，再到距离滤波器的应用，文章全面解析了角点检测的算法步骤。

摘要由CSDN通过智能技术生成

紧接着“Harris角点检测及代码分析”，这里主要分析OpenCV的cvGoodFeaturesToTrack()函数，这才是角点提取的真正代码。
Jianbo Shi, Carlo Tomasi. Good Features to Track. CVPR94
时间：2015-05-25 16:12

算法步骤

cvCornerHarris()函数仅计算出每个像素单元的 $R$ 值，然后再对 $R$ 进行过滤，仅留下 $\ge \beta \times max{R}$ 的像素点；
非极大值抑制（Non-Maximum Suppression），我看到网上很多人直接Dilate-Image && Origin-Image，但这样会形成角点块，除非有下面的步骤；
对该阶段的角点 $R$ 值进行降序排序；
距离滤波器，过滤点靠太近的点，留下 $R$ 值高的角点。

代码分析

便于理解先看OpenCV1.0的cvGoodFeaturesToTrack()源码，已删除了与理解算法无关的代码。

void cvGoodFeaturesToTrack(...)
{
    // 是否使用harris角点响应器
    if(use_harris)
    {
        // harris角点响应值，det(M)-a*(trace(M))^2
        cvCornerHarris(img, eig, block_size, 3, harris_k);
    }
    else
    {
        // 用最小特征值作为响应值
        cvCornerMinEigenVal(img, eig, block_size, 3);
    }
    cvMinMaxLoc(eig, 0, &max_val, 0, 0, mask);
    cvThreshold(eig, eig, max_val * 0.01, 0, CV_THRESH_TOZERO);
    cvDilate(eig, tmp);

    min_dist = cvRound(min_distance * min_distance);

    // 检查原响应图与dilate图的每个像素，如果不为0且相等，则为候选角点
    for(y = 1, k = 0; y < size.height - 1; y++）
    {
        for(x = 1; x < size.width - 1; x++)
        {
            if(eig_data[x] != 0 && eig_data[x] == tmp_data[x])
                ptr_data[k++] = eig_data + x;
        }
    }

    // 对候选角点按响应值降序排序
    icvSortFeatures( ptr_data, k, 0 );

    // 遍历所有候选角点，如果当前角点与前面所有角点的欧式距离小于min_dist，
    // 则标记为真正角点
    for(i = 0; i < k; i++)
    {
        for(j = 0; j < count; j++)
        {
            int dx = x - ptr[j].x;
            int dy = y - ptr[j].y;
            int dist = dx * dx + dy * dy;

            if(dist < min_dist)
                break;
        }

        if(j == count)
        {
            ptr[count].x = x;
            ptr[count].y = y;
            if(++count >= max_count)
                break;
        }
    }
}