Matlab基于腐蚀和膨胀的边缘检测

1、Matlab基于腐蚀和膨胀的边缘检测2010年07月22日 星期四 16:15腐蚀：删除对象边界某些像素。膨胀：给图像中的对象边界添加像素。在操作中，输出图像中所有给定像素的状态都是通过对输入图像的相应像素及邻域使用一定的规则进行确定。在膨胀操作时，输出像素值是输入图像相应像素邻域内所有像素的最大值。在二进制图像中，如果任何像素值为1，那么对应的输出像素值为1；而在腐蚀操作中，输出像素值是输入图像相应像素邻域内所有像素的最小值。在二进制图像中，如果任何一个像素值为0，那么对应的输出像素值为0。结构元素的原点定义在对输入图像感兴趣的位置。对于图像边缘的像素，由结构元素定义的邻域将会有一部分位于图

2、像边界之外。为了有效处理边界像素，进行形态学运算的函数通常都会给出超出图像、未指定数值的像素指定一个数值，这样就类似于函数给图像填充了额外的行和列。对于膨胀和腐蚀操作，它们对像素进行填充的值是不同的。对于二进制图像和灰度图像，膨胀和腐蚀操作使用的填充方法如下表：腐蚀和膨胀填充图像规则表规 则腐蚀超出图像边界的像素值定义为该数据类型允许的最大值，对于二进制图像，这些像素值设置为1；对于灰度图像，unit8类型的最小值也为255。膨胀超出图像边界的像素值定义为该数据类型允许的最小值，对于二进制图像，这些像素值设置为0；对于灰度图像，unit8类型的最小值也为0。通过对膨胀操作使用最小值填充和对腐蚀

3、操作使用最大值填充，可以有效地消除边界效应（输出图像靠近边界处的区域与图像其它部分不连续）。否则，如果腐蚀操作使用最小值进行填充，则进行腐蚀操作后，输出图像会围绕着一个黑色边框。结构元素：膨胀和腐蚀操作的最基本组成部分，用于测试输出图像，通常要比待处理的图像小的多。二维平面结构元素由一个数值为0或1的矩阵组成。结构元素的原点指定了图像中需要处理的像素范围，结构元素中数值为1的点决定结构元素的邻域像素在进行膨胀或腐蚀操作时是否需要参与计算。三维或非平面的结构元素使用0，1定义结构元素在x和y平面上的范围，第三维z定义高度。(1)任意大小和维数的结构元素B原点坐标的获取： origin = flo

4、or(size(nhood)+1)/2) 其中nhood 是指结构元素定义的邻域（STREL对象的属性nhood）(2)创建结构元素：（strel函数来创建任意大小和形状的STREL 对象，支持如线形line、钻石形diamond、圆盘形disk、球形ball等许多种常用的形状） se = strel (diamond ,3)se = Flat STREL object containing 25 neighbors.Decomposition: 3 STREL objects containing a total of 13 neighborsNeighborhood:0 0 0 1 0 0

5、 00 0 1 1 1 0 00 1 1 1 1 1 01 1 1 1 1 1 10 1 1 1 1 1 00 0 1 1 1 0 00 0 0 1 0 0 0% se返回了结构元素的有关信息。(3) 结构元素的分解为了提高执行效率，stel函数可能会将结构元素拆为较小的块，这种技术称为结构元素的分解。例如要对一个1111的正方形结构元素进行膨胀操作，可以首先对111的结构元素进行膨胀操作，然后再对111的结构元素进行膨胀，通过这样的分解，在理论上可以使执行速度提高6.5倍。对圆盘形和球形结构元素进行分解，其结构是近似的，而对于其他形状的分解，得到的分解结果是精确的。可以调用getsequen

6、ce函数来查看分解所得的结构元素序列。 seq=getsequence(sel)seq =4x1 array of STREL objects seq(1)ans =Flat STREL object containing 5 neighbors.Neighborhood:0 1 01 1 10 1 0 seq(2)ans =Flat STREL object containing 4 neighbors.Neighborhood:0 1 01 0 10 1 0 seq(3)ans =Flat STREL object containing 4 neighbors.Neighborhood:0

7、 0 1 0 00 0 0 0 01 0 0 0 10 0 0 0 00 0 1 0 0 seq(4)ans =Flat STREL object containing 4 neighbors.Neighborhood:0 1 01 0 10 1 01.图像膨胀的Matlab实现：可以使用imdilate函数进行图像膨胀,imdilate函数需要两个基本输入参数，即待处理的输入图像和结构元素对象。结构元素对象可以是strel函数返回的对象，也可以是一个自己定义的表示结构元素邻域的二进制矩阵。此外，imdilate还可以接受两个可选参数：PADOPT(padopt) 影响输出图片的大小、PACK

