基于分水岭和形态学的图像分割算法研究毕业设计.doc

《基于分水岭和形态学的图像分割算法研究毕业设计.doc》由会员分享，可在线阅读，更多相关《基于分水岭和形态学的图像分割算法研究毕业设计.doc（20页珍藏版）》请在三一办公上搜索。

1、本 科 生 毕 业 设 计（申请学士学位）论文题目基于分水岭和形态学的图像分割算法研究 作者姓名 专业名称 指导教师 2013年6月学 生： （签字）学 号： 答 辩 日 期：2013年6月15日指 导 教 师： （签字）目录摘要1Abstract.11 绪论21.1 研究目的和意义21.2 图像分割的研究进展21.3 论文主要内容和组织结构32 数学形态学32.1 膨胀与腐蚀32.1.1 灰度膨胀42.1.2 灰度腐蚀42.2 形态学的开运算和闭运算42.3 形态学重建53 基于分水岭和形态学的图像分割算法53.1 分水岭算法原理63.2 形态学算子的改进73.3 改进的图像分割算法描述83

2、.4 实验结果与数据分析84 结束语10参考文献11附录11致谢14基于分水岭和形态学的图像分割算法研究 摘要：图像分割是图像分析和处理中一个重要的研究方向，也是目标的检测和识别的重要步骤。而且由于图像的多样性和复杂性，很难用统一的方法来描述感兴趣的对象，因而在实际应用中只能根据各种领域的需求来选择合适的分割方法，导致各种图像分割方法具有特定的局限性和针对性。目前还没有一种通用的方法，能使各种类型的图象达到最优分割质量。传统的图像分割算法中对图像噪声敏感，噪声会直接恶化图像的梯度图像，计算量大，分割过程耗时长，分割效率低，从而产生过分割问题。为了降低过分割现象对图像分割的影响，提出了一种改进的

3、分水岭算法的图像分割方法。本文先进行分水岭变换，再利用数学形态学的方法，使用多尺度形态梯度算子，利用结构元素度优点以达到改善图像分割中的过分割现象。实验结果表明，改进的算法有效地改善了过分割现象，具有较好的图像分割效果。关键字：图像分割；分水岭算法；形态学算子Image segmentation algorithm based on watershed and morphological studyAbstract:Image analysis and image segmentation is an important research direction, also is the impo

4、rtant process of target detection and recognition. And because of the diversity and complexity of the image, it is difficult to use uniform method to describe the object of interest, and therefore can only according to the various fields in the practical application needs to choose the appropriate m

5、ethod, lead to all kinds of image segmentation method has certain limitations and pertinence. There is no a common method, can make various types of image to achieve the optimal segmentation quality.The traditional image segmentation algorithm is sensitive to image noise, the noise will deteriorate

6、image gradient image directly, large amount of calculation, the segmentation process takes long, segmentation efficiency is low, resulting in a over-segmentation problem. In order to reduce the over-segmentation phenomena influence on image segmentation, an improved watershed algorithm is proposed f

7、or image segmentation method. Watershed transform first in this paper, and then the mathematical morphology method, using multi-scale morphological gradient operator, using the structure elements of advantage to improve the image segmentation of over-segmentation phenomenon. The experimental results

8、 show that the improved algorithm improved the over-segmentation phenomena effectively, has the good image segmentation effect.Key words: Image segmentation; Watershed algorithm; Morphological operator1 绪论1.1 研究目的和意义图像分割是一种重要的图像技术，在理论研究和实际应用中都得到了人们的广泛重视。图像分割的方法和种类有很多，有些分割运算可直接应用于任何图像，而另一些只能适用于特殊类别的图

9、像。许多不同种类的图像或景物都可作为待分割的图像数据，不同类型的图像，已经有相对应的分割方法对其分割;但某些分割方法只是适合于某些特殊类型的图像分割，所以分割结果的好坏需要根据具体的场合及要求衡量。图像分割是从图像处理到图像分析的关键步骤，可以说，图像分割结果的好坏直接影响对图像的理解。图像分割是由图像处理到图像分析的关键步骤，在图像工程中占有重要位置。一方面，它是目标表达的基础，对特征测量有重要的影响。另一方面，因为图像分割及其基于分割的目标表达、特征提取和参数测量等将原始图像转化为更抽象、更紧凑的表达形式，使得更高层的图像分析和理解成为可能。因此在实际应用中，图像分割不仅仅要把一幅图像分成

