OpenCV2.3的cvCalcHist函数有问题?255级值总为0,索性自己写一个直方图计算函数,附源码!

欢迎大家加入图像识别技术交流群：271891601，另外，特别欢迎成都从事图像识别工作的朋友交流，我的QQ号248787278

-------------------------------------------

我在写直方图规定化的代码过程中，发现OpenCV自带的cvCalcHist函数计算出的直方图的第255分量总是为0，测试了几张图都是这样，代码如下：

#include <opencv2/opencv.hpp>
#include <opencv2/legacy/compat.hpp>
#include <fstream>
using namespace std;
#pragma comment(linker, "/subsystem:\"windows\" /entry:\"mainCRTStartup\"")  

CvHistogram* CreateGrayImageHist(IplImage **ppImage)
{
    int nHistSize = 256;  //这句话表明有256个灰度级
    float fRange[] = {0, 255};  //灰度级的范围是0到255
    float *pfRanges[] = {fRange};   //下一句使用的cvCreateHist函数要求第四个参数是指针的指针,
								    //所以要作这样的处理,至于为什么要用指针的指针,我也不清楚
    CvHistogram *pcvHistogram = cvCreateHist(1, &nHistSize, CV_HIST_ARRAY, pfRanges); //在使用直方图计算函数cvCalcHist前,
																					  //需要先创建cvCalcHist的第二个参数需要数据类型
    cvCalcHist(ppImage, pcvHistogram);//直方图计算函数cvCalcHist,注意第一个参数要求是指针的指针,具体为什么要这样,我也不清楚
				       //至于怎么取用存放在pcvHistogram中的直方图数据,看完下面的main函数就知道了
    return pcvHistogram;
}  

int main()

{
	IplImage *a = cvLoadImage("coins.png", CV_LOAD_IMAGE_UNCHANGED);

	CvHistogram *a_hist = CreateGrayImageHist(&a);// CreateGrayImageHist是自己写的函数哦

	int index;
	double watch_a_hist[256];
	for (index = 0; index<256; index++)
	{
		 watch_a_hist[index] = cvQueryHistValue_1D(a_hist,index);
	}

return 0;
}

和MATLAB运行结果的对比截图如下：

于是自己写了一个直方图计算函数，经测试，无误 ，代码如下：

void mycvCalcHist(IplImage *img,double out_hist[256])
{
	int i=0, j=0;
	double temp1=0;
	int temp2=0;
    const int hist_sz = 256;//0到255，一共256个灰度值
    double hist[hist_sz];
    memset(hist, 0, sizeof(hist));
    for( i = 0; i < img->height; i++ )
    {
        for( j = 0; j < img->width; j++ )
		{	temp1=cvGet2D(img,i,j).val[0];
			temp2=int(temp1);//作类型转换
            hist[temp2]++; //这里实现了hist中存储各灰度值出现的次数
		}
    }  

	memcpy(out_hist,hist, sizeof(hist)); //肯定有人要问,为啥不用数组名作为参数传递从而改变实参数组的值
										 //这种方法一般情况下都可以,我也测试了,然而这里就是不行,我估计与
										 //memset(hist, 0, sizeof(hist));这句语句有关

}

int main()
{
	IplImage *a = cvLoadImage("coins.png", CV_LOAD_IMAGE_UNCHANGED);  

	double a_hist[256];

	mycvCalcHist(a,a_hist);

return 0;
}

补充说明一下为啥用double类型的数组，因为直方图数据经常要参与浮点数运算，所以干脆用这个！