Matlab基于腐蚀和膨胀的边缘检测
文/天神
http://blog.sina.com.cn/s/blog_4a93ccea0100d25q.html
一.课题背景:
形态学运算只针对二值图像(二进制图像),并依据数学形态学(Mathermatical Morphogy)集合论方法发展起来的图像处理方法,起源于岩相对岩石结构的定量描述工作,在数字图像处理和机器视觉领域中得到了广泛的应用,形成了一种独特的数字图像分析方法和理论。数学形态学是图像处理和模式识领域的新方法,其基本思想是:用具有一定形态的结构元素去量度和提取图像中的对应形状,以达到图像分析和识别的目的。优势有以下几点:有效滤除噪声,保留图像中原有信息,算法易于用并行处理方法有效实现(包括硬件实现),基于数学形态学的边缘信息提取处理优于基于微分运算的边缘提取算法,提取的边缘比较平滑,提取的图像骨架也比较连续,断点少。
二、课题相关原理:
形态学基本运算:
特殊领域运算形式——结构元素(Structure Element),在每个像素位置上与二值图像对应的区域进行特定的逻辑运算。运算结果是输出图像的相应像素。运算效果取决于结构元素大小内容以及逻辑运算性质。
常见形态学运算有腐蚀(Erosion)和膨胀(Dilation)两种。
集合论是数学形态学的基础。有集合、元素、子集、并集、补集、位移、映像(镜像对称)、差集等集合的基本概念。
对象和结构元素的3种关系:『对象X(Object)、结构元素B(Structure Element)』
B include in X 包含于 、B hit X 击中(不全包含) 、B miss X
平移、对称集:Bx=Uy{x+y}
腐蚀:一种消除边界点,使边界向内部收缩的过程。利用它可以消除小而且无意义的物体。B对X腐蚀所产生的二值图像E是满足以下条件的点(x,y)的集合:如果B的原点平移到点(x,y),那么B将完全包含于X中。
膨胀:将与物体接触的所有背景点合并到该物体中,使边界向外部扩张的过程。利用它可以填补物体中的空洞。B对X膨胀所产生的二值图像D是满足以下条件的点(x,y)的集合:如果B的原点平移到点(x,y),那么它与X的交集非空。
腐蚀和膨胀运算中存在对偶原理:X⊕B,它是所有满足以下条件的点X'的集合:在B中存在一点y,而且在X中存在一点x,使得x'=x+y。
基本运算:1.开运算(先腐蚀后膨胀的过程):利用它可以消除小物体,在纤细点处分离物体,平滑较大物体边界,但同时并不明显改变原来物体的面积。OPEN(X,B)
2.闭运算(先膨胀后腐蚀的过程):利用它可以填充物体内细小空洞,连接临近物体、平滑其边界,但同时并不明显改变原来物体的面积。CLOSE(X,B)
通常由于噪声的影响,图像在阈值化后所得到的边界通常都很不平滑,物体区域具有一些噪声孔,而背景区域上散布着一些小的噪声物体,连续的开和闭运算可以有效的改善这种情况,而有时,我们需要经过多次腐蚀之,后再加上相同次数的膨胀,才能产生比较好的处理效果。
另外两种是 3.击中,击不中变换HMT(模板严格匹配) 以及 4.边缘和骨架(Boundary and Skeleton)
三、腐蚀和膨胀的Matlab实现:
腐蚀:删除对象边界某些像素。
膨胀:给图像中的对象边界添加像素。
在操作中,输出图像中所有给定像素的状态都是通过对输入图像的相应像素及邻域使用一定的规则进行确定。在膨胀操作时,输出像素值是输入图像相应像素邻域内所有像素的最大值。在二进制图像中,如果任何像素值为1,那么对应的输出像素值为1;而在腐蚀操作中,输出像素值是输入图像相应像素邻域内所有像素的最小值。在二进制图像中,如果任何一个像素值为0,那么对应的输出像素值为0。
结构元素的原点定义在对输入图像感兴趣的位置。对于图像边缘的像素,由结构元素定义的邻域将会有一部分位于图像边界之外。为了有效处理边界像素,进行形态学运算的函数通常都会给出超出图像、未指定数值的像素指定一个数值,这样就类似于函数给图像填充了额外的行和列。对于膨胀和腐蚀操作,它们对像素进行填充的值是不同的。
对于二进制图像和灰度图像,膨胀和腐蚀操作使用的填充方法如下表:
腐蚀和膨胀填充图像规则表
| 规 |
腐蚀 | |
膨胀 | |
通过对膨胀操作使用最小值填充和对腐蚀操作使用最大值填充,可以有效地消除边界效应(输出图像靠近边界处的区域与图像其它部分不连续)。否则,如果腐蚀操作使用最小值进行填充,则进行腐蚀操作后,输出图像会围绕着一个黑色边框。
结构元素:膨胀和腐蚀操作的最基本组成部分,用于测试输出图像,通常要比待处理的图像小的多。二维平面结构元素由一个数值为0或1的矩阵组成。结构元素的原点指定了图像中需要处理的像素范围,结构元素中数值为1的点决定结构元素的邻域像素在进行膨胀或腐蚀操作时是否需要参与计算。三维或非平面的结构元素使用0,1定义结构元素在x和y平面上的范围,第三维z定义高度。
(1)任意大小和维数的结构元素B原点坐标的获取:
>> origin = floor((size(nhood)+1)/2)
其中nhood 是指结构元素定义的邻域(STREL对象的属性nhood)
(2)创建结构元素:(strel函数来创建任意大小和形状的STREL 对象,支持如线形line、钻石形diamond、圆盘形disk、球形ball等许多种常用的形状)
>> se = strel ('diamond'
se =
Flat STREL object containing 25 neighbors.
Decomposition: 3 STREL objects containing a total of 13 neighbors
Neighborhood:
(3) 结构元素的分解
为了提高执行效率,stel函数可能会将结构元素拆为较小的块,这种技术称为结构元素的分解。例如要对一个11×11的正方形结构元素进行膨胀操作,可以首先对1×11的结构元素进行膨胀操作,然后再对11×1的结构元素进行膨胀,通过这样的分解,在理论上可以使执行速度提高6.5倍。
对圆盘形和球形结构元素进行分解,其结构是近似的,而对于其他形状的分解,得到的分解结果是精确的。可以调用getsequence函数来查看分解所得的结构元素序列。
>> seq=getsequence(sel)
seq =
4x1 array of STREL objects
>> seq(1)
ans =
Flat STREL object containing 5 neighbors.
Neighborhood:
>> seq(2)
ans =
Flat STREL object containing 4 neighbors.
Neighborhood:
>> seq(3)
ans =
Flat STREL object containing 4 neighbors.
Neighborhood:
>> seq(4)
ans =
Flat STREL object containing 4 neighbors.
Neighborhood:
1.图像膨胀的Matlab实现:
可以使用imdilate函数进行图像膨胀,imdilate函数需要两个基本输入参数,即待处理的输入图像和结构元素对象。结构元素对象可以是strel函数返回的对象,也可以是一个自己定义的表示结构元素邻域的二进制矩阵。此外,imdilate还可以接受两个可选参数:PADOPT(padopt) ——影响输出图片的大小、PACKOPT(packopt).——说明输入图像是否为打包的二值图像(二进制图像)。举个实例如下:
步骤1,首先创建一个包含矩形对象的二值图像矩阵。
>> BW=zeros(9,10);
>> BW(4:6,4:7) =1
BW =
步骤2,使用一个3×3的正方形结构元素对象对创建的图像进行膨胀。
>> SE=strel('square',3)
SE =
Flat STREL object containing 9 neighbors.
Neighborhood:
步骤3,将图像BW和结构元素SE传递给imdilate函数。
>> BW2=imdilate(BW,SE)
BW2 =
步骤4,显示结果。
>> imshow(BW,'notruesize')
>> imshow(BW2,'notruesize')
膨胀前后效果图:
2.图像腐蚀的Matlab实现:
可以使用imerode函数进行图像腐蚀。imerode函数需要两个基本输入参数:待处理的输入图像以及结构元素对象。此外,imerode函数还可以接受3个可选参数:PADOPT(padopt) ——影响输出图片的大小、PACKOPT(packopt).——说明输入图像是否为打包的二值图像(二进制图像)。M——指定原始图像的行数。
以下程序示例说明了如何对某一副具体图像进行腐蚀操作,腐蚀前后的效果对比如图末。
步骤1,读取图像cameraman.tif (该图像是Matlab当前目录下自带的图片)
>> BW1=imread('cameraman.tif');
步骤2,创建一个任意形状的结构元素对象
>> SE=strel('arbitrary',eye(5));
步骤3,以图像BW1和结构元素SE为参数调用imerode函数进行腐蚀操作。
>> BW2=imerode(BW1,SE);
步骤4,显示操作结果
>> imshow(BW1)
>> figure,imshow(BW2)
图像cameraman.tif 腐蚀前后的效果对比:
3.膨胀和腐蚀联合操作(图像开运算操作):
下面以图像开启为例,说明如何综合使用imdilate和imerode这两个函数,实现图像处理操作。
步骤1,创建结构元素:
>> clear;close all
>> SE = strel('rectangle',[40 30]);
步骤2,使用结构元素腐蚀图像:
>> BW1=imread('circbw.tif');
>> BW2=imerode(BW1,SE);
>> imshow(BW2)
>> figure,imshow(BW1)
步骤3,恢复矩形为原有大小,使用相同的结构元素对腐蚀过的图像进行膨胀.
>> BW3=imdilate(BW2,SE);
>> figure,imshow(BW3)
最终效果如下图:
a.原始图像->b.腐蚀后的图像->c.膨胀后的图像
4.基于膨胀与腐蚀的形态操作——骨架化和边缘检测
(1)骨架化:
某些应用中,针对一副图像,希望对图像中所有对象简化为线条,但不修改图像的基本结构,保留图像基本轮廓,这个过程就是所谓的骨架化。提供了专门的函数bwmorph,可以实现骨架化操作。
>> clear;close all
>> BW1=imread('circbw.tif');
>> BW2=bwmorph(BW1,'skel',Inf);
>> imshow(BW1)
>>
(2)边缘检测
对于一副灰度二进制图像,如果图像像素值为1,则该像素的状态为ON,如果其像素值为0,则该像素的状态为OFF。在一副图像中,如果图像某个像素满足以下两个条件:
1.该像素状态为ON,
2.该像素邻域中有一个或多个像素状态为OFF。
则认为该像素为边缘像素。
Matlab中提供了专门的函数bwperim,可以用于判断一副二进制图像中的哪些像素为边缘像素。
以下程序代码示例就是利用bwperim函数,对图像circbw.tif进行边缘检测,其边缘像素检测效果如尾图。
>> clear;close all
>> BW1=imread('circbw.tif');
>> BW2=bwperim(BW1);
>> imshow(BW1)
>> figure,imshow(BW2)
>>
基于腐蚀和膨胀的形态操作函数如下:
bwhitmiss 图像逻辑"与"操作,该函数使用一个结构元素对图像进行腐蚀操作后,再使用第二个结构元素对图像进行腐蚀操作
imbothat 从原始图像中减去经过形态关闭后的图像,该函数可用来寻找图像中的灰度槽
imclose 闭合操作.首先对图像进行膨胀,然后再对膨胀后的图像进行腐蚀,两个操作使用同样的结构元素
imopen 开启操作,首先对图像进行腐蚀,然后再对腐蚀后的图像进行膨胀,两个操作使用同样的结构元素
imtophat 从原始图像中减去形态开启后的图像,可以用来增强图像的对比度
边缘检测算子:Roberts(2×2)、Prewitt(3×3)、Sobel(3×3)、Isotropic Sobel (3×3)、log、candy、zerocross。Matlab提供了专门的边缘检测edge函数。(用法详请参见Matlab图像处理函数大全)
四、所需用到的图像基本操作
1.图像类型转换:mat2gray(数据矩阵转灰度图像)、rgb2gray(彩色转灰度)、im2bw(将真彩色、索引色、灰度图转换成二值图)
I=mat2gray(A,[amin amax])
I=mat2gray(A)
I=rgb2gray(RGB)
newmap=rgb2gray(map)
BW=im2bw(I,level)
BW=im2bw(X,map,level)
BW=im2bw(RGB,level)
2.图像的块边缘操作:(对图像某一部分块进行操作,一次只处理一个图像像素)
像素的边缘是个像素集,由该像素(中心像素)的相关位置确定。像素的边缘就是一个矩形块,当操作从图像矩阵的一个元素移动到另一个元素时,该矩形块也以相向的方向同时移动。
通常,对于m×n的边缘来说,中心像素的计算方法:floor([(m+1 n+1)]/2)
Matlab对图像边缘操作的过程如下:
步骤1,选择像素,即读入所需要的图片像素。
>>I=imread('tire.tif');
步骤2,确定该像素的边缘块。
>>f=inline('max(x(:))');
步骤3,调用适当的函数对边缘中的元素进行操作
>>I2=nlfilter(I,[3 3],f);
步骤4,查找对应于输出图像中心像素的每个像素点,重复(1)~(4)步骤。
>>imshow(I);
>>figure,imshow(I2);
注意:本文是我完成本学期毕业设计所摘录的一些资料包含一些原理加个人思考总结的一些设计思路,不代表毕业设计的全部,希望对那些正在研究Matlab图像处理的朋友能够带来一些启示。