10、满足上面五个条件的各具特性的区域，而且要把其中感兴趣的目标区域提取出来。只有这样才算真正完成了图像分割的任务，为下一步的图像分析做好准备，使更高层的图像分析和理解成为可能。图像分割在很多方面，如医学图像分析，交通监控等，都有着非常广泛的应用，具有重要的意义。(1)分割的结果常用于图像分析，如不同形式图像的配准与融合，结构的测量，图像重建以及运动跟踪等。(2)在系统仿真，效果评估，图像的3D重建以及三维定位等可视化系统中，图像分割都是预处理的重要步骤。(3)图像分割可在不丢失有用信息的前提下进行数据压缩，这就降低了传输的带宽，对提高图像在因特网上的传输速度至关重要。(4)分割后的图像与噪声的关系

11、减弱，具有降噪功能，便于图像的理解。1.2 图像分割的研究进展图像分割是图像处理中的一项关键技术，至今已提出上千种分割算法。但因尚无通用的分割理论，现提出的分割算法大都是针对具体问题的，并没有一种适合所有图像的通用分割算法。前人的方法主要有三大类：阈值分割方法、边缘检测方法和区域提取方法。（1）阈值分割方法阈值分割法分为全局阈值法和局部阈值分割法。所谓局部阈值分割法是将原始图像划分成较小的图像，并对每个子图像选取相应的阈值。在阈值分割后，相邻子图像之间的边界处可能产生灰度级的不连续性，因此需用平滑技术进行排除。局部阈值法常用的方法有灰度差直方图法、微分直方图法。局部阈值分割法虽然能改善分割效果

12、，但存在以下几个缺点：每幅子图像的尺寸不能太小，否则统计出的结果无意义。每幅图像的分割是任意的，如果有一幅子图像正好落在目标区域或背景区域，而根据统计结果对其进行分割，也许会产生更差的结果。局部阈值法对每一幅子图像都要进行统计，速度慢，难以适应实时性的要求。全局阈值分割方法在图像处理中应用比较多，它在整幅图像内采用固定的阈值分割图像。经典的阈值选取以灰度直方图为处理对象。根据阈值选择方法的不同，可以分为模态方法、迭代式阈值选择等方法。这些方法都是以图像的直方图为研究对象来确定分割的阈值的。另外还有类间方差阈值分割法、二维最大熵分割法、模糊阈值分割法以及共生矩阵分割法等等。 （2）基于边缘检测法

13、边缘检测法是一种处理不连续性图像的分割技术。图像的大部分信息不会只存于某个特定的区域，而是存于不同区域的边缘上，而且人的视觉系统在很大程度上都是根据图像边缘差异对图像进行识别分析的。所以通过对图像的边缘信息检测，可以实现对图像的分割。按照处理技术可以分为并行边缘检测技术和串行边缘检测技术。检测过程中可以通过空域微分算子来完成卷积。这些微分算子包括、Sobel梯度算子、Prewitt梯度算子、综合正交算子等。这些方法主要是对检测图像中灰度的变化，图像边缘是灰度突变的地方。在有噪声时，得到的边缘常是孤立不连续的，为了得到完整的边缘信息，还需进行边界闭合处理。边界闭合是根据像素梯度的幅度及梯度方向满

14、足规定的条件将边缘素连接起来，就有可能得到闭合的边界。（3）基于区域分割法区域分割法主要包括：区域生长和分裂合并法，其分割过程后续步骤的处理要根据前面步骤的结果进行判断而确定。区域生长法是根据预先规定好的指标，提取图像中相互连接区域的方法，它是利用区域一致性准则对目标进行分割。规定的指标包括图像的灰度信息，边缘，某种特性。区域生长法一般都会放在一系列过程中使用，不会单独使用。它主要的缺陷是，每一个需要提取的区域，都必须先给出种子点，然后提取出和种子一样，符合规定的指标的区域，这样有多少区域就必须给出多少个种子数。这种法对噪声也很敏感，会造成分割区域不连续。相反的，局部且大量的噪声会使影响会使原

