图像的直方图显示与直方图均衡处理

最新推荐文章于 2023-02-21 17:52:15 发布

Emiyasstar__

最新推荐文章于 2023-02-21 17:52:15 发布

阅读量5k

点赞数

分类专栏：图像处理文章标签： arrays float matlab 图像处理文档图形

本文链接：https://blog.csdn.net/Emiyasstar__/article/details/6918024

版权

图像处理专栏收录该内容

3 篇文章 0 订阅

订阅专栏

一.直方图概念

直方图其实就是一副图像的灰度级统计，设一副图像的灰度级为[0,L-1],则可以建立一个二维坐标系（x,y),x取值[0,L-1],y取值为相应灰度级的像素数量。如果是归一化直方图，则y代表该灰度级在图像中的出现概率，且。

二.基于opencv的直方图显示

相较于matlab ,opencv的直方图显示要麻烦的多，因为它生成的只是一个数组，需要自己手动画出来，这里介绍基于c++ API的显示方法:

核心函数:

void calcHist(const Mat*arrays, intnarrays, const int*channels, InputArraymask, OutputArrayhist, intdims, const int*histSize, const float**ranges, booluniform=true, boolaccumulate=false )

常用参数简介，详细请参考opencv文档：

arrays:需要显示的图像数组；

narrays:输入的图像数目，一般都是取1吧;

channels:输入图像的通道数目，这里比较麻烦，详细参见官方文档，一般取1，也就是单通道灰度图；

mask:不明，通常取mat();

hist:输出的直方图数组；

dims:直方图维度；

histSize：直方图尺寸，如果是8位深就取256;

ranges;显示范围

这里给出一个简单的直方图显示类，将常用参数封装了一下：

class Histogram
{
public:
	Histogram()
	{
       histSize[0] = 256;
	   hrangee[0] = 0.0;
       hrangee[1] = 255.0;
	   ranges[0]= hrangee;//8位深度，256灰度级
	   channels=0;
	}
	Mat calcHistogram(const Mat& imageGray)//创建直方图显示
	{   
		Mat hist,histImage;	
	    calcHist(&imageGray, 1, 0, Mat(), hist, 1, histSize, ranges);//构建直方图矩阵
		histImage=showHistogarm(hist);//图形显示
		return histImage;

	}
	Mat calcEqualizeHistogram(const Mat& imageGray,Mat& imageGrayEqualize)//直方图均衡化
	{
		Mat hist,histImage;	
		equalizeHist(imageGray,imageGrayEqualize);//直方图均衡
	    calcHist(&imageGrayEqualize, 1, 0, Mat(), hist, 1, histSize, ranges);
		histImage=showHistogarm(hist);
		return histImage;

	}
	~Histogram()
	{

	}
private:
	Mat showHistogarm( Mat& hist)//直方图的图形显示
	{
		Mat histImage = Mat::ones(200, 320, CV_8U)*255;//初始化图像
		normalize(hist, hist, 0, histImage.rows, CV_MINMAX, CV_32F);//归一化
		histImage = Scalar::all(255);

		int binW = cvRound((double)histImage.cols/histSize[0]);//各灰度级显示的宽度

		for( int i = 0; i < histSize[0]; i++ )//灰度256级，依次画出每一级的矩形表示
            rectangle( histImage, Point(i*binW, histImage.rows),
                   Point((i+1)*binW, histImage.rows - cvRound(hist.at<float>(i))),
                   Scalar::all(0), -1, 8, 0 );
		return histImage;
	}
private:
	int histSize[1]; //直方图尺寸
	float hrangee[2];
    const float *ranges[1];//灰度级
	int channels;//通道数
};

三.直方图的简单应用，直方图均衡:

直方图均衡可以说是最常用也是最基础的图像增强方法了,大致原理既是将图像的灰度范围动态展宽，是所有图像灰度级实现等概率分布。

假设r为输入图像灰度，s为增强图像灰度值，P(r),P(s)分别为相应的概率密度，设s,r之间存在灰度变换公式：