《山东一医大医学数字图像处理习题集.docx》由会员分享，可在线阅读，更多相关《山东一医大医学数字图像处理习题集.docx（11页珍藏版）》请在课桌文档上搜索。
1、医学数字图像处理习题集第1章习题1.1 连续图象f(x,y)与数字图象I(r,c)中各量的含义是什么？它们有什么联系和区别？它们的取值各在什么范围？1.2 图象处理、图象分析和图象理解各有什么特点？它们之间有哪些联系和区别？1.3 图象工程主要与哪些学科相关？它们之间如何互相影响？1.4 近年来有哪些新的数学工具在图象工程中得到应用？1.5 简要叙述图象采集、图象显示、图象存储以及图像通信与图象处理和分析的联系。1.6 近年来在图象采集、图象显示、图象存储各方面最突出的进展是什么？对图象处理和分析产生了哪些影响？1.7 本书主要可分成哪些层次？互相之间有什么关系？第2章习题2.1 人类视觉中最
2、基本的几个要素是什么？2.2 马赫带效应和同时对比度反映了什么共同问题？2.3 试说明舌形图里连接红、绿、蓝构成和三角形与HIS颜色三角形有何不同？2.4 证明计算S的公式对所有HlS颜色三角形里的点都成立。1.1 5为什么从RGB空间向HIS空间转换时分2段计算H,而从HIS空间向RGB空间转换时分3段计算H?2.6 发光强度及亮度与照度各有什么不同？2.7 空间点(1,2,3)经人=0.5的镜头透视后的摄象机坐标和图象平面坐标各是什么？2.8 波特率(baudrate)是一种常用的离散数据传输量度。当采用二进制时，它等于每秒所传输的比特数。现设每次先传输1个起始比特，再传输8个比特的信息，
3、最后传输1个终止比特，试计算(1)以9600波特传输1幅102425610=655360,则以9600波特的速率传输需要时间为655360/9600=68.27s。(2)传输的比特数为10241024X30=31457280,则以384000波特的速率传输需要时间为31457280/384000=81.92s。2.9 设有2个图象子集S和T,如图题2.9所示。rir0000000】I01 O0I00IO0Ioooooo00IrIOOOOO001IIOOlII图收2.9(1)如果v=l,试指出它们是否4-连通，-8连通，in-连通；(2)子集S和子集T以外的所有象素看成另一个子集U,试指出子集S
4、和子集T是否子集U4-毗邻，8-毗邻，m-毗邻。(3)如果将子集S和子集T以外的所有象素看成另一个子集U,试指出子集S和子集T是否与子集U-4毗邻，8-毗邻，m-毗邻。2 .10试画出(1)将单象素宽的8-通路转换为4-通路的程序流程；(2)将单象素宽的m-通路转换为4-通路的程序流程。2.11讨论画出由矩阵B代表的各元素间所有可能的空间排列。abcl0laB=011)Ll1U42.12证明2.7.3小节中计算B+的算法给出与式(2.7.3)完全相同的结果。2.13图题2.13给出了一个图象子集。312l()2202)211(A1O2困收2.13(1)令V=0,1,计算P和q之间的D4,D8和
5、Dm距离；(2)令V=1,2,仍计算上述3个距离。2. 14(1)证明P和q之间的D4,距离等于它们之间最短的4-通路的长度；(2)上述通路是惟一的吗？2.15如何仅利用逻辑运算从所给图象(图题2.15)中检测出二值圆环？2. 16(1)给出将图象顺时针旋转45的变换矩阵。(2)如何利用上述矩形阵实现图象旋转？(3)利用(1)中得到的矩阵转图象点(x,y)二(1,0)o2.17设给定如下平移变换矩阵T和尺度变换矩阵S,分别计算对空间点(1,2,3)先平移变换后尺度变换和先尺度变换后平移所得到的结果，并进行比较讨论。POO2 Io 1 0 4 （0016U） 0 0 1.40 0 01 0 3
6、0 0! 0 0 2 0 JO 0 0 1.2.18已知空间1个点成象在图象平面(a,b)处，现要将其移到(c,d)处，试分别写出用以下变换进行移动所需的变换矩阵：(1)只用平移变换；(2)只用尺度变换；(3)只用旋转变换。第3章习题3. 12-D傅里叶变换的分离性有什么实际意义？3. 2证明式(3.2.18)和式(3.2.19)成立。3. 3试举例说明2-D傅里叶变换的周期性和共枕对称性的用途。3.4 设x,y都是连续变量，计算下列各式的傅里叶变换。3.5 证明f(X)的自相关函数的傅里叶变换就是f(X)的功率谱|尸()。3.6 证明离散傅里叶变换和反变换都是周期函数(为简例可以用1-D函数
7、为例)。3.7 试讨论连续卷积和离散卷积的不同。3.8 证明2个函数卷积的傅里叶变换是这2个函数傅里叶变换的乘积(为简便可以用单变量函数为例)。3.9 1个实函数f(X)可分解成为1个奇函数fdd(X)和1个偶函数feven(x)之和。证明fevenx)=f(-x)/2Jodd(X)=f(x)-f(-x)/2;(2)证明Ffeven(x)=ReFf(x),Ffodd(x)ylmF(x).3. 10讨论证明；W滋=Wbl(2) Wf=Wh(3) W肪M=-W3.11在3.3.2小节中指出,为计算N点的FFT需要NIogzN次加法和(1/2)NlOg2N次乘法。如计算1幅NXN图的2-DFFT需要
8、多少次加法和乘法？3.12根据位对换规则对输和数据进行排序以计算1个16点的快速傅里叶变换。3.13(1)证明式(3.4.15)和式(3.4.17)组成1个变换对；(2)以24为例，验证式(3.4.23)和式(3.4.24)；(3)试根据式(3.4.34)和式(3.4.35),写出NE阶的哈达玛矩阵。3.14证明式(3.4.34)代入式(3.4.46)组成1个变换对。3.15写出N=2时的2-DDCT的正反变换形核的值；讨论反变换核的可分离性和对称性。3.16不考虑a(u),算出对应u=v=l块中的各值。3.17对N=8,计算哈尔变换矩阵。3.18对N=8,计算斯拉特变换矩阵。3. 19设有1
9、组随机矢量x=xlx2x3,其中xl=0Ol,x2=010,x3=lO0,请分别给出X的协方差矩阵和经霍特林变换所得到的矢量y的协方差矩阵。3. 20设有1组64X64的图象，算出它们的协方差矩阵是单位矩阵。如果只使用一半的原始特征值计算重建图象，那么原始图和重建图之间的均方误差是多少？第4章习题4. 1用灰度剪切变换函数T(r)=0,rE0,127;T(r)=255,rE128,255能从1幅8位面图象中提取出第7位面。据此给出1组变换函数以提取该图象的其他各个位面。4. 2为什么一般情况下对离散图象的直方图均衡化并不能产生完全平坦的直方图？4 .3设已用直方图均衡化技术对1幅数字图象进行了
10、增强，试证明再用这个方法对所得结果增强并不会改变其结果。5 .4设1幅图象有如图题4.4(a)所示直方图，拟对其进行规定直方图变换，所需规定直方图如图题4.4(b)所示。请参照上册表4.2.2列表给出直方图规定化计算结果(不必统计直方图各灰度级象素），并比较SML方法和GML方法的误差情况。4. 5证明式（4.2.12）和式（4.2.13）成立。4.6 设工业检测中工件的图象受到零均值不相关噪声的影响，如果图象采集装置每秒可采集30幅图，要采用图象平均方法将噪声的方差减少到1/10,那么工件需保持多长时间固定在采集装置前？4.7 空间滤波器在具体实现时需要让模板的中心移过图象中的每个位置，计算
11、模板系数和对应象素的乘积并对它们求和。对所有模板系数均为1的低通滤波器，可使用称为盒滤波器或移动平均的算法程序，即每次只计算在模板移动中其值变化的部分。4.8 编1个程序实现nXn中值滤波器，当模板中心移动图象中每个位置时，设计一种简便地更新中值的方法。4. 9利用33滤波器，证明在熔域中利用Highpass=Original-Lowpass可得到高通滤波的图象。4.10 画出高频增强滤波器的转移函数和脉冲响应曲线。4.11 设仅利用象素点（x,y）的4-近邻象素（不用点（x,y）本身）组成1个低通滤波器。（1）给出它在频域的等价滤波器H（u,v）;（2）证明所得结果确实是1个低通滤波器。4.
12、12 有一种计算梯度的基本步骤是计算f（x,y）和f（x+l,y）的差。（1）给出在频域进行等价计算所用的滤波器转移函数H（u,v）;（2）证明这个运算相当于1个高通滤波器的功能。4.13 讨论用于空间滤波的平滑滤波器和锐化滤波器的相同点，不同点以及联系。4. 14原什么条件下巴特沃斯低通滤波器变成理想低通滤波器？4.15 证明可以通过在频域内用原始图减去低通滤波图得到高通滤波的结果。4.16 从巴特沃斯低通滤波器出发推导它对应的高通滤波器。4.17 有1种常用的图象增强技术是将高频增强和直方图均衡化结合起来以达到使边缘锐化的反差增强效果，以上2个操作的先后次序对增强效果有影响吗？为什么？4.
13、18 在天体研究所获图象中有一些相距很远的对应恒星的亮点。由于大气散射原因而迭加的照度常使得这些亮点很难看清楚。如果对这类图象模型化为恒定亮度的背景与1组脉冲的乘积，根据同态滤波的概念设计1种增强方法将对应恒星的亮点提取出来。4.19试证明：如果H(u,v)是实对称的，则h(x,y)一定也得实对称的。4.20在1条自动装配线上，有3类形状相同的工件。为了方便检测，将工件用不同颜色标注。现只有1个单色摄影机，请提出1种用这个摄影机检测3种颜色的方法。4.21设有1个能输出RGB模拟信号的彩色摄影机，1个能将这些模拟信号转化为以(l30s)的视频速度输出RGB或HIS图明的数字化器，3块能以视频速
14、度接受图象的帧缓存卡，以及1个以视频速度计算直方图的硬件。所有这些都可与I台微机组合在一起。现要解决如下问题：生产线上有一系列形状相同但颜色不同的工件，它们按红、黄、绿、蓝的次序排列。请借助以上硬件设计1个图象处理软件系统将工件的颜色检测出来。这里假设工件移动速度相当慢，所以可忽略由此产生的图象模糊问题。请画出系统的流程图并对每个模块及所选处理技术进行讨论。第5章习题5.1 设有1个线性位移不变的图象退化系统，它的脉冲响应可以写成h(x-a,y-b)=exp-x-a2+(y-b)20如果给系统输入1个位于x=a处的无穷长细直线信号，它可用f(x,y)=6(-a)模型化，求系统的输出。5.2 设
15、1幅图象的模糊是由物体在X方向的匀加速运动产生的。当t=0时物体静止，在t=0到t=T间物体加速度是Xo(t)=at22,求转移函数H(u,v),并讨论匀速运动和匀加速运动所造成的模糊的不同特点。5.3 成象时由于长时间曝光受到大气干扰而产生图象模糊可以用转移函数H(u,v)=exp-(u2+v2)202表示，设噪声可忽略，求恢复这类模糊的维纳滤波器的方程。5.4 设恢复滤波器R(u,v)满足(#)=H(n,u)g(,并假设强制恢复图象的功率谱?(#)等于原始图象的功率谱倒(1)根据I产(外训2,H(%咪和N(比,训2求出R(u)o(2)用(1)中所得的结果以类似于式(5.4.8)的形式写出F
16、(u,v).5.5 设1台X射线成象设备所产生的模糊可模型化为1个卷积过程，且卷积函数为循环对称的，即11(!卜+2。2)/。2刖-2_2耳其中於.2+丫2。为恢复这类图象要设计1个有约束最小平方米恢复滤波器，请推导它的转移函数。5.6 左下角为原点，求表示几何失真过程的一对双线性等式和校正公式。如设f(l,1)=1,f(7,1)=7,f(l,7)=7,f(7,7)=14,求点f(2,4)的灰度值。5.7 证明(1)如果f(x,y)是旋转对称的，那么它可以由单个投影重建；(2)如果f(x,y)是可以分解成g(x)和h(y)的乘积，那么它可以由2个坐标轴垂直的投影重建。5.8 试证明傅里叶变换投
17、影定理。5.9 设在点(x,y)处的线性衰减系数为d(x,y),根据X射线衰减规律，有PIrLjay)=/(7)其中Ir为射线入射强度，I为射线透过物体后的强度，(S,o)代表从发射源到接受器的直线与物体相交的部分。现定义CT值为c=kd(K)-4,*，)d，x,y)其中dw(x,y)为水的线性衰性系数，k为归一化系数。当使用卷积逆投影重建法明，已知重建滤波器中激响应为h(s),求计算CT值的表达公式。第6章习题6.1 当对数底数为e时，得到的信息单位称为奈特(nat)；当对数底数为10时，得到的信息单位称为哈特利(hareley)。推导它们与比特的换算关系。6.2 对1个具有q个符号的零记忆
18、信源，证明它的端的最大值为Iogq,这个值当且仅当所有源符号出现概率相同时达到。提示：考虑logqf(u),并利用不等式InXWX6.3客观保真度准则和主观保真度准则各有什么特点？6.4 除书中介绍的保真度准则外，还有什么方法可以描述解码图象相对于原始图象的偏离程度？231106.5 设A=0,1,B=0,1,u=3414T,Q=Liz39zi0J,计算与信道有关的各个概率，包括P(a=0),P(a=1),P(b=0),P(b=1),P(b=Oa=O),P(b=Oa=1),P(b=la=o),P(b=la=l),P(a=Ob=O),P(a=Ob=l),P(a=lb=O),P(a=l,b=l),
19、P(a=0,b=0),P(a=0,b=l),P(a=l,b=0),P(a=l,b=l).6.6 考虑例6.2.2中的信源和二元对称信道，令pbs=34,pe=l3,试问(1)信源的蜡是多少？(2)当接受到输出时，关于输报告文学的不确定性减少了多少？(3)不确定性的改变与信道容量在数值上是什么关系？6.7 (1)请说明是否能用变长编友法压缩1幅已直方图均衡化的具有2n级灰度的图？(2)这样的图象中包含象素间冗余吗？6.8 (1)对1个具有3个符号的信源，有多少种惟一的哈夫曼码？(2)构造这些码。6.9 (1)计算在6.2.2中给出符号概率的信源的燃；(2)对信源符号构造哈夫曼码，解释这样构造的码
20、与6.2.2中第2种码的区别；(3)构造最优的Bl码；(4)构造最优的2bit二元平移码；(5)将所有符号分成2组，每组4个，然后构造最优的哈夫曼平移码；(6)对每个码计算平均字长，并将它们与(1)中算得的蜡进行比较。6.10 对6.3.3中所用信源的符号进行哈夫曼编码，给出码字、码字的平均长度和编码效率(并与算术编码进行比较)。6.11已知符号a,e,i,o,u,X的出现概率分别是0.2,0.3,0.1,0.2,0.1,0.1,对0.23355进行算术解码。6. 12（1）构造完整的4bit灰度码；（2）设计1个通用的方法将已灰度编码的数字转换成它的二值对应物，并用此法对0110101001
21、11进行解码。6. 13对1幅64X64二值图用每块有4个象素的I-DWBS法（见6.3.2小节）进行编码，已知对图象某一行的编码为OIlOOlOOOOOOlOOOOloOlOOOOO00,其中0代表黑色象素。（1）将这行码解出来；（2）设计1个1-D的迭代WBS方法，开始先搜索所有白线（64个象素的块），然后连续地将非白色区间二等分，直至达到4个象素的块；（3）用（2）中的方法对（1）中解出来的行编码，并与l-DWBS法比较所用的比特数。7. 14将给定图象（图题6.14）分解成3个位平面，然后用洲程编程编码方法逐行编码,给出码字，计算编码效率。8. 15（1）为什么在6.3.5中要将上1进
22、中对应e之后的第1个类似过渡点当作c，?（2）你能设计1个其他的方法吗？6 .16将6.4.1填充满。7 .17设1幅图象具有如（6.4.18）形式的自相关函数，其中6h=0,它已用1个二阶预测器按DPCM方式编码。（1）组成自相关矩阵R和矢量r；（2）找出最优的预测系数；（3）计算利用上述最优预测系数时的预测误差的方差。6. 18对LE和均匀概率密度函数1（2A）-A5AP（S）=1O其它推导Lloyd-Max的判断和重建值。第7章习题6.1 设1幅NXN二值图象中心有1个值为1的nXn的正方形区域，此处的象素值均为0。（1）根据式（4,3,5）,使用索贝宁算子计算这幅图的梯度，并画出梯度幅
23、度图（给出梯度幅度图中所有象素的值;(2)画出根据式(7,2,3)得到的梯度方向的直方图，并标出直方图每个峰的高度；(3)画出根据上册图7,2,4(a)的模板算出的拉普拉斯图，给出图中所有象素的值。7. 2如上册4.3.3小节中指出，从原始图中减去模糊图称为非锐化掩模。证明用上册图7.2.4(a)模板算得的拉普拉斯值正比于(只差1个1/4的系数)从原始图中减去4-邻域平均图所得的结果。7.1 (1)证明式(7.2.6)怕给算子的平均值为零；(2)证明任意图象与上述算子卷积后其平均值也为零；(3)说明当使用如图7.2.4中的模板似计算式(7.2.6)时，上述(2)中的结论是否仍能成立？7.2 给
24、出能使图7.2.7中直线子空间基的各个模板产生最强响应的直线的方向。里设直线的宽度为1个象素。1.1 5对某些2-D模板在图象中漫游1遍进行卷积的结果也可用相应的1-D模板在图象中分别漫游2遍并进行卷积得到。例如4.3.2小节中的3X3平滑滤波器的2-D模块对图象的处理效果可通过先用模板111在图象中漫游卷积1遍，并将这个1-D模板旋转90。后再漫游卷积前次结果而得到。现请证明图7.2.2中的索贝尔模板的效果可通过先用差分模板T01在图象中漫游卷积1遍，再用平滑模板121旋转90后再漫游卷积前次结果而得到。7.6 设1幅二值图象中含有水平、垂直、倾斜45和倾斜135。的各种直线，请设计1组可以
25、用来检测这些直线中单象素宽间断长度为多个象素的3X3模板。7.7 提出一种在梯度图象中检测直线段中的长度为1到L象素之间断的方法。这里可设背景灰度为常数，所有线段都是单象素宽且已标为同样的，与背景灰度值不同的灰度值。提示：不要试图构造检测断裂的模板，而要根据2.7.2小节中关于8-连通的概念进行。7.8 (1)解释为什么图7.2.13(a)中点1在图7.2.13(b)中对应的是1条直线？(2)点1是能产生如此结果的惟一一个点吗？(3)解释图7.2.13(1)中的反射相连关系。7.9 将图7.3.3中指出的所有边缘元素都标在图7.3.2(a)的图象上，并计算最小代价通路的代价。7.10 设有1幅
26、无噪声的NXN图象，其中左半边象素的灰度值为I,右半边象素的灰度值为J,且JI。现设有另1幅无噪声的NXN的图象，其灰度值从最左1列的0线性增加到最后1列的K,KJ,将2幅图相乘，得到1幅图象，春直方图是怎样的？7.11 11从式(7.4.9)出发推导式(7.4.10)。7.12 12从式(7.4.5)出发推导式(7.4.11)。7.13设1幅图象具有如图题7.13所示的灰度分布，其中Pl(z)对应目标，P2(z)对应背景。如果P1=P2,求分割目标和背景的最佳阈值。图即7.137.14一幅图象背景部分的均值为25、方差为625,在背景上分布着一些互不重叠的均值为150、方差为400的小目标。设所有目标合起来约占图象总面积的20%,提出1个基于取阈值的分别割算法将这些目标分割出来。7.15提出1种基于区域生长的方法解决上册题7.14中的问题。7.16用分裂合并法分割如图题7.16的图象，并给出对应分割结果的四叉树。田眄7.167.17参照表7.1.1将自己所知道的图象分割算法分别归到PB,PR,SB,SR这4类中去。7.18试分析为什么典型的EAGh泌(L)和EAGe(L)曲线都是单峰曲线。
