OpenCV1.x下的灰度直方图绘制及显示的源码,带详细注释

图像处理开发需求、图像处理接私活挣零花钱，请加微信/QQ 2487872782
图像处理开发资料、图像处理技术交流请加QQ群，群号 271891601

注：本文的源码参考下面这篇博文：

https://blog.csdn.net/morewindows/article/details/8364656

下面的代码是“OpenCV1.x下的灰度直方图绘制及显示的源码”，自己写上了详细注释，以供大家阅读时参考。

代码中用到的图片的下载链接：

https://pan.baidu.com/s/1mJK80q1i4MyukKOt05ysYQ  提取码：shdw 

//写这段代码注释者的微信/QQ 2487872782
//有问题可以联系作者交流
//图像处理技术交流QQ群 271891601

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


void FillWhite(IplImage *pImage)  
{  
    cvRectangle(pImage, cvPoint(0, 0), cvPoint(pImage->width, pImage->height), CV_RGB(255, 255, 255), CV_FILLED);  
}  

// 计算绘制直方图所需要的数据
CvHistogram* CreateGrayImageHist(IplImage **ppImage)  
{  
    int nHistSize = 256;  //这句话表明有256个统计区间
    float fRange[] = {0, 255};  //灰度级的范围是0到255    
    float *pfRanges[] = {fRange};   //下一句使用的cvCreateHist函数要求第四个参数是指针的指针,之所以要这样定义,是因为函数cvCreateHist()可计算多个通道的直方图数据
    CvHistogram *pcvHistogram = cvCreateHist(1, &nHistSize, CV_HIST_ARRAY, pfRanges); //调用函数cvCalcHist()前要先定义类型为CvHistogram的变量
    cvCalcHist(ppImage, pcvHistogram);//调用OpenCV1.X的直方图数据统计函数cvCreateHist()
    return pcvHistogram;  
}  

// 绘制直方图
IplImage* CreateHisogramImage(int nImageWidth, int nScale, int nImageHeight, CvHistogram *pcvHistogram)  
{  
    IplImage *pHistImage = cvCreateImage(cvSize(nImageWidth * nScale, nImageHeight), IPL_DEPTH_8U, 1);  //创建直方图图像
    FillWhite(pHistImage);  //把创建好的直方图图像全部填充为白色,FillWhite也是自己写的函数哦
  
    //统计直方图中的最大直方块，以便于绘图时进行归一化处理
    float fMaxHistValue = 0;  
    cvGetMinMaxHistValue(pcvHistogram, NULL, &fMaxHistValue, NULL, NULL);  
  
    //分别将每个直方块的值绘制到图中  
    int i;  
    for(i = 0; i < nImageWidth; i++)  
    {  
        float fHistValue = cvQueryHistValue_1D(pcvHistogram, i); //灰度级i的直方图统计值 
        int nRealHeight = cvRound((fHistValue / fMaxHistValue) * nImageHeight);  //灰度级i在直方图中的实际高度 
        cvRectangle(pHistImage,  
            cvPoint(i * nScale, nImageHeight - 1),  //需要绘制的直方图矩形的一个顶点
													//要明白这几行代码,自己举一个实际的例子并仔细分析就清楚了，我的博文中也配有一个例子。
													//为什么这里nImageHeight要减去1，是因为为了保证即使某个灰度值区间的像素个数为0，也能绘制出相应的矩形，如果不减1，当某个灰度值区间的像素个数为0，则其柱状矩形的两个纵坐标都为nImageHeight，函数cvRectangle()就会报错了。
            cvPoint((i + 1) * nScale - 1, nImageHeight - nRealHeight),  //需要绘制的直方图矩形的另一个顶点
            cvScalar(i, 0, 0, 0),   //这是填充的颜色
            CV_FILLED //这个参数表示矩形需要填充，即矩行是实心的
        );   
    }  
    return pHistImage;  

}  
int main( int argc, char** argv )  
{     
    const char *pstrWindowsSrcTitle = "原图";  
    const char *pstrWindowsGrayTitle = "灰度图";  
    const char *pstrWindowsHistTitle = "直方图";  
  
    // 从文件中加载原图  
    IplImage *pSrcImage = cvLoadImage("yang_mi.jpg", CV_LOAD_IMAGE_UNCHANGED);  
    IplImage *pGrayImage = cvCreateImage(cvGetSize(pSrcImage), IPL_DEPTH_8U, 1);  
    
	// 转换为灰度图
    cvCvtColor(pSrcImage, pGrayImage, CV_BGR2GRAY);  
  
    // 计算绘制直方图所需要的数据
    CvHistogram *pcvHistogram = CreateGrayImageHist(&pGrayImage);  // CreateGrayImageHist是自己写的函数
      
    // 绘制直方图
    int nHistImageWidth = 255;   //直方图的宽度
    int nHistImageHeight = 150;  //直方图的高度
    int nScale = 4;  //这个值的意义你要看了函数CreateHisogramImage()的源码后才能了解:
	                 //先透露一下,实际上是把最后Windows绘制出的直方图图像的窗口的横轴按比例拉伸了
					 //避免Windows绘制出的直方图图像窗口看着太小,可以自己改下值看下结果,或者看我下面给的截图
    IplImage *pHistImage = CreateHisogramImage(nHistImageWidth, nScale, nHistImageHeight, pcvHistogram);  //调用自己写的绘制直方图的函数
  
    // 显示  
    cvNamedWindow(pstrWindowsSrcTitle, CV_WINDOW_AUTOSIZE);  
    cvNamedWindow(pstrWindowsGrayTitle, CV_WINDOW_AUTOSIZE);  
    cvNamedWindow(pstrWindowsHistTitle, CV_WINDOW_AUTOSIZE);  
    cvShowImage(pstrWindowsSrcTitle, pSrcImage);  
    cvShowImage(pstrWindowsGrayTitle, pGrayImage);  
    cvShowImage(pstrWindowsHistTitle, pHistImage);  
  
  
    cvWaitKey(0);  
  
    cvReleaseHist(&pcvHistogram);  
  
    cvDestroyWindow(pstrWindowsSrcTitle);  
    cvDestroyWindow(pstrWindowsGrayTitle);  
    cvDestroyWindow(pstrWindowsHistTitle);  
    cvReleaseImage(&pSrcImage);  
    cvReleaseImage(&pGrayImage);  
    cvReleaseImage(&pHistImage);  
    return 0;  
}  

下面的截图体现nScale取不同值时的效果:

注释中绘制直方图时说自己可以举例理解的地方举例如下:

i=0、1、2时，nRealHeight的值分别为130、110、90。

设nScale=3，则有：

当i=0时，需要绘制的矩形(柱状图)的两个坐标分别是：

(0*3,150-1)和((0+1)*3-1,150-130)

↓↓↓↓↓↓↓

(0,149)和(2,20)

当i=1时，需要绘制的矩形(柱状图)的两个坐标分别是：

(1*3,150-1)和((1+1)*3-1,150-110)

↓↓↓↓↓↓↓

(3,149)和(5,40)

当i=2时，需要绘制的矩形(柱状图)的两个坐标分别是：

(2*3,150-1)，((2+1)*3-1,150-90)

↓↓↓↓↓↓↓

(6,149)和(8,60)

手绘上述过程如下图所示：

从图中，我们可以看出：2和3、5和6、之间是空白的，这样柱状图的效果就出来了。即每两根柱状之间通过空白进行了分隔和区别。