图像局部描述符-Harris角点检测

Harris角点检测

图像特征类型

图像特征类型可以分为三种：

1、边缘

2、角点（感兴趣关键点）

3、斑点（感兴趣区域）

角点是个很特殊的存在。它们在图像中可以轻易定位，同时它们在人造物体中场景中，如门窗，桌子等中随处可见。因为角点位于两条边缘的交点处，代表了两个边缘变化的方向上的点，所以它们是可以精确定位的二维特征，甚至可以达到亚像素的精度。且其图像梯度有很高的变化，这种变化是可以用来帮助检测角点的。需要注意的是，角点与位于相同强度区域上的点不同，与物体轮廓上的点也不同，因为轮廓点难以在相同的其他物体上精确定位

角点的定义

‘如果某一点在任意方向的一个微小变动都会引起灰度很大的变化，那么我们就把它称之为角点’
关于角点的具体描述可以有几种：
1、一阶导数（即灰度的梯度）的局部最大所对应的像素点
2、两条及两条以上边缘的点
3、图像中的梯度值和梯度方向的变化速率都很高的点
4、角点处的一阶导数最大，二阶导数为0，指示物体边缘变化不连续的方向

角点特征

轮廓之间的交点。
对于同一场景，即使视角发生变化，通常具备稳定性质的特征。
该点附近区域的像素点无论在梯度方向上还是其梯度幅值上有着较大变化。
角点处的一阶导数最大，二阶导数为0。
角点指示了物体边缘变化不连续的方向。

Harris角点算法

Harris角点检测是特征点检测的基础，提出了应用邻近像素点灰度差值概念，从而进行判断是否为角点、边缘、平滑区域。Harris角点检测原理是利用移动的窗口在图像中计算灰度变化值，其中关键流程包括转化为灰度图像、计算差分图像、高斯平滑、计算局部极值、确认角点。

算法思想

利用矩形窗在图像上移动，若窗内包含有角点，则窗口向各个方向移动时，窗内的灰度值都会发生变化。从而达到检测图像角点的目的。如果像素周围显示存在多于一个方向的边，我们认为该点为兴趣点。该点就称为角点。

harris角点检测实例代码

# -*- coding: utf-8 -*-
from pylab import *
from PIL import Image
from PCV.localdescriptors import harris

"""
Example of detecting Harris corner points (Figure 2-1 in the book).
"""

# 读入图像
im = array(Image.open('../data/empire.jpg').convert('L'))

# 检测harris角点
harrisim = harris.compute_harris_response(im)

# Harris响应函数
harrisim1 = 255 - harrisim

figure()
gray()

#画出Harris响应图
subplot(141)
imshow(harrisim1)
print harrisim1.shape
axis('off')
axis('equal')

threshold = [0.01, 0.05, 0.1]
for i, thres in enumerate(threshold):
    filtered_coords = harris.get_harris_points(harrisim, 6, thres)
    subplot(1, 4, i+2)
    imshow(im)
    print im.shape
    plot([p[1] for p in filtered_coords], [p[0] for p in filtered_coords], '*')
    axis('off')

#原书采用的PCV中PCV harris模块
#harris.plot_harris_points(im, filtered_coords)

# plot only 200 strongest
# harris.plot_harris_points(im, filtered_coords[:200])

show()

测试场景1-奖杯（边缘丰富区域）

原图

正面图片

侧面图片

稍近距离图片

照明情况下图片

threshold = [0.01, 0.05, 0.1]

随着 α的值增大，角点检测的数量随之减少，也就是角点的检测灵敏度降低了。

同一事物的不同角度也会影响角点的检测，比如侧面的照片在α值为0.1的情况下仍然在奖杯的底部检测到了大量的角点，相较于正面图片的检测数量增加了很多，我觉得这是因为侧面图片让本来从正面看不明显的边缘显得更加立体明显，从而导致即使α的值很大，仍然能够检测到大量的角点。

照明条件更好的一组中，角点检测数量并没有和没有照明的一组出现太大的差别，这是因为Harris角点检测时，使用了微分算子对图像进行微分运算，而微分运算对图像密度的拉升或收缩和对亮度的抬高或下降不敏感。换言之，对亮度和对比度的仿射变换并不改变Harris响应的极值点出现的位置，但是，由于阈值的选择，可能会一定程度影响角点检测的数量