实验三 图像频域变换.doc

实验三、图像的频域变换 实验目的 1了解图像变换的意义和手段； 2熟悉傅里叶变换的基本性质； 3热练掌握FFT方法及应用； 4通过实验了解二维频谱的分布特点； 5通过本实验掌握利用MATLAB编程实现数字图像的傅立叶变换及滤波锐化和复原处理； 6 了解理想、巴特沃兹、高斯等不同滤波器的结构及滤波效果。 实验原理 1应用傅立叶变换进行图像处理 傅里叶变换是线性系统分析的一个有力工具，它能够定量地分析诸如数字化系统、采样点、电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能，将有助于解决大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说，把时间用在学习和掌握博里叶变换上是很有必要的。 2傅立叶（Fourier）变换的定义 对于二维信号，二维Fourier变换定义为：　　 　　二维离散傅立叶变换为： 图像的傅立叶变换与一维信号的傅立叶变换变换一样，有快速算法，具体参见参考书目，有关傅立叶变换的快速算法的程序不难找到。实际上，现在有实现傅立叶变换的芯片，可以实时实现傅立叶变换。 3利用MATLAB软件实现数字图像傅立叶变换的程序： I=imread(‘原图像名.gif’); %读入原图像文件 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; %设定窗口 imshow(A); %显示原图像的频谱 实验内容 （一）研究以下程序，分析程序功能；输入执行各命令行，认真观察命令执行的结果。熟悉程序中所使用函数的调用方法，改变有关参数，观察试验结果。 1．傅立叶变换 (1) 简单人工二值图像 clear all, close all f = zeros(50,50); f(15:35,23:28) = 1; figure(1), imshow(f,notruesize) F = fft2(f,128,128); F1 = fftshift(F); figure(2), imshow(log(abs(F1)), [-1 5]); colormap(gray); colorbar figure(3), mesh(1:128,1:128, abs(F1)); colormap(gray); colorbar F2 = fft2(imrotate(f,90),128,128); F3 = fftshift(F2); figure(4), imshow(imrotate(f,90),notruesize) figure(5), imshow(log(abs(F3)), [-1 5]); colormap(gray); colorbar figure(6), mesh(1:128,1:128, abs(F3)); colormap(gray); colorbar （2）实际图像傅立叶变换 clear all, close all A=imread(rice.tif); %读一幅灰度值图象.tif imshow(A); FT=fft2(A); %傅立叶变换 AbsFT=abs(FT); %求FT的绝对值，AbsFT为A的振幅谱 %可用以下方式观看振幅谱: F2=log(AbsFT); figure, imshow(F2,[-1 10],notruesize); FinvT=ifft2(FT); %傅立叶反变换 AbsFinvT=abs(FinvT); F3=log(AbsFinvT); figure, imshow(mat2gray(FinvT); (3) 频域滤波 close all, clear all bw = imread(text.tif); imshow(bw); C = real(ifft2(fft2(bw) .* fft2( fspecial(average),256,256))); figure, imshow(C,[]); %Display, scaling data to appropriate range. 2．离散余弦变换（DCT） (1) 对cameraman.tif图像的