资源描述：
《图形与图像处理实验七》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、%读取图片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文n處雇①)BO(W)Mqd□dJ嘛QQBMtftM%产生空间“sober模板h=fspecialfsober);%查看相应频域滤波器
2、图像figure(2)，freqz2(h);titlef相应频域滤波器图像•);642O1%产生滤波时所需大小的矩阵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
3、*AQ)IIKDIlClttO贈助(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('ro
4、se.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(•滤波后图像’);
5、7Fijure5—□X文件(D病禳(DH(Y)踵入(DUKD蠲面(D)宙口QflO格助(tD□UJkaDQ■oKilHft濾波gri值
