Matlab圖像處理(四)——常用濾波的實現

2 [height, width]=size(A_salt_pepper);   %輸入圖像是p×q的,且p>n,q>n  
3 A_salt_pepper2=zeros(height+n-1,width+n-1); %存放延拓後的圖片
4 for i=1+(n-1)/2:1+(n-1)/2+height-1
5     for j=1+(n-1)/2:1+(n-1)/2+width-1
6         A_salt_pepper2(i,j)=A_salt_pepper(i-(n-1)/2,j-(n-1)/2);
7     end
8 end
9 A_salt_pepper3=A_salt_pepper2;   %灰度值更改後會對後續的濾波產生影響，因此需要一個不變的圖片
10a(1:n,1:n)=1;
11for i=1:height  
12    for j=1:width  
13        c=A_salt_pepper3(i:i+(n-1),j:j+(n-1)).*a; %取出x1中從(i,j)開始的n行n列元素與模板相乘  
14        s=sum(sum(c));                 %求c矩陣中各元素之和  
15        A_salt_pepper2(i+(n-1)/2,j+(n-1)/2)=s/(n*n); %將與模板運算後的各元素的均值賦給模板中心位置的元素  
16    end  
17end  
18A_salt_pepper2=uint8(A_salt_pepper2((n-1)/2+1:height+(n-1)/2,(n-1)/2+1:width+(n-1)/2));%通過計算，此處只是一個矩陣，用uint8將矩陣變成圖像
19imshow(A_salt_pepper2);
20title('自己編寫均值濾波圖')