直方图和模块匹配

#include "cv.h"
#include "highgui.h"
#include <cxcore.h>

int main(int agrc, char** agrv)
{
 IplImage* src_image = cvLoadImage("2.jpg");
 IplImage* dst_image = cvCloneImage(src_image);
 IplImage* gray_src_image = cvCreateImage(cvGetSize(src_image), IPL_DEPTH_8U, 1);
 IplImage* gray_dst_image = cvCloneImage(gray_src_image);

 cvCvtColor(src_image, gray_src_image, CV_RGB2GRAY);
 cvNamedWindow("src", 1);
 cvNamedWindow("graysrc", 1);
 cvShowImage("src", src_image);
 cvShowImage("graysrc", gray_src_image);
//
//
//计算和显示直方图
/* 
 IplImage* hsv = cvCreateImage(cvGetSize(src_image), 8, 3);
 cvCvtColor(src_image, hsv, CV_RGB2HSV);

 IplImage* h_plane = cvCreateImage(cvGetSize(src_image), 8, 1);
 IplImage* s_plane = cvCreateImage(cvGetSize(src_image), 8, 1);
 IplImage* v_plane = cvCreateImage(cvGetSize(src_image), 8, 1);
 IplImage* planes[] = {h_plane, s_plane};
 cvCvtPixToPlane(hsv, h_plane, s_plane, v_plane, 0);

 int h_bins = 30, s_bins = 32;
 CvHistogram* hist;

 int hist_size[] = {h_bins, s_bins};
 float h_ranges[] = {0, 180};
 float s_rangs[] = {0, 255};
 float* ranges[] = {h_ranges, s_rangs};

 hist = cvCreateHist(2, hist_size, CV_HIST_ARRAY, ranges, 1);

 cvCalcHist(planes, hist, 0, 0);
 cvNormalizeHist(hist, 1.0);
 
 int scale = 10;
 IplImage* hist_img = cvCreateImage(cvSize(h_bins * scale, s_bins * scale), 8, 3);
 cvZero(hist_img);

 float max_value = 0;
 cvGetMinMaxHistValue(hist, 0, &max_value, 0, 0);
 for (int h = 0; h < h_bins; h++)
 {
  for (int s = 0; s < s_bins; s++)
  {
   float bin_val =cvQueryHistValue_2D(hist, h, s);
   int intensity = cvRound(bin_val * 255 / max_value);
    cvRectangle(
     hist_img,
     cvPoint(h * scale, s * scale),
     cvPoint((h + 1) * scale - 1, (s + 1) * scale - 1),
     CV_RGB(intensity, intensity, intensity),
     CV_FILLED);
  }
 }

 cvNamedWindow("H-S Histogram", 1);
 cvShowImage("H-S Histogram", hist_img);
*/
//
//
//模板匹配
/*
 IplImage* temp = cvLoadImage("2.2.jpg");
 IplImage* result = cvCreateImage(cvSize(src_image->width - temp->width + 1, src_image->height - temp->height + 1), IPL_DEPTH_32F, 1);

 cvMatchTemplate(src_image, temp, result, CV_TM_SQDIFF);

 cvNamedWindow("SQDIFF", 1);
 cvShowImage("SQDIFF", result);

 
 cvMatchTemplate(src_image, temp, result, CV_TM_CCOEFF);
 
 cvNamedWindow("CCOEFF", 1);
 cvShowImage("CCOEFF", result);

 
 cvMatchTemplate(src_image, temp, result, CV_TM_CCORR);
 
 cvNamedWindow("CCORR", 1);
 cvShowImage("CCORR", result);

 
 cvMatchTemplate(src_image, temp, result, CV_TM_SQDIFF_NORMED);
 
 cvNamedWindow("SQDIFF_NORMED", 1);
 cvShowImage("SQDIFF_NORMED", result);

 
 cvMatchTemplate(src_image, temp, result, CV_TM_CCORR_NORMED);
 
 cvNamedWindow("CCORR_NORMED", 1);
 cvShowImage("CCORR_NORMED", result);

 
 cvMatchTemplate(src_image, temp, result, CV_TM_CCOEFF_NORMED);
 
 cvNamedWindow("CCOEFF_NORMED", 1);
 cvShowImage("CCOEFF_NORMED", result);
*/
//
//

 cvReleaseImage(&src_image);
 cvReleaseImage(&dst_image);
 cvReleaseImage(&gray_src_image);
 cvReleaseImage(&gray_dst_image);
 int c = cvWaitKey(0);
 if (c == 27)
 {
  return -1;
 }

 return 0;

}