www.pudn.com > Gauss_Process.rar > num23.m

clear; 
A=imread('num23','bmp'); 
subplot(2,2,1); 
B=rgb2gray(A); 
subimage(B); 
title('处理前的图'); 
%求高斯阵 
G33=zeros(3,3);%3*3的 
G77=zeros(7,7);%7*7的 
G91111=zeros(11,11);%11*11的 
%求G33 
G(1,1)=exp(-(1+1)/2); 
for j=-1:1 
    for i=-1:1 
        G33(j+2,i+2)=exp(-(j^2+i^2)/2)/G(1,1); 
    end 
end 
%求G77 
G(1,1)=exp(-(9+9)/2); 
for j=-3:3 
    for i=-3:3 
        G77(j+4,i+4)=exp(-(j^2+i^2)/2)/G(1,1); 
    end 
end 
%求G1111 
G(1,1)=exp(-(25+25)/2); 
for j=-5:5 
    for i=-5:5 
        G1111(j+6,i+6)=exp(-(j^2+i^2)/2)/G(1,1); 
    end 
end 
B=double(B); 
C=B; 
xsize=size(B); 
%用G33处理 
gg=sum(G33); 
gg=sum(gg); 
for k=2:(xsize(1)-1) 
    for j=2:(xsize(2)-1) 
        t=sum(B(k-1:k+1,j-1:j+1).*G33);         
        C(k,j)=sum(t)/gg; 
    end 
end 
C=uint8(C); 
subplot(2,2,2); 
subimage(C); 
title('用3*3高斯矩阵处理后的图'); 
%用G77处理 
gg=sum(G77); 
gg=sum(gg); 
for k=4:(xsize(1)-3) 
    for j=4:(xsize(2)-3) 
        t=sum(B(k-3:k+3,j-3:j+3).*G77);         
        C(k,j)=sum(t)/gg; 
    end 
end 
C=uint8(C); 
subplot(2,2,3); 
subimage(C); 
title('用7*7高斯矩阵处理后的图'); 
%用G1111处理 
gg=sum(G1111); 
gg=sum(gg); 
for k=6:(xsize(1)-5) 
    for j=6:(xsize(2)-5) 
        t=sum(B(k-5:k+5,j-5:j+5).*G1111);         
        C(k,j)=sum(t)/gg; 
    end 
end 
C=uint8(C); 
subplot(2,2,4); 
subimage(C); 
title('用11*11高斯矩阵处理后的图');