图像直方图的原理及代码实现

目录

一、图像直方图的基本概念

二、代码部分

1.代码解读

2.完整代码


一、图像直方图的基本概念

图像直方图(Histogram)是一种统计报告图,它通过一系列高度不等的纵向条纹或线段表示图像数据的分布情况。通常,横轴表示数据类型,纵轴表示分布情况。对于灰度图像,直方图显示的是不同灰度级别的像素数量分布;对于彩色图像,直方图可以按通道分别显示R、G、B三个颜色通道的像素数量分布,也可以合并显示彩色图像的像素数量分布。

构建直方图需要将图像数据的范围分段,然后计算每个段落的像素数量。对于灰度图像,通常将像素灰度值从0到255分成256个段,然后统计每个段落的像素数量;对于彩色图像,可以按通道(R、G、B)分别进行统计,也可以将三个通道的像素数量合并后再进行统计。

二、代码部分

1.代码解读

cv2.calcHist(images,channels,mask,histSize,ranges)

images: 原图像图像格式为 uint8 或 float32。当传入函数时应 用中括号 [] 括来例如[img]
channels: 表示传入的图像通道数。如果输入图像是灰度图它的值就是 [0]。如果是彩色图像 的传入的参数可以是 [0][1][2] 它们分别对应着 BGR。
mask: 掩模图像。统计整幅图像的直方图就把它为None。但是如果你想统计图像某一部分的直方图,你就制作一个掩模图像并使用它。
histSize:BINS的数目。也需用中括号括来10
BINS :上面的直方图显示了每个像素值的像素数,即从0到255。即您需要256个值才能显示上述直方图。
但是请考虑一下,如果您不需要单独查找所有像素值的像素数,而是在像素值间隔内查找像素数,
该怎么办?例如,您需要找到介于 0 到 15 之间的像素数,然后是 16 到 31、...、240 到 255。
您只需要 16 个值来表示直方图。
因此,只需将整个直方图拆分为 16 个子部分,每个子部分的值就是其中所有像素计数的总和。这每个子部分都称为"BIN"。在第一种情况下,条柱数为256(每个像素一个),而在第二种情况下,它只有16。BINS 在 OpenCV 文档中由术语histSize表示。
ranges: 像素值范围常为 [0 256]

2.完整代码
phone = cv2.imread('mx.png',cv2.IMREAD_GRAYSCALE)
phone_hist = cv2.calcHist([phone],[0],None,[16],[0,256])
a = phone.ravel()
plt.hist(a, bins=256)#numpy中的ravel将数组多维度拉成一维数组
plt.show()

plt.plot(phone_hist)#使用calcHist的值绘制曲线图
plt.show()


#绘制三个通道的直方图曲线
img = cv2.imread('phone.png')
color = ('b', 'g', 'r')
for i, col in enumerate(color):
    histr = cv2.calcHist([img], [i], None, [256], [0, 256])
    plt.plot(histr, color=col)
plt.show()

运行结果:

绘制了mx.png的一维灰度图直方图

绘制了之前计算的灰度图像的直方图:

定义了颜色通道后绘制直方图

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
图像直方图均衡化是一种用于增强图像对比度的技术,它通过对图像的像素值进行调整,使得图像中的亮度分布更均匀。下面是一个简单的C++实现原理: 1. 首先,计算图像的灰度直方图,即统计每个灰度级别的像素数量。 2. 根据直方图,计算每个灰度级别的累积分布函数(CDF),表示灰度级别小于等于当前级别的像素所占的比例。 3. 根据CDF,计算每个灰度级别的映射表,将原始图像的每个像素值映射到新的像素值上。映射表的计算公式为: new_value = round((L-1) * CDF(old_value)) 其中,L表示灰度级别的数量,CDF(old_value)表示旧像素值对应的累积分布函数值。 4. 遍历原始图像的每个像素,根据映射表将其像素值替换为新的像素值。 下面是一个简单的示例代码: ```cpp #include <iostream> #include <opencv2/opencv.hpp> using namespace cv; using namespace std; void histogramEqualization(Mat& image) { // 计算直方图 Mat hist; int histSize = 256; float range[] = {0, 256}; const float* histRange = {range}; calcHist(&image, 1, 0, Mat(), hist, 1, &histSize, &histRange, true, false); // 计算累积分布函数 Mat cdf; hist.copyTo(cdf); for (int i = 1; i < histSize; i++) { cdf.at<float>(i) += cdf.at<float>(i - 1); } cdf /= image.total(); // 计算映射表 Mat lut(1, histSize, CV_8UC1); for (int i = 0; i < histSize; i++) { lut.at<uchar>(i) = saturate_cast<uchar>((histSize - 1) * cdf.at<float>(i)); } // 应用映射表 LUT(image, lut, image); } int main() { // 读取图像 Mat image = imread("input.jpg", IMREAD_GRAYSCALE); if (image.empty()) { cout << "Failed to read image" << endl; return -1; } // 直方图均衡化 histogramEqualization(image); // 显示结果 imshow("Equalized Image", image); waitKey(0); return 0; } ``` 请注意,这只是一个简单的示例代码,实际应用中可能需要处理更多的细节,例如处理彩色图像、处理边界情况等。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值