第9章 形态学处理
腐蚀操作 imerode
- 通过腐蚀操作(imerode)去掉图像的某些部分:
%% 第9章 形态学处理
%% imerode 腐蚀操作
clc
clear
A1=imread('./1.jpg');
subplot(221),imshow(A1);title('原始图像');
%strel函数的功能是运用各种形状和大小构造结构元素
se1=strel('disk',5);%这是创建一个半径为5的平坦圆盘结构元素
A2=imerode(A1,se1);
subplot(222),imshow(A2);title('使用disk(5,5)腐蚀后的图像');
se2=strel('rectangle',[3,3]);
A3=imerode(A1,se2);
subplot(223),imshow(A3),title('使用rectangle[5,5]腐蚀操作后');
se3=strel('disk',10);%这是创建一个半径为5的平坦圆盘结构元素
A3=imerode(A1,se3);
subplot(224),imshow(A3);title('使用disk(10,10)腐蚀后的图像');
figure
se4=strel('line',10,0);
A4=imerode(A1,se4);
subplot(221),imshow(A4);title('使用line(10,0)腐蚀后的图像');
se5=strel('line',10,45);
A5=imerode(A1,se5);
subplot(222),imshow(A5);title('使用line(10,45)腐蚀后的图像');
se6=strel('line',10,90);
A6=imerode(A1,se6);
subplot(223),imshow(A6);title('使用line(10,90)腐蚀后的图像');
通过观察,发现结构元的尺寸越大,腐蚀的效果越明显;
当用45°的线进行操作时,留下了该方向的细线和竖直的较宽的细线;
当用90°的结构元进行操作时,对中间较粗的细线几乎没有进行任何腐蚀。
规律似乎是这样的:对图像中与结构元相似的形状进行腐蚀效果不明显,比结构元规格大的形状腐蚀完仍会有保留。
膨胀操作 imdilate
- 膨胀操作用于桥接裂缝
%% 第9章 形态学处理
%% imdilate 膨胀处理
clc
clear
A1=imread('./2.jpg');
subplot(221);imshow(A1);title('原始图像');
info=imfinfo('./2.jpg');
%se1=strel('disk',1);
se1=strel('line',10,0);
A2=imdilate(A1,se1);
%subplot(222);imshow(A2);title('disk(1)-1次膨胀');
subplot(222);imshow(A2);title('line(10,0)-1次膨胀');
A3=imdilate(A2,se1);
%subplot(223);imshow(A3);title('disk(1)-2次膨胀');
subplot(223);imshow(A3);title('line(10,0)-2次膨胀');
A4=imdilate(A3,se1);
%subplot(224);imshow(A4);title('disk(1)-13膨胀');
subplot(224);imshow(A4);title('line(10,0)-3次膨胀');
使用disk(1)进行膨胀操作,多次操作,效果加深。
使用line(10,0)进行多次膨胀,随着膨胀次数的增多,对垂直方向的效果比水平方向的效果明显。
规律:膨胀操作时只要是有背景与结构元的反射有一点重叠,即为新集合的元素,因此用横向的结构元膨胀时,对横向目标不如对纵向的目标效果明显。腐蚀是以结构元的中心遍历图像的每一个像素进行访问,如果结构元完全包含在图像中,则该中心位置为腐蚀后的新集合中的元素;
膨胀是以结构元的中心遍历图像的每一个像素进行访问,如果结构元与图像的交集不为空,则该结构元所占区域为新结合中的元素;
开运算 &闭操作
开运算一般会平滑物体的轮廓,断开较窄的狭颈并消除细的突出物。
闭操作同样也会平滑物体的轮廓,但是与开操作相反,它通常会弥合较窄的间断和细长的沟壑,消除小的空洞,填补轮廓线中的断裂。
%% 第9章 形态学处理
%% 开运算
clc
clear
%% 开运算
A1=imread('./5.jpg');
subplot(221),imshow(A1);title('原始图像');
se1=strel('disk',20);%这是创建一个半径为5的平坦圆盘结构元素
A2=imopen(A1,se1);
subplot(222),imshow(A2);title('开运算');
%%自带闭运算
A3=imclose(A1,se1);
subplot(223),imshow(A3);title('闭运算');
%% 第2副图像
A1=imread('./5.jpg');
%subplot(234),imshow(A1);title('原始图像');
se1=strel('disk',20);%这是创建一个半径为5的平坦圆盘结构元素
A2=imclose(A1,se1);
%subplot(235),imshow(A2);title('闭运算');
%%自带闭运算
A3=imopen(A1,se1);
%subplot(236),imshow(A3);title('开运算');
开运算:将较细的线断开,消除突出物
闭运算:弥合较窄的沟壑和细长的线段
对比上图中两幅原始图像(反相)的开运算和闭运算,可以发现,对背景(黑色区域)的开运算相当于对目标(白色区域)的闭运算;对目标(白色区域)的开运算相当于对背景(黑色区域)的闭运算。
击中或击不中变换
击中或击不中变换是形状检测的一个基本工具。假设集合B由B1和B2组合成,B1对图像A腐蚀,B2对图像的补集腐蚀,两者得到的结果相减即为B在图像A中的匹配。
可以参考知乎中的解释:https://www.zhihu.com/question/23669012
%% 击中击不中变换
clc
clear
A=im2bw(rgb2gray(imread('./AT9A3130.JPG')));
subplot(161);imshow(A);title('原始图像');
B=~A;
se1= strel('disk',3);
se2= strel('disk',5);
A11=imerode(A,se1);
A1=mat2gray(A11);
subplot(162);imshow(A1);title('se1腐蚀');
A22=imerode(B,se2);
A2=mat2gray(A22);
subplot(163);imshow(A2);title('se2腐蚀');
subplot(164);imshow(mat2gray(A-A1));title('原图减去腐蚀1的结果');
subplot(165);imshow(mat2gray(A1-A2));title('击中击不中结果');
subplot(166);imshow(mat2gray(~(abs(A1-A2))));title('上图结果取绝对值再取反');
上图中原图减去se1腐蚀结果的图像,可以用来做边界提取。
空洞修复
%% 空洞修复
A=im2bw(imread('./AT9A3130.JPG'));
subplot(221);imshow(A);title('原始图像');
B=imfill(A,'holes');
subplot(222);imshow(B);title('填充后图像');
%% 注解
BW4 = im2bw(imread('kongdong.jpg'));
subplot(223);imshow(BW4);title('原始图像');
BW5 = imfill(BW4,'holes');
subplot(224);imshow(BW5);title('填充后');
如何控制空洞的大小呢???
- 未完待续…