实验三图像频域变换.docxVIP

实验三、图像的频域变换 一、 实验目的 1 了解图像变换的意义和手段； 2熟悉傅里叶变换的基本性质； 3热练掌握FFT方法及应用； 4通过实验了解二维频谱的分布特点； 5通过本实验掌握利用MATLAB编程实现数字图像的傅立叶变换及滤波锐化和复原处 理； 6 了解理想、巴特沃兹、高斯等不同滤波器的结构及滤波效果。 二、 实验原理 1应用傅立叶变换进行图像处理 傅里叶变换是线性系统分析的一个有力工具，它能够定量地分析诸如数字化系统、采样 点、电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能，将有助于 解决大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说，把时间 用在学习和掌握博里叶变换上是很有必要的。 2傅立叶(Fourier)变换的定义 对于二维信号，二维Fourier变换定义为： F(仏 V)=「「/u, y)^720+vv)^ J—J—8 ej0 =cos + j sin 二维离散傅立叶变换为: F(u,v) =M-l N-lx=0 v=0 F(u,v) = M-l N-l x=0 v=0 for u = 0J2…,M -1, v = 0,12…，N 一 1 M-1N-1 /(x, y)二 Y工 F(ti, for x = 0,1,2,...,M-\,y = 0,1,2,...,N-1 z/=0 v=0 图像的傅立叶变换与一维信号的傅立叶变换变换一样，有快速算法，具体参见参考书目， 有关傅立叶变换的快速算法的程序不难找到。实际上，现在有实现傅立叶变换的芯片，可以 实吋实现傅立叶变换。 3利用MATLAB软件实现数字图像傅立叶变换的程序： I=imread(4原图像名.gif); %读入原图像文件 imshow(I); fftl=fft2(l); sfftl=fftshift(fftl); RR=real(sfftl); II=imag(sfftI);%显示原图像 imshow(I); fftl=fft2(l); sfftl=fftshift(fftl); RR=real(sfftl); II=imag(sfftI); %二维离散傅立叶变换 %直流分量移到频谱中心 %収傅立叶变换的实部 %取傅立叶变换的虚部 A=sqrt(RR.A2+n.A2);% 计算频谱幅值 A= (A-min(min(A))) /(max(max(A))-min(min(A)))*225; %归一-化 figure;%设定窗I I figure; 三、 实验内容 （一）研究以下程序，分析程序功能：输入执行各命令行，认真观察命令执行的结果。 熟悉程序中所使用函数的调用方法，改变有关参数，观察试验结果。 1.傅立叶变换 简单人工二值图像 clear all, close all f = zeros(50,50); f(15:35,23:28)= 1; figure(l), imshow(f,notruesize*) F = fft2(f,128,128); Fl = fftshift(F); figure(2), imshow(log(abs(Fl)), [-1 5]); colormap(gray); colorbar figure(3), mesh(l:128,l:128, abs(Fl)); 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(l:128,1:128, abs(F3)); colormap(gray); colorbar  100120100 -8080606040、2040201001000J15050iiiiiiiiiihiHiiiiilllhiin140120100806040 100 120 100 - 80 80 60 60 40、 20 40 20 100 100 0J 150 50 iiiiiiiiiih iHiiiiilllhiin 140 120 100 80 60 40 20 0 150 120 100 实际图像傅立叶变换 clear all, close all A=imread(rice.tif); %读一幅灰度值图象 rice.tif imshow(A);FT imshow(A); FT 二 fft2(A); AbsFT=abs(FT); %求FT的绝对值，AbsFT为A的振幅谱 %可用以下方式观看振幅谱：