医学图像真皮区域分割算法

《医学图像真皮区域分割算法》由会员分享，可在线阅读，更多相关《医学图像真皮区域分割算法（27页珍藏版）》请在装配图网上搜索。

1、 课 程 现代机器学习 题 目 医学图像真皮区域分割算法 学 号 1401120522 专 业 电子与通信工程 姓 名 董文欣 医学图像真皮区域分割算法研究摘要图像分割是一种重要的和关键的图像分析技术，也是当今医学领域备受研究人员关注的热点问题。其目的是把图像分成各具特性的区域并提取出感兴趣的部分，图像分割的结果是图像特征提取和识别等图像理解的基础。本文在前人研究成果的基础上，简述医学病理切片图像中真皮区域分割的实现流程，研究并实现基于K-means算法和live wire算法的医学病理切片图像中真皮区域分割算法，引用了重合比率参数比较分析两种分类算法的优缺点。关键词：医学图像区域分割，K-均

2、值聚类，重合率 目录1. 绪论11.1 引言11.2 课题研究的背景及意义11.3 医学图像分割的研究现状21.4 论文的主要工作及结构安排32. 颜色模型基本概念42.1 引言42.2 颜色模型定义42.3 颜色模型分类42.4 颜色模型的选择和转换62.5本章小结73. 医学图像分割概述83.1 引言83.1.1 医学图像分割特点83.1.2 医学图像分割定义93.2 医学图像分割典型算法103.3 本章小结134. 医学病理切片中真皮组织区域分割算法144.1引言144.2 基于K-means医学图像分割算法144.3 算法评估与分析184.4 本章小结215. 总结与展望225.1 总

3、结225.2 展望22参考文献24I 1. 绪论1.1 引言图像处理的目的应是符合对图像的正确理解，即对图像中特定物体的正确认识，以指导下一步的研究工作。如图1-1所示，依据研究方法的不同，图像技术具体可以分为三个层次1：处理、分析和理解。在整个过程中，图像分割（Image Segmentation）是最关键的一步。图1-1 图像技术的层次结构图像分割是一种关键的图像技术，它不仅得到人们的广泛重视，也在实际中得到大量的应用。图像分割能够自动或半自动描绘出医学图像中的解剖结构和其它感兴趣的区域, 从而有助于医学诊断。图像分割在不同领域中有时也用其它名称，如阈值化技术、目标轮廓技术、图像区分或求差

4、技术、目标跟踪技术、目标识别技术、目标检测技术等，这些技术本身或核心实际上也是图像分割技术。1.2 课题研究的背景及意义医学图像分割技术是医学图像处理与分析领域的重要课题之一，也是近些年来备受研究人员高度关注的热点问题。医学图像分割是病变组织或者正常组织的三维可视化、图形引导手术、手术模拟等后续处理工作的基础，分割的准确性对诊断患者疾病的实际情况并给出相应的治疗方案是至关重要的。所以图像分割在医学领域中具有特殊的重要意义，特别是对人类生命健康息息相关皮肤病变图像中真皮区域图像分割研究成为医学图像可视化研究中的一个非常重要的部分，是科学计算可视化中研究的热点难题。医学图像分割2的目的是把图像分成

5、各具特性的区域并提取出感兴趣的部分，并使分割结果尽可能的接近解剖结构，从而为临床诊疗和病理学研究提供可靠依据。由于人体皮肤解剖结构的组织器官形状的不规则性及个体之间的差异性，一般将单一的图像分割方法应用于医学图像并不能得到理想分割效果，为此必须针对某一部分寻找一种特定有效的医学图像分割方法。1.3 医学图像分割的研究现状医学图像分割技术一直是医学图像领域研究的一个热点课题。如图1-2所示其发展过程经历了一个从人工分割到半自动分割到现在的自动分割。完全的人工分割方法是由医务工作者或者专家根据经验在原始图像上直接标出期望的边界2，这个方法的特点是工作量非常大，而且容易让入厌烦，容易出错。它的精确度

6、也不容易保证。半自动分割法是凭借医生的工作经验以及对图像知识的理解，通过一定的人机交互，由计算机完成分割，具有较快的分割速度和比较高的分割精度，与完全人工分割相比，半自动分割方法大大减少了人为因素的影响，而且分割速度快、分割精度高，但是操作者的知识和经验仍然是影响图像分割质量的一个重要因素。而且人机交互的半自动方法费时费力、可重复性低，这就要求有一种可重复性好、高效率、高精度的自动的图像分割方法。有鉴于此，一种可重复性好、效率高、精度高的自动图像分割方法已成为当前应用中的亟需。图1-2 图像分割技术发展过程自从上世纪六七十年代以来，人们研制开发出了众多的医学图像分割法，这些方法主要可分为基于边

