OpenCV:实现灰度直方图和单通道直方图拉伸
原文链接:http://blog.csdn.net/xiaowei_cqu/article/details/7600666
本文略有修改,如有疑问或者版权问题,请移步原作者或者告知本人。
直接贴代码:
//计算直方图cv::MatND colorWish::getHist(const cv::Mat &image, int histSize){cv::Mat im;if (image.channels() == 3)cv::cvtColor(image, im, CV_RGB2GRAY, 0);else im = image;float r[2]; r[0] = 0; r[1] = 255;const float *ranges[1];ranges[0] = r;cv::MatND hist;//float range[] = { 0,255 }; //灰度级的范围 //float* ranges[] = { range };const int channels[1] = { 0 };cv::calcHist(&im, 1, channels, cv::Mat(), hist, 1, &histSize, ranges);return hist;}//直方图拉伸cv::Mat colorWish::stretch(const cv::Mat &image, int minValue /* = 0 */){int histSize = 256;const cv::Mat hist = getHist(image, histSize);int mini = 0;for (; mini < histSize; ++mini) {if (hist.at<float>(mini) > minValue)break;}int maxi = histSize - 1;for (; maxi >= 0; --maxi) {if (hist.at<float>(maxi) > minValue)break;}cv::Mat lookup(1, 256, CV_8U);for (int i = 0; i < histSize; ++i) {if (i < mini)lookup.at<uchar>(i) = 0;else if (i > maxi)lookup.at<uchar>(i) = 255;elselookup.at<uchar>(i) = static_cast<uchar>((i - mini) * 255 / (maxi - mini) + 0.5);}cv::Mat result;cv::LUT(image, lookup, result);return result;}
灰度直方图是数字图像中最简单且有用的工具,这一篇主要总结OpenCV中直方图CvHistogram的结构和应用。
灰度直方图的定义
灰度直方图是灰度级的函数,描述图像中该灰度级的像素个数(或该灰度级像素出现的频率):其横坐标是灰度级,纵坐标表示图像中该灰度级出现的个数(频率)。
一维直方图的结构表示为
高维直方图可以理解为图像在每个维度上灰度级分布的直方图。常见的是二维直方图。如红-蓝直方图的两个分量分别表示红光图像的灰度值和蓝光图像灰度值的函数。其图像坐标(Dr,Db)处对应在红光图像中具有灰度级Dr同时在蓝光图像中具有灰度级Db的像素个数。这是基于多光谱——每个像素有多个变量——的数字图像,二维中对应每个像素统计个变量。
OpenCV中的直方图CvHistogram
注意我们在上面理解直方图的意义时更多把他想象成一幅“图”,继而理解图中横坐标,纵坐标的意义。而在OpenCV中,应该更多把直方图看做“数据结构”来理解。
OpenCV中用CvHistogram表示多维直方图(CvHistogram):
typedef struct CvHistogram { int type; CvArr* bins; //存放每个灰度级数目的数组指针 float thresh[CV_MAX_DIM][2]; //均匀直方图 float** thresh2; //非均匀直方图 CvMatND mat; //直方图数组的内部数据结构 } CvHistogram;
这个结构看起来简单(比IplImage*元素少多了。。。)其实并不太好理解。
第一个成员type用来指定第二个成员bins的类型。OpenCv中常见到CvArr*的接口,可以用以指定诸如CvMat、CvMatND、IplImage的类型,其实CvArr*的是一个指向void的指针。在函数内部有时需要得到确切的指向类型,这就需要type来指定。
thresh用来指定统计直方图分布的上下界。比如[0 255]表示用来统计图像中像素分别在灰度级[0 255]区间的分布情况,CV_MAX_DIM对应直方图的维数,假如设定二维红-蓝直方图的thresh为[0 255;100 200],就是分别统计红色图像灰度级在[0 255]以及蓝色图像在灰度级[100 200]的分布情况。
thresh用以指定均匀直方图的分布,我们按每个像素理解自然是“均匀分布”,其实也可以统计像素在几个区间的分布。如果统计像素在2个区间的分布,则对应[0 255]的上下界,均匀分布统计的区间即[0 127] [127 255]分布的概率,这也是为什么thresh第二个维数默认为2——会自动均分上下界;而thresh2指定非均匀的分布,这就需要指定每个区间的上下界,如果要统计直方图在区间(0,10,100,255)的分布,那需要指定thresh2的一个维度为[0 10 100 255],所以用float**形式表示。
mat简单说就是存储了直方图的信息,即我们统计的直方图分布概率。
创建直方图 cvCreateHist()
OpenCV中用cvCreateHist()创建一个直方图:
CvHistogram* cvCreateHist( int dims, //直方图维数 int* sizes,//直翻图维数尺寸 int type, //直方图的表示格式 float** ranges=NULL, //图中方块范围的数组 int uniform=1 //归一化标识 );
size数组的长度为dims,每个数表示分配给对应维数的bin的个数。如dims=3,则size中用[s1,s2,s3]分别指定每维bin的个数。
type有两种:CV_HIST_ARRAY 意味着直方图数据表示为多维密集数组 CvMatND; CV_HIST_TREE 意味着直方图数据表示为多维稀疏数组 CvSparseMat。
ranges就是那个复杂的不好理解的thresh的范围,他的内容取决于uniform的值。uniform为0是均匀的,非0时不均匀。
计算图像直方图的函数为CalcHist():
void cvCalcHist( IplImage** image, //输入图像(也可用CvMat**) CvHistogram* hist, //直方图指针 int accumulate=0, //累计标识。如果设置,则直方图在开始时不被清零。 const CvArr* mask=NULL //操作 mask, 确定输入图像的哪个象素被计数 );
要注意的是这个函数用来计算一张(或多张)单通道图像的直方图,如果要计算多通道,则用这个函数分别计算图像每个单通道。
实践:一维直方图
下面实践一下用OpenCV生成图像的一维直方图
int main( ) { IplImage * src= cvLoadImage("baboon.jpg"); IplImage* gray_plane = cvCreateImage(cvGetSize(src),8,1); cvCvtColor(src,gray_plane,CV_BGR2GRAY); int hist_size = 256; //直方图尺寸 int hist_height = 256; float range[] = {0,255}; //灰度级的范围 float* ranges[]={range}; //创建一维直方图,统计图像在[0 255]像素的均匀分布 CvHistogram* gray_hist = cvCreateHist(1,&hist_size,CV_HIST_ARRAY,ranges,1); //计算灰度图像的一维直方图 cvCalcHist(&gray_plane,gray_hist,0,0); //归一化直方图 cvNormalizeHist(gray_hist,1.0); int scale = 2; //创建一张一维直方图的“图”,横坐标为灰度级,纵坐标为像素个数(*scale) IplImage* hist_image = cvCreateImage(cvSize(hist_size*scale,hist_height),8,3); cvZero(hist_image); //统计直方图中的最大直方块 float max_value = 0; cvGetMinMaxHistValue(gray_hist, 0,&max_value,0,0); //分别将每个直方块的值绘制到图中 for(int i=0;i<hist_size;i++) { float bin_val = cvQueryHistValue_1D(gray_hist,i); //像素i的概率 int intensity = cvRound(bin_val*hist_height/max_value); //要绘制的高度 cvRectangle(hist_image, cvPoint(i*scale,hist_height-1), cvPoint((i+1)*scale - 1, hist_height - intensity), CV_RGB(255,255,255)); } cvNamedWindow( "GraySource", 1 ); cvShowImage("GraySource",gray_plane); cvNamedWindow( "H-S Histogram", 1 ); cvShowImage( "H-S Histogram", hist_image ); cvWaitKey(0); }
试验结果:
对应的,我们可以用一样的思路统计每个通道的直方图,并绘制图像每个通道像素的分布:
实践:二维直方图
我们也可以结合OpenCV的例子生成二维直方图:
IplImage* r_plane = cvCreateImage( cvGetSize(src), 8, 1 ); IplImage* g_plane = cvCreateImage( cvGetSize(src), 8, 1 ); IplImage* b_plane = cvCreateImage( cvGetSize(src), 8, 1 ); IplImage* planes[] = { r_plane, g_plane }; //将HSV图像分离到不同的通道中 cvCvtPixToPlane( src, b_plane, g_plane, r_plane, 0 ); // 生成二维直方图数据结构 int r_bins =256, b_bins = 256; CvHistogram* hist; { int hist_size[] = { r_bins, b_bins }; float r_ranges[] = { 0, 255 }; // hue is [0,180] float b_ranges[] = { 0, 255 }; float* ranges[] = { r_ranges,b_ranges }; hist = cvCreateHist( 2, hist_size, CV_HIST_ARRAY, ranges, 1); } //计算一张或多张单通道图像image(s) 的直方图 cvCalcHist( planes, hist, 0, 0 );
刚才的图我们是对应每个横坐标绘制纵坐标的直方块,二维的图需要绘制每个点:
for( int h = 0; h < r_bins; h++ ) { for( int s = 0; s < b_bins; s++ ) { float bin_val = cvQueryHistValue_2D( hist, h, s ); //查询直方块的值 int intensity = cvRound( bin_val * 255 / max_value ); cvRectangle( hist_img, cvPoint( h*scale, s*scale ), cvPoint( (h+1)*scale - 1, (s+1)*scale - 1), CV_RGB(intensity,intensity,intensity), CV_FILLED); } }
最终生成二维直方图:
直方图的应用以后再讨论。
转载请注明出处:http://blog.csdn.net/xiaowei_cqu/article/details/7600666
实验代码下载:http://download.csdn.net/detail/xiaowei_cqu/4328881
Mat格式的参考这里:http://blog.csdn.net/xiaowei_cqu/article/details/8833799
OpenCV:实现灰度直方图和单通道直方图拉伸相关推荐
- vs+opencv创建空白图片(单通道,三通道)
// 从Mat类的众多重载构造函数Mat中,选择合适的 Mat img = imread("../image/1.png", IMREAD_GRAYSCALE); if (img. ...
- opencv python 灰度图转三通道彩色图
灰度图转RGB IMG_OUT = cv2.cvtColor(IMG_IN, cv2.COLOR_GRAY2RGB)
- 将图像转换为8位单通道_数字图像存储
微信公众号:枫叶AI,专注计算机视觉,机器学习,人工智能等 数字化图像数据有两种存储方式:位图存储(Bitmap)和矢量存储(Vector).位图图像又称作点阵图像.位映射图像,它是由一系列像素组成的 ...
- C#+Opencv图像处理01:16位影像转8位、单通道转三通道、彩色影像增强(直方图均衡化、对数变换和伽马变换)
目录 1 读取不同类型的影像 2 转化16位影像为8位影像 3 将单通道影像转化为多通道影像
- 图像灰度图,直方图,像素通道问题
1.图像直方图概述 直方图广泛运用于很多计算机视觉运用当中,通过标记帧与帧之间显著的边缘和颜色的统计变化,来检测视频中场景的变化.在每个兴趣点设置一个有相近特征的直方图所构成 "标签&q ...
- 【Opencv入门】RGB三通道直方图的计算与绘制
文章目录 一. 直方图概述 Overview of histogram 二.直方图的建立 Establishment of histogram 三.直方图的作用 The function of his ...
- 【 OpenCV】——灰度直方图
[ OpenCV]--灰度直方图 前言 本文介绍了灰度直方图的基础内容. 对函数参数的解释 void calcHist { const Mat*images,//输入的矩阵数组或数据集 int nim ...
- python 灰度直方图_python3+opencv 使用灰度直方图来判断图片的亮暗操作
1.如何让计算机自动判断一张图是否偏暗?或是判断一张图是否是处于夜晚?我们可以先把图片转换为灰度图,然后根据灰度值的分布来判断,如: 我们可以从上图看到,晚上的图片的灰度值是集中在前段的,如0~30多 ...
- python将图像转换为8位单通道_Python OpenCV读取16位单通道图像并转换为8位灰度图显示...
语义.实例分割数据集的标注图像以及一些深度图像等都是由单通道16位整型图像存储的,我们通常需要读取这种图像并显示出来,由于OpenCV一般只能够对8位图像进行显示,也就是像素范围在0-255的图像,而 ...
最新文章
- Openfiler 排错
- codeforces round721 div2. E
- java访问jar中的资源问题代码
- 以太网物理层(PHY)
- mysql8 安装_MySQL8.x安装使用
- conn.setAutoCommit(true) and conn.close() 关系
- GTJ2018如何导出全部工程量_新清单计量规范征求意见稿第二期来啦!来看看那些让你烦恼的操作如何解决...
- Android CircleMenu:旋转转盘选择Menu
- java 合并两个有序数组
- 20幻读是什么,幻读有什么问题
- 卸载掉WPS后安装Office文档图标显示异常
- 软件资格证考试——初级程序员
- 传统行业如何在互联网时代转型
- LZ77压缩算法原理剖析
- 零基础入门推荐系统 - 新闻推荐(一)
- easyuefi只能在基于uefi启动的_只需2个命令,就能将win10 BIOS启动方式转换为UEFI,你见过吗!...
- seo关键字优化条例
- SNMP TRAP发包工具
- css样式实现td显示字体超长显示...鼠标悬浮显示全部
- https+ip实现摄像头访问浏览器
热门文章
- KVM 虚拟化环境安装
- Linux运维工程师:30道面试题整理
- linux svn 常用命令
- ubuntu 下telnet 操纵memcache 实现
- 基础概念总结(spring security、Quartz、JUnit测试)
- 一个民工的数字化生活
- 用vmware-converter4把linux 迁移到ESX4.1中
- 查看、修改linux系统的最大链接数限制、文件描述符限制、端口范围限制、虚拟内存等...
- source ~/.bashrc 什么意思
- jquery.cookie 使用方法