圖像處理中消除噪聲的方法研究——學(xué)士學(xué)位畢業(yè)論文(編輯修改稿)
2024-07-19 20:50 本頁面
【文章內(nèi)容簡介】 利用Matlab軟件對含噪圖像的去噪算法進(jìn)行仿真，將應(yīng)用中值濾波法、均值濾波法和小波變換法對含有高斯噪聲和椒鹽噪聲圖像的去噪效果進(jìn)行比較，從而得到相應(yīng)結(jié)論。 中值濾波的仿真 本節(jié)選用中值濾波法對含有高斯噪聲和椒鹽噪聲的圖像進(jìn)行去噪，并用Matlab軟件仿真。1）給圖像加入均值為0，分別選擇33模板、55模板和77模板進(jìn)行去噪Matlab部分代碼：j=imnoise(I,39。gaussian39。,0,)。x=j(:,:,1)。subplot(221)。imshow(x)。title(39。高斯噪聲圖片39。)。k1=medfilt2(x,[3 3])。k2=medfilt2(x,[5 5])。k3=medfilt2(x,[7 7])。仿真結(jié)果如圖1所示:圖1 中值濾波法對高斯噪聲去噪的仿真結(jié)果2），分別選擇33模板、55模板和77模板進(jìn)行去噪。Matlab部分代碼：i=imread(39。39。)。j=imnoise(I,39。salt amp。 pepper39。,)。x=j(:,:,1)。subplot(221)。imshow(x)。title(39。椒鹽噪聲圖片39。)。k1=medfilt2(x,[3 3])。k2=medfilt2(x,[5 5])。k3=medfilt2(x,[7 7])。仿真結(jié)果如圖2所示:圖2 中值濾波法對椒鹽噪聲去噪的仿真結(jié)果從仿真結(jié)果可以看出：對圖像加入椒鹽噪聲后，應(yīng)用中值濾波，如圖2所示，噪聲的斑點(diǎn)幾乎全部被濾去，它對濾除圖像的椒鹽噪聲非常有效。而對于高斯噪聲來說，如圖1所示，雖然也有一些去噪效果，但效果不佳。由此可知，中值濾波法運(yùn)算簡單，易于實(shí)現(xiàn)，而且能較好地保護(hù)邊界，但有時(shí)會失掉圖像中的細(xì)線和小塊區(qū)域。并且采用窗口的大小對濾波效果影響很大，窗口越大，圖像去噪效果越好，但代價(jià)是模糊的程度越大。 均值濾波的仿真本節(jié)選用均值濾波法對含有高斯噪聲和椒鹽噪聲的圖片進(jìn)行去噪，并用Matlab軟件仿真1）給圖像加入均值為0，選擇33模板去噪Matlab部分代碼：j=imnoise(x,39。gaussian39。,0,)。h=ones(3,3)。 h=h/9。k=conv2(j,h)。 仿真結(jié)果如圖3所示：圖3 均值濾波法對高斯噪聲去噪的仿真結(jié)果2），選擇33模板去噪Matlab部分代碼：j=imnoise(x,39。salt amp。 pepper39。,)。h=ones(3 3)。h=h/9。k=conv2(j,h)。仿真結(jié)果如圖4所示：圖4 均值濾波對椒鹽噪聲去噪的仿真結(jié)果從仿真結(jié)果可以看出：均值濾波實(shí)現(xiàn)起來很方便，適用于消除圖像中的顆粒噪聲，但需要指出這種方法既平滑了圖像信號，同時(shí)使圖像的細(xì)節(jié)部分變得模糊。由以上處理后的圖像可以看到：均值濾波法消弱了圖像的邊緣，使圖像變得有些模糊。如圖3所示，均值濾波對高斯噪聲的抑制是比較好的，但對椒鹽噪聲的抑制作用不好，如圖4所示，椒鹽噪聲仍然存在，只不過被削弱了而已。為了改善均值濾波細(xì)節(jié)對比度不好、區(qū)域邊界模糊的缺陷，常用門限法來抑制椒鹽噪聲和保護(hù)細(xì)小紋理，用加權(quán)法來改善圖像的邊界模糊，用選擇平均的自適應(yīng)技術(shù)來保持圖像的邊界。 基于小波變換法的仿真本節(jié)選用模糊小波變換法對含有高斯噪聲和椒鹽噪聲的圖像進(jìn)行去噪，并用Matlab軟件仿真。1）給圖像加入均值為0，用小波函數(shù)coif2對圖象進(jìn)行2層分解，選擇33模板去噪Matlab部分代碼：function y=zishiying(x)x11=medfilt2(x,[3 3])。 x12=double(x11)。[a,b]=size(x12)。[c,s]=wavedec2(x12,3,39。coif239。)。 n=[1,2,3]。 p4=*(sqrt(2*log(a*b)))。 size(detcoef2(39。h39。,c,s,1))。 size(detcoef2(39。v39。,c,s,1))。size(detcoef2(39。d39。,c,s,1))。p1(1)=detcoef2(39。h39。,c,s,1)。p2(1)=detcoef2(39。v39。,c,s,1)。p3(1)=detcoef2(39。d39。,c,s,1)。p1(2)=detcoef2(39。h39。,c,s,2)。 p2(2)=detcoef2(39。v39。,c,s,2)。p3(2)=detcoef2(39。d39。,c,s,2)。for i=1:1:2 p1(i)=1/((p1(i)p4)^2+1)。 if p1(i)=p4
電大資料相關(guān)推薦
鄂ICP備17016276號-1