距离变换

#include <cv.h>#include <highgui.h>int main(){    char* filename = "cs.jpg";    IplImage* src_image = cvLoadImage(filename,1);    if(!src_image)        return -1;    cvNamedWindow("src");    CvSize size = cvGetSize(src_image);    IplImage* gray_image = cvCreateImage(size,8,1);    cvCvtColor(src_image,gray_image,CV_BGR2GRAY);    IplImage* dist_image = cvCreateImage(size,32,1);    IplImage* bi_src = cvCreateImage(size,8,1);    IplImage* dist8u_image = cvCreateImage(size,8,1);    IplImage* bi_dist = cvCreateImage(size,8,1);    //原图像二值化    cvThreshold(gray_image,bi_src,100,255,CV_THRESH_BINARY);    //距离变换    cvDistTransform(bi_src,dist_image,CV_DIST_L2,3,0,0);    //找最大值    double max;    cvMinMaxLoc(dist_image,0,&max,0,0);    cvCvtScale(dist_image,dist8u_image,255./max);    cvShowImage("dist8u_image",dist8u_image);//这一句我自己加的，显示出我们想要看到的距离变换之后的图像    //对距离图像二值化，去除手指部分   cvThreshold(dist8u_image,bi_dist,80,255,CV_THRESH_BINARY);    //求重心    float s=0.0, x=0.0, y=0.0;     uchar* data = (uchar*)bi_dist->imageData;    int step = bi_dist->widthStep;    for(int h=0;h<bi_dist->height;h++)        for(int w=0;w<bi_dist->width;w++)            if(255 == data[step*h+w])            {                x += w;                y += h;                s++;            }    if(s>0)    {        x = x/s;        y = y/s;    }    CvPoint pos = cvPoint((int)x,(int)y);    cvCircle(src_image,pos,3,CV_RGB(255,0,0),1,CV_AA);    cvShowImage("src",src_image);    cvWaitKey(-1);    cvDestroyWindow("src");    cvReleaseImage(&src_image);    cvReleaseImage(&gray_image);    cvReleaseImage(&bi_src);    cvReleaseImage(&dist8u_image);    cvReleaseImage(&bi_dist);    return 0;}