實(shí)驗(yàn)四-圖像的傅立葉變換與頻域?yàn)V波11頁

1、實(shí)驗(yàn)四 圖像的傅立葉變換與頻域?yàn)V波一、 實(shí)驗(yàn)?zāi)康?了解圖像變換的意義和手段；2熟悉傅里葉變換的基本性質(zhì)；3熟練掌握FFT方法的應(yīng)用；4通過實(shí)驗(yàn)了解二維頻譜的分布特點(diǎn)；5通過本實(shí)驗(yàn)掌握利用MATLAB編程實(shí)現(xiàn)數(shù)字圖像的傅立葉變換。6、掌握怎樣利用傅立葉變換進(jìn)行頻域?yàn)V波7、掌握頻域?yàn)V波的概念及方法8、熟練掌握頻域空間的各類濾波器9、利用MATLAB程序進(jìn)行頻域?yàn)V波二、 實(shí)驗(yàn)原理1應(yīng)用傅立葉變換進(jìn)行圖像處理傅里葉變換是線性系統(tǒng)分析的一個(gè)有力工具，它能夠定量地分析諸如數(shù)字化系統(tǒng)、采樣點(diǎn)、電子放大器、卷積濾波器、噪音和顯示點(diǎn)等的作用。通過實(shí)驗(yàn)培養(yǎng)這項(xiàng)技能，將有助于解決大多數(shù)圖像處理問題。對任何想在工作中

2、有效應(yīng)用數(shù)字圖像處理技術(shù)的人來說，把時(shí)間用在學(xué)習(xí)和掌握博里葉變換上是很有必要的。2傅立葉（Fourier）變換的定義對于二維信號(hào)，二維Fourier變換定義為： 二維離散傅立葉變換為： 圖像的傅立葉變換與一維信號(hào)的傅立葉變換變換一樣，有快速算法，具體參見參考書目，有關(guān)傅立葉變換的快速算法的程序不難找到。實(shí)際上，現(xiàn)在有實(shí)現(xiàn)傅立葉變換的芯片，可以實(shí)時(shí)實(shí)現(xiàn)傅立葉變換。3利用MATLAB軟件實(shí)現(xiàn)數(shù)字圖像傅立葉變換的程序：I=imread(原圖像名.gif); %讀入原圖像文件imshow(I); %顯示原圖像fftI=fft2(I); %二維離散傅立葉變換sfftI=fftshift(fftI); %

