wbc特征提取

　　/* CvBox2D skin_rect;
　　CvMemStorage* storage = cvCreateMemStorage(0);
　　CvSeq* contour = 0;
　　CvMemStorage* storage2 = cvCreateMemStorage(0);
　　CvSeq* contour2 = 0;
　　cvFindContours( fore, storage, &contour, sizeof(CvContour),CV_RETR_LIST, CV_CHAIN_APPROX_NONE, cvPoint(0,0)); 
　　CvSeq* contour_max=0;
　　for( ; contour != 0; contour = contour->h_next )
　　{
　　int tmparea=abs(cvContourArea(contour,CV_WHOLE_SEQ));/ /轮廓区域的面积
　　if (tmparea maxarea)
　　{
　　maxarea = tmparea;//最大面积连通区域
　　contour_max=cvCloneSeq(contour);
　　}
　　if(tmpcount > count)
　　{
　　count = tmpcount;//最大轮廓的像素数
　　}
　　if(tmpbox > axisratio)
　　{
　　axisratio = tmpbox;//最大长短轴的比值
　　}
　　}
　　if (contour_max!=0)
　　{
　　CvScalar color = CV_RGB( rand()&255, 0, 0 );
　　cvDrawContours(dst, contour_max, color, color, -1, -1,;
　　ConcavityRepair(dst,3);//对fore图像进行修补
　　cvFindContours( dst, storage2, &contour2, sizeof(CvContour),CV_RETR_LIST, CV_CHAIN_APPROX_NONE, cvPoint(0,0)); 
　　cvDrawContours(src, contour2, color, color, -1, 1,;/*
　　/************************************************* ***********************/
　　/*                            图像轮廓的HU矩                            */
　　/************************************************* ***********************/
　　/*
　　CvMoments m;
　　CvHuMoments hu;
　　//cvMoments(cvGetSubRect(fore,&mat,rct),&m,0);
　　cvMoments(contour2,&m,0);
　　cvGetHuMoments(&m,&hu);
　　huju[0]=(float)fabs(log10(fabs(hu.hu1)));
　　huju[1]=(float)fabs(log10(fabs(hu.hu2)));
　　huju[2]=(float)fabs(log10(fabs(hu.hu3)));
　　huju[3]=(float)fabs(log10(fabs(hu.hu4)));
　　huju[4]=(float)fabs(log10(fabs(hu.hu5)));
　　huju[5]=(float)fabs(log10(fabs(hu.hu6)));
　　huju[6]=(float)fabs(log10(fabs(hu.hu7)));*/
　　/************************************************* ***********************/
　　/*                            图像重心及半径分布                        */
　　/************************************************* ***********************/
　　/*  double m00,x,y;
　　double r_long=0;
　　double r_short=100;
　　int num_long,num_short;//长、短轴对应的弧度
　　float dist_sum=0;//距离和
　　m00=cvGetSpatialMoment(&m,0,0);
　　x=cvGetSpatialMoment(&m,1,0)/m00;  //重心坐标
　　y=cvGetSpatialMoment(&m,0,1)/m00;
　　int contour_max_count=contour2->total;   //this is number point in contour
　　CvPoint* pointarray;
　　float* dist_center=new float[contour_max_count];
　　pointarray=(CvPoint*)malloc(contour_max_count*size of(CvPoint));
　　cvCvtSeqToArray(contour2,pointarray,CV_WHOLE_SEQ);
　　//查找最长和最短轴
　　for (i=0;idist_center[i])
　　{
　　r_short=dist_center[i];
　　num_short=i;
　　}
　　}
　　axisratio=r_long/r_short;//计算长短轴的比值
　　axis_angle=(float)abs(num_long-num_short)/contour_ max_count;//计算长短轴对应的角度
　　dist_avar=dist_sum/contour_max_count;//均值
　　float vari_temp=0;
　　for (i=0;iheight;i++)
　　{
　　for (int j=0;jwidth;j++)
　　{
　　if (((uchar*)(dst->imageData+i*dst->widthStep))[j]==255 )
　　{
　　sum_pix++;
　　int temp=((uchar*)(src->imageData+i*src->widthStep))[j];
　　sum_gray+=temp;
　　ffreq[temp]++;
　　}
　　}
　　}
　　for (i=0;i 0)
　　{
　　// 计算图像熵
　　shan_maxcontour -= ffreq[i] * log(ffreq[i]) / log(2.0);
　　}
　　}
　　averagegray=(double)sum_gray/sum_pix;//均值
　　int grey_s=0;
　　for (i=0;iheight;i++)
　　{
　　for (int j=0;jwidth;j++)
　　{
　　if (((uchar*)(dst->imageData+i*dst->widthStep))[j]==255 )
　　{
　　grey_s+=(((uchar*)(src->imageData+i*src->widthStep)) [j]-averagegray)*(((uchar*)(src->imageData+i*src->wid thStep))[j]-averagegray);
　　}
　　}
　　}
　　variance_gray=sqrt(grey_s/(sum_pix+eps)); //方差
　　}
　　else
　　{
　　return 0;
　　}
　　CString aa(Path);
　　aa+="-1.jpg";//保存图像为新的路径，但不影响原先的bmp图像
　　cvSaveImage(aa,src);
　　cvReleaseMemStorage(&storage2);
　　cvReleaseMemStorage(&storage);*/