图像算法：图像阈值分割

代码

/*===============================图像分割=====================================*/
/*---------------------------------------------------------------------------*/
/*手动设置阀值*/
    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);

转自：http://www.cnblogs.com/skyseraph/archive/2010/12/21/1913058.html