《opencv实用探索·十六》opencv直方图计算calcHist函数解析

直方图理解:
(对于8位灰度图像亮度/灰度为(0-255),12位灰度图像亮度/灰度为(0-4095))
以8位图像为例,亮度分为0到255共256个数值,数值越大,代表的亮度越高。其中0代表纯黑色的最暗区域,255表示最亮的纯白色。亮度分为了5个区域,分别是黑色,阴影,中间调,高光和白色。直方图用来统计一幅图像某个亮度像素数量。
在这里插入图片描述

直方图函数解析:

void calcHist( const Mat* images, int nimages,const int* channels, InputArray mask,OutputArray hist, int dims, const int* histSize,const float** ranges, bool uniform=true, bool accumulate=false );

images:输入的图像的指针,可以是多幅图像,所有的图像必须有同样的深度(CV_8U or CV_32F)。同时一副图像可以有多个channes。
nimages:输入图像的个数
channels:需要统计直方图的第几通道。用来计算直方图的channes的数组。比如输入是2副图像,第一副图像有0,1,2共三个channel,第二幅图像只有0一个channel,那么输入就一共有4个channes,如果int channels[3] = {3, 2, 0},那么就表示是使用第二副图像的第一个通道和第一副图像的第2和第0个通道来计算直方图。
mask:掩膜,如果是空矩阵则表示图像中所有位置的像素都计入直方图中,如果mask不为空,那么它必须是一个8位(CV_8U)的数组,并且与输入图像尺寸相同,找到mask非0像素区域,把这些区域对应到输入图像区域,然后把这些区域参与直方图计算。
hist:输出的直方图数组
dims:需要统计直方图通道的个数,通常为1
histSize:指的是直方图分成多少个区间,就是 bin的个数。
const float** ranges: 统计像素值得区间。
uniform=true::是否对得到的直方图数组进行归一化处理
accumulate=false:在多个图像时,是否累计计算像素值得个数

几个例子理解直方图参数含义:
首先我们对一个mat图像做几个设定如下:
(1)8位单通道灰度图,那么这个图像灰度范围为0-255
(2)mat图像总像素个数为10000
(3)灰度范围在0-127时像素个数为8000, 灰度范围128-255时像素个数为2000
(3)灰度范围在0-84时像素个数为3000, 灰度范围85-170时像素个数为6000,灰度范围171-255时像素个数为1000

int main()
{Mat mat = cv::imread("1.tif", -1);int max_value = 0;int channels[] = { 0 };const int histSize[] = { 256 };float range[] = { 0, 256};const float* ranges[] = { range };cv::Mat hist;cv::calcHist(&mat, 1, channels, cv::Mat(), hist, 1, histSize, ranges, true, false);int sum = 0;for (int i = 0; i < 255; i++){sum += hist.at<float>(i);}return 0;
}

上面代码中可知像素值范围range为0-255(不包含256),histSize为256相当于把rang(0-255)分为256份。for循环对histSize总份数进行遍历,相当于遍历0-255每个灰度,这样可以拿到每个灰度下的像素个数,最终sum总和等于10000。

int main()
{Mat mat = cv::imread("1.tif", -1);int max_value = 0;int channels[] = { 0 };const int histSize[] = { 2 };float range[] = { 0, 256};const float* ranges[] = { range };cv::Mat hist;cv::calcHist(&mat, 1, channels, cv::Mat(), hist, 1, histSize, ranges, true, false);int sum = 0;for (int i = 0; i < 2; i++){sum += hist.at<float>(i);}return 0;
}

上面代码中range灰度范围为0-255,histSize把灰度范围分为两份,第一份灰度范围0-127,第二份灰度范围128-255
for循环中遍历总份数,即遍历两份,第一份可以拿到灰度范围0-127的所有像素之和,即8000,第二份可以拿到128-255范围所有像素之和,即2000

int main()
{Mat mat = cv::imread("1.tif", -1);int max_value = 0;int channels[] = { 0 };const int histSize[] = { 3 };float range[] = { 0, 256};const float* ranges[] = { range };cv::Mat hist;cv::calcHist(&mat, 1, channels, cv::Mat(), hist, 1, histSize, ranges, true, false);int sum = 0;for (int i = 0; i < 3; i++){sum += hist.at<float>(i);}return 0;
}

