灰度图形态学运算知识点总结 gray-scale morphology

最新推荐文章于 2022-03-31 21:41:41 发布

Full_Speed_Turbo

最新推荐文章于 2022-03-31 21:41:41 发布

阅读量5.5k

点赞数 1

分类专栏：图像处理数字图像处理文章标签：灰度图形态学

本文链接：https://blog.csdn.net/full_speed_turbo/article/details/70048715

版权

图像处理同时被 2 个专栏收录

26 篇文章 6 订阅

订阅专栏

数字图像处理

22 篇文章 7 订阅

订阅专栏

这里写图片描述

一般, 结构元素B是flat.

1. 腐蚀

f以正无穷大扩边, 选取B覆盖f区域中的最小值.

2. 膨胀

f以负无穷大扩边, 选取B覆盖f区域中的最大值.

3. 开运算

先腐蚀后膨胀.
开运算会抑制小size的白点(高灰度对象).

4. 闭运算

先膨胀后腐蚀.
闭运算会填充小size的黑洞(低灰度对象).

5. 腐蚀,膨胀,开,闭的一维仿真

四种运算的仿真代码:

%% 灰度图形态学运算
clc;close all;clear;
x=1:600;
y=2*sin(0.01*x)+sin(0.02*x)+sin(0.04*x);
SELen=51; % odd number
SE=ones(1,SELen);
%% 腐蚀运算 
y_erosion=nan(size(y));
for i=1:length(y)
    h=max(1,i-floor(SELen/2));
    t=min(length(y),i+floor(SELen/2));
    y_erosion(i)=min(y(h:t));
end
%% 膨胀
y_dilation=nan(size(y));
for i=1:length(y)
    h=max(1,i-floor(SELen/2));
    t=min(length(y),i+floor(SELen/2));
    y_dilation(i)=max(y(h:t));
end
%% 开运算
y_open=y_erosion; % 先腐蚀后膨胀
for i=1:length(y_erosion)
    h=max(1,i-floor(SELen/2));
    t=min(length(y_erosion),i+floor(SELen/2));
    y_open(i)=max(y_erosion(h:t));
end
%% 闭运算
y_close=y_dilation; % 先膨胀后腐蚀
for i=1:length(y_dilation)
    h=max(1,i-floor(SELen/2));
    t=min(length(y_dilation),i+floor(SELen/2));
    y_close(i)=min(y_dilation(h:t));
end
figure;plot(x,y,'r-',x,y_erosion,'g-',x,y_dilation,'b-');
xlabel('点序号');ylabel('灰度值');legend('原曲线','腐蚀','膨胀');
figure;plot(x,y_open,'r-',x,y_close,'b-');
xlabel('点序号');ylabel('灰度值');
legend('开','闭');

输出:
这里写图片描述

这里写图片描述
注意: 闭运算曲线开头处, 是一段水平线, 没有和原曲线重合. 和DIP书中不一致.

使用Matlab的内置函数验证一下:

%% 内置灰度形态运算
clc;close all;clear;
x=1:600;
y=2*sin(0.01*x)+sin(0.02*x)+sin(0.04*x);

SELen=51; % odd number
SE = strel('line', 50, 0);
%% imerode imdilate
y_imerode=imerode(y,SE);
y_imdilate=imdilate(y,SE);
%% imclose 闭运算
y_imclose = imclose(y,SE);
%% imopen 开运算
y_imopen=imopen(y,SE);
%% plot
figure;plot(x,y,'r-',x,y_imerode,'g-',x,y_imdilate,'b-');
xlabel('点序号');ylabel('灰度值');legend('原曲线','imerode','imdilate');
figure;plot(x,y_imopen,'r-',x,y_imclose,'b-');
legend('imopen','imclose');