利用霍夫梯度法进行圆检测的原理概要及利用OpenCV的函数HoughCircles()实现霍夫梯度法圆检测的代码_霍夫梯度法检测圆的原理-程序员宅基地

技术标签: 霍夫梯度法  HoughCircles  图像处理原理、工具、代码  

霍夫圆变换的基本原理与霍夫线变换原理类似。利用霍夫变换检测直线的相关知识可以参考我的另一篇博文,链接 https://blog.csdn.net/wenhao_ir/article/details/51774444

对直线检测而言,一条直线可由极坐标参数r和θ确定,对于圆来说,则需要三个参数来确定一个圆(为什么是三个,圆心坐标需要两个参数,还要加上一个半径)。标准霍夫圆变换就依据这个思路将直角坐标转换到描述圆的三维度空间中,然后使用这三个维度进行累加度量(投票),根据投票的结果判断是否为圆。

由于是三个参数,所以相当于转换到了三维空间,这样带来的问题就是计算量是几何级增大的,所以我们不用这个思路来实现圆检测法。为了克服这个问题,便有了霍夫梯度法来检测圆的算法。

霍夫梯度法的思路:第一步根据每个点的模向量来找到圆心,这样三维的累加平面就转化为二维累加平面;第二步根据所有候选中心的边缘非零像素对其的支持程度来确定半径。

如果想了解霍夫梯度法的算法原理及其优点和缺点,可以参考下面这两篇文章:

http://t.zoukankan.com/jsxyhelu-p-7192711.html

https://blog.csdn.net/DIPDWC/article/details/117513601

由于上面的第一篇文章并不是CSDN的,容易出现网页丢失的情况,所以博主把网页保存到了本地,大家如果打不开了,可以下载下来浏览,百度网盘下载链接如下:https://pan.baidu.com/s/1qSXnIkLGq3kSsePebqnWUw?pwd=t0w6 

OpenCV中提供了函数HoughCircles()来实现霍夫梯度法,其函数原型如下

void HoughCircles(	InputArray image, 
					OutputArray circles,
					int method, 
					double dp, 
					double minDist,
					double param1 = 100, 
					double param2 = 100,
					int minRadius = 0, 
					int maxRadius = 0 );

官方文档对参数的解释如下:

image – 8-bit, single-channel, grayscale input image. 
circles – Output vector of found circles. Each vector is encoded as a 3-element floating-point vector (x, y, radius) . 
method – Detection method to use. Currently, the only implemented method is CV_HOUGH_GRADIENT , which is basically 21HT , described in [Yuen90]. 
dp – Inverse ratio of the accumulator resolution to the image resolution. For example, if dp=1 , the accumulator has the same resolution as the input image. If dp=2 , the accumulator has half as big width and height. 
minDist – Minimum distance between the centers of the detected circles. If the parameter is too small, multiple neighbor circles may be falsely detected in addition to a true one. If it is too large, some circles may be missed. 
param1 – First method-specific parameter. In case of CV_HOUGH_GRADIENT , it is the higher threshold of the two passed to the Canny() edge detector (the lower one is twice smaller). 
param2 – Second method-specific parameter. In case of CV_HOUGH_GRADIENT , it is the accumulator threshold for the circle centers at the detection stage. The smaller it is, the more false circles may be detected. Circles, corresponding to the larger accumulator values, will be returned first. 
minRadius – Minimum circle radius. 
maxRadius – Maximum circle radius. 

我的翻译如下

image-源图像,要求是8 位、单通道、灰度输入图像。

circles-检测到的圆的三个参数就存在这个动态数组中,数组的每个成员包含三个参数,分别为圆心的x坐标、圆心的y坐标和圆的半径。

method-检测方法选择参数,目前只支持 CV_HOUGH_GRADIENT这种方法,即霍夫梯度法。

dp-累加器图像分辨率的反比。要理解这个,必须举例子,比如若dp==1,那么累加器的分辨率与原图像一样,如果dp==2,那么累加器的分辨率只有原图像的一半。可见,这个值越大,累加器的分辨率反而越低,所以这是个反比值。

minDist-检测到圆的圆心之间的最小距离,显然,这个值越小,伪圆可能越多,而这个值越大,则有越多的圆被漏检。

param1-函数HoughCircles是用到了Canny作边缘检测的,而Canny算法中的滞后阈值法要求设定高低两个阈值,这个参数就是设定这个高阈值的,而低阈值为这个高阈值的二分之一。

