数字图像处理笔记（五）：直方图均衡与直方图匹配（规定化） |
您所在的位置：网站首页 › 直方图匹配的基本原理 › 数字图像处理笔记（五）：直方图均衡与直方图匹配（规定化） |
数字图像处理笔记（五）：直方图均衡与直方图匹配（规定化）
|
1 - 直方图处理
直方图是多种空间域处理技术的基础，可以用于图像增强。 2 - 直方图均衡考虑连续灰度值，并用变量r表示待处理图像的灰度。设r的取值区间为[0,L-1],且r = 0表示黑色，r = L - 1表示白色。令 p r ( r ) p_r(r) pr(r)和 p s ( s ) p_s(s) ps(s)分别表示随机变量r和s的概率密度函数，概率密度函数归一化后可以画出图像的直方图 import numpy as np import cv2 as cv import matplotlib.pyplot as plt image = cv.imread('images/test_1.jpg', 0) hist = cv.calcHist([image],[0],None,[256],[0,255]) plt.plot(hist) plt.show()
可以看到较暗的图片灰度直方图分布的非常不均衡，那么我们构造一个出变换 s = T ( r ) = ( L − 1 ) ∫ 0 r p r ( w ) d w , 0 ≤ r ≤ L − 1 s=T(r)=(L-1)\int_{0}^{r}p_r(w)dw,0\leq r \leq L-1 s=T(r)=(L−1)∫0rpr(w)dw,0≤r≤L−1 使得对应的输出图像概率密度函数变成均匀的 对于离散值，我们处理其概率（直方图值）与求和来替代处理概率密度函数与积分。 一副数字图像中灰度级 r k r_k rk出现的概率近似为： p r ( r k ) = n k M N , k = 0 , 1 , 2 , … , L − 1 p_r(r_k)=\frac{n_k}{MN}, k=0,1,2,\dots,L-1 pr(rk)=MNnk,k=0,1,2,…,L−1 其中，MN是图像中像素的综述， n k n_k nk是灰度为 r k r_k rk像素的个数。 则离散形式的变化公式为： s k = T ( r k ) = ( L − 1 ) ∑ j = 0 k p r ( r j ) = ( L − 1 ) M N ∑ j = 0 k n j , k = 0 , 1 , 2 , … , L − 1 s_k=T(r_k)=(L-1)\sum_{j=0}^kp_r(r_j)=\frac{(L-1)}{MN}\sum_{j=0}^kn_j, k=0,1,2,\dots,L-1 sk=T(rk)=(L−1)j=0∑kpr(rj)=MN(L−1)j=0∑knj,k=0,1,2,…,L−1 import numpy as np import cv2 as cv import matplotlib.pyplot as plt image = cv.imread('images/test_1.jpg', 0) image2 = cv.equalizeHist(image) cv.imshow("image_2",image2) hist = cv.calcHist([image2],[0],None,[256],[0,255]) plt.plot(hist) plt.show() cv.waitKey(0)
可以看到直方图的分布更加均匀了，图像的对比度和相关细节显示的也更加明显了。 3 - 直方图匹配直方图均衡能自动地确定变换函数，当需要自动增强时，这是一种好方法，因为由于这种技术得到的结果可以预知，并且这种方法实现也很简单，但是对于某些应用，采用均匀直方图的基本增强并不是最好的方法，有时我们希望处理后的图像具有规定的直方图形状可能更有用。这种采用处理后有特殊直方图的方法称为直方图匹配或直方图规定划。 例：直方图规定划 设一副图片的灰度概率密度函数（PDF）为 p r ( r ) = 2 r / ( L − 1 ) 2 , 0 ≤ r ≤ ( L − 1 ) p_r(r)=2r/(L-1)^2, 0\leq r\leq (L-1) pr(r)=2r/(L−1)2,0≤r≤(L−1) 我们想要寻找一个变换函数，使得产生的图像的灰度PDF是 p z ( z ) = 3 z 2 / ( L − 1 ) 3 , 0 ≤ r ≤ ( L − 1 ) p_z(z)=3z^2/(L-1)^3,0\leq r\leq (L-1) pz(z)=3z2/(L−1)3,0≤r≤(L−1) 求出原图像的直方图均衡变换： s = T ( r ) = ( L − 1 ) ∫ 0 r p r ( w ) d w = 2 ( L − 1 ) ∫ o r w d w = r 2 ( L − 1 ) , 0 ≤ r ≤ L − 1 s=T(r)=(L-1)\int_{0}^{r}p_r(w)dw=\frac{2}{(L-1)}\int_{o}^rwdw=\frac{r^2}{(L-1)},0\leq r \leq L-1 s=T(r)=(L−1)∫0rpr(w)dw=(L−1)2∫orwdw=(L−1)r2,0≤r≤L−1 求出规定图像的直方图均衡变换： G ( z ) = ( L − 1 ) ∫ 0 z p z ( w ) d w = 3 ( L − 1 ) 2 ∫ o z w 2 d w = z 3 ( L − 1 ) 2 , 0 ≤ r ≤ L − 1 G(z)=(L-1)\int_{0}^{z}p_z(w)dw=\frac{3}{(L-1)^2}\int_{o}^zw^2dw=\frac{z^3}{(L-1)^2},0\leq r \leq L-1 G(z)=(L−1)∫0zpz(w)dw=(L−1)23∫ozw2dw=(L−1)2z3,0≤r≤L−1 最后使其相等 G ( z ) = s G(z)=s G(z)=s 利用反变换得到 z = G − 1 [ T ( r ) ] = G − 1 ( s ) z=G^-1[T(r)]=G^-1(s) z=G−1[T(r)]=G−1(s) 画出函数的关系图就很好理解了 在离散形式下 s k = T ( r k ) = ( L − 1 ) ∑ j = 0 k p r ( r j ) = ( L − 1 ) M N ∑ j = 0 k n j , k = 0 , 1 , 2 , … , L − 1 s_k=T(r_k)=(L-1)\sum_{j=0}^kp_r(r_j)=\frac{(L-1)}{MN}\sum_{j=0}^kn_j,k=0,1,2,\dots,L-1 sk=T(rk)=(L−1)j=0∑kpr(rj)=MN(L−1)j=0∑knj,k=0,1,2,…,L−1 G ( z q ) = ( L − 1 ) ∑ i = 0 q p z ( z i ) G(z_q)=(L-1)\sum_{i=0}^qp_z(z_i) G(zq)=(L−1)i=0∑qpz(zi) 对于一个q值，有 G ( z q ) = s k G(z_q)=s_k G(zq)=sk 再利用反变换，有 Z q = G − 1 ( s k ) Z_q=G^-1(s_k) Zq=G−1(sk) 那么我们总结直方图规定划过程如下： 计算给定图像的直方图 p r ( r ) p_r(r) pr(r)，并用它寻找直方图均衡变换，把 s k s_k sk四舍五入为范围[0,L-1]内的整数。计算变换函数G的所有值，将G的值储存在一个表中对于每个值，使用储存的G值寻找相应的 z Q z_Q zQ值，使 G （ z q ） G（z_q） G（zq）最接近 s k s_k sk。首先对输入图像进行均衡，然后使用步骤3找到的映射把该图像中的每个均衡后的像素值 s k s_k sk映射为直方图规定划后的图像相应 z q z_q zq值具体python实现： 将一张较为黑暗的原始图像 实验结果： |
与一张拍摄光线较好的匹配图像
进行直方图匹配，目的是为了把原始图像的直方图匹配成匹配图像的直方图分布，使得原始图像能够呈现出匹配图像的光线效果。
【本文地址】