1 /* This is a standalone program. Pass an image name as a first parameter of the program.
2 Switch between standard and probabilistic Hough transform by changing "#if 1" to "#if 0" and back */
7 int main(int argc, char** argv)
9 const char* filename = argc >= 2 ? argv[1] : "pic1.png";
10 IplImage* src = cvLoadImage( filename, 0 );
13 CvMemStorage* storage = cvCreateMemStorage(0);
20 dst = cvCreateImage( cvGetSize(src), 8, 1 );
21 color_dst = cvCreateImage( cvGetSize(src), 8, 3 );
23 cvCanny( src, dst, 50, 200, 3 );
24 cvCvtColor( dst, color_dst, CV_GRAY2BGR );
26 lines = cvHoughLines2( dst, storage, CV_HOUGH_STANDARD, 1, CV_PI/180, 100, 0, 0 );
28 for( i = 0; i < MIN(lines->total,100); i++ )
30 float* line = (float*)cvGetSeqElem(lines,i);
32 float theta = line[1];
34 double a = cos(theta), b = sin(theta);
35 double x0 = a*rho, y0 = b*rho;
36 pt1.x = cvRound(x0 + 1000*(-b));
37 pt1.y = cvRound(y0 + 1000*(a));
38 pt2.x = cvRound(x0 - 1000*(-b));
39 pt2.y = cvRound(y0 - 1000*(a));
40 cvLine( color_dst, pt1, pt2, CV_RGB(255,0,0), 3, CV_AA, 0 );
43 lines = cvHoughLines2( dst, storage, CV_HOUGH_PROBABILISTIC, 1, CV_PI/180, 50, 50, 10 );
44 for( i = 0; i < lines->total; i++ )
46 CvPoint* line = (CvPoint*)cvGetSeqElem(lines,i);
47 cvLine( color_dst, line[0], line[1], CV_RGB(255,0,0), 3, CV_AA, 0 );
50 cvNamedWindow( "Source", 1 );
51 cvShowImage( "Source", src );
53 cvNamedWindow( "Hough", 1 );
54 cvShowImage( "Hough", color_dst );