15、本来分开的区域连接起来。分裂合并法是从整个图像出发，将图像分割成各个子区域，再把前景的区域合并起来，这样就实现了目标的提取。分裂合并法的目标区域由一些相互连通的像素组成的，如果把图像分割到像素级的话，就可以判断该像素是否为目标像素。当所有的分割的子区域都判断完，把目标区域就可得到前景目标。这种方法处理复杂图像时效果较好，但算法比较复杂，计算量也比较大，在分裂过程中可能会破坏目标区域的边界。1.3 论文主要内容和组织结构本文详细介绍了数学形态学和分水岭算法原理，对基于改进分水岭的图像分割算法进行了较为详细的探讨和研究，突出了分水岭图像分割算法在图像处理中的优点；重点研究了形态学算子的改进应用于分

16、水岭算法的图像分割，并通过不同分割方法得到的分割结果，对分割结果进行了分析。本文共分为3章，各章内容安排如下：第1章绪论部分，简要介绍了本课题研究的目的和意义，重点分析了图像分割的研究进展。第2章详细介绍了数学形态学，了解了数学形态学中的膨胀与腐蚀以及形态学的开闭运算。第3章详细分析了分水岭分割原理，针对基于改进分水岭的图像分割算法进行了探讨，利用改进的形态学算子应用分水岭算法的图像分割，并进行了实验验证：针对研究图像，利用不同的算法做比较实验，并对实验结果及数据进行了综合分析。第4章对本文的主要工作进行了小结，并对图像分割的前景进行了展望。2 数学形态学数学形态学是数学基础和所用语言是集合论

17、。数学形态学的算法具有天然的并行实现的结构，实现了形态学分析和处理算法的并行，大大提高了图像分析和处理的速度。在有关数学形态学的许多研究成果中，其已经影响到包括图像增强、分割、恢复、边缘检测、纹理分析、颗粒分析、特征生成、骨架化、形状分析、压缩、成分分析及细化等诸多领域。2.1 膨胀与腐蚀将二值数学形态学推广至灰度形态学,则在以下讨论中将处理的是数学图像函数而不是集合。设是输入图像,是结构元素,它可被看作是一个子图像函数。如果z表示的是实整数集合,同样假设是来自的整数,I和b是对坐标为象素灰度值的函数(来自实数集R的实数)。如果灰度也是整数,则Z可由整数R所代替。最基本的数学形态学算子有膨胀和

18、腐蚀或者由这两种基本算子复合而成的其它的数学形态学算子,如形态开运算算子、形态闭运算算子等。下面列举一些基本的用于灰度图像处理的数学形态学运算。设I表示灰度图像,则它可定义为根据组成的的子集,其中N是之间的正整数。2.1.1 灰度膨胀灰度膨胀：用b对灰度图像I进行的灰度膨胀表示为， （3-1） 其中，和分别是b和I的定义域。 根据在每个结构元素的位置,这一点的膨胀值是在跨度为b的区间内I和b之和的最大值,是局部最大值滤波。2.1.2 灰度腐蚀 灰度腐蚀：用b对灰度图像I进行的灰度腐蚀表示为， （3-2） 其中, 是分别是b和I的定义域。 腐蚀操作是以在结构元素形状定义的区间内最小值为基础的,是

19、局部最小值滤波。 源图像a a膨胀后图像 a腐蚀后图像图2-1细胞图像a经过膨胀和腐蚀后的结果2.2 形态学的开运算和闭运算 开运算一般能平滑图像的轮廓,削弱狭窄的部分,去掉细的突出。闭运算也是平滑图像的轮廓，与开运算相反,它常常能融合窄的缺口和细长的弯口,去掉小洞,填补轮廓上的缝隙。设I是原始图像,b是结构元素图像,则集合I被结构元素b作开运算,记为,其定义为: （3-3）就是I被b腐蚀后的结果再被b膨胀。设I是原始图像,b是结构元素图像,则集合I被结构元素b作闭运算,记为,其定义为: （3-4）就是I被b膨胀后的结果再被b腐蚀。 源图像b b开运算后图像 b闭运算后图像图2-2细胞图像b经

20、过开闭运算后的结果2.3 形态学重建形态学重建是形态学图像处理的另一个重要的内容。形态学重建具有下面特性：（1）处理基于两幅图像，而不是一幅图像和一个结构元素；（2）进行重复处理直到稳定(即图像不再改变)；（3）处理基于连通性，而不是基于结构元素；所谓形态学重建，就是根据一幅图像(称之为掩模图像)的特征对另一幅图像(称之为标记图像)进行重复的膨胀，直到该图像的像素值不再变化，图像稳定为止。本文引入测地膨胀的概念来理解灰度图像形态学重建。根据阈值分解原理，灰度图像的基本测地膨胀定义为： （3-5）其中代表逐点比较取最小值，B为尺寸为1的结构元素，是表示J被B进行膨胀运算。因此灰度图像形态学重建可

