《数字图像处理第三章答案.doc》由会员分享,可在线阅读,更多相关《数字图像处理第三章答案.doc(20页珍藏版)》请在三一办公上搜索。
1、3.1 a为正常数的指数式对于构造灰度平滑变换函数是非常有用的。由这个基本函数开始,构造具有下图形状的变换函数。所示的常数是输入参数,并且提出的变换必须包含这些参数的特定形式(为了答案曲线中的不是所要求的参数)。解:由(a)图所示,设,则在r=0时,T(r)=A在r=时,T(r)=A/2联立,解得则由(b)图所示,可以由(a)图翻转得到,所以(b)图的表达式s=(c)图是(b)图沿y轴平移得到,所以(c)图的表达式3.19 (a)在3.6.2节中谈到,分布在图像背景上的孤立的亮和暗的像素团块,当它们小于中值滤波器区域的一半时,经过中值滤波器处理后会被滤除(被其邻值同化)。假定滤波器尺寸为,n为
2、奇数,解释这种现象的原因?答:在的滤波器中有个像素,n为奇数,中值为,则有个像素小于或者等于,其它的大于或等于。当其中孤立的亮或者暗的像素A在像素团块中小于中值滤波器的一半时,即使在当所有群集点包含过滤屏蔽的极端情况下,没有足够的在其中任何一个集群点等于中值。如果在区域的中心点是一个群集点,它将被设置为中位数值,而背景的阴影将“淘汰”出集群。这一结论适用于当集群区域包含积分少集群的最大规模的较极端情况下。(b)考虑一副有不同像素团块的图像,假设在一个团块的所有点都比背景凉或者暗(但不是同时既比背景亮又比背景暗),并且每个团块的尺寸不大于。试求当n符合什么条件时,有一个或多个这样的团块像(a)中
3、所说的那样被分离出来?答:在A的结论下,我们考虑的团块的像素个数不可能超过,两个相近的或亮或暗的团块不可能同时出现在相邻的位置。在这个的网格里,两个团块的最小距离至少大于,也就是说至少在对角线的区域分开跨越(n-1)个像素在对角线上。3.29 CCD电视摄像机用于每天24小时,每月30天对同一区域进行长期观测研究。5分钟拍一次数字图像并传送到中心场所。场景的照明,白天为自然光,晚上为人造光,没有无照明的时间,因此摄像机本身并不需要使用任何补偿装置。另外,使用数字技术对图像进行后处理并归一化,这样就使图像与恒定照明是等效的。对此,设计一种方法。可以在实验室内使用希望的任何方法,但要在设计中明确列
4、出所做的所有假设。答:本题是考虑到范围的照明停留在线性部分的相机的反应范围,但是没有值范围内给出。图像待在线性的范围。唯一的方式,建立一个基准值的照明就是当变量(日光照明已不存在)。让F0(X,Y)表示图像只有采取人工光照条件下,没有移动的物体(如人或车辆)在现场。这成为标准的所有其他图像将正常化。选择各代表子区域(x,y)不可能被遮蔽移动物体,计算他们的平均强度。然后,选择所有个人的平均值,fmin表示最低目标,fmax最高目标,处理任何输入图像,f(x,y),因此,其最低和最大的分别等于fmin和fmax。线性变换函数的形式:其中f是输出图像。这是很容易验证,输出的图像将有(如果选择所需的
5、最低和最高值):在fmax和fmin是输入图像的最高值与最低值。假设所有的图像的线性范围之内。相机的工作范围内,从而饱和等非线性不是一个问题。3.3提出一组能够产生8比特单色图像所有独立位平面的灰度分成变换(例如,变换函数T(r)=255,当r在0,127范围内时,T(r)=0,而当r在128,255范围内,T(r)=255,此时的函数可以产生一幅8比特图像的第7位平面图像)解: 3.10一幅图像的灰度PDF, 示于下图。现在对比此图像进行灰度变换,使其灰度表达式为下面右图的。假设灰度值连续,求完成这一操作的变换(r到z)。解:由左图可知由右图得到:即:由图可知:故: 3.20(a)提出一种过
6、程来求一个l领域中值?(b)试提出一种技术,逐像素地移动邻域的中心来更新中值。解:(a)设的中值为m,其中最大值设为a则(b)一旦值已经被分类一次,我们仅仅是删除在缓慢移动向领域的值,插入首要领域的值到分类排列的最恰当的位置。2.18在下一章中我们将讨论算子,其函数在一个很小的子图像区S计算像素总数。说明这些都是线性算子。答:让H表示领域的求和运算符,让f和g表示两个不同的小子图像领域,让f+g表示f图像和g图像里的相应像素值的总和,H是在给定一个领域里计算像素值总和的算子,将f和g分别乘以两个常量a、b,所以表示f图像的像素值得a倍加上g图像的像素值得b倍,所以我们可以推导: = = =正如
7、式(2.6.1)所示,所以这些计算图像区域像素总数的算子都为线性算子。Prob4: (a)通常,如果将低阶比特面设为零值,对一幅图像的直方图有何影响?答:如果将低阶比特面设为零值,该图像会丢失细节。即不同灰度值的像素个数将减少,这会导致直方图的成分数减少。由于像素个数不会改变,这将在总体上导致直方图峰值高度上升。通常,较低的灰度值变化将减少对比度。 (b)如果将高阶比特面设为零值,对直方图有何影响?答:如果将高阶比特面设为零值,该图像会丢失轮廓,即丢手视觉上的很多数据。最明显的影响是使图像非常模糊,根据灰度变换函数,将0127之间的所有灰度映射为0,下降的最高位将限制到127的8位图像中最亮的
8、水平。由于像素数将保持不变,一些直方图峰值的高度会增加。一般直方图的形状将更高更窄,过去127没有直方图组件。Prob21: (a)在识别的应用领域,文本页通过图3.2(b)所示的阈值变换函数简化为二值图像。这遵循如下过程,即细化字符直到它们成为全“0”背景上的一串“1”。由于有噪声存在,故二值化和细化处理时,导致在连1处有缝隙存在,缝隙有13个像素宽。修复缝隙的一种方法是,对二值图像使用均值掩膜来模糊它,这样会在缝隙间桥连非零像素。试求出能执行该任务的均值掩膜的最小尺寸?答:最极端的情况是,当面具被定位在沿着一条薄段离中心像素3像素的差距,在这种情况下,一个33掩模将包括一个完全空白的领域。
9、因为这是最大的差距,下一个(奇数)面具的大小应包括一些在薄段的像素。因此,能执行该任务的均值掩膜的最小尺寸为55。 (b)桥连缝隙后,为了转换回二值形式,要进行阈值处理。在(a)中得出的答案里,完成这一任务且不产生断线所要求的最小阈值是什么?答:当面具包含了只有两个像素的片段时将产生最小的平均值。该均值是一个灰色的刻度值,而不是二进制,像片段像素值其余的部分。用指代最小平均值,用指代薄段的像素的二进制值。很明显,比小,那么略高于的二值化阈值设置将在掩模中心创建一个二进制像素值B。3.7 在实际应用中,将输人图像的直方图模型化为高斯概率密度函数,其概率密度函数形 式为:其中m和分别是高斯PDF的
10、平均值与标准差。具体处理方法是将m和看做给定图像的平均灰度级和对比度,试求出直方图均衡化的变换函数。解:一般直方图均衡化的变换函数为:高斯密度函数一般-到+,实际不可能实现。一是假设标准偏差足够小,r对P(r)影响可忽略不计。;二是比例放大值直到区域到尾部可忽略不计,变换函数: 一般范围取0 255。3.22 以下的三幅图像是分别通过n=23,25和45的方形均值掩模处理后的模糊图像。图(a)和(c)中左下角的垂直竖条被模糊了,但竖条与竖条之间的分割仍然很清楚。但图(b)中的竖条却已经融人了整幅图像,尽管产生这幅图像的掩模要比处理图像(c)的小得多,请解释这一现象。解:从图可知,垂直线有5个像
11、素宽,100像素高,他们的间隔是20像素。问题是相关的现象与水平之间的间隔线有关,所以我们可以简化问题,考虑一个单一的扫描行通过线的图像。回答这个问题的关键在于实际之间的距离(无像素)开始的线条,下一个(其右面)是25个像素。考虑扫描线,如图,同样显示是一个断面25 x25掩膜。掩膜反应包括的像素是平均的。我们注意到,当一个像素掩膜移动右面,它失去了左边竖线的价值,可是它捡起一个相同的一个在右边,所以反应不会改变。事实上,多少像素属于垂直线和包含在掩膜并不会改变,无论在掩膜的任何地方(只要是包含在线内,而不是在边缘附近线)。这一事实的线像素数量低于掩膜并不会改变是由于特有的线条和分隔线之间的宽
12、度的相当于25像素。这个常数宽度的反应是没有看到白色的差距在问题的声明中图像显示的理由。注意这个常数不发生在23 x23或45 x45的掩膜,因为他们不是同步与线条宽度和将它们分开的距离。MATLAB实现:i = imread(prob3_22.tif);% subplot(2,2,1); % imshow(i); subplot(2,2,2);%imshow(i,30 200); w1 = fspecial(average,23 23);% J=imadd(w,100);% g1=imfilter(i,w1,replicate); imshow(g1); subplot(2,2,3);%im
13、hist(i) w2 = fspecial(average,25 25); g2=imfilter(i,w2,replicate); imshow(g2); subplot(2,2,4);%imhist(i) w2 = fspecial(average,45 45); g3=imfilter(i,w2,replicate); imshow(g3);效果如下:3.12有两幅图像f(x,y)和g(x,y),其直方图分别为和。给出能确定直方图的条件,并简述在每种情况下如何的得到直方图(a)(b)(c)(d)解:对图像进行重新规整,设作四则运算后图像的像素值为f,像素值的最大值和最小值,规整后的像素值
14、为f(a)(b)(c)(d)根据f的分布情况即可得出直方图3.23考虑图3.26所示的应用,即消除图像中比q*q像素大小的方形包围的物体小的目标。假设想要将目标的平均灰度值减少为原来平均灰度值的1/10.用这种方法,那些目标可以接近背景灰度并用门限法消除。给出平均掩模的最小尺寸(奇数),该掩模仅对整幅图像处理一次就将平均灰度级减少到所希望的程度.解:由掩模定义的邻域像素灰度的平均值为 设一个平均掩膜尺寸为,则对于大小的物体来说,经过滤波后的灰度,要等于原来的1/10,那么,边长约为物体的3倍长,这就是滤波器最小所需的尺寸。3.13 考虑两幅8比特图像,它们的灰度覆盖了整个0 ,255的范围。(
15、e) 讨论反复将图像(b)从图像(a)中减去的最终效果。(f) 如果把图像的顺序颠倒会得到不同的结果吗?解:(a)做反复减法的过程可用下式表示: 最大数和最小数实行减法的次数以及灰度范围取值可分为两种情况:(1) 设整个图像中至少有一个像素值为255,所以,进行减法不超过-511的最大取值K为3,此时,对坐标(s,t)有,a(s,t)=b(s,t)=255.在这种情况下,差分图像中的值得可能范围是-511到255。取255的情况为,对坐标(s,t)有,a(s,t)=255,b(s,t)=0。(2) K能取的最小值为2,对坐标(s,t)有,a(s,t)=0,b(s,t)=255。在这种情况下,差
16、分图像中的值得可能范围是-510到255。取255的情况为,对坐标(s,t)有,a(s,t)=255,b(s,t)=0。(b)做反复减法的过程可用下式表示: 最大数和最小数实行减法的次数以及灰度范围取值可分为两种情况:(1)设整个图像中至少有一个像素值为255,所以,进行减法不超过-511的最大取值K为3,此时,对坐标(s,t)有,b(s,t)=a(s,t)=255.在这种情况下,差分图像中的值得可能范围是-511到255。取255的情况为,对坐标(s,t)有,b(s,t)=255,a(s,t)=0。(2)K能取的最小值为2,对坐标(s,t)有,b(s,t)=0,a(s,t)=255。在这种情
17、况下,差分图像中的值得可能范围是-510到255。取255的情况为,对坐标(s,t)有,b(s,t)=255,a(s,t)=0。3.24 在给定的应用中,一个均值掩模被用于输入图像以减少噪声,然后再用一个拉普拉斯掩模来增强图像中的小细节,如果将这两个步骤交换一下,结果是否会相同?解:交换步骤结果相同。1)因为,所以,2)由1)中知道先用拉普拉斯掩模结果为 再用均值掩模为:因为以及等的均值都为零,所以进行均值掩模后原式15证明式3.4.4和3.4.5的正确性。证明:一副将噪声加入到原始图像f(x,y)中所形成的带有噪声的图像g(x,y),即g(x,y)=f(x,y)=(x,y), 我们通过累加一
18、组噪声图像来减少噪声, 即有g(x,y)= 那么Eg(x,y)=f(x,y) 其中,在所有坐标点(x,y)上,Eg(x,y)是g的期望值,与分别是g与n的方差,在平均图像中,任何一点的标准差为: 故原式得以证明。25 证明拉普拉斯变换是各向同性的(旋转不变),需要下列轴旋转角的坐标方程: X=x -y Y=x +y 其中(x,y)为非旋转坐标,而(x,y)为旋转坐标。解:拉普拉斯变换:f= 旋转变换:f= 而已知:X=x -y Y=x +y 代入计算得, 而 可见拉普拉斯变换是旋转不变的。3.16在工业应用中,X光摄影法用来检查组合铸件的内部结构。其目的是发现铸件内部是否有缝隙.这些缝隙在图像
19、中一般表现为小气泡。但是,由于铸件材料的性质以及我们使用的X线能量,较高的噪声常会导致检查困难,所以我们决定使用图像均值处理来减小噪声和改进视觉对比度。在计算均值的过程中,为减少时间,在成像中保留固定分块以尽可能保持小的图像数量非常重要。经过数次实践,发现块数为10时噪声的方差明显地减小了。如果成像装置每秒产生30幅图像,那么在成像时保持固定块数而想达到预先要求的铸件噪声方差的减小量需要多少时间?假设图像中的噪声为非相关的且其均值为O。解:当i= 1(不平均),则有 g(1)=g1 且2 g(1)= 2 当i=k时, 有 g(k)=1/k*gi 且2g(k)=1/k*2 依题意可知要使:2 g
20、(1)=10*2g(k) 即有 k=10成像装置每秒产生30幅图像,所以在成像时保持固定块数而又想达到预期的铸件噪声方差的减小量需要1/3 秒。 3.26试求一个3*3的反锐化掩模来对一幅图像进行一次掩模处理。解:由式3.7.11可得如下两种掩模:Prob17.线性空间滤波处理要求在整幅图像中移动掩模的中心点,在每个处理区域中,计算掩模系数与该区域相应像素值乘积的总和。在低通滤波器中,所有系数和为1,我们使用所谓的盒式滤波法或移动均值算法(这种方法一次只更新掩模计算的一部分,并从一部分到另一部分顺次进行)。(g) 对一个nn的滤波器公式化这样一个算法,说明涉及的计算规律,以及围绕图像移动掩模时
21、所用的扫描序列。(h) 使用大规模处理的运算量与盒式滤波算法执行的运算量的比值称为“计算优势”。在本例中求出计算优势并作为n的函数,其中n1。由于这两种处理方法的标定系数都是1/,因此在获得计算优势时根本不用考虑他的影响。假设图像有为0的外部轮廓,这个轮廓足够厚,完全可以忽略它在处理过程中的影响。解:(a)如图所示33的掩模:R=w(-1,-1)f(x-1,y-1)+w(-1,0)f(x-1,y)+.+w(0,0)f(x,y)+.+w(1,0)f(x+1,y)+w(1,1)f(x+1,y+1)以33的掩模可推得:由于所有系数为1,则该低通滤波器掩模响应是将模板下的每个像素的灰度值相加所得。当模
22、板向右移动一个像素的时候获得新的一列计算得:其中,C1是模板移动之前所对应的列像素总和,C3是移动之后的列像素总和。C3是由2次相加所得,然后一次相加和相减得到,移动一次进行4次运算。扫描图像需要进行递归,即当一排的所有像素扫描完成后将向下移动一个像素然后又向相反的方向扫描。同理:对于nn的模板,C3需要由(n-1)次相加所得,那么每次移动需要(n+1)次运算得到。所以每次移动需要(-1)次加法运算。(b)令“计算优势”为A则A为n的线性函数,且n=2时,A=1。Prob27.使用式(3.7.4)给出的拉普拉斯变换的定义,证明将一幅图像减去其相应的拉普拉斯图像等同于对图像进行反锐化掩模处理。解
23、:式3.7.4为所以 其中表示(x,y)点所在像素和其上下左右4个像素的平均值。 与反锐化掩模处理表示相同,即证得将一幅图像减去相应的拉普拉斯图像等同于对图像进行反锐化掩模处理。3.18 讨论用一个3*3低通空间滤波器反复对一幅数字图像处理的结果,可以不考虑边界的影响。看一个空间滤波器的反复应用的最简单的方法之一是使用superProblem位置,让f(x,y)和h(x,y)分别表示该图象和滤波功能。假设正方形图像大小为NxN,为方便起见,我们可以表达F(X,Y)为在大多数幅图像的总和,其中每个只有一个非零像素(最初我们假设N可以是无限的)。然后,可以运行H(X,Y)的过程中,超过F(X,Y)
24、下面的卷积表示为:假设用于说明目的,连接(X,Y)在其中心具有价值1,而其他像素值0,如上所述(见图。P3.18a)。如果h(X,Y)是3掩模1/9zs(图P3.18b),然后卷积H(X,Y)连接(X,Y),将产生一个图像3x3,在其中心和0zs别处1/9zs阵列,如图P3.18(c)所示。如果h(X,Y)现在这个形象,将产生的图像如图P3.18(D)所示,需要注意的是两个非零像素总和为零。P3.18(c)及(d)是相同的,等于原像素值。因此,它是直观明显,连续应用的H(X,Y)会扩散fi(X,Y)的非零值(这一点也不意外因为H(X,Y)是一个模糊的过滤器)。因为结果保持不变,非零元素的值会变
25、得越来越小,作为申请过滤器的数目增加。总的结果是加入卷积 因此连续应用低通空间的净效应滤波器h(X,Y),其结果越来越模糊,每个像素值之间重新分配。因此模糊图像的平均价值将与F(X,Y)的平均值相同。据悉,每一个迭代模糊进一步从起点向外扩散。在极限情况,其值将得到无限小,但因为平均值保持不变,这将需要一个无限的空间比例的图象。这个交界处就是重要的边界条件,虽然它不是必需的问题的陈述,它是指导课堂上讨论H(X,Y)连续应用的影响有限比例的图像。净效应定义是图像的边界值无法向外扩散,分母平均值最终击败像素值的连续应用,使得图像极限为零。一个简单的例子,这是图P3.18(e),显示了一个大小为1x7
26、的阵列,连续应用1x3掩模,我们可以看到,只要模糊值可以化解,总和S所产生的像素是1。然而边界满足时,则必须在如何处理边境上的掩模操作。在这里,我们使用常用假设,即过去的边界像素值是0,但是掩模操作不超越边界。在这个例子中,我们看到连续应用时像素值的总和开始下降,。在极限情况下,周期压倒的像素值的总和,产生的阵列。3.28(a)证明式(3.7.13)给出的梯度值是一种各项同性的操作过程(见习题3.25.(b)证明如果梯度用式(3.7.14)进行计算,将使其失去各向同性的性质。(a) 从题3.25可知: 并且: 从中可知: 或者 所以,我们可以的吃(3.7.13)给出的梯度值是一种各项同性的操作过程。(b)从式(3.7.12)(3.7.14)以及前面推导的结果可得: 并且: 很明显 所以将失去各向同性的性质。