山东一医大医学数字图像处理习题集.docx

《山东一医大医学数字图像处理习题集.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(u,v).5.5 设1台X射线成象设备所产生的模糊可模型化为1个卷积过程，且卷积函数为循环对称的，即11(!卜+2。2)/。2刖-2_2耳其中於.2+丫2。为恢复这类图象要设计1个有约束最小平方米恢复滤波器，请推导它的转移函数。5.6 左下角为原点，求表