7、界的分割法和基于区域的分割法两大类别。基于边界的分割方法主要是利用梯度信息确定目标的边界，包括局部微分算子，如 Rboert 算子，Sobel算子，Laplacian 二阶差分算子等。这些方法不依赖于已处理像素的结果，适合于并行化，但缺点是对噪声敏感，而且当边缘像素值变化不明显时，容易产生假边界或不连续的边界。基于区域的分割法依赖于图像空间的特征，如灰度、纹理及其它类型统计特征等，典型的基于区域的分割法有阈值分割、区域生长与分裂、分类器和聚类等。1.4 本文的主要工作及结构安排本文的主要研究内容如下：首先介绍图像分割的背景及发展现状，颜色模型。分析了图像分割典型分割算法和评估方法，为后续的算法

8、研究提供理论基础。其次概述了基于K-means 的医学图像分割算法，分析了算法原理及实现流程。再次，K-means分割算法进行Matlab仿真实验，和半自动分割算法live wire比较并对结果进行分析。 2. 颜色模型基本概念2.1 引言自然界中的光根据颜色可分为彩色光和非变化光，按频率成分可分为单色光和复合光，按频率和颜色综合考虑可分为谱色光（Spectrum Color）和非谱色光，按发光方式可分为直射光，反射光和透射光。无论是什么光，它的颜色都是取决于主观与客观两方面的因素。主观因素是因为人的眼睛有视觉特性，不同的人对于同一功率谱的光的色感可能是不相同的，光给人的颜色感觉与人眼的视觉特

9、性有关，客观因素是它的功率波谱分布，即光源的颜色，直接取决于它的功率谱，因此物体的颜色不仅取决于它的反射特性和透射特性，而且还与照射光源的功率谱有密切关系。2.2 颜色模型定义颜色模型3就是指某个三维颜色空间中的一个可见光子集，它包含某个颜色域的所有颜色。例如，RGB颜色模型就是三维直角坐标颜色系统的一个单位正方体。颜色模型的用途是在某个颜色域内方便的指定颜色，由于每一个颜色域都是可见光的子集，所以任何一个颜色模型都无法包含所有的可见光。在大多数的彩色图形显示设备一般都是使用红、绿、蓝三原色，我们的真实感图形学中的主要的颜色模型也是RGB模型，但是红、绿、蓝颜色模型用起来不太方便，它与直观的颜

10、色概念如色调、饱和度和亮度等没有直接的联系。2.3 颜色模型分类颜色模型主要有RGB、CMYK、LAB等。它们在不同的行业各有所指，但在计算机技术方面运用最为广泛。2.3.1 RGB颜色模型RGB（Red，Green，Blue）颜色模型通常使用于彩色阴极射线等彩色光栅图形显示设备中，彩色光栅图形的显示器都使用R、G、B数值来驱动R、G、B 电子枪发射电子，并分别激发荧光屏上的R、G、B三种颜色的荧光粉，发出不同亮度的光线，并通过相加混合产生各种颜色；扫描仪也是通过吸收原稿经反射或透射而发送来的光线中的R、G、B成分，并用它来表示原稿的颜色。RGB颜色模型称为与设备相关的颜色模型，RGB颜色模型

11、所覆盖的颜色域取决于显示设备荧光点的颜色特性，是与硬件相关的。它是我们使用最频繁，最熟悉的颜色模型。它采用三维直角坐标系。红、绿、蓝原色是加性原色，各个原色混合在一起可以产生复合色。可以用单位立方体来表示。在正方体的主对角线上，各原色的强度相等，产生由暗到明的白色，也就是不同的灰度值。（0，0，0）为黑色，（1，1，1）为白色。正方体的其他六个角点分别为红、黄、绿、青、蓝和品红。2.3.2 CMYK颜色模型CMYK（Cyan，Magenta，Yellow）颜色模型应用于印刷工业，印刷业通过青（C）、品（M）、黄（Y）三原色油墨的不同网点面积率的叠印来表现丰富多彩的颜色和阶调，这便是三原色的CM

