图像分割(阈值分割和区域生长)vc++代码.doc

实验6 图像分割 阈值分割法 实验代码及结果： #include cv.h #include highgui.h #include stdio.h #include stdlib.h #include math.h int main( int argc, char** argv ) { IplImage* pImg=0; //声明IplImage指针 int height,width,step,channels; uchar *data; int i,j; //载入图像 pImg = cvLoadImage(2.jpg,0); if( !pImg ) { printf(Image was not loaded.\n); return -1; } //获取图像信息 height=pImg-height; width=pImg-width; step=pImg-widthStep; channels=pImg-nChannels; data=(uchar *)pImg-imageData; printf(Processing a %d*%d with %d channels\n,height,width,channels); int MaxP=0; float t,T,m=0,n=0; float w0,w1,u0,u1,u; double q2=0,temp=0; cvNamedWindow( win1, 1);//创建窗口 cvShowImage( win1, pImg );//显示图像 //计算最大灰度值MaxP for(i=0;iheight;i++) for(j=0;jwidth;j++) { if(data[i*width+j]MaxP) {MaxP=data[i*step+j];} } printf(最大灰度值：%d\n,MaxP); //计算最大方差时的阈值t for(T=0;TMaxP;T++) { w0=0,w1=0,u0=0,u1=0,u=0,q2=0,m=0,n=0; for(i=0;iheight;i++) { for(j=0;jwidth;j++) { if(data[i*width+j]T) {u0+=data[i*width+j];m++;} else {u1+=data[i*width+j];n++;} } } w0=m/(width*height); w1=1-w0; u0/=m; u1/=n; u=w0*u0+w1*u1; q2=w0*w1*pow(u0-u1,2); if(q2temp) { t=T;temp=q2;} } printf(最佳阈值：%lf\n,t); for(i=0;iheight;i++) for(j=0;jwidth;j++) { if(data[i*width+j]t) data[i*width+j]=0; else data[i*width+j]=255; } cvNamedWindow(win2, 1); cvShowImage( win2, pImg );//显示图像 cvWaitKey(0); cvDestroyWindow( win1 );//销毁窗口 cvDestroyWindow( win2 ); cvReleaseImage( pImg ); //释放图像 return 0; } 区域生长法 实验代码及结果： #include cv.h #include highgui.h #include stdio.h #include stdlib.h #include math.h #includemalloc.h #define MaxSize 1000000 struct queue{ int qu[MaxSize]; int front; int rear; int tag; //front==rear时表示队列满或空的标志：tag==1表示满，tag==0表示空 }; struct queue *InitQu() {struct queue *q; q=(struct queue *)malloc(sizeof(struct queue));//分配空间 q-front=0; q-rea