上面代码中range范围0-255,histSize把range总共分为三份,即0-84,85-170,171-255
for循环遍历histSize,在第一份中可以拿到灰度范围在0-84的总像素个数,即3000,在第二份中可以拿到灰度范围在85-170的总像素个数,即6000,第三份可以拿到171-255范围的总像素个数,即1000

下面是直方图的一个具体应用,8位三通道计算直方图案例:

#include <opencv2/opencv.hpp>
#include <iostream>using namespace cv;
using namespace std;int main() {// 读取图像Mat image = imread("your_image.jpg");if (image.empty()) {cerr << "Error: Could not read the image." << endl;return -1;}// 设置直方图参数int histSize = 256;  // 直方图的 bin 数量float range[] = {0, 256};  // 像素值范围const float* histRange = {range};bool uniform = true;  // 直方图是否均匀bool accumulate = false;  // 是否累积直方图// 计算蓝色通道的直方图int channelB = 0;  // 蓝色通道索引Mat histB;calcHist(&image, 1, &channelB, Mat(), histB, 1, &histSize, &histRange, uniform, accumulate);// 计算绿色通道的直方图int channelG = 1;  // 绿色通道索引Mat histG;calcHist(&image, 1, &channelG, Mat(), histG, 1, &histSize, &histRange, uniform, accumulate);// 计算红色通道的直方图int channelR = 2;  // 红色通道索引Mat histR;calcHist(&image, 1, &channelR, Mat(), histR, 1, &histSize, &histRange, uniform, accumulate);// 绘制直方图int histWidth = 512;int histHeight = 400;int binWidth = cvRound((double) histWidth / histSize);Mat histImage(histHeight, histWidth, CV_8UC3, Scalar(0, 0, 0));// 归一化直方图数据normalize(histB, histB, 0, histImage.rows, NORM_MINMAX, -1, Mat());normalize(histG, histG, 0, histImage.rows, NORM_MINMAX, -1, Mat());normalize(histR, histR, 0, histImage.rows, NORM_MINMAX, -1, Mat());// 绘制直方图for (int i = 1; i < histSize; ++i) {line(histImage, Point(binWidth * (i - 1), cvRound(histB.at<float>(i - 1))),Point(binWidth * (i), cvRound(histB.at<float>(i))),Scalar(255, 0, 0), 2, 8, 0);line(histImage, Point(binWidth * (i - 1), cvRound(histG.at<float>(i - 1))),Point(binWidth * (i), cvRound(histG.at<float>(i))),Scalar(0, 255, 0), 2, 8, 0);line(histImage, Point(binWidth * (i - 1),cvRound(histR.at<float>(i - 1))),Point(binWidth * (i), cvRound(histR.at<float>(i))),Scalar(0, 0, 255), 2, 8, 0);}// 显示原始图像和直方图imshow("Original Image", image);imshow("Histogram", histImage);waitKey(0);return 0;
}

在这里插入图片描述

也可也遍历拿到BGR每个通道下0-255每个灰度的像素个数
比如遍历B通道

