數(shù)字圖像處理-加噪與濾波(含MATLAB代碼).docx

《數(shù)字圖像處理-加噪與濾波(含MATLAB代碼).docx》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在工程資料-天天文庫。

1、數(shù)字圖像處理實驗三15生醫(yī)一、實驗內(nèi)容產(chǎn)生教材104頁題圖4.16所示的灰度圖像（白為255，黑為0），分別加入高斯白噪聲和椒鹽噪聲，再分別進(jìn)行3x3的平均濾波和中值濾波，顯示原圖像、加噪圖像和濾波結(jié)果圖像，并比較四種濾波結(jié)果。二、運行環(huán)境MATLABR2014a三、運行結(jié)果及分析①由實驗結(jié)果可以看出：平均濾波比中指濾波的效果差一點，平均濾波對原圖像的邊緣保留得不好，交界處較為模糊。同時，均值濾波器是用鄰域內(nèi)各像素的灰度平均值代替該像素原來的灰度值，實現(xiàn)圖像的平滑，它對高斯噪聲的抑制效果好。中值濾波器是用像素鄰域內(nèi)的中間值代替該像素值。它對脈沖干擾及椒鹽噪

2、聲的抑制效果好，在抑制隨機(jī)噪聲的同時能有效保護(hù)邊緣少受模糊。②原因分析：平均濾波屬于局部平均的線性濾波，所用的線性濾波器的工作原理可以比喻為用水沖澆桌面上的污點，沖洗的結(jié)果是污點并沒有消失，只是被淡化，也有可能是整個桌面都被污點所影響。而中值濾波的基本思想是將污點直接去除。在含噪圖像中,噪聲往往以孤立點的形式出現(xiàn)，尤其是干擾脈沖和椒鹽類噪聲。這些干擾或噪聲所占的像素數(shù)很少，而圖像則是由像素數(shù)目較多、面積較大的小塊組成?；谶@一事實，可采用中值濾波方法來消除噪聲，同時避免采用平均濾波時帶來的目標(biāo)物邊緣模糊的問題。③放大右邊結(jié)果部分的截圖：四、心得體會通過MA

3、TLAB編程更加直觀地了解了平均濾波與中值濾波的效果差異，使課本理論得到驗證，另一方面也對噪聲模型及其適用的濾波模板有了一定的認(rèn)識。同時，對MATLAB產(chǎn)生噪聲的函數(shù)調(diào)用及實現(xiàn)像素點灰度值的替換機(jī)理也有所掌握，比如后邊附的程序中會提到的“%”標(biāo)注的思考。五、具體程序%創(chuàng)建原圖像f=zeros(256,256);unit=255*ones(210,7);fori=0:8f(24:233,(29+24*i):(29+24*i+6))=unit;%坐標(biāo)從1開始（不是0），因此白色區(qū)域的縱向起始像素點為(256-210)/2+1=24，%止點為256-24+1=2

4、33，同理橫向起始像素點為(256-7x9-17x8)/2=28.5，取28%每個白條寬為7，間隔為7+17=24，因此下一個白條始于28+24+1，以此類推endf1=imnoise(f,'gaussian',0,0.1);%均值為零方差為0.1的高斯噪聲f2=imnoise(f,'salt&pepper',0.1);%噪聲密度為0.1的椒鹽噪聲fori=2:1:255forj=2:1:255%3x3的8領(lǐng)域平均濾波f1_mean(i,j)=(f1(i-1,j-1)+f1(i-1,j)+f1(i-1,j+1)+f1(i,j-1)+f1(i,j+1)+f1

5、(i+1,j-1)+f1(i+1,j)+f1(i+1,j+1))/8;f2_mean(i,j)=(f2(i-1,j-1)+f2(i-1,j)+f2(i-1,j+1)+f2(i,j-1)+f2(i,j+1)+f2(i+1,j-1)+f2(i+1,j)+f2(i+1,j+1))/8;%3x3的方形中值濾波a1=[f1(i-1,j-1),f1(i-1,j),f1(i-1,j+1),f1(i,j-1),f1(i,j),f1(i,j+1),f1(i+1,j-1),f1(i+1,j),f1(i+1,j+1)];a2=[f2(i-1,j-1),f2(i-1,j),f2(

6、i-1,j+1),f2(i,j-1),f2(i,j),f2(i,j+1),f2(i+1,j-1),f2(i+1,j),f2(i+1,j+1)];b1=sort(a1);b2=sort(a2);f1_mid(i,j)=b1(5);f2_mid(i,j)=b2(5);endendsubplot(2,4,1);imshow(f);title('原圖像');subplot(2,4,2);imshow(f1);title('加入高斯噪聲');subplot(2,4,3);imshow(f1_mean);title('進(jìn)行平均濾波');subplot(2,4,4);i

7、mshow(f1_mid);title('進(jìn)行中值濾波');subplot(2,4,5);imshow(f);title('原圖像');subplot(2,4,6);imshow(f2);title('加入椒鹽噪聲');subplot(2,4,7);imshow(f2_mean);title('進(jìn)行平均濾波');subplot(2,4,8);imshow(f1_mid);title('進(jìn)行中值濾波');