模式识别报告

1、基于CV-SVM的医学图像分割摘要：随着计算机成像技术的不断完善，医学图像在疾病判别中发挥着越来越 重要的作用。近年来，针对医学图像尤其是脑部医学图像的研究，成为一门新兴 学科。核磁共振脑图像中的肿瘤识别及提取是医学图像处理中的难点，本文主要 使用改进的支持向量机方法研究肿瘤化组织的自动检测和提取，以提高工作效率 和疾病诊断的准确率。受脑图像复杂程度的限制，单独应用支持向量机和CV模 型水平集方法都无法准确有效的将脑肿瘤分割出来，本文提出了一种CV模型 和支持向量机结合的CV-SVM方法。该方法在脑部图像四分类(背景、脑脊液、 白质和灰质)中具有良好的分割效果，能够在相似灰度区域中提取出肿瘤化

2、组织 图像特征，便于进行脑瘤医学诊断，具有应用价值。关键字:支持向量机、关键字:支持向量机、CV、医学图像分割Abstract: This work is presentKeywords： SVM、CV、medical image segmentation研究背景与目的近20年来,医学影像已成为医学技术中发展最快的领域之一，其结果使临床医 生对人体内部病变部位的观察更直接、更清晰，确诊率也更高。20世纪70年代初,X -CT的发明曾引发了医学影像领域的一场革命，与此同时,核磁共振成像 (MRI:Magnetic Resonance Imaging)超声成像、数字射线照相术、发射型计算机 成像和

3、核素成像等也逐步发展。计算机和医学图像处理技术作为这些成像技术的 发展基础,带动着现代医学诊断正产生着深刻的变革。各种新的医学成像方法的 临床应用,使医学诊断和治疗技术取得了很大的进展,同时将各种成像技术得到的 信息进行互补，也为临床诊断及生物医学研究提供了有力的科学依据。在目前的影像医疗诊断中,主要是通过观察一组二维切片图象去发现病变体, 往往需要借助医生的经验来判定。至于准确地确定病变体的空间位置、大小、几 何形状及与周围生物组织的空间关系，仅通过观察二维切片图象是很难实现的。 因此,利用计算机图象处理技术对二维切片图象进行分析和处理,实现对人体器 官、软组织和病变体的分割提取、三维重建和

4、三维显示，可以辅助医生对病变体 及其它感兴趣的区域进行定性甚至定量的分析，可以大大提高医疗诊断的准确性 和可靠性。此外，它在医疗教学、手术规划、手术仿真及各种医学研究中也能起 重要的辅助作用。图1三维可视化大体步骤医学图像的三维可视化的方法很多，但基本步骤大体相同，如图1。从CTM/ R 或超声等成像系统获得二维断层图像,然后需要将图像格式转化成计算机方便处 理的格式。通过二维滤波,减少图像的噪声影响，提高信噪比和消除图像的尾迹。 采取图像插值方法,对医学关键部位进行各向同性处理，获得体数据。经过三维滤 波后,不同组织器官需要进行分割和归类，对同一部位的不同图像进行配准和融合, 以利于进一步对

5、某感兴趣部位的操作。根据不同的三维可视化要求和系统平台的 能力，选择不同的方法进行三维体绘制，实现三维重构。最重要的环节便是医学图像分割。理想的图像分割技术能将医学图像中各个 部位和组织完整的分离开来，使得医学影像的成像结果更加直观，也使病变部位 的识别判断更加准确。现有的图像分割技术，无论是单一的分割方法，还是基于 单一的特征信息进行分割的方法，在对分割精确度和效率要求较高的临床医学领 域都难以取得令人满意的分割效果，因而人们开始寻找新的概念和方法并将其引 入图像分割领域。近几年来提出了几种新方法，其中基于多特征融合和多谱图像 分析相结合的方法应用最为广泛，然而，这些方法多数是基于渐近假设理

