fssss.m

自己编的腐蚀源程序

clear all;
clc;
yuantu=imread('kids.jpg');
figure,subplot(2,2,1),imshow(yuantu);
[height,width]=size(yuantu);
s=ones(5,5);
a=100;
b=130;
[sh,sw]=size(s);
sum=length(find(s==1));
sorigin=[3 3];

imresult=zeros(height,width);
for i=3:1:height-2
    for j=3:1:width-2
        count=1;
        for k=1:1:sh
            for w=1:1:sw
                if yuantu(i-3+k,j-3+w)>=a & yuantu(i-3+k,j-3+w)<=b & s(k,w)==1
                    if count>=25
                    imresult(i,j)=yuantu(i,j); 
                    else
                        count=count+1;
                    end
                else
                    break;
                end
            end
        end
    end
end

subplot(2,2,2),imshow(imresult);