[高中教育]实验三 频域增强.pptVIP

实验三 傅里叶变换及频域增强 一．实验内容 二．实验目的 三．实验步骤 3.1傅里叶变换性质 3.2低通滤波 3.3高通滤波 二．实验目的 1、理解傅里叶变换的原理，掌握傅里叶变换的性质 2、掌握频域平滑原理，学会用理想低通滤波器、Butterworth低通滤波器、高斯低通滤波器进行图像处理。 3、掌握频域锐化原理，学会用理想高通滤波器、Butterworth高通滤波器、高斯高通滤波器进行图像处理。 三．实验步骤 3.1傅里叶变换性质 一个图像尺寸为 的函数 的离散傅里叶变换有下列等式给出： 其中， 相应地，可以通过傅里叶反变换求的 ，由表达式： 其中， 通常在进行傅里叶变换之前用 乘以输入的 图像函数。由于指数的性质，有以下表达式： 该表达式说明，用 乘以 将 原点即 被设置在 和 。 实验步骤：1.首先构造一幅黑白二值图像，在128×128 的黑色背景中心产生一个4×4的白色方块，对其进行傅里 叶变换；(Matlab中用fft2实现2D傅里叶变换) 2.把低频分量移到图象中心，而把高频分量移到四 个角上；(方法有两种：其一，在FT以前对测试图 象逐点加权(-1)^(i+j);其二，利用FFTSHIFT数）; 3.利用图象增强中动态范围压缩的方法增强2DFT；（Y＝C*log（1＋abs（X))）； 4 构造一幅黑白二值图像，在128×128的黑色背景中令第32行至36行、第32列至第36列的值为1（即产生一个4×4的白色方块），对其进行傅里叶变换； 5 将上图旋转300，再进行傅里叶变换 (imrotate) 构造二幅黑白二值图像，在128×128的黑色背景中分别令第60行至68行、第60列至第68列的值为1，第64行至65行、第64列至第65列的值为1产生两幅图像，分别对这两幅图像进行傅里叶变换. 用matlab编程，其程序如下： clear all; close all; clc; f=zeros(128); f(63:66,63:66)=1; g=fft2(f); m=fftshift(g); y=log(1+abs(m)); f1=zeros(128); f1(32:36,32:36)=1; g2=imrotate(f1,30,bilinear,crop); g3=fft2(g2); g1=fft2(f1); f3=zeros(128);f4=zeros(128); f3(60:68,60:68)=1; f4(64:65,64:65)=1; g4=fft2(f3);g5=fft2(f4); subplot(2,2,1),imshow(f);title(原始图像); subplot(2,2,2),imshow(g);title(未平移傅里叶变换); subplot(2,2,3),imshow(m);title(低频及高频转移); subplot(2,2,4),imshow(y);title(对数变换); figure,subplot(2,2,1),imshow(f);title(原始图像); subplot(2,2,2),imshow(fftshift(abs(g)),[-1,5]);title(平移后傅里叶变换); subplot(2,2,3),imshow(f1);title(原始图像); subplot(2,2,4),imshow(fftshift(abs(g1)),[-1,5]);title(平移后傅里叶变换); figure,subplot(2,2,1),imshow(f);title(原始图像); subplot(2,2,2),imshow(fftshift(abs(g)),[-1,5]);title(平移后傅里叶变换); subplot(2,2,3),imshow(g2);title(旋转30); subplot(2,2,4),imshow(fftshift(abs(g3)),[-1,5]);title(平移后傅里叶变换); figure,subplot(2,3,1),imshow(f4);title(原始图像); subplot(2,3,2),imshow(fftshift(abs(g5)),[-1,5]);t