6、论发展 而来的，且都默认样本数目趋于无穷大，但实际问题中样本数目往往是有限的， 所以它们在很多方面都出现了问题，有些推广能力较差，有些学习速度较慢，还 有些难以收敛，在高维度特征和小样本学习方面都难以取得良好的效果，无法进 行分割。为了解决上述问题，支持向量机(Support Vector Machine，SVM)方法 应运而生。支持向量机理论最初就是基于样本线性可分时的最优分类超平面理论发展 起来的，它的主要目的是寻找一个规则能把Rn空间内的样本点分成两部分，这 个规则就是能使两类样本离超平面的距离最大化的最优分类超平面。最优分类面 的目的就是能将两类样本正确分开且使分类间隔最大，这实际上也

7、是对推广能力 的控制，是SVM的核心思想之一。SVM是以统计学习理论为基础发展起来的 一种新的模式识别分类方法，依据Vapnic结构风险最小化原则和VC维理论 在小样本、高维度、非线性和局部极小点数据空间中获得了良好的推广能力，被 看作是对传统分类器的发展。尽管支持向量机在图像分割方面获得了较好的分类效果，但在图像中提取的 分类样本中含有噪声或野值样本时，可能在构造最优分类面时出现偏差，得到的 结果并不是真正的最优分类面，导致分类误差较大而失去作用。近年来针对医学 图像的特点人们又提出了许多关于SVM算法的改进和应用研究。经典SVM 方法通过参数C对目标函数中的分类误差项进行惩罚，因而参数C的

8、选择对 图像分类结果至关重要，针对此问题，Scholkoph等人提出了v-SVM方法。针 对SVM方法对噪声和野值比较敏感的问题，Suykens等人提出了最小二乘支持 向量机(LS-SVM)，简化了运算复杂性，却由于每个样本对分类器的贡献失去了 稀疏性的优点。后来Lin等人提出了模糊支持向量机方法(Fuzzy-SVM)，这种方 法针对不同的样本数据采用不同的惩罚系数，使不同的样本在目标函数中的贡献 不同，并且使含有噪声或野值的样本的权值较小，达到消除噪声或野值影响的目 的。医学图像分割中经常会出现不同类样本数量差异较大的情况，使用SVM经 常会出现分类结果偏向于样本数较多的类的问题，而医学图像

9、中某些重要的类 (如肿瘤部位)和其他部位相比样本数量比较小，但又对分类精度要求较高，这 就需要更多新的解决方法。另外一种比较受公众关注的图像分割方法进入人类的视野一一基于水平集 方法求解主动轮廓模型的曲线演化方法。这种方法的原理是利用轮廓曲线的几何 特性建立能够驱动轮廓曲线演化的能量泛函，然后通过最小化能量泛函驱动轮廓 曲线逐渐逼近图像中的目标边界得到演化方程，最后利用水平集函数将轮廓曲线 的演化方程转化为求解数值化偏微分方程的问题。CV-SVM新方法针对医学图像的复杂多样性对支持向量机算法造成的运算 速度较慢，识别时间较长的影响，在支持向量机进行特征提取后，使用PCA方 法对特征向量进行降维

10、运算，提取有价值的主成分进行分析，使运算量大大下降， 识别过程得到简化，提高了运算速度。再通过两种方法的结合有效地减弱了各自 的缺点对分割结果的影响，使医学图像的分割结果在相似灰度及噪声影响下，细 节方面表现突出，在肿瘤的自动提取过程中也达到了较为满意的结果。本文提出利用CV模型和支持向量机相结合的CV-SVM新方法，分别对经 典脑图像和医院获取的脑部医学肿瘤图像进行分割和肿瘤提取，通过Matlab平 台给出仿真实验结果，通过实验数据和结果分析了新方法对图像分割和提取的影 响，证明了利用新方法对医学图像处理的可行性和有效性。方案设计针对MRI图像，支持向量机方法对噪声和野值比较敏感，所以在图像

