用OpenCv找出最大连通域，并对其进行连通域进行ROI

用OpenCv找出最大连通域，并对其进行连通域进行ROI

(2011-08-21 21:53:33)

转载▼

标签： 连通域 标定 opencv it

分类：OpenCV

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

#include "cv.h"
#include "cxcore.h"
#include "highgui.h"
 
int main( int argc, char** argv )
{
 //声明IplImage指针
 IplImage* pImg = cvLoadImage("e:/black.jpg",0);
 IplImage* pContourImg = NULL;
 CvMemStorage * storage = cvCreateMemStorage(0);
 CvSeq * contour = 0;
 CvSeq *contmax = 0;
 int mode = CV_RETR_EXTERNAL;
 cvShowImage( "src", pImg );
 //为轮廓显示图像申请空间
 //3通道图像，以便用彩色显示
 pContourImg = cvCreateImage(cvGetSize(pImg),
  IPL_DEPTH_8U,
  3);
 //copy source image and convert it to BGR image
 cvCvtColor(pImg, pContourImg, CV_GRAY2BGR);
 //查找contour
 cvFindContours( pImg, storage, &contour, sizeof(CvContour),
  mode, CV_CHAIN_APPROX_SIMPLE, cvPoint(0,0));
 //将轮廓画出   
 cvDrawContours(pContourImg, contour,
  CV_RGB(255,0,0), CV_RGB(255, 0, 0),
  2, 2, 8, cvPoint(0,0));
 int area,maxArea = 10;//设面积最大值大于10Pixel
 for(;contour;contour = contour->h_next)
 {
  area = fabs(cvContourArea( contour, CV_WHOLE_SEQ )); //获取当前轮廓面积
  printf("area == %lf\n", area);
  if(area > maxArea)
  {
   contmax = contour;
   maxArea = area;
  }
 }
 CvRect aRect = cvBoundingRect( contmax, 0 );
 cvSetImageROI( pContourImg,aRect);
 //显示图像
 cvShowImage( "contour", pContourImg );
 cvSaveImage("e:/contour.jpg",pContourImg);
 cvWaitKey(0);

 //销毁窗口
 cvDestroyWindow( "src" );
 cvDestroyWindow( "contour" );
 //释放图像
 cvReleaseImage( &pImg );
 cvReleaseImage( &pContourImg );
 cvReleaseMemStorage(&storage);

 return 0;
}

处理前的连通域

处理后的连通域

 

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

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

[cpp] view plain copy print ?

1.    

 

[cpp] view plain copy print ?

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