21、定义为： （3-6） （3-7）其中代表逐点比较取最大值。用J对I的灰度图像形态学重建可通过在不超过I的前提下反复迭代I的基本灰度测地膨胀直至稳定。形态学重建是图像形态处理的重要操作之一。通常用来强调图像中与掩模图像指定对象相一致的部分，同时忽略图像中的其他对象。通过形态学重建可以有效地去除由噪声和量化误差所引起的由少量象素或者邻近象素值差别很小的局部“谷底”，进一步抑制噪声点。3 基于分水岭和形态学的图像分割算法分水岭分割算法是一种常见的算法，其优点在于它可以得到封闭、连续的物体轮廓线，而且当图像的边缘信息不清或微弱时也依然有很好的分割效果。算法中对图像噪声敏感，噪声会直接恶化图像的梯度图像

22、，计算量大，分割过程耗时长，分割效率低，从而产生过分割问题。为此，需进一步研究图像分割算法，以期望获得更好的分割结果。3.1 分水岭算法原理分水岭算法最早是由Beucher和Lantuejoul 在20世纪70年代末提出的“浸没”模型并引入到图像处理领域。算法在已知图像每个区域极小值影响的范围内，通过形态学厚化运算，逐渐扩展所影响的区域范围得到分水线。分水岭算法的直观概念来源于地理学，其基本思想是把图像看作地貌。对分水岭算法的一种理解是模拟降水的过程：不断下落的雨水降沿着山地表面的地势向下流动，如果雨水降落在2个不同的点上，最终都会流到相同的局部最低点，则这2个点属于同一区域。只有山脊上的雨水

23、流入周围的区域的概率是相同的。因此，水降落在山地表面能流入到相同局部最低点的表面形成一个集水盆地，即分割区域，而分割集水盆地的山脊就是区域边界，即分水岭。其工作原理如图3-1所示。图3-1 模拟降水的过程 另一种理解方法是想象泛洪的过程：该地貌正被浸入一个大水池中，它局部有些极小值处有供水流入的小孔，水流从局部极小值处逐渐填满集水盆地，当不同集水盆地的集水汇合在一起的时候，就会在这些相遇点建立起水坝。而当水平线达到该地貌最高点时，整个地貌的浸入过程结束。最后，该地貌就被整个水坝分成了不同的集水区域或集水盆地，这些水坝就被称之为分水岭。其工作原理如图3-2所示。图3-2模拟泛洪的过程可以用两种最

24、基本的方法来模拟该过程以进行图像分割：一种是通过取积水盆地的补集得到分水岭；另一种则是把一个完整的将图像分成集水盆地的分割计算出来，再利用边缘检测得到分水岭。在实际操作中，通常不会将分水岭变换应用于原始图像，而是把它应用于原始图像的梯度或是数学形态学梯度图像。这样做会在灰度值不连续的像素点处得到分水岭，而这正是图像分割想要得到的。3.2 形态学算子的改进经典的形态学算子在图像处理中有很多的不足之处，本节讨论其中的改进办法，以适应分水岭算法。形态学算子定义为： （3-8） 其中 B 是以原点为中心所命名的扁平结构元素。因为利用扁平结构元素作腐蚀和膨胀可以得到极大极小滤波器的效果，故在每一个点，形

25、态学算子都可以得到由扁平结构元素所确定的领域上的极大值和极小值的差值。其中，“模糊”边缘可以用斜坡边缘来建模，并将边缘两边之间的像素灰度差可称之为“边缘高度”。对于斜坡边缘，由于传统梯度算子的输出是边缘的斜率，因此，如果边缘的斜率不大，则不能用取阈值的方法来区分斜坡边缘与噪声或量化误差。简而言之，其性能取决于结构元素 B 的大小。即如果 B 足够大，则对斜坡边缘来说，这个梯度算子的输出即等于边缘高度，可遗憾的是，大的结构元素会造成边缘间严重的互相影响，这将导致梯度极大值与边缘的不一致；然而，若结构元素过小，则梯度算子虽有高的空间分辨率，但对斜坡边缘会产生一个很小的输出结果。为了达到利用大结构元

