数字图像处理腐蚀与膨胀小程序

%%%%%%%%%%%%%%%%下面为用于处理的样本‘圆’%%%%%%%%%%%%%%%%%%%%%%%%%
R=8;%%半径R
S=zeros(2*R);
for x1=-R:R-1;
for y1=-R:R-1;
if round(sqrt(x1^2+y1^2))<=R%%圆心在原点的圆的表达函数
%%对圆心在原点的圆进行平移，使x,y都为非负数
x=R+x1+1;
y=R+y1+1;
S(x,y)=1;
end
end
end
figure(1),imshow(S); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%下面为要处理的图形%%%%%%%%%%%%%%%%%%%%%%%%%
W=256;L=512; %%定义边界
X=zeros(W,L);
for i=1:128
for j=1:256
if i<=15|i>=114
if j<103|j>153
X(i+W/4,j+L/4)=1;
end
else if (i>15&i<=56)|(i>=72&i<114)
if (j>50&j<103)|j>153
X(i+W/4,j+L/4)=1;
end
else if i>56&i<72
if j>50&j<206
X(i+W/4,j+L/4)=1;
end
end
end
end
end
end
figure(2),imshow(X); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%下面用S对X进行腐蚀操
作%%%%%%%%%%%%%%%%%%%%%%%%
X1=zeros(W,L);%X1用于存放处理后的图像
M=zeros(2*R);
for m=0:W-2*R
for n=0:L-2*R
%%下面用S对X区域进行‘与’操作%%
for i=1:2*R
for j=1:2*R
M(i,j)=X(i+m,j+n)&S(i,j);
end
end
if M==S %%判断S是否处于X中，等则在S中
X1(m+R,n+R)=1;
end
M=zeros(2*R);
end
end
figure(3),imshow(X1); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%下面用S对X1进行膨胀操作%%%%%%%%%%%%%%%%%%%%%%%%%%
X2=zeros(W,L);%X2用于存放处理后的图像
cnt=0;
for m=0:W-2*R
for n=0:L-2*R
%%下面用S对X区域进行‘与’操作%%
for i=1:2*R
for j=1:2*R
if X1(i+m,j+n)&S(i,j)==1 %%判断S是否与X相交
cnt=1;
break;
end
end
if cnt==1 %%判断是否是由于cnt=1而跳出循环
break;
end
end
if cnt==1 %%判断是否是由于cnt=1而跳出循环
X2(m+R,n+R)=1;
cnt=0; %%将cnt的值置0
end
end
end
figure(4),imshow(X2); ( S )
( X )
( X1 )
( X2 )。