1、130452 计算题复习一、 直方图均衡化（P68）对已知图像进行直方图均衡化修正。例：表 1 为已知一幅总像素为 n=6464 的 8bit 数字图像（即灰度 级数为 8），各灰度 级（出现的频率）分布列于表中。要求将此幅图像进行均衡化修正（ 变换），并画出修正（变换）前后的直方图。表 1原图像灰度级 rk原各灰度级像素个数nk原分布概率 pr(rk)r0=0 790 0.19r1=1 1023 0.25r2=2 850 0.21r3=3 656 0.16r4=4 329 0.08r5=5 245 0.06r6=6 122 0.03r7=7 81 0.02解：对已知图像均衡化过程见下表：原图
2、像灰度级 rk原各灰度级像素个数nk原分布概率 pr(rk)累积分布函数 sk 计取整扩展sk 并确定映射对应关系rks k新图像灰度级 sk新图像各灰度级像素个数 nsk新图像分布概率ps(sk)r0=0 790 0.19 0.19 1 01 1 790 0.19r1=1 1023 0.25 0.44 3 13 3 1023 0.25r2=2 850 0.21 0.65 5 25 5 850 0.21r3=3 656 0.16 0.81 6 36r4=4 329 0.08 0.89 6 46 6 985 0.24r5=5 245 0.06 0.95 7 57r6=6 122 0.03 0.9
3、8 7 67r7=7 81 0.02 1.00 7 777 448 0.11画出直方图如下：（a）原始图 像直方图 （b）均衡化后直方图*以下部分不用写在答题中。其中： rk、nk 中 k = 0，1，72 pr（rk）= nkn，即计算各灰度级像素个数占所有像素个数的百分比，其中 ，在此 题中 n=6464。kjn0 ，即计算在本灰度级之前（包含本灰度级）所有百分比之和。kjjrs0)(计 ，其中 L 为图像的灰度级数（本 题中 L = 8），int 表示对方括号中的数字取整。5.01int计并 kksL 并ks nsk 为映射对应关系 rks k 中 rk 所对应的 nk 之和。 ，或为映
4、射对应关系 rks k 中 rk 所对应的 pr(rk)之和。npsks/)(3二、 模板运算 使用空间低通滤波法对图像进行平滑操作（P80）空间低通滤波法是应用模板卷积方法对图像每一个像素进行局部处理。模板（或称掩模）就是一个滤波器，它的响应为 H(r,s)，于是滤波输出的数字图像 g(x,y)用离散卷积表示为 )6.24(,(,llskr srHrfyxg式中：x，y = 0 ，1，2，N -1；k 、l 根据所选邻域大小来决定。具体过程如下：（1）将模板在图像中按从左到右、从上到下的顺序移动，将模板中心与每个像素依次重合（边缘像素除外） ；（2）将模板中的各个系数与其对应的像素一一相乘，
5、并将所有的结果相加；（3）将（2）中的结果赋给图像中对应模板中心位置的像素。对于空间低通滤波器而言，采用的是低通滤波器。由于模板尺寸小，因此具有计算量小、使用灵活、适于并行计算等优点。常用的 3*3 低通滤波器（模板）有：191H1202 12463H1084 01245H模板不同，邻域内各像素重要程度也就不同。但无论怎样的模板，必须保证全部权系数之和为 1，这样可保证输出图像灰度值在许可范围内，不会产生灰度“溢出”现象。例：以 为模板，对下图做低通滤波处理，写出处理结果。1202H1 7 1 8 1 7 1 11 1 1 5 1 1 1 11 1 5 5 5 1 1 71 1 5 5 5 1
6、 8 18 1 1 5 1 1 1 18 1 1 5 1 1 8 11 1 1 5 1 1 1 11 7 1 8 1 7 1 1解：低通滤波的步骤为：（1）将模板在图像中按从左到右、从上到下的顺序移动，将模板中心与每个像素依次重合（边缘像素除外）；（2）将模板中的各个系数与其对应的像素一一相乘，并将所有的 结果相加；（3）将（2）中的结果赋给图像中对应模板中心位置的像素。如图中第 2 行第 2 列处的值 = （1*1+1*7+1*1+1*1+2*1+1*1+1*1+1*1+1*5）/10 = 2（其他位置同样方法计算可得）由此步骤可得处理结果为（空白处自己计算后填入）1 7 1 8 1 7 1
7、 11 2 11 71 18 18 11 11 7 1 8 1 7 1 14三、 中值滤波与邻域平均中值滤波（P81）中值滤波是对一个滑动窗口内的诸像素灰度值排序，用中间值代替窗口中心像素的原来灰度值，它是一种非线性的图像平滑法。它对脉冲干扰及椒盐噪声的抑制效果好，在抑制随机噪声的同时能有效保护边缘少受模糊。但它对点、线等细节较多的图像却不太合适。局部平滑法（邻域平均法 或 移动平均法） （P76）局部平滑法是一种直接在空间域上进行平滑处理的技术。用邻域内各像素的灰度平均值代替该像素原来的灰度值，实现图像的平滑。设有一幅 NN 的图像 f(x,y)，若平滑图像为 g(x,y)，则有 )1.24
8、(,1,(,sjifM式中 x,y = 0,1,N-1；s 为（x,y）邻域内像素坐标的集合；M 表示集合 s 内像素的总数。可见邻域平均法就是将当前像素邻域内各像素的灰度平均值作为其输出值的去噪方法。设图像中的噪声是随机不相关的加性噪声，窗口内各点噪声是独立同分布的，经过（4.2.1）平滑后，信号与噪声的方差比可望提高 M 倍。这种算法简单，但它的主要缺点是在降低噪声的同时使图像产生模糊，特别在边缘和细节处。而且邻域越大，在去噪能力增强的同时模糊程度越严重。例：对下图做 3*3 中值滤波处理和 3*3 邻域平均处理，写出处理结果，并比 较邻域平均与中值滤波的差异。1 7 1 8 1 7 1
9、11 1 1 5 1 1 1 11 1 5 5 5 1 1 71 1 5 5 5 1 8 18 1 1 5 1 1 1 18 1 1 5 1 1 8 11 1 1 5 1 1 1 11 7 1 8 1 7 1 1解：（1）中值滤波是对一个滑动窗口内的诸像素灰度值排序，用中间值代替窗口中心像素的原来灰度值，它是一种非线性的图像平滑法。题目中的图像经 3*3 中值滤波后的结果为（忽略边界）：1 7 1 8 1 7 1 11 1 5 5 5 1 1 11 1 5 5 5 1 1 71 1 5 5 5 1 1 18 1 1 5 1 1 1 18 1 1 1 1 1 1 11 1 1 1 1 1 1 1
10、1 7 1 8 1 7 1 1（2）局部平滑法（邻域平均法 或 移动平均法）是用邻域内各像素的灰度平均值代替该像素原来的灰度值， 实5现图像的平滑。题目中的图像经 3*3 局部平滑法（邻域平均法 或 移动平均法）后的结果为（忽略边界）：1 7 1 8 1 7 1 11 19/9 38/9 40/9 38/9 23/9 21/9 11 71 18 18 11 11 7 1 8 1 7 1 1中值滤波法和局部平滑法（邻域平均法 或 移动平均法）均能有效削弱椒盐噪声，但中 值滤波法比邻域平均法更有效，且滤波后图像中的 轮廓比较清晰。6四、 霍夫曼编码（P124）例：设有一信源 A=a1, a2, a
11、3, a4, a5, a6，对应概率 P=0.1, 0.4, 0.06, 0.1, 0.04, 0.3.（1）进行霍夫曼编码（要求大概率的赋码字 0，小概率的 赋码 字 1），给出码字；（2）计算平均码长，信源熵和编码 效率。解：（1）编码步骤1）缩减信源符号数量将信源符号按出现概率从大到小排列，然后 结合2）对每个信源符号赋值从（消减到）最小的信源开始，逐步回到初始信源由此可得哈夫曼编码结果见下表符号 a1 a2 a3 a4 a5 a6概率 0.1 0.4 0.06 0.1 0.04 0.06编码结果 011 1 01010 010001011 00平均码长 2.504.6.04.3.02.
12、4.010 LiipB（其中， 是灰度值为 i 的编码长度， 为灰度值为 的概率，L 为灰度级数）i ipia信源熵 14.2log10iLiipH编码效率 973.B0.40.30.10.10.60.40.40.30.10.10.10.40.30.10.2 0.40.30.3 0.40.61 2 3 4a12a5a3a6a4a信信 信信信信 信0.40.30.10.10.60.40.40.30.10.10.10.40.30.10.2 0.40.30.3 0.40.61 2 3 4a12a5a3a6a4a信信 信信信信 信 100101001 100101 1001 0110010100100
13、11信7五、 费诺仙农编码（P126）费诺仙农编码与 Huffman 编码相反，采用从上到下的方法。香农-范诺编码算法步骤：（1）按照符号出现的概率减少的顺序将待编码的符号排成序列。（2）将符号分成两组，使这两组符号概率和相等或几乎相等。（3）将第一组赋值为 0，第二组赋值为 1。（4）对每一组，重复步骤 2 的操作。例：设一副灰度级为 8 的图象中，各灰度所 对应的概率分别为 0.04，0.05，0.06，0.07，0.10，0.10，0.18，0.40，要求对其进行费诺.仙侬编码。灰度值 S0 S1 S2 S3 S4 S5 S6 S7出现频率 0.40 0.18 0.10 0.10 0.0
14、7 0.06 0.05 0.04解：根据费诺仙农编码的方法进行分组和赋值如下图所示s 0 , s 1 , s 2 , s 3 , s 4 , s 5 , s 6 , s 7s 2 , s 3 , s 4 , s 5 , s 6 , s 7s 0 , s 10 . 5 80 . 4 2s 2 , s 3 s 4 , s 5 , s 6 , s 7s 0 s 1s 4 , s 5 s 6 , s 7s 2 s 3s 4 s 5 s 60 . 2 00 . 2 20 . 1 3 0 . 0 90 10 1 01010 1010 s1所得编码结果如下表灰度值 S0 S1 S2 S3 S4 S5 S6
15、S7费诺仙农码 00 01 100 101 11001101111011118六、 算术编码（P127）例：编码来自 1 个 4-符号信源 a1, a2, a3, a4的由 5 个符号组成的符号序列：b 1b2b3b4b5 = a1a2a3a3a412a34信源符号 概率 初始子区间0 . 20 . 20 . 40 . 2 0 , 0 . 2 0 . 8 , 1 . 0 0 . 4 , 0 . 8 0 . 2 , 0 . 4 解：由（新子区间的起始位置=前子区间的起始位置+当前符号的区间左端*前子区间长度）LCFNls（新子区间的结束位置=前子区间的起始位置+当前符号的区间右端*前子区间长度）
16、rse可得，对于a 1，a2，a3，a3，a4，有a1 0，0.2a1a2 0.2*0.2，0.2*0.4=0.04，0.08a1a2a3 0.04+0.04*0.4，0.04+0.04*0.8=0.056，0.072a1a2a3a3 0.056+0.016*0.4，0.056+0.016*0.8=0.0624，0.0688a1a2a3a3a4 0.0624+0.0064*0.8，0.056+0.0064*1=0.06752，0.0688解码过程 0.06870234.( 1 ) 0 . 0 6 8 在 区 间 0 , 0 . 2 ， 可 知 第 一 个 源 符 号 为 a1( 2 )在 区
17、间 0 . 2 - 0 . 4 中 ， 第 二 个 为 a2( 3 )在 区 间 0 . 4 - 0 . 8 中 ， 第 三 个 为 a3( 4 )在 区 间 0 . 4 - 0 . 8 中 ， 第 四 个 为 a3( 5 )在 区 间 0 . 8 - 1 中 ， 第 五 个 为 a45.9七、 区域分割状态法（峰谷法、灰度阈值法） （P155）基本思想是，确定一个合适的阈值 T。将大于等于阈值的像素作为物体或背景，生成一个二值图像。阈值的选定可以通过如下图中灰度直方图确定。方法：首先统计最简单图像的灰度直方图，若直方图呈双峰且有明显的谷，则将谷所对应的灰度值 T 作为阈值，按图右侧的等式进行
18、二值化，就可将目标从图像中分割出来。这种方法适用于目标和背景的灰度差较大、有明显谷的情况。Tyxfyxg),(10),(在四邻域中有背景的像素，既是边界像素。例：对下面的图像用状态法进行二值化，并 计算二值图像的欧拉数。0 1 3 2 1 3 2 10 5 7 6 2 5 6 71 6 0 6 1 6 3 42 6 7 5 3 5 6 53 2 2 7 2 6 1 62 6 5 0 2 7 5 01 2 3 2 1 2 1 23 1 2 3 1 2 2 1解：（1）首先根据已知列出灰度级分布表灰度级 0 1 2 3 4 5 6 7像素个数 5 12 16 8 1 7 10 5（2）画出图像的直
19、方图02468101214160 1 2 3 4 5 6 710（3）由此可确定阈值 T = 4，根据 ，二值化的图像如下：Tyxfyxg),(10),(0 0 0 0 0 0 0 00 1 1 1 0 1 1 10 1 0 1 0 1 0 00 1 1 1 0 1 1 10 0 0 1 0 1 0 10 1 1 0 0 1 1 00 0 0 0 0 0 0 00 0 0 0 0 0 0 0（4）在二值图像中，1 像素连接成分数 C 减去孔数 H 的差值 叫做这幅图像的欧拉数。本题从图中可以看出，其取值为 1 的像素的连接成分数 C = 2，孔数 H = 2，所以这幅图像的欧拉数为E = C H = 2 2 = 0
