OpenCV图像处理中“找圆技术”的使用

一、为什么“找圆”

     圆是基本图形的一种，更为重要的是，自然情况下采集的图像，很少大量存在“圆”；但凡存在的，大都是人工的，那么就必然代表特定的意义，从而方便定位、分割和识别。

    OpenCV现有代码中能够直接“找圆”，主要有2个，一个是“HoughCircle ”，另一个是“BlobDetector ”，此外基本的轮廓分析也能够用于圆的寻找。但是这些基础的方法，涉及到的参数比较多，一方面我们需要深入理解、一方面需要融合运用，才能够有效提高识别准确率。因此结合实践，整理相关内容如下：

1. “找圆”在图像处理中的价值和应用案例；

2. 深入理解“HoughCircle ”的参数设置和优缺点；

3. 深入理解 “BlobDetector”的参数设置和应用实践；

4. 进一步理解”阈值-轮廓-分割“的分割方法和在“找圆”上的运用；

5. 融合目前技术，提出”找圆算法链“，提高识别准确率。

6. 对圆度、凸性、惯性比等基础知识的进一步认识。

希望能够为图像处理工程师、爱好者提供一些启发。

二、有效“找圆”的方法

OpenCV现有代码中，设计“找圆”算法的，主要有2个，一个是“HoughCircle ”，另一 个是“ BlobDetector   ”，此外基本的轮廓分析也能够用于圆的寻找。

 2.1HoughCircle   霍夫圆变换

  

代码： 

#include <opencv2/imgproc.hpp>
#include <opencv2/highgui.hpp>
#include <math.h>
using namespace cv;
using namespace std;
int main(int argc, char** argv)
{
    Mat img, gray;
    if( argc != 2 || !(img=imread(argv[1], 1)).data)
        return -1;
    cvtColor(img, gray, COLOR_BGR2GRAY);
    // smooth it, otherwise a lot of false circles may be detected
    GaussianBlur( gray, gray, Size(9, 9), 2, 2 );
    vector<Vec3f> circles;
    HoughCircles(gray, circles, HOUGH_GRADIENT,2, gray.rows/4, 200, 100 );
    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]);
         // draw the circle center
         circle( img, center, 3, Scalar(0,255,0), -1, 8, 0 );
         // draw the circle outline
         circle( img, center, radius, Scalar(0,0,255), 3, 8, 0 );
    }
    namedWindow( "circles", 1 );
    imshow( "circles", img );
    waitKey(0);
    return 0;
}

特别需要注意的是，目前版本出现新参数“HOUGH_GRADIENT_ALT”，在 默认参数下比以前有很大程度的精度提升：

但是HoughCircle的缺点也是显而易见的，简单来说，在默认参数下，它非常容易丢目标。

2.2BlobDetector 

所谓Blob就是图像中一组具有某些共同属性（例如，灰度值）的连接像素。OpenCV提供了一种方便的方法来检测斑点并根据不同的特征对其进行过滤。

// Setup SimpleBlobDetector parameters.
SimpleBlobDetector : :Params params;
// Change thresholds
params.minThreshold  =  10;
params.maxThreshold  =  200;
// Filter by Area.
params.filterByArea  =  true;
params.minArea  =  1500;
// Filter by Circularity
params.filterByCircularity  =  true;
params.minCircularity  =  0. 1;
// Filter by Convexity
params.filterByConvexity  =  true;
params.minConvexity  =  0. 87;
// Filter by Inertia
params.filterByInertia  =  true;
params.minInertiaRatio  =  0. 01;
# if CV_MAJOR_VERSION  <  3    // If you are using OpenCV 2
   // Set up detector with params
  SimpleBlobDetector detector(params);
   // You can use the detector this way
   // detector.detect( im, keypoints);
# else
 

        // Set up detector with params

        cv::Ptr<cv::SimpleBlobDetector> detector = cv::SimpleBlobDetector::create(params);

        vector<KeyPoint> keypoints;

        detector->detect(screw1, keypoints);  
# endif

在OpenCV中实现的叫做SimpleBlobDetector，它基于以下描述的相当简单的算法，并且进一步由参数控制，具有以下步骤。

