数字图像处理频域增强供参习.docVIP

中国地质大学（武汉） 数字图像处理上机实习 （第三专题） 学生姓名： 班 级： 学 号： 指导老师： 实验内容 一 图计算图象的傅氏变换频谱函数 要求（1-6）： 设计图象f6(x,y) 为3*30*30/256*256，居中垂直排列，选用Matlab函数直接调用实现，重点观察空域图象和频域频谱的对应关系； 补充完成：设计120*30/256*256，观察空域图象和频域频谱的对应关系。 算法设计 2.程序代码 %观察空域图象和频域频谱的对应关系 %设计图象f6(x,y) 为3*30*30/256*256 f=zeros(256,256); f([30:60],[113:143])=1; f([90:120],[113:143])=1; f([150:180],[113:143])=1; subplot(221);imshow(f); % 设计图象f2(x,y)为120*30/256*256,并作fft变换 f2 = zeros(256,256); f2(114:143,69:188) = ones(30,120); subplot(223);imshow(f2); %二维傅里叶变换 F=fft2(f); F2 = fft2(f2); %绘制fft图 subplot(222);imshow(fftshift(log(abs(F)))); %title(频谱图) subplot(224);imshow(fftshift(log(abs(F2)))); %title(频谱图（量化）) figure subplot(121);mesh(fftshift(abs(F))); subplot(122);mesh(fftshift(abs(F2))); 结果分析 （1）空域图象和频域频谱对比 频谱图（量化）对比 clear; f=zeros(256,256); f([30:60],[113:143])=1; f([90:120],[113:143])=1; f([150:180],[113:143])=1; subplot(232);imshow(f); title(原图像) f1 = double(f); F1 = fft2(f1); %fft F2 = dct2(f1); %DCT F3 = (1/256) * hadamard(256)*f1*hadamard(256); %hadamard F4 = (1/256) * walsh(256)*f1*walsh(256); %walsh subplot(2,4,5);imshow(uint8(abs(F1)));title(fft变换图像) subplot(2,4,6);imshow(abs(F2));title(dct变换图像) subplot(2,4,7);imshow(abs(F3));title(hadamard变换图像) subplot(2,4,8);imshow(abs(F4));title(walsh变换图像) 离散余弦变换函数： function b=dct2(arg1,mrows,ncols) [m, n] = size(arg1); if (nargin == 1), if (m 1) (n 1), b = dct(dct(arg1).).; return; else mrows = m; ncols = n; end end % Padding for vector input. a = arg1; if nargin==2, ncols = mrows(2); mrows = mrows(1); end mpad = mrows; npad = ncols; if m == 1 mpad m, a(2, 1) = 0; m = 2; end if n == 1 npad n, a(1, 2) = 0; n = 2; end if m == 1, mpad = npad; npad = 1; end % For row vector. b = dct(a, mpad); if m 1 n 1, b = dct(b., npad).; end 运行结果 （1）离散余弦变换，显示其频谱函数 （2）实现离散傅立叶变换、离散余弦变换、Walsh变换和Hadama