傅立叶变换

  • 格式:doc
  • 大小:1.45 MB
  • 文档页数:6

实验六——图像的傅立叶变换
一、实验目的
1、了解图像变换的意义和手段;
2、熟悉傅立叶变换的基本性质;
3、熟悉掌握FFT变换方法及应用;
4、通过实验了解二维频谱的分布方法;
5、通过本实验掌握利用MA TLAB编程实现数字图像的傅立叶变换;
6、评价人眼对图像幅频特性和相频特性的敏感度。

二、实验内容
I=imread ('原图像名.gif');
I=rgb2gray(I);
imshow(I);
fftI=fft2(I);
sfftI=fftshift(fftI);
RR=real(sfftI);
II=imag(sfftI);
A=sqrt(RR.^2+II.^2);
A=(A-min(min(A)))/(max(max(A))-min(min(A)))*225
figure(5111),imshow(A);
f=imread('2.jpg');
f=rgb2gray(f);
[M,N]=size(f);
F=fft2(f);
sig=10;
H=lpfilter('gaussian',M,N,sig);
G=H.*F;
g=real(ifft2(G));
figure(5112),imshow(g,[])
PQ=paddedsize(size(f));
Fp=fft2(f,PQ(1),PQ(2));
Hp=lpfilter('gaussian',PQ(1),PQ(2),2*sig); Gp=Hp.*Fp;
gp=real(ifft2(Gp));
gpc=gp(1:size(f,2),1:size(f,2));
figure(5113),imshow(gp,[])
f=imread('3.jpg');
f=rgb2gray(f);
F=fft2(f);
S=fftshift(log(1+abs(F)));
S=gscale(S);
figure(5114),imshow(S);
h=fspecial('sobel')'
h =
1 0 -1
2 0 -2
1 0 -1
freqz2(h)
PQ=paddedsize(size(f));
H=freqz2(h,PQ(1),PQ(2));
H1=ifftshift(H);
imshow(abs(H),[])
figure(5115),imshow(abs(H1),[])
gs=imfilter(double(f),h);
gf=dftfilt(f,H1);
imshow(gs,[])
figure(5116),imshow(gf,[])
figure(5117),imshow(abs(gs),[])
figure(5118),imshow(abs(gf),[])
figure(5119),imshow(abs(gs)>0.2*abs(max(gs(:)))) figure(51110),imshow(abs(gf)>0.2*abs(max(gf(:))))
三、 实验小结
通过这次学习,对图像变换的意义和手段、傅立叶变换的基本性质有了一定了解 ,基本掌握了FFT 变换方法及应用,也通过实验了解了二维频谱的分布方法,以及学会掌握利用MATLAB 编程实现数字图像的傅立叶变换,评价人眼对图像幅频特性和相频特性的敏感度。