26、素和小结构元素的各自优点的目的，提出以下一种改进的算子即多尺度的形态梯度算子。设为一组正方形的结构元素，的大小为，则多尺度梯度定义为 （3-9）对于阶跃边缘，运算即产生一条与边缘重合的2个像素宽的直线。由于该直线的高度等于边缘高度，因此，在这种情形下，多尺度的形态梯度算子即相当于单尺度的形态梯度算子。在实践中，由于多尺度的形态学梯度算子适用了取平均运算，所以抗噪声的能力更强。对于斜坡边缘，用h和w分别表示边缘的高度和宽度，而运算则产生一条对应于边缘的线。当时，则该线为梯形；当时，则该线为三角形，梯形或三角形的底边宽度为个像素，梯形的高度为，而三角形的高度为，且它们都大于边缘斜率。如果n的取值足

27、够大，则的值近似等于。此外，多尺度梯度算法有很强的抗边缘间相互影响的能力，如对应于某一边缘梯度极大值的位置不会因其他边缘的存在而受到影响。两个相邻边缘间的距离用d来表示，当时，则运算则产生与边缘一致的梯度极大值；当时，则运算的结果是在两个相邻梯度极大值之间的间隙用较小的极大值填充，然后通过对所有i的值取平均值，则在恰当的位置保持梯度极大值不变。对于斜坡边缘间的相互影响，情况与上述类似。3.3 改进的图像分割算法描述分水岭分割算法的优点在于它可以得到封闭、连续的物体轮廓线，而且当图像的边缘信息不清或微弱时也依然有很好的分割效果。然而算法中存在着过分割现象，难以根据分割出来的轮廓识别实际的目标。使

28、用改进的形态学算子形态学在其中可以改善这一现象，即在利用分水岭算法分割图像前,对图像进行有限次的膨胀、腐蚀、开闭运算和重建等形态学操作,在一定程度上降低了噪声和灰度不连续的影响。结合了形态学和分水岭算法各自的优势，本文提出了基于形态学和分水岭的图像分割改进算法，算法步骤如下描述。（1）首先计算图像的梯度图，对目标图像进行灰度化处理。（2）使用多尺度算子计算形态学梯度,利用迭代法得到合适的阈值。（3）求滤波后图像的局部极大值二值图像；对局部极大值二值图像进行开闭运算平滑边缘，并去掉图像中像素个数小于 20 的局部极小值区域，得到前景标记图；对第三步的滤波图像用迭代法求阈值，然后进行二值化，得到二

29、值图像。对该二值图像进行距离变换，分水岭分割，得到的分水线图像作为背景标记；（4）利用形态学极小值标定技术修改第二步得到的梯度图像，在根据前景标记和背景标记修改梯度图像，进行分水岭分割；（5）利用评价指标分割结果进行质量评价。3.4 实验结果与数据分析 为了客观地验证本文算法的有效性，使用迭代法、分水岭分割方法与本文方法进行对比实验。即A为迭代法，B为分水岭分割方法，以及本文方法。实验软件采用Matlab平台进行验证，所有实验都在主频为2.0G的CPU（双核)和内存1GB DDRRAM的计算机上进行，操作系统使用Windows XP（SP3）。本文实验图像选取三组图像，尺寸均为300300。如

30、图3-3所示，其中图1是梨图像，图中梨表面有虫眼，需要将其分割出来，从而有利于梨表面缺陷检测；图2是血细胞图像，其细胞具有可变形性，形状并不固定，需要对血细胞的细胞以及细胞内的物质进行分割，对医学临床检测具有着重要意义；图3是白细胞图像，细胞存在着细胞浆与细胞核，需要将其细胞核分割出来，对医学的后续研究有着重要意义。 （d）图像1 （e）图像2 (f)图像3图3-3 实验选取的三组图像将实验图像进行对比实验，得到三组结果。如图3-4所示，即A为迭代法，B为分水岭分割方法，以及本文方法。 图像1 A 方法分割结果 B方法分割结果 本文方法分割结果 图像2 A 方法分割结果 B 方法分割结果 本文

