MATLAB 二值图像腐蚀

clear
clc
I=imread('fushi.png');
[L,R,D]=size(I);
IR=I(:,:,1);
IG=I(:,:,2);
IB=I(:,:,3);
Y=0.299.*IR+0.587.*IG+0.114.*IB;
num=110;
for i=1:L
    for j=1:R
        if Y(i,j)>num
            F(i,j)=1;
        else
            F(i,j)=0;
        end
    end
end
figure(1)
subplot(1,2,1);
imshow(F);
title('原图二值化（１１０）')

[L,R]=size(F);
s=[1 1 1 1 1
    1 1 1 1 1
    1 1 1 1 1
    1 1 1 1 1
    1 1 1 1 1]    %形态学为5*5个1
[s_L,s_R]=size(s);
num1=s_L-1;
num2=s_R-1;
for i=1:L
    for j=1:R
        if F(i,j) == 1 & i <= L-num1 & j <= R-num2
            E1 = (s(1,1)&F(i,j)) & (s(1,2)&F(i,j+1)) & (s(1,3)&F(i,j+2)) & (s(1,4)&F(i,j+3)) & (s(1,5)&F(i,j+4));
            E2 = (s(2,1)&F(i+1,j)) & (s(2,2)&F(i+1,j+1)) & (s(2,3)&F(i+1,j+2)) & (s(2,4)&F(i+1,j+3)) & (s(2,5)&F(i+1,j+4));
            E3 = (s(3,1)&F(i+2,j)) & (s(3,2)&F(i+2,j+1)) & (s(3,3)&F(i+2,j+2)) & (s(3,4)&F(i+2,j+3)) & (s(3,5)&F(i+2,j+4));
            E4 = (s(4,1)&F(i+3,j)) & (s(4,2)&F(i+3,j+1)) & (s(4,3)&F(i+3,j+2)) & (s(4,4)&F(i+3,j+3)) & (s(4,5)&F(i+3,j+4));
            E5 = (s(5,1)&F(i+4,j)) & (s(5,2)&F(i+4,j+1)) & (s(5,3)&F(i+4,j+2)) & (s(5,4)&F(i+4,j+3)) & (s(5,5)&F(i+4,j+4));
            E(i,j) =E1 & E2 & E3 & E4 & E5;
        else
            E(i,j)=0;
        end
    end
end
figure(1)
subplot(1,2,2);
imshow(E);
title('腐蚀')

效果：细线被腐蚀了

可以改变s的结构元素或改变s的大小，看看效果
膨胀也是差不多一样的道理，有空再写写