資源描述：
《數(shù)字圖像處理代碼Ch3《頻率域?yàn)V波》.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、來自《數(shù)字圖像處理MATLAB版》書中示例，代碼復(fù)制與運(yùn)行結(jié)果。例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)無填充時(shí)頻率域中低通濾波后的圖像');%=======有填充的濾波==========PQ=paddedsize(size(f));%f是浮點(diǎn)數(shù)Fp=fft2(f,PQ(1),PQ(2));%計(jì)算填充的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)有填充時(shí)頻率域中低通濾波后的圖像');運(yùn)行結(jié)果：例3.2空間濾波和頻率域?yàn)V波的比較。clc,clear;closeall;%空間濾波和頻率域?yàn)V波比較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)相應(yīng)于垂直sobel空間濾波器的頻率域?yàn)V波器的絕對(duì)值');PQ=paddedsize(size(f));H=freqz2(h,PQ(1),PQ(2));H1=ifftshift(H);subplot222;mesh(abs(H1)');t
5、itle('(b)經(jīng)函數(shù)ifftshift處理后的同一濾波器');subplot223;imshow(abs(H),[]);title('以圖像方式展示兩個(gè)濾波器');subplot224;imshow(abs(H1),[]);title('以圖像方式展示兩個(gè)濾波器');gs=imfilter(f,h);%使用0填充圖像邊界gf=dftfilt(f,H1);figure;subplot221;imshow(gs,[]);title('(a)用垂直Sobel模板在空間域?qū)υ瓐D像濾波結(jié)果');subplot222;imshow(gf,[]);title('(b)用
6、濾波器b在頻率域中得到的結(jié)果');subplot223;imshow(abs(gs),[]);title('(c)a的絕對(duì)值');subplot224;imshow(abs(gf),[]);title('(d)b的絕對(duì)值');figure;%閾值處理，使主要邊緣顯示更為清晰subplot121;imshow(abs(gs)>0.2*abs(max(gs(:))));title('(a)a圖經(jīng)閾值處理的二值圖像');subplot122;imshow(abs(gf)>0.2*abs(max(gf(:))));title('(b)b圖經(jīng)閾值處理的二值圖像');運(yùn)行
7、結(jié)果：d=abs(gs-gf);max(d(:))min(d(:))例3.3函數(shù)dftuv的使用。clc,clear,closeall;%函數(shù)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))運(yùn)行結(jié)果：例3.4低通濾波器。clc,clear,closeall;f=imread('Fig0413(a)(original_test_pattern).tif');[f,revertclass]=tofloat(
8、f);PQ=paddedsize(si
