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

news2024/10/6 12:26:19

直方图理解:
(对于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.coloradmin.cn/o/1302115.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

XCube——用于超高分辨率 3D 形状和场景的生成模型!

他们的方法在稀疏体素网格的层次结构上训练潜在扩散模型的层次结构。他们在稀疏结构 VAE 的潜在空间上进行扩散&#xff0c;它为层次结构的每个级别学习紧凑的潜在表示。 XCube 是稀疏体素层次上的分层潜在扩散模型&#xff0c;即从粗到细的 3D 稀疏体素网格序列&#xff0c;使…

OpenWRT搭建本地web站点并结合内网穿透实现公网远程访问

文章目录 前言1. 检查uhttpd安装2. 部署web站点3. 安装cpolar内网穿透4. 配置远程访问地址5. 配置固定远程地址 前言 uhttpd 是 OpenWrt/LuCI 开发者从零开始编写的 Web 服务器&#xff0c;目的是成为优秀稳定的、适合嵌入式设备的轻量级任务的 HTTP 服务器&#xff0c;并且和…

【高数:3 无穷小与无穷大】

【高数&#xff1a;3 无穷小与无穷大】 1 无穷小与无穷大2 极限运算法则3 极限存在原则4 趋于无穷小的比较 参考书籍&#xff1a;毕文斌, 毛悦悦. Python漫游数学王国[M]. 北京&#xff1a;清华大学出版社&#xff0c;2022. 1 无穷小与无穷大 无穷大在sympy中用两个字母o表示无…

数据结构--二叉树

目录 1.二叉树链式结构的实现 1.1 前置说明 1.2 二叉树的遍历 1.2.1 前序、中序以及后序遍历 1.2.2 层序遍历及判断是否为完全二叉树 1.3 节点个数&#xff0c;叶子节点个数&#xff0c;第k层节点个数以及高度等 1.4 二叉树的创建和销毁 1.二叉树链式结构的实现 1.1 前置说…

使用人工智能优化常见业务流程

在现代商业环境中&#xff0c;人工智能(AI)正在改变企业的运营方式。将人工智能集成到业务流程中可以提高效率和准确性&#xff0c;从而节省大量时间和成本。 这使员工能够专注于更具战略性的任务。人工智能在商业中的应用范围从自动化日常任务到提供高级分析&#xff0c;以做…

Citespace、vosviewer、R语言的文献计量学可视化分析

文献计量学是指用数学和统计学的方法&#xff0c;定量地分析一切知识载体的交叉科学。它是集数学、统计学、文献学为一体&#xff0c;注重量化的综合性知识体系。特别是&#xff0c;信息可视化技术手段和方法的运用&#xff0c;可直观的展示主题的研究发展历程、研究现状、研究…

《杂文选刊》明年休刊之随笔

笔者今天偶然发现&#xff0c;网民对于近期登上社交网站热搜榜的一条新闻&#xff0c;既感兴趣又觉迷惑&#xff0c;因此关注度显得较高。 荣登社交网站热搜榜的这条新闻是&#xff1a;12月4日发出“休刊启事 ”&#xff0c;在宣布“《 杂文选刊》2024年1月1日起休刊”的同时&…

和鲸科技荣获第三届光合组织解决方案大赛集智赛道优秀奖

2023年11月28日&#xff0c;历经数月的“第三届光合组织解决方案大赛”落下帷幕&#xff0c;获奖榜单正式出炉。 本次大赛中&#xff0c;上海和今信息科技有限公司&#xff08;简称“和鲸科技”&#xff09;凭借多年深耕数据智能领域&#xff0c;提供关键基础设施催生人工智能…

解决HTTP 429错误的Scrapy中间件配置

引言 在进行网络数据抓取时&#xff0c;经常会遇到HTTP 429错误&#xff0c;表示请求速率已超出API限制。为避免封禁或限制访问&#xff0c;需要调整Scrapy的请求速率&#xff0c;以在不触发HTTP 429错误的情况下完成数据抓取。针对这一问题&#xff0c;可使用Scrapy的AutoThr…

Win11预览体验计划显示Your PC does not meet the minimum hardware requirements...的解决方案

某一天你心血来潮&#xff0c;打算参与Win11 预览体验计划&#xff0c;但体验计划页面却显示“Your PC does not meet the minimum hardware requirements for Windows11…”。 一种解决思路&#xff1a; 去以下网页下载Offline Insider Enroll软件&#xff0c;管理员权限运行…

安全开发:身份认证方案之 Google 身份验证器和基于时间的一次性密码 TOTP 算法

参考资料在文末注明&#xff0c;如本文有错漏欢迎评论区指出&#x1f44f; 目前很多应用都逐步采用了双因子认证或者说MFA认证方案&#xff0c;因此本文介绍一下背后的机制和TOTP算法原理。使用TOTP算法&#xff0c;只要满足两个条件&#xff1a;1&#xff09;基于相同的密钥&…

RocketMQ源码

RocketMQ的核心三流程 启动流程 RocketMQ服务端由两部分组成NameServer和Broker&#xff0c;NameServer是服务的注册中心&#xff0c;Broker会把自己的地址注册到NameServer&#xff0c;生产者和消费者启动的时候会先从NameServer获取Broker的地址&#xff0c;再去从Broker发…

BEUEC品牌比例阀放大器厂家

HE-SP1/HE-SP2/HE-SP2-U/HE-MPS1/HE-MPS2/HE-MAPQ-V/HE-MAPQ-K/HE-MPT2-Y/HE-MPT2-U比例放大器适配控制各种不带位置反馈比例阀&#xff1b; BEUEC品牌比例放大器控制如博世力士乐&#xff08;Bosch Rexroth&#xff09;、伊顿威格士&#xff08;EATON Vickers&#xff09;、油…

shiro入门demo

搭建springboot项目&#xff0c;引入以下依赖&#xff1a; <dependencies><dependency><groupId>org.springframework.boot</groupId><artifactId>spring-boot-starter-web</artifactId></dependency><!--单元测试--><depe…

[CIKM 2023] 基于会话推荐的双通道多稀疏图注意网络

Bi-channel Multiple Sparse Graph Attention Networks for Session-based Recommendation GitHub - QEpiphany/MSGAT: Multiple Sparse Graph Attention Networks for Session-based Recommendation 摘要 基于会话的推荐&#xff08;SBR&#xff09;最近受到了极大的关注&…

安装程序无法自动安装Virtual Machine Communication Interface Sockets(VSock)驱动程序

环境情况&#xff1a; 物理机win10系统 虚拟机windowserver08系统 vmware 16.0的版本 问题触发&#xff1a; 在虚拟机win7系统上安装vmware tools出现提示&#xff0c;报错信息“安装程序无法自动安装Virtual Machine Communication Interface Sockets&#xff08;VSock&a…

iPhone 数据恢复:iMyFone D-Back iOS

iMyFone D-Back iOS 最佳 iPhone 数据恢复&#xff0c;最好的 iPhone 数据恢复软件&#xff0c;恢复成功率最高。 直接从iOS设备、iTunes/iCloud/第三方程序备份快速恢复数据。 有选择地恢复已删除的照片、WhatsApp、消息和 18 多种其他数据类型。 仅通过 iCloud 帐户访问即可从…

jmeter 压测需要的部分配置

修改jmeter 目录的bin目录下的jmeter.properties文件 解除KeepAlive设置 修改接口的高级中的实现和超时 解除httpclient4.retrycount前的注释符并将0修改为1 即修改为&#xff1a;httpclient4.retrycount1 解除httpclient4.idletimeout前的注释符并修改为合适间隔 即修改为…

java--Arrays类、自定义排序规则Comparable、自定义比较器Comparator

1.Arrays 用来操作数组的一个工具类。 2.Arrays类提供的常见方法 3.如果数组中存储的是对象&#xff0c;如何排序 方法一&#xff1a;让该对象的类实现Comparable(比较规则)接口&#xff0c;然后重写CompareTo方法&#xff0c;自己来定制比较规则。 方法二&#xff1a;使用下…

​ 轻量应用服务器:亚马逊云科技打造全球领先的云计算解决方案

随着“第四次工业革命”的爆炸式发展&#xff0c;众多企业都将自己的业务与迅速发展的应用开发和网站建设领域高度绑定。而对于众多有上云需求的企业和个人用户来说&#xff0c;选择一款自己的服务器配置就成为了一项至关重要的任务。而随着需求端的不断扩大&#xff0c;云服务…