12、Y颜色模型。实际印刷中，一般采用青（C）、品（M）、黄（Y）、黑（BK）四色印刷，在印刷的中间调至暗调增加黑版。当红绿蓝三原色被混合时，会产生白色，但是当混合蓝绿色、紫红色和黄色三原色时会产生黑色。既然实际用的墨水并不会产生纯正的颜色，黑色是包括在分开的颜色，而这模型称之为CMYK。CMYK颜色空间是和设备或者是印刷过程相关的，则工艺方法、油墨的特性、纸张的特性等，不同的条件有不同的印刷结果。所以CMYK颜色空间称为与设备有关的表色空间。而且CMYK具有多值性，也就是说对同一种具有相同绝对色度的颜色，在相同的印刷过程前提下，可以用分种CMYK数字组合来表示和印刷出来。这种特性给颜色管理带来了很

13、多麻烦，同样也给控制带来了很多的灵活性。在印刷过程中，必然要经过一个分色的过程，所谓分色就是将计算机中使用的RGB颜色转换成印刷使用的CMYK颜色。在转换过程中存在着两个复杂的问题，其一是这两个颜色模型在表现颜色的范围上不完全一样，RGB的色域较大而CMYK则较小，因此就要进行色域压缩；其二是这两个颜色都是和具体的设备相关的，颜色本身没有绝对性。因此就需要通过一个与设备无关的颜色模型来进行转换，即可以通过介绍的Lab色空间来进行转换。2.3.3 Lab颜色模型Lab颜色模型是由CIE（国际照明委员会）制定的一种色彩模式。自然界中任何一点的颜色都可以在Lab空间中表达出来，它的色彩空间比RGB空

14、间还要大。另外，这种模式是以数字化方式来描述人的视觉感应，与设备无关，所以它弥补了RGB必须依赖于设备色彩特性的不足。由于Lab的色彩空间要比RGB模式的色彩空间大。这就意味着RGB所能描述的色彩信息在Lab空间中都能得以影射。Lab颜色模型中L表示照度（ Luminance），a 表示从红色到绿色的范围， b 表示从蓝色到黄色的范围。在Lab 颜色模式下，图像的色彩信息和亮度信息被分开保存，调整颜色通道时亮度通道将保持不变。这样L通道可以看作是一影像的灰度版，其中保存了图像的细节信息，因此利用L通道可容易区分自然图像中的明暗细节。此外Lab模型具有宽阔的色域，不仅包含了RGB 的所有色域,

15、而且弥补了RGB色彩模型色彩分布不均的问题。2.4 颜色模型的选择和转换随着印刷行业从模拟到数字的变化，色彩的准确再现问题已经变得十分重要。因此我们需要更好、更快、更准确地获得彩色图像。上述三种颜色模型各有特点，而Lab颜色模型相比之下有以下优点4：1）Lab模式既不依赖光线，也不依赖于颜料。2）处理速度与RGB模式同样快，比CMYK模式快很多。3）Lab模式在转换成CMYK模式时色彩没有丢失或被替换。因此为了避免色彩损失的方法是：应用Lab模式编辑图像，再转换为CMYK模式打印输出。4）Lab模式还弥补了前面两种色彩模式的不足。RGB在蓝色与绿色之间的过渡色太多，绿色与红色之间的过渡色又太少

16、，CMYK模式在编辑处理图片的过程中损失的色彩则更多，而Lab模式在这些方面都有所补偿。在表达色彩范围上，处于第一位的是Lab模式，第二位的是RGB模式，第三位是CMYK。颜色特征是在图像分割中应用最为广泛的视觉特征, 在本文的医学图像分割算法中, 我们采用了Lab 色彩模型。颜色模型转换5是指将一个颜色模型中的颜色数据转换或表示成另一个颜色模型中的相应数据，即用不同的颜色模型中的数据表示同一颜色。在本文中要将与设备相关的RGB颜色模型转换到与设备无关的Lab颜色模型。颜色模型转换的方法有很多种，本文主要介绍三维查表插值法。三维查找表法是当今颜色模型转换较为常用的一种算法。三维查找表算法的基本

17、思想是，将源颜色模型进行分割，划分成一个个规则的立方体，而每个立方体的八个顶点的数据都是已知的，然后将所有源模型的已知点组成一张三维查找表。当给定源模型中任意一个点时，就能够通过查表找到与之相邻的八个数据点构成一个小立方体格子的节点，通过对这个小立方体的八个顶点进行插值，得到目标模型对应的数据。一般情况下都是将查找表法和插值法结合起来使用，使之成为带有插值法的三维查找表法，此方法步骤如图2-1所示：图2-1 颜色模型转换 1) 分割：将源色彩空间按一定的采样间隔分区，建立三维查找表；2) 查找：对于一个已知的输入点，搜索源空间，找出包含它的由八个栅格点构成的立方体；3) 插值：在一个立方体的栅

