圖像處理作業第三次

學號:2019E8013261007
班級:705
姓名:蔡少斐

圖像處理作業第三次

1.根據書中對傅立葉變換的定義,證明課本165頁上有關傅立葉變換的平移性質。

F(uu0,vv0)F(u-u_0,v-v_0)

=x=0M1y=0N1f(x,y)ej2π((uu0)x/M+(vv0)y/N)=\sum_{x=0}^{M-1}\sum_{y=0}^{N-1}f(x,y)e^{-j2\pi((u-u_0)x/M+(v-v_0)y/N)}

=x=0M1y=0N1f(x,y)ej2π(u0x/M+v0y/N)ej2π(ux/M+vy/N)=\sum_{x=0}^{M-1}\sum_{y=0}^{N-1}f(x,y)e^{j2\pi(u_0x/M+v_0y/N)}e^{-j2\pi(ux/M+vy/N)}

=DFT(f(x,y)ej2π(u0x/M+v0y/N))=DFT(f(x,y)e^{j2\pi(u_0x/M+v_0y/N)})

f(xx0,yy0)f(x-x_0,y-y_0)

=u=0M1v=0N1F(u,v)ej2π((xx0)u/M+(yy0)v/N)=\sum_{u=0}^{M-1}\sum_{v=0}^{N-1}F(u,v)e^{j2\pi((x-x_0)u/M+(y-y_0)v/N)}

=u=0M1v=0N1f(x,y)ej2π(xu0/M+yv0/N)ej2π(xu/M+yv/N)=\sum_{u=0}^{M-1}\sum_{v=0}^{N-1}f(x,y)e^{-j2\pi(xu_0/M+yv_0/N)}e^{j2\pi(xu/M+yv/N)}

=IDFT(F(u,v)ej2π(xu0/M+yv0/N))=IDFT(F(u,v)e^{j2\pi(xu_0/M+yv_0/N)})

2. 課本171頁上習題4.9。

f(x,y)f(x,y)(1)(x+y)f(x,y) \rightarrow f(x,y)(-1)^{(x+y)}

則根據如下公式,可得

F(uM2,vN2)=f(x,y)(1)(x+y)ej2π(xu/M+yv/N)F(u-\frac{M}{2},v-\frac{N}{2}) = \sum\sum f(x,y)(-1)^{(x+y)} e^{-j2\pi(xu/M+yv/N)}

F(uM2,vN2)=DFT(f(x,y)(1)(x+y))....(a)\rightarrow F(u-\frac{M}{2},v-\frac{N}{2})=DFT(f(x,y)(-1)^{(x+y)}) ....(a)

進行共軛變換,可得

F(uM2,vN2)=F(M2u,N2v)F^*(u-\frac{M}{2},v-\frac{N}{2}) = F(\frac{M}{2}-u,\frac{N}{2}-v)

根據比例性 ,那麼(a)(a)式可以寫成:

F(M2u,N2v)=DFT(f(x,y)(1)(x+y))F(\frac{M}{2}-u,\frac{N}{2}-v) = DFT(f(-x,-y)(-1)^{(x+y)})

最後每個像素乘以(1)(x+y)(-1)^{(x+y)}得到

f(x,y)=f(Mx,Ny)f(-x,-y) = f(M-x,N-y)

由此,關於中心對稱。

3.證明高斯的傅立葉變換還是高斯函數。

已知:

eπ(x2+y2)=IDFT(eπ(u2+v2))e^{-\pi (x^2+y^2)} = IDFT(e^{-\pi (u^2+v^2)})

根據比例性

u12πσu,v12πσvu \leftarrow \frac{1}{\sqrt{2\pi}\sigma}u,v \leftarrow \frac{1}{\sqrt{2\pi}\sigma}v

那麼比例係數 a=b=12πσa = b = \frac{1}{\sqrt{2\pi}\sigma}

f(ax,by)=1abF(u/a,v/b)f(ax,by) = \frac{1}{|ab|}F(u/a,v/b)

即有

F(u,v)=IDFT(Ae(u2+v2)/2σ2)=A2πσ2eπ2σ2(x2+y2)F(u,v) =IDFT(Ae^{- (u^2+v^2)/2\sigma^2}) = A2\pi\sigma^2e^{- \pi2\sigma^2(x^2+y^2)}

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章