[数字图像处理]3*3模板灰度图像空间域变换

function [  ] = picSmooth0(A ,s)
%picSmooth0 此处显示有关此函数的摘要
%picSmooth0(A,s)  A为图像矩阵    s为3*3变换矩阵  
%   此处显示详细说明
[M,N,a]=size(A);
if a == 3
    B=rgb2gray(A);
else
    B=A;
end
B=double(B);    %防止溢出
P=zeros(M,N);
k=0;
sum=0;
[m,n]=size(s);
if( m~=3 || n~=3)
       fprintf('请输入正确的3*3空间域平滑模板');
else
        for i=1:3
            for j=1:3
                sum=sum+s(i,j);
            end
        end
        
        for i=1:M
            for j = 1:N
                k=0;
                if((0<i-1&&i-1<=M) && (0<j-1&&j-1<=N))
                    k=k+s(1,1)*B(i-1,j-1);
                end
                if(0<i&&i<=M && 0<j-1&&j-1<=N)
                    k=k+s(2,1)*B(i,j-1);
                end
                if(0<i+1&&i+1<=M && 0<j-1&&j-1<=N)
                    k=k+s(3,1)*B(i+1,j-1);
                end
                    
                if(0<i-1&&i-1<=M && 0<j&&j<=N)
                    k=k+s(1,2)*B(i-1,j);
                end
                if(0<i&&i<=M && 0<j&&j<=N)
                    k=k+s(2,2)*B(i,j);
                end
                if(0<i+1&&i+1<=M && 0<j&&j<=N)
                    k=k+s(3,2)*B(i+1,j);
                end
                    
                if(0<i-1&&i-1<=M && 0<j+1&&j+1<=N)
                    k=k+s(1,3)*B(i-1,j+1);
                end
                if(0<i&&i<=M && 0<j+1&&j+1<=N)
                    k=k+s(2,3)*B(i,j+1);
                end
                if(0<i+1&&i+1<=M && 0<j+1&&j+1<=N)
                    k=k+s(3,3)*B(i+1,j+1);
                end
                P(i,j)=k;
            end
        end            
end
P=P/sum;
P=uint8(P);
B=uint8(B);

figure;
imshow(B);
figure;
imshow(P);

end

 转载请注明来源