hough变换 正在的c代码风格实现

图像中直线的检测 -- Hough变换原理  

一、简单介绍

Hough变换是图像处理中从图像中识别几何形状的基本方法之一。Hough变换的基本原理在于利用点与线的对偶性，将原始图像空间的给定的曲线通过曲线表达形式变为参数空间的一个点。这样就把原始图像中给定曲线的检测问题转化为寻找参数空间中的峰值问题。也即把检测整体特性转化为检测局部特性。比如直线、椭圆、圆、弧线等。

二、Hough变换的基本思想

设已知一黑白图像上画了一条直线，要求出这条直线所在的位置。我们知道，直线的方程可以用y=k*x+b 来表示，其中k和b是参数，分别是斜率和截距。过某一点(x0,y0)的所有直线的参数都会满足方程y0=kx0+b。即点(x0,y0)确定了一族直线。方程y0=kx0+b在参数k--b平面上是一条直线，(你也可以是方程b=-x0*k+y0对应的直线)。这样，图像x--y平面上的一个前景像素点就对应到参数平面上的一条直线。我们举个例子说明解决前面那个问题的原理。设图像上的直线是y=x, 我们先取上面的三个点：A(0,0), B(1,1), C(22)。可以求出，过A点的直线的参数要满足方程b=0, 过B点的直线的参数要满足方程1=k+b, 过C点的直线的参数要满足方程2=2k+b, 这三个方程就对应着参数平面上的三条直线，而这三条直线会相交于一点(k=1,b=0)。　同理，原图像上直线y=x上的其它点(如(3,3),(4,4)等)　对应参数平面上的直线也会通过点(k=1,b=0)。这个性质就为我们解决问题提供了方法，就是把图像平面上的点对应到参数平面上的线，最后通过统计特性来解决问题。假如图像平面上有两条直线，那么最终在参数平面上就会看到两个峰值点，依此类推。

简而言之，Hough变换思想为：在原始图像坐标系下的一个点对应了参数坐标系中的一条直线，同样参数坐标系的一条直线对应了原始坐标系下的一个点，然后，原始坐标系下呈现直线的所有点，它们的斜率和截距是相同的，所以它们在参数坐标系下对应于同一个点。这样在将原始坐标系下的各个点投影到参数坐标系下之后，看参数坐标系下有没有聚集点，这样的聚集点就对应了原始坐标系下的直线。

在实际应用中， y=k*x+b 形式的直线方程没有办法表示 x=c 形式的直线 ( 这时候，直线的斜率为无穷大 ) 。所以实际应用中，是采用参数方程 p=x*cos(theta)+y*sin(theta) 。这样，图像平面上的一个点就对应到参数 p---theta 平面上的一条曲线上，其它的还是一样。

      以下代码实现了最简单的Hough变换直线检测，输入为width*height的二值图（背景为0，前景为255），存放在矩阵src中，iThreshold为判断为直线的域值。输出pR为原点到直线的距离，pTh为直线的角度
void Hough(BYTE *src,int width,int height, int *pR, int *pTh, int iThreshold)
{
        int *pArray;
        int iRMax = (int)sqrt(width * width + height * height) + 1;
        int iThMax = 361;
        int iTh = 0;
        int iR;
        int iMax = -1;
        int iThMaxIndex = -1;
        int iRMaxIndex = -1;

        pArray = new int[iRMax * iThMax];
        memset(pArray, 0, sizeof(int) * iRMax * iThMax);

        float fRate = (float)(PI/180);

        for (int y = 0; y < height; y++)
        {
                for (int x = 0; x < width; x++)
                {
                        if(*src == 255)
                        {
                                for(iTh = 0; iTh < iThMax; iTh += 1)
                                {
                                        iR = (int)(x * cos(iTh * fRate) + y * sin(iTh * fRate));

                                        if(iR > 0)
                                        {
                                                pArray[iR/1 * iThMax + iTh]++;
                                        }
                                }
                        }

                        src++;
                } // x
        } // y

        for(iR = 0; iR < iRMax; iR++)
        {
                for(iTh = 0; iTh < iThMax; iTh++)
                {
                        int iCount = pArray[iR * iThMax + iTh];
                        if(iCount > iMax)
                        {
                                iMax = iCount;
                                iRMaxIndex = iR;
                                iThMaxIndex = iTh;
                        }
                }
        }

        if(iMax >= iThreshold)
        {
                *pR = iRMaxIndex;
                *pTh = iThMaxIndex;
        }

        delete []pArray;

        return;
} // end of Hough

摘自 <<利用Hough变换实现直线检测的代码>> 和 Hough变换原理
参考 Hough transform