for(int i = 0; i < 256; i++)
{float count = cvRound(histB.at<float>(i - 1);
}

0-256为每个灰度,for循环遍历每个灰度,并拿到该灰度下像素的个数,比如i=100,count=20表示,灰度为100的像素有20个。

在这里插入图片描述

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/265605.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Linux_CentOS_7.9 VNC安装卸载以及相关配置开机自启动服务简易记录

VNC安装卸载以及相关配置开机自启动服务&#xff1a; 查看环境&#xff1a;&#xff08;yum镜像源配置可以参考我之前文章里面有详细参考http://t.csdnimg.cn/mzGoI&#xff09; [rootorcl238 ~]# rpm -qa | grep vnc ##查看系统现有VNC软件版本 gtk-vnc2-0.7.0-3.el7.x86…

〖Python网络爬虫实战㊷〗- 极验滑块介绍(四)

订阅&#xff1a;新手可以订阅我的其他专栏。免费阶段订阅量1000 python项目实战 Python编程基础教程系列&#xff08;零基础小白搬砖逆袭) 说明&#xff1a;本专栏持续更新中&#xff0c;订阅本专栏前必读关于专栏〖Python网络爬虫实战〗转为付费专栏的订阅说明作者&#xff1…

在前端开发中,什么是SEO(Search Engine Optimization)?如何优化网站的SEO?

聚沙成塔每天进步一点点 ⭐ 专栏简介 前端入门之旅&#xff1a;探索Web开发的奇妙世界 欢迎来到前端入门之旅&#xff01;感兴趣的可以订阅本专栏哦&#xff01;这个专栏是为那些对Web开发感兴趣、刚刚踏入前端领域的朋友们量身打造的。无论你是完全的新手还是有一些基础的开发…

泰裤辣!!!我又学到了监控服务器的各种方法!

好耶&#xff01;又能学习到监控服务器的各种方法了。像是利用linux命令去监控的方法就有好几种哦 方法一&#xff1a;动态监控top命令 方法二&#xff1a;内存监控free -m 方法三&#xff1a;磁盘iostat -d 2 3 还有利用jmeter插件去监控服务器哦我们可以使用ServerAgent插件…

为什么字符串一旦创建就不可以改变了么?

今天看到一个String和StringBuilder的题&#xff0c;很是恼怒&#xff01; 总结下吧。 判断string类创建后可以被修改&#xff08;&#xff09; 理解&#xff1a; 字符串的值一旦被创建&#xff0c;就不能改变指的是字符串的内容不能发生改变。而不是字符串的引用不能发生改…

数据结构与算法—递归及其应用(八皇后,小球出迷宫)

递归 文章目录 递归1.递归应用场景2.递归的概念3.递归调用机制4.递归能解决什么问题5.递归需要遵守的重要规则6.递归-迷宫、八皇后问题6.1 迷宫问题6.2 八皇后问题 1.递归应用场景 迷宫问题(回溯)&#xff0c;递归(Recursion) 2.递归的概念 简单的说&#xff1a;递归就是方法自…

AMEYA360:大唐恩智浦荣获 2023芯向亦庄 “汽车芯片50强”

2023年11月28日&#xff0c;由北京市科学技术委员会和北京市经济和信息化局指导、北京经济技术开发区管理委员会主办、盖世汽车协办的“芯向亦庄”汽车芯片大赛在北京亦庄成功闭幕。 在本次大赛中 大唐恩智浦的 电池管理芯片DNB1168 (应用于新能源汽车BMS系统) 凭卓越的性能及高…

AidLux:手机/平板上的Linux环境与AI开发利器

AidLux是一个基于ARM架构的跨生态&#xff08;鸿蒙/AndroidLinux&#xff09;一站式智能物联网&#xff08;AIoT&#xff09;应用开发和部署平台&#xff0c;正受到越来越多开发者和用户的青睐。既可以作为手机/平板上的一个Linux环境使用&#xff0c;也可以作为AI开发利器以发…

HCIA-H12-811题目解析(5)

1、【单选题】 以下关于Hybrid端口说法正确的有&#xff1f; 2、【单选题】使用命令"vlan batch 10 20"和"valn batch 10 to 20"&#xff0c;分别能创建的vlan数量是&#xff1f;&#xff08;&#xff09; 3、【单选题】二层ACL的编号范围是&#xff1f;…

厘米级高精度定位系统为什么更倾向于UWB技术?

超宽带&#xff08;Ultra Wide-Band&#xff0c;UWB&#xff09;是一种新型的无线通信技术&#xff0c;根据通信委员会的规范&#xff0c;UWB的工作频带为3.1~10.6GHz&#xff0c;系统-10dB带宽与系统中心频率之比大于20%或系统带宽至少为500MHz。 UWB信号的发生可通过发射时间…

Java预科知识

以下内容是根据狂神的Java说、chatgpt和csdn相关博客&#xff0c;结合自己的理解完成的。 Java了解 基于Java 开发了巨多的平台&#xff0c;系统&#xff0c;工具 构建工具&#xff1a; Ant, Maven, Jekins应用服务器&#xff1a;Tomcat, Jetty, Jboss, Websphere, weblogic…

网络攻击(二)--情报搜集阶段

4.1. 概述 在情报收集阶段&#xff0c;你需要采用各种可能的方法来收集将要攻击的客户组织的所有信息&#xff0c;包括使用社交网络、Google Hacking技术、目标系统踩点等等。 而作为渗透测试者&#xff0c;你最为重要的一项技能就是对目标系统的探查能力&#xff0c;包括获知…