数字图像处理上机实验

数字图像处理上机作业 班级：020812 班 学号：******** 姓名：*** 2011/6/28 数字图像处理上机作业 1. 产生右图所示亮块图像 f1(x,y)（128×128 大小，暗处=0，亮处=255），对其进行FFT： （1） 同屏显示原图f1和 FFT(f1)的幅度谱图； 图像： （2） 若令 f2(x,y)=（-1）x+y f1(x,y)，重复以上过程，比较二者幅度谱的异同，简述理由； （3）若将 f2(x,y)顺时针旋转 45 度得到 f3(x,y)，试显示FFT(f3)的幅度谱，并与 FFT(f2)的幅度谱进行比较。 结论：不同点：f2 的频谱是对 f1 频谱的移位，它时 f1 的频谱从原点（0，0）移到了中心（64，64），而得到了一个完整的频谱。 相同点：频谱的实质没有改变，幅度等都没有发生变化。 （3）若将 f2(x,y)顺时针旋转 45 度得到 f3(x,y)，试显示FFT(f3)的幅度谱，并与 FFT(f2)的幅度谱进行比较。 源程序：f1=zeros(128,128); for i=38:1:90 for j=58:1:70 f1(i,j)=255; end end figure(1) subplot(1,2,1); imshow (f1); subplot(1,2,2); imshow (fft2(f1)); % f2(x,y)=（-1）^(x+y)* f1(x,y) for i=1:1:128 for j=1:1:128 f2(i,j)=(-1)^(i+j)*f1(i,j); end end figure(2); subplot(1,3,1); imshow (f2); f3=imrotate(f2,-45,'bilinear');%将f2 顺时针旋转45 度 subplot(1,3,2); imshow(fft2(f2));%显示f2 的频谱 subplot(1,3,3); imshow(fft2(f3));%显示f3 的频谱 结论：均衡化后的直方图并非完全均匀分布的原因：因为图像的像素个数和灰度等级均为离散值，而且均衡化后使灰度级并归。 2. 对256256 大小、256 级灰度的数字图像len a 进行频域的理想低通、高通滤波，同屏显示原图、幅度谱图和低通、高通滤波的结果图。 源代码： figure(1); fid=fopen('d:\img\lena.img','r'); data=(fread(fid,[256,256],'uint8'))'; subplot(2,2,1) imagesc(data); colormap(gray); title('LENA','Color','r'); subplot(2,2,2); imshow(fft2(data)); s=fftshift(fft2(data)); [M,N]=size(s); %分别返回s 的行数到M 中，列数到N 中 n=2; %对 n 赋初值 %GLPF 滤波，d0=5，15，30（程序中以 d0=30 为例） d0=30; %初始化 d0 n1=floor(M/2); %对 M/2 进行取整 n2=floor(N/2); %对 N/2 进行取整 for i=1:M for j=1:N d=sqrt((i-n1)^2+(j-n2)^2); %点（i,j）到傅立叶变换中心的距离 h=1*exp(-1/2*(d...