代码
/*===============================图像分割=====================================*//*---------------------------------------------------------------------------*//*手动设置阀值*/ IplImage* binaryImg = cvCreateImage(cvSize(w, h),IPL_DEPTH_8U, 1); cvThreshold(smoothImgGauss,binaryImg,71,255,CV_THRESH_BINARY); cvNamedWindow("cvThreshold", CV_WINDOW_AUTOSIZE ); cvShowImage( "cvThreshold", binaryImg );//cvReleaseImage(&binaryImg); /*---------------------------------------------------------------------------*//*自适应阀值 //计算像域邻域的平均灰度，来决定二值化的值*/ IplImage* adThresImg = cvCreateImage(cvSize(w, h),IPL_DEPTH_8U, 1);double max_value=255;int adpative_method=CV_ADAPTIVE_THRESH_GAUSSIAN_C;//CV_ADAPTIVE_THRESH_MEAN_C int threshold_type=CV_THRESH_BINARY;int block_size=3;//阈值的象素邻域大小 int offset=5;//窗口尺寸 cvAdaptiveThreshold(smoothImgGauss,adThresImg,max_value,adpative_method,threshold_type,block_size,offset); cvNamedWindow("cvAdaptiveThreshold", CV_WINDOW_AUTOSIZE ); cvShowImage( "cvAdaptiveThreshold", adThresImg ); cvReleaseImage(&adThresImg);/*---------------------------------------------------------------------------*//*最大熵阀值分割法*/ IplImage* imgMaxEntropy = cvCreateImage(cvGetSize(imgGrey),IPL_DEPTH_8U,1); MaxEntropy(smoothImgGauss,imgMaxEntropy); cvNamedWindow("MaxEntroyThreshold", CV_WINDOW_AUTOSIZE ); cvShowImage( "MaxEntroyThreshold", imgMaxEntropy );//显示图像 cvReleaseImage(&imgMaxEntropy ); /*---------------------------------------------------------------------------*//*基本全局阀值法*/ IplImage* imgBasicGlobalThreshold = cvCreateImage(cvGetSize(imgGrey),IPL_DEPTH_8U,1); cvCopyImage(srcImgGrey,imgBasicGlobalThreshold);int pg[256],i,thre; for (i=0;i<256;i++) pg[i]=0;for (i=0;i<imgBasicGlobalThreshold->imageSize;i++) // 直方图统计 pg[(BYTE)imgBasicGlobalThreshold->imageData[i]]++; thre = BasicGlobalThreshold(pg,0,256); // 确定阈值 cout<<"The Threshold of this Image in BasicGlobalThreshold is:"<<thre<<endl;//输出显示阀值 cvThreshold(imgBasicGlobalThreshold,imgBasicGlobalThreshold,thre,255,CV_THRESH_BINARY); // 二值化 cvNamedWindow("BasicGlobalThreshold", CV_WINDOW_AUTOSIZE ); cvShowImage( "BasicGlobalThreshold", imgBasicGlobalThreshold);//显示图像 cvReleaseImage(&imgBasicGlobalThreshold);/*---------------------------------------------------------------------------*//*OTSU*/ IplImage* imgOtsu = cvCreateImage(cvGetSize(imgGrey),IPL_DEPTH_8U,1); cvCopyImage(srcImgGrey,imgOtsu);int thre2; thre2 = otsu2(imgOtsu); cout<<"The Threshold of this Image in Otsu is:"<<thre2<<endl;//输出显示阀值 cvThreshold(imgOtsu,imgOtsu,thre2,255,CV_THRESH_BINARY); // 二值化 cvNamedWindow("imgOtsu", CV_WINDOW_AUTOSIZE ); cvShowImage( "imgOtsu", imgOtsu);//显示图像 cvReleaseImage(&imgOtsu);/*---------------------------------------------------------------------------*//*上下阀值法：利用正态分布求可信区间*/ IplImage* imgTopDown = cvCreateImage( cvGetSize(imgGrey), IPL_DEPTH_8U, 1 ); cvCopyImage(srcImgGrey,imgTopDown); CvScalar mean ,std_dev;//平均值、 标准差double u_threshold,d_threshold; cvAvgSdv(imgTopDown,&mean,&std_dev,NULL); u_threshold = mean.val[0] +2.5* std_dev.val[0];//上阀值 d_threshold = mean.val[0] -2.5* std_dev.val[0];//下阀值//u_threshold = mean + 2.5 * std_dev; //错误//d_threshold = mean - 2.5 * std_dev; cout<<"The TopThreshold of this Image in TopDown is:"<<d_threshold<<endl;//输出显示阀值 cout<<"The DownThreshold of this Image in TopDown is:"<<u_threshold<<endl; cvThreshold(imgTopDown,imgTopDown,d_threshold,u_threshold,CV_THRESH_BINARY_INV);//上下阀值 cvNamedWindow("imgTopDown", CV_WINDOW_AUTOSIZE ); cvShowImage( "imgTopDown", imgTopDown);//显示图像 cvReleaseImage(&imgTopDown);/*---------------------------------------------------------------------------*//*迭代法*/ IplImage* imgIteration = cvCreateImage( cvGetSize(imgGrey), IPL_DEPTH_8U, 1 ); cvCopyImage(srcImgGrey,imgIteration);int thre3,nDiffRec; thre3 =DetectThreshold(imgIteration, 100, nDiffRec); cout<<"The Threshold of this Image in imgIteration is:"<<thre3<<endl;//输出显示阀值 cvThreshold(imgIteration,imgIteration,thre3,255,CV_THRESH_BINARY_INV);//上下阀值 cvNamedWindow("imgIteration", CV_WINDOW_AUTOSIZE ); cvShowImage( "imgIteration", imgIteration); cvReleaseImage(&imgIteration);/*===============================图像分割== 