18、格内，计算出非栅格点上的颜色值。2.5本章小结本章主要介绍各类颜色模型，Lab色彩模型它不仅包含了RGB，CMYK的所有色域，还能表现它们不能表现的色彩。人的肉眼能感知的色彩，都能通过Lab模型表现出来。另外，Lab色彩模型的绝妙之处还在于它弥补了RGB色彩模型色彩分布不均的不足，因为RGB模型在蓝色到绿色之间的过渡色彩过多，而在绿色到红色之间又缺少黄色和其他色彩 。本文是针对医学真皮图像，我们想在数字图像的处理中保留尽量宽阔的色域和丰富和色彩，所以选择Lab颜色模型。3. 医学图像分割概述3.1 引言随着影像医学在临床医学的成功应用，图像分割在影像医学中发挥着越来越大的作用。如图3-1所示，

19、进行图像分割首先选择合适的颜色模型，形成特征向量空间；然后，在此特征空间中，应用提出的算法进行区域分割；最后抽取图像区域的特征。图像分割是提取影像图像中特殊组织的定量信息所不可缺少的手段，同时也是可视化实现的预处理步骤和前提。分割后的图像正被广泛应用于各种场合，如组织容积的定量分析，诊断，病变组织的定位，解剖结构的学习和计算机指导手术6。医学图像分割到今天仍然没有获得圆满解决，一个重要的原因是医学图像的复杂性和多样性。图3-1 图像分割过程框图3.1.1 医学图像分割特点随着数字图像技术的不断发展，医学图像分割方法也层出不穷，但仍然无法满足临床医学的实际要求。与一般的自然图像比较，医学图像有以

20、下特点7：1) 模糊性和不均匀性，主要来源于图像灰度上的含糊性，在同一组织中灰度值会出现大幅度的变化，即使密度值相同也会出现分布不均匀的现象。2) 不确定性，病态会导致出现一些正常组织或部位没有的结构。3) 多模态性，不同的图像模态反映了不同表现的图像。原因是不同的医学成像设备的成像原理各不相同。上述众多因素给医学图像分割带来了不可预知的难度，这就要求对医学图像的分割有较强的针对性。因此，对特定的医学图像分割任务提出切实可行的分割方法是当前医学图像分割邻域的主要策略。3.1.2 医学图像分割定义图像分割8是指将图像中具有特殊涵义的不同区域区分开来；这些区域是互相不交叉的，每一个区域都满足特定区

21、域的一致性。定义将一幅图像，其中 x ， y ,进行分割就是将图像划分为满足如下条件的子区域.：(1)，即所有子区域组成了整幅图像；(2)是连通的区域；(3)（k=1，2，.，N）；J=1，2，N；k7)，即任意两个子区域不存在公共元素；(4) 区域 m 满足一定的均一性条件，均一性(或相似性)一般指同一区域内的像素点之间的灰度值差异较小或灰度值的变化较缓慢。如果连通性的约束被取消,那么对像素集的划分就称为分类（pixel classification），每一个像素集称为类（class）。在下面的叙述中，为了简单，我们将经典的分割和像素分类通称为分割。3.1.3 图像分割算法评估标准医学图像分

22、割质量的高低直接影响着临床的诊断结果；分割速度的快慢直接影响着分割算法是否符合临床诊断对实时性的要求，从而能为临床医疗工作者提供需要的实时交互。医学图像分割的评估标准的种类非常多，算法评价9对指导图像分割任务意义重大。通过对图像分割算法的对比评价，可以针对特定的分割任务选择最优的分割算法；通过对图像分割结果进行评价，可以进一步分析所选算法的有效性，以及评判新的图像分割算法的分割性能。一般来说，对图像分割的结果以人的主观判断作为评价准则，尽管对大多数图像处理问题而言，最后的住宿是人的视觉，但对不同图像分割方法的处理结果作一定量的比较、评价也是必需的，如何对现有的图像分割方法做一个全面系统的比较，

23、并给出公正的评价，这是一个有意义但是比较困难的工作。到目前为止，已经有研究专家在这方面做出了一些工作，他们认为对图像分割方法的评价可以通过分析和实验两种方式来进行。分析法是直接分析分割算法本身的原理及性能，在实践中很少得到应用。实验法是利用由分割算法得到的实际分割图像与参考图像进行比较，根据实际分割图像与参考图像的符合程度来评价分割算法的优劣，这是一个比较实际的方法。实验法已经提出了不少，目前定性实验准则以及常用的定量实验准则；根据临床中的需求评估标准可以分为以下两类：1) 客观评价图像分割结束后算法花费的时间，算法的鲁棒性等等关键性指标都是可以得到的，可以通过定性或者定量的实验给予评价，使用

