OpenCV计算连通区域数目与最大连通区域并标示出

申明，本文非笔者原创，原文转载自：http://blog.sina.com.cn/s/blog_5ca0198e0101hrt8.html

[cpp] view plaincopy

1. #include   
2. #include   
3. #include   
4.   
5. #pragma comment(lib, "cv.lib")  
6. #pragma comment(lib, "cxcore.lib")  
7. #pragma comment(lib, "highgui.lib")  
8.   
9. int main( int argc, char** argv )    
10. {  
11.     IplImage* src = cvLoadImage(".\\wind.png", CV_LOAD_IMAGE_GRAYSCALE);  
12.     IplImage* dst = cvCreateImage(cvGetSize(src), 8, 3);  
13.     CvMemStorage* storage = cvCreateMemStorage(0);  
14.     CvSeq* contour = 0;  
15.     cvThreshold(src, src,120, 255, CV_THRESH_BINARY);   // 二值化  
16.     cvNamedWindow("Source", 1);  
17.     cvShowImage("Source", src);  
18.     // 提取轮廓  
19.     int contour_num = cvFindContours(src, storage, &contour, sizeof(CvContour), CV_RETR_CCOMP, CV_CHAIN_APPROX_SIMPLE);  
20.     cvZero(dst);        // 清空数组  
21.     CvSeq *_contour = contour;   
22.     double maxarea = 0;  
23.     double minarea = 100;  
24.     int m = 0;  
25.     for( ; contour != 0; contour = contour->h_next )    
26.     {    
27.   
28.         double tmparea = fabs(cvContourArea(contour));  
29.         if(tmparea < minarea)     
30.         {    
31.             cvSeqRemove(contour, 0); // 删除面积小于设定值的轮廓  
32.             continue;  
33.         }    
34.         CvRect aRect = cvBoundingRect( contour, 0 );   
35.         if ((aRect.width/aRect.height)<1)    
36.         {    
37.             cvSeqRemove(contour, 0); //删除宽高比例小于设定值的轮廓  
38.             continue;  
39.         }    
40.         if(tmparea > maxarea)    
41.         {    
42.             maxarea = tmparea;  
43.         }    
44.         m++;  
45.         // 创建一个色彩值  
46.         CvScalar color = CV_RGB( 0, 255, 255 );  
47.   
48.         //max_level 绘制轮廓的最大等级。如果等级为0，绘制单独的轮廓。如果为1，绘制轮廓及在其后的相同的级别下轮廓  
49.         //如果值为2，所有的轮廓。如果等级为2，绘制所有同级轮廓及所有低一级轮廓，诸此种种  
50.         //如果值为负数，函数不绘制同级轮廓，但会升序绘制直到级别为abs(max_level)-1的子轮廓  
51.         cvDrawContours(dst, contour, color, color, -1, 1, 8);   //绘制外部和内部的轮廓  
52.     }    
53.     contour = _contour;  
54.     int count = 0;  
55.     for(; contour != 0; contour = contour->h_next)  
56.     {    
57.         count++;  
58.         double tmparea = fabs(cvContourArea(contour));  
59.         if (tmparea == maxarea)    
60.         {    
61.             CvScalar color = CV_RGB( 255, 0, 0);  
62.             cvDrawContours(dst, contour, color, color, -1, 1, 8);  
63.         }    
64.     }    
65.     printf("The total number of contours is:%d", count);  
66.     cvNamedWindow("Components", 1);  
67.     cvShowImage("Components", dst);  
68.     cvWaitKey(0);  
69.     cvDestroyWindow("Source");  
70.     cvReleaseImage(&src);  
71.     cvDestroyWindow("Components");  
72.     cvReleaseImage(&dst);  
73.   
74.     return 0;  
75. }