8、OPT(packopt).说明输入图像是否为打包的二值图像(二进制图像)。举个实例如下：步骤1，首先创建一个包含矩形对象的二值图像矩阵。 BW=zeros(9,10); BW(4:6,4:7) =1BW =0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 1 1 1 1 0 0 00 0 0 1 1 1 1 0 0 00 0 0 1 1 1 1 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0步骤2，使用一个33的正方形结构元素对象对创建的图像

9、进行膨胀。 SE=strel(square,3)SE =Flat STREL object containing 9 neighbors.Neighborhood:1 1 11 1 11 1 1步骤3，将图像BW和结构元素SE传递给imdilate函数。 BW2=imdilate(BW,SE)BW2 =0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 00 0 1 1 1 1 1 1 0 00 0 1 1 1 1 1 1 0 00 0 1 1 1 1 1 1 0 00 0 1 1 1 1 1 1 0 00 0 1 1 1 1 1 1 0 00 0 0 0 0 0 0 0

10、 0 00 0 0 0 0 0 0 0 0 0步骤4，显示结果。 imshow(BW,notruesize) imshow(BW2,notruesize)膨胀前后效果图：2.图像腐蚀的Matlab实现：可以使用imerode函数进行图像腐蚀。imerode函数需要两个基本输入参数：待处理的输入图像以及结构元素对象。此外，imerode函数还可以接受3个可选参数：PADOPT(padopt) 影响输出图片的大小、PACKOPT(packopt).说明输入图像是否为打包的二值图像(二进制图像)。M指定原始图像的行数。以下程序示例说明了如何对某一副具体图像进行腐蚀操作，腐蚀前后的效果对比如图末。步骤

11、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这两个函

12、数，实现图像处理操作。步骤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.原始图像

13、-b.腐蚀后的图像-c.膨胀后的图像4.基于膨胀与腐蚀的形态操作骨架化和边缘检测（1）骨架化：某些应用中，针对一副图像，希望对图像中所有对象简化为线条，但不修改图像的基本结构，保留图像基本轮廓，这个过程就是所谓的骨架化。提供了专门的函数bwmorph，可以实现骨架化操作。 clear;close all BW1=imread(circbw.tif); BW2=bwmorph(BW1,skel,Inf); imshow(BW1) figure,imshow(BW2)（2）边缘检测对于一副灰度二进制图像,如果图像像素值为1,则该像素的状态为ON,如果其像素值为0,则该像素的状态为OFF。在一副图像

14、中，如果图像某个像素满足以下两个条件：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 图像逻辑与操作,该函数使用一个结

15、构元素对图像进行腐蚀操作后,再使用第二个结构元素对图像进行腐蚀操作imbothat 从原始图像中减去经过形态关闭后的图像,该函数可用来寻找图像中的灰度槽imclose 闭合操作.首先对图像进行膨胀,然后再对膨胀后的图像进行腐蚀,两个操作使用同样的结构元素imopen 开启操作,首先对图像进行腐蚀,然后再对腐蚀后的图像进行膨胀,两个操作使用同样的结构元素imtophat 从原始图像中减去形态开启后的图像,可以用来增强图像的对比度边缘检测算子：Roberts（22）、Prewitt（33）、Sobel（33）、Isotropic Sobel (33)、log、candy、zerocross。Mat

16、lab提供了专门的边缘检测edge函数。（用法详请参见Matlab图像处理函数大全）四、所需用到的图像基本操作1.图像类型转换：mat2gray(数据矩阵转灰度图像)、rgb2gray（彩色转灰度）、im2bw(将真彩色、索引色、灰度图转换成二值图)I=mat2gray(A,amin amax) %其中amin amax 是0,1间的任意值I=mat2gray(A) %数据矩阵转换为可视化灰度图像I=rgb2gray(RGB) %输入是真彩色，则输出图像类型可以是unit8或doublenewmap=rgb2gray(map) %输入输出图像都是double型的。BW=im2bw(I,leve

17、l) %灰度图像二值化，level是归一化的阈值，取值范围在0,1之间BW=im2bw(X,map,level) %索引色图像二值化(map为调色板)BW=im2bw(RGB,level) %彩色图像二值化2.图像的块边缘操作：（对图像某一部分块进行操作，一次只处理一个图像像素）像素的边缘是个像素集，由该像素（中心像素）的相关位置确定。像素的边缘就是一个矩形块，当操作从图像矩阵的一个元素移动到另一个元素时，该矩形块也以相向的方向同时移动。通常，对于mn的边缘来说，中心像素的计算方法：floor(m+1 n+1)/2)Matlab对图像边缘操作的过程如下：步骤1，选择像素，即读入所需要的图片像素

