直方图均衡化推导整理

这篇博客详细介绍了图像均衡化的数学推导过程,通过概率论的知识,从原始像素的概率分布出发,目标是得到理想情况下均匀分布的概率。作者简化了连续图像的假设,并利用概率密度函数和反函数的关系,推导出像素值的转换函数T(x)。最终得出均衡化公式,展示了如何从原始像素分布p(x)计算得到新的分布p(y)。博客适合对图像处理和概率论有一定基础的读者。
摘要由CSDN通过智能技术生成

看了几篇别人的博客,不是图像有问题,就是讲的不清晰(自己知识面不够),所以整理一下。
话不多说,直接进行推导,图像什么的我就不放了。

设随机变量 x , y x,y x,y为均衡化前后的值,由 y = T ( x ) y = T(x) y=T(x)
p ( x ) , p ( y ) p(x),p(y) p(x),p(y) x , y x,y x,y的概率密度函数。其中, p ( x ) = n k N p(x) = \frac{n_k}{N} p(x)=Nnk n k 代 表 某 一 像 素 值 的 个 数 , N 代 表 图 片 的 像 素 个 数 n_k代表某一像素值的个数,N代表图片的像素个数 nkN,显然 p ( x ) p(x) p(x)是已知的。
由统计观察可以知道,要使图像的对比度强烈,感官更好,其直方图应该是更加均匀,理想条件下应该是一条直线。那么其各级概率分布就是 1 L − 1 , L = 256 \frac{1}{L-1},L = 256 L11,L=256,所以 p ( y ) = 1 L − 1 p(y) = \frac{1}{L-1} p(y)=L11
那么这个问题就变成了已知 p ( x ) , p ( y ) , 求 T ( x ) p(x),p(y),求T(x) p(x),p(y)T(x)

由于图像的分布函数是离散的,我们暂不做考虑,假设其是连续的。
下面复习一下概率论的知识:

已知 p ( x ) , x 和 y 的 函 数 为 T ( x ) , 求 p ( y ) p(x),x和y的函数为T(x),求p(y) p(x),xyT(x),p(y)
设X有概率密度函数 f ( x ) f(x) f(x) Y = g ( x ) Y = g(x) Y=g(x) g g g单调增加,且可到,则其反函数 X = h ( y ) X = h(y) X=h(y)存在。
任取 y y y
P ( Y ≤ y ) = P ( g ( x ) ≤ y ) = P ( x ≤ h ( y ) ) = ∫ − ∞ h ( y ) f ( t ) d t P(Y≤y) = P(g(x)≤y) = P(x≤h(y)) = \int_{-∞}^{h(y)}f(t)dt P(Yy)=P(g(x)y)=P(xh(y))=h(y)f(t)dt
(这一串要好好理解,不懂的翻翻书,做两道题就懂了)

然后做一下变换,该替换的替换掉。
F ( y ) = P ( g ( x ) ≤ y ) = F ( x ) = ∫ − ∞ T − 1 ( y ) f ( t ) d t F(y) = P(g(x)≤y) =F(x)= \int_{-∞}^{T^{-1}(y)}f(t)dt F(y)=P(g(x)y)=F(x)=T1(y)f(t)dt

对上式进行求导得:
f ( y ) = f ( T − 1 ( y ) ∗ d [ T − 1 ( y ) ] d y ) = f ( x ) ∗ d x d y f(y) = f(T^{-1}(y)*\frac{d[T^{-1}(y)]}{dy})=f(x)*\frac{dx}{dy} f(y)=f(T1(y)dyd[T1(y)])=f(x)dydx
(这里的变换也要好好理解下,就是把符号变了变)

再继续变换符号
p ( y ) = p ( x ) ∗ d x d y p(y) = p(x)*\frac{dx}{dy} p(y)=p(x)dydx

带入已知的值
1 L − 1 = n k N ∗ d x d [ T ( x ) ] \frac{1}{L-1} = \frac{n_k}{N}*\frac{dx}{d[T(x)]} L11=Nnkd[T(x)]dx

整理后可得:
d [ T ( x ) ] = n k ( L − 1 ) N ∗ d x d[T(x)] = \frac{n_k(L-1)}{N}*dx d[T(x)]=Nnk(L1)dx
对两边进行积分,因为图像函数并不是连续的,所以积分只是各级求和

T ( x ) = L − 1 N ∑ i = 0 k n i T(x) = \frac{L-1}{N}\sum_{i=0}^kn_i T(x)=NL1i=0kni
到此,推导结束。主要还是几个公式符号的变化可能会绕一下。

借鉴:
https://blog.csdn.net/schwein_van/article/details/84336633
https://blog.csdn.net/spongebob1234/article/details/77778709

  • 1
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值