图像处理与分析-王伟强-作业题及标准答案汇总-.docx
- 文档编号:150632
- 上传时间:2023-04-28
- 格式:DOCX
- 页数:28
- 大小:1.53MB
图像处理与分析-王伟强-作业题及标准答案汇总-.docx
《图像处理与分析-王伟强-作业题及标准答案汇总-.docx》由会员分享,可在线阅读,更多相关《图像处理与分析-王伟强-作业题及标准答案汇总-.docx(28页珍藏版)》请在冰点文库上搜索。
【作业1】
1、完成课本习题3.2(a)(b),课本中文版《处理》第二版的113页。
可以通过matlab帮助你分析理解。
a:
s=Tr=11+(mr)E
b:
E控制函数的斜坡,也就是函数的倾斜程度,E越大,函数倾斜程度越大,如下图1,图2所示:
图1:
E=5
图2:
E=20
2、一幅8灰度级图像具有如下所示的直方图,求直方图均衡后的灰度级和对应概率,并画出均衡后的直方图的示意图。
(计算中采用向上取整方法,图中的8个不同灰度级对应的归一化直方图为[0.170.250.210.160.070.080.040.02])
【解答】直方图均衡采用公式
Sr=Gw=0rPrw-1
式中,G为灰度级数,取8,pr(w)为灰度级w的概率,Sr为变换后的灰度,计算过程如下表所示:
灰度级r
各级概率Pr(r)
累积概率w=0rPrw
累积概率×8-1
向上取整sr
0
0.17
0.17
0.36
1
1
0.25
0.42
2.36
3
2
0.21
0.63
4.04
5
3
0.16
0.79
5.32
6
4
0.07
0.86
5.88
6
5
0.08
0.94
6.52
7
6
0.04
0.98
6.84
7
7
0.02
1
7
7
则新灰度级的概率分别是:
Ps(0)=0
Ps
(1)=Pr(0)=0.17
Ps
(2)=0
Ps(3)=Pr
(1)=0.25
Ps(4)=0
Ps(5)=Pr
(2)=0.21
Ps(6)=Pr(3)+Pr(4)=0.23
Ps(7)=Pr(5)=Pr(6)=Pr(7)=0.14
编写matlab程序并绘制直方图:
s=0:
1:
7;
p=[00.1700.2500.210.230.14];
bar(s,p);
axis([-1800.3]);
可以看出,此图较题目原图更加“均匀”。
【作业2】
1、完成课本数字图像处理第二版114页,习题3.10。
【解答】
由图可知
prr=-2r+2,0≤r≤1
pzz=2z,0≤z≤1
将两图做直方图均衡变换
s1=T1r=0rprwdw=0r-2w+2dw=-r2+2r
s2=T2z=0zpzwdw=0z2wdw=z2
令上面两式相等,则
z2=-r2+2r
因为灰度级非负,所以
z=-r2+2r
2、请计算如下两个向量与矩阵的卷积计算结果。
(1)[123454321]*[20-2]
(2)-101-202-101*1320410323041052321431042
【解答】
(1)设向量a=[123454321],下标从-4到4,即a(-4)=1,a(-3)=2……a(4)=1;设向量b=[20-2],下标从-1到1,即b(-1)=2,b(0)=0,b
(1)=-2;设向量c=a*b,下标从-5到5。
根据卷积公式可知
cx=t=-∞∞atbx-t=t=-44atbx-t
其中,-5≤x≤5,则
c(-5)=a(-4)b(-1)=1*2=2
c(-4)=a(-4)b(0)+a(-3)b(-1)=1*0+2*2=4
c(-3)=a(-4)b
(1)+a(-3)b(0)+a(-2)b(-1)=1*(-2)+2*0+3*2=4
c(-2)=a(-3)b
(1)+a(-2)b(0)+a(-1)b(-1)=2*(-2)+3*0+4*2=4
c(-1)=a(-2)b
(1)+a(-1)b(0)+a(0)b(-1)=3*(-2)+4*0+5*2=4
c(0)=a(-1)b
(1)+a(0)b(0)+a
(1)b(-1)=4*(-2)+5*0+4*2=0
c
(1)=a(0)b
(1)+a
(1)b(0)+a
(2)b(-1)=5*(-2)+4*0+3*2=-4
c
(2)=a
(1)b
(1)+a
(2)b(0)+a(3)b(-1)=4*(-2)+3*0+2*2=-4
c(3)=a
(2)b
(1)+a(3)b(0)+a(4)b(-1)=3*(-2)+2*0+1*2=-4
c(4)=a(3)b
(1)+a(4)b(0)=2*(-2)+1*0=-4
c(5)=a(4)b
(1)=1*(-2)=-2
所以卷积结果为:
[244440-4-4-4-4-2]
(2)设矩阵
b=-101-202-101
下标从(-1,-1)到(1,1),即b(-1,-1)=-1,b(-1,0)=0……b(1,1)=1;
设矩阵
a=1320410323041052321431042
下标从(-2,-2)到(2,2),即a(-2,-2)=3,a(-2,-1)=2……a(2,2)=4;
设矩阵c=a*b=b*a,下标从(-3,-3)到(3,3)。
根据卷积公式可知
cx,y=s=-∞∞t=-∞∞as,tbx-s,y-t=s=-22t=-22as,tbx-s,y-t
其中,-3≤x≤3,-3≤y≤3,则
c(-3,-3)=a(-2,-2)b(-1,-1)=3*(-1)=-3
……
c(0,0)=a(-1,-1)b(1,1)+a(-1,0)b(1,0)+a(-1,1)b(1,-1)
+a(0,-1)b(0,1)+a(0,0)b(0,0)+a(0,1)b(0,-1)
+a(1,-1)b(-1,1)+a(1,0)b(-1,0)+a(1,1)b(-1,-1)
=3*1+4*2+0*1+2*0+1*0+3*0+1*(-1)+0*(-2)+2*(-1)
=8
……
c(3,3)=a(2,2)b(1,1)=4*1=4
所以卷积结果为:
-1-3-13-204
-3-6-44-4211
-3-7-63-6415
-3-11-48-10317
-7-1125-10615
-8-56-4-698
-3-13-3-242
【作业3】
1、高斯型低通滤波器在频域中的传递函数是
Hu,v=Ae-u2+v22σ2
根据二维傅里叶性质,证明空间域的相应滤波器形式为
hx,y=A2πσ2e-2π2σ2x2+y2
这些闭合形式只适用于连续变量情况。
在证明中假设已经知道如下结论:
函数e-πx2+y2的傅立叶变换为e-πu2+v2
【解答】
IDFTHu,v=-∞∞-∞∞Ae-u2+v22σ2ej2πux+vydudv
=-∞∞-∞∞Ae-u22σ2+j2πuxe-v22σ2+j2πvydudv
=-∞∞-∞∞Ae-12σ2u2-j4πσ2uxe-12σ2v2-j4πσ2vydudv
=-∞∞-∞∞Ae-12σ2u2-j4πσ2ux-4π2σ4x2+4π2σ4x2e-12σ2v2-j4πσ2vy-4π2σ4y2+4π2σ4y2dudv
=-∞∞-∞∞Ae-u-j2πσ2x22σ2e-2π2σ2x2e-v-j2πσ2y22σ2e-2π2σ2y2dudv
令r=u-j2πσ2x,s=v-j2πσ2y,则du=dr,dv=ds,上式写成:
=-∞∞-∞∞Ae-r22σ2e-2π2σ2x2e-s22σ2e-2π2σ2y2drds
=Ae-2π2σ2x2+y2-∞∞e-r22σ2dr-∞∞e-s22σ2ds
=Ae-2π2σ2x2+y22πσ2πσ12πσ-∞∞e-r22σ2dr12πσ-∞∞e-s22σ2ds
因为后两项是高斯分布,在-∞到∞积分为1,故上式等于:
=A2πσ2e-2π2σ2x2+y2=hx,y
命题得证。
2、第二版课本习题4.6(a)
【解答】
先来证明结论-1x+y=ejπx+y
根据欧拉公式展开等式右边
ejπx+y=cosx+yπ+jsinx+yπ
因为x,y均为整数,故sinx+yπ=0,
当x+y为奇数时,cosx+yπ=-1
当x+y为偶数时,cosx+yπ=1
故-1x+y=ejπx+y
再来证明题中等式
DFTfx,y-1x+y=DFTfx,yejπx+y
=1MNx=0M-1y=0N-1fx,yejπx+ye-j2πuxM+vyN
=1MNx=0M-1y=0N-1fx,ye-j2π-xM2M-yN2Ne-j2πuxM+vyN
=1MNx=0M-1y=0N-1fx,ye-j2πxu-M2M+yv-N2N
=Fu-M2,v-N2
3、观察如下所示图像。
右边的图像这样得到:
(a)用左侧图像乘以-1x+y;(b)计算离散傅里叶变换(DFT);(c)对变换取复共轭;(d)计算离散傅里叶反变换;(e)结果的实部再乘以-1x+y。
用数学方法解释为什么会产生右图的效果。
(忽略中间和右侧的黑白条纹,原题没有)
【解答】
已知
IDFTFu,v=1MNu=0M-1v=0N-1Fu,vej2πuxM+vyN=fx,y
则傅里叶变换的共轭复数进行傅里叶反变换的结果如下:
IDFTF*u,v=1MNu=0M-1v=0N-1Fu,ve-j2πuxM+vyN
=1MNu=0M-1v=0N-1Fu,vej2πu-xM+v-yN=f-x,-y
设原始图像为fx,y,经过(a)变换后得到
-1x+yfx,y
经过(b)变换后得到
Fu,v=1MNu=0M-1v=0N-1-1x+yfx,ye-j2πuxM+vyN
经过(c)变换后得到
F*u,v=1MNu=0M-1v=0N-1-1x+yfx,yej2πuxM+vyN
经过(d)变换后得到
IDFTF*u,v=1MNu=0M-1v=0N-11MNu=0M-1v=0N-1-1x+yfx,yej2πuxM+vyNej2πuxM+vyN
其实部为-1x+yf-x,-y,经过(e)变换后得到
-1x+y-1x+yf-x,-y=f-x,-y
最终效果是将原图像上下颠倒,左右颠倒,实现了旋转180度的效果。
【作业4】
1、请用公式列举并描述出你所知道的有关傅里叶变换的性质。
【解答】
1、时移性
2、频移性
3、均值
4、共轭对称性
5、对称性
6、周期性
7、线性
8、微分特性
9、卷积定理
10、相关定理
11、相似性
12、几种特殊函数的傅里叶变换
2、中文课本173页习题4.21
【解答】
没有区别。
补0延拓的目的是在DFT相邻隐藏周期之间建立一个“缓冲区”。
如果把左边的图像无限复制多次,以覆盖整个平面,那么将形成一个棋盘,棋盘中的每个方格都是本图片和黑色的扩展部分。
假如将右边的图片做同样的处理,所得结果也是一样的。
因此,无论哪种形式的延拓,都能达到相同的分离图像的效果。
3、
(1)假设我们有一个[0,1]上的均匀分布随机数发生器U(0,1),请基于它构造指数分布的随机数发生器,推导出随机数生成方程。
(2)若我们有一个标准正态分布的随机数发生器N(0,1),请推导出对数正态分布的随机数生成方程。
【解答】
(1)设U(0,1)可生成随机数w,用它来生成具有指数CDF的随机数z,其CDF具有下面的形式
Fz=1-e-az0,z≥0,z<0
令F(z)=w,当z≥0时,解方程1-e-az=w
得z=-1aln1-w
由于w∈[0,1],由上式可知z≥0,因而不存在z<0的情况。
所以随机数生成方程为:
z=-1aln1-U0,1
(2)设N(0,1)可生成随机数w,用它来生成具有对数正态分布CDF的随机数z,其CDF具有下面的形式
Fz=0z12πbe-lnv-a22b2dv
令F(z)=w,解方程
0z12πbe-lnv-a22b2dv=w
得z=ebw+a
所以随机数生成方程为:
z=ebN0,1+a
当b=1,a=0时,标准对数正态分布的随机数生成方程为:
z=eN0,1
4、对于公式
fx,y=s,t∈Sxygs,tQ+1s,t∈Sxygs,tQ
给出的逆谐波滤波回答下列问题:
(a)解释为什么当Q是正值时滤波对去除“胡椒”噪声有效?
(b)解释为什么当Q是负值时滤波对去除“盐”噪声有效?
【解答】
将原公式变形
fx,y=s,t∈Sxygs,tQ+1s,t∈Sxygs,tQ
=s,t∈Sxygs,tQgs,ts,t∈Sxygs,tQ
=s,t∈Sxygs,tQs,t∈Sxygs,tQgs,t
上式可看做求(x,y)邻域内所有(s,t)点的加权平均值,权重的分母s,t∈Sxygs,tQ是个常数,只需考虑分子gs,tQ的大小。
(a)当Q>0时,gs,tQ对gs,t有增强作用,由于“胡椒”噪声值较小,对加权平均结果的影响较小,所以滤波后噪声点处(x,y)取值和周围其他值更接近,有利于消除“胡椒”噪声。
(b)当Q<0时,gs,tQ对gs,t有削弱作用,由于“盐”噪声值较大,取倒数后较小,对加权平均结果的影响较小,所以滤波后噪声点处(x,y)取值和周围其他值更接近,有利于消除“盐”噪声。
【作业5】
1、请证明带通与带阻的频域关系公式,即课本中的关系公式
Hbpu,v=1-Hbru,v
【解答】证法一
设D0是带宽的径向中心,W是带宽,D是D(u,v)距滤波器中心的距离。
理想带通滤波器的频域公式为:
Hbpu,v=1,D0-W2≤D≤D0+W20,others
理想带阻滤波器的频域公式为:
Hbru,v=0,D0-W2≤D≤D0+W21,others
由此可知
Hbpu,v=1-Hbru,v
图形化表示即为:
图中黑色代表0,白色代表1,左图是带阻滤波器,右图是带通滤波器。
可以发现两图是“互补”的,若将对应点数值相加,则全为1。
证法二
一张图像f(x,y)可拆分为带阻部分fr(x,y)和带通部分fp(x,y),即
fx,y=frx,y+fpx,y
也可看做由带阻滤波器和带通滤波器分别卷积后叠加所得:
fx,y=fx,y*hbrx,y+fx,y*hbpx,y
等式两边取傅里叶变换得
Fu,v=Fu,vHbru,v+Fu,vHbpu,v
Fu,v=Fu,vHbru,v+Hbpu,v
Hbru,v+Hbpu,v=1
Hbpu,v=1-Hbru,v
2、复习理解课本中最佳陷波滤波器进行图像恢复的过程,请推导出w(x,y)最优解的计算过程,即从公式到
的推导过程。
【解答】
因为
σ2x,y=12a+12b+1s=-aat=-bbgx+s,y+t-wx,yηx+s,y+t-gx,y-wx,yηx,y2
所以对w(x,y)求偏导数,并令其等于0
∂σ2x,y∂wx,y=∂∂wx,y12a+12b+1s=-aat=-bbgx+s,y+t-gx,y+wx,yηx,y-ηx+s,y+t2=0
令a=ηx,y-ηx+s,y+t,b=gx+s,y+t-gx,y,N=2a+12b+1,则上式可简写成
∂σ2x,y∂wx,y=∂∂wx,y1Ns=-aat=-bbawx,y+b2
=1Ns=-aat=-bb∂awx,y+b2∂wx,y=0
已知
∂awx,y+b2∂wx,y=2awx,y+ba=0
的解为
wx,y=-ba
所以
wx,y=1Ns=-aat=-bbgx+s,y+t-gx,yηx+s,y+t-ηx,y
=1Ns=-aat=-bbgx+s,y+t-gx,yηx+s,y+t+ηx,yηx+s,y+t-ηx,yηx+s,y+t+ηx,y
令g=gx+s,y+t,η=ηx+s,y+t,则上式等于
=1Ns=-aat=-bbgη-gx,yη+gηx,y-gx,yηx,yη2x+s,y+t-η2x,y
=gx,yηx,y-gx,yηx,y+gx,yηx,y-gx,yηx,yη2x,y-η2x,y
=gx,yηx,y-gx,yηx,yη2x,y-η2x,y
推导成立。
3、考虑在x方向均匀加速导致的图像模糊问题。
如果图像在t=0静止,并用均匀加速x0(t)=at2/2加速,对于时间T,找出模糊函数H(u,v),可以假设快门开关时间忽略不计。
【解答】
由定义可知
Hu,v=0Te-j2πux0t+vy0tdt
将x0t=at22,y0t=0代入得
Hu,v=0Te-j2πuat22dt=0Te-jπuat2dt
上述积分难以化简,故不再推导。
4、已知一个退化系统的退化函数H(u,v),以及噪声的均值与方差,请描述如何利用约束最小二乘方算法计算出原图像的估计。
【解答】
频域中原图像的估计由下式给出
Fu,v=H*u,vHu,v2+γPu,v2Gu,v
其中P(u,v)是拉普拉斯算子的傅里叶变换。
定义“残差”向量r=g-Hf,由于Fu,v是γ的函数,则f和r都是γ的函数。
令ϕγ=rTr=r2,则它是γ的单调递增函数。
再调整γ使r2=η2±a,a是一个精确度因子。
已知噪声的均值为mu,方差ση2,和Hu,v、Pu,v,
(1)设定一个γ的初始值,
(2)计算r2,
(3)若满足r2=η2±a则执行第4步,若不满足,则调整γ大小,然后返回第2步。
(4)使用最新的γ,计算
Fu,v=H*u,vHu,v2+γPu,v2Gu,v
(5)再通过傅里叶反变换即可得到估计图像。
【作业6】
1、r,g,b是RGB彩色空间沿R,G,B轴的单位向量,定义向量
u=∂R∂xr+∂G∂xg+∂B∂xb
v=∂R∂yr+∂G∂yg+∂B∂yb
将gxx,gyy,gxy定义为这两个向量的点乘:
gxx=u⋅u=uTu=∂R∂x2+∂G∂x2+∂B∂x2
gyy=v⋅v=vTv=∂R∂y2+∂G∂y2+∂B∂y2
gxy=u⋅v=uTv=∂R∂x∂R∂y+∂G∂x∂G∂y+∂B∂x∂B∂y
推导出最大变换率方向θ和点(x,y)在θ方向上变化率的值F(θ)。
【解答】
要求最大变换率方向,即求使ucosθ+vsinθ2取最大值的θ。
ucosθ+vsinθ2=u2cos2θ+2uvsinθcosθ+v2sin2θ
=12gxx1+cos2θ+gxysin2θ+12gyy1-cos2θ
=12gxx+gyy+12gxx-gyycos2θ+gxysin2θ
将上式对θ求偏导得
∂∂θ12gxx+gyy+12gxx-gyycos2θ+gxysin2θ
=-gxx-gyysin2θ+2gxycos2θ
令上式等于0,解方程求出θ的值
-gxx-gyysin2θ+2gxycos2θ=0
gxx-gyysin2θ=2gxycos2θ
tan2θ=2gxygxx-gyy
θ=12tan-12gxygxx-gyy
此即是最大变换率方向,对应的变化率值为
Fθ=12gxx+gyy+12gxx-gyycos2θ+gxysin2θ
【作业7】
1、请根据课本中Z变换的定义,证明如下结论。
(1)若xn的Z变换为Xz,则-1nxn的Z变换为X-z。
(2)若xn的Z变换为Xz,则x-n的Z变换为X(1/z)。
【解答】
(1)由Z变换定义可知,xn的Z变换是:
Xz=-∞∞xnz-n
则-1nxn的Z变换为:
-∞∞-1nxnz-n=-∞∞-1-nxnz-n=-∞∞xn-z-n=X-z
(2)x-n的Z变换为:
-∞∞x-nz-n
令n'=-n,则上式可写成:
-∞∞xn'zn'=-∞∞xn'1z-n'=X1z
2、若G1z=-z-2K+1G0-z-1成立,请证明g1n=-1n+1g02K-1-n。
【解答】
已知g1n的Z变换是G1z,g0n的Z变换是G0z,
根据Z变换的时间翻转性x-n⟺Xz-1可得
g0-n⟺G0z-1
根据Z变换的平移性xn-k⟺z-kXz可得
g0-n-2K+1⟺z-2K+1G0z-1
根据Z变换的Z域翻转性-1nxn⟺X-z可得
-12K-1-ng02K-1-n⟺-z-2K+1G0-z-1
-12K-1-1-ng02K-1-n⟺-1-2K+1z-2K+1G0-1-1z-1
因为2K-1为奇数,所以-12K-1=-1,-1-n=-1n,-2K+1为奇数,故-1-2K+1=-1,上式继续推导得
-1n+1g02K-1-n⟺-z-2K+1G0-z-1=G1z
所以G1z的Z反变换是-1n+1g02K-1-n
又因为g1n也是G1z的Z反变换,所以
g1n=-1n+1g02K-1-n
3、假设课本中给出完美重建滤波器的正交族对应的三个滤波器间的关系式是正确的,请以此为基础,推导h0,h1的关系。
【解答】
①已知
g0n=-1n+1h1n
g1n=-1nh0n
g1n=-1n+1g02K-1-n
将第1、2个式子代入第3个式子得到
-1nh0n=-1n+1-12K-1-n+1h12K-1-n
h0n=-12K-n+1h12K-1-n
h0n=-1n+1h12K-1-n
②若
g0n=-1nh1n
g1n=-1n+1h0n
g1n=-1n+1g02K-1-n
则将第1、2个式子代入第3个式子得到
-1n+1h0n=-1n+1-12K-1-nh12K-1-n
h0n=-12K-1-nh12K-1-n
h0n=-1n+1h12K-1-n
因此,两种情况的解一致,都是:
h0n=-1n+1h12K-1-n
4、请证明完美重建滤波器组的双正交性质,即课本282页的公式7.1.20。
【解答】
①
②
③
已知书上公式7.1.9,7.1.10,7.1.18
H0zG0z+H1zG1z=2
H0-zG0z+H1-zG1z=0
H0zG0z+H0-zG0-z=2
调制矩阵
Hmz=H0zH0-zH1zH1-z
④
⑤
用调制矩阵的行列式表示G0z和G1z
G0z=2detHmzH1-z
G1z=-2detHmzH0-z
(1)先来证第1个公式:
g1k,h12n-k=δn
令Pz=G1zH1z,将⑤式代入得
Pz=G1zH1z=-2detHmzH0-zH1z
由于detHm-z=-detHmz,则将④式两端同时乘以H0z得
G0zH0z=2detHmzH0zH1-z=P-z
⑥
因此G0zH0z=G1-zH1-z,代入①式得
G1zH1z+G1-zH1-z=2
根据Z变换的Z域翻转性-1nxn⟺X-z,对上式反Z变换得到
kg1kh1n-k+-1nkg1kh1n-k=2δn
由于冲激函数δn在n=0时等于1,其他情况等于0,且n为奇数时,上式奇次方项可以相互抵消,因此只取n为偶数的情况,用2n代替上式的n可得
2kg1kh12n-k=2δ2n=2δn
kg1kh12n-k=g1k,h12n-k=δn
(2)再来证第2个公式:
g0k,h12n-k=0
由②式可得
H0-z=-H1-zG1zG0z
⑦
③-①得
H0-zG0-z-H1zG1z=0
将上面两式合并得
-H1-zG1zG0zG0-z-H1zG1z=0
H1-zG1zG0-z+H1zG1zG0z=0
由于G1z是滤波器,不为0,故可将上式两端同时除以G1z得
H1-zG0-z+H1z
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 图像 处理 分析 王伟强 作业题 标准答案 汇总