24、重合率以及错分率来评价评价图像分割的精度。2）目标检测目标检测是主观评价方法，可以使用一些相关领域的专家和学者进行目测检验。这种评估方法也是很重要的。因为做分割的最终目标就是得到他们认可的分割效果，从而在临床中使用。3.2 医学图像分割典型算法图像分割算法的研究一直受到人们的高度重视，到目前为止，提出的分割算法已经多达上千种，由于现有的分割算法非常多，所以将它们进行分类的方法也提出了不少。图像分割通常利用不同对象之间特征的不连续性，以及同一对象内部的特征相似性。而基于区域的算法侧重于利用区域内特征的相似性。例如有把分割算法10分成3类的：区域生长和分裂合并法、阈值分割法、分类器和聚类。本章从实

25、际应用的角度考虑，详细介绍了图像分割的如下算法：区域生长和分裂合并法、阈值分割法、分类器和聚类等。3.2.1 区域生长和分裂合并法基于区域的图像分割算法11是以区域为对象依照图像纹理特征、灰度值、换图像的特征等共同属性来划分图像区域。该方法给图像中的每一个像素都赋予一个类别属性，使具有预定义属性的像素聚集在同类别中。最常见的基于区域分割法是区域生长和分裂合并法。区域生长的基本思想是将具有相似特性的像素合并到某个区域。首先为每个需要分割的区域确定一个种子像素作为生长起点（见图3-2），然后按一定的生长准则把它周围与其特性相同或相似的像素合并到种子像素所在的区域中，把这些新像素作为种子继续生长（见

26、图3-3），直到没有满足条件的像素可被包括，这时一个区域就形成了。区域生长法所要解决的主要问题是区域生长准则的设计。生长准则往往和具体问题有关，直接影响最后形成的区域，如果选取不当，就会造成过分割或欠分割的现象。区域生长法的优点是计算简单，对于较均匀的连通目标有较好的分割效果。它的缺点是需要人为确定种子点，对噪声敏感，可能导致区域内有卒洞；另外，它是一种串行算法，当图像目标较大时，分割速度较慢，因此在设计算法时，要考虑尽量提高效率。 图3-2 设定种子点 图3-3 生长过程分裂合并法的基本思想是从整幅图像开始通过小断分裂合并得到各个邻域。这种方法对较复杂图像的分割，效果比较好，但算法较复杂、计

27、算量大，同时分裂还有可能破坏区域的边界。3.2.2 阈值分割法阈值法11是一种较常见的、并行的、直接检测区域的分割方法。该方法是基于对灰度图像的一种假设，即目标或背景内相邻像素间的灰度值是相似的，不同的目标或背景的像素在灰度上有差异的，反映在直方图上，不同的目标或背景对应不同的峰值。阈值是指用于区分不同目标图像的灰度值。如果图像中只有目标和背景两类，那么只需要选取一个阈值，称为单阈值分割（如图3-4），这种方法是将图像中每个像素的灰度值与该阈值进行比较，灰度值大于该阈值的像素为一类，而小于阈值的像素为另一类。如果图像中有多个目标，就需要选取多个阈值将每个目标分开，这种方法称为多阈值分割（如图3

28、-5）。阈值法基于对灰度图像的一种假设：目标或背景内的相邻像素间的灰度值是相似的，但不同目标或背景的像素在灰度上有差异，反映在图像直方图上，不同目标和背景则对应不同的峰。选取的阈值应位于两个峰之间的谷，从而将各个峰分开。阈值分割的优点有计算简单、算法效率较高、速度快等，当不同类物体的灰度值或其它特征值相差很大时，能有效地对图像进行分割。该方法经常用于医学图像中目标区域和背景区域的分割，如骨骼和皮肤等不同组织的分割。论文就是针对真皮区域图像的分割法。 图3-4 单阈值分割 图3-5 多阈值分割3.2.3 分类器和聚类分类2是模式识别领域中一种基本的统计分析方法。分类的目的是利用已知的训练样本集在

