1. 简介

利用OpenCV读取图像,转换为灰度图像,绘制该灰度图像直方图。点击直方图,控制台输出该灰度级像素个数。

2. 原理

(1) 实现原理较为简单,主要利用了OpenCV读取图像,并转换为灰度图像;

    srcImg = imread(" ......");    // “....” 代表图像地址
if (srcImg.empty()) {
return -;
}
imshow(WINDOW_SRCIMG, srcImg);
Mat grayImg;
cvtColor(srcImg, grayImg, CV_BGR2GRAY);
imshow(WINDOW_GRAYIMG, grayImg);

(2) 利用Mat类新建一个固定分辨率的画布,统计处于每一灰度级像素个数,在该画布上绘制灰度直方图。

    int nRows = ,nCols=;
Mat g_dstImg(nRows,nCols, CV_8UC1, Scalar::all()); // 新建画布

同时避免画布中该灰度级太高而超出画布范围,在本程序中采用了等比例缩小的方法。

    int MaxCount = arrayMax(grayLevel,);//寻找在处于某一灰度级中个数最多的像素个数
yscaleRate = double(nRows)/MaxCount ;//y缩放比例
double xscaleRate = nCols / ;//x缩放比例

3. 实施细节

// 灰度直方图.cpp : 定义控制台应用程序的入口点。
// Topic: 绘制灰度图的直方图; 转载请注明出处:Chen_HW (https://www.cnblogs.com/chen-hw/p/11668119.html)
// Env: VS2015 + Debug x64 + OpenCV3.4.1
// Date:2018.11.22 by Chen_HW
#include "stdafx.h"
#include <opencv2/opencv.hpp>
#include <iostream>
#define WINDOW_SRCIMG "【源图】"
#define WINDOW_GRAYIMG "【灰度图】"
#define WINDOW_HIST "【直方图】"
using namespace std;
using namespace cv;
int arrayMax(int g_arr[], int num);
Mat drawHist(Mat &g_srcImg);
void MouseHandle(int event, int x, int y, int flags, void *param);
Point clickPoint,displayPoint;
bool downFlag = false;
Mat srcImg;
double yscaleRate = ; int main()
{
system("color 3f");
srcImg = imread(" ......"); // “....” 代表图像地址
if (srcImg.empty()) {
return -;
}
imshow(WINDOW_SRCIMG, srcImg);
Mat grayImg;
cvtColor(srcImg, grayImg, CV_BGR2GRAY);
imshow(WINDOW_GRAYIMG, grayImg);
Mat histImg = drawHist(grayImg);
imshow(WINDOW_HIST, histImg);
setMouseCallback(WINDOW_HIST, MouseHandle, (void *)&histImg);//(void *)&srcImg传递给void *param
waitKey();
return ;
}
Mat drawHist(Mat &g_srcImg) {
int nRows = ,nCols=;
Mat g_dstImg(nRows,nCols, CV_8UC1, Scalar::all()); // 新建画布
int grayLevel[] = {};
for (int i = ; i < g_srcImg.rows; ++i) {
for (int j = ; j < g_srcImg.cols; ++j) {
grayLevel[(int)g_srcImg.at<uchar>(i, j)]++;
}
}
int MaxCount = arrayMax(grayLevel,);//寻找在处于某一灰度级中个数最多的像素个数
yscaleRate = double(nRows)/MaxCount ;//y缩放比例
double xscaleRate = nCols / ;//x缩放比例
int yAxis[], xAxis[];
for (int m = ; m < ; m++) {
yAxis[m] = int(grayLevel[m]*yscaleRate);
}
//绘制直方图
for (int n = ; n < ; n++) {
if (n == ) {
xAxis[n + ] = xAxis[n];
}
rectangle(g_dstImg, Point(xAxis[n], yAxis[n]), Point( xAxis[n+], ), Scalar(), -);
//-1表示填充矩形框;正值表示不填充矩形框,更方便观察灰度级像素个数的分布;
//但因为后面需要用到填充的情况,故设置成填充状态 } return g_dstImg;
}
//num:数组元素个数;g_max:返回最大值
int arrayMax(int g_arr[],int num) {
int g_max = ;
int i = ;
while (i < num) {
if ( g_arr[i]>= g_max) {
g_max = g_arr[i];
i++;
}
else {
i++;
}
} return g_max;
}
//event 鼠标事件(如按下鼠标左键、左键抬起、鼠标移动等) x、y 鼠标坐标
void MouseHandle(int event, int x, int y, int flags, void *param) {
Mat &g_srcImg = *(Mat *)(param);
Mat g_tempImg = g_srcImg.clone();
int nCount=;//该列白色像素点个数
int nLevelCount = ;
char text[];//存储文本信息
float g_rate;//该列像素点占总像素点个数的比例
imshow(WINDOW_HIST, g_srcImg);
switch (event) {
case EVENT_LBUTTONDOWN:
clickPoint.x = x;
clickPoint.y = y;
downFlag = true;
break;
default:
break;
}
if (downFlag) {
//displayPoint.x = clickPoint.x;
for (int i = ; i < g_srcImg.rows; ++i) {
if (int(g_srcImg.at<uchar>(i, x)) == ) {
++nCount;
}
}
nLevelCount = (nCount / yscaleRate);
g_rate = float(nCount/ yscaleRate) / (srcImg.rows*srcImg.cols);//此处计算的比例不是精确的
cout << "该灰度级像素点个数: " << nLevelCount << ";占总像素个数的比例: "<<g_rate << endl;
sprintf_s(text, "Rate:%f",g_rate );//该灰度级像素个数占总个数的比例;
putText(g_tempImg, text, clickPoint,FONT_HERSHEY_PLAIN,,Scalar(,,));
imshow(WINDOW_HIST, g_tempImg);
downFlag = false;
}
}

4. 结果

  结果如下图所示,点击右侧每一级灰度直方图,在控制台中会输出该灰度级像素个数,并显示占总像素比例。

【图像处理】利用C++编写函数,绘制灰度图像直方图的更多相关文章

  1. 使用matplotlib中的bar函数绘制柱状图

    使用柱状图显示三日电影的票房信息 要显示的数据为2018年12月7日-9日四场电影的票房信息 四场电影分别为:无名之辈,狗十三,毒液:知名守卫者,憨豆特工3 2018年12月7日四场电影票房分别为:[ ...

  2. C++:利用如下公式,编写函数计算∏的值,直到最后一项的绝对值小于e,主程序接收从键盘输入的e,输出∏的值(保留5位小数)。 ∏/4 = 1-1/3+1/5-1/7...

    利用如下公式,编写函数计算∏的值,直到最后一项的绝对值小于e,主程序接收从键盘输入的e,输出∏的值(保留5位小数). ∏/4 = 1-1/3+1/5-1/7... #include <iostr ...

  3. Python 图像处理 OpenCV (16):图像直方图

    前文传送门: 「Python 图像处理 OpenCV (1):入门」 「Python 图像处理 OpenCV (2):像素处理与 Numpy 操作以及 Matplotlib 显示图像」 「Python ...

  4. paper 2:图像处理常用的Matlab函数汇总

    一 图像的读写 1 imread imread函数用于读入各种图像文件,如:a=imread('e:\w01.tif') 注:计算机E盘上要有w01相应的.tif文件. 2 imwrite imwri ...

  5. 利用java编写的盲注脚本

    之前在网上见到一个盲注的题目,正好闲来无事,便用java写了个盲注脚本,并记录下过程中的坑 题目源码: <?php header("Content-Type: text/html;ch ...

  6. 品味性能之道<九>:利用Loadrunner编写socket性能测试脚本简述

            一.概述         Loadrunner拥有极为丰富的工具箱,供予我们制造出各种奇妙魔法的能力.其中就有此次要讨论的socket套接字操作.     二.socket概述     ...

  7. 向日葵远程RCE漏洞分析及漏洞利用脚本编写

    0x00 漏洞概述 向日葵是一款免费的,集远程控制电脑.手机.远程桌面连接.远程开机.远程管理.支持内网穿透等功能的一体化远程控制管理软件.如果想要手机远控电脑,或者电脑远控手机可以利用向日葵:如果是 ...

  8. 利用c++编写bp神经网络实现手写数字识别详解

    利用c++编写bp神经网络实现手写数字识别 写在前面 从大一入学开始,本菜菜就一直想学习一下神经网络算法,但由于时间和资源所限,一直未展开比较透彻的学习.大二下人工智能课的修习,给了我一个学习的契机. ...

  9. iOS基本动画/关键帧动画/利用缓动函数实现物理动画效果

    先说下基本动画部分 基本动画部分比较简单, 但能实现的动画效果也很局限 使用方法大致为: #1. 创建原始UI或者画面 #2. 创建CABasicAnimation实例, 并设置keypart/dur ...

随机推荐

  1. Entity Framework使用EntityState和Attach来保存数据变化以及更新实体的个别字段

    在使用Entity Framework作为ORM来存取数据的过程中,最常规的操作就是对数据对象的更新.本文将会包含如何Attach Entity到一个数据Context中,以及如何使用EntitySt ...

  2. CF #636 (Div. 3) 对应题号CF1343

    unrated 选手悠闲做题,然后只做出四个滚蛋了 符合 div3 一贯风格,没啥难算法 E最后就要调出来了,但还是赛后才A的 CF1343A Candies 传送门 找到一个 \(x\),使得存在一 ...

  3. CF1324D Pair of Topics

    好像题解里都是树状数组(起码我翻到的是 说一种cdq分治的(这应该算是cdq分治了 用cdq比较简单,所以可以作为一个练手题 cdq分治其实是一种模糊的思想,处理\([l,r]\)区间内,有多少\(( ...

  4. windows中配置安装mysql数据库

    MySql 是一种免费的关系型数据库,相较于 MsSqlServer 和 Oracle 比较轻量化,安装也很简单,而且免费不需要的版权费用,个人认为一般的小项目采用还是比较合适的,当然也有部分数据量很 ...

  5. 《C程序设计语言》 练习1-22

    问题描述 练习1-22 编写一个程序,把较长的输入行“折”成短一些的两行或者多行,折行的位置在输入行的第N列之前的最后一个非空格之后.要保持程序能够智能地处理输入行很长以及在制定的列前没有空格或者制表 ...

  6. ipc$链接批量爆破

    写了一个bat bat code: @echo off cls echo Useage: % ip.txt pass.txt ) do ( FOR /F ) do ( echo net use \\% ...

  7. POJ3734(矩阵快速幂)

    \(假设现在到第i个积木\) \(红绿恰都是偶数a种方案,恰都是奇数为b种方案,一奇一偶为c种方案\) \(由此考虑i+1个积木的情况\) Ⅰ.一奇一偶的方案 \(如果第i层恰是奇数的情况,那么本次只 ...

  8. EChart将timeline图标设置为加号

    一 要实现的效果 二 官方文档 看一下官方文档上,替换symbol的描述: timeline.symbol string [ default: 'emptyCircle' ] timeline标记的图 ...

  9. ASP.Net Core中使用Swagger

    我们先简单介绍下什么是Swagger,主要是用来干嘛?? 在Swagger诞生之前,我们通常在开发接口的过程中,需要前后端共同维护一个接口文档,然后大家按照接口文档的规范进行对接.接口文档俨然成了接口 ...

  10. kubernetes部署redis主从高可用集群

    1.redis主从高可用集群结构 2.k8s部署有状态的服务选择 对于K8S集群有状态的服务,我们可以选择deployment和statefulset statefulset service& ...