【OpenCV文档】用于角点检测的Fast算法

来源：互联网发布：深圳关键词排名软件编辑：程序博客网时间：2024/04/27 20:41

原文地址：http://docs.opencv.org/trunk/doc/py_tutorials/py_feature2d/py_fast/py_fast.html#fast-algorithm-for-corner-detection

目标

理解FAST算法的基本原理
使用OpenCV的FAST函数进行角点（corners）检测

原理

我们已知很多种特征检测的方法，而且它们其中很多效果都非常不错。但是，当从一个实时运行的程序角度出发，它们还不够快。一个最好的例子就是SLAM（Simultaneous Localization and Mapping）移动机器人，它的可计算资源是受限的。

作为上述问题的一个解决方法，FAST（Features from Accelerated Segment Test）算法由Edward Rosten 和 Tom Drummond在他们2006年的论文“Machine learning for high-speed corner detection”（在2010年再次被修订）中被提出。下面是该算法的一个基本总结。你可以参见原论文去学习更多的细节（下面所有的图片是从原论文中提取的）。

使用FAST进行特征检测

从图片中选取一个像素点 $p$ ，下面我们将判断它是否是一个特征点。我们首先把它的密度（即像素值）设为 $I_p$ 。
设定一个合适的阙值 $t$ 。
考虑该像素点周围的16个像素。（见下图）
现在，如果在这个大小为16个像素的圆上有 $n$ 个连续的像素点，它们的像素值要么都比 $I_p + t$ 大，要么都比 $I_p − t$ 小，那么它就是一个角点。（如上图中白色虚线所示）。 $n$ 这里被设定为12。
我们现在应该提出一个高效的测试，来快速排除一大部分是非角点的点。该测试仅仅检查在位置1、9、5和13四个位置的像素（首先检查1和9，看它们是否亮于或暗于阙值。如果是，再检查5和13）。如果 $p$ 是一个角点，那么上述四个像素点中至少有3个应该必须都大于 $I_p + t$ 或者小于 $I_p − t$ （因为若是一个角点，超过四分之三个圆的部分应该满足判断条件，半圆比包含上述某三个点）。如果都不满足，那么 $p$ 不可能是一个角点。完整的分段测试可以被用于接受的所有候选点，通过检测圆上的所有点。这种检测有很好的性能，但是有一些缺点：
1. 当n < 12时不能拒绝许多候选点。
2. 检测出来的角点不是最优的，这是因为它的效率是依靠角点外形的排列和分布的。
3.
4. 相邻的多个特征点会被检测到。

前三个问题可以使用机器学习的方法解决。最后一个可以使用non-maximal suppression。

机器学习

选择一个图片集合进行学习（最好是来自于目标应用定义域）。
在每一张图上运行FAST算法，找到特征点。
对于每个特征点，存储它周围的16个像素点到一个vector中。为所有的图片做同样的事，得到它们所有特征点的vector $P$ 。
这16个像素中的每一个像素（假设为 $x$ ），可以有下面三种状态中的一种：
依靠这些状态，特征向量 $P$ 被划分为3个子集， $P_d$ , $P_s$ , $P_b$ 。
定义一个新的布尔变量， $K_p$ 。如果 $p$ 是一个角点，那些 $K_p$ 为真；否则为假。
使用ID3算法（决策树分类器）来查询每一个子集。
递归计算所有的子集直到它的嫡为0。
被创建的决策树就被用于其他图片的fast检测。

non-maximal suppression

从邻近的位置选取了多个特征点是另一个问题。我们可以使用non-maximal suppression来解决。

为每一个检测到的特征点计算它的分数函数（score function）， $V$ 。 $V$ 是 $p$ 和它周围16个像素点的绝对偏差的和。
考虑两个相邻的特征点，并比较它们的 $V$ 值。
$V$ 值较低的点将会被剔除。

总结

FAST算法比其他已知的角点检测法要快很多倍。

但是当图片的噪点较多时，它的健壮性并不好。这依靠一个阙值。

OpenCV中的FAST特征检测

和OpenCV中其他特征检测器的调用相同。如果你想，你可以指定一个阙值，或者决定是否使用non-maximal suppression来判断邻近特征点。

对于相邻特征点，OpenCV定义了三个flags：

cv2.FAST_FEATURE_DETECTOR_TYPE_5_8
cv2.FAST_FEATURE_DETECTOR_TYPE_7_12
cv2.FAST_FEATURE_DETECTOR_TYPE_9_16

下面是一个简单的示例代码。

import numpy as npimport cv2from matplotlib import pyplot as pltimg = cv2.imread('simple.jpg',0)# Initiate FAST object with default valuesfast = cv2.FastFeatureDetector()# find and draw the keypointskp = fast.detect(img,None)img2 = cv2.drawKeypoints(img, kp, color=(255,0,0))# Print all default paramsprint "Threshold: ", fast.getInt('threshold')print "nonmaxSuppression: ", fast.getBool('nonmaxSuppression')print "neighborhood: ", fast.getInt('type')print "Total Keypoints with nonmaxSuppression: ", len(kp)cv2.imwrite('fast_true.png',img2)# Disable nonmaxSuppressionfast.setBool('nonmaxSuppression',0)kp = fast.detect(img,None)print "Total Keypoints without nonmaxSuppression: ", len(kp)img3 = cv2.drawKeypoints(img, kp, color=(255,0,0))cv2.imwrite('fast_false.png',img3)

下面是测试结果。左边是使用了nonmaxSuppression的FAST，右边则没有使用。

0 0