OpenCV 文档：霍夫圆变换

参考：https://blog.csdn.net/dcrmg/article/details/52506538

霍夫圆变换的基本思路是认为图像上每一个非零像素点都有可能是一个潜在的圆上的一点，跟霍夫线变换一样，也是通过投票，生成累积坐标平面，设置一个累积权重来定位圆。

在笛卡尔坐标系中圆的方程为：

           

其中（a,b）是圆心，r是半径，也可以表述为：

即

所以在abr组成的三维坐标系中，一个点可以唯一确定一个圆。

而在笛卡尔的xy坐标系中经过某一点的所有圆映射到abr坐标系中就是一条三维的曲线：

经过xy坐标系中所有的非零像素点的所有圆就构成了abr坐标系中很多条三维的曲线。

在xy坐标系中同一个圆上的所有点的圆方程是一样的，它们映射到abr坐标系中的是同一个点，所以在abr坐标系中该点就应该有圆的总像素N0个曲线相交。通过判断abr中每一点的相交（累积）数量，大于一定阈值的点就认为是圆。

以上是标准霍夫圆变换实现算法，问题是它的累加面试一个三维的空间，意味着比霍夫线变换需要更多的计算消耗。Opencv霍夫圆变换对标准霍夫圆变换做了运算上的优化。它采用的是“霍夫梯度法”。它的检测思路是去遍历累加所有非零点对应的圆心，对圆心进行考量。如何定位圆心呢？圆心一定是在圆上的每个点的模向量上，即在垂直于该点并且经过该点的切线的垂直线上，这些圆上的模向量的交点就是圆心。

霍夫梯度法就是要去查找这些圆心，根据该“圆心”上模向量相交数量的多少，根据阈值进行最终的判断。

函数原型：

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是输入图像矩阵，要求是灰度图像；转换类型为：COLOR_BGR2GRAY

第二个参数 circles是一个包含检测到的圆的信息的向量，向量内第一个元素是圆的横坐标，第二个是纵坐标，第三个是半径大小；

第三个参数 methodmethod是所使用的圆检测算法，目前只有CV_HOUGH_GRADIENT一个可选；

第四个参数 dp是累加面与原始图像相比的分辨率的反比参数，dp=2时累计面分辨率是元素图像的一半，宽高都缩减为原来的一半，dp=1时，两者相同。（关于这个分辨率的概念没有理解透，按道理低分辨率应该意味着更快的检测速度，然而实测恰恰相反）

第五个参数 minDist定义了两个圆心之间的最小距离；

第六个参数param1是Canny边缘检测的高阈值，低阈值被自动置为高阈值的一半；

第七个参数param2是累加平面对是否是圆的判定阈值；

第八和第九个参数定义了检测到的圆的半径的最大值和最小值；
 

注意：

因为霍夫圆检测对噪声比较敏感，所以使用霍夫圆变换时首先要对图像做滤波：中值滤波、高斯滤波等。

示例程序：

int main(int argc, char** argv) {
    Mat src, dst;
    src = imread(STRPAHT3);
    if (!src.data) {
        printf("could not load image...\n");
        return -1;
    }
    
    Mat moutput;
    //中值滤波medianBlur
    medianBlur(src, moutput, 3);
    cvtColor(moutput, moutput, CV_BGR2GRAY);

    //基于灰度空间
    vector<Vec3f> pcircles;
    HoughCircles(moutput, pcircles, CV_HOUGH_GRADIENT, 1, 10, 100, 30, 5, 50);
    src.copyTo(dst);
    //重新传回rgb色彩空间
    cvtColor(dst, dst, CV_BGR2GRAY);
    for (size_t i = 0; i < pcircles.size(); i++) {
        Vec3f cc = pcircles[i];
        circle(dst, Point(cc[0], cc[1]), cc[2], Scalar(0, 0, 255), 2, LINE_AA);
        circle(dst, Point(cc[0], cc[1]), 2, Scalar(198, 23, 155), 2, LINE_AA);
    }
    imshow("OUTPUT_TITLE", dst);

    waitKey(0);
}