“數(shù)字圖像處理與識別”實驗報告-實驗6-圖像帶阻濾波

1、“數(shù)字圖像處理與識別”實驗報告實驗六圖像帶阻濾波一、實驗目的掌握巴特沃斯帶阻濾波器的設計。進一步加深理解和掌握圖象頻譜的特點和帶阻濾波的原理。理解圖象帶阻濾波的處理過程和特點。二、實驗內容設計程序，實現(xiàn)二階巴特沃斯帶阻濾波器對圖像特定干擾頻率的濾波處理。觀察處理前后圖象效果，分析實驗結果和算法特點。(帶阻濾波器半徑和寬度根據(jù)圖像干擾頻率的特點由實驗選定)三、實驗原理N階Butterworth帶阻濾波器傳遞函數(shù)：H(u,v)=皐計丄十02(“)_詁其中D(u,v)是圖像頻譜圖中的點到中心的距離，W是阻帶寬度，6是阻帶的中心半徑。帶阻濾波就是阻止特定范圍頻率的信號通過，而允許其他范圍的信號通過。在

2、徑向對稱的濾波器的情況下，這種濾波器可以由一個剖面完全確定。在點(uO,vO)處的一個圓形鄰域上產生的帶阻濾波器是用高通濾波器的平移實現(xiàn)的。四、算法設計(含程序設計流程圖)獲取圖像的頻域圖像，函數(shù)觀察干擾噪聲的頻率成分，找出所有的干擾分量，由于干擾分量在大致在一個圓上，所以可以用一個頻譜特性為一個圓環(huán)的帶阻濾波器來濾除干擾分量。濾波完成后將頻域圖像實施逆變換，即可顯示處理后的圖像。五.實驗結果及分析（需要給出原始圖像和處理后的圖像）,*、，幾,flx:;後題冷.-w-5;:cwrws:;:?-.:;s.2o:v.:;.;?;v-H.?:=:&-:-x;4:c:八oc;c*cov歹0wcv*g;

3、Q4Q.rwc-ok:tg浜bwHJl靈:lvicc:;$sd;:.:v皺s-HKdn-Q2:c“3Hogsvxou八介s“s9?-.-關二=:?,0匕&66空?;仆-Z:;0苫卞03亦:eft9$!靈：*。：3乂暑!1來兌0洽介9;:;.歪一0?.0:0左,仝ccnmo芒xl-4:仝;i/.;c.My-.H-cx.2a3mov32遼亠(pcaoxHWC/e仝:.:-:-s;v-i;:os.fiM審至氓!Kv.k:u-8;u.fty亞淇Ag:s9譽.cwb;:.亠;:-.0:55空空出出亠歩疋ro:g溟比法.5;.訂0療3麥se2xvHvg艾京oomsx*Q背占X仝;:VV:O:EV2:久ui

4、wviwuxlanoiwowci.:7:HF1V“矚空:nf-K.!vf4I9.-fi;v-ff:s0:s3-x:x:tv?rNm.k2R52?:ro;YC3;wx:v;:x:iMGBM?-a:3:?i-65i:da-a:&r-;u:6:r:;R:;?;念cgfi:3:Jo;s?$:c復Iro*世n:K*sH.u;5:&rci.3:sv5wa8:CG婁!皿：卞*上.JHA4qzAy*pv*wA攻、V丿，MA*-tn3*53c.$3CW0vw:p*Nr&x色啜&總其06錢oxigigojxaexaov渋営八。荷:gxxaoflQ4waGxagKi;p:D&gac0sKcy*y:6屈oKv!rtv

5、Av:o2KHwwxv5x:nc:0.:G-3;_;H0;-r4*tc7:.-:nc-:/;d:k:crX八c$繼Mo;:H5icrji祜y/:$:Atr-:g用324vMC*5QPP-n0161:八x:o:v:c:2:ivo:!-:一後謖*.vxo:oi囁圣x.vpHtwo奏QH瑟主備the溼4s5OHLbT-!A*$;xHsio、：A工RawtM姜,垓-ii?2;-5浴*rt:;沖6上心.xsbulM裁4汛點2:tQig七vi:9QmYg1:帶阻濾波實驗原圖帶阻濾波結果帶阻濾波原圖的頻譜圖分析：此實驗利用阻帶濾波器阻止特定范圍內的頻率通過，而其他頻率則被濾掉。當設置截止久=175,而帶寬w

6、=150,利用此阻帶進行濾波消除周期性噪聲，且細小的細節(jié)和紋理也被修復。但是效果不是很好，有可能是混有其他的干擾，或是需要用其他的濾波器過濾噪聲。附：程序源代碼clear;data=xmread(fjpg)；numset=175;w=150;data=double(data);row=size(data,1);col=size(data,2);f=zeros(row,col);fori=1:rowforj=1:colf(i,j)=(-DA(i+j)*data(i,j);endendf=fft2(f);R=real(f);I=xmag(f);fl=sqrt(R人2+1.2);f1=(f1-min(min(fl)/(max(max(fl)-min(min(f1)*225;iimshow(fl);imwrite(fl,fbmp*);hh=zeros(row,col);fori=l:rowforj=1:coltemp=sqrt(i-row/2)2+(j-col/2)A2);hh(izj)=1/(1+(temp*w/(temp2-numsetA2);endendout=f*hh;fr=real(ifft2(out);fori=1:rowf