模式识别报告基于CVSVM的医学图像分割.docx
《模式识别报告基于CVSVM的医学图像分割.docx》由会员分享，可在线阅读，更多相关《模式识别报告基于CVSVM的医学图像分割.docx（14页珍藏版）》请在冰豆网上搜索。
模式识别报告基于CVSVM的医学图像分割
基于CV-SVM的医学图像分割
摘要：
随着计算机成像技术的不断完善，医学图像在疾病判别中发挥着越来越重要的作用。
近年来，针对医学图像尤其是脑部医学图像的研究，成为一门新兴学科。
核磁共振脑图像中的肿瘤识别及提取是医学图像处理中的难点，本文主要使用改进的支持向量机方法研究肿瘤化组织的自动检测和提取，以提高工作效率和疾病诊断的准确率。
受脑图像复杂程度的限制，单独应用支持向量机和cv模型水平集方法都无法准确有效的将脑肿瘤分割出来，本文提出了一种CV模型和支持向量机结合的CV-SVM方法。
该方法在脑部图像四分类（背景、脑脊液、白质和灰质）中具有良好的分割效果，能够在相似灰度区域中提取出肿瘤化组织图像特征，便于进行脑瘤医学诊断，具有应用价值。
关键字：
支持向量机、CV、医学图像分割
Abstract:
Thisworkispresent
Keywords:
SVM、CV、medicalimagesegmentation
1.研究背景与目的
近20年来,医学影像已成为医学技术中发展最快的领域之一，其结果使临床医生对人体内部病变部位的观察更直接、更清晰，确诊率也更高。
20世纪70年代初,X-CT的发明曾引发了医学影像领域的一场革命，与此同时,核磁共振成像（MRI：
MagneticResonanceImaging）、超声成像、数字射线照相术、发射型计算机成像和核素成像等也逐步发展。
计算机和医学图像处理技术作为这些成像技术的发展基础，带动着现代医学诊断正产生着深刻的变革。
各种新的医学成像方法的临床应用，使医学诊断和治疗技术取得了很大的进展，同时将各种成像技术得到的信息进行互补，也为临床诊断及生物医学研究提供了有力的科学依据。
在口前的影像医疗诊断中，主要是通过观察一组二维切片图象去发现病变体，往往需要借助医生的经验来判定。
至于准确地确定病变体的空间位置、大小、儿何形状及与周圉生物组织的空间关系，仅通过观察二维切片图象是很难实现的。
因此,利用计算机图象处理技术对二维切片图象进行分析和处理，实现对人体器官、软组织和病变体的分割提取、三维重建和三维显示，可以辅助医生对病变体及其它感兴趣的区域进行定性其至定量的分析，可以大大提高医疗诊断的准确性和可靠性。
此外，它在医疗教学、手术规划、手术仿真及各种医学研究中也能起重要的辅助作用。
图1三维可视化大体步骤
医学图像的三维可视化的方法很多，但基本步骤大体相同，如图1。
从CTM/R或超声等成像系统获得二维断层图像，然后需要将图像格式转化成讣算机方便处理的格式。
通过二维滤波，减少图像的噪声影响，提高信噪比和消除图像的尾迹。
采取图像插值方法，对医学关键部位进行各向同性处理,获得体数据。
经过三维滤波后，不同组织器官需要进行分割和归类，对同一部位的不同图像进行配准和融合，以利于进一步对某感兴趣部位的操作。
根据不同的三维可视化要求和系统平台的能力，选择不同的方法进行三维体绘制，实现三维重构。
最重要的环节便是医学图像分割。
理想的图像分割技术能将医学图像中各个部位和组织完整的分离开来，使得医学影像的成像结果更加直观，也使病变部位的识别判断更加准确。
现有的图像分割技术，无论是单一的分割方法，还是基于单一的特征信息进行分割的方法，在对分割精确度和效率要求较高的临床医学领域都难以取得令人满意的分割效果，因而人们开始寻找新的概念和方法并将其引入图像分割领域。
近儿年来提出了儿种新方法，其中基于多特征融合和多谱图像分析相结合的方法应用最为广泛，然而，这些方法多数是基于渐近假设理论发展而来的，且都默认样本数LI趋于无穷大，但实际问题中样本数LI往往是有限的，所以它们在很多方面都出现了问题，有些推广能力较差，有些学习速度较慢，还有些难以收敛，在高维度特征和小样本学习方面都难以取得良好的效果，无法进行分割。
为了解决上述问题，支持向量机（SupportVectorMachine,SVM）方法应运而生。
支持向量机理论最初就是基于样本线性可分时的最优分类超平面理论发展起来的，它的主要LI的是寻找一个规则能把Rn空间内的样本点分成两部分，这个规则就是能使两类样本离超平面的距离最大化的最优分类超平面。
最优分类面的U的就是能将两类样本正确分开且使分类间隔最大，这实际上也是对推广能力的控制，是SVM的核心思想之一。
SVM是以统计学习理论为基础发展起来的一种新的模式识别分类方法，依据Vapnic结构风险最小化原则和VC维理论在小样本、高维度、非线性和局部极小点数据空间中获得了良好的推广能力，被看作是对传统分类器的发展。
尽管支持向量机在图像分割方面获得了较好的分类效果，但在图像中提取的分类样本中含有噪声或野值样本时，可能在构造最优分类面时出现偏差，得到的结果并不是真正的最优分类面，导致分类误差较大而失去作用。
近年来针对医学图像的特点人们乂提出了许多关于SVM算法的改进和应用研究。
经典SVM方法通过参数C对LI标函数中的分类误差项进行惩罚，因而参数C的选择对图像分类结果至关重要，针对此问题，Scholkoph等人提出了v-SVM方法。
针对SVM方法对噪声和野值比较敏感的问题，Suykens等人提出了最小二乘支持向量机（LS-SVM）,简化了运算复杂性，却由于每个样本对分类器的贡献失去了稀疏性的优点。
后来Lin等人提岀了模糊支持向量机方法（Fuzzy-SVM）,这种方法针对不同的样本数据采用不同的惩罚系数，使不同的样本在口标函数中的贡献不同，并且使含有噪声或野值的样本的权值较小，达到消除噪声或野值影响的L1的。
医学图像分割中经常会出现不同类样本数量差异较大的情况，使用SVM经常会出现分类结果偏向于样本数较多的类的问题，而医学图像中某些重要的类（如肿瘤部位）和其他部位相比样本数量比较小，但乂对分类精度要求较高，这就需要更多新的解决方法。
另外一种比较受公众关注的图像分割方法进入人类的视野一一基于水平集方法求解主动轮廓模型的曲线演化方法。
这种方法的原理是利用轮廓曲线的儿何特性建立能够驱动轮廓曲线演化的能量泛函，然后通过最小化能量泛函驱动轮廓曲线逐渐逼近图像中的LI标边界得到演化方程，最后利用水平集函数将轮廓曲线的演化方程转化为求解数值化偏微分方程的问题。
CV-SVM新方法针对医学图像的复杂多样性对支持向量机算法造成的运算速度较慢，识别时间较长的影响，在支持向量机进行特征提取后，使用PCA方法对特征向量进行降维运算，提取有价值的主成分进行分析，使运算量大大下降，识别过程得到简化，提高了运算速度。
再通过两种方法的结合有效地减弱了各自的缺点对分割结果的影响，使医学图像的分割结果在相似灰度及噪声影响下，细节方面表现突出，在肿瘤的自动提取过程中也达到了较为满意的结果。
本文提出利用CV模型和支持向量机相结合的CV-SVM新方法，分别对经典脑图像和医院获取的脑部医学肿瘤图像进行分割和肿瘤提取，通过Matlab平台给出仿真实验结果，通过实验数据和结果分析了新方法对图像分割和提取的影响，证明了利用新方法对医学图像处理的可行性和有效性。
2.方案设计
针对MRI图像，支持向量机方法对噪声和野值比较敏感，所以在图像分割时容易出现边缘分类不清的问题；若图像中出现不同类的样本数LI差异较大时，分类结果还会向数LI较多的类偏移，造成分类偏向问题。
而CV模型水平集方法则需要对整个待分割LI标图像中所有点的水平集函数值进行讣算和更新，工作效率比较低；且在受到偏移场影响后的图像中，演化曲线无法演化到相应的边界上造成较大的分割误差其至分割错误。
这些缺点一旦在医学图像分割领域得到体现就会造成严重后果。
山于支持向量机对噪声和野值的敬感度较高，CV模型受偏移场影响较大，在较复杂的医学图像领域都不能独自胜任，本文提出利用CV模型和支持向量机相结合的CV-SVM新方法。
图2基于CV-SVY的医学图像分割流程
图2给出了CV-SVM的流程图。
利用CV模型的水平集对噪声具有较强的鲁棒性先对图像进行粗略的分割，使用图像的边缘轮廓提取作为初始轮廓线，进行迭代后将图像大致分割为不同的部分，然后使用支持向量机方法对提取的分类的内部区域进行细节修正，修正过程使用分类后的LI标区域进行特征提取，使用主成分分析法对提取出的特征进行降维，使用简化的主成分学习后并细致的分类，最后得到较为准确的目标区域。
3.算法描述
（-）CV模型水平集粗略分割
设Q为图像区域，
r（0为演化曲线，水平集函数定义为：
卩（儿y,Z）＞0
，如果（x,y）在
r（r）内部
0（禺几Z）=0,如果（x,y）在
rfe）
外部
设
c（£）=Ca，y,ij表示t时刻的二维闭合曲面，在超平面
0
上，把图像上的点X到初始曲线C（0）的符号距离
d
定义为水平集函数在X处的函数值
（P（乙t=0）=+d
（P（才'£=0）＞0
0（才,f=0）＜0
点X在C（0）内部,
;点X在C（0）外部，
O初始状态t二0时，有
水平集方法的主要U的是为水平集函数
0（兀t）
构造一个能使它在演化界处的值为零的方程。
设X（t）为演化界面上点的演化路径，由于函数在演化界面处的值为零，即
根据符合函数求导法则，得
©+工冬r=0
Z=1
式中
为X的第i个分量。
定义
="（孔））
为垂直于演化界面的速度，因为
£他耳=（%%…％）（“％•••宀）=“（了⑴）131
用
1w1
代表函数
卩（兀t）
的梯度，F代表
讥儿t）
的速度函数。
水平集函数
讥儿t）
的演化方程就可以写成
0（z）+z7vy|=0
多相位水平集分割n个复朵程度不同的拓扑结构相位只需要
log,n
个水平集函数，且能够很好地解决空白与重叠问题。
针对文中的四分类问题我们只
需考虑
m=log2n
（m二4）个单相位水平集函数
则
的所有零水平集集合构成的分割边界线就是多LI标分割的分界线。
我们用c代表演化的曲线，用
C\
C；
分别代表曲线c内部区域和外部区域的均值。
假设某
均值。
假设某图
ft
山两块近似常量区域组成，其像素值分别为
、
。
其中要检测的U标区域由是像素值
此
其边界为c。
那么
此
、X
分别代表曲线内外的区域了。
令
尸（C）=片（C）+E（C）=JI4/0-c,fdxdy+|丨4一幺「弘力
那么U标的边界就是上式中的能量函数达到最小时的边界，即
F©=幷“）+巧（C）«0
（二）特征提取和归一化
在
（1）中使用图像的边缘轮廓提取作为初始轮廓线，进行迭代后将图像大致分割为不同的部分。
然后使用支持向量机方法对提取的分类的内部区域进行细节修正，具体如下：
针对医学MR图像的复杂特点主要采用纹理特征与灰度特征（灰度共生矩阵）组合的特征提取方法来进行处理。
角二阶矩：
a=｛兀;力『
iJ
对比度：
4=丿"（：
，丿）
/J
相关性：
EE（y＞丿"（/，丿）一匕竹
4=
5
1
2
3
4
5
6
7
8
9
反差距:
。
（厶丿）
1+（/-J）2
和均值:
4=W，）
/=2
差均值:
人=Z讥川）
2=2
方差:
-uYp（2,J）
和方差:
/J
2Nr2
mwa
Z=2
差方差:
—”
4=工G-4）"
7-0
上述儿种特征能够反映图像灰度的变化特性。
10和爛：
血一£匕卫）1昭（C））
/=2
11爛：
心=-SE^Aj）iog（PUj））
iJ
12差爛：
坨一£7"1。
心丿））
/=0
特征⑩到？
的各种爛则主要用来测量共生矩阵的系数变化复杂度和纹理变化复杂度。
脑部医学MRI图像的窗口选择为7*7大小。
在纹理特征提取时使用灰度共生矩阵在四个方向的12种纹理特征共48种特征，再加上象素灰度，7*7窗口内的灰度中值及平均值等特征共51种特征作为图像的特征矢量。
为避免特征的偏移，同时在计算特征向量的内积时避免大数汁算引发的计算溢出的困难，对数据进行归一化，采用如下公式
片_仃戶的二特征向量-特征向量最小值
忆川谊M•-特征向量最大值一特征向虽最小値
（三）主成分分析
归一化后的特征向量的范围被限制在［0,1］之间，简化了运算过程，提高了运算速度。
采用主成分分析法对繁杂的特征向量进行提取，以降低分类代价，提高分类性能。
PCA为了获得数口和汁算量最小的主成分需要从信息压缩最大方向寻求图像的低维信息表达，所以对应的K-L变换的信息燔最小，获得的识别空间也是对原图像的最优逼近。
步骤如下：
（1）对各维数据进行无量纲化；
设
那么得到
其中
是各维特征向量的均值,
是各
Xq.gX,i=1,2,…丿,
维特征向量的标准差。
（2）计算Y的协方差矩阵S:
S=—[K-KZ]x[r-i7]rN
（3）解特征值和特征向量
根据特征方程
uzmo
求解S的特征值
和特征向量u。
而后将
兄
从大到小排列：
人A心A…A1A：
并找到与之对应的°
Uj=［勺.,％…，％］,/=1,2,…川
（4）计算主成分
//=儿旳+力灼+…+y丹
这里
J-1,2,•••，卫MN
（5）选取主成分
P的取值原则要让主成分满足最大程度包含原信息量最小值的条件。
这里依
据主成分贡献率。
at=孕£&
选取主成分。
_
选用Gauss径向基核函数（高斯核函数）
K（x,x'）=exp（-卜一”『/2a2）
，惩罚
因子C取1000,高斯核宽度参数CT2取］,窗口大小为7*7。
单独便fflSVM分劇后的图像
使用CV3M分割后的图像
图3SVM与CV-SVM对比
山图3中可以看岀，CV-SVM的分割结果比经典支持向量机在细节和边缘处的优势更加明显，且对噪声具有较强的鲁棒性。
实验证明，CV-SVM可以应用在复杂度高，精确度要求也较高的MRI脑部医学图像的分割领域。
4.实验结果
为了测试基于CV-SVM的医学图像分割方法的效果，采用青岛大学医学院附属医院的带肿瘤组织的脑部图像作为实验对象。
下面给出CV-SVM分类后提取岀的肿瘤化区域和专家手工标注的肿瘤化区域对比图。
PE标岀肿瘤部位
图4专家手工标注肿瘤区域的脑补图像
手工提収肝瘤部位
SVM提HZ肿瘤鄂位
CV・SVM提取肿瘤部位
图5三种方法提取出的肿瘤区域对比图
表1三种方法分割效果评价
方法
评价标准
手工分割
使用经典SVM
使用CV-SVM
像素个数
2333
2042
2356
时间（S）
7.8
6.73
5.96
与手工分割对比
漏检率
0.1466
0.03129
误检率
0.02186
0.02143
从上图5中我们可以直观的看到三种方法提取出的肿瘤组织的形状差异，111于肿瘤疾病的复杂性和技术条件的局限性，医学图像分割的模板或者说准则还没有形成，我们暂时使用医学专家手工提取的肿瘤组织区域作为比对模板。
山表1可以看岀经典SVM和CV-SVM方法与手工分割方法相比，CV-SVM的像素个数跟接近专家给出方案中的像素个数，而且误检率、漏检率也比经典SVM方法要低，CV-SVM方法可以应用在肿瘤组织区域的提取领域，并可以达到较好的效果。
5.结论
本文总结支持向量机方法对推广能力控制的优点，针对它的一些缺点和不足比如在噪声和野值影响下的误差问题和样本数LI差异带来的分类偏向问题，提出CV模型水平集与支持向量机结合的新方法CV-SVM,并对同一对象分割对比进行验证，随后乂通过对同一对象进行肿瘤提取，并与专家意见进行比较，实验表明：
CV-SVM方法能有效地分割图像并能较准确地识别提取肿瘤部位。
6.参考文献：
[1]张翔•支持向量机及其在医学图像分割中的应用[D].华中科技大学,2004
⑵林瑶，田捷.医学图像分割方法综述[J].模式识别与人工智
能.2002,Vol.15
（2）192-204
[3]吴恩惠.医学影像学[Ml.北京：
人民卫生出版社,2001:
33-49
[4]卢健，胡志忠,杨如乃.医学图像融合技术的研究[J].上海生物医学工程.2006（03）
[5]王栋.基于加权多宽度高斯核函数的支持向量机聚类算法研究[D].中国海洋大学，2011
[6]马波.支持向量机多类分类算法的分析与设计[D].扬州大学,2008
[7]BurgesC.J.C.ATutorialonSupportVectorMachinesforPatternRecognitionEj].
KnowledgeDiscoveryandDataMining,1998,2
（2）:
121-167
[8]边肇棋，张学工•模式识别[M]•北京:
清华大学出版社,2000
[9]郭磊，武优西，刘学娜，颜威利，沈雪勤.基于主成份分析和支持向量机的
MRI图像多目标分割[J].中国生物医学工程学报,2007,26（4）
[10]李俊•基于曲线演化的图像分割方法及应用[D]・上海交通大学,2001
升级会员 