opencv函数源码之直方图规定化（1）

由于在opencv中没有找到相关的源代码，于是网上搜索到别人的写的代码，将他拿来研究一下，地址：点击打开链接

// do the hist match histv[] is the vector that you need match （直方图与想要匹配的直方图进行匹配）
void HistMatch(IplImage *histimg,double histv[])  //histv[]需要匹配的直方图
{
	int bins = 256;
	int sizes[] = {bins};
	CvHistogram *hist = cvCreateHist(1,sizes,CV_HIST_ARRAY);
	cvCalcHist(&histimg,hist);//计算直方图
	
	cvNormalizeHist(hist,1);//直方图归一化
	
	double val_1 = 0.0;
	double val_2 = 0.0;
	uchar T[256] = {0};
	double S[256] = {0};
    double G[256] = {0};

	for (int index = 0; index<256; ++index)
	{
		val_1 += cvQueryHistValue_1D(hist,index);//返回一个特定元素的单通道阵列；可以安全的使用单通道；见附录一
    val_2 += histv[index];

		G[index] = val_2;
		S[index] = val_1;
	}
	
	double min_val = 0.0;
	int PG = 0;
//比较原图的累计直方图与需要匹配的累计直方图，如果之间的差值小于min_val，
//则认定两者匹配，将匹配的值保存存到T[ ]中
	for ( int i = 0; i<256; ++i)
	{
		min_val = 1.0;
		for(int j = 0;j<256; ++j)
		{
		    if( (G[j] - S[i]) < min_val && (G[j] - S[i]) >= 0)
			{
               min_val = (G[j] - S[i]);
			   PG = j;
			}
		
		}
		T[i] = (uchar)PG;
	}
    
	uchar *p = NULL;
	for (int x = 0; x<histimg->height;++x)
	{ 
		p = (uchar*)(histimg->imageData + histimg->widthStep*x);
		for (int y = 0; y<histimg->width;++y)
		{
			p[y] = T[p[y]];
		}
		
   }
 //===================附录一：cvQueryHistValue_1D=====================//
 #define cvQueryHistValue_1D( hist, idx0 ) \
  ((float)cvGetReal1D( (hist)->bins, (idx0)))