31、方法分割结果 图像3 A 方法分割结果 B 方法分割结果 本文方法分割结果图3-4 不同算法分割的结果为了评价融合图像的质量，本文从主观视觉并结合了算法时间（Time,T）、区域内部均匀性(Uniformity Measure，UM)和区域对比度(Gray-Level Contrast，GC)等客观因素的不同方法对图像分割进行评价。从主观地评价分割图像的质量：主要是从图像分割后的分割图像边缘的完整程度来评价。T是指图像分割前后的时间差，T越小表明算法分割效率越快。从客观地评价融合图像的质量：其中UM是指把原始图像分成具有相似特性的区域的内部均匀的程度，用表示第个区域，表示其面积，则分割图像的均

32、匀性UM可表示为： (3-10) 其中C为归一化系数，在二值图像中即为2，UM的值越大。分割效果就越好，所用的算法分割效果也就越好。GC是指图像分割的目标和背景两个区域的灰度分布特性的对比度。对图像分割中的目标和背景这两个区域来说，如果用和分别表示它们各自的平均灰度，则它们之间的灰度对比度GC可表示为： (3-11)由上述分析可知区域对比度GC越大，则分割效果就越好表3-1给出了上述三组实验图像分割结果的性能评价数据。表3-1实验数据统计结果分割方法算法时间T区域对比度GC区域内部均匀性UM边缘保持图像1A方法2.3415.81667.1578欠佳B方法2.1416.55777.5634一般本

33、文方法1.7619.23467.9878良好图像2A方法1.9116.22347.3456欠佳B方法1.4817.66017.5676一般本文方法1.1219.65817.8997良好图像3A方法2.0317.81007.1231欠佳B方法1.8318.11337.4632一般本文方法1.6719.34347.8567良好从上述图3-1实验结果和数据可以看出，本文所采用的改进分割方法不论是在主观视觉还是客观质量评价指标上，都明显优于迭代法和分水岭分割方法。从三组数据上可以很直观的看到本文的分割速度较其他算法有很大的提高，能够较快的将目标对象分割出来，分割效率比其他方法提高了很多，而且分割图像内

34、图像边缘较完整，对于图像1实验结果，迭代法分割方法并没有很好将表面有虫分割出来，不利于梨表面的缺陷检测，相比于本文图像分割方法还是有较大的差距，本文的分割方法结果能够准确的分割出表面有虫的部位，有利于有利于梨表面缺陷检测。其中本文算法时间比迭代法在算法时间上降低了24.8%，区域对比度上提高了21.6%，本文方法比分水岭分割方法在算法时间上降低了17.8%，而在区域对比度上提高了16.2%，在区域内部均匀性上本文方法明显比迭代法提高了11.6%，较分水岭分割方法法提高了5.6%，这说明分割效率有了显著提高；对于图像2，为研究血细胞的细胞以及细胞内的物质进行分割，将细胞以及细胞内的物质进行分割出

35、来，迭代法和分水岭分割方法分割效果差，区域对比度不高，缺少了细胞内的一些物质，导致血细胞产生去分割现象，本文分割方法很好的保留了血细胞图像分割的完整性，从数据上可以看出本文图像分割方法的区域对比度和区域内部均匀性较迭代法分别提高了21.2%和7.5%，比分水岭分割方法提高了11.3%和4.4%，本文图像分割方法的算法时间比迭代法降低了41.4%，较分水岭图像分割方法降低了24.3%；对于图像3，本文方法区域对比度较迭代法提高了8.6%，较分水岭图像分割方法提高了6.8%，区域内部均匀性较迭代法提高了10.3%，较分水岭图像分割方法提高了5.3%，本文方法的算法时间比迭代法和分水岭图像分割方法分

36、别降低了17.7%和8.7%，从图像上看迭代法和分水岭图像分割方法有明显的过分割现象，而本文提出的改进的图像分割方法可以得到完整并且连续的图像边缘且图像边缘完整、区域对比度高，其分割结果更具有实用性。4 结束语本文首先对图像分割的定义和研究发展进行了概述，然后叙述了数学形态学膨胀和腐蚀以及开闭运算，在以数学形态学的基础上引入分水岭变换。分水岭变换具有计算数度较快，边界定位精确，在图像分割中得到物体轮廓线的优点是具有封闭性，但分水岭变换对灰度变化有很大的敏感性，常常使分割结果出现过分割现象，难以得到很满意的分割效果。如何克服分水岭变换的过分割现象的是图像分割领域的一个研究热点，至今未有方法很好的

