Question

帮忙调个OPENCV程序，这个程序是用直方图的方法求两张图片的相似度， cvCalcHist过程总有问题

代码如下
#include <highgui.h>
#include <cv.h>
#include <iostream>
double compare_image(IplImage *image1,IplImage *image2)
{
CvHistogram *hist1=0;
CvHistogram *hist2=0;//声明两张直方图
int dim=3;//直方图维数
int size[]={50,50,50};//直方图的维数尺寸
float r_a[]={0,255};
float g_a[]={0,255};
float b_a[]={0,255};
float *range[]={r_a,g_a,b_a};
hist1=cvCreateHist(dim,size,CV_HIST_ARRAY,range,1);
hist2=cvCreateHist(dim,size,CV_HIST_ARRAY,range,1);
cvCalcHist(&image1,hist1);
cvCalcHist(&image2,hist2);
double temp=cvCompareHist(hist1,hist2,CV_COMP_INTERSECT);
cvReleaseHist(&hist1);
cvReleaseHist(&hist2);
return temp;
}
int main()
{
IplImage *image1,*image2;
image1=cvLoadImage("D:/1.jpg");
image2=cvLoadImage("D:/2.jpg");
std::cout<<compare_image(image1,image2)<<"\n";
cvNamedWindow("image1");
cvNamedWindow("image2");
cvShowImage("image1",image1);
cvShowImage("image2",image2);
cvWaitKey();
cvDestroyWindow("image1");
cvDestroyWindow("image2");
cvReleaseImage(&image1);
cvReleaseImage(&image2);
return 0;
}
我主要这个程序编译没问题，运行就会有内存出错，我需要比较两张彩色图片，想法是分别建立两张图的直方图，用直方图的相似度cvCompareHist(hist1,hist2,CV_COMP_INTERSECT)来代表两帧画面的相似度。。。

Answer 1

//计算直方图，并用图片显示出来直方图结果图像；
double compare_image(IplImage *image1,IplImage *image2)
{
int size=256;
float range[]={0,255};
float* ranges[]={range};
double binnum[256]={0};

CvHistogram* hist1=cvCreateHist(1,&size,CV_HIST_ARRAY,ranges,1);
CvHistogram* hist2=cvCreateHist(1,&size,CV_HIST_ARRAY,ranges,1);
cvCalcHist(&image1,hist1,0,NULL);
cvCalcHist(&image2,hist2,0,NULL);

float max1=0;
float max2=0;
cvGetMinMaxHistValue(hist1,NULL,&max1,NULL,NULL);
cvGetMinMaxHistValue(hist2,NULL,&max2,NULL,NULL);

IplImage* dst=cvCreateImage(cvSize(400,300),8,3);
cvSet(dst,cvScalarAll(255),0);

double bin_width=(double)dst->width/size;
double bin_unith1=(double)dst->height/max1;
double bin_unith2=(double)dst->height/max2;
for(int i=0;i<size;i++)
{
CvPoint p0=cvPoint(i*bin_width,dst->height);
CvPoint p1=cvPoint((i+1)*bin_width,dst->height-cvGetReal1D(hist1->bins,i)*bin_unith1);
cvRectangle(dst,p0,p1,cvScalar(0,0,255),-1,8,0); //用红色显示直方图

CvPoint p2=cvPoint((i+1)*bin_width,dst->height-cvGetReal1D(hist2->bins,i)*bin_unith2);
cvRectangle(dst,p0,p2,cvScalar(255,0,0),-1,8,0); //用绿色显示直方图
}

double temp=cvCompareHist(hist1,hist2,CV_COMP_INTERSECT);//比较直方图
cvShowImage("gray image histogram",dst);
cvReleaseHist(&hist1);
cvReleaseHist(&hist2);

return temp;
}

追问

我不要显示直方图 我只要能比较就行了。。。

追答

double temp=cvCompareHist(hist1,hist2,CV_COMP_INTERSECT);//比较直方图
没出结果吗？

追问

你给我一个 你能最终运行出结果版本看看。。。