《数字图像处理计算题复习精华版要点.pdf》由会员分享,可在线阅读,更多相关《数字图像处理计算题复习精华版要点.pdf(19页珍藏版)》请在三一文库上搜索。
1、1 30452 计算题复习 一、直方图均衡化( P68) 对已知图像进行直方图均衡化修正。 例:表 1 为已知一幅总像素为n=64 64 的 8bit 数字图像(即灰度级数为8) ,各灰度级(出现的频率)分布 列于表中。要求将此幅图像进行均衡化修正(变换),并画出修正(变换)前后的直方图。 表 1 原图像灰 度级rk 原各灰度级 像素个数nk 原分布概率 pr(rk) r0=0 790 0.19 r1=1 1023 0.25 r2=2 850 0.21 r3=3 656 0.16 r4=4 329 0.08 r5=5 245 0.06 r6=6 122 0.03 r7=7 81 0.02 解:
2、对已知图像均衡化过程见下表: 原图像灰 度级rk 原各灰度级 像素个数nk 原分布概率 pr(rk) 累积分布函 数sk计 取整扩展 sk并 确定映射 对应关系 rksk 新图像灰 度级sk 新图像各灰 度级像素个 数nsk 新图像分 布概率 ps(sk) r0=0 790 0.19 0.19 1 011 790 0.19 r1=1 1023 0.25 0.44 3 133 1023 0.25 r2=2 850 0.21 0.65 5 255 850 0.21 r3=3 656 0.16 0.81 6 36 6 985 0.24 r4=4 329 0.08 0.89 6 46 r5=5 245
3、 0.06 0.95 7 57 7 448 0.11 r6=6 122 0.03 0.98 7 67 r7=7 81 0.02 1.00 7 77 画出直方图如下: 2 (a)原始图像直方图(b)均衡化后直方图 * 以下部分不用写在答题中。 其中: rk 、n k中 k = 0,1,7 pr(rk)= nkn, 即计算各灰度级像素个数占所有像素个数的百分比,其中 k j j nn 0 ,在此题中 n=64 64。 k j jrk rps 0 )( 计 ,即计算在本灰度级之前(包含本灰度级)所有百分比之和。 5.0)1int( 计并kk sLs ,其中 L 为图像的灰度级数(本题中L = 8)
4、, int 表示对方括号中的数字取 整。 并kk ss nsk为映射对应关系 rksk 中 rk所对应的nk之和。 nnsp skks /)( ,或为映射对应关系 rksk中rk所对应的pr(rk)之和。 3 二、模板运算使用空间低通滤波法对图像进行平滑操作(P80) 空间低通滤波法是应用模板卷积方法对图像每一个像素进行局部处理。模板(或称掩模)就是一个滤 波器,它的响应为H(r,s),于是滤波输出的数字图像g(x,y)用离散卷积表示为 )6.2.4(),(),(),( l ls k kr srHsyrxfyxg 式中: x,y = 0,1,2,N-1;k、 l 根据所选邻域大小来决定。 具体
5、过程如下: (1)将模板在图像中按从左到右、从上到下的顺序移动,将模板中心与每个像素依次重合(边缘像素除外) ; (2)将模板中的各个系数与其对应的像素一一相乘,并将所有的结果相加; (3)将( 2)中的结果赋给图像中对应模板中心位置的像素。 对于空间低通滤波器而言,采用的是低通滤波器。由于模板尺寸小,因此具有计算量小、使用灵活、 适于并行计算等优点。常用的3*3 低通滤波器(模板)有: 111 111 111 9 1 1 H 111 121 111 10 1 2 H 121 242 121 16 1 3H 111 101 111 8 1 4 H 00 1 00 2 1 4 1 4 1 4 1
6、 4 1 5 H 模板不同,邻域内各像素重要程度也就不同。但无论怎样的模板,必须保证全部权系数之和为1,这样可保 证输出图像灰度值在许可范围内,不会产生灰度“溢出”现象。 例:以 111 121 111 10 1 2 H为模板,对下图做低通滤波处理,写出处理结果。 1 7 1 8 1 7 1 1 1 1 1 5 1 1 1 1 1 1 5 5 5 1 1 7 1 1 5 5 5 1 8 1 8 1 1 5 1 1 1 1 8 1 1 5 1 1 8 1 1 1 1 5 1 1 1 1 1 7 1 8 1 7 1 1 解:低通滤波的步骤为: (1)将模板在图像中按从左到右、从上到下的顺序移动,将
7、模板中心与每个像素依次重合(边缘像素除外) ; (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 1 1 2 1 4 1 7 1 1 8 1 8 1 1 1 1 7 1 8 1 7 1 1 5 三、中值滤波与邻域平均 中值滤波( P81) 中值滤波是对一个滑动窗口内的诸像素灰度值排序,用中
8、间值代替窗口中心像素的原来灰度值,它是 一种非线性的图像平滑法。 它对脉冲干扰及椒盐噪声的抑制效果好,在抑制随机噪声的同时能有效保护边缘少受模糊。但它对点、 线等细节较多的图像却不太合适。 局部平滑法(邻域平均法或 移动平均法) (P76) 局部平滑法是一种直接在空间域上进行平滑处理的技术。用邻域内各像素的灰度平均值代替该像素原来的 灰度值,实现图像的平滑。 设有一幅 NN 的图像 f(x,y),若平滑图像为g(x,y),则有 )1.2 .4(),( 1 ),( ,sji jif M yxg 式中 x,y = 0,1, ,N-1; s为( x,y)邻域内像素坐标的集合; M 表示集合s内像素的
9、总数。 可见邻域平均法就是将当前像素邻域内各像素的灰度平均值作为其输出值的去噪方法。 设图像中的噪声是随机不相关的加性噪声,窗口内各点噪声是独立同分布的,经过(4.2.1)平滑后, 信号与噪声的方差比可望提高M 倍。 这种算法简单,但它的主要缺点是在降低噪声的同时使图像产生模糊,特别在边缘和细节处。而且邻 域越大,在去噪能力增强的同时模糊程度越严重。 例:对下图做3*3 中值滤波处理和3*3 邻域平均处理,写出处理结果,并比较邻域平均与中值滤波的差异。 1 7 1 8 1 7 1 1 1 1 1 5 1 1 1 1 1 1 5 5 5 1 1 7 1 1 5 5 5 1 8 1 8 1 1 5
10、 1 1 1 1 8 1 1 5 1 1 8 1 1 1 1 5 1 1 1 1 1 7 1 8 1 7 1 1 解: (1)中值滤波是对一个滑动窗口内的诸像素灰度值排序,用中间值代替窗口中心像素的原来灰度值,它是 一种非线性的图像平滑法。 题目中的图像经3*3 中值滤波后的结果为(忽略边界): 1 7 1 8 1 7 1 1 1 1 5 5 5 1 1 1 1 1 5 5 5 1 1 7 1 1 5 5 5 1 1 1 6 8 1 1 5 1 1 1 1 8 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 7 1 8 1 7 1 1 (2)局部平滑法 (邻域平均法或 移动平均
11、法) 是用邻域内各像素的灰度平均值代替该像素原来的灰度值, 实现图像的平滑。 题目中的图像经3*3 局部平滑法(邻域平均法或 移动平均法)后的结果为(忽略边界): 1 7 1 8 1 7 1 1 1 19/9 38/9 40/9 38/9 23/9 21/9 1 1 7 1 1 8 1 8 1 1 1 1 7 1 8 1 7 1 1 中值滤波法和局部平滑法(邻域平均法或 移动平均法)均能有效削弱椒盐噪声,但中值滤波法比邻 域平均法更有效,且滤波后图像中的轮廓比较清晰。 7 四、霍夫曼编码( P124) 例:设有一信源A=a1, a2, a3, a4, a5, a6 ,对应概率P=0.1, 0.
12、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 0100 01011 00 0.4 0.3 0.1 0.1 0.06 0.04 0.4 0.3 0.1 0.1 0
13、.1 0.4 0.3 0.1 0.2 0.4 0.3 0.3 0.4 0.6 1234 a1 2 a 5 a 3 a 6 a 4 a 初始信源信源的消减步骤 符号概率 0.4 0.3 0.1 0.1 0.06 0.04 0.4 0.3 0.1 0.1 0.1 0.4 0.3 0.1 0.2 0.4 0.3 0.3 0.4 0.6 1234 a1 2 a 5 a 3 a 6 a 4 a 初始信源对消减信源的赋值 符号概率 1 00 011 0100 0101 1 00 010 011 1 00 01 0 1 1 00 011 0100 01010 01011 码字 8 平均码长2.2504.0.
14、0506.0.041 .031 .023 .014.0 1 0 L i iip B (其中, i是灰度值为 i 的编码长度, i p为灰度值为 i a的概率, L 为灰度级数) 信源熵14.2log 2 1 0 i L i i ppH 编码效率973.0 2 . 2 14.2 B H 9 五、费诺仙农编码( P126) 费诺仙农编码与Huffman 编码相反,采用从上到下的方法。香农-范诺编码算法步骤: (1)按照符号出现的概率减少的顺序将待编码的符号排成序列。 (2)将符号分成两组,使这两组符号概率和相等或几乎相等。 (3)将第一组赋值为0,第二组赋值为1。 (4)对每一组,重复步骤2 的操
15、作。 例:设一副灰度级为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.07 0.06 0.05 0.04 解:根据费诺 仙农编码的方法进行分组和赋值如下图所示 s0,s1,s2,s3,s4,s5,s6,s7 s2,s3,s4,s5,s6,s7s0,s1 0.58 0.42 s2,s3s4,s5,s6,s7s0s1 s4,s5s6,s7 s2s3 s4s5s6 0.20 0.22 0
16、.130.09 01 010 1 0 1 01 01 0 s7 1 所得编码结果如下表 灰度值 S0 S1 S2 S3 S4 S5 S6 S7 费诺 仙农码 00 01 100 101 1100 1101 1110 1111 10 六、算术编码( P127) 例:编码来自1 个 4-符号信源 a1, a2, a3, a4的由 5 个符号组成的符号序列:b1b2b3b4b5 = a1a2a3a3a4 1 a 2 a 3 a 4 a 信源符号概率初始子区间 0.2 0.2 0.4 0.2 0 , 0.2 0.8 , 1.0 0.4 , 0.8 0.2 , 0.4 解:由 LCFN lss *(新子
17、区间的起始位置=前子区间的起始位置+当前符号的区间左端*前子区间长度) LCFN rse *(新子区间的结束位置=前子区间的起始位置+当前符号的区间右端*前子区间长度) 可得,对于 a1,a2,a3, a3,a4 ,有 a10 , 0.2 a1a2 0.2*0.2, 0.2*0.4=0.04, 0.08 a1a2a3 0.04+0.04*0.4, 0.04+0.04*0.8=0.056, 0.072 a1a2a3a3 0.056+0.016*0.4, 0.056+0.016*0.8=0.0624, 0.0688 a1a2a3a3a4 0.0624+0.0064*0.8, 0.056+0.006
18、4*1=0.06752, 0.0688 解码过程0.068 70 20 20340 . . (1)0.068 在区间0 ,0.2 ,可知第一个源符号为a1 (2)在区间 0.2-0.4 中,第二个为a2 (3)在区间 0.4-0.8 中,第三个为a3 (4) 在区间 0.4-0.8 中,第四个为a3 (5)在区间 0.8-1 中,第五个为a4 340 20 00680 . . . 750 40 4070 . . 8750 40 40750 . . 11 七、区域分割 状态法(峰谷法、灰度阈值法)(P155) 基本思想是,确定一个合适的阈值T。将大于等于阈值的像素作为物体或背景,生成一个二值图像
19、。阈值的 选定可以通过如下图中灰度直方图确定。 方法:首先统计最简单图像的灰度直方图,若直方图呈双峰且有明显的谷,则将谷所对应的灰度值T 作为 阈值,按图右侧的等式进行二值化,就可将目标从图像中分割出来。这种方法适用于目标和背景的灰度差 较大、有明显谷的情况。 Tyxf Tyxf yxg ),(1 ),(0 ),( 在四邻域中有背景的像素,既是边界像素。 例:对下面的图像用状态法进行二值化,并计算二值图像的欧拉数。 0 1 3 2 1 3 2 1 0 5 7 6 2 5 6 7 1 6 0 6 1 6 3 4 2 6 7 5 3 5 6 5 3 2 2 7 2 6 1 6 2 6 5 0 2
20、7 5 0 1 2 3 2 1 2 1 2 3 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)画出图像的直方图 12 0 2 4 6 8 10 12 14 16 01234567 (3)由此可确定阈值T = 4,根据 Tyxf Tyxf yxg ),(1 ),(0 ),(,二值化的图像如下: 0 0 0 0 0 0 0 0 0 1 1 1 0 1 1 1 0 1 0 1 0 1 0 0 0 1 1 1 0 1 1 1 0 0 0 1 0 1 0 1 0 1 1 0 0 1 1 0
21、0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 (4)在二值图像中,1 像素连接成分数C 减去孔数 H 的差值叫做这幅图像的欧拉数。 本题从图中可以看出,其取值为1 的像素的连接成分数C = 2,孔数 H = 2,所以这幅图像的欧拉数为 E = C H = 2 2 = 0 13 八、区域增长 简单区域扩张法(P159) 步骤:以图像的某个像素为生长点,比较相邻像素的特征,将特征相似的相邻像素合并为同一区域;以合 并的像素为生长点,继续重复以上的操作,最终形成具有相似特征的像素是最大连通集合。这种方法称简 单(单一型)区域扩张法。 步骤: (1)从图像最左上角开始,对图像进行光栅
22、扫描,找到不属于任何的像素。 (2)把这个像素灰度同其周围(4 邻域或 8 邻域)不属于其他区域的像素的灰度值和已存在区域的像素灰 度平均值进行比较,若灰度差值小于阈值,则合并到同一区域,并对合并的像素赋予标记。 (3)从新合并的像素开始,反复进行(2)的操作。 (4)反复进行(2) 、 (3)的操作,直至不能再合并。 (5)返回( 1)操作,寻找新区域出发点的像素。 例:对下面的图像采用简单区域生长法进行区域生长,给出灰度差值 T = 1;T = 3 ;T = 8 三种情况 下的分割图像。 1 0 4 7 5 1 0 4 7 7 0 1 5 5 5 2 0 5 6 5 0 2 5 6 4 解
23、:以 8 邻域为比较范围, 当阈值 T = 1 时,图像分割结果如下图所示6 个区域: 1 0 4 7 5 1 0 4 7 7 0 1 5 5 5 2 0 5 6 5 0 2 5 6 4 当阈值T = 3 时,图像分割结果如下图所示为2 个区域 1 0 4 7 5 1 0 4 7 7 0 1 5 5 5 2 0 5 6 5 0 2 5 6 4 当阈值T = 8 时,图像分割结果如下图所示为1 个区域 1 0 4 7 5 1 0 4 7 7 0 1 5 5 5 2 0 5 6 5 14 0 2 5 6 4 15 九、分裂合并 分裂合并法(基于四叉树思想的方法)(P161) 算法实现: 1)对于图
24、像中灰度级不同的区域,均分为四个子区域。 2)如果相邻的子区域所有像素的灰度级相同,则将其合并。 3)反复进行上两步操作,直至不再有新的分裂与合并为止。 例:用分裂合并法分割图像,并给出对应分割结果的四叉树。 解:四叉树法分裂结果如下: 00 01 1 030 031 02 033 032 300 301 31 20 21 303 302 33 32 23 22 对应四叉树为: node 0123 00010203 300301302303 2021222330313233 030031032033 合并后结果为: 16 其中背景区域 包括四叉树中00、01、030、033、1、 21、22、
25、23、32、33、300、303 子块 图像区域 包括四叉树中031、032、02、20、 301、302、31 子块 17 十、膨胀和腐蚀、开运算和闭运算(P172) 1)膨胀 膨胀就是把二值图像各1 像素连接成分的边界扩大一层的处理。 膨胀的原理 :设二值图像为F,结构元素为B,Bs代表 B 关于原点对称的结构元素。当结构元素B s的原点 移到( x,y)处时,结构元素用 s xy B表示。则图像F 被结构元素B 膨胀的定义式为: FBBF s xy yx|,(即 B s 击中 F,均指各自的非零像素) 其含义是 :当结构元素Bs的原点移动到 (x,y)位置时, 如果 s xy B所覆盖范
26、围内的F 的子图像与结构元素 s xy B 相应位置上至少有一个元素相同且不为0,则把该子图像中与 s xy B的原点位置对应的(x,y)点的那个像素位 置标注为1,否则为0。图像 F 上标注出的所有这样的像素组成的集合,即为膨胀运算的结果。 膨胀运算的基本过程是: (a)求结构元素B 关于其原点的反射集合Bs; (b)每当结构元素在目标图像F 上平移后,结构元素B s 与其覆盖的子图像中至少有一个元素相交时,就 将目标图像中与结构元素Bs的原点对应的那个位置的像素值置为“1”,否则置为0。 注意 : (a)结构元素中原点位置所对应的目标图像子图像位置处的值是0 时,仍可进行膨胀运算,无需强求
27、是1。 (b)当结构元素在目标图像上平移时,允许结构元素中的非原点像素超出目标图像范围。 结构元素形状对膨胀运算结果的影响:当目标图像不变,但所给的结构元素的形状改变时;或结构元素 的形状不变,而其原点位置改变时,膨胀运算的结果会发生改变。 2)腐蚀(或收缩) 腐蚀是把二值图像各1 像素连接成分的边界点去掉从而缩小一层的处理。 腐蚀的原理 :设 F 为目标图像, B 为结构元素,则目标图像F 被结构元素B 腐蚀可定义为 FBBF xy yx|,(即 B 包含于 F,均指各自的非零像素) 其含义是 :当结构元素B 的原点移动到目标图像F 中的( x,y)位置时,如果(x,y)处像素值为1,并且
28、Bxy 所覆盖范围内的F 的子图像的其他像素能够包含Bxy的其他像素或与Bxy的其他像素完全相同,则保留该子 图像中与 Bxy的原点位置对应的(x,y)点的像素值1,否则均为0。图像 F 上保留的所有这样值为1 的像素 组成的集合,即为腐蚀运算的结果。 这里的 “ 包含 ” 是指结构元素B 和目标图像F 的子图像中值为1 的像素两两之间的对应关系。 腐蚀运算的基本过程是: 把结构元素B 看作为一个卷积模板,每当结构元素平移到其原点位置与目标图像F 中那些像素值为“1” 的位置重合时,就判断被结构元素覆盖的子图像的其它像素的值是否都与结构元素相应位置的像素值相同; 只有当其都相同时,就将结果图像
29、中的那个与原点位置对应的像素位置的值置为“1”,否则置为0。 注意 : (a)结构元素中的原点位置处的像素值可以不为1,但要求目标图像中的子图像与结构元素B 的原点对应 的那个位置的像素值是1。 (b)当结构元素在目标图像上平移时,结构元素中的任何元素不能超出目标图像的范围。 腐蚀运算的结果不仅与结构元素的形状(矩形、圆形、菱形等 )选取有关, 而且还与原点位置的选取有关。 3)膨胀和腐蚀的作用: 膨胀的作用是使孔洞收缩,目标扩大。对消除图像目标中的小颗粒噪声和填补凹陷非常有效。 腐蚀的左右是使目标收缩,孔洞扩大。对去除图像小颗粒噪声和目标之间的粘连非常有效。 4)开运算V.S.闭运算 A、开运算(先腐蚀再膨胀):BBFBF)( 作用:光滑目标轮廓、消除小目标(如去掉毛刺和孤立点等),在纤细点处分离物体,同时并不明显改变目 18 标面积; B、闭运算(先膨胀再腐蚀):BBFBF)( 作用:在保持原目标的大小与形态的同时,填充凹陷、弥合孔洞和裂缝。 例 1:用结构元素B 对目标图像F 进行腐蚀运算。结构元素B 中红色为原点。 例 2:结构元素不同时的腐蚀运算 例 3:结构元素原点不同时的腐蚀运算 例 4:用结构元素B 对目标图像F 进行膨胀运算。结构元素B 中红色为原点。 19 图 8.10 膨胀运算实例 例 5 例 6
链接地址:https://www.31doc.com/p-5208699.html