37、解决这个问题。本文先进性分水岭变换，再利用数学形态学的方法，得到最终分割结果。用灰度图像进行测试，测试结果表明改进的算法能够有效改善过分割现象，分割结果更有意义。基于分水岭和形态学的图像分割方法具有广泛的应用前景，虽然到目前为止，这种技术还不够成熟，相信通过人们坚持不懈的努力探索，分水岭变换在图像分割领域必定获得更大的发展。参考文献1Liang K H, Tjahajadi T, Yang Y H. Roof Edge Detection Using Regularised CubicB-spelins Filting. Pattern Recognition, 1997, 30(5): 71

38、9-7282Michael Kass, Andrew Witkin, Demetri Terzopoulos.Snake:Active Contour Models.Internatinal Journal of Computer Vision,321-331(1988).3Gashtashy A.Design and Recovery of 2-D and 3-D Shapes Using Rational GanssianCurves and Surfaces.International Journal of Computer Vision,1993,10(3):233-256. 4崔屹.

39、图像处理与分析数学形态学方法及应用M.北京:科学出版社,2000.5宋利伟,宋朝昀,庄天戈.基于分水岭算法的磁共振脑图像自动分割J.上海交通大学学报,2003,37(11):1754-1756.6王蜀,李永宁,陈楷民,黄戈.基于数学形态学的医学图像分割J.计算机应用，2005,25(10):2381-2386.7王鹏伟，吴秀清，张名成.基于多尺度形态学融合的分水岭图像分割方法J.数据采集与处理，2006,21(4):398-402.8冯林,管慧娟,孙焘,滕弘飞.基于分水岭变换和核聚类算法的图像分割J.大连理工大学学报,2006,46(6):851-856.9白相志,周付根.基于改进形态学算子的

40、多尺度边缘检测J.中国图像图形学报，2007,12(9):1610-1613.10梁萍,程伟.基于分水岭的粘连雾滴图像分割方法J.重庆工学院学报,2008,22(10):171-178.11杨洪薇,肖志涛,翁秀梅.基于分水岭和模糊C均值聚类的图像分割方法J.天津工业大学学报,2008,27(1):53-55.12冯林,孙焘,吴振宇,管慧娟.基于分水岭变换和图论的图像分割方法J.仪器仪表学报,2008,29(3):649-653.13孙敬飞,杨红卫.形态学分水岭算法在粘连大米图像分割中的应用J.河南工业大学学报,2009,30(6):84-87.14陈洁,胡永,刘泽国.基于标记的分水岭图像分割算

41、法研究J.软件杂志,2012,33(9):115-117.15纪凤仙.分水岭算法在医学图像分割中的改进方法研究J.技术应用,2012,1(4):89-91.附录1.程序源代码（1）经典迭代法源程序clear;close all; I=imread(d1.jpg); if length(size(I)=3 I=rgb2gray(I);endX,Y=size(I);T=0.5*(double(min(I(:)+double(max(I(:);%初始分割阈值done=false;whiledone%迭代分割 g=I=T; Tnext=0.5*(mean(I(g)+mean(I(g); done=ab

42、s(T-Tnext)0.5; T=Tnext;endfor i=1:X for j=1:Y if(I(i,j)=T) I(i,j)=0; else I(i,j)=1; end endendI=double(I);imshow(I);（2）分水岭图像分割算法源程序clear;close all; I=imread( pears.png ); if length(size(I)=3 I=rgb2gray(I);endX,Y=size(I);V_max=max(max(I);V_min=min(min(I);T=graythresh(I)%运用graythresh函数计算阈值TmaxTmax=T*2

43、55for i=1:X for j=1:Y if(I(i,j)=Tmax) I(i,j)=0; else I(i,j)=1; end endendI=double(I);imshow(I);（3）本文算法源程序clear alla=imread(10,jpg);figure,imshow(a)count=imhist(a);m,n=size(a);N=m*n;L=256;count=count/N; countfor i=1:L if count(i)=0 st=i-1; break; endendstfor i=L:-1:1 if count(i)=0 nd=i-1; break; endendndf=count(st+1:nd+1); %f是每个灰度出现的概率size(f)E=;for Th=st:nd-1 %设定初始分割阈值为Thav1=0;av2=0;Pth=sum(count(1:Th+1);%第一类的平均