图像形态学处理——灰度平滑（Matlab）

clc;

clear all;

close all;

%灰度平滑测试图像

I=im2double(imread('D:Gray Files9-38.tif'));

%==============================灰度平滑====================================

%半径为3的磁盘（disk）结构元素

% r=3;

% 半径为5的磁盘（disk）结构元素

r=5;

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_Smoothing=zeros(M,N);

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

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_Opening_pad=padarray(J_Opening,[n_l,n_l],'symmetric');

%膨胀图像

J_Dilation=zeros(M,N);

for i=1:M

for j=1:N

%获得图像子块区域

Block=J_Opening_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_Smoothing(i,j)=min(C);

end

end

subplot(121)

imshow(I)

subplot(122)

imshow(J_Smoothing)