param2 -圆心阈值参数。圆心是通过投票得出的,那么多少票才算是圆心呢?这个值就是确定这个问题的。

minRadius -检测到的圆的最小半径

maxRadius-检测到的圆的最大半径

接下来上示例代码:

代码中用到的图像的下载链接:https://pan.baidu.com/s/1bo6hEft

//博主微信/QQ 2487872782
//有问题可以联系博主交流
//有图像处理开发需求请联系博主
//图像处理技术交流QQ群 271891601

//OpenCV版本:3.0
//VS版本:2013

#include <opencv2/core/core.hpp>
#include <opencv2/highgui/highgui.hpp>
#include<opencv2/imgcodecs/imgcodecs.hpp>
#include <opencv2/imgproc/imgproc.hpp>


using namespace cv;
using namespace std;

int main(int argc, char** argv)
{
	cv::Mat srcImage = imread("F:/material/images/P0043-circle.jpg");
	if (!srcImage.data)
		return -1;
	cv::imshow("srcImage", srcImage);
	// 转换为灰度图像
	cv::Mat src_gray;
	cvtColor(srcImage, src_gray, CV_BGR2GRAY);
	// 高斯平滑滤波
	GaussianBlur(src_gray, src_gray, Size(9, 9), 2, 2);
	vector<Vec3f> circles;
	// 霍夫圆检测
	HoughCircles(src_gray, circles, CV_HOUGH_GRADIENT,
		1, src_gray.rows / 8, 200, 100, 0, 0);
	// 将得到的结果绘图
	for (size_t i = 0; i < circles.size(); i++)
	{
		Point center(cvRound(circles[i][0]), cvRound(circles[i][1]));
		int radius = cvRound(circles[i][2]);
		// 绘制圆中心
		circle(srcImage, center, 3, Scalar(0, 255, 0), -1, 8, 0);
		// 绘制圆轮廓
		circle(srcImage, center, radius, Scalar(120, 120, 120), 3, 8, 0);
	}
	cv::imshow("HoughResult", srcImage);
	cv::waitKey(0);
	return 0;
}

运行结果如下图所示

注意:利用OpenCV的函数HoughCircles()实现霍夫梯度法圆检测时参数dp的值对于最终结果影响是挺大的,详情见博文 https://blog.csdn.net/wenhao_ir/article/details/125140993

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/wenhao_ir/article/details/51783566

智能推荐

class和struct的区别-程序员宅基地

文章浏览阅读101次。4.class可以有⽆参的构造函数,struct不可以,必须是有参的构造函数,⽽且在有参的构造函数必须初始。2.Struct适⽤于作为经常使⽤的⼀些数据组合成的新类型,表示诸如点、矩形等主要⽤来存储数据的轻量。1.Class⽐较适合⼤的和复杂的数据,表现抽象和多级别的对象层次时。2.class允许继承、被继承,struct不允许,只能继承接⼝。3.Struct有性能优势,Class有⾯向对象的扩展优势。3.class可以初始化变量,struct不可以。1.class是引⽤类型,struct是值类型。

android使用json后闪退,应用闪退问题:从json信息的解析开始就会闪退-程序员宅基地

