利用霍夫梯度法进行圆检测的原理概要及利用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

智能推荐

如何确定IP在同一网段_判断ip是否在一个网段-程序员宅基地

文章浏览阅读2.6w次,点赞20次,收藏84次。一 。IP地址分类我们将IP地址分为了【A B C D E】五大类,但是D和E不对民用组织开放,所以我们能使用IP地址的只有ABC三类。判断四组IP的第一个数组来确定是哪类。【注意:127.0.0.1不再里面,是因为它代表当前计算机自己】二。判断每一类IP是否属于一个网段(一)A类1. A类IP用第一个数字来表示不同网段。例如:1.0.0.0和2.0.0.0是不同的网络..._判断ip是否在一个网段

html最基础的几个标签_html中最基本的标签aheadbhtmlctitle-程序员宅基地

文章浏览阅读313次。基本网页标签将一个网页看成人的话,那么一个网页标签也分别有表示皮肤、头、脸、身体的部位<html> <head> <!--head中一般写入一些浏览器中的配置标签--> <meta charset='utf-8' /> <title> 这是我的第一个网页</title> </head> &..._html中最基本的标签aheadbhtmlctitle

电子类经典电子书分享(信号系统 模电 ADI)_放大器设计实践125问全解-程序员宅基地

文章浏览阅读668次。最近阿里冲容量,拿些压箱底的东西出来吧。没有注册的,可以用我的链接注册一下,再用手机app登一下,咱们各享500G空间。我在使用不限速「阿里云盘」,赠送你 500GB 快来试试吧------------我在使用不限速「阿里云盘」,赠送你 500GB 快来试试吧 ------------点此链接领取福利:https://pages.aliyundrive.com/mobile-page/web/beinvited.html?code=21ed59d1 信号与系统包括以下内._放大器设计实践125问全解

Windows域控禁用U盘和光盘的策略 【全域策略生效】_域控关闭自动打开u盘策略-程序员宅基地

文章浏览阅读2.6k次,点赞5次,收藏2次。目录开始--管理工具--组策略管理进入组策略界面在主域下编辑DefaultDomainPolicy按照以下路径打开策略编辑器可移动存储访问U盘策略和光盘策略分开特殊情况测试域策略立即生效客户端生效生效效果开始--管理工具--组策略管理进入组策略界面在主域下编辑DefaultDomainPolicy按照以下路径打开策略编辑器可移动存储访问可移动光盘:控制U盘的读取、执行、写入权限CD和DVD:控制光盘的读取、执行、写入权._域控关闭自动打开u盘策略

通信原理板块——数字数据压缩编码之霍夫曼编码_霍夫曼压缩编码-程序员宅基地

文章浏览阅读206次。信源字符表可以是二进制的,也可使多字符的,非二进制字符可以通过一个字符编码表映射为二进制码字;③x7和x8合并后的复合字符的概率为P(x7)+P(x8)=1/8,并将新得到一组字符按照概率不增大的次序排列,注意:新复合字符与x3和x4概率相同,可放置在x2和x5之间的任何位置,此例子放置在x4之后,替换x5。等长码中表示每个字符的码字长度是相同的,但是各字符所含有的信息量是不同的。等长码可以通过计数的方法确定字符的分界,但变长码则不可以,接收端收到一长串变长码,不一定能确定每个字符的分界。_霍夫曼压缩编码

centos/ubuntu—yum/apt-get软件安装_yum安装apt-get-程序员宅基地

文章浏览阅读6.7k次。这里注意,info后面必须是软件完整名称,例如查询mysql-connector-java.noarch,名称应该是mysql-connector-java,写mysql会报错。或者写作mysql*,查询所有含有mysql名称的软件包参考centos-vault | 镜像站使用帮助 | 清华大学开源软件镜像站 | Tsinghua Open Source Mirror鸟哥私房菜 - 第二十二章、软件安装 RPM, SRPM 与 YUM (vbird.org)/var/l_yum安装apt-get

随便推点

linux编程常用指令_linux实现代码的指令-程序员宅基地

文章浏览阅读154次。一、网络函数(1)htons(2)inet_addr 点格式转换为无符号长整型 ina.sin_addr.s_addr = inet_addr("132.241.5.10");(3)inet_ntoa 无符号长整型转换为点格式 struct in_addr addr1; ulong l1; l1= inet_addr("192.168.0.74"..._linux实现代码的指令

memmove函数_void *memmove(void *dest, const void *src, size_t -程序员宅基地

文章浏览阅读1.1k次。在这个例子中,字符串"hello world!"被移动,使得输出为 “hello hello world!是一个标准库函数,用于 C++ 中的内存操作。它主要用于在内存中移动或复制字节。可以处理源内存区和目标内存区重叠的情况。如果源内存区和目标内存区重叠,此函数会返回指向目标内存区的指针。仍然可以正确地复制字节,而。_void *memmove(void *dest, const void *src, size_t n);

2013最新版Subversion 1.7.10 for Windows x86 + Apache 2.4.4 x64 安装配置教程+错误解决方案...-程序员宅基地

文章浏览阅读140次。一 、工作环境  操作系统:Windows Server 2008 R2 SP1 x64  Apache版本:2.4.4  Subversion版本: Setup-Subversion-1.7.10.msi  TortoiseSVN版本:TortoiseSVN-1.7.13.24257-x64-svn-1.7.10.msi +LanguagePack_1.7.13.24257..._windows subversion x64

python拟合直线的斜率_线性曲线拟合总是得到斜率和y的中间值为1-程序员宅基地

文章浏览阅读1.1k次。我得到了一些东西作为Excel的线性拟合,使用scipy basinhopping而不是曲线拟合和大量迭代。运行迭代需要一点时间,而且还需要一个错误函数,但它是在没有缩放原始数据的情况下完成的。Basinhopping docs.import numpy as npimport matplotlib.pyplot as pltfrom scipy.optimize import basinhopp..._一阶线性拟合斜率

广西计算机应用基础试题及答案,修订版成人大专考试试卷(有答案).-程序员宅基地

文章浏览阅读376次。广西医科大学融水函授站成人大专考试试卷课程:计算机应用基础专业:班级:姓名:学号:一、Word与Excel部分1.编辑排版一个文档完毕后,若要知道其打印效果,可选择( A )功能。A.打印预览 B.模拟打印 C.提前打印 D。屏幕打印2.Word窗口中打开文档MWA,修改后另存为MWC,则文档( B )。A.MWA是当前文档 B.MWC是当前文档C.MWC和MWA都是当前文档 D.MWC和MWA均...

T-sql语句修改数据库逻辑名、数据库名、物理名-程序员宅基地

文章浏览阅读104次。--更改MSSQL数据库物理文件名Sql语句的写法--注意:要在活动监视器里面确保没有进程连接你要改名的数据库!!!!!!!!!!!!!!!!!!!!-- Sql语句如下USE master --改逻辑名ALTER DATABASE YQBlog MODIFY FILE(NAME='YQBlogAA',NEWNAME='YQBlog') -- GOALTER..._tdsql修改数据库名称

推荐文章

热门文章

相关标签