OpenCV实现分水岭算法

OpenCV实现分水岭算法:

  1. //  分水岭算法原理   
  2. //     
  3.   
  4. IplImage* marker_mask = 0;  
  5. IplImage* markers = 0;  
  6. //IplImage* img0 = 0, *img = 0, *img_gray = 0, *wshed = 0;   
  7. IplImage  *img_gray = 0, *wshed = 0;  
  8. CvPoint prev_pt = {-1,-1};  
  9.   
  10. void on_mouse( int event, int x, int y, int flags, void* param )  
  11. {  
  12.     if( !img )  
  13.         return;  
  14.   
  15.     if( event == CV_EVENT_LBUTTONUP || !(flags & CV_EVENT_FLAG_LBUTTON) )  
  16.         prev_pt = cvPoint(-1,-1);  
  17.     else if( event == CV_EVENT_LBUTTONDOWN )  
  18.         prev_pt = cvPoint(x,y);  
  19.     else if( event == CV_EVENT_MOUSEMOVE && (flags & CV_EVENT_FLAG_LBUTTON) )  
  20.     {  
  21.         CvPoint pt = cvPoint(x,y);  
  22.         if( prev_pt.x < 0 )  
  23.             prev_pt = pt;  
  24.         cvLine( marker_mask, prev_pt, pt, cvScalarAll(255), 5, 8, 0 );  
  25.         cvLine( img, prev_pt, pt, cvScalarAll(255), 5, 8, 0 );  
  26.         prev_pt = pt;  
  27.         cvShowImage( "image", img );  
  28.     }  
  29. }  
  30.   
  31. void CCVMFCView::OnWatershed()//分水岭   
  32. {  
  33.     int   flag=0;  
  34.     CvRNG rng = cvRNG(-1);  
  35.   
  36.     img0 = cvCloneImage( workImg );         //  建立工作位图   
  37.     cvFlip(img0);  
  38.   
  39.     cvNamedWindow( "image", 1 );  
  40.     //  cvNamedWindow( "watershed transform", 1 );   
  41.   
  42.     img = cvCloneImage( img0 );  
  43.     img_gray = cvCloneImage( img0 );  
  44.     wshed = cvCloneImage( img0 );  
  45.     marker_mask = cvCreateImage( cvGetSize(img), 8, 1 );  
  46.     markers = cvCreateImage( cvGetSize(img), IPL_DEPTH_32S, 1 );  
  47.     cvCvtColor( img, marker_mask, CV_BGR2GRAY );  
  48.     cvCvtColor( marker_mask, img_gray, CV_GRAY2BGR );  
  49.   
  50.     cvZero( marker_mask );  
  51.     cvZero( wshed );  
  52.     cvShowImage( "image", img );  
  53.     //  cvShowImage( "watershed transform", wshed );   
  54.     cvSetMouseCallback( "image", on_mouse, 0 );  
  55.   
  56.     m_ImageType=-3;  
  57.     for(;;)  
  58.     {  
  59.         int c = cvWaitKey(0);  
  60.   
  61.         if( c == 27 ) {  
  62.             if (!flag) {                    //  未加标记   
  63.                 wshed = cvCloneImage( img0 );  
  64.             }  
  65.             break;  
  66.         }  
  67.   
  68.         if( c == 'r' )  
  69.         {  
  70.             cvZero( marker_mask );  
  71.             cvCopy( img0, img );  
  72.             cvShowImage( "image", img );  
  73.         }  
  74.   
  75.         if( c == 'w' || c == '\r' )  
  76.         {  
  77.             CvMemStorage* storage = cvCreateMemStorage(0);  
  78.             CvSeq* contours = 0;  
  79.             CvMat* color_tab;  
  80.             int i, j, comp_count = 0;  
  81.             //cvSaveImage( "wshed_mask.png", marker_mask );   
  82.             //marker_mask = cvLoadImage( "wshed_mask.png", 0 );   
  83.             cvFindContours( marker_mask, storage, &contours, sizeof(CvContour),  
  84.                 CV_RETR_CCOMP, CV_CHAIN_APPROX_SIMPLE );  
  85.   
  86.             CvSeq* contourn = contours;  
  87.             int n;  
  88.             for (n=0; contourn != 0; contourn = contourn->h_next,n++) {  
  89.             }                               //  检查边界数   
  90.             if (n) {                        //  已作标记才进行处理                
  91.                 cvZero( markers );  
  92.                 for( ; contours != 0; contours = contours->h_next, comp_count++ )  
  93.                 {  
  94.                     cvDrawContours( markers, contours, cvScalarAll(comp_count+1),  
  95.                         cvScalarAll(comp_count+1), -1, -1, 8, cvPoint(0,0) );  
  96.                 }  
  97.   
  98.                 color_tab = cvCreateMat( 1, comp_count, CV_8UC3 );  
  99.                 for( i = 0; i < comp_count; i++ )  
  100.                 {  
  101.                     uchar* ptr = color_tab->data.ptr + i*3;  
  102.                     ptr[0] = (uchar)(cvRandInt(&rng)%180 + 50);  
  103.                     ptr[1] = (uchar)(cvRandInt(&rng)%180 + 50);  
  104.                     ptr[2] = (uchar)(cvRandInt(&rng)%180 + 50);  
  105.                 }  
  106.   
  107.                 {  
  108.                     double t = (double)cvGetTickCount();  
  109.                     cvWatershed( img0, markers );  //  分水岭算法处理   
  110.                     t = (double)cvGetTickCount() - t;  
  111.                     //                  printf( "exec time = %gms\n", t/(cvGetTickFrequency()*1000.) );   
  112.                 }  
  113.   
  114.                 // paint the watershed image   
  115.                 for( i = 0; i < markers->height; i++ ) {  
  116.                     for( j = 0; j < markers->width; j++ )  
  117.                     {  
  118.                         int idx = CV_IMAGE_ELEM( markers, int, i, j );  
  119.                         uchar* dst = &CV_IMAGE_ELEM( wshed, uchar, i, j*3 );  
  120.                         if( idx == -1 )  
  121.                             dst[0] = dst[1] = dst[2] = (uchar)255;  
  122.                         else if( idx <= 0 || idx > comp_count )  
  123.                             dst[0] = dst[1] = dst[2] = (uchar)0; // should not get here   
  124.                         else  
  125.                         {  
  126.                             uchar* ptr = color_tab->data.ptr + (idx-1)*3;  
  127.                             dst[0] = ptr[0]; dst[1] = ptr[1]; dst[2] = ptr[2];  
  128.                         }  
  129.                     }  
  130.                 }  
  131.   
  132.                 cvAddWeighted( wshed, 0.5, img_gray, 0.5, 0, wshed );  //  图像合成   
  133.                 //              cvShowImage( "watershed transform", wshed );   
  134.                 cvReleaseMemStorage( &storage );  
  135.                 cvReleaseMat( &color_tab );  
  136.             }  
  137.             else {                          //  未加标记   
  138.                 wshed = cvCloneImage( img0 );  
  139.             }  
  140.             cvCopy(wshed,workImg);  
  141.             cvFlip(workImg);  
  142.   
  143.             CClientDC dc(this);  
  144.             StretchDIBits(dc.m_hDC,         //  刷新主窗口   
  145.                 0,0,workImg->width,workImg->height,  
  146.                 0,0,workImg->width,workImg->height,  
  147.                 workImg->imageData,m_lpBmi,DIB_RGB_COLORS,SRCCOPY);  
  148.             flag=1;  
  149.         }  
  150.     }  
  151.   
  152.     cvDestroyWindow( "image" );  
  153.     cvReleaseImage(&img0);  
  154.     cvReleaseImage(&img);  
  155.     cvReleaseImage(&img_gray);  
  156.     cvReleaseImage(&marker_mask);  
  157.     cvReleaseImage(&markers);  
  158.   
  159.     cvFlip(wshed);  
  160.     m_dibFlag=imageReplace(wshed,&workImg);  
  161.   
  162.     Invalidate();  
  163. }  

相关推荐