OpenCv的连通域操作

由于项目需要，要对图像中的最大连通域进行标定，并且存储。首先需要使用cvFindCountour对边缘进行标定，其实它的原理就是连通域的边缘提取；其次就是对连通域进行大小判断找出最大的连通域；最后当然就是进行Rect并且ROI了。如果有需要可以进行存储。直接上源码吧。

[cpp] view plain copy

1. #include "cv.h"  
2. #include "cxcore.h"  
3. #include "highgui.h"  
4.    
5. int main( int argc, char** argv )  
6. {  
7.  //声明IplImage指针  
8.  IplImage* pImg = cvLoadImage("black.bmp",0);  
9.  IplImage* pContourImg = NULL;  
10.  CvMemStorage * storage = cvCreateMemStorage(0);  
11.  CvSeq * contour = 0;  
12.  CvSeq *contmax = 0;  
13.  int mode = CV_RETR_EXTERNAL;  
14.  cvShowImage( "src", pImg );  
15.  //为轮廓显示图像申请空间  
16.  //3通道图像，以便用彩色显示  
17.  pContourImg = cvCreateImage(cvGetSize(pImg),  
18.   IPL_DEPTH_8U,  
19.   3);  
20.  //copy source image and convert it to BGR image  
21.  cvCvtColor(pImg, pContourImg, CV_GRAY2BGR);  
22.  //查找contour  
23.  cvFindContours( pImg, storage, &contour, sizeof(CvContour),  
24.   mode, CV_CHAIN_APPROX_SIMPLE, cvPoint(0,0));  
25.  //将轮廓画出     
26.  cvDrawContours(pContourImg, contour,  
27.   CV_RGB(255,0,0), CV_RGB(255, 0, 0),  
28.   2, 2, 8, cvPoint(0,0));  
29.  int area,maxArea = 10;//设面积最大值大于10Pixel  
30.  for(;contour;contour = contour->h_next)  
31.  {  
32.   area = fabs(cvContourArea( contour, CV_WHOLE_SEQ )); //获取当前轮廓面积  
33.   printf("area == %lf\n", area);  
34.   if(area > maxArea)  
35.   {  
36.    contmax = contour;  
37.    maxArea = area;  
38.   }  
39.  }  
40.  CvRect aRect = cvBoundingRect( contmax, 0 );  
41.  cvSetImageROI( pContourImg,aRect);  
42.  //显示图像  
43.  cvShowImage( "contour", pContourImg );  
44.  cvSaveImage("contour.bmp",pContourImg);  
45.  cvWaitKey(0);  
46.   
47.  //销毁窗口  
48.  cvDestroyWindow( "src" );  
49.  cvDestroyWindow( "contour" );  
50.  //释放图像  
51.  cvReleaseImage( &pImg );  
52.  cvReleaseImage( &pContourImg );  
53.  cvReleaseMemStorage(&storage);  
54.   
55.  return 0;  
56. }  

处理前的连通域

处理后的连通域

1. a. 二值化  
2. b. 得到轮廓的个数  
3. c. 将面积小于100的轮廓删除  
4. d. 将宽、高 比例小于1的轮廓删除  
5. e. 把面积最大的米粒用红色框框画出来  

[cpp] view plain copy

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

以下是结果：

//-----------------------------------------------

转自：http://blog.csdn.net/augusdi/article/details/8864748

//------------------------------------------------