直方图均衡化

dwm88888888

已于 2023-07-23 15:44:50 修改

阅读量66

点赞数

分类专栏： Opencv 文章标签： opencv 人工智能计算机视觉

于 2023-07-21 10:00:36 首次发布

本文链接：https://blog.csdn.net/dwm88888888/article/details/131845548

版权

Opencv 专栏收录该内容

9 篇文章 0 订阅

订阅专栏

该代码示例展示了如何使用OpenCV库在C++中对灰度图像进行直方图均衡化处理，然后计算并绘制直方图。通过`cvtColor`函数将BGR图像转换为灰度，接着使用`equalizeHist`进行直方图均衡化，最后利用`calcHist`计算直方图并用自定义函数`drawHist`绘制直方图。

摘要由CSDN通过智能技术生成

#include "iostream"
#include "opencv2/opencv.hpp"
using namespace std;
using namespace cv;


void drowHist(Mat &mat, string name)
{
	//归一化
	Mat dst;
	double alpha = 1;
	double beta = 0;
	int norm_type = NORM_INF;
	int  dtype = -1;
	normalize(mat, dst, alpha, beta, norm_type, dtype, Mat());

	//绘制直方图
	int hist_height = 300;
	int hist_width = 512;
	int unit_width = 2;
	float rowvalue = 0;
	Point zs, yx, ys;//直方图左上方的点，右上方的点，右下方的点
	Mat picture = Mat::zeros(Size(hist_width, hist_height), CV_8UC3);
	for (int i = 0; i < dst.rows; i++)
	{
		rowvalue = dst.at<float>(i);
		int heighttemp = hist_height - cvRound(hist_height * rowvalue);//翻转
		zs = Point(i*unit_width, heighttemp);
		ys = Point(i*unit_width + 1, heighttemp);
		yx = Point(i*unit_width + 1, hist_height);
		//直方图
		rectangle(picture, zs, yx, Scalar(255, 0, 0));
		//线
		line(picture, zs, ys, Scalar(0, 0, 255));
	}

	imshow(name, picture);
}

int main()
{
	Mat  img, gray, equalizehist;

	img = imread("bao.jpg");
	cvtColor(img, gray, COLOR_BGR2GRAY);
	namedWindow("gray", WINDOW_NORMAL);
	namedWindow("equalizehist", WINDOW_NORMAL);
	equalizeHist(gray, equalizehist);

	imshow("gray", gray);
	imshow("equalizehist", equalizehist);

	int nimages = 1;//图片数量
	const int channels[] = { 0 };//需要统计直方图的第几通道
	Mat hist ;//直方图计算的输出值 Mat
	Mat hist2;//直方图计算的输出值 Mat
	int dims = 1;
	const int histSize[] = { 256 }; //直方图每一个维度划分的柱条的数目
	float pranges[] = { 0, 255 };//取值区间
	const float* ranges[] = { pranges };
	calcHist(&gray, nimages, channels, Mat(), hist, dims, histSize, ranges);
	calcHist(&equalizehist, nimages, channels, Mat(), hist2, dims, histSize, ranges);

	drowHist(hist, "temp1");
	drowHist(hist2, "temp2");

	waitKey(0);

	return 1;
}