SimpleBlobDetector : :Params : :Params()
{
    thresholdStep  =  10;     //二值化的阈值步长，即公式1的t
    minThreshold  =  50;    //二值化的起始阈值，即公式1的T1
    maxThreshold  =  220;     //二值化的终止阈值，即公式1的T2
     //重复的最小次数，只有属于灰度图像斑点的那些二值图像斑点数量大于该值时，该灰度图像斑点才被认为是特征点
    minRepeatability  =  2;   
     //最小的斑点距离，不同二值图像的斑点间距离小于该值时，被认为是同一个位置的斑点，否则是不同位置上的斑点
    minDistBetweenBlobs  =  10;
 
    filterByColor  =  true;     //斑点颜色的限制变量
    blobColor  =  0;     //表示只提取黑色斑点；如果该变量为255，表示只提取白色斑点
 
    filterByArea  =  true;     //斑点面积的限制变量
    minArea  =  25;     //斑点的最小面积
    maxArea  =  5000;     //斑点的最大面积
 
    filterByCircularity  =  false;     //斑点圆度的限制变量，默认是不限制
    minCircularity  =  0. 8f;     //斑点的最小圆度
     //斑点的最大圆度，所能表示的float类型的最大值
    maxCircularity  = std : :numeric_limits < float > : :max();
 
    filterByInertia  =  true;     //斑点惯性率的限制变量
    minInertiaRatio  =  0. 1f;     //斑点的最小惯性率
    maxInertiaRatio  = std : :numeric_limits < float > : :max();     //斑点的最大惯性率
 
    filterByConvexity  =  true;     //斑点凸度的限制变量
    minConvexity  =  0. 95f;     //斑点的最小凸度
    maxConvexity  = std : :numeric_limits < float > : :max();     //斑点的最大凸度
}

• 阈值：通过使用以minThreshold开始的阈值对源图像进行阈值处理，将源图像转换为多个二进制图像。这些阈值以thresholdStep递增，直到maxThreshold。因此，第一个阈值为minThreshold，第二个阈值为minThreshold + thresholdStep，第三个阈值为minThreshold + 2 x thresholdStep，依此类推；
• 分组：在每个二进制图像中，连接的白色像素被分组在一起。我们称这些二进制blob；
• 合并：计算二进制图像中二进制斑点的中心，并合并比minDistBetweenBlob更近的斑点；
• 中心和半径计算：计算并返回新合并的Blob的中心和半径。

并且可以进一步设置SimpleBlobDetector的参数来过滤所需的Blob类型。

• 按颜色：首先需要设置filterByColor =True。设置blobColor = 0可选择较暗的blob，blobColor = 255可以选择较浅的blob。
• 按大小：可以通过设置参数filterByArea = 1以及minArea和maxArea的适当值来基于大小过滤blob。例如。设置minArea = 100将滤除所有少于100个像素的斑点。
• 按圆度：这只是测量斑点距圆的距离。例如。正六边形的圆度比正方形高。要按圆度过滤，请设置filterByCircularity =1。然后为minCircularity和maxCircularity设置适当的值。圆度定义为（）。圆的为圆度为1，正方形的圆度为PI/4，依此类推。
• 
• 按凸性：凸度定义为（斑点的面积/凸包的面积）。现在，形状的“凸包”是最紧密的凸形，它完全包围了该形状，用不严谨的话来讲，给定二维平面上的点集，凸包就是将最外层的点连接起来构成的凸多边形，它能包含点集中所有的点。直观感受上，凸性越高则里面“奇怪的部分”越少。要按凸度过滤，需设置filterByConvexity = true，minConvexity、maxConvexity应该属于[0,1]，而且maxConvexity> minConvexity。
• 按惯性比：这个词汇比较抽象。我们需要知道Ratio可以衡量形状的伸长程度。简单来说。对于圆，此值是1，对于椭圆，它在0到1之间，对于直线，它是0。按惯性比过滤，设置filterByInertia = true，并设置minInertiaRatio、maxInertiaRatio同样属于[0,1]并且maxConvexity> minConvexity。

  按凸性（左低右高）	按惯性比（左低右高）

这里的基础知识可能比较复杂，关键是默认参数下，识别的效果应该说出奇的好。

cv : :Ptr <cv : :SimpleBlobDetector > detector  = cv : :SimpleBlobDetector : :create();

但是存在的主要问题是由于blob分析的配置参数太多，优化起来存在困难。同时对于某些情况，明显识别错误。

2.3 基本轮廓分析

更普通的情况下，我们还是需要从轮廓分析开始，通过上面提出的“圆度”来寻找圆，主要是用来“查漏补缺”，或者是用于特殊情况的查找。

三、算法融合、协作增效

在前面已经详细分析3种主要算法的基础上，本文的重点创造一个“算法链”找到的目标有效地融合起来，并且进一步横向分析研究算法间的关系，希望多少能够给关注这个方向、有类似需求的创作者一些思考。

算法流程

首先对于自然图片，通过blod detection获得准确的半径； 而后基于准确的半径，分别调用HoughCircle以查漏补缺； 最后，以上获得的结果，需要进行融合筛选。这个方法，我在“钢管识别”项目上得到了突出的成功应用，最终能够实现非常高的准确识别。主要是基于以下几点：

1、blobdetector能够找到准确的圆的半径，但是会找错、找漏；

2、HoughCircle在有“准确的圆的半径”的加持下，能够很大程度上提高准确识别效率；

3、目标物体是有“固有特征”的，比如这里需要寻找的钢管，他们的“半径”基本上是一致的。

如果对于这个项目感兴趣，可以在51cto课程上搜索，感谢阅读至此，希望有所帮助。

来自为知笔记(Wiz)