灰度直方图巴士距离排序.docVIP

灰度直方图巴士距离排序 #includecvaux.h #includecv.h #includehighgui.h #includeiostream #includemath.h #includevector #define SAMPLES 200//样本个数 using namespace std; using namespace cv; vectorint bin; vectorint bin_cmp;//存放 每个图像的直方图信息 void covar(Mat frame_get);//求取已知图片的协方差矩阵 void covar_cmp(Mat frame_get);//求取待选图片的协方差矩阵 double compare_hist();//求取直方图的相似性 int main() { vectordouble bc_dis; vectorint flag; vectorvector int flaga; vectorvector double bc_disa; Mat img,img_cmp; int sum=0; double sim; char imgname[100]; char imgname_cmp[100]; for(int num=0;numSAMPLES;num++) { sprintf_s(imgname,D:\\viper\\a\\%d.bmp,num); img=imread(imgname); for(int numc=0;numcSAMPLES;numc++) { sprintf_s(imgname_cmp,D:\\viper\\b\\%d.bmp,numc); img_cmp=imread(imgname_cmp); covar(img);//求原图的协方差矩阵 covar_cmp(img_cmp);//求比较图像的协方差矩阵 sim=compare_hist(); bc_dis.push_back(sim); } for(int i=0;ibc_dis.size();i++) { flag.push_back(i); } for(int i=0;ibc_dis.size();i++) { double temp; int tem; for(int j=0;jbc_dis.size();j++) { if (bc_dis[i]bc_dis[j]) { temp=bc_dis[i]; bc_dis[i]=bc_dis[j]; bc_dis[j]=temp; tem=flag[i]; flag[i]=flag[j]; flag[j]=tem; } } } coutflag[0] flag[1] flag[2] flag[3] flag[4] flag[5] flag[6] flag[6] endl; //coutbc_dis[0] bc_dis[1] bc_dis[2]endl; bc_disa.push_back(bc_dis); bc_dis.clear(); flaga.push_back(flag); flag.clear(); } for(int i=0;iflaga.size();i++) { if(flaga[i][0]==i) sum++; } coutsumendl; imshow(img,img); waitKey(); return 0; } void covar(Mat frame_get) { Mat img; Size dsize=Size(64,128); Mat frame=Mat(dsize,CV_8UC3); resize(frame_get,frame,dsize); cvtColor(frame,img,CV_BGR2GRAY); int sum[16]={0}; // 数组sum是每个量化等级的像素点的个数统计之和 int x[16]={0},y[16]={0}; // 数组x,y存放的是每个量化等级的像素坐标均值(x,y) double cov1[16]={0.0},cov2[16]={0.0