《数字图像处理计算题复习精华版要点》由会员分享,可在线阅读,更多相关《数字图像处理计算题复习精华版要点(15页珍藏版)》请在金锄头文库上搜索。
1、30452计算题复习一、 直方图均衡化(P68)对已知图像进行直方图均衡化修正。例:表1为已知一幅总像素为n=6464的8bit数字图像(即灰度级数为8),各灰度级(出现的频率)分布列于表中。要求将此幅图像进行均衡化修正(变换),并画出修正(变换)前后的直方图。表1原图像灰度级rk原各灰度级像素个数nk原分布概率pr(rk)r0=07900.19r1=110230.25r2=28500.21r3=36560.16r4=43290.08r5=52450.06r6=61220.03r7=7810.02解:对已知图像均衡化过程见下表:原图像灰度级rk原各灰度级像素个数nk原分布概率pr(rk)累积分
2、布函数sk计取整扩展sk并确定映射对应关系rksk新图像灰度级sk新图像各灰度级像素个数nsk新图像分布概率ps(sk)r0=07900.190.1910117900.19r1=110230.250.44313310230.25r2=28500.210.6552558500.21r3=36560.160.8163669850.24r4=43290.080.89646r5=52450.060.9575774480.11r6=61220.030.98767r7=7810.021.00777画出直方图如下: (a)原始图像直方图 (b)均衡化后直方图*以下部分不用写在答题中。其中: rk、nk中k
3、= 0,1,7 pr(rk)= nkn,即计算各灰度级像素个数占所有像素个数的百分比,其中,在此题中n=6464。 ,即计算在本灰度级之前(包含本灰度级)所有百分比之和。 ,其中L为图像的灰度级数(本题中L = 8),int 表示对方括号中的数字取整。 nsk 为映射对应关系rksk 中rk所对应的nk之和。 ,或为映射对应关系rksk 中rk所对应的pr(rk)之和。二、 模板运算 使用空间低通滤波法对图像进行平滑操作(P80)空间低通滤波法是应用模板卷积方法对图像每一个像素进行局部处理。模板(或称掩模)就是一个滤波器,它的响应为H(r,s),于是滤波输出的数字图像g(x,y)用离散卷积表示
4、为式中:x,y = 0,1,2,N-1;k、l根据所选邻域大小来决定。具体过程如下:(1)将模板在图像中按从左到右、从上到下的顺序移动,将模板中心与每个像素依次重合(边缘像素除外);(2)将模板中的各个系数与其对应的像素一一相乘,并将所有的结果相加;(3)将(2)中的结果赋给图像中对应模板中心位置的像素。对于空间低通滤波器而言,采用的是低通滤波器。由于模板尺寸小,因此具有计算量小、使用灵活、适于并行计算等优点。常用的3*3低通滤波器(模板)有: 模板不同,邻域内各像素重要程度也就不同。但无论怎样的模板,必须保证全部权系数之和为1,这样可保证输出图像灰度值在许可范围内,不会产生灰度“溢出”现象。
5、例:以为模板,对下图做低通滤波处理,写出处理结果。1718171111151111115551171155518181151111811511811115111117181711解:低通滤波的步骤为:(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(其他位置同样方法计算可得)由此步骤可得处理结果为(空白处自
6、己计算后填入)17181711121171181811117181711三、 中值滤波与邻域平均中值滤波(P81)中值滤波是对一个滑动窗口内的诸像素灰度值排序,用中间值代替窗口中心像素的原来灰度值,它是一种非线性的图像平滑法。它对脉冲干扰及椒盐噪声的抑制效果好,在抑制随机噪声的同时能有效保护边缘少受模糊。但它对点、线等细节较多的图像却不太合适。局部平滑法(邻域平均法 或 移动平均法)(P76)局部平滑法是一种直接在空间域上进行平滑处理的技术。用邻域内各像素的灰度平均值代替该像素原来的灰度值,实现图像的平滑。设有一幅NN的图像f(x,y),若平滑图像为g(x,y),则有 式中x,y = 0,1,
7、N-1; s为(x,y)邻域内像素坐标的集合; M表示集合s内像素的总数。可见邻域平均法就是将当前像素邻域内各像素的灰度平均值作为其输出值的去噪方法。设图像中的噪声是随机不相关的加性噪声,窗口内各点噪声是独立同分布的,经过(4.2.1)平滑后,信号与噪声的方差比可望提高M倍。这种算法简单,但它的主要缺点是在降低噪声的同时使图像产生模糊,特别在边缘和细节处。而且邻域越大,在去噪能力增强的同时模糊程度越严重。例:对下图做3*3中值滤波处理和3*3邻域平均处理,写出处理结果,并比较邻域平均与中值滤波的差异。17181711111511111155511711555181811511118115118
8、11115111117181711解:(1)中值滤波是对一个滑动窗口内的诸像素灰度值排序,用中间值代替窗口中心像素的原来灰度值,它是一种非线性的图像平滑法。题目中的图像经3*3中值滤波后的结果为(忽略边界):1718171111555111115551171155511181151111811111111111111117181711(2)局部平滑法(邻域平均法 或 移动平均法)是用邻域内各像素的灰度平均值代替该像素原来的灰度值,实现图像的平滑。题目中的图像经3*3局部平滑法(邻域平均法 或 移动平均法)后的结果为(忽略边界):17181711119/938/940/938/923/921/9
9、1171181811117181711中值滤波法和局部平滑法(邻域平均法 或 移动平均法)均能有效削弱椒盐噪声,但中值滤波法比邻域平均法更有效,且滤波后图像中的轮廓比较清晰。四、 霍夫曼编码(P124)例:设有一信源A=a1, a2, a3, a4, a5, a6,对应概率P=0.1, 0.4, 0.06, 0.1, 0.04, 0.3.(1)进行霍夫曼编码(要求大概率的赋码字0,小概率的赋码字1),给出码字;(2)计算平均码长,信源熵和编码效率。解:(1)编码步骤1)缩减信源符号数量将信源符号按出现概率从大到小排列,然后结合2)对每个信源符号赋值从(消减到)最小的信源开始,逐步回到初始信源由
10、此可得哈夫曼编码结果见下表符号a1a2a3a4a5a6概率0.10.40.060.10.040.06编码结果01110101001000101100平均码长(其中,是灰度值为i的编码长度,为灰度值为的概率,L为灰度级数)信源熵编码效率五、 费诺仙农编码(P126)费诺仙农编码与Huffman编码相反,采用从上到下的方法。香农-范诺编码算法步骤:(1)按照符号出现的概率减少的顺序将待编码的符号排成序列。(2)将符号分成两组,使这两组符号概率和相等或几乎相等。(3)将第一组赋值为0,第二组赋值为1。(4)对每一组,重复步骤2的操作。例:设一副灰度级为8的图象中,各灰度所对应的概率分别为0.04,0
11、.05,0.06,0.07,0.10,0.10,0.18,0.40,要求对其进行费诺.仙侬编码。灰度值S0S1S2S3S4S5S6S7出现频率0.400.180.100.100.070.060.050.04解:根据费诺仙农编码的方法进行分组和赋值如下图所示所得编码结果如下表灰度值S0S1S2S3S4S5S6S7费诺仙农码00011001011100110111101111六、 算术编码(P127)例:编码来自1个4-符号信源a1, a2, a3, a4的由5个符号组成的符号序列:b1b2b3b4b5 = a1a2a3a3a4解:由(新子区间的起始位置=前子区间的起始位置+当前符号的区间左端*前
12、子区间长度)(新子区间的结束位置=前子区间的起始位置+当前符号的区间右端*前子区间长度)可得,对于a1,a2,a3,a3,a4,有a10,0.2a1a20.2*0.2,0.2*0.4=0.04,0.08a1a2a30.04+0.04*0.4,0.04+0.04*0.8=0.056,0.072a1a2a3a30.056+0.016*0.4,0.056+0.016*0.8=0.0624,0.0688a1a2a3a3a40.0624+0.0064*0.8,0.056+0.0064*1=0.06752,0.0688解码过程 0.068七、 区域分割状态法(峰谷法、灰度阈值法)(P155)基本思想是,确定一个合适的阈值T。将大于等于阈值的像素作为物体或背景,生成一个二值图像。阈值的选定可以通过如下图中灰度直方图确定。方法:首先统计最简单图像的