距离变换

#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;
}