图形图像处理实验二.docVIP

甘肃政法学院 本科生实验报告 （二） 姓名:柳文会 学院:计算机科学学院 专业:计算机科学与技术 班级:2010级计本班 实验课程名称:图像变换及频域滤波 实验日期:2010年 4 月 12日 指导教师及职称:王秋云 实验成绩: 开课时间：2011-2012 学年 第二 学期 甘肃政法学院实验管理中心印制 实验题目 图像变换及频率滤波 小组合作 否 姓名 柳文会 班级 10计本班 学 号 201081010124 一、实验目的 （1）编写快速傅里叶变换算法程序，验证二维傅里叶变换的平移性和旋转不变性； （2）实现图像频域滤波，加深对频域图像增强的理解。 二．实验环境 Windws XP、MATALAB7.0 三、实验内容与步骤 实验内容：1．验证二维傅里叶变换的平移性和旋转不变性； a．要验证证其平移特性,就先建立一个二维图象,然后再对其平移,通过观察两者的频谱图来观察平移特性,为了方便起见,我们选择特殊情况来分析,令u0=v0=N/2,使??F(u-N/2,v-N/2),达到将原始F(U,V)四周频谱移到中心的效果,及达到频谱中心化。 b．验证旋转不变性可以通过将原始数组的通过移动45度,然后再比较旋转后与旋转前的频谱,得出频谱旋转不变性的结论。 具体步骤： 1）产生如图1所示图像（128×128大小，暗处=0，亮处=255） 2）同屏显示原图和的幅度谱图。 3）若令，重复以上过程，比较二者幅度谱的异同。 4）将顺时针旋转45度得到，显示的幅度谱，并与的幅度谱进行比较。 图1实验图象 图1实验图象f1(x, y) 实验内容：2.实现图像频域滤波，加深对频域图像增强的理解。频率域中进行增强是相当直观的,主要步骤有: 1）计算需要增强的图象的傅立叶变换; 2）将其与一个(根据需要设计的)转移的函数相乘; 3）再将结果反傅立叶变换以得到增强的图象. 为了直观的展示频域增强,可以通过下面任务来展现: 对如图2所示的数字图像lena.img（256×256大小、256级灰度）进行频域的理想低通、高通滤波，同屏显示原图、幅度谱图和低通、高通滤波的结果图。 图2 实验图象lena.img 图2 实验图象lena.img 四、实验过程与分析 1．验证二维傅里叶变换的平移性和旋转不变性； 1.1 产生的图像与fft2(f1)幅度谱的程序如下： 1.2 令，则图像与的幅度谱图如下： 1.3若将f1(x,y)顺时针旋转45°得到，则及的幅度谱图如下 a=zeros(128,128) for y=54:74 for x=34:94 a(x,y)=1; end end figure(1) a1=fft2(a); subplot(3,2,1); imshow(a); subplot(3,2,2); a2=abs(a1); mesh(a2); for x=1:128 for y=1:128 b(x,y)=(-1).^(x+y).*a(x,y); end end b1=fft2(b); subplot(3,2,3); imshow(b); subplot(3,2,4); b2=abs(b1); mesh(b2); t=imrotate(a,45,nearest,crop) t1=fft2(t); subplot(3,2,5); imshow(t); subplot(3,2,6); t2=abs(t1); surf(t2); 结果分析：对比两图可以得到,1.2得的图就是将1.1得的图的频谱往中心移.即称为频谱中心化,将能量集中的4个角往中心靠拢。 结果分析：将1.3得的图与1.1得的图比较可知,将原图移动旋转45度以后,幅度谱图仍然没的改变,图象能量依然集中在4个角。 2.实现图像频域滤波 2.1 对数字图像lena.img 进行频域的理想低通，同屏显示原图、幅度谱图和低通滤波的结果图。其中，取理想低通滤波的半径R分别为88、24、11和5。 程序代码如下（取r=8时） a=fopen(D:\图像实验 \img\lena.img,r); b=fread(a,[256,256],uchar); fclose(a); figure(1) subplot(1,2,1) imshow(b,[0,255]); b=fft2(b) m=abs(b); subplot(1,2,2) m0=15*log(m+1.001) surf(m) q=b; t=fftshift(q) r=8; for x=1:256 for y=1:256 if (x-1