cvWatershed例子:
#include<cv.h>
#include<highgui.h>
#include<iostream>
using namespace std;
IplImage* marker_mask = 0;
IplImage* markers = 0;
IplImage* img0 = 0, *img = 0, *img_gray = 0, *wshed = 0;
CvPoint prev_pt = {-1,-1};
void on_mouse( int event, int x, int y, int flags, void* param )//opencv 會自動給函數傳入合適的值
{
if( !img )
return;
if( event == CV_EVENT_LBUTTONUP || !(flags & CV_EVENT_FLAG_LBUTTON) )
prev_pt = cvPoint(-1,-1);
else if( event == CV_EVENT_LBUTTONDOWN )
prev_pt = cvPoint(x,y);
else if( event == CV_EVENT_MOUSEMOVE && (flags & CV_EVENT_FLAG_LBUTTON) )
{
CvPoint pt = cvPoint(x,y);
if( prev_pt.x < 0 )
prev_pt = pt;
cvLine( marker_mask, prev_pt, pt, cvScalarAll(255), 5, 8, 0 );//CvScalar 成員:double val[4] RGBA值A=alpha
cvLine( img, prev_pt, pt, cvScalarAll(255), 5, 8, 0 );
prev_pt = pt;
cvShowImage( "image", img);
}
}
int main( int argc, char** argv )
{
char* filename = argc >= 2 ? argv[1] : (char*)"fruits.jpg";
CvMemStorage* storage = cvCreateMemStorage(0);
CvRNG rng = cvRNG(-1);
if( (img0 = cvLoadImage(filename,1)) == 0 )
return 0;
printf( "Hot keys: \n"
"\tESC - quit the program\n"
"\tr - restore the original image\n"
"\tw or SPACE - run watershed algorithm\n"
"\t\t(before running it, roughly mark the areas on the image)\n"
"\t (before that, roughly outline several markers on the image)\n" );
cvNamedWindow( "image", 1 );
cvNamedWindow( "watershed transform", 1 );
img = cvCloneImage( img0 );
img_gray = cvCloneImage( img0 );
wshed = cvCloneImage( img0 );
marker_mask = cvCreateImage( cvGetSize(img), 8, 1 );
markers = cvCreateImage( cvGetSize(img), IPL_DEPTH_32S, 1 );
cvCvtColor( img, marker_mask, CV_BGR2GRAY );
cvCvtColor( marker_mask, img_gray, CV_GRAY2BGR );//這兩句只用將RGB轉成3通道的灰度圖即R=G=B,用來顯示用
cvZero( marker_mask );
cvZero( wshed );
cvShowImage( "image", img );
cvShowImage( "watershed transform", wshed );
cvSetMouseCallback( "image", on_mouse, 0 );
for(;;)
{
int c = cvWaitKey(0);
if( (char)c == 27 )
break;
if( (char)c == 'r' )
{
cvZero( marker_mask );
cvCopy( img0, img );//cvCopy()也可以這樣用,不影響原img0圖像,也隨時更新
cvShowImage( "image", img );
}
if( (char)c == 'w' || (char)c == ' ' )
{
CvSeq* contours = 0;
CvMat* color_tab = 0;
int i, j, comp_count = 0;
//下面選將標記的圖像取得其輪廓, 將每種輪廓用不同的整數表示
//不同的整數使用分水嶺算法時,就成爲不同的種子點
//算法本來就是以各個不同的種子點爲中心擴張
cvClearMemStorage(storage);
cvFindContours( marker_mask, storage, &contours, sizeof(CvContour),
CV_RETR_CCOMP, CV_CHAIN_APPROX_SIMPLE );
cvZero( markers );
for( ; contours != 0; contours = contours->h_next, comp_count++ )
{
cvDrawContours(markers, contours, cvScalarAll(comp_count+1),
cvScalarAll(comp_count+1), -1, -1, 8, cvPoint(0,0) );
}
//cvShowImage("image",markers);
if( comp_count == 0 )
continue;
color_tab = cvCreateMat( 1, comp_count, CV_8UC3 );//創建隨機顏色列表
for( i = 0; i < comp_count; i++ ) //不同的整數標記
{
uchar* ptr = color_tab->data.ptr + i*3;
ptr[0] = (uchar)(cvRandInt(&rng)%180 + 50);
ptr[1] = (uchar)(cvRandInt(&rng)%180 + 50);
ptr[2] = (uchar)(cvRandInt(&rng)%180 + 50);
}
{
double t = (double)cvGetTickCount();
cvWatershed( img0, markers );
cvSave("img0.xml",markers);
t = (double)cvGetTickCount() - t;
printf( "exec time = %gms\n", t/(cvGetTickFrequency()*1000.) );
}
// paint the watershed image
for( i = 0; i < markers->height; i++ )
for( j = 0; j < markers->width; j++ )
{
int idx = CV_IMAGE_ELEM( markers, int, i, j );//markers的數據類型爲IPL_DEPTH_32S
uchar* dst = &CV_IMAGE_ELEM( wshed, uchar, i, j*3 );//BGR三個通道的數是一起的,故要j*3
if( idx == -1 ) //輸出時若爲-1,表示各個部分的邊界
dst[0] = dst[1] = dst[2] = (uchar)255;
else if( idx <= 0 || idx > comp_count ) //異常情況
dst[0] = dst[1] = dst[2] = (uchar)0; // should not get here
else //正常情況
{
uchar* ptr = color_tab->data.ptr + (idx-1)*3;
dst[0] = ptr[0]; dst[1] = ptr[1]; dst[2] = ptr[2];
}
}
cvAddWeighted( wshed, 0.5, img_gray, 0.5, 0, wshed );//wshed.x.y=0.5*wshed.x.y+0.5*img_gray+0加權融合圖像
cvShowImage( "watershed transform", wshed );
cvReleaseMat( &color_tab );
}
}
return 1;
}
運行情況:同上面分析可看出,因爲不相連的標記後,程序在masker中的數值不同
經過分水嶺算法後,
不同的標記肯定會在不同的區域中,
例如頭髮部分,我畫了一條線標記 ,, 處理後就把頭髮部分分割了出來
還比如胳膊那一塊,正好也分割出來了
opencv中的算法是先把輸入圖像轉化成梯度圖(標量)
如果把梯度圖看成是一個地形的話,就會發現,梯度高的地方就成了山脈,梯度低的地方就是山谷
我們經過標記爲不同的區域後,
就從各個標記的地方注水進去,注入的水越來越多的時候,就會出現把流過低些的山脈,從而流到別的山谷中,那麼他們就連一了一片區域。
區域分割的要求是把不同的標記分割成不同的地方。所以如果一直注水,可能就會覆蓋別的區域了。這時算法就採取某種方法,修大壩使標記的不同區域不會因爲注水而相連
他們會互不相干的擴張領地,直到把整個領地都擴張完爲止。
函數輸出時,自己的標記擴張出來的區域都用之前標記的值表示,代表一個區域。 不同的值代表不同的區域
區域與區域之間的邊界由由值-1表示.