18、。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画二次曲面一、螺旋线1.静态螺旋线a=0:0.1:20*pi;h=plot3(a.*cos(a),a.*sin(a),2.*a,b,linewidth,2);axis(-50,50,-50,50,0,150);grid onset(h,

19、erasemode,none,markersize,22);xlabel(x轴);ylabel(y轴);zlabel(z轴);title(静态螺旋线); 2.动态螺旋线t=0:0.1:10*pi;i=1;h=plot3(sin(t(i),cos(t(i),t(i),*,erasemode,none);grid onaxis(-2 2 -2 2 0 35)for i=2:length(t) set(h,xdata,sin(t(i),ydata,cos(t(i),zdata,t(i); drawnow pause(0.01)endtitle(动态螺旋线);(图略) 3.圆柱螺旋线t=0:0.1:1

20、0*pi;x=r.*cos(t);y=r.*sin(t);z=t;plot3(x,y,z,h,linewidth,2);grid onaxis(square)xlabel(x轴);ylabel(y轴);zlabel(z轴);title(圆柱螺旋线) 二、旋转抛物面b=0:0.2:2*pi;X,Y=meshgrid(-6:0.1:6);Z=(X.2+Y.2)./4;meshc(X,Y,Z);axis(square)xlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(旋转抛物面)或直接用：ezsurfc(X.2+Y.2)./4) 三、椭圆柱面lo

21、ad clownezsurf(2*cos(u),4*sin(u),v,0,2*pi,0,2*pi)view(-105,40) %视角处理shading interp %灯光处理colormap(map) %颜色处理grid on %添加网格线axis equal %使x,y轴比例一致xlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(椭圆柱面) %添加标题四、椭圆抛物面b=0:0.2:2*pi;X,Y=meshgrid(-6:0.1:6);Z=X.2./9+Y.2./4;meshc(X,Y,Z);axis(square)xlabel(x轴);

22、ylabel(y轴);zlabel(z轴);shading flat;title(椭圆抛物面)或直接用：ezsurfc(X.2./9+Y.2./4)五、双叶双曲面ezsurf(8*tan(u)*cos(v),8.*tan(u)*sin(v),2.*sec(u),-pi./2,3*pi./2,0,2*pi)axis equalgrid onaxis squarexlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(双叶双曲面)六、双曲柱面load clownezsurf(2*sec(u),2*tan(u),v,-pi/2,pi/2,-3*pi,3

23、*pi)hold on %在原来的图上继续作图ezsurf(2*sec(u),2*tan(u),v,pi/2,3*pi/2,-3*pi,3*pi)colormap(map)shading interpview(-15,30)axis equalgrid onaxis equalxlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(双曲柱面)七、双曲抛物面（马鞍面）X,Y=meshgrid(-7:0.1:7);Z=X.2./8-Y.2./6;meshc(X,Y,Z);view(85,20)axis(square)xlabel(x轴);ylabel

24、(y轴);zlabel(z轴);shading flat;title(双曲抛物面)或直接用：ezsurfc(X.2./8-Y.2./6) 八、抛物柱面X,Y=meshgrid(-7:0.1:7);Z=Y.2./8;h=mesh(Z);rotate(h,1 0 1,180) %旋转处理%axis(-8,8,-8,8,-2,6);axis(square)xlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(抛物柱面)或直接用：ezsurfc(Y.2./8) 九、环面ezmesh(5+2*cos(u)*cos(v),(5+2*cos(u)*sin(v

25、),2*sin(u),0,2*pi,0,2*pi)axis equalgrid onxlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(环面)十、椭球ezsurfc(5*cos(u)*sin(v),(3*sin(u)*sin(v),4*cos(v),0,2*pi,0,2*pi)axis equalgrid onxlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(椭球)十一、单叶双曲面ezsurf(4*sec(u)*cos(v),2.*sec(u)*sin(v),3.*tan(u),-pi./2,pi./2,0,2*pi)axis equalgrid onxlabel(x轴);ylabel(y轴);zlabel(z轴);shadi