【数字图像处理】直方图均衡化与规定化
目录
直方图处理技术概述
直方图均衡化
公式推导
Matlab代码实现
图像的规定化
数学推导
单映射
Matlab代码
效果展示编辑
组映射
Matlab代码
直方图处理技术概述
灰度级范围[0,L-1]的数字图像，在实际使用中，用哥灰度级的像素总数除以MN的整幅图像的像素数量，归一化到规律中进行计算，这样可以解除面积对哥灰度级在所有橡树中分布的影响，仅仅靠概率进行分析个灰度级的分布。
在较暗的图像中，直方图的分量集中在灰度级的低段。
低对比度图像具有较窄的直方图，且都集中于[0,L-1]的灰度级的某一部分，较为集中。
高对比度图像的直方图的分量覆盖了很宽的灰度范围，且像素数量的分布也没有太不均匀，只有少数灰度级的像素数量要高出许多。
可以直观的得出结论：若以一幅图像的像素占据整个可能的灰度级并且均与分布，则该图像会有高对比度的图像外观并展示灰色调的较大变化，最终效果会得到一幅灰度细节丰富且动态范围较大的图像。
直方图均衡化
直方图均衡化，其目的是为了将图像中的一些动态范围较小、较集中的的灰度级改变，让图像的灰度级分布的更为均与一些。可以依据输入的图像就能得到一个变换函数来自动的实现直方图分布均匀的效果。
我对于直方图的推到论文并没有看过，所有最底层的细节不太了解。全基于我对于均衡化的计算理解。
从实际情况分析
观察这张图的直方图，在整个[0,255]的灰度级中，大部分的灰度级都是集中在较高的灰度部分。整个图片都是偏亮。我们希望提高图像的对比度，让图像的更多细节变的突出，同时希望图像不要太亮，让灰度级像素分布均匀一些。
所以，我们的目的是希望挪动直方图的灰度级，而不是修改直方图的灰度级的数量，大规模修改灰度级的像素的数量会改变整幅图像。且会影响灰度级的分布，有需要重新取计算，所有我们的目的是，不修改各自灰度级的像素数量，而是仅仅移动灰度级。
r是指原始图像的灰度级，s是目的变换后的灰度级。
可以列出某种函数关系
变换前与变换后都是要满足同样的灰度级范围。
在输入一个灰度级的像素，经过变换后就是输出后的一个新的灰度级。
关于变换函数，可以得到其反变换，也就是反函数
而函数存在反函数的情况是，函数是一个严格单调递增函数，这样才不会出现一对多的二义性矛盾。
但是在于实际的图片的变换情况下是，任然存在多个灰度级变换后变成一个灰度级的情况，也就是多对一情况。这是因为，作用于图像的整数，是任然存在多对一的情况。
对于图像的不同的灰度级，对于不同的灰度级的灰度数量，在整幅图像的分布，可以用概率来替代表示，这是属于归一化的思想。
公式推导
得到r的s映射关系后，可以通过相应的映射可以进行灰度级变换。
Matlab代码实现
mat1=imread("Image\gray.jpg");
mat2=mat1;
L=256;%统计一副图像的不同灰度级的在图像中的概率密度（该灰度级的像素数量/整个图像的像素总数量）
[M,N]=size(mat1);
probaDensity=zeros(1,L);%从1开始到256-[0,255]
for a=1:Mfor b=1:NprobaDensity(mat1(a,b)+1)=probaDensity(mat1(a,b)+1) + 1;end
end
for a=1:LprobaDensity(a)=probaDensity(a)/(M*N);
end
%bar(1:L,probaDensity)
%计算累计概率密度-s=T(r)=积分p(r)*(L-1)*dr
cumprobaDensity=zeros(1,L);
cumprobaDensity(1)=probaDensity(1);
for a=2:LcumprobaDensity(a)=cumprobaDensity(a-1)+probaDensity(a);
end%完成变换，乘以（L-1）取整，完成s的映射，可以构建s的映射表
S=zeros(1,L);
for a=1:LS(a)=uint8(cumprobaDensity(a)*(L-1)+0.5);
end
%bar(1:L,S)%完成了r-s的映射表，然后遍历图像，将图像的灰度级(r)对应到变换后的sfor a=1:Mfor b=1:Nmat2(a,b)=S(mat2(a,b)+1);end
endfigure,
subplot(221),imhist(mat1),title("原直方图")
subplot(222),imhist(mat2),title("均衡化直方图")
subplot(223),imshow(mat1),title("原图")
subplot(224),imshow(mat2),title("均衡化");
clf;效果展示
这样一来，灰度级分布的较为均匀，基本上在[0,L-1]的灰度级都覆盖到了。均衡化后的图像对比度提升，细节增加。
图像的规定化
规定化，其本身与均衡化一样，不过均衡化的变换函数/增强函数是由图像自动决定的，让图像的灰度级变的更为均匀。而规定化可以让我们自己指定希望变换的图像的增强函数。让我们可以达到自己想要的效果。可以选择让灰度级集中到低灰度级区域，让阴影的细节更为丰富。
我们提供一个图像，以其的直方图作为变换的模板。
我们提供了两张图，两张图都进行变换，然后得到两个的映射函数表，再进行一次映射，可以得到新的映射表。
数学推导
关键是在于怎么完成Sr到z的反映射。
有两种方法-单映射与组映射，也就是SML与GML算法。
单映射
根据SML的思想，以为基准，让
去查找，差值的绝对值最小的哪个灰度级就是新映射的灰度级。
Matlab代码
mat1=imread("Image\原图.jpg");
tmat1=mat1;
mat0=imread("Image\规定图.jpg");
tmat0=mat0;L=256;[M1,N1]=size(tmat1);
[M0,N0]=size(tmat0);%%统计r与z的概率密度
density0=zeros(1,L);
density1=zeros(1,L);for a=1:M0for b=1:N0density0(tmat0(a,b)+1)=density0(tmat0(a,b)+1)+1;end
end
for a=1:M1for b=1:N1density1(tmat1(a,b)+1)=density1(tmat1(a,b)+1)+1;end
endfor a=1:Ldensity1(a)=density1(a)/(M1*N1);density0(a)=density0(a)/(M0*N0);end%%分别得到关于r的T(r)-累计概率密度*（L-1）
%%关于z的G(z)-累计概率密度*（l-1）cumdensity0=zeros(1,L);
cumdensity1=zeros(1,L);cumdensity0(1)=density0(1);
cumdensity1(1)=density1(1);for a=2:Lcumdensity0(a)=cumdensity0(a-1)+density0(a);cumdensity1(a)=cumdensity1(a-1)+density1(a);end% bar(1:L,cumdensity0)
% bar(1:L,cumdensity1)%%完成了量组对于S的映射，还需要关于Sr的反映射到z
%找Sr与Sk中差值最小的
Rz1=zeros(1,L);%为避免取整误差，直接用累计概率密度进行计算%%单映射以z为基准，进行r中的查找最小的开始映射
temp=0.0;
for a=1:Ltemp=cumdensity1(a);MinDen=1.0;minIndex=1;for b=1:Lif abs(temp-cumdensity0(b)) <= MinDenminIndex=b;MinDen=abs(temp-cumdensity0(b));endend%可以找出最小的差值和下标Rz1(a)=minIndex;
endfor a=1:M1for b=1:N1tmat1(a,b)=Rz1(tmat1(a,b));end
endfigure,
subplot(231),imhist(mat1),title("原直方图")
subplot(232),imhist(mat0),title("规定直方图")
subplot(233),imhist(tmat1),title("单映射直方图")
subplot(234),imshow(mat1),title("原图");
subplot(235),imshow(mat0),title("规定图");
subplot(236),imshow(tmat1),title("单映射后图");
clf;%%组映射以r为基准，查找r中的累计概率密度的最小差值效果展示
组映射
都说单映射的查找方法是存在误差的，所以提供了组映射。
组映射与单映射反过来，以原图图的为基准，让规定图的
去查找原图中累计概率密度的差值的绝对值最小的哪个灰度级，依次为一个分别，让从0到这个灰度级之间的所有灰度级都映射到当前的规定图的灰度级，然后更新边界，直到下一次找到最小的边界。
具体组映射和单映射的理解可以参考直方图规定化：单映射规则、组映射规则的手动求解过程_Nefu_lyh的博客-CSDN博客_直方图规定化例题
这个博主讲解的还是很清晰的。
Matlab代码
mat1=imread("Image\原图.jpg");
tmat1=mat1;
mat0=imread("Image\规定图.jpg");
tmat0=mat0;L=256;[M1,N1]=size(tmat1);
[M0,N0]=size(tmat0);%%统计r与z的概率密度
density0=zeros(1,L);
density1=zeros(1,L);for a=1:M0for b=1:N0density0(tmat0(a,b)+1)=density0(tmat0(a,b)+1)+1;end
end
for a=1:M1for b=1:N1density1(tmat1(a,b)+1)=density1(tmat1(a,b)+1)+1;end
endfor a=1:Ldensity1(a)=density1(a)/(M1*N1);density0(a)=density0(a)/(M0*N0);end%%分别得到关于r的T(r)-累计概率密度*（L-1）
%%关于z的G(z)-累计概率密度*（l-1）cumdensity0=zeros(1,L);
cumdensity1=zeros(1,L);cumdensity0(1)=density0(1);
cumdensity1(1)=density1(1);for a=2:Lcumdensity0(a)=cumdensity0(a-1)+density0(a);cumdensity1(a)=cumdensity1(a-1)+density1(a);end
%%完成了量组对于S的映射，还需要关于Sr的反映射到z
%找Sr与Sk中差值最小的
Rz2=zeros(1,L);%为避免取整误差，直接用累计概率密度进行计算
%%组映射以r为基准，查找r中的累计概率密度的最小差值，以组为单位进行映射
%设置一个边界，在查找完一轮后，映射区域内的r，然后更新边界
bound=1;
temp=0.0;
for a=1:Ltemp=cumdensity0(a);MinDen=1.0;minIndex=1;for b=1:Lif abs(temp-cumdensity1(b))<=MinDenminIndex=b;MinDen=abs(temp-cumdensity1(b));endend%查找完，开始映射[bound,minIndex]for cd=bound:minIndexRz2(cd)=a;end%更新边界bound=minIndex+1;
endfor a=1:M1for b=1:N1tmat1(a,b)=Rz2(tmat1(a,b));end
endfigure,
subplot(231),imhist(mat1),title("原直方图")
subplot(232),imhist(mat0),title("规定直方图")
subplot(233),imhist(tmat1),title("组映射直方图")
subplot(234),imshow(mat1),title("原图");
subplot(235),imshow(mat0),title("规定图");
subplot(236),imshow(tmat1),title("组映射后图");
clf;效果展示
【数字图像处理】直方图均衡化与规定化相关推荐
- 图像直方图均衡化c语言程序,数字图像处理——直方图均衡化
1.相关了解 直方图均衡化又称为直方图平坦化.直方图均衡化的基本思想是将原始图像的不均衡的直方图变换为均匀分布的形式.即将输入图像转化为在每一灰度级上都有相同的像素点数(即输出的直方图是平坦的,其分布 ...
- 数字图像处理 直方图均衡化 MATLAB实验
一.原理 直方图均衡化是通过对原图像进行某种变换使原图像的灰度直方图修正为均匀的直方图的一种方法. 二.步骤 ①读入原图像huafen.jpg,并显示图像及其直方图: ②对图像进行直方图均衡化处理,并 ...
- 数字图像处理——直方图均衡化
直方图均衡化 直方图均衡化(Histogram equalization)是一种常用的灰度变换方法. 基本原理 直方图均衡化的基本原理是:对在图像中像素个数多的灰度值(即对画面起主要作用的灰度值)进行 ...
- 数字图像处理实验之对比度拉伸、直方图均衡化和规定化
[实验内容]: 一.运用Matlab对Lena进行对比度拉伸增强(利用 imadjust ) 二.运用Matlab对图像进行直方图均衡化和规定化 Matlab版本:R2015a 一.运用Matl ...
- 数字图像处理——直方图累计
数字图像处理--直方图累计 %函数名称为Image_Hist,输入参数Image,输出参数Hist function [Hist] = Image_Hist(Image) %获取矩阵的行.列.波段数 ...
- [图像处理] 直方图均衡化原理 - 数学推导
直方图均衡化 效果 代码 import cv2 as cv import numpy as np import matplotlib.pyplot as pltsrc = cv.imread(&quo ...
- 图像处理-直方图均衡化
个人博客:http://www.chenjianqu.com/ 原文链接:http://www.chenjianqu.com/show-11.html 概念 在统计学中,直方图(英语:Histogra ...
- 直方图均衡化取整怎么计算_玩转直方图处理之直方图均衡化、规定化
冒泡....双十一刚过~购物狂欢完还是要收心学习鸭!今天来说一说直方图. 直方图 定义: 直方图是一种统计报告图,由一系列高度不等的纵向条纹或线段表示数据分布的情况. 一般用横轴表示数据类型,纵轴表示 ...
- 详细解析图像处理 直方图均衡化计算
什么是直方图均衡化 看一个图了解一下什么是直方图均衡化: 第一个图灰度都集中在左边,整体图像较暗 第二个图灰度都集中在右边,整体图像较亮 第三个图灰度都集中在中间,整体图像适中,但是雾蒙蒙的并不清晰 ...
最新文章
- oracle如何复制表的索引,Oracle表与索引管理
- Java 序列化 之 Serializable
- js原生实现点击事件只创建一个元素
- GitHub 创建项目
- [CityLife]“背后的故事”---贫嘴曾志伟
- ElasticSearch前缀匹配查询和范围查询（中文检索）
- Cordova 快速入门记录
- 【动态规划】路径dp：蓝桥2020： 数字三角形
- C++/Python/Java/C，四大语言对比，你更喜欢哪种？
- 05、ADS使用记录之集总参数匹配
- 钉钉自动打卡上下班微信自动打卡（原理及源码）
- Redis五大数据类型及常见命令
- fastadmin-addons使用
- 失业在家做什么赚钱好？失业在家怎么赚钱？
- [入门篇]初识操作系统OS(这篇博客将帮你打开操作系统OS大门的第一步 超超超超超详细)
- java计算机毕业设计web企业档案管理系统MyBatis+系统+LW文档+源码+调试部署
- 升级iOS8和iOS9系统后，保险箱Pro、私人保险箱、私密相冊打开就闪退的官方解决方式...
- Docker Swarm nginx 集群搭建
- 按键精灵生成小精灵获取服务器信息失败,{为什么我在按键精灵调试是正常的，可是一生成小精灵就提示错误了呢} _ 会员专区 - 按键精灵论坛...
- word按空格键时显示很多小点，怎么办？