数字图像处理第三专题报告剖析.doc

中国地质大学（武汉） 数字图像处理上机实习 （第三专题） 姓 名： 班级序号： 学 号： 指导老师：傅华明 一 、实验题目 因为我座次号为13，所以做第4小题 一、计算图象的傅氏变换频谱函数 4，设计图象f4(x,y) 为30*120/256*256，居中垂直排列； 二、计算显示图象的频谱函数 4，对p03-04的离散哈达玛变换，显示其频谱函数 三、设计图象的频域滤波 4，根据频率采样法设计一个带阻滤波器，对图象f4(x,y)进行带阻滤波，观察分析空域图象和频谱分布的变化（设计2个阻带）。 四、实际图象的频域滤波 4，根据频率采样法设计一个带阻滤波器，对图象p03-04进行带阻滤波，观察分析空域图象和频谱分布的变化（设计2个阻带）。 实验内容 一、计算图象的傅氏变换频谱函数：设计图象f4(x,y) 为30*120/256*256，居中垂直排列 算法设计 程序代码 f = zeros(256,256); f(68:188,113:143) = 1; subplot(121); imshow(f); title(原图); F = fft2(f); F2=fftshift(F); F2 = log(abs(F2)); subplot(122); imshow(F2); title(频谱图); 运行结果 二、计算显示图象的频谱函数：对p03-04的离散哈达玛变换，显示其频谱函数 （1）程序代码 I=imread(E:\第03专题-频域增强\P03-04.tif); H=hadamard(256); I=double(I)/256; hI=H*I*H; hI=hI/256; subplot(121);imshow(I); title(原图); subplot(122);imshow(hI); title(图像的离散哈达玛变换); 运行结果 三、设计图象的频域滤波：根据频率采样法设计一个带阻滤波器，对图象f4(x,y)进行带阻滤波，观察分析空域图象和频谱分布的变化（设计2个阻带） 算法设计 程序代码 f = zeros(256,256); f(68:188,113:143) = 1; subplot(221); imshow(f); title(原图); F = fft2(f); F2=fftshift(F); subplot(222); imshow(log(abs(F2))); title(频谱图); [M,N]=size(F2); D0=50; D1=90; m=fix(M/2); n=fix(N/2); for i=1:M for j=1:N D=sqrt((i-m)^2+(j-n)^2); if(D=D0) h=1; else if(D=D1D=D0) h=0; else h=1; end end Q(i,j)=h*F2(i,j); end end subplot(223),imshow(Q); title(理想带阻滤波器); Q=ifftshift(Q); J1=ifft2(Q); J2=uint8(real(J1)); subplot(224),imshow(J2); title(复原图); （3）运行结果 实际图象的频域滤波：根据频率采样法设计一个带阻滤波器，对图象p03-04进行带阻滤波，观察分析空域图象和频谱分布的变化（设计2个阻带） 算法设计 同上一题 （2）程序代码 f = imread(E:\第03专题-频域增强\P03-04.tif); subplot(221); imshow(f); title(原图); f1=double(f); F = fft2(f1); F2=fftshift(F); subplot(222); imshow(log(abs(F2)),[-1 10]); title(频谱图); [M,N]=size(F2); D0=50; D1=90; m=fix(M/2); n=fix(N/2); for i=1:M for j=1:N D=sqrt((i-m)^2+(j-n)^2); if(D=D0) h=1; else if(D=D1D=D0)