分水岭算法VC实现.docVIP

函数名：?????????????? 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())???????? {?????????