OpenCV-Python:Harris角点检测与Shi-Tomasi角点检测
一、Harris角点检测
原理:
角点特性:向任何方向移动变换都很大。
Chris_Harris 和 Mike_Stephens 早在 1988 年的文章《A CombinedCorner and Edge Detector》中就已经提出了焦点检测的方法,被称为Harris 角点检测。将窗口向各个方向移动(u,v)然后计算所有差异的总合:表达式如下:
角点检测中要使E(u,v)的值最大。这就是说必须使方程的第二项的取值最大。对上面的等式进行泰勒级数展开,然后再通过几步数学换算(参考其他标准教材),我们得到下面的等式:
其中
这里 Ix 和 Iy 是图像在 x 和 y 方向的导数(可以使用cv2.Sobel()计算得到)
根据一个用来判定窗口内是否包含角点的等式进行打分
其中
- λ1 和 λ 2 是矩阵M的特征值所以根据这些特征中我们可以判断一个区域是否是一个角点、边界或是平面
- 当λ1 和 λ 2都小时,|R| 也小,这个区域就是一个平坦区域
- 当 λ 1 ≫ λ 2 或者 λ 1 ≪ λ 2时,R小于0,这个区域是边缘
- 当 λ 1 和 λ 2 都很大,并且 λ 1 ~λ 2 中的时,R 也很大,(λ 1 和 λ 2 中最小值都大于阈值),说明这个区域是角点
用下图来表示我们的结论:
所以Harris 角点检测的结果是一个由角点分数构成的灰度图像。选取适当的阈值对结果进行二值化,我们就检测到了图像中的角点。
OpenCV中的角点检测
使用 cv2.cornerHarris(),参数如下:
img:数据类型为float32的输入图像(灰度图)
blockSize:角点检测中要考虑的领域大小
ksize:Sobel求导中使用的窗口大小
k:Harris角点检测方程中的自由参数,取值参数为[0.04,0.06]
import cv2
import numpy as np img = cv2.imread('blox.jpg') # . Harris角点检测基于灰度图
gray = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY) # . Harris角点检测
dst = cv2.cornerHarris(gray, , , 0.04)
# 腐蚀一下,便于标记
dst = cv2.dilate(dst, None) # 角点标记为红色
img[dst > 0.01 * dst.max()] = [, , ]
cv2.imwrite('blox-RedPoint.png', img)
cv2.imshow('dst', img)
cv2.waitKey()
二、Shi-Tomasi角点检测 & 适合于跟踪的图像特征
Shi-Tomasi 角点检测是在Harris角点检测上进行了修改。我们知道Harris角点检测的打分公式为:
但 Shi-Tomasi 使用的打分函数为:
如果打分查过阈值,我们就认为它是一个角点,我们可以把它绘制到 λ 1 ~λ 2 空间中,就会得到下图:
从这幅图中,我们可以看出只有当λ 1 和 λ 2 都大于最小值时,才被认为是角点(绿色区域)
OpenCV使用cv2.goodFeatureToTrack()。这个函数可以帮我们使用Shi-Tomasi方法获取图像中的N个最好的角点。
image : 输入灰度图
maxCorners:想要检测的角点数目
qualityLevel:设置角点的最低质量,0 到 1之间。低于这个数的所有角点都会被忽略
minDistance:两个角点之间的最短欧式距离
根据这些信息,函数就能在图像上找到角点。所有低于质量水平的角点都会被忽略。然后再把合格角点按角点质量进行降序排列。函数会采用角点质量最高的那个角点(排序后的第一个),然后将它附近(最小距离之内)的角点都删掉。按着这样的方式最后返回 N 个最佳角点。
import numpy as np
import cv2
from matplotlib import pyplot as plt img = cv2.imread('blox.jpg')
gray = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY) # Shi-Tomasi角点检测
corners = cv2.goodFeaturesToTrack(gray, , 0.01, )
corners = np.int0(corners) # 个角点坐标 for i in corners:
# 压缩至一维:[[, ]] -> [, ]
x, y = i.ravel()
cv2.circle(img, (x, y), , (, , ), -) cv2.imwrite('Shi-Tomasi-corner.jpg', img)
cv2.imshow('dst', img)
cv2.waitKey()
这个函数很适合在目标跟踪中使用
三、Harris和Shi-Tomasi的比较
import numpy as np
import cv2
from matplotlib import pyplot as plt img = cv2.imread('blox.jpg')
imgShi, imgHarris = np.copy(img), np.copy(img)
gray = cv2.cvtColor(img, cv2.COLOR_BGR2GRAY) # . Shi-Tomasi角点检测
# 优点:速度相比Harris有所提升,可以直接得到角点坐标
corners = cv2.goodFeaturesToTrack(gray, , 0.01, )
corners = np.int0(corners) for i in corners:
# 压缩至一维:[[, ]] -> [, ]
x, y = i.ravel()
cv2.circle(imgShi, (x, y), , (, , ), -) # . Harris角点检测
dst = cv2.cornerHarris(gray, , , 0.04)
# 腐蚀一下,便于标记
dst = cv2.dilate(dst, None)
# 角点标记为红色
imgHarris[dst > 0.01 * dst.max()] = [, , ] cv2.imwrite('compare.png', np.hstack((imgHarris, imgShi)))
cv2.imshow('compare', np.hstack((imgHarris, imgShi)))
cv2.waitKey()
左图为Harris角点检测结果,右图为Shi-Tomasi结果。可以看出,Shi-Tomasi角点检测的质量更高,数量也相对较少。一般情况下,Shi-Tomasi相比Harris速度有所提升,并且可以直接得到角点坐标。
OpenCV-Python:Harris角点检测与Shi-Tomasi角点检测的更多相关文章
- OpenCV + python 实现人脸检测(基于照片和视频进行检测)
OpenCV + python 实现人脸检测(基于照片和视频进行检测) Haar-like 通俗的来讲,就是作为人脸特征即可. Haar特征值反映了图像的灰度变化情况.例如:脸部的一些特征能由矩形特征 ...
- OpenCV Python教程(3、直方图的计算与显示)
转载请详细注明原作者及出处,谢谢! 本篇文章介绍如何用OpenCV Python来计算直方图,并简略介绍用NumPy和Matplotlib计算和绘制直方图 直方图的背景知识.用途什么的就直接略过去了. ...
- opencv+python视频实时质心显示
利用opencv+python实现以下功能: 1)获取实时视频,分解帧频: 2)将视频做二值化处理: 3) 将视频做滤波处理(去除噪点,获取准确轮廓个数): 4)识别图像轮廓: 5)计算质心: 6)描 ...
- Erlang/Elixir: 使用 OpenCV, Python 搭建图片缩略图服务器
这篇文章是在OSX上测试和运行的的, Ubuntu下的安装和配置请移步到这里 应用程序进程树, 默认 Poolboy 中初始化10个用于处理图片的 Python 工作进程(Worker) 首先安装Op ...
- opencv python实用操作
画多边形 fillConvexPloy与fillConvexPloy的区别 fillConvexPloy 用来画单个凸多边形: 如果点集的连线不是凹多边形,则会找一个最小的凸多边形把该凹多边形包住画出 ...
- Win7 64bit+Anaconda(3-5.0.1,Python3.6)+Pycharm(community-2017.3.3)+OpenCV(python‑3.4.0‑cp36‑cp36m)(转载)
Anaconda(3-5.0.1,Python3.6)下载链接:https://pan.baidu.com/s/1bqFwLMB 密码:37ih Pycharm(community-2017.3.3) ...
- 解决Windows Server 2012 R2 Datacenter云服务器无法运行opencv python程序的问题
写了个基于opencv的python程序,pyinstaller 32位机打包后在win7/win10 32/64正常运行,在Windows Server 2012 R2 Datacenter云服务器 ...
- OpenCV Python 4.0安装
1.安装OpenCV-Python 如果你是第一次使用OpenCV Python开发包,想要安装OpenCV Python4.0只要执行如下命令行即可: pip install opencv-pyth ...
- Ubuntu系统---安装Caffe (+OpenCV+Python+CPU-only)
安装配置Ubuntu14.04+Caffe (+OpenCV+Python+CPU-only) 记录 [作者:Wu Ping.时间:20180428.] 本人已经安装很多次的Caffe了:从开始的初探 ...
- 使用Python基于VGG/CTPN/CRNN的自然场景文字方向检测/区域检测/不定长OCR识别
GitHub:https://github.com/pengcao/chinese_ocr https://github.com/xiaofengShi/CHINESE-OCR |-angle 基于V ...
随机推荐
- jatoolsprinter html实现每隔几秒获取数据直接后台打印不弹窗
1.流程说明 jatoolspringter 必须要能在html代码里面看到 id =page1 page2 page3..... 才能打印,所以无法动态打印,必须先把要打印的内容放到页面某个地方隐 ...
- Python_001_开始学习的一些准备
从Hello World!开始,希望你能坚持一下,然后每一天有一点收获. ♥ 安装了pip和setuotools库: pip用于安装管理Python扩展包,setuptools用于分布Python ...
- [PL]如果天空是黑暗的,那就摸黑生存
“如果天空是黑暗的,那就摸黑生存:如果发出声音是危险的,那就保持沉默:如果自觉无力发光的,那就蜷缩于墙角.但不要习惯了黑暗,就为黑暗辩护:不要为自己的苟且而得意:不要嘲讽那些比自己更勇敢更热情的人们. ...
- Linux-Jenkins安装部署
Jenkins 安装及插件安装 Jenkins简介: Jenkins只是一个平台,真正运作的都是插件.这就是jenkins流行的原因,因为jenkins什么插件都有 Hudson是Jenkins的前身 ...
- CC++语法::数组名退化(array decaying)
参考: CSDN::C/C++中数组名退化为指针的情况 stackoverflow::What is array decaying? 起因 笔者在写memset的时候总想偷一点懒(因为我们一般都是为了 ...
- elasticsearch中head插件中的定制增加用户名密码范例
在head插件目录下一般 在 elasticsearch目录下的 plugins\head目录 下 在 或 plugins\head\site目录下 有 一个index.html文件.把这个文件用下面 ...
- 《Linux下cp XXX1 XXX2的功能》的实现
<Linux下cp XXX1 XXX2的功能>的实现 一.题目要求 编写MyCP.java 实现类似Linux下cp XXX1 XXX2的功能,要求MyCP支持两个参数: java MyC ...
- java第一课 面向对象的编程概念
一.什么是对象(object)? 对象是相关状态和行为的软件包. 1.现实社会的对象都有两个共同特征:状态和行为.如:狗有状态(名称,颜色,品种,饥饿)和行为(吠叫,取出,摇尾). 2.软件对象在概念 ...
- pt-show-grants的用法
pt-show-grants的用法 1.先查找所有用户和Host 2.然后逐个执行show grants pt-show-grants的功能是格式化打印输出MySQL上的赋权,以便你可以有效地复制.比 ...
- (二) V4L2引入(含浅析UVC)
title: V4L2引入(含浅析UVC) date: 2019/4/23 19:00:00 toc: true --- V4L2引入(含浅析UVC) 基本框架 V4L2全名是video for li ...