最新实验一图像的二维离散傅立叶变换.docxVIP

实验一图像的二维离散傅立叶变换 实验目的 掌握图像的二维离散傅立叶变换以及性质 实验要求 实验要求 64的黑色1）建立输入图像，在 64的黑色 图像矩阵的中心建立16 16的 白色矩形图像点阵，形成图像文件。对输入图像进行二维傅立叶变换， 将原始图像及变换图像（三维、中心化）都显示于屏幕上。 2） 调整输入图像中白色矩形的位置，再进行变换，将原始图像及变换图像（ 维、中心化）都显示于屏幕上，比较变换结果。 3） 调整输入图像中白色矩形的尺寸 40， 4 4），再进行变换，将原始图 像及变换图像（三维、中心化）都显示于屏幕上，比较变换结果。 三、实验仪器设备及软件 HP D538、MATLAB 四、实验原理 在二维情况下，定义f(x,y)的傅立叶变换F (u,v): F(u v) f (x, y)ej2 (ux vy)dxdy f(x,y) F(u,v)ej2 (uxvy)dudv 它表明了空间频率成分与二维图像信号之间的相互关系 对于我们要处理的实际二维图像，其傅氏变换一般是在频率域上有界的， 亦即有 用成分总是落在一定的频率域范围之内 上述的频率域性质的依据在于： 一是图像中景物的复杂性具有一定的限度，其中大部分内容是变化不大的区域 完全像“雪花”点似的图像没有任何实际意义。 二是人眼对空间复杂性(频率)的分辨率以及显示器的分辨能力都是具有一定 限度。 若实变量函数f(x)是绝对可积的，即： 且F(u)是可积的，则傅立叶变换对一定存在。 F f x F u f x exp j 2 ux dx F F u f x F u exp j 2 ux du 如果f(x)考虑为实函数，它的傅立叶变换通常是复数形式，即： Fu Ru jI u也可表为： F u |F u |ej u 若二变量函数f(x,y)是绝对可积的，即： 且F(u,v是可积的，则傅立叶变换对一定存在。 F f x,y F u, v f x, y exp j 2 ux vy F -1 F u , v f x, y F u , v exp j 2 ux vy dudv 二维函数的傅立叶谱，振幅谱 相位谱 和能量谱分别为： 2 2 1/2 Fu,v| R u,v I u , v I u , v u, v arctg R u, v Eu,v R 2 u,v I 2 u , v 五、实验步骤及程序 %clear 9原始图象 f=zeros(64,64); %俞入64*64的黑色图像矩阵 f(25:40,25:40)=1; %建立16*16的白色矩行图像点阵 figure(1); subplot(231),imshow(f); title(原始图像’)％显示原图像 F=fft2(f); %傅立叶变换 subplot(232) imshow(abs(F));title(傅里叶变换图像’);％显示傅里叶变换图像 F2=fftshift(abs(F)); %频谱中心化 subplot(233); imshow(abs(F2));title(中心化傅里叶频谱图’);％显示中心化傅里叶频谱图 x=1:64; y=1:64; subplot(234); mesh(abs(real(F)));title(三维频谱图’);％显示三维频谱图 subplot(235) mesh(x,y,F2(x,y)); title( FFT) 2、调整输入图像中白色矩形的位置，再进行变换后的程序 clear %原始图象 f=zeros(64,64); %俞入64*64的黑色图像矩阵 f(47:63,47:63)=1; %建立16*16的白色矩行图像点阵 figure(1); subplot(231),imshow(f); title(原始图像’)％显示原图像 F=fft2(f); %傅立叶变换 subplot(232) imshow(abs(F));title(傅里叶变换图像’);％显示傅里叶变换图像 F2=fftshift(abs(F)); %频谱中心化 subplot(233); imshow(abs(F2));title(中心化傅里叶频谱图’);％显示中心化傅里叶频谱图 x=1:64; y=1:64; subplot(234); mesh(abs(real(F)));title(三维频谱图’);％显示三维频谱图 subplot(235) mesh(x,y,F2(x,y)); title( FFT) 3 、 整 输 入 图 像 中 白 色 矩 形 的 尺 40( 40 40 440 X 4 40 X40 clear %原始图象 Irzeros(64b4)-雪 64*64sllB??wa再 f(」3ul2L3ul2llr 禧日 16U6 孚 口?answw再 figures八 subpof(231) jm