資源描述:
《matlab傅立葉變換及圖像的頻域處理》由會員上傳分享,免費在線閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、實驗六傅立葉變換及圖像的頻域處理一、實驗?zāi)康?、理解離散傅立葉變換的基本原理;2、掌握應(yīng)用MATLAB語言進行FFT及逆變換的方法;3、熟悉圖像在頻域中處理方法,應(yīng)用MATLAB語言作簡單的低通濾波器。二、實驗原理1、傅立葉變換的基本知識2、MATLAB提供的快速傅立葉變換函數(shù)三、實驗要求1、讀取圖像girl.bmp,顯示這幅圖像,對圖像作傅立葉變換,顯示頻域振幅圖像。作傅立葉逆變換,顯示圖像,看是否與原圖像相同。A=imread('F:1girl.bmp','bmp');subplot(1,3,1),imshow(A);title('原圖
2、像');B=fftshift(fft2(A));subplot(1,3,2),imshow(log(abs(B)),[],'notruesize');title('二維傅立葉變換');C=ifft2(B);subplot(1,3,3),imshow(log(abs(C)),[],'notruesize');title('逆變換后圖像');2、設(shè)計一個簡單的理想低通濾波器(截止頻率自選),對圖像作頻域低通濾波,再作反變換,觀察不同的截止頻率下反變換后的圖像與原圖像的區(qū)別。A=imread('F:1girl.bmp','bmp');subplo
3、t(3,2,1),imshow(A);title('原圖像');B=fftshift(fft2(A));6E=ifft2(B);subplot(3,2,2),imshow(log(abs(E)),[],'notruesize');title('未濾波的逆變換圖像');[X,MAP]=imread('F:1girl.bmp','bmp');[m,n]=size(X);H1=zeros(m,n);H1(m/2:m,n/2:n)=1;LOWPASS1=B.*H1;C=ifft2(LOWPASS1);subplot(3,2,3),imshow(lo
4、g(abs(LOWPASS1)),[],'notruesize');title('低通1');subplot(3,2,4),imshow(log(abs(C)),[],'notruesize');title('低通1逆變換圖像');H2=zeros(m,n);H2(m/10:m,n/10:n)=1;LOWPASS2=B.*H2;D=ifft2(LOWPASS2);subplot(3,2,5),imshow(log(abs(LOWPASS2)),[],'notruesize');title('低通2');subplot(3,2,6),imshow
5、(log(abs(D)),[],'notruesize');title('低通2逆變換圖像');1、設(shè)計一個簡單的理想高通濾波器(截止頻率自選),對圖像作頻域高通濾波,再作反變換,觀察不同的截止頻率下反變換后的圖像與原圖像的區(qū)別。A=imread('F:1girl.bmp','bmp');subplot(3,2,1),imshow(A);title('原圖像');B=fftshift(fft2(A));E=ifft2(B);subplot(3,2,2),imshow(log(abs(E)),[],'notruesize');title('未
6、濾波的逆變換圖像');[X,MAP]=imread('F:1girl.bmp','bmp');6[m,n]=size(X);H1=ones(m,n);H1(m/3:m*2/3,n/3:n*2/3)=0;LOWPASS1=B.*H1;C=ifft2(LOWPASS1);subplot(3,2,3),imshow(log(abs(LOWPASS1)),[],'notruesize');title('高通1');subplot(3,2,4),imshow(log(abs(C)),[],'notruesize');title('高通1逆變換圖像')
7、;H2=ones(m,n);H2(m*3/5:m,n*3/5:n)=0;LOWPASS2=B.*H2;D=ifft2(LOWPASS2);subplot(3,2,5),imshow(log(abs(LOWPASS2)),[],'notruesize');title('高通2');subplot(3,2,6),imshow(log(abs(D)),[],'notruesize');title('高通2逆變換圖像');1、(選做內(nèi)容):對一幅圖像作傅立葉變換,顯示一幅頻域圖像的振幅分布圖和相位分布圖,分別對振幅分布和相位分布作傅立葉逆變換,觀察兩幅
8、逆變換后的圖像,體會頻域圖像中振幅與位相的作用。I=imread('girl.bmp');A=fftshift(fft2(I));A1=angle(A