数字图像处理上机实习-频域处理.docVIP

数字图像处理上机实习 （第三专题——频域处理） 学生姓名： XX 班 级： XXX 学 号： XXXXXXX 指导老师： XXXXX 一、计算图象的傅氏变换频谱函数 要求：设计图象f2(x,y) 为30*256/256*256，居中垂直排列； 算法设计： 程序设计： f=zeros(256,256); f(1:256,113:143)=1; subplot(211);imshow(f);title(图像f2(x,y)); F=fft2(f); subplot(212);imshow(fftshift(uint8(F))); title(频谱图); figure; mesh(fftshift(abs(F))); title(傅氏变换频谱函数); 程序运行结果： 二、计算显示图象的频谱函数 要求：对p03-02的离散余弦变换，显示其频谱函数 算法设计： 程序设计： M=imread(D:\matlab2011\图片\p03-02.tif); subplot(211);imshow(M);title(原图像); M1=dct2(M); subplot(212);imshow(abs(M1));title(二维离散余弦变换结果); 程序运行结果： 设计图象的频域滤波 要求：根据频率采样法设计一个高通滤波器，对图象f2(x,y)进行高通滤波，观察分析空域图象和频谱分布的变化（设计2个截止频率）。 算法设计： 程序设计： %%%%%%%%%%%%%%%%%%%子函数%%%%%%%%%%%%%%%%%%%%%%%%% function [G,g] = imhp(f,D1) % 高通滤波 % f为输入图像,G为频谱图，g为复原图像 F=fftshift(fft2(f)); [m,n] = size(F); G=F; u0=floor(m/2); v0=floor(n/2); for u = 1:m for v = 1:n if ( sqrt((u-u0).^2+(v-v0).^2))= D1 H = 1; else H = 0; end G(u,v) = F(u,v)*H; end End g = ifft2(G); %%%%%%%%%%%%%%%%%%%主函数%%%%%%%%%%%%%%%%%%%%%%%%% f=zeros(256,256); f(80:173,113:143)=1; subplot(231);imshow(f);title(原图像); F1 = fft2(f); F1_shift = fftshift(F1); subplot(2,3,4);imshow(uint8(F1_shift));title(原频谱) %高通滤波 [F2,j] = imhp(f,50); g = abs(j); subplot(2,3,2);imshow(g);title(高通滤波后图像 D1=50) subplot(2,3,5);imshow(uint8(F2));title(高通后的频谱) [F2,j] = imhp(f,100); g = abs(j); subplot(2,3,3);imshow(g);title(高通滤波后图像 D1=100) subplot(2,3,6);imshow(uint8(F2));title(高通后的频谱) 程序运行结果： 四、实际图象的频域滤波 要求：根据频率采样法设计一个高通滤波器，对图象p03-02进行高通滤波，观察分析空域图象和频谱分布的变化（设计2个截止频率）。 算法设计： 程序设计： %%%%%%%%%%%%%%%%%%%子函数%%%%%%%%%%%%%%%%%%%%%%%%% function [G,g] = imhp(f,D1) % 高通滤波 % f为输入图像,G为频谱图，g为复原图像 F=fftshift(fft2(f)); [m,n] = size(F); G=F; u0=floor(m/2); v0=floor(n/2); %找到中心点 for u = 1:m for v = 1:n if ( sqrt((u-u0).^2+(v-v0).^2))= D1 H = 1; else H = 0; end G(u,v) = F(u,v)*H; end end g = ifft2