分水岭算法源代码及其分析和注释.doc

/*====================================================================函数名：??????????? Watershed功能：????????????? 用标记-分水岭算法对输入图像进行分割算法实现：????????? 无输入参数说明：????? OriginalImage --输入图像（灰度图，0~255）??????????????????? SeedImage?????? --标记图像（二值图，0-非标记，1-标记）??????????????????? LabelImage????? --输出图像（1-第一个分割区域，2-第二个分割区域，...）??????????????????? row???????????? --图像行数??????????????????? col???????????? --图像列数返回值说明：??????? 无????????====================================================================*/void WINAPI CDib::Watershed(unsigned char **OriginalImage, char** SeedImage, int **LabelImage, int row, int col){// using namespace std;//标记区域标识号，从1开始int Num=0;int i,j;//保存每个队列种子个数的数组vectorint* SeedCounts;//临时种子队列queuePOINT quetem;//保存所有标记区域种子队列的数组，里面放的是种子队列的指针vectorqueuePOINT* vque;int* array;//指向种子队列的指针queuePOINT *pque;POINT temp;for(i=0;irow;i++){??? for(j=0;jcol;j++)???? LabelImage[i][j]=0;}int m,n,k=0;BOOL up,down,right,left,upleft,upright,downleft,downright;//8 directions...//预处理,提取区分每个标记区域，并初始化每个标记的种子队列//种子是指标记区域边缘的点，他们可以在水位上升时向外淹没（或者说生长）//pans words：我的理解是梯度值较小的象素点，或者是极小灰度值的点。for(i=0;irow;i++){??? for(j=0;jcol;j++)??? {???? //如果找到一个标记区域???? if(SeedImage[i][j]==1)???? {????? //区域的标识号加一????? Num++;????? //分配数组并初始化为零，表示可有256个灰阶????? array=new int[256];????? ZeroMemory(array,256*sizeof(int));????? //种子个数数组进vector，每次扫描则生成一个数组，并用区域标识号来做第一维。灰度级做第二维。????? //表示某个盆地区域中某灰阶所对应的点的数目。????? SeedCounts.push_back(array);????? //分配本标记号的优先队列，256个种子队列，????? //表示对应一个灰阶有一个队列，并且每个队列可以存储一个集合的点信息????? pque=new queuePOINT[256];????? //加入到队列数组中，对应的是本标记号Num的????? vque.push_back(pque);????? //当前点放入本标记区域的临时种子队列中????? temp.x=i;????? temp.y=j;????? quetem.push(temp);????? //当前点标记为已处理????? LabelImage[i][j]=Num;????? SeedImage[i][j]=127;//表示已经处理过????????? //让临时种子队列中的种子进行生长直到所有的种子都生长完毕????? //生长完毕后的队列信息保存在vque中，包括区域号和灰阶，对应点数存储在seedcounts中????? while(!quetem.empty())????? {?????? up=down=right=left=FALSE;?????? upleft=upright=downleft=d