imageerode.m

上个形态学的基本方法的程序的

%--------------------------------------------------------------------------
%  image erosion function
%
%  result: 腐蚀图像中灰度值为黑色的像素
%
%  principle:分析原始图像当前点与结构元素结构相同的子图像，
%            当结构元素中某一像素为0，而子图像相应位置处的灰度值不为0，
%            则将目标图像中当前点的灰度值置为1，否则置为0      
% 
%  input parameters:
%    imoriginal--原始图像
%    s--结构元素
%    sorigin--结构元素原点
%--------------------------------------------------------------------------
function imresult = imageerode(imoriginal,s,sorigin);

[imh,imw] = size(imoriginal);   % 读取原始图像和结构元素的尺寸
[sh,sw] = size(s);

imresult = ones(imh,imw);   % 初始化目标图像

% 腐蚀
for i = sorigin(1,1):1:imh-(sh-sorigin(1,1))
    for j = sorigin(1,2):1:imw-(sw-sorigin(1,2))
        flag = 0;
        
        for si = 1:1:sh;
            for sj = 1:1:sw
                if s(si,sj)==0 & imoriginal(i-sorigin(1,1)+si,j-sorigin(1,2)+sj)~=0
                    flag = 1;
                    break;
                end
            end
        end

        if flag == 0
            imresult(i,j) = 0;
        end
    end
end