霍夫圆变换

opencv中实现霍夫圆变换是通过霍夫梯度法。该原理是：首先对图像应用边沿检测（cvCanny）。然后，对边缘图像中每一个非0点，考虑其局部梯度（Sobel()函数计算X、Y方向的sobel一阶导数得到梯度）。利用得到的梯度，由斜率指定的直线上的每一个点在累加器呗累计，这里斜率是从一个指定的最小值到指定的最大值的距离。同时表姐边缘图像中的每一个非零像素的位置。然后从二维累加器中这些点中选择候选中心，这些中心大于给定阈值并且大于其所有近邻。候选中心在累加器中式降序排列的，以便于最支持像素中心的首先出现。接下来对每一个中心，考虑所有非零像素与中心的距离，选择最支持的一条半径。

opencv中霍夫圆检测函数：

CvSeq * cvHoughCircles (
CvArr *image,
void * Circle_storage,
int method,
double dp,
double min_dist,
double param1 = 100,
double param2 = 300,
int min_radius = 0,
int max_radius = 0,
);

circle_strorage 可以是CV_32FC3的单列数组，三个通道分别存储圆的位置和半径。或者内存存储器（memory storage），圆将变成一个序列CvSeq,由cvHoughCircle()返回一个指向这个序列的指针，同时将method 参数设为：CV_HOUGH_GRADIENT.

dp 是累加器图像的分辨率，如果为1 ，则分辨率是相同的，如果设置更大的值（例如2），累加器的分辨率受此影响会变小（此情况为一半）。dp值大于等于1.

min_dist可以区分两个不同圆之间的最小距离。

param1,param2是边缘阈值、累加器阈值。

实验代码：

#include <stdio.h>
#include <cv.h>
#include "highgui.h"



int main()
{
	cvNamedWindow("grayimg",1);
	//cvNamedWindow("edges",1);
	cvNamedWindow("hough",1);
	cvMoveWindow("grayimg",0,0);
	//cvMoveWindow("edges",300,60);
	cvMoveWindow("hough",30,0);
	IplImage * img = cvLoadImage("e:/circle.jpg");
	IplImage * gray = NULL,*edges = NULL, *circleimg = NULL;
	CvSize size;
	size.width = img->width;
	size.height = img->height;
	CvMemStorage * storage = cvCreateMemStorage(0);

	CvSeq *result = 0;

	gray = cvCreateImage(size,img->depth,1);
	edges = cvCreateImage(size,img->depth,1);
	circleimg = cvCreateImage(cvGetSize(img),img->depth,img->nChannels);
	cvCvtColor(img,gray,CV_RGB2GRAY);
	//cvCanny(gray,edges,50,100,3);
	
	cvCopy(img,circleimg);
	
	result = cvHoughCircles(gray,storage, CV_HOUGH_GRADIENT, 2, gray->width/10 );
	

	int index; // index 为直线索引
	for (index = 0; index < result->total; index ++)
	{
		float *p = (float *)cvGetSeqElem(result,index);
		CvPoint pt = cvPoint(cvRound( p[0]), cvRound(p[1]));
		cvCircle(circleimg, pt, cvRound(p[2]), CV_RGB(255,0,0));

		
	}


	cvShowImage("grayimg",gray);
	cvShowImage("edges",edges);
	cvShowImage("circleimg",circleimg);
	cvWaitKey();
	cvDestroyAllWindows();
	cvReleaseImage( &gray );
	cvReleaseImage(&img);
	//cvReleaseImage( &edges );
	system("pause");
	return 0;
	
}

result: