数字图像处理：部分课后习题参考答案

第PAGE\*MERGEFORMAT#页（共13页）第PAGE\*MERGEFORMAT#页（共13页）第1页(共13页)第一章1.连续图像中，图像为一个二维平面，(x,y)图像中的任意一点，f(x,y)为图像于(x,y)于处的值连续图像中，(x,y)的取值是连续的，f(x,y)也是连续的数字图像中，图像为一个由有限行有限列组成的二维平面，(i,j)为平面中的任意一点，g(i,j)则为图像在(i,j)处的灰度值，数字图像中，(i,j)的取值是不连续的，只能取整数，对应第i行j列,g(i,j)也是不连续的，表示图像i行j列处图像灰度值。联系：数字图像g(i,j)是对连续图像f(x,y)经过采样和量化这两个 步骤 新产品开发流程的步骤课题研究的五个步骤成本核算步骤微型课题研究步骤数控铣床操作步骤 得到的。其中g(i,j)=f(x,y)|x=i,y=j2.图像工程的内容可分为图像处理、图像分析和图像理解三个层次，这三个层次既有联系又有区别，如下图所示。抽象程度数据量图像处理的重点是图像之间进行的变换。尽管人们常用图像处理泛指各种图像技术，但比较狭义的图像处理主要是对图像进行各种加工，以改善图像的视觉效果并为自动识别奠定基础，或对图像进行压缩编码以减少所需存储空间图像分析主要是对图像中感兴趣的目标进行检测和测量，以获得它们的客观信息，从而建立对图像的描述。如果说图像处理是一个从图像到图像的过程，则图像分析是一个从图像到数据的过程。这里的数据可以是目标特征的测量结果，或是基于测量的符号表示，它们描述了目标的特点和性质。图像理解的重点是在图像分析的基础上，进一步研究图像中各目标的性质和它们之间的相互联系，并得出对图像内容含义的理解以及对原来客观场景的解释，从而指导和规划行动。如果说图像分析主要以观察者为中心来研究客观世界，那么图像理解在一定程度上是以客观世界为中心，借助知识、经验等来把握整个客观世界(包括没有直接观察到的事物)的。联系：图像处理、图像分析和图像理解处在三个抽象程度和数据量各有特点的不同层次上。图像处理是比较低层的操作，它主要在图像像素级上进行处理，处理的数据量非常大。图像分析则进入了中层，分割和特征提取把原来以像素描述的图像转变成比较简洁的非图形式的描述。图像理解主要是高层操作，基本上是对从描述抽象出来的符号进行运算，其处理过程和 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 与人类的思维推理有许多类似之处。第二章：1.参见第一章：第1题jpg（jpeg）JPEG联合摄影专家组JPEG图片以24位颜色存储单个光栅图像。JPEG是与平台无关的格式，支持最高级别的压缩，不过，这种压缩是有损耗的。渐近式JPEG文件支持交错。可以提高或降低JPEG文件压缩的级别。但是，文件大小是以图像质量为代价的。压缩比率可以高达100:1。（JPEG格式可在10:1到20:1的比率下轻松地压缩文件，而图片质量不会下降。）JPEG压缩可以很好地处理写实摄影作品。但是，对于颜色较少、对比级别强烈、实心边框或纯色区域大的较简单的作品，JPEG压缩无法提供理想的结果。有时，压缩比率会低到5:1，严重损失了图片完整性。这一损失产生的原因是，JPEG压缩 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 可以很好地压缩类似的色调，但是JPEG压缩方案不能很好地处理亮度的强烈差异或处理纯色区域。优点：摄影作品或写实作品支持高级压缩，利用可变的压缩比可以控制文件大小。支持交错（对于渐近式JPEG文件）。JPEG广泛支持Internet 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 。缺点：有损耗压缩会使原始图片数据质量下降。当您编辑和重新保存JPEG文件时，JPEG会混合原始图片数据的质量下降。这种下降是累积性的。不适用于所含颜色很少、具有大块颜色相近的区域或亮度差异十分明显的较简单的图片。是最常见的格式之一。BMP：Windows位图Windows位图可以用任何颜色深度（从黑白到24位颜色）存储单个光栅图像。Windows位图文件格式与其他MicrosoftWindows程序兼容。它不支持文件压缩，也不适用于Web页。从总体上看，Windows位图文件格式的缺点超过了它的优点。为了保证照片图像的质量，请使用PNG文件、JPEG文件或TIFF文件。BMP文件适用于Windows中的墙纸。优点：BMP支持1位到24位颜色深度。BMP格式与现有Windows程序（尤其是较旧的程序）广泛兼容。缺点：BMP不支持压缩，这会造成文件非常大，BMP文件不受Web浏览器支持。GIF:图形交换格式；GIF图片以8位颜色或256色存储单个光栅图像数据或多个光栅图像数据。GIF图片支持透明度、压缩、交错和多图像图片（动画GIF）。PGIF透明度不是alpha通道透明度，不能支持半透明效果。GIF压缩是LZW压缩，压缩比大概为3:1。GIF文件 规范 编程规范下载gsp规范下载钢格栅规范下载警徽规范下载建设厅规范下载 的GIF89a版本中支持动画GIF。优点：GIF广泛支持Internet标准。支持无损耗压缩和透明度。动画GIF很流行，易于使用许多GIF动画程序创建。很多QC表情都是GIF的〜缺点：GIF只支持256色调色板，因此，详细的图片和写实摄影图像会丢失颜色信息第四章1.如下表：Pr（rk）、Pz（Zl）分别表示原直方图与规定直方图—、对原直方图进行均衡化处理，得到映射关系rk—Sk（第四列）二、对规定直方图进行均衡化处理，得到映射关系乙—Vl（第五列）:2.第PAGE\*MERGEFORMAT#页（共13页）2.第PAGE\*MERGEFORMAT#页（共13页）rkPr(rk)Pz(Zi)Sk=T(rk)vk=G(zi)rk=>zi规定化后的直方图r00.1400.141/70r0=>Z30n0.2200.363/70r1=>Z40r20.2600.624/70r2=>Z50ra0.170.190.796/70.19r3=>Z60.14r40.090.250.886/70.44r4=>Z60.22r50.060.210.947/70.65r5=>Z70.26r60.040.240.987/70.89r6=>Z70.26r70.020.111.007/71r7=>Z70.12三、对于每个Sk,迭代计算出满足下式的最小V1，得到映射关系Sk—V1,再由九―Sk得到「k—V1最后由Zi—V1的逆变换V1—乙求出rk—乙的变换：lkV|SkPz(Zj)Pr(rJ0j0对k=0,l=3时,v3-so=O.19-O.14>=O,对k=1,l=3时,v4-si=0.44-0.36>=0,对k=2,l=3时,v5-S2=0.65-0.62>=0,对k=3,l=3时,v6-s3=0.89-0.79>=0,对k=4,l=3时,v6-s4=0.89-0.88>=0,对k=5,l=3时,v7-s5=1.00-0.94>=0,对k=6,l=3时,v7-s6=1.00-0.98>=0,对k=6,l=3时,v7-s7=1.00-1.00>=0,i0开始满足上式，于是有r0=>Z3开始满足上式，于是有r1=>Z4开始满足上式，于是有r2=>Z5开始满足上式，于是有r3=>Z6开始满足上式，于是有r4=>Z6开始满足上式，于是有r5=>Z7开始满足上式，于是有r6=>Z7开始满足上式，于是有r7=>Z7直方图规定化结果如最后一列所示，规定化后的直方图与规定直方图基本一致第PAGE\*MERGEFORMAT#页（共13页）第PAGE\*MERGEFORMAT#页（共13页）rknkPr(「k)Sk=T(rk)rk=>SkPs(Sk)ro5600.140.141/70=>10ri9200.220.363/71=>30.140.143r210460.260.624/72=>40.286ra7050.170.796/73=>60.220.429r43560.090.886/74=>60.260.571r52670.070.947/70.714r61700.040.997/75,6,7->70.260.857r7720.021.007/70.13111111225523333325543333325546334525582346781111112585830233598759323331598860633326062823467811111122552333332554333332554633452558234678第4页原图111111223333333333334446334578页2页）346783*3中值滤波原图3*3领域平均法32*m+15•均值滤波：把每个像素都用周围的8个像素来做均值操作。可以平滑图像，速度快，算法简单。但是无法去掉噪声，这能微弱的减弱它。中值滤波：常用的非线性滤波方法，也是图像处理技术中最常用的预处理技术。它在平滑脉冲噪声方面非常有效，同时它可以保护图像尖锐的边缘。为什么一般情况下对离散图象的直方图均衡化并不能产生完全平坦的直方图？由于离散图象的直方图也是离散的，其灰度的累积分布函数是一个不减的阶梯函数。如果映射后的图象仍能取到所有256级灰度，那一定是原图象没有任何改变，这种情况只可能发生在原图象的直方图已经是一条水平线的情况下。一般情况下映射后所得到的图象只能取到少于256级灰度，这样在变换后的直方图中会有某些灰度级空缺，当然这些空缺应该均匀分布在0到255之间。于是问题就变成了将原有的256个值，即各灰度的概率，按顺序分成n（n<256）份，每份的概率总和应该相等。显然这个问题是不一定有解的，因此我们只能找到一个近似解。其结果就是最后得到一幅有空缺且不太平坦的直方图不会发生变化，因为再次均衡化，所用的变换函数是首次均衡化后得到的增强图像的累积直方图，不会改变其结果。10.相同点：都能减弱或消除傅立叶空间的某些分量，而不影响或较少影响其他分量，从而达到了增强某些频率分量的效果。不同点：平滑滤波器减弱或消除了傅立叶空间的高频分量，所以达到了增强低频分量，平滑图像中细节的效果。锐化滤波器减弱或消除了傅立叶空间的低频分量，所以达到了增强高频分量，锐化图像中细节的效果。:两者效果相反，互为补充，从原始图像中减去平滑滤波器的结果得到锐化滤波器的效果，而从原始图像中减去锐化滤波器的结果则可得到平滑滤波器的结果。第五章：P1051.0000100011001-21%-20-2$0-2(>()21010(10(J1参见教材P83第六章P1412•—般地，按比例将原图像放大k倍时，如果按照最近邻域法则需要将一个像素值添在新图像的kxk的子块，如果放大倍数太大，按照这种方法处理会出现马赛克效应。为了提高几何变换后的图像质量，常采用线性插值法。该方法的原理是，当求出的分数地址与像素点不一致时，求出周围四个像素点的距离比，根据该比率，由四个邻域的像素灰度值进行线性插值图像旋转之后，会出现许多的空洞点。对这些空洞点必须进行填充处理，否则画面效果不好一般也称这种操作为插值处理。最简单的方法是行插值方法或列插值方法：列插值算法如下：找出当前列的最小和最大的非白点的坐标，记作（k1，j）、（k2,j）。在（k1，k2）范围内进行插值，插值的方法是：空点的像素值等于上一点的像素值。同样的操作重复到所有列。经如上的插值处理之后，图像效果就变得自然了变换矩阵：T=0.70700.70700-0.70700.70700001.0000本题图像共有16个像素，变换前的3X16矩阵如下：P0=123412334414111222111111111变换后的3X16矩阵：P=T*P0516234234100-1-1-210-1111111111变换结果如下图所示4123412233334411111114344544-1110-1211111111Xx坐标012340123456第7页(共13页)第7页(共13页)6.f(221,396)=18,f(221,397)=45,f(222,396)=52,f(222,397)=36,试分别用最邻近插值法和双线性插值法，分别计算f(221.3,396.7)的值.解：设已知点(221.3,396.7)的周围像素的灰度值，用最邻近插值法，求点(221.3,396.7)的灰度值，•••221.3-221<222-221.3且396.7-396>397-396.7即所求点离点(221.397)最近•••f(221.3，396.7)=f(221，397)=457.2.双线性插值法，设x,y为所求点至点（221,396）的x,y坐标增量，如图所示:f(x,396)=f(221,396)+x*(f(222,396)-f(221,396))=18+34*xf(x,397)=f(221,397)+x*(f(222,397)-f(221,397))=45-9*xf(x,y)=f(x,396)+y*(f(x,397)-f(x,396))=18+34*x+y*(45-9*x-18-34*x)=18+34x+27y-43xy•••f(0.3,0.7)=38(221,396)(221(221,397)插值点•7)(222,397)首先将原点平移到（100,260)A=-100-260然后旋转B=cos60-sin60sin60cos60然后在平移回来100260以上变换为复合变换矩阵T=C*B*A注意是用的齐次坐标[x,y,1]'=T[X0Y01]'第七章：1、图像的频域处理就是把图像从空间域变换到频域，分析图像的频谱特性，据此进行图像处理；它的理论基础是：“任何波形都可以用单纯的正弦波的加权和表示”第8页（共13页）第PAGE\*MERGEFORMAT#页(共13页)第PAGE\*MERGEFORMAT#页（共13页）第PAGE\*MERGEFORMAT#页（共13页）常用变换：①傅里叶变换：它是应用最广泛和最重要的变换。它的变换核是复指数函数，转换域图像是原空间域图像的二维频谱，其“直流”项与原图像亮度的平均值成比例，高频项表征图像中边缘变化的强度和方向。为了提高运算速度，计算机中多采用傅里叶快速算法。②沃尔什-哈达玛变换：它是一种便于运算的变换。变换核是值+1或-1的有序序列。这种变换只需要作加法或减法运算,不需要象傅里叶变换那样作复数乘法运算,所以能提高计算机的运算速度，减少存储容量。其他还有余弦变换、正弦变换等也在图像处理中得到应用3、不管是连续傅立叶变换还是离散傅立叶变换，变换域均反映了被变换域的频谱。不同：1）连续傅立叶变换，信号量和自变量均是连续的，而离散傅立叶变换，自变量和信号量均是离散的。2）连续傅立叶变换，信号可以是无限长的，信号量也可以是无穷大；而对于离散傅立叶变换，信号应该是有限长的，信号量也应该是有限值，才能用计算机进行处理。5、根据二维离散傅立叶变换的公式，有：33j2(ux/Mvy/N)F(u,v)f(x,y)ex0y033j2ux/Mj2vy/N=ef(x,y)ex0y0=PfQj2ux/Mej2vy/Nex,y,u,v0,1,2,3M,N4j2/NeW0W0W0W0W1W2W0W2W4W0W3W6W0W3W6W9利用w的周期性，得：W2=-W0,W4=W0,W6=-W0,和w的对称性，得：w3=—w1,w2=—w0,则有：F(u,.v)F(u,.v)W0W0W0W0W0W1W0WW0W0W0WW0W1W0W111101P=Q=1111010211111j1j03041j1j1111050611111j1j07081j1j4j364j088j0080088j011jjPfQ=（写到上式就可以了）3688j88j8f=[0102;0304;0506;0708]p=[1111;1-j-1j;1-11-1;1j-1-j]第八章11）膨胀的结果为半径为5*r/4的圆（图略）2）膨胀的结果为边长为为3*r/2的正方形，在直角处用半径为r的内切圆弧连接（图略）3）为简化问题,设图像为等边三角形,则膨胀结果是边长为3*r/2和等边三角形，在夹角处以半径为r/4的圆弧连接(图略)4)腐蚀的结果为半径为3*r/4的圆(图略)腐蚀的结果为边长为为r/2的正方形(图略)为简化问题，设图像为等边三角形,则腐蚀结果是边长为0.134r的等边三角形(图略)第九章1、x=imread('baboon.bmp');[mnk]=size(x);y=uint8(zeros(m,n,k));y(:,:,1)=x(:,:,1);y(:,:,2)=x(:,:,3);y(:,:,3)=x(:,:,2);figure,imshow(x),title('原图像');figure,imshow(y),title('交换绿蓝通道');z=x;z(:,:,3)=2*x(:,:,3)fori=1:1:mforj=1:1:nif(z(i,j,3)>255)z(i,j,3)=255endendendfigure,imshow(z),title('蓝色通道加倍');2、依题意,R=200,G=50,B=150AA"I第十一章：1.现有8个待编码符号M0,M1,……,M7,它们的概率分别为0.40,0.25,0.11,0.09,0.06,0.04,0.03,0.01，试求这一组符号的信号熵，利用哈夫曼编码求出这一组符号的编码，画出哈夫曼树，并计算平均码长和编码效率M4:1M7:01M0:0001M5:00001M6:00000M2:0010M1:00111M3:00110平均码长：R=1*0.39+2*0.25+4*0.11+5*0.05+5*0.06+4*0.08+5*0.02+5*0.04=2.5图像的熵H为：nHPkgPk2.43k1编码效率：n=H/R=2.43/2.5=97.2%2.二分法香农-范诺编码方法。其步骤如下：1）首先统计出每个符号出现的概率；2）从左到右对上述概率从大到小排序；3）从这个概率集合中的某个位置将其分为两个子集合，并尽量使两个子集合的概率和近似相等，给前面一个子集合赋值为0，后面一个子集合赋值为1；4）重复步骤3,直到各个子集合中只有一个元素为止；5）将每个元素所属的子集合的值依次串起来，即可得到各个元素的香农-范诺编码。码字符号出现概率0M40.390.39(0)第PAGE\*MERGEFORMAT#页（共13页）第PAGE\*MERGEFORMAT#页（共13页）100M70.250.61(1)0.36(0)0.25(0)101M00.110.11(1)1000M20.080.25(0)0.14(0)0.08(0)1001M60.060.06(1)1010M50.050.11(1)0.05(0)10110M30.040.06(1)0.04(0)10111M10.020.02(1)平均码长R=2.53图像熵H=2.43编码效率n=2.43/2.53=96.04%参考答案二:符号出现概率-log2Pi-log2Pi+1码字长度累加概率转为二进制M4:0.391.362.3620r00M70.252.003.0020.3901M00.113.184.1840.641010M2:0.083.644.6440.751100M60.064.065.0650.8311010M50.054.325.3250.8911100M3:0.044.645.6450.9411110M10.025.646.6460.98111110平均码长R2.91熵2.43编码效率83.50515