数字图像处理实验4空间域和频率域图像处理.doc

电子信息学院 实验报告书 课 程 名： 《数字图像处理》课程实验 1、实验内容或题目 题目 对于二维信号，二维傅立叶变换定义为： 其逆变换： 二维离散傅立叶变换为： 其逆变换： 2、实验目的与要求 熟悉图像空间域滤波的MATLAB函数； 熟练掌握应用MATLAB软件编程进行空间域滤波的方法。 了解频域图像处理的意义和手段； 熟悉离散傅里叶变换的基本性质； 熟练掌握图像傅里叶变换的方法及应用； 了解二维频谱的分布特点； 掌握利用MATLAB工具箱实现频域图像处理 要求 将直方图增强和空间域滤波结果打印出来提交； 提交程序源代码M文件； 比较直方图增强前后图像的视觉效果有何不同。 3、实验仪器 MATLAB软件一套； 移动存储器（U盘等）； 实验图像若干。 4、实验图像 5、实验步骤 打开计算机并启动MATLAB程序，在“work”文件夹下找到待处理的图像文件； 统计原的灰度直方图，并利用直方图均衡方法进行增强； 记录和整理实验报告。 打开计算机并启动MATLAB程序； 读入test.bmp图像，并计算其傅立叶变换的模值； 分别构造gaussian滤波器（,）和motion滤波器（，）； 利用上一步骤构造的两个滤波器分别对test.bmp的频域模值进行滤波处理； 将原图像和频域滤波的结果显示出来； 重新构造一幅黑白二值图像，在的黑色背景中心一个的白色方块 了解变换为可分离核变换的性质，利用1D FFT子程序实现2D变换； 把低频分量移到中心，而把高频分量移到四个角上（利用函数） 利用增强中动态范围压缩的方法增强2DFT（） 研究变换幅度的平移不变性； 研究变换幅度的旋转性质； 研究变换的尺度变换性质 记录和整理实验报告。 6、程序代码 clear clc a=imread(girl2.jpg);%读取图像 b=rgb2gray(a);%转换成灰度图像 c=histeq(b,256); figure subplot(221),imshow(a);title(原图); subplot(222),imshow(c);title(均衡化图); subplot(223),imhist(b);title(原图灰度直方图); subplot(224),imhist(c);title(均衡化直方图); d=fspecial(gaussian,3,1);%构造gaussian滤波器（hsize=3,sigma=1） e=fspecial(motion,9,45);%构造motion滤波器（len=9，theta=45） f=imfilter(a,d);%gaussian滤波处理 g=imfilter(a,e);%motion滤波处理 figure subplot(221),imshow(a);title(原图); subplot(222),imshow(f);title(高斯滤波); subplot(223),imshow(g);title(运动平衡滤波); %频率域图像处理 A=zeros(128); A(62:64,62:66)=1;%构造一幅黑白二值图像 figure F1=fft2(A,128,128);%2D傅里叶变换 F2=fftshift(F1);%将原点移至中心 F3=log(1+abs(F2)); A1=zeros(128); A1(22:26,22:26)=1; F4=fft2(A1,128,128); A2=imrotate(A,-45);%旋转45度 F5=fft2(A2,128,128);%旋转45度后的傅里叶变换 A3=zeros(128); A3(60:68,60:68)=1; F6=fft2(A3,128,128); subplot(341),imshow(A);title(原图); subplot(342),imshow(F1);title(2D傅里叶变换图); subplot(343),imshow(F2);title(fftshift); subplot(344),imshow(F3);title(增强图); subplot(345),imshow(A1);title(平移图); subplot(346),imshow(F4);title(平移后的傅里叶变换图); subplot(347),imshow(A2);title(旋转45度图); subplot(348),imshow(F5);title(旋转45度后的傅里叶变换图); subplot(349),imshow(A3);title(尺度变换图); subplot(3,4,10),imshow(F6);title(尺度变换后的傅里叶变换图); 7、结果