文章浏览阅读586次。想实现的功能是点击顶部按钮之后按关键字进行搜索,已经可以从服务器收到反馈的json信息,但从json信息的解析开始就会闪退,加载listview也不知道行不行public abstract class loadlistview{public ListView plv;public String js;public int listlength;public int listvisit;public..._rton转json为什么会闪退

如何使用wordnet词典,得到英文句子的同义句_get_synonyms wordnet-程序员宅基地

文章浏览阅读219次。如何使用wordnet词典,得到英文句子的同义句_get_synonyms wordnet

系统项目报表导出功能开发_积木报表 多线程-程序员宅基地

文章浏览阅读521次。系统项目报表导出 导出任务队列表 + 定时扫描 + 多线程_积木报表 多线程

ajax 如何从服务器上获取数据?_ajax 获取http数据-程序员宅基地

文章浏览阅读1.1k次,点赞9次,收藏9次。使用AJAX技术的好处之一是它能够提供更好的用户体验,因为它允许在不重新加载整个页面的情况下更新网页的某一部分。另外,AJAX还使得开发人员能够创建更复杂、更动态的Web应用程序,因为它们可以在后台与服务器进行通信,而不需要打断用户的浏览体验。在Web开发中,AJAX(Asynchronous JavaScript and XML)是一种常用的技术,用于在不重新加载整个页面的情况下,从服务器获取数据并更新网页的某一部分。使用AJAX,你可以创建异步请求,从而提供更快的响应和更好的用户体验。_ajax 获取http数据

Linux图形终端与字符终端-程序员宅基地

文章浏览阅读2.8k次。登录退出、修改密码、关机重启_字符终端

随便推点

Python与Arduino绘制超声波雷达扫描_超声波扫描建模 python库-程序员宅基地

文章浏览阅读3.8k次,点赞3次,收藏51次。前段时间看到一位发烧友制作的超声波雷达扫描神器,用到了Arduino和Processing,可惜啊,我不会Processing更看不懂人家的程序,咋办呢?嘿嘿,所以我就换了个思路解决,因为我会一点Python啊,那就动手吧!在做这个案例之前先要搞明白一个问题:怎么将Arduino通过超声波检测到的距离反馈到Python端?这个嘛,我首先想到了串行通信接口。没错!就是串口。只要Arduino将数据发送给COM口,然后Python能从COM口读取到这个数据就可以啦!我先写了一个测试程序试了一下,OK!搞定_超声波扫描建模 python库

凯撒加密方法介绍及实例说明-程序员宅基地

文章浏览阅读4.2k次。端—端加密指信息由发送端自动加密,并且由TCP/IP进行数据包封装,然后作为不可阅读和不可识别的数据穿过互联网,当这些信息到达目的地,将被自动重组、解密,而成为可读的数据。不可逆加密算法的特征是加密过程中不需要使用密钥,输入明文后由系统直接经过加密算法处理成密文,这种加密后的数据是无法被解密的,只有重新输入明文,并再次经过同样不可逆的加密算法处理,得到相同的加密密文并被系统重新识别后,才能真正解密。2.使用时,加密者查找明文字母表中需要加密的消息中的每一个字母所在位置,并且写下密文字母表中对应的字母。_凯撒加密

工控协议--cip--协议解析基本记录_cip协议embedded_service_error-程序员宅基地

文章浏览阅读5.7k次。CIP报文解析常用到的几个字段:普通类型服务类型:[0x00], CIP对象:[0x02 Message Router], ioi segments:[XX]PCCC(带cmd和func)服务类型:[0x00], CIP对象:[0x02 Message Router], cmd:[0x101], fnc:[0x101]..._cip协议embedded_service_error

如何在vs2019及以后版本(如vs2022)上添加 添加ActiveX控件中的MFC类_vs添加mfc库-程序员宅基地

文章浏览阅读2.4k次,点赞9次,收藏13次。有时候我们在MFC项目开发过程中,需要用到一些微软已经提供的功能,如VC++使用EXCEL功能,这时候我们就能直接通过VS2019到如EXCEL.EXE方式,生成对应的OLE头文件,然后直接使用功能,那么,我们上篇文章中介绍了vs2017及以前的版本如何来添加。但由于微软某些方面考虑,这种方式已被放弃。从上图中可以看出,这一功能,在从vs2017版本15.9开始,后续版本已经删除了此功能。那么我们如果仍需要此功能,我们如何在新版本中添加呢。_vs添加mfc库

frame_size (1536) was not respected for a non-last frame_frame_size (1024) was not respected for a non-last-程序员宅基地

文章浏览阅读785次。用ac3编码,执行编码函数时报错入如下:[ac3 @ 0x7fed7800f200] frame_size (1536) was not respected for anon-last frame (avcodec_encode_audio2)用ac3编码时每次送入编码器的音频采样数应该是1536个采样,不然就会报上述错误。这个数字并非刻意固定,而是跟ac3内部的编码算法原理相关。全网找不到,国内音视频之路还有很长的路,音视频人一起加油吧~......_frame_size (1024) was not respected for a non-last frame

Android移动应用开发入门_在安卓移动应用开发中要在活动类文件中声迷你一个复选框变量-程序员宅基地

文章浏览阅读230次,点赞2次,收藏2次。创建Android应用程序一个项目里面可以有很多模块,而每一个模块就对应了一个应用程序。项目结构介绍_在安卓移动应用开发中要在活动类文件中声迷你一个复选框变量

推荐文章

热门文章

相关标签