實驗的3圖像變換

1、實用標準文案實驗三圖像變換一、實驗目的1 、熟悉及掌握圖像的變換原理及性質(zhì)，實現(xiàn)圖像的傅里葉變換。二、實驗環(huán)境MATLAB 6.5 以上版本 、WIN XP或 WIN2000計算機三、相關(guān)函數(shù)圖像的變換1 fft2fft2函數(shù)用于數(shù)字圖像的二維傅立葉變換，如：i=imread('e:w01.tif');j=fft2(i);由于 MATLAB無法顯示復數(shù)圖像，因此變換后的結(jié)果還需進行求模運算，即調(diào)用 abs 函數(shù)。之后常常還進行對數(shù)變換，即調(diào)用 log 函數(shù)，以減緩傅里葉譜的快速衰減，更好地顯示高頻信息。2 ifft2ifft2函數(shù)用于數(shù)字圖像的二維傅立葉反變換，如：i=imre

2、ad('e:w01.tif');j=fft2(i);k=ifft2(j);3 fftshift用于將變換后圖像頻譜中心從矩陣的原點移動到矩陣的中心B=fftshift（i ）4利用 fft2計算二維卷積精彩文檔實用標準文案利用 fft2函數(shù)可以計算二維卷積, 如:a=8,1,6;3,5,7;4,9,2;b=1,1,1;1,1,1;1,1,1;a(8,8)=0;b(8,8)=0;c=ifft2(fft2(a).*fft2(b);c=c(1:5,1:5);利用 conv2( 二維卷積函數(shù) ) 校驗 ,如：a=8,1,6;3,5,7;4,9,2;b=1,1,1;1,1,1;1,1,1

3、;c=conv2(a,b);四、實驗內(nèi)容1、對一幅圖像進行平移， 顯示原始圖像與處理后圖像， 分別對其進行傅里葉變換，顯示變換后結(jié)果，分析原圖的傅里葉譜與平移后傅里葉頻譜的對應關(guān)系。s=imread('f:1.jpg');%讀入原圖像i=rgb2gray(s)i=double(i)j=fft2(i);%傅里葉變換k=fftshift(j);%直流分量移到頻譜中心l=log(abs(k);%對數(shù)變換m=fftshift(j);%直流分量移到頻譜中心RR=real(m);%取傅里葉變換的實部II=imag(m);%取傅里葉變換的虛部A=sqrt(RR.2+II.2);%計算頻譜府幅

4、值精彩文檔實用標準文案A=(A-min(min(A)/(max(max(A)*255; %歸一化b=circshift(s,800 450); %對圖像矩陣 im 中的數(shù)據(jù)進行移位操作b=rgb2gray(b)b=double(b)c=fft2(b);%傅里葉變換e=fftshift(c);%直流分量移到頻譜中心l=log(abs(e);%對數(shù)變換f=fftshift(c);%直流分量移到頻譜中心WW=real(f);%取傅里葉變換的實部 BZZ=imag(f);%取傅里葉變換的虛部B=sqrt(WW.2+ZZ.2);%計算頻譜府幅值B=(B-min(min(B)/(max(max(B)*25

5、5;%歸一化subplot(2,2,1);imshow(s);title('原圖像 ')subplot(2,2,2);imshow(uint8(b);title('平移圖像 ')subplot(2,2,3);imshow(A);title('離散傅里葉頻譜 ');subplot(2,2,4);imshow(B);title('平移圖像離散傅里葉頻譜')精彩文檔實用標準文案2、對一幅圖像進行旋轉(zhuǎn)， 顯示原始圖像與處理后圖像， 分別對其進行傅里葉變換，顯示變換后結(jié)果，分析原圖的傅里葉譜與旋轉(zhuǎn)后傅里葉頻譜的對應關(guān)系。s=imread(&

6、#39;f:1.jpg');%讀入原圖像i=rgb2gray(s)i=double(i)j=fft2(i);%傅里葉變換k=fftshift(j);%直流分量移到頻譜中心l=log(abs(k);%對數(shù)變換m=fftshift(j);%直流分量移到頻譜中心RR=real(m);%取傅里葉變換的實部II=imag(m);%取傅里葉變換的虛部A=sqrt(RR.2+II.2);%計算頻譜府幅值A(chǔ)=(A-min(min(A)/(max(max(A)*255;%歸一化精彩文檔實用標準文案b=imrotate(s, -90); %對圖像矩陣 im 中的數(shù)據(jù)進行移位操作b=rgb2gray(b)b

7、=double(b)c=fft2(b);%傅里葉變換e=fftshift(c);%直流分量移到頻譜中心l=log(abs(e);%對數(shù)變換f=fftshift(c);%直流分量移到頻譜中心WW=real(f);%取傅里葉變換的實部 BZZ=imag(f);%取傅里葉變換的虛部B=sqrt(WW.2+ZZ.2);%計算頻譜府幅值B=(B-min(min(B)/(max(max(B)*255;%歸一化subplot(2,2,1);imshow(s);title('原圖像 ')subplot(2,2,2);imshow(uint8(b);title('平移圖像 ')s

8、ubplot(2,2,3);imshow(A);title('離散傅里葉頻譜 ');subplot(2,2,4);imshow(B);title('平移圖像離散傅里葉頻譜')精彩文檔實用標準文案傅里葉變換、反變換、I=imread('19.jpg');A=rgb2gray(I);x1=fft2(A);x2=fftshift(x1);x3=ifft(x1)/10;figure,subplot(1,4,1);imshow(A)title('原圖 ');subplot(1,4,2);imshow(x1)title('頻譜圖 &#

9、39;);subplot(1,4,3);imshow(log(abs(x2)+1),0 10);title('直流分量移至頻譜圖中心');subplot(1,4,4);imshow(x3,0 10)title('傅里葉反變換');DCT變換、反變換X=imread( '19.jpg');I=rgb2gray(X);subplot(1,3,1);精彩文檔實用標準文案imshow(I);title('原圖' );subplot(1,3,2);J=dct2(I);imshow(log(abs(J),0 20);title(' 二維離散余弦變換' );subplot(1,3,3);K=idct2(J)/20;imshow(K,0 20);title(' 二維離散反余弦變換' );利用 DCT變換壓縮圖像I=imread('19.jpg');A=rgb2gray(I);B=DCT2(A);B(abs(B)<0.1)=0;C=idct2(B)/255;figure,subp