cvCalcBackProjectPatch() 基于块的反向投影

使用在输入图像平面上的滑动块来设置目标图像上相应的像素（块中心）

测试图像：

find.jpg

test.jpg

程序：

#include <iostream>
#include <cv.h>
#include <highgui.h>
using namespace std;

void GetHSV(const IplImage *image , IplImage **h,IplImage **s,IplImage **v)
{
	IplImage *hsv = cvCreateImage(cvGetSize(image),8,3);
	cvCvtColor(image,hsv,CV_BGR2HSV);

	if((h != NULL) && (*h == NULL))
		*h = cvCreateImage(cvGetSize(image),8,1);
	if((s != NULL) && (*s == NULL))
		*s = cvCreateImage(cvGetSize(image),8,1);
	if((v != NULL) && (*v == NULL))
		*v = cvCreateImage(cvGetSize(image),8,1);

	cvSplit(hsv,*h,(s == NULL)?NULL:*s,(v==NULL)?NULL:*v,NULL);
	cvReleaseImage(&hsv);
}

int main()
{
	IplImage *src = cvLoadImage("find.jpg");//输入待寻找块
	IplImage *h_src = NULL ,*s_src = NULL;
	GetHSV(src,&h_src,&s_src,NULL);//获取H、S通道
	IplImage *images[] = {h_src,s_src};//生成H、S双通道图像


	CvHistogram *hist_src;//创建直方图
	int dims = 2;//直方图维数为2( H和S )
	int size[] = {10,20};//H通道(0-180)和S(0-256)通道的直方图横轴单位矩形的大小(时间效率和精确度之间的取舍)
	float range_h[]={0,180},range_s[]={0,256};//H通道灰度范围和S通道的灰度范围
	float *ranges[]={range_h,range_s};//二维直方图的灰度范围
	hist_src = cvCreateHist(dims,size,CV_HIST_ARRAY,ranges);
	cvCalcHist(images,hist_src);//计算二维直方图
	cvNormalizeHist(hist_src,1);//归一化

	IplImage *dst=cvLoadImage("test.jpg");//待搜索图像
	IplImage *h_dst=NULL,*s_dst =NULL;
	GetHSV(dst,&h_dst,&s_dst,NULL);
	images[0]=h_dst,images[1]=s_dst; 

	CvSize patch_size=cvSize(src->width,src->height);
	IplImage *result=cvCreateImage(cvSize(h_dst->width-patch_size.width+1,h_dst->height-patch_size.height+1),IPL_DEPTH_32F,1);
	cvCalcBackProjectPatch(images,result,patch_size,hist_src,CV_COMP_CORREL,1);
	cvShowImage("result",result);

	CvPoint max_location;
	cvMinMaxLoc(result,NULL,NULL,NULL,&max_location,NULL);
	max_location.x += cvRound(patch_size.width/2);
	max_location.y += cvRound(patch_size.height/2);

	cvCircle(dst,max_location,3,CV_RGB(255,0,0),-1);
	cvShowImage("dst",dst);

	cvWaitKey();
	cvReleaseImage(&src);
	cvReleaseImage(&dst);
	cvReleaseImage(&h_src);
	cvReleaseImage(&h_dst);
	cvReleaseImage(&s_dst);
	cvReleaseImage(&s_src);
	cvReleaseHist(&hist_src);
	cvReleaseImage(&result);
	cvDestroyAllWindows();
}

运行结果（大概运行了10秒）：