图形与图像处理实验七.doc

%读取图片 I = imread(’oldman.jpg); I = rgb2gray(l); %对图像进行傅里叶变换 F = fft2(l); %对变换后图像进行队数变化，并对其坐标平移，使其中心化 S = fftshift(log(1+abs(F))); %将频谱图象标度在0-256的范围内 S = gscale(S); %显示频谱图 figure(1)，imshow(S);titleC 频谱图?); 7. f 吻的 1 - □ X 文 nD MkD MOD BM 處雇①)BO(W) Mqd □ d J 嘛 Q QB MtftM %产生空间“sober模板 h = fspecialfsober); %查看相应频域滤波器图像 figure(2)，freqz2(h);titlef 相应频域滤波器图像?); 6 4 2O 1 6 4 2 O 1 %产生滤波时所需大小的矩阵 Pa = paddedsize(size(l)); %产生频域中的‘sobel’滤波器 H = freqz2(h，Pa(1)，Pa(2)); %重排数据序列，使得原点位于频率矩阵的左上角 H1=ifftshift(H); %以图形形式显示滤波器 figure(3)，subplot(1,2,1);imshow(abs(H},口) figure(3),subplot(1,2,2);imshow(abs(H1),0) ** Rgtre3 文邮 MkD .馨00 *AQ) IIKD IlClttO 贈助(M) □ X * □ Q %用模版h进行空域滤波 gs=imfilter(double(l),h); %用滤波器对图像进行频域滤波 gf=dftfilt(l,H1); figure(4),subplot(2,2,1),imshow(gs,[]) figure(4),subplot(2,2,2),imshow(gf,Q) figure(4),subplot(2,2,3),imshow(abs(gs),[]) figure(4),subplot(2,2,4),imshow(abs(gf),[]) %  %读取图片 f=imread(rose.jpg,); f=rgb2gray(f)； %产生滤波时所需大小的矩阵 Pa=paddedsize(size(f)); %设定高斯高通滤波器的阈值 D0=0.05*Pa(1); %产生高斯高通滤波器 H=hpfilter(gaussian,Pa(1),Pa(2),D0); %对图像进行滤波 g=dftfilt(f,H); %显不原图像 figure(5),subplot(1,2,1);imshow(f); titlef原图像); %显不滤波后图像 figure(5),subplot(1,2,2);imshow(g,[]); title figure(5),subplot(1,2,2);imshow(g,[]); title(?滤波后图像’); 7 Fijure 5 — □ X 文件(D病禳(D H(Y)踵入(D UKD蠲面(D)宙口QflO格助(tD □ U J k \ a DQ ■ o KilHft 濾波gri值