图象增强对比度.pptVIP

第三章 图象增强 对比度增强 图象平滑 图象锐化 同态滤波 伪彩色与假彩色处理 （一）直方图均衡化 Histogram Equalization 直方图:表示数字图象中的每一灰度级与其出现的频率(该灰度级的象素数目)间的统计关系,用横坐标表示灰度级, 纵坐标表示频数(也可用概率表示)。 直方图调整法 直方图 I=imread(rice.tif); [M,N]=size(I); s = 256; %灰度级 H = zeros(1,s); for i=1:M for j=1:N g=I(i,j); H(1,g+1) = H(1,g+1) + 1; end; end; 计算直方图 灰度直方图 彩色直方图 直方图均衡化是将原图象的直方图通过变换函数修正为均匀的直方图，然后按均衡直方图修正原图象。 图象均衡化处理后，图象的直方图是平直的，即各灰度级具有相同的出现频数，那么由于灰度级具有均匀的概率分布，图象看起来就更清晰了。 直方图均衡化 直方图均衡化 首先假定连续灰度级的情况，推导直方图均衡化变换公式，令r代表灰度级，P ( r ) 为概率密度函数。 r值已归一化，最大灰度值为1。 灰度级数为0,1,…,L-1共L级。 ? ? 连续灰度的直方图非均匀分布 连续灰度的直方图均匀分布 直方图均衡化目标 直方图均衡化 直方图均衡化 要找到一种变换 S=T ( r ) 使直方图变平直，为使变换后的灰度仍保持从黑到白的单一变化顺序，且变换范围与原先一致，以避免整体变亮或变暗。必须规定： （1）在0≤r≤1中，T(r)是单调递增函数，且0≤T(r)≤1； （2）反变换r=T-1(s),T-1(s)也为单调递增函数，0≤s≤1。 ? rj rj+?r sj sj+?s 直方图均衡化变换公式推导图示 直方图均衡化 考虑到灰度变换不影响象素的位置分布，也不会增减象素数目。所以有 ? 直方图均衡化 应用到离散灰度级，设一幅图象的象素总数为n，分L个灰度级。 nk: 第k个灰度级出现的频数。 第k个灰度级出现的概率 P(rk)=nk/n 其中0≤rk≤1，k=0,1,2,...,L-1 形式为： ? 例：设图象有64*64=4096个象素，有8个灰度级，灰度分布如表所示。进行直方图均衡化。 rk r0=0 r1=1/7 r2=2/7 r3=3/7 r4=4/7 r5=5/7 r6=6/7 r7=1 nk 790 1023 850 656 329 245 122 81 p(rk) 0.19 0.25 0.21 0.16 0.08 0.06 0.03 0.02 1. 由（3-2）式计算sk rk r0=0 r1=1/7 r2=2/7 r3=3/7 r4=4/7 r5=5/7 r6=6/7 r7=1 nk 790 1023 850 656 329 245 122 81 p(rk) 0.19 0.25 0.21 0.16 0.08 0.06 0.03 0.02 sk计算 0.19 0.44 0.65 0.81 0.89 0.95 0.98 1.00 rk r0=0 r1=1/7 r2=2/7 r3=3/7 r4=4/7 r5=5/7 r6=6/7 r7=1 nk 790 1023 850 656 329 245 122 81 p(rk) 0.19 0.25 0.21 0.16 0.08 0.06 0.03 0.02 sk计算 0.19 0.44 0.65 0.81 0.89 0.95 0.98 1.00 sk舍入 1/7 3/7 5/7 6/7 6/7 1 1 1 2. 把计算的sk就近安排到8个灰度级中 rk r0=0 r1=1/7 r2=2/7 r3=3/7 r4=4/7 r5=5/7 r6=6/7 r7=1 nk 790 1023 850 656 329 245 122 81 p(rk) 0.19 0.25 0.21 0.16 0.08 0.06 0.03 0.02 sk计算 0.19 0.44 0.65 0.81 0.89 0.95 0.98 1.00 sk舍入 1/7 3/7 5/7 6/7 6/7 1 1 1 sk s0 s1 s2 s3 s4 nsk 790 1023 850 985 448 p(sk) 0.19 0.