圖像的傅立葉變換數(shù)字圖像處理實(shí)驗(yàn)報(bào)告

1、精選優(yōu)質(zhì)文檔-傾情為你奉上計(jì)算機(jī)與信息工程學(xué)院驗(yàn)證性實(shí)驗(yàn)報(bào)告專業(yè)：通信工程 年級/班級：2011級 20132014學(xué)年第一學(xué)期課程名稱數(shù)字圖像處理指導(dǎo)教師段新濤本組成員學(xué)號姓名實(shí)驗(yàn)地點(diǎn)計(jì)科樓111實(shí)驗(yàn)時(shí)間周五5-6節(jié)項(xiàng)目名稱圖像的傅立葉變換實(shí)驗(yàn)類型驗(yàn)證性一、實(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評價(jià)人眼對圖像幅頻特性和相頻特性的敏感度。二、實(shí)驗(yàn)原理1 應(yīng)用傅立葉變換進(jìn)行圖像處理傅里葉變換是線性系統(tǒng)分析的一個(gè)有力工具，它能夠定量地分析諸如數(shù)

2、字化系統(tǒng)、采樣點(diǎn)、電子放大器、卷積濾波器、噪音和顯示點(diǎn)等的作用。通過實(shí)驗(yàn)培養(yǎng)這項(xiàng)技能，將有助于解決大多數(shù)圖像處理問題。對任何想在工作中有效應(yīng)用數(shù)字圖像處理技術(shù)的人來說，把時(shí)間用在學(xué)習(xí)和掌握博里葉變換上是很有必要的。2 傅立葉（Fourier）變換的定義對于二維信號，二維Fourier變換定義為：逆變換：二維離散傅立葉變換為：逆變換：圖像的傅立葉變換與一維信號的傅立葉變換變換一樣，有快速算法，具體參見參考書目，有關(guān)傅立葉變換的快速算法的程序不難找到。實(shí)際上，現(xiàn)在有實(shí)現(xiàn)傅立葉變換的芯片，可以實(shí)時(shí)實(shí)現(xiàn)傅立葉變換。3利用MATLAB軟件實(shí)現(xiàn)數(shù)字圖像傅立葉變換的程序：I=imread(原圖像名.gif)

3、;%讀入原圖像文件imshow(I); %顯示原圖像fftI=fft2(I); %二維離散傅立葉變換sfftI=fftshift(fftI); %直流分量移到頻譜中心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); %顯示原圖像的頻譜三、實(shí)驗(yàn)步驟1 將圖像內(nèi)容讀入內(nèi)存；2 用Fourier變換算法，對圖像作二維Fourier變換；3 將其幅度譜進(jìn)行

4、搬移，在圖像中心顯示；4 用Fourier系數(shù)的幅度進(jìn)行Fourier反變換；5 用Fourier系數(shù)的相位進(jìn)行Fourier反變換；6 比較4、5的結(jié)果，評價(jià)人眼對圖像幅頻特性和相頻特性的敏感度。7 記錄和整理實(shí)驗(yàn)報(bào)告。四、實(shí)驗(yàn)儀器1計(jì)算機(jī)；2 MATLAB程序；3移動(dòng)式存儲(chǔ)器（軟盤、U盤等）。4 記錄用的筆、紙。五、實(shí)驗(yàn)報(bào)告內(nèi)容1敘述實(shí)驗(yàn)過程；2 提交實(shí)驗(yàn)的原始圖像和結(jié)果圖像。六、實(shí)驗(yàn)圖片number.tif七、實(shí)驗(yàn)程序及結(jié)果clc;clear allI=imread('Fig0707(a)(Original).tif.tif'); %讀入原圖像文件imshow(I); %

5、顯示原圖像title('原始圖像')fftI=fft2(I); %二維離散傅立葉變換sfftI=fftshift(fftI); %直流分量移到頻譜中心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); %顯示原圖像的頻譜title('原始圖像的頻譜')f1=ifft2(A); %用Fourier系數(shù)的幅度進(jìn)行Fourier反變換f2=ifft2(angle(fftI); %用Fourier系數(shù)的相位進(jìn)行Fourier反變換；figuresubplot 121;imshow(f1,)title('幅度進(jìn)行Fouri