形态学重建 matlab,图像形态学处理——灰度开操作和闭操作（Matlab）

clc;

clear all;

close all;

%灰度膨胀、腐蚀、开操作、闭操作测试图像

I=im2double(imread('D:\Gray Files\9-37.tif'));

%========================灰度开操作和闭操作================================

%半径为3的磁盘(disk)结构元素

r=3;

B=ones(2*r-1,2*r-1);

%将距离大于40的点置零

for i=1:2*r-1

for j=1:2*r-1

d=sqrt((i-r)^2+(j-r)^2);

if d>r

B(i,j)=0;

end

end

end

n=2*r-1;

ind=find(B==0);

n_l=floor(n/2);

I_pad=padarray(I,[n_l,n_l],'symmetric');

[M,N]=size(I);

%-------------------------------灰度开操作---------------------------------

J_Opening=zeros(M,N);

%腐蚀操作

J_Erosion=zeros(M,N);

for i=1:M

for j=1:N

%获得图像子块区域

Block=I_pad(i:i+2*n_l,j:j+2*n_l);

C=Block.*B;

%删除0值，保留4连通数值

C=C(:);

C(ind)=[];

J_Erosion(i,j)=min(C);

end

end

%对腐蚀图像进行扩展

J_Erosion_pad=padarray(J_Erosion,[n_l,n_l],'symmetric');

%膨胀图像

for i=1:M

for j=1:N

%获得图像子块区域

Block=J_Erosion_pad(i:i+2*n_l,j:j+2*n_l);

C=Block.*B;

%删除0值，保留4连通数值

C=C(:);

C(ind)=[];

J_Opening(i,j)=max(C);

end

end

%-------------------------------灰度闭操作---------------------------------

J_Closing=zeros(M,N);

%膨胀图像

J_Dilation=zeros(M,N);

for i=1:M

for j=1:N

%获得图像子块区域

Block=I_pad(i:i+2*n_l,j:j+2*n_l);

C=Block.*B;

%删除0值，保留4连通数值

C=C(:);

C(ind)=[];

J_Dilation(i,j)=max(C);

end

end

%对膨胀图像进行扩展

J_Dilation_pad=padarray(J_Dilation,[n_l,n_l],'symmetric');

%腐蚀操作

for i=1:M

for j=1:N

%获得图像子块区域

Block=J_Dilation_pad(i:i+2*n_l,j:j+2*n_l);

C=Block.*B;

%删除0值，保留4连通数值

C=C(:);

C(ind)=[];

J_Closing(i,j)=min(C);

end

end

subplot(131)

imshow(I)

subplot(132)

imshow(J_Opening)

subplot(133)

imshow(J_Closing)