29、图像的特征空间找到点（一维）、曲线（二维）、曲面（三维）或超曲而（高维），从而实现对图像的划分。用分类器进行分割是一种有监督的统计方法，它需要用手工分割得到的样本集作为对新的图像进行自动分割的参考。分类器有两个优点：1)不需要迭代运算，故计算量相对较小；2)能应用于多通道图像。但是分类器同样没有考虑空间信息，因此对灰度不均匀的图像分割效果不好，分类器还要求由手工分类生成训练集，而手工分类的工作量很大，同时用小量的训练集训练的分类器对大量的样本空间进行分类时会产生误差。聚类法是将图像空间中的像素用对应的特征空间点表示，通过将相似的特征点聚集成团，然后再将它们映射到原图像空间，得到分割结果。而图像

30、分割可看成是目标对象的分类问题，可以借助于模式识别中的模式分类技术。聚类算法与分类器算法极为类似，只是它不需要训练样本，因此聚类是一种无监督的（unsupervised）统计方法，聚类算法迭代地执行图像分类和提取各类的特征值，是一种自我训练的分类。聚类算法包括K均值法、模糊C均值（FCM）法。K均值算法先对当前的每一类求均值，然后按均值对像素进行重新分类（将像素归入均值最近的类），对新生成的类再迭代执行前面的步骤。聚类分析不需要训练集，但是聚类时需要有一个初始分割提供初始参数，即聚类的数目，初始参数对最终分类结果影响较大。此外，聚类也没有考虑空间关联信息，因此也对噪声和灰度不均匀敏感。3.3

31、本章小结 本章对医学应用领域中一些有代表性的图像分割方法做了一个比较全面的综述，并分析了几类方法的特点和局限性。近几年来虽然仍然有很多研究人员致力于医学图像分割的研究，发表了很多的研究成果，但由于问题本身的困难性，目前的方法多是针对某个具体任务而言的，还没有一个通用的解决办法。4. 医学病理切片中真皮组织区域分割算法4.1引言图像分割即图像分类，图像分类学是人类认识世界的基础科学。聚类分析和判别分析是研究事物分类的基本方法，广泛地应用于自然科学、社会科学、工农业生产的各个领域。本文是对医学病理切片中的真皮组织区域进行分割，若事先已经建立类别，则使用线性鉴别分析法（Linear Discrimi

32、nant Analysis ，LDA）；若事先没有建立类别，则使用K-means聚类分析。4.2 基于K-means医学图像分割算法K-means12分割算法在图像分割中得到广泛应用，属于无监督分类的一个分支。其主要目的是对具有相同数据类型的样本数据按距离最短规则进行集合的划分,最终获取各等价类。K-means算法是典型的基于原型的目标函数聚类方法的代表，它是数据点到原型的某种距离作为优化的目标函数，利用函数求极值的方法得到迭代运算的调整规则。K-means算法以欧式距离作为相似度测度，它是求对应某一初始聚类中心向量V最优分类，使得评价指标J最小。算法采用误差平方和准则函数作为聚类准则函数。4

33、.2.1 K-means算法简介假设有一组包含K个聚类的数据,其中第k个聚类可以用集合来表示,假设包含个资料（，. ），此聚类中心为，则该聚类的平方差可以定义为: (4-1) 其中是属于第k类的资料点。而这K个聚类的总和平方差E便是每个聚类的平方差总和： (4-2)K-means算法聚类其算法步骤如图4-1所示，通过迭代的方式, 设法降低E 的值，以使得各聚类本身尽可能的紧凑，而各聚类之间尽可能的分开。图4-1算法基本流程4.2.2 仿真结果本文随机选取医学真皮区域切片图作为实验的材料，在实现彩色图像聚类分割后进一步基于纹理特征实现区域描述，根据真皮切片选定聚类的种类k=3,如图4-2，4-3

34、中给出了两组实验结果。其中图(a)是原始真皮病理切片图像,可以看出该图像可以分为白色背景区域，紫色部分为真皮区域，其他浅粉色部分为组织区域，图（b）是基于空间、颜色和纹理特征进行聚类分割的结果，灰色部分为背景区域，黑色区域为真皮区域，白色区域为组织区域。为了标记分割区域，分区域显示，图（c)的紫色部分为真皮区域，也是本文的感兴趣区域。图（d）中的浅粉色区域为组织区域，图（e）黑色部分为去除背景的区域。 （a)原始图像 (b)聚类结果图 (c)分割区域一 （d）分割区域二(e)分割区域三图4-2 仿真实验组1 （a)原始图像 (b)聚类结果图 (c)分割区域一 （d）分割区域二(e)分割区域三图

