图像形态学处理——灰度开操作和闭操作（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)