11、分割 时容易出现边缘分类不清的问题；若图像中出现不同类的样本数目差异较大时， 分类结果还会向数目较多的类偏移，造成分类偏向问题。而CV模型水平集方法 则需要对整个待分割目标图像中所有点的水平集函数值进行计算和更新，工作效 率比较低；且在受到偏移场影响后的图像中，演化曲线无法演化到相应的边界上 造成较大的分割误差甚至分割错误。这些缺点一旦在医学图像分割领域得到体现 就会造成严重后果。由于支持向量机对噪声和野值的敏感度较高，CV模型受偏 移场影响较大，在较复杂的医学图像领域都不能独自胜任，本文提出利用CV模 型和支持向量机相结合的CV-SVM新方法。图2基于CV-SVN的医学图像分割流程图2给出了

12、 CV-SVM的流程图。利用CV模型的水平集对噪声具有较强 的鲁棒性先对图像进行粗略的分割，使用图像的边缘轮廓提取作为初始轮廓线， 进行迭代后将图像大致分割为不同的部分，然后使用支持向量机方法对提取的分 类的内部区域进行细节修正，修正过程使用分类后的目标区域进行特征提取，使 用主成分分析法对提取出的特征进行降维，使用简化的主成分学习后并细致的分 类，最后得到较为准确的目标区域。算法描述(一)cv模型水平集粗略分割设Q为图像区域，Rt)为演化曲线，水平集函数定义为：(P (x, y,t) 0 ,如果(x,y)在Rt)内部(p(x,y,t)=O,如果(x,y)在Rt)上(p (x,y,t) 0 ；

13、点 X在 C(O)外部，(P(X,t = 0) 0。初 始状态t=0时，有c G = o) = * | cp (x,t = o) = o水平集方法的主要目的是为水平集函数中G, t)构造一个能使它在演化界处的值为零的方程。设X(t)为演化界面上点的演化路径，由于函数在演化界面处的值为零， 即(P(X (t), t) = 0根据符合函数求导法则，得 TOC o 1-5 h z 甲 + 甲x = 0tx ii=l i t式中X为X的第i个分量。定义X = F(X(t)为垂直于演化界面的速度，因为 it5L xx i . i t 1=1 , cp , (p )(x , x , ., x ) = F(

14、X(t) 5L xx i . i t 1=1X! X2XN 12N用I V(p I代表函数(p(x,t)的梯度，F代表(p(x,t)的速度函数。水平集函数(p(x,t)的演化方程就可以写成多相位水平集分割n个复杂程度不同的拓扑结构相位只需要10气n个水平集函数，且 能够很好地解决空白与重叠问题。针对文中的四分类问题我们只需考虑m = 1og2 n(m=4) 个单相位水平集函数q.，则甲的所有零水平集集合构成的分割边界线就是多目标分割的分 界线。我们用C代表演化的曲线，用C1、C2分别代表曲线C内部区域和外部区域的均值。假 设某R 1均值。假设某图R0由两块近似常量区域组成，其像素值分别为片、牟

15、。其中要检测的目标区域由是像素值片，其边界为C。那么R.、*分别代表曲线内外的区域了。令F(C) = F(C) + F (C) = j | u - c |2dxdy + j | u - c |2dxdyinsede(C )outside(C )那么目标的边界就是上式中的能量函数达到最小时的边界，即F(C) = F(C) + F(C) = 0(二)特征提取和归一化在(1)中使用图像的边缘轮廓提取作为初始轮廓线，进行迭代后将图像大 致分割为不同的部分。然后使用支持向量机方法对提取的分类的内部区域进行细 节修正，具体如下：针对医学MR图像的复杂特点主要采用纹理特征与灰度特征(灰度共生矩阵)组合的 特

16、征提取方法来进行处理。角二阶矩:角二阶矩:f = E排1.=EZ (i-j)2P(i,j)j)P(i,j)-uuf =二3b b对比度:相关性:反差距:f =提54. . 1 + (i j )2i J和均值：f5 = 2np (i)i = 2i =2差均值：f6 =为iP (1)i =2方差：f7 = (口u)2p()和方差：f8 =云(i - f5) P (i)i = 2W12差方差：f9 =也(i - ?七y(i)i = 0上述几种特征能够反映图像灰度的变化特性。和熵：fi0 =-翅七(i)log(P (i)i = 2 熵：fii = -EEP(i,j)iog(P(i,j) 差熵：f* =

17、 -E P (i)log(P (i)i = 0特征到的各种熵则主要用来测量共生矩阵的系数变化复杂度和纹理变化复杂度。脑部医学MRI图像的窗口选择为7*7大小。在纹理特征提取时使用灰度共生矩阵在四个方向的12种纹理特征共48种特征，再加 上象素灰度，7*7窗口内的灰度中值及平均值等特征共51种特征作为图像的特征矢量。为避免特征的偏移，同时在计算特征向量的内积时避免大数计算引发的计算溢出的困 难，对数据进行归一化，采用如下公式归一化后特征特征向量-特征向量最小值归化后特征特征向量最大值-特征向量最小值(三)主成分分析归一化后的特征向量的范围被限制在0,1之间，简化了运算过程，提高了运算速度。采 用

18、主成分分析法对繁杂的特征向量进行提取，以降低分类代价，提高分类性能。PCA为了获得数目和计算量最小的主成分需要从信息压缩最大方向寻求图像的低维信 息表达，所以对应的K-L变换的信息熵最小，获得的识别空间也是对原图像的最优逼近。步骤如下：(1)对各维数据进行无量纲化；x x设x.j e X,i = 1,2,N,那么得到yij = -j G Yj其中xj = M E xkj是各维特征向量的均值，Sj = JmJ(xij - xj)2是各k=1i =1维特征向量的标准差。(2)计算Y的协方差矩阵S:n = Nn = NE yik y =ItJ ,Tl = |1,1,，1|1X NS = ；Y - Y

19、l X Y - YlT(3)解特征值和特征向量根据特征方程QI - S)U = 0求解S的特征值人和特征向量u。而后将人从大到小排列：七2.法2七,并找到与之对应的舟:u = u , u ，,u ,j = 1,2, , Nj 1j 2jNj计算主成分Z = y u + y u + + y uji1 1ji2 2jiN Nj这里 j = 1,2,p N选取主成分P的取值原则要让主成分满足最大程度包含原信息量最小值的条件。这里依据主成分贡 献率。a = .I丈.选取主成分。i =1选用Gauss径向基核函数(高斯核函数)K(x,x，) = exp(卜x2 / 2。2)，惩罚因子C取1000，高斯核

20、宽度参数。2取1，窗口大小为7*7。单-独使用SVM分割后的图单-独使用SVM分割后的图像使用CV-SVM分割后的图像图3 SVM与CV-SVM对比由图3中可以看出，CV-SVM的分割结果比经典支持向量机在细节和边缘处的优势更加 明显，且对噪声具有较强的鲁棒性。实验证明，CV-SVM可以应用在复杂度高，精确度要 求也较高的MRI脑部医学图像的分割领域。实验结果为了测试基于CV-SVM的医学图像分割方法的效果，采用青岛大学医学院附 属医院的带肿瘤组织的脑部图像作为实验对象。下面给出CV-SVM分类后提取 出的肿瘤化区域和专家手工标注的肿瘤化区域对比图。手工标出肿瘤部位图4专家手工标注肿瘤区域的脑

21、补图像手工提取肿瘤部位 SVM提取肿瘤部住CV-SVM提职肿瘤部位图5三种方法提取出的肿瘤区域对比图表1三种方法分割效果评价方法评价标准-手工分割使用经典SVM使用CV-SVM像素个数233320422356时间（S）7.86.735.96与手工分割对比漏检率0.14660.03129误检率-一0.021860.02143从上图5中我们可以直观的看到三种方法提取出的肿瘤组织的形状差异，由 于肿瘤疾病的复杂性和技术条件的局限性，医学图像分割的模板或者说准则还没 有形成，我们暂时使用医学专家手工提取的肿瘤组织区域作为比对模板。由表1 可以看出经典SVM和CV-SVM方法与手工分割方法相比，CV-SVM的像素个数 跟接近专家给出方案中的像素个数，而且误检率、漏检率也比经典SVM方法要 低，CV-SVM方