35、4-3 仿真实验组2实验表明，不同的医学病理图片，分割的效果不同，聚类 K 的选择对K-means算法分类结果影响很大，若K取小了就会有的区域不能分割出来，若K取大了就会出现过分割。所以在分割图像前K值的大小一般都由先验知识判断得出，当需要处理大量图像时，由人为判断给出K值是不现实的。K-means算法计算复杂度低，仍作为一种广泛应用的分割算法所分割的结果仍无法满足有待今后进一步研究和改进。4.3 算法评估与分析K-means聚类法在图像分割中得到广泛应用，基于K-means聚类法进行图像分割是通过找到特征空间中像素值的空间聚类，并把每一个像素划分到不同的聚类中来实现图像的分割。作为一种经典的

36、聚类算法，K-means算法有其优越性, 当然也有其局限性，如聚类数K 值必须依据先验知识提前设定、最初的类中心是随机选取的（这常常导致每次聚类产生不同的分类结果）、聚类效果的好坏依赖于距离判定公式、图像特征值等条件的选取。为了比较分析分类算法的优缺点，在这里作为对比选用了由Artur Chodorowskia提出的Color Lesion Boundary Detection Using Live Wire13 的半自动分割算法进行了仿真实验。如图4-7示该文献所提出的live wire算法主要应用于口腔病变区域的划定，用live wire算法进行了检测和比较，将其用于医学病理真皮区域图像边

37、界划定。其算法思想是把包括括色渐变分析大小，颜色渐变方向，Canny边缘检测，拉普拉斯算子和零交叉等多重成本方面纳入最后的总成本函数，使它们直接作用于三个分量的彩色图像，而不是使用单一的派生色带。live wire算法被证明是相比于自动分割更精确。用于评估K-means聚类算法图4-7 口腔病变区域划定如图4-8，4-9示对仿真实验组1分别用K-means和live wire分割后的真皮区域分进行位置标注，并坐标归一化为500800mm。 （a）K-means区域标注 (d) K-means坐标标注 (c) 评估live wire区域标注 （f）live wire坐标标注 图4-8 分割区域标

38、注 图4-9 区域坐标归一如图4-10，4-11所示对仿真实验组2分别用K-means和标准live wire分割后的真皮区域进行位置信息标注，并坐标归一550500mm。 （a）K-means区域标注 (d) K-means坐标标注 (c) 评估live wire区域标注 （f）live wire坐标标注 图4-10 分割区域标注 图4-11 区域坐标归一评估基本思想是如式（4-8）所示，通过计算分割区域的重合比例评估K-means 和LDA算法精度。 （4-8)其中A分别代表K-means算法分割区域，live wire算法分割区域 。B代表医生手动分割区域，通过实验组1，实验组2仿真实验

39、结果如表5所示：表5 算法评估算法实验组1实验组2K-means0.65850.7182Live wire0.68710.7664实验数据如表5所示，可以看出在同一实验组live wire算法相比K-means算法具有更高的精确度，K-means分割算法是一种无监督的学习方法，从K-means聚类的算法看，K-means聚类是个持续迭代的过程，影响聚类结果的主要因素就是划分的类数K、初始聚类中心C 的选择和对象间的距离d（X ，C )。实际上从人眼的视觉特性出发，人眼能够直接分辨的色阶一般不会超过60 个，对医学彩色图像分割处理而言，分类的类数会远远小于此数，实际类数的选择可以依赖具体的应用环

40、境和图像分割的目的人为确定。由于K-means聚类采用目标函数最小化作为聚类结束的条件，并且初始聚类中心采用随机选择,，所以比较容易陷入局部最优，K-means对初始聚类中心有很强的敏感性，不同的初始聚类中心可能导致完全不同的聚类效果。K-means聚类算法简洁，具有很强的搜索力，适合处理数据量大的情况，在数据挖掘领域中得到广泛的应用。而live wire是一种半自动分割算法，算法精度相比自动分割算法在特殊的人体真皮组织要高。4.4 本章小结本章主要通过在对两种分割算法原理的深入学习后，在Matlab平台下仿真，将一种应用于口腔病变区域划定的live wire半自动分割算法用于实现医学病理切片

