傅里叶变换

gaohaikuo1

1. clc;clear all
   
2. I=zeros(400,400);
   
3. for i=1:400
   
4.     for j=1:400
   
5.         I(i,j)=200*cos(j/10);
   
6.     end
   
7. end
   
8. subplot(1,3,1);imshow(I,[]);
   
9. F=fft2(I);
   
10. subplot(1,3,2);imshow(F,[]);
    
11. Fc=fftshift(F);
    
12. subplot(1,3,3);imshow(abs(Fc),[]);
    

复制代码

gaohaikuo1

1. clc;clear all
   
2. 
   
3. I= imread('lena.jpg')
   
4. I_gray = rgb2gray(I);
   
5. subplot(2,3,1);imshow(I_gray,[]);
   
6. F=fft2(I_gray); % fast floriur transform 快速傅里叶变换，
   
7. %空域变为频域，转换完之后，以原点（图像左上角对称）
   
8. subplot(2,3,2);imshow(uint8(log(abs(F)+1)),[]);
   
9. Fc=fftshift(F);%坐标原点移动到图像的中心
   
10. subplot(2,3,3);imshow(uint8(log(abs(Fc)+1)),[]);
    
11. F1=ifftshift(Fc);%坐标原点移动到图像的中心
    
12. subplot(2,3,4);imshow(uint8(log(abs(F1)+1)),[]);
    
13. I1=ifft2(F1);
    
14. subplot(2,3,5);imshow(abs(I1),[]);

复制代码