图像处理开发需求、图像处理接私活挣零花钱,请加微信/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、之间是空白的,这样柱状图的效果就出来了。即每两根柱状之间通过空白进行了分隔和区别。