實驗六---圖像陷波處理3頁

1、數(shù)字圖像處理實驗報告姓名：田蕾 學號：20091202098 專業(yè)：信號與信息處理 年級：09 實驗六 圖像陷波處理一、 實驗目的掌握巴特沃斯陷波濾波器的設計。進一步加深理解和掌握圖像頻譜的特點和陷波濾波的原理。理解圖像陷波濾波的處理過程和特點。二、 實驗內容設計程序，實現(xiàn)二階巴特沃斯陷波濾波器對圖像特定干擾頻率的濾波處理。觀察處理前后圖像效果，分析實驗結果和算法特點。（陷波濾波器半徑根據圖像干擾頻率的特點由實驗選定）三、 實驗原理n階 Butterworth 陷波器的傳遞函數(shù)為： 其中D0為陷波器半徑， D1， D2分別為點 （u,v） 到陷波器中心 （u0,v0） 和 （-u0,-v0）的

2、距離。該傳遞函數(shù)使得以陷波器中心為圓心，D0為半徑的圓內的所有u,v對應的濾波函數(shù)H (u,v)輸出為 0，也即將位于（u0,v0）和（-u0,-v0）處的頻率成分去除。四、 算法設計（含程序設計流程圖）獲取圖像的頻域圖像，使用Matlab函數(shù)pixval觀察干擾噪聲的頻率成分，找出所有的（u0,v0）對，針對每一個噪聲源分別設計相應的陷波器。濾波完成后將頻域圖像實施逆變換，即可顯示處理后的圖像。 五、 實驗結果及分析（需要給出原始圖像和處理后的圖像） 實驗結果分析：陷波濾波器阻止事先定義的中心頻率鄰域內的頻率，本實驗利用pixval函數(shù)得出相應地（u0，v0），設置相應地陷波半徑，以合理的程

3、度屏蔽噪聲，但不會帶來模糊，因此濾波效果比較的好。 附：程序源代碼r=imread('D:實驗6原圖.tif');m,n=size(r);imshow (r);title('原圖象');y=fft2(r);y0=abs(y);figure,imshow(y0,);title('圖像頻譜');g=fftshift(y);figure,imshow(abs(g),);title('中心移位后頻譜');y1=log(1+abs(g);figure,imshow(y1,);title('取對數(shù)變換后頻譜');pixval;

4、for u=1:m for v=1:nd11=(u-257)2+(v-267)2)(1/2);d12=(u-257)2+(v-247)2)(1/2);d10=3;H1(u,v)=1/(1+(d10*d10/(eps+d11*d12)2); endendg1=g.*H1;for u=1:m for v=1:nd21=(u-237)2+(v-257)2)(1/2);d22=(u-277)2+(v-257)2)(1/2);d20=3;H2(u,v)=1/(1+(d20*d20/(eps+d21*d22)2); endendg2=g1.*H2;for u=1:m for v=1:nd31=(u-217)2+(v-217)2)(1/2);d32=(u-297)2+(v-297)2)(1/2);d30=3;H3(u,v)=1/(1+(d30*d30/(eps+d31*d32)2); endendg3=g