Hough变换原理与Opencv实现

转自：https://blog.csdn.net/ycj9090900/article/details/52944708

霍夫变换(Hough Transform)是图像处理中的一种特征提取技术，它通过一种投票算法检测具有特定形状的物体。该过程在一个参数空间中通过计算累计结果的局部最大值得到一个符合该特定形状的集合作为霍夫变换结果。霍夫变换于1962年由Paul Hough 首次提出[53]，后于1972年由Richard Duda和Peter Hart推广使用[54]，经典霍夫变换用来检测图像中的直线，后来霍夫变换扩展到任意形状物体的识别，多为圆和椭圆.

经过几天的学习，发现各位大牛的理解方式之前都是有一些区别的，但是核心的思想没有变化，因此记录一下自己对霍夫变换直线检测的认识。

一、原理介绍：

1、对于直角坐标系中的任意一点A(x0,y0)，经过点A的直线满足Y0=k*X0+b.(k是斜率，b是截距)

2、那么在X-Y平面过点A(x0,y0)的直线簇可以用Y0=k*X0+b表示，但对于垂直于X轴的直线斜率是无穷大的则无法表示。因此将直角坐标系转换到极坐标系就能解决该特殊情况。

3、在极坐标系中表示直线的方程为ρ=xCosθ+ySinθ（ρ为原点到直线的距离）,如图所示：

4、如上图，假定在一个8*8的平面像素中有一条直线，并且从左上角（1,8）像素点开始分别计算θ为0°、45°、90°、135°、180°时的ρ，图中可以看出ρ分别为1、(9√2)/2、8、(7√2)/2、-1，并给这5个值分别记一票，同理计算像素点(3,6)点θ为0°、45°、90°、135°、180°时的ρ，再给计算出来的5个ρ值分别记一票，此时就会发现ρ = (9√2)/2的这个值已经记了两票了，以此类推，遍历完整个8*8的像素空间的时候ρ = (9√2)/2就记了5票， 别的ρ值的票数均小于5票，所以得到该直线在这个8*8的像素坐标中的极坐标方程为 (9√2)/2=x*Cos45°+y*Sin45°，到此该直线方程就求出来了。（PS：但实际中θ的取值不会跨度这么大，一般是PI/180）。

二、Opencv实现直线检测：

1、Opencv1.0版本：

#include
#include

int main()
{
    IplImage* pImgSrc = NULL;    //源图像
    IplImage* pImg8u = NULL;     //灰度图
    IplImage* pImgCanny = NULL;  //边缘检测后的图
    IplImage* pImgDst = NULL;    //在图像上画上检测到的直线后的图像
    CvSeq* lines = NULL;
    CvMemStorage* storage = NULL;

    /*边缘检测*/
    pImgSrc = cvLoadImage(".\\res\\street.jpg", 1);
    pImg8u = cvCreateImage(cvGetSize(pImgSrc), IPL_DEPTH_8U, 1);
    pImgCanny = cvCreateImage(cvGetSize(pImgSrc), IPL_DEPTH_8U, 1);
    pImgDst = cvCreateImage(cvGetSize(pImgSrc), IPL_DEPTH_8U, 1);
    cvCvtColor(pImgSrc, pImg8u, CV_BGR2GRAY);
    cvCanny(pImg8u, pImgCanny, 20, 200, 3);

    /*检测直线*/
    storage = cvCreateMemStorage(0);
    lines = cvHoughLines2(pImgCanny, storage, CV_HOUGH_PROBABILISTIC, 1, CV_PI / 180, 80, 200, 10);
    pImgDst = cvCreateImage(cvGetSize(pImgSrc), IPL_DEPTH_8U, 3);
    cvCvtColor(pImg8u, pImgDst, CV_GRAY2BGR);

    /*在pImgDst上画出检测到的直线*/
    for (int i = 0; i < lines->total; i++)
    {
        CvPoint* line = (CvPoint*)cvGetSeqElem(lines, i);
        cvLine(pImgDst, line[0], line[1], CV_RGB(255, 0, 0), 3, 8);
    }

    cvNamedWindow("src", 1);
    cvNamedWindow("canny", 1);
    cvNamedWindow("hough", 1);
    cvShowImage("src", pImgSrc);
    cvShowImage("canny", pImgCanny);
    cvShowImage("hough", pImgDst);

    cvWaitKey(0);

    cvReleaseImage(&pImgSrc);
    cvReleaseImage(&pImg8u);
    cvReleaseImage(&pImgCanny);
    cvReleaseImage(&pImgDst);
    cvReleaseMemStorage(&storage);

    return 0;
}