3、直流分量移到頻譜中心RR=real(sfftI); %取傅立葉變換的實(shí)部II=imag(sfftI); %取傅立葉變換的虛部A=sqrt(RR.2+II.2);%計(jì)算頻譜幅值A(chǔ)=（A-min(min(A)）/(max(max(A)-min(min(A)*225; %歸一化figure; %設(shè)定窗口imshow(A); %顯示原圖像的頻譜域?yàn)V波分為低通濾波和高通濾波兩類，對應(yīng)的濾波器分別為低通濾波器和高通濾波器。頻域低通過濾的基本思想：G(u,v)=F(u,v)H(u,v)F(u,v)是需要鈍化圖像的傅立葉變換形式，H(u,v)是選取的一個(gè)低通過濾器變換函數(shù)，G(u,v)是通過H(u,v)減少F

4、(u,v)的高頻部分來得到的結(jié)果，運(yùn)用傅立葉逆變換得到鈍化后的圖像。理想地通濾波器(ILPF)具有傳遞函數(shù)：其中，為指定的非負(fù)數(shù)，為(u,v)到濾波器的中心的距離。的點(diǎn)的軌跡為一個(gè)圓。n階巴特沃茲低通濾波器(BLPF)(在距離原點(diǎn)處出現(xiàn)截至頻率)的傳遞函數(shù)為與理想地通濾波器不同的是，巴特沃茲率通濾波器的傳遞函數(shù)并不是在處突然不連續(xù)。高斯低通濾波器(GLPF)的傳遞函數(shù)為其中，為標(biāo)準(zhǔn)差。相應(yīng)的高通濾波器也包括：理想高通濾波器、n階巴特沃茲高通濾波器、高斯高通濾波器。給定一個(gè)低通濾波器的傳遞函數(shù)，通過使用如下的簡單關(guān)系，可以獲得相應(yīng)高通濾波器的傳遞函數(shù)：利用MATLAB實(shí)現(xiàn)頻域?yàn)V波的程序f=imr

5、ead('room.tif');F=fft2(f); %對圖像進(jìn)行傅立葉變換S=fftshift(log(1+abs(F);%對變換后圖像進(jìn)行隊(duì)數(shù)變化，并對其坐標(biāo)平移，使其中心化S=gscale(S); %將頻譜圖像標(biāo)度在0-256的范圍內(nèi)imshow(S) %顯示頻譜圖像h=special('sobel'); %產(chǎn)生空間sobel模版freqz2(h) %查看相應(yīng)頻域?yàn)V波器的圖像PQ=paddedsize(size(f);%產(chǎn)生濾波時(shí)所需大小的矩陣H=freqz2(h,PQ(1),PQ(2);%產(chǎn)生頻域中的sobel濾波器H1=ifftshift(H); %重

6、排數(shù)據(jù)序列，使得原點(diǎn)位于頻率矩陣的左上角imshow(abs(H),) %以圖形形式顯示濾波器figure,imshow(abs(H1),)gs=imfilter(double(f),h); %用模版h進(jìn)行空域?yàn)V波gf=dftfilt(f,H1); %用濾波器對圖像進(jìn)行頻域?yàn)V波figure,imshow(gs,) figure,imshow(gf,)figure,imshow(abs(gs),)figure,imshow(abs(gf),)f=imread('number.tif');%讀取圖片PQ=paddedsize(size(f);%產(chǎn)生濾波時(shí)所需大小的矩陣D0=0.05

7、*PQ(1); %設(shè)定高斯高通濾波器的閾值H=hpfilter('gaussian',PQ(1),PQ(2),D0);%產(chǎn)生高斯高通濾波器g=dftfilt(f,H); %對圖像進(jìn)行濾波figure,imshow(f) %顯示原圖像figure,imshow(g,) %顯示濾波后圖像三、 實(shí)驗(yàn)步驟1 生成如下圖所示的一個(gè)二維矩形信號(hào)。黑白2 利用一維FFT計(jì)算二維付里葉變換。分別顯示行計(jì)算結(jié)果和列變換結(jié)果。（立體結(jié)果，用mesh(F)顯示）3 利用MatLab工具箱中的函數(shù)編制FFT頻譜顯示的函數(shù);4 a). 調(diào)入、顯示“實(shí)驗(yàn)一”獲得的圖像；圖像存儲(chǔ)格式應(yīng)為“.gif”; b)

8、 對這三幅圖像做FFT并利用自編的函數(shù)顯示其頻譜;c) 討論不同的圖像內(nèi)容與FFT頻譜之間的對應(yīng)關(guān)系。5 利用MATLAB提供的低通濾波器實(shí)現(xiàn)圖像信號(hào)的濾波運(yùn)算，并與空間濾波進(jìn)行比較。6利用MATLAB提供的高通濾波器對圖像進(jìn)行處理。7 記錄和整理實(shí)驗(yàn)報(bào)告。四、 實(shí)驗(yàn)報(bào)告內(nèi)容1敘述實(shí)驗(yàn)過程；2提交實(shí)驗(yàn)的原始圖像和結(jié)果圖像。五、 思考題1傅里葉變換有哪些重要的性質(zhì)?2圖像的二維頻譜在顯示和處理時(shí)應(yīng)注意什么?3.用數(shù)據(jù)和圖片給出各個(gè)步驟中取得的實(shí)驗(yàn)結(jié)果，并進(jìn)行必要的討論，必須包括原始圖像及其計(jì)算/處理后的圖像。4.結(jié)合實(shí)驗(yàn)，評(píng)價(jià)頻域?yàn)V波有哪些優(yōu)點(diǎn)？5.在頻域?yàn)V波過程中需要注意哪些事項(xiàng)？ clc;c

9、lear;I=zeros(256,256);for i=118:138 for j=118:138 I(i,j)=1; endend% imshow(I); m,n=size(I);X,Y=meshgrid(1:m,1:n);J=I.*(-1).(X+Y);% H=fft(J);figure,mesh(log(1+abs(H);title('行變換'); L=fft(J')'figure,mesh(log(1+abs(L);title(列變換);clc;clear;I=zeros(256,256);for i=118:138 for j=118:138 I(i,j)=1; endend% imshow(I);I=fft(I);J=fft(I')'J=fftshift(J);imshow(log(1+abs(J),) clc;clear;I=zeros(256,256);for i=118:138 for j=118:138 I(i,j)=1