资源描述：
《数字图像处理代码Ch3《频率域滤波》.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、来自《数字图像处理MATLAB版》书中示例，代码复制与运行结果。例3.1有和没有填充的滤波效果。clc,clear;closeall;%有和没有填充的滤波效果f=imread('Fig0405(a)(square_original).tif');%========无填充的滤波========[M,N]=size(f);[f,revertclass]=tofloat(f);F=fft2(f);sig=10;H=lpfilter('gaussian',M,N,sig);G=H.*F;g=ifft2(G);g=revertclass(g);subplot221;im
2、show(f);title('(a)尺寸为256×256像素的一幅图像');subplot222;imshow(g);title('(b)无填充时频率域中低通滤波后的图像');%=======有填充的滤波==========PQ=paddedsize(size(f));%f是浮点数Fp=fft2(f,PQ(1),PQ(2));%计算填充的FFTHp=lpfilter('gaussian',PQ(1),PQ(2),2*sig);Gp=Hp.*Fp;gp=ifft2(Gp);gpc=gp(1:size(f,1),1:size(f,2));gpc=revertcla
3、ss(gpc);subplot223;imshow(gp);title('(c)图像滤波后用ifft2得到全填充图像');h=fspecial('gaussian',15,7);gs=imfilter(f,h);subplot224;imshow(gs);title('(d)有填充时频率域中低通滤波后的图像');运行结果：例3.2空间滤波和频率域滤波的比较。clc,clear;closeall;%空间滤波和频率域滤波比较f=imread('Fig0409(a)(bld).tif');f=tofloat(f);F=fft2(f);S=fftshift(log(1
4、+abs(F)));subplot121;imshow(f);title('(a)一幅灰度图像');subplot122;imshow(S,[]);title('(b)该图像的傅里叶谱');h=fspecial('sobel');%生成sobel滤波器figure;subplot221;freqz2(h);title('(a)相应于垂直sobel空间滤波器的频率域滤波器的绝对值');PQ=paddedsize(size(f));H=freqz2(h,PQ(1),PQ(2));H1=ifftshift(H);subplot222;mesh(abs(H1)');t
5、itle('(b)经函数ifftshift处理后的同一滤波器');subplot223;imshow(abs(H),[]);title('以图像方式展示两个滤波器');subplot224;imshow(abs(H1),[]);title('以图像方式展示两个滤波器');gs=imfilter(f,h);%使用0填充图像边界gf=dftfilt(f,H1);figure;subplot221;imshow(gs,[]);title('(a)用垂直Sobel模板在空间域对原图像滤波结果');subplot222;imshow(gf,[]);title('(b)用
6、滤波器b在频率域中得到的结果');subplot223;imshow(abs(gs),[]);title('(c)a的绝对值');subplot224;imshow(abs(gf),[]);title('(d)b的绝对值');figure;%阈值处理，使主要边缘显示更为清晰subplot121;imshow(abs(gs)>0.2*abs(max(gs(:))));title('(a)a图经阈值处理的二值图像');subplot122;imshow(abs(gf)>0.2*abs(max(gf(:))));title('(b)b图经阈值处理的二值图像');运行
7、结果：d=abs(gs-gf);max(d(:))min(d(:))例3.3函数dftuv的使用。clc,clear,closeall;%函数dftuv的使用[U,V]=dftuv(8,5);DSQ=U.^2+V.^2fftshift(DSQ)D=hypot(U,V)timeit(@()U.^2+V.^2)timeit(@()hypot(U,V))运行结果：例3.4低通滤波器。clc,clear,closeall;f=imread('Fig0413(a)(original_test_pattern).tif');[f,revertclass]=tofloat(
8、f);PQ=paddedsize(si
