OpenCV Hough Circle 變換

目標(biāo)

在本教程中，您將學(xué)習(xí)如何：

• 使用OpenCV函數(shù)cv :: HoughCircles來(lái)檢測(cè)圖像中的圓。

理論

Hough Circle變換

• Hough Circle變換的工作原理與上一個(gè)教程中解釋的Hough Line變換大致相似。
  
• 在行檢測(cè)案例中，一行由兩個(gè)參數(shù)(r,θ)。在圓圈中，我們需要三個(gè)參數(shù)來(lái)定義一個(gè)圓：
  

定義中心位置（綠點(diǎn)），r是半徑，這允許我們完全定義一個(gè)圓，如下所示：

• 為了實(shí)現(xiàn)效率，OpenCV實(shí)現(xiàn)了一個(gè)比標(biāo)準(zhǔn)霍夫變換（Hough Transform）稍微棘手的檢測(cè)方法：霍夫梯度法，由兩個(gè)主要階段組成。第一階段涉及邊緣檢測(cè)并找到可能的圓心，第二階段為每個(gè)候選中心找到最佳的半徑。有關(guān)更多詳細(xì)信息，請(qǐng)查看學(xué)習(xí)OpenCV或您最喜歡的計(jì)算機(jī)視覺(jué)參考書(shū)

Code

• 這個(gè)程序是做什么的？
  

1. 加載圖像并使其模糊，
   
2. 以減少噪音將Hough Circle變換應(yīng)用于模糊圖像。
   
3. 在窗口中顯示檢測(cè)到的圓。
   

• 我們將從這里下載我們將要解釋的示例代碼。這里可以找到一個(gè)稍微推薦的版本（顯示用于更改閾值的軌跡條）。
  

#include "opencv2/imgcodecs.hpp"
#include "opencv2/highgui.hpp"
#include "opencv2/imgproc.hpp"
#include <iostream>
using namespace cv;
using namespace std;
static void help()
{
    cout << "\nThis program demonstrates circle finding with the Hough transform.\n"
            "Usage:\n"
            "./houghcircles <image_name>, Default is ../data/board.jpg\n" << endl;
}
int main(int argc, char** argv)
{
    cv::CommandLineParser parser(argc, argv,
        "{help h ||}{@image|../data/board.jpg|}"
    );
    if (parser.has("help"))
    {
        help();
        return 0;
    }
    string filename = parser.get<string>("@image");
    Mat img = imread(filename, IMREAD_COLOR);
    if(img.empty())
    {
        help();
        cout << "can not open " << filename << endl;
        return -1;
    }
    Mat gray;
    cvtColor(img, gray, COLOR_BGR2GRAY);
    medianBlur(gray, gray, 5);
    vector<Vec3f> circles;
    HoughCircles(gray, circles, HOUGH_GRADIENT, 1,
                 gray.rows/16, // change this value to detect circles with different distances to each other
                 100, 30, 1, 30 // change the last two parameters
                                // (min_radius & max_radius) to detect larger circles
                 );
    for( size_t i = 0; i < circles.size(); i++ )
    {
        Vec3i c = circles[i];
        circle( img, Point(c[0], c[1]), c[2], Scalar(0,0,255), 3, LINE_AA);
        circle( img, Point(c[0], c[1]), 2, Scalar(0,255,0), 3, LINE_AA);
    }
    imshow("detected circles", img);
    waitKey();
    return 0;
}

說(shuō)明

• 加載圖像
  

    string filename = parser.get<string>("@image");
    Mat img = imread(filename, IMREAD_COLOR);
    if(img.empty())
    {
        help();
        cout << "can not open " << filename << endl;
        return -1;
    }

• 將其轉(zhuǎn)換為灰度：
  

    Mat gray;
    cvtColor(img, gray, COLOR_BGR2GRAY);

• 應(yīng)用中值模糊以減少噪音，避免假圈檢測(cè)：
  

    medianBlur(gray, gray, 5);

繼續(xù)應(yīng)用Hough Circle變換：

    vector<Vec3f> circles;
    HoughCircles(gray, circles, HOUGH_GRADIENT, 1,
                 gray.rows/16, // change this value to detect circles with different distances to each other
                 100, 30, 1, 30 // change the last two parameters
                                // (min_radius & max_radius) to detect larger circles
                 );

有論據(jù)：

1. gray：輸入圖像（灰度）。
2. circles：存儲(chǔ)套3個(gè)的值的矢量：為每個(gè)檢測(cè)到的圓。xc,yc,r
3. HOUGH_GRADIENT：定義檢測(cè)方法。目前這是OpenCV中唯一可用的。
4. dp = 1：分辨率的反比。
5. min_dist = gray.rows / 16：檢測(cè)到的中心之間的最小距離。
6. param_1 = 200：內(nèi)部Canny邊緣檢測(cè)器的上限閾值。
7. param_2 = 100 *：中心檢測(cè)閾值。
8. min_radius = 0：要檢測(cè)的最小半徑。如果未知，則將零置為默認(rèn)值。
9. max_radius = 0：要檢測(cè)的最大半徑。如果未知，則將零置為默認(rèn)值。

• 繪制檢測(cè)到的圓：
  

    for( size_t i = 0; i < circles.size(); i++ )
    {
        Vec3i c = circles[i];
        circle( img, Point(c[0], c[1]), c[2], Scalar(0,0,255), 3, LINE_AA);
        circle( img, Point(c[0], c[1]), 2, Scalar(0,255,0), 3, LINE_AA);
    }

你可以看到，我們將用紅色圓圈和中心繪制一個(gè)小綠點(diǎn)

• 顯示檢測(cè)到的圓，并等待用戶(hù)退出程序：
  

    imshow("detected circles", img);
    waitKey();

結(jié)果

使用測(cè)試圖像運(yùn)行代碼的結(jié)果如下所示：