41、真皮区域分割，将其作为对比K-means，对于两种分割算法，首先运用分割算法进行分割并提取出感兴趣区域，对区域进行描述。利用重合率参数fraction，对比K-means和live wire算法，实验取得了比较好的效果。5. 总结与展望5.1 总结图像分割是图像处理技术中的关键，对于图像分割的研究，现已提出了数以千计的不同算法。虽然这些算法都在不同程度上取得了一定的分割效果，但图像分割研究仍然面临着极大的挑战。至今，不但没有一个算法能够对所有图像都可以给出很好的分割结果，而且没有一个方法能够对由不同算法得到的分割结果给出与主观感觉一致的评价，也没有一个统一的理论标准来指导如何根据不同的图像选择

42、合适的分割算法。所以，在解决一些实际的图像分割问题时，往往还是依赖于自己的知识和经验，这必将限制了图像分析和理解研究的发展。医学图像由于其与人类生命健康息息相关，对其进行图像分割研究更加引起了人们的广泛关注。医学图像分割是人体组织及器官的三维可视化、手术模拟、图形引导手术等后续处理的基础，分割的准确性对医生判断患者疾病的真实情况并做出相应的诊断计划至关重要。因此，图像分割在医学应用中具有特殊的重要意义。虽然对于复杂的医学图像分割，目前已有学者研究出不少边缘检测、区域分割及目标提取的方法，但也没有一种普遍适用图像分割方法，只能针对某一具体应用而采用相对行之有效的算法。由于人体组织及其分割技术的复

43、杂性，图像分割成为医学图像处理技术发展的瓶颈。本文针对医学病理真皮区域图像，对图像分割问题进行了深入研究及详细讨论，而且进行了相应的分割算法仿真并给出实验结果分析，并对下一步工作做出展望。5.2 展望通常将医学图像分割面向具体的临床应用，因此图像分割的准确性更为重要。图像分割一直是一个很困难的问题，目前的自动分割方法虽然在一些方面取得了一定的成功，但还远远不能满足医学图像处理的实践中对分割结果准确性的要求因而，近年来由用户参与控制、引导的交互式分割方法在医学图像分割中正受到越来越多的关注。目前人们仍在继续研究更先进的成像技术和更复杂的图像处理算法图像分割方法的研究与分析与影像数据的物理成像原理

44、、图像形成和重构算法的关系更为密切，而且图像分割与其它图像处理分析任务（ 如图像增强、匹配、可视化）在识别对象结果和功能上是相关的，因此将它们结合起来若同研究是未来研究的一种趋势。图像分割结束后算法花费的时间，算法的鲁棒性等等关键性指标都是可以通过定性或者定量的实验给予评价，由于时间的关系本文只针对分割评价指标中的重合率对K-means和live wire进行评价，在以后的研究中将继续拓展利用错分率来评价评价图像分割的精度。医学图像的分割算法的验证和评估方法仍需进一步研究。总之，医学图像的分割算法的研究道路曲折而充满挑战。第 22 页 参考文献1 Zhang Yu Jing. Image Se

45、gmentation M, Beijing, China: Science Publisher, 2001, PP. 2-3.2 刘茂奇. 基于人体切片的图像分割技术研究 D, 广州: 广东工业大学, 2006.3 李新友, 张云松. 颜色模型浅析J, 理论探讨, 8（27), 2010: 132-133.4 王易, 赵勋杰. 基于K均值聚类分割彩色图像算法的改进J, 计算机应用与软件, 8（27）, 2010: 127-130. 5 古晶. RGB到Lab颜色空间转换N, 印刷技术, 2009: 11-13 .6 Duncan J S, Ayache N. Medical image Ana

46、lysis: Progress over Two Decades and the Challenges Abead IEEE Trans on Pattern Analysis and Machine Intelligence. 2010. 22(1): 181204.7 黄琴波. 医学图像分割算法研究及其在骨分割中的应用D, 西安: 西安电子科技大学, 2011. 8 Gonzalez R C, Wood R E. Digital Image Processing.3nJ.eds Addison wesley , 1992，40-48.9 韦明祥. 基于水平集方法的医学图像分割算法研究D,

47、西安: 西安电子科技大学, 2012. 10 冈萨雷斯 等著, 阮秋琦 等译. 数字图像处理的MATLAB实现（第二版）M, 北京: 清华大学出版社, 2013. 11 黄文博, 燕杨, 王云吉. 医学图像分割方法综述J, 长春师范学院学报（自然科学版）, 2013, 32(2): 22-25.12 李翠, 冯冬青. 基于改进k-均值聚类的图像分割算法研究D, 郑州大学学报, 1（43), 2011: 109-113.13 Artur Chodorowski , Ulf Mattsson . Color Lesion Boundary Detection Using Live Wire Chalmers University of Technology Central Hospital, Karlstad, Sweden.第 23 页