[图像识别]11.OpenCV的角点检测 FAST和OBR算法

回顾，上节课你学到了什么？

#实例化sift

sift=cv.SIFT_create()  

#检测关键点并计算

kp,des=sift.detectAndCompute(gray,None)

#将关键点绘制在图像上

cv.drawKeypoints(image,keypoints,outputimage,color,flags)

目录

FAST

1.原理

2.步骤

3.代码API

 OBR

1.原理

2.步骤

3.代码API

 总结，这节课你学到了什么？

---

FAST

1.原理

前面讲述的SIFT和SURF算法效率还是太低了，为解决这个问题提出了FAST算法

FAST算法通俗的将就是检测一个像素周围有一定数量的像素与该点像素值不同，就认为它是角点。

2.步骤

【1】在图像中选取一个像素点p,来判断它是不是关键点。Ip等于像素点p的灰度值。
【2】以r为半径画圆，覆盖p点周围的M个像素，通常情狂下，设置r=3,则M=16

【3】设置一个阈值t,如果在这16个像素点中存在n个连续像素点的灰度值都在（Ip- t，Ip+ t）之外，那么像素点p就被认为是一个角点。如上图中的虚线所示，n一般取值为12。
【4】为了快速筛选，首先对候选点的周围每个90度的点:上述图片的 1, 9, 5,13进行测试(先测试1和19,如果它们符合阈值要求再测试5和13)。如果p是角点，那么这四个点中至少有3个要符合阈值要求，否则直接剔除。对保留下来的点再继续进行测试(是否有12的点符合阈值要求)。

缺点:●获得的候选点比较多●特征点的选取不是最优的●大量的点被丢弃●很多特征点都是相邻的
前3个问题可以通过机器学习的方法解决，最后一个问题可以使用非最大值抑制的方法解决。

3.代码API

#实例化fast

fast=cv.FastFeatureDetector_create(threshold,nonMaxSuppression)

threshold:阈值默认10

nonMaxSuppression：是否进行非极大抑制

#检测关键点，没有关键点描述

kp=fast.detect(gray,None)

kp:关键点信息，位置，尺度，方向

gray：灰度图像

#将观测点绘制到图像上

cv.drawKeyPoints(img,keypoints,outputimg,color,flags)

举例：

#10.1

import cv2 as cv
import numpy as np
import matplotlib.pyplot as plt
#解决中文显示问题,固定格式，直接复制下面俩行代码就行
plt.rcParams['font.sans-serif']=['SimHei']
plt.rcParams['axes.unicode_minus']=False

#1.读取灰度图像
pic1 =cv.imread("10.jpg")
gray=cv.cvtColor(pic1,cv.COLOR_BGR2GRAY)
#2.FAST实例化
fast=cv.FastFeatureDetector_create(15,True)
#3.检测关键点
kp=fast.detect(gray,None)
#4.绘制关键点
pic2=pic1.copy()
cv.drawKeypoints(pic2,kp,pic2,(0,0,255),flags=4)
#绘制图像
fig,axes=plt.subplots(nrows=2,ncols=2,figsize=(10,10))
axes[0,0].set_title("原图")
axes[0,0].imshow(pic1[:,:,::-1])
axes[0,1].set_title("灰度图")
axes[0,1].imshow(gray,plt.cm.gray)
axes[1,1].set_title("FAST特征点检测")
axes[1,1].imshow(pic2[:,:,::-1])

plt.show()
cv.waitKey(0)

 

 结果：

 OBR

1.原理

ORB算法结合了Fast和Brief算法，提出了构造金字塔，为Fast特征点添加了 方向，从而使得关键点具有了尺度不变性和旋转不变性。

2.步骤

【1】构造尺度金字塔，金字塔共有n层，每一层仅有一副幅图像。第0层为初始图像σ0尺寸默认1.2，第s层的尺度为:

 第s层图像大小

【2】在不同的尺度上利用Fast算法检测特征点，采用Harris角点响应函数，根据角点响应值R排序，选取前N个特征点，作为本尺度的特征点。
【3】计算特征点的主方向，计算以特征点为圆心半径为r的圆形邻域内的灰度质心位置，将从特征点位置到质心位置的方向做特征点的主方向。

质心位置

主方向位置

【4】为了解决旋转不变性，将特征点的邻域旋转到主方向.上利用Brief算法构建特征描述符，至此就得到了ORB的特征描述向量。

Brief算法：BRIEF是一种特征描述子提取算法，生成二值化描述子的算法，不提取代价低，匹配只需要使用简单的汉明距离利用比特之间的异或操作就可以完成。因此，时间代价低，空间代价低,效果还挺好。

（1）图像滤波，去除噪声

（2）选取点对:以特征点为中心,取S"S的邻域窗口，在窗口内随机选取N组点对，一般N=128,256,512,默认是256，有多种形式

（3）构建描述符，根据公式二值化：

 p(x)和p(y)是点对的俩个端点的像素值。

形成BRIEF的关键 点的描述特征向量，该向量一般为128-512位的字符串。

3.代码API

#实例化orb

orb=cv.ORB_creat(nfeatures)

nfeatures:特征点最大数量

#检测关键点并计算

kp,des=orb.detectAndCompute(gray,None)

kp:关键点集合

des:关键点描述符

gray:灰度图

#关键点检测结果绘制到图像上面

cv.drawKeypoints(img,keypoints,outputimg,flags)

举例：

#10.2

import cv2 as cv
import numpy as np
import matplotlib.pyplot as plt
#解决中文显示问题,固定格式，直接复制下面俩行代码就行
plt.rcParams['font.sans-serif']=['SimHei']
plt.rcParams['axes.unicode_minus']=False

#1.读取灰度图像
pic1 =cv.imread("10.jpg")
gray=cv.cvtColor(pic1,cv.COLOR_BGR2GRAY)
#2.FAST实例化
orb=cv.ORB_create(10000)
#3.检测关键点
kp,des=orb.detectAndCompute(gray,None)
#4.绘制关键点
pic2=pic1.copy()
cv.drawKeypoints(pic2,kp,pic2,(0,0,255),flags=0)
#绘制图像
fig,axes=plt.subplots(nrows=2,ncols=2,figsize=(10,10))
axes[0,0].set_title("原图")
axes[0,0].imshow(pic1[:,:,::-1])
axes[0,1].set_title("灰度图")
axes[0,1].imshow(gray,plt.cm.gray)
axes[1,1].set_title("ORB特征点检测")
axes[1,1].imshow(pic2[:,:,::-1])

plt.show()
cv.waitKey(0)

结果：

 总结，这节课你学到了什么？

#实例化fast

fast=cv.FastFeatureDetector_create(threshold,nonMaxSuppression)

#检测关键点，没有关键点描述

kp=fast.detect(gray,None)

#关键点检测结果绘制到图像上面

cv.drawKeypoints(img,keypoints,outputimg,flags)

#实例化orb

orb=cv.ORB_creat(nfeatures)

#检测关键点并计算

kp,des=orb.detectAndCompute(gray,None)

#关键点检测结果绘制到图像上面

cv.drawKeypoints(img,keypoints,outputimg,flags)