从空间滤波器获得频率域滤波器

3-2从空间滤波器获得频率域滤波器一、实验目的掌握从空间滤波器获得频率域滤波器的方法二、实验内容1 .从Sobel空间滤波器用freqz2获得频率域滤波器2 .观察空间滤波器和对应的频率域滤波器的滤波结果三、实验步骤1 .执行如下代码，实现空间域滤波和频域滤波的比较1读入原始图像：clcclearf=imread(lena.jpg;imshow(ftitle(原始图像2进行傅里叶变换：F=fft2(f;S=fftshift(log(1+abs(F;S=gscale(S;imshow(Stitle(傅立叶频谱图像原始图像3使用f=im2double(f之后再进行处理的傅立叶频谱图像：f=im2double(f;%转换为F=fft2(f;S=fftshift(log(1+abs(F;S=gscale(S;imshow(Stitle(使用f=im2double(f之后再进行处理的傅立叶频谱图像2.freqz2P90增强垂直边缘sobelH=freqz2(h,PQ(1,PQ(2:1产生的滤波器原点在矩阵中心处读入图像：clcclearf=imread(lena.jpg;imshow(f进行傅里叶变换：F=fft2(f;S=fftshift(log(1+abs(F;S=gscale(S;Imshow(S增强垂直边缘：h=fspecial(sobel”;增强垂直边缘figure,freqz2(h;%usesn2n1=6464.PQ=paddedsize(size(f;H=freqz2(h,PQ(1,PQ(2;%产生的滤波器原点在矩阵中心处H1=ifftshift(H;%迁移原点到左上角figure,mesh(abs(H1(1:20:400,1:20:400imshow(abs(H,口imshow(abs(H1,口gs=imfilter(double(f,h;gf=dftfilt(f,H1;imshow(gs,口imshow(gf,口imshow(abs(gs,口imshow(abs(gf,口imshow(abs(gs0.2*abs(max(gs(:imshow(abs(gf0.2*abs(max(gf(:d=abs(gs-gf;max(d(:min(d(:3.reqz2P90增强水平边缘sobel%fft2(f产生的频域F的原点在左上角读入图像clcclearf=imread(lena.jpg;傅里叶变换：F=fft2(f;S=fftshift(log(1+abs(F;S=gscale(S;h=fspecial(sobel;%增强水平边缘figure,freqz2(h;%size_h=size(tempPQ=paddedsize(size(f;H=freqz2(h,PQ(1,PQ(2;H1=ifftshift(H;figure,mesh(abs(H1(1:20:400,1:20:400imshow(abs(H,口imshow(abs(H1,口gs=imfilter(double(f,h;gf=dffilt(f,H1;d=abs(gs-gf;max(d(:min(d(:四、实验总结通过本次实验，我了解了从空间滤波器获得频率域滤波器的方法，对于这部分的理论知识也有了更进一步的认识，通过在这几次实训I,我对MATLAB的操作更加熟悉了，整个实验进行的也更加顺利了