对cvMatchContourTree的测试代码(还是没找出关系) ， 第八章11题

代码如下：

#include <cv.h>#include <cxcore.h>#include <highgui.h>#include <cvaux.h>#include <iostream>#include <string>#include <strstream>#include <iomanip>#include <cmath>using namespace std;static void imRotate(IplImage *src,IplImage *dst,double angle , CvPoint2D32f center,double scale = 1 );//缩放图像int main(){IplImage *src = cvLoadImage("f:\\images\\test3.bmp",CV_LOAD_IMAGE_GRAYSCALE);IplImage *src_cpy = cvCloneImage(src);CvMemStorage *storage = cvCreateMemStorage();cvThreshold(src,src,50,255,CV_THRESH_BINARY);CvSeq *contour_src = NULL;cvFindContours(src_cpy,storage,&contour_src);//此函数会改变输入图像CvContourTree *tree_src = cvCreateContourTree(contour_src,storage,0);cvReleaseImage(&src_cpy);CvSeq *contour = NULL;CvRNG rng = cvRNG(-1);bool flag = false; //控制scale和angle哪个变量变化while(1){//控制其中一个变化，另一个固定//与原图像比较double scale = 1 , angle = 0;if(flag){scale = cvRandReal(&rng)*1.3;if(scale < 0.5)continue;}elseangle = cvRandReal(&rng)*360;time_t StartTime = clock();IplImage *dst = cvCreateImage(cvSize(src->width,src->height),8,1);imRotate(src,dst,angle,cvPoint2D32f(src->width/2,src->height/2),scale);cvShowImage("dst",dst);cvThreshold(dst,dst,100,255,CV_THRESH_BINARY);cvFindContours(dst,storage,&contour);CvContourTree *tree = cvCreateContourTree(contour,storage,0);double result = cvMatchContourTrees(tree_src,tree,CV_CONTOUR_TREES_MATCH_I1,0);time_t EndTime = clock();cout<<"angle = "<<setprecision(3)<<angle<<" scale = "<<scale<<"\t";cout<<"result = "<<setprecision(3)<<result<<"\ttime = "<<double(EndTime - StartTime)/CLOCKS_PER_SEC<<endl;cvReleaseImage(&dst);cvClearSeq(contour);int key = cvWaitKey();if(key == 27)break;else if(key == 'n' || key == 'N')flag = !flag;}cvReleaseMemStorage(&storage);cvReleaseImage(&src);}void imRotate(IplImage *src, IplImage *dst,double angle , CvPoint2D32f center,double scale){assert(src->width == dst->width && src->height == dst->height &&src->depth == dst->depth &&src->nChannels == dst->nChannels);CvMat *mapMatrix = cvCreateMat(2,3,CV_32FC1);IplImage *tmp = cvCreateImage(cvGetSize(src),src->depth,src->nChannels);cv2DRotationMatrix(center,angle,scale,mapMatrix); //旋转缩放为仿射变换，此处求变换矩阵cvWarpAffine(src,(src->imageData == dst->imageData) ? tmp : dst,mapMatrix); //此函数不能src == dstif(src->imageData == dst->imageData)cvCopyImage(tmp,dst);cvReleaseImage(&tmp);}

使用图像为:

在固定其中一个变量的情况下做匹配，发现得到的值差异很大，不知道这个轮廓树比较得到的结果有什么用处，请各位路过的大侠指教!

试验结果：

从结果中看几乎无关联 , 我又修改了程序

使无旋转时，图像的大小随着scale变化（上面的程序图像大小不变)，即用cvResize缩放,仍然是一样的情况　，看不出关联