图像变换及频域滤波.docVIP

实 验 报 告 实验名称 图像变换及频域滤波 课程名称 数字图像处理 姓名 成绩 班级 学号 日期 2011/10/30 地点 综合实验楼 备注： 1.实验目的 3 2.实验环境（软件、硬件及条件） 3 3. 实验方法 3 （1）傅立叶变换 3 （2）频域滤波 5 4. 实验分析 9 5. 实验结论 9 附件 9 1.实验目的 编写快速傅立叶变换算法程序，验证二维傅立叶变换的平移和旋转不变性； 实现图像频域滤波，加深对图像频域增强的理解。 2.实验环境（软件、硬件及条件） Windows2000/XP MATLAB 6.x or above Visual C++、Visual Basic 或其他 实验方法 （1）傅立叶变换 产生如图1所示图像f1 (x, y)（128×128 大小，黑色区域灰度值为0，中心白色区域大小为16×64，灰度值等于255： 图1 %源程序1 a=zeros(128,128); a(33:97,57:73)=255; imshow(a); （2）用MATLAB 中的fft2 函数对其进行FFT： ① 同屏显示原图f1和FFT( f1 )的幅度谱图； 图2 %源程序2 a=zeros(128,128); a(33:97,57:73)=255; F1=fft2(a); subplot(121); imshow(a); subplot(122); mesh(1:128,1:128,abs(F1)); ② 若令，重复以上过程，比较二者幅度谱的异同，简述理由； 图3 %源程序3 for x=1:128; for y=1:128; b(x,y)=(-1)^(x+y)*a(x,y); end end F2=fft2(b); subplot(121); imshow(b); subplot(122); mesh(1:128,1:128,abs(F2)); ③ 若将 f2 (x, y) 顺时针旋转45 度得到f3 (x, y)，试显示FFT( f3 ) 的幅度谱，并与FFT( f2 )的幅度谱进行比较。 图4 %源程序4 c=imrotate(a,-45,nearest,crop); for x=1:128; for y=1:128; c1(x,y)=(-1)^(x+y)*c(x,y); end end F=fft2(c1); F3=fftshift(F); subplot(121); imshow(c1); subplot(122); mesh(1:128,1:128,abs(F3)); （2）频域滤波 对如 图5 所示的数字图像lena.img（256×256 大小、256 级灰度）进行频域的 理想低通、高通滤波，同屏显示原图、幅度谱图和低通、高通滤波的结果图。 图5 ①原图频谱图 图6 %源程序6 file=fopen(C:\Users\thebluedb\Desktop\images\lena.img,r); lena=fread(file,[256,256],uint8); F=fft2(lena); F1=fftshift(F); F1temp=log(F1+1); mesh(0:255,0:255,abs(F1)); ②理想低通滤波 图7 %源程序7 figure; r = 88; %另取24、11、5 m=256;n=256; for i=1:m for j=1:n if (((i-m/2)^2+(j-n/2)^2) r^2 ) mask(i,j) = 1; else mask(i,j) = 0; end end end F2=F1.*mask; F2temp=log(F2+1); mesh(1:256,1:256,abs(F2)) ③低通滤波结果图 图8 %源程序8 Figure; ilena1=ifft2(ifftshift(F2)); imshow(ilena1,[0,255]); ④理想高通滤波 图9 %源程序9 figure; r = 2; %另取8、24 m=256;n=256; for i=1:m for j=1:n if (((i-m/2)^2+(j-n/2)^2) r^2 ) mask(i,j) = 0; else