HoughCircles函数

HoughCircles函数用于在灰度图像中使用霍夫变换查找圆。该函数通过修改霍夫变换来实现，通常可以很好地检测出圆的中心，但可能无法找到正确的半径。可以通过指定半径范围(minRadius和maxRadius)来协助该函数，或者在#HOUGH_GRADIENT方法中将maxRadius设置为负数以仅返回圆心而不进行半径搜索，并使用其他过程找到正确的半径。此外，还可以对图像进行一定程度的平滑处理，除非它已经很软。例如，可以使用7x7内核和1.5x1.5 sigma或类似的模糊处理来平滑图像。

函数原型：

<code>CV_EXPORTS_W 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位单通道图像，如果使用彩色图像，需要先转换成灰度图像；

circles：用来存储HoughCircles的结果，类型为list，list中对象格式为x,y,r；

method：定义检测图像中圆的方法。目前唯一实现的方法是cv::HOUGH_GRADIENT；

dp：图像像素分辨率与参数空间分辨率的比值（官方文档上写的是图像分辨率与累加器分辨率的比值，它把参数空间认为是一个累加器，毕竟里面存储的都是经过的像素点的数量），dp=1，则参数空间与图像像素空间（分辨率）一样大，dp=2，参数空间的分辨率只有像素空间的一半大；#通过设置dp可以减少计算量

minDist：检测到的圆中心（x,y）坐标之间的最小距离。如果minDist太小，则会保留大部分圆心相近的圆。如果minDist太大，则会将圆心相近的圆进行合并（若两圆心距离 < minDist，则认为是同一个圆）。

param1：canny 边缘检测的高阈值，低阈值被自动置为高阈值的一半，默认为 100。；

param2：累加平面某点是否是圆心的判定阈值。大于该阈值才判断为圆。当值设置的很小是，检测到的圆越多。默认值为 100；

minRadius：半径的最小大小（以像素为单位）默认为 0；

maxRadius：半径的最大大小（以像素为单位）默认为 0。

使用场景

HoughCircles 是一种图像处理算法，用于检测图像中的圆。霍夫圆检测能检测出目标图像中存在的圆，但在实际使用中，参数调节存在很大的困难。

参数设置不合理会不仅导致计算量非常大，并且可能无法找不到你想保留的圆。

1）边缘噪声滤除

对图像进行 中值滤波 或者 均值滤波 或者 高斯滤波 优化，减少图像中存在的边缘噪声。

2）检测参数如何设定

a. 检索半径设定

通过观察原图初步确定的圆的像素半径范围(对应minRadius和maxRadius）。

b. 圆心累加值设定（对应参数param2）

因为累加圆心(圆弧上的在圆心上的累加值，累加值超过该阈值则被认为是一个圆)

c. 圆心距离设定（对应参数miniDist）

通过观察原图所有希望找出的圆，然后确定2个圆心之间最小的距离。

使用案例

#include <iostream>

#include <opencv2/opencv.hpp>

using namespace cv;

using namespace std;

int main()

{

// 读取原始图像

Mat src = imread("3480e7ce_mk_tk_a_c2_6400_4480.png", IMREAD_COLOR);

if (src.empty())

{

cout << "无法读取图像" << endl;

return -1;

}

// 转换为灰度图像

Mat grayImage;

cvtColor(src, grayImage, COLOR_BGR2GRAY);

imshow("灰度图", grayImage);

waitKey(0);

// 使用滤波器去噪

Mat blurImage;

// GaussianBlur(grayImage, blurImage, Size(9, 9), 2, 2); // 高斯滤波

// medianBlur(grayImage, blurImage, 7); // 中值滤波

blur(grayImage, blurImage, Size(7,7)); // 均值滤波

imshow("模糊滤波", blurImage);

waitKey(0);

// 使用 Canny 边缘检测

Mat edgesImage;

Canny(blurImage, edgesImage, 50, 100);

imshow("边缘检测", edgesImage);

waitKey(0);

clock_t t1 = clock();

// 使用 HoughCircles 检测圆

vector<Vec3f> circles;

// 设置Hough变换参数

int minDist = 30; // 最小距离

int param1 = 100; // Canny边缘检测高阈值 # 保持和和前面Canny的阈值一致

int param2 = 15; // Hough变换高阈值

int minRadius = 10; // 最小半径

int maxRadius = 100; // 最大半径

HoughCircles(blurImage, circles, HOUGH_GRADIENT, 1, minDist, param1, param2, minRadius, maxRadius);

clock_t t2 = clock();

// 在原图上绘制检测到的圆

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(src, center, 3, Scalar(0, 255, 0), -1, 8, 0); // 绘制圆心

circle(src, center, radius, Scalar(0, 0, 255), 3, 8, 0); //绘制空心圆

}

clock_t t3 = clock();

double timeForHoughCircle = (double) (t2 - t1) / CLOCKS_PER_SEC * 1000;

double timeForDrawCircle = (double) (t3 - t2) / CLOCKS_PER_SEC * 1000;

cout << "HoughCircles:" << timeForHoughCircle << "ms, DrawCircle:" << timeForDrawCircle << "ms" <<endl;

// 显示结果

imshow("霍夫圆检测", src);

waitKey(0);

return 0;

}

灰度图：

 模糊滤波：

 边缘检测：

霍夫圆检测结果：

参考资料

opencv 十一 霍夫圆检测原理及高级使用案例（含优化步骤）