直方图均衡化原理
2015-07-27 15:54
106 查看
项目中遇到要直方图均衡化,在网上找到如下介绍:
直方图均衡化算法分为三个步骤,第一步是统计直方图每个灰度级出现的次数,第二步是累计归一化的直方图,第三步是计算新的像素值。
第一步:
for(i=0;i<height;i++){
for(j=0;j<width;j++){
n[s[i][j]]++;
}
}
for(i=0;i<L;i++){
p[i]=n[i]/(width*height);
}
这里,n[i]表示的是灰度级为i的像素的个数,L表示的是最大灰度级,width和height分别表示的是原始图像的宽度和高度,所以,p[i]表示的就是灰度级为i的像素在整幅图像中出现的概率(其实就是p[]这个数组存储的就是这幅图像的归一化之后的直方图)。
第二步:
for(i=0;i<=L;i++){
for(j=0;j<=i;j++){
c[i]+=p[j];
}
}
c[]这个数组存储的就是累计的归一化直方图。
第三步:
max=min=s[0][0];
for(i=0;i<height;i++){
for(j=0;j<width;j++){
if(max<s[i][j]){
max=s[i][j];
}else if(min>s[i][j]){
min=s[i][j];
}
}
}
找出像素的最大值和最小值。
for(i=0;i<height;i++){
for(j=0;j<width;j++){
t[i][j]=c[s[i][j]]*(max-min)+min;
}
}
t[][]就是最终直方图均衡化之后的结果。
处理前的图片:
![](http://p.blog.csdn.net/images/p_blog_csdn_net/vincentzhao2009/EntryImages/20090918/source.jpg)
处理后的图片:
![](http://p.blog.csdn.net/images/p_blog_csdn_net/vincentzhao2009/EntryImages/20090918/destination.jpg)
对于彩色的图片来说,直方图均衡化一般不能直接对R、G、B三个分量分别进行上述的操作,而要将RGB转换成HSV来对V分量进行直方图均衡化的操作。
直方图均衡化算法分为三个步骤,第一步是统计直方图每个灰度级出现的次数,第二步是累计归一化的直方图,第三步是计算新的像素值。
第一步:
for(i=0;i<height;i++){
for(j=0;j<width;j++){
n[s[i][j]]++;
}
}
for(i=0;i<L;i++){
p[i]=n[i]/(width*height);
}
这里,n[i]表示的是灰度级为i的像素的个数,L表示的是最大灰度级,width和height分别表示的是原始图像的宽度和高度,所以,p[i]表示的就是灰度级为i的像素在整幅图像中出现的概率(其实就是p[]这个数组存储的就是这幅图像的归一化之后的直方图)。
第二步:
for(i=0;i<=L;i++){
for(j=0;j<=i;j++){
c[i]+=p[j];
}
}
c[]这个数组存储的就是累计的归一化直方图。
第三步:
max=min=s[0][0];
for(i=0;i<height;i++){
for(j=0;j<width;j++){
if(max<s[i][j]){
max=s[i][j];
}else if(min>s[i][j]){
min=s[i][j];
}
}
}
找出像素的最大值和最小值。
for(i=0;i<height;i++){
for(j=0;j<width;j++){
t[i][j]=c[s[i][j]]*(max-min)+min;
}
}
t[][]就是最终直方图均衡化之后的结果。
处理前的图片:
![](http://p.blog.csdn.net/images/p_blog_csdn_net/vincentzhao2009/EntryImages/20090918/source.jpg)
处理后的图片:
![](http://p.blog.csdn.net/images/p_blog_csdn_net/vincentzhao2009/EntryImages/20090918/destination.jpg)
对于彩色的图片来说,直方图均衡化一般不能直接对R、G、B三个分量分别进行上述的操作,而要将RGB转换成HSV来对V分量进行直方图均衡化的操作。
相关文章推荐
- 【HDU】5298 Solid Geometry Homework【暴力+二分图】
- RHEL 6.5安装DRBD
- Linux_LVM Couldn't find device with uuid
- JSVC启动Tomcat
- javaWeb-mvc之利用c3p0写入数据库出现乱码
- python基础教程总结15——3 XML构建网址
- AndroidManifest.xml -【manifest】
- Android New Relic 4.244 Agent 安装问题
- 四层负载均衡和七层负载均衡的区别
- [Tue, 21 Jul 2015 ~ Mon, 27 Jul 2015] Deep Learning in arxiv
- 浅谈Linux进程
- mysql 分组后取每个组内最新的一条数据
- 在服务器上运行Jar包
- 百度地图API(二)
- 浅谈Linux进程
- hdu 1398 Square Coins (母函数)
- Oracle中的延迟块的清除方法
- orachk现在软件里面都已经内置了
- iCloudManager部署---iCloudManager的安装及设置(三)
- html5之indexdb(nosql存储)