一种基于SVM的医疗图像血管识别方法技术

本发明专利技术公开一种基于SVM的医疗图像血管识别方法，包括：首先SVM分类器是通过专家手工分割后的样本训练出来的SVM模型对FCM自动选取的样本训练出的SVM模型进行交叉验证、优化后的结果，使分割结果更加精确；然后采用SVM分割血管，其实就是将像素点分为前景(即血管)和背景(即除血管外其他部分)两类，然后把血管部分提取出来，最后通过形态学处理以及阈值化可以达到增强血管网络、保留血管分叉以及交汇部分，最终转化为二值图像可以更直接地反映血管分布。本发明专利技术采用了FCM、SVM和形态学图像处理相结合，可使识别效果更好。

【技术实现步骤摘要】

本专利技术属于计算机领域，提出一种基于SVM的医疗图像血管识别方法，能很好的提取出医疗图像血管网络，尤其是眼底图像血管。
技术介绍
在医疗诊断和医学研究中，血管是非常重要的生物组织，许多器官发生病变时，其血管会出现异常。如眼睛中的视网膜血管、心脏中的冠状动脉血管、肺部小血管等。其中，视网膜是人体中唯一可以以非创伤方式直接观察到的较为深层的血管，而且视网膜血管是眼底图像中可见的最主要解剖结构，其结构特征变化直接反映身体的各种病变。视网膜眼底图像还可以作为判断白内障等多种疾病的重要依据，计算机的自动眼底图像处理与分析在辅助医疗诊断中有重要的价值，而且所提取的血管可以作为一种特征用于白内障眼底图像分类。此外，考虑到眼底筛查结果客观性、可重复性、准确性以及大批量的需求，临床上对视网膜血管的自动检测以及定量分析提出来迫切要求。因此，采取有效的方法自动提取视网膜眼底图像中的血管结构，有助于对白内障等眼底疾病的早期诊断、治疗跟踪，在临床上具有重要的应用价值。目前，针对视网膜眼底图像血管的提取，国内外已经做了较多的研究，不断有新方法涌现。根据分割时是否有采用标准图像提供的特征数据规则，现有的视网膜血管分割方法分为非监督和监督方法两大类。非监督方法主要包括局部滤波法、局部自适应阈值法和血管跟踪法等。有监督的方法包括基于脊波的血管分割法和基于神经网络的像素分类法。作为图像分析和模式识别的重要问题，图像分割的方法也越来越被众多人关注，诸如阈值分割法、边缘检测分割法、基于聚类分割法等。支持向量机(SVM)在图像分割领域也占有一席之地。2012年，由NilanjanDey，AnamitraBardhanRoy等发表的InternationalJournalofComputerScienceandNetwork论文“FCMBasedBloodVesselSegmentationMethodforRetinalImages”中，提出了一种基于基础图像处理技术和FCM技术结合的分割眼底图像血管的方法。首先提取眼底图像的G分量通道(即绿色通道)，接着进行直方图均衡化对图像进行增强。直方图均衡化是一种利用灰度变换自动调节图像对比度质量的方法，基本思想是通过灰度级的概率密度函数求出灰度变换函数，它是一种以累计分布函数变换方法为基础的直方图修正法。直方图均衡化后紧接着对眼底图像的前景与背景相减，然后进行中值滤波，完成图像的初步处理。其中，中值滤波可以去除图像中的椒盐噪声，平滑效果优于均值滤波，在抑制噪声的同时还能够保持图像的边缘清晰。中值滤波后通过FCM对眼底图像进行分割，即通过迭代求值更新隶属度矩阵和聚类中心，最终将图像分为前景背景两类(即血管和非血管)，提取出血管网络。最后经过二值化处理提取出最终的血管图像。2002年，由ConorHeneghan,JohnFlynn等发表的MedicalImageAnalysis论文“Characterizationofchangesinbloodvesselwidthandtortuosityinretinopathyofprematurityusingimageanalysis”提出一种利用形态学技术提取眼底血管的方法。形态学是一种应用于图像处理和模式识别领域的方法，具有完备的数学基础，这为形态学图像分析和处理等奠定了坚实的基础。经过形态学的处理以及滤波可以达到增强血管网络、保留血管分叉以及交汇部分。此论文首先用形态学预处理来突出血管结构，然后用几个算子进一步强调细血管结构、过滤血管，最后使用形态学滤波以及阈值化，得到最终的血管图像。然而只通过上面介绍的已有方法不能太精确的提取到血管，尤其是患有白内障等模糊图像的眼底血管，使得一般的提取方法所达到的效果并不理想，采用SVM，并通过专家手动分割后训练的SVM进行优化，可以更有效精确地提取血管网络。
技术实现思路
本专利技术提出一种基于SVM的医疗图像血管识别方法。结合了FCM自动选取训练样本和专家手工标记相结合，使血管的提取效果更好，尤其是眼底图像的血管提取效果，患有白内障的模糊血管也可以提取出部分血管。由于视网膜是人体中唯一可以以非创伤方式直接观察到的较为深层的血管，且眼底图像在医疗图像里运用最广泛，所以本专利技术均以眼底图像为例介绍。为实现上述目的，本专利技术采用如下的技术方案：一种基于SVM的医疗图像血管识别方法包括以下步骤：S1、获取眼底图像；S2、对所述眼底图像进行血管的绿色通道分量提取，得到待处理眼底图像；S3、通过FCM自动选取训练样本，提取样本颜色和纹理特征，即分别提取血管部分和非血管部分的特征，然后训练SVM模型；S4、通过专家标记后的图像训练出的SVM模型来交叉验证FCM自动选取的样本所训练的SVM模型，优化SVM参数；S5、采用优化好的SVM分类器对待处理眼底图像进行分割，即初步提取眼底图像血管；S6、对所述初步提取眼底图像血管进行形态学操作及阈值化处理，实现最终提取眼底血管。作为优选，基于FCM训练样本自动选取以及训练SVM模型包括以下步骤：(1)、灰度化以及初始化参数：像素点的灰度值作为像素点的属性值，初始化参数，设定聚类个数，取c＝2(即分为血管和背景)，以及收敛误差ε；(2)、更新迭代求聚类中心、隶属度矩阵：对图像所有领域窗口计算局部空间信息、局部灰度信息、局部相似性度量，然后运用上述参数值计算线性加权求和图像，更新隶属度矩阵，更新聚类中心，用新的聚类中心计算方法更新迭代，直到灰度之差小于ε为止；(3)、对像素点分类：通过隶属度函数uk(xi,yi)(k＝1,2,…,c；i＝1,2,…,n)对像素点进行分类,uJ(xi,yi)＝max(u1(xi,yi),u2(xi,yi),…,uc(xi,yi)),即第i个像素点(xi,yi)属于第J类，所述隶属度函数uk(xi,yi)为SVM参数；(4)、训练样本的自动选取和标注：对于属于第J类的像素点，随机选取个数为uJ/t的像素点作为第J类的训练样本，其中，uJ是第J类像素点的个数，t值的选取依赖于样本的个数。作为优选，步骤S6中的形态学操作及阈值化处理包括：灰度取反、高帽滤波、阈值化处理。本专利技术的基于SVM的医疗图像血管识别方法，包括：首先SVM分类器是通过专家手工分割后的样本训练出来的SVM模型对FCM自动选取的样本训练出的SVM模型进行交叉验证、优化后的结果，使分割结果更加精确；然后采用SVM分割血管，其实就是将像素点分为前景(即血管)和背景(即除血管外其他部分)两类，然后把血管部分提取出来，最后通过形态学处理以及阈值化可以达到增强血管网络、保留血管分叉以及交汇部分，最终转化为二值图像可以更直接地反映血管分布。本专利技术采用了FCM、SVM和形态学图像处理相结合，可使识别效果更好。与现有技术相比，本专利技术具有以下优点：本专利技术以大量眼底图像数据为条件，结合了模式识别和图像处理技术，运用FCM、SVM和形态学图像处理作为基础算法。同时，训练集通过由FCM处理结果以及人工处理结果两部分组成，通过他们的交叉验证与优化后得到分类器，这样既考虑了分类的准确性，也考虑了数据集的多样性，能更好的提高血管识别的准确率。附图说明结合附图，从下面对本专利技术实施例的详细描述，将更好地理解本专利技术，其中：图1为本文档来自技高网...

【技术保护点】
一种基于SVM的医疗图像血管识别方法，其特征在于，包括以下步骤：S1、获取眼底图像；S2、对所述眼底图像进行血管的绿色通道分量提取，得到待处理眼底图像；S3、通过FCM自动选取训练样本，提取样本颜色和纹理特征，即分别提取血管部分和非血管部分的特征，然后训练SVM模型；S4、通过专家标记后的图像训练出的SVM模型来交叉验证FCM自动选取的样本所训练的SVM模型，优化SVM参数；S5、采用优化好的SVM分类器对待处理眼底图像进行分割，即初步提取眼底图像血管；S6、对所述初步提取眼底图像血管进行形态学操作及阈值化处理，实现最终提取眼底血管。

【技术特征摘要】
1.一种基于SVM的医疗图像血管识别方法，其特征在于，包括以下步骤：S1、获取眼底图像；S2、对所述眼底图像进行血管的绿色通道分量提取，得到待处理眼底图像；S3、通过FCM自动选取训练样本，提取样本颜色和纹理特征，即分别提取血管部分和非血管部分的特征，然后训练SVM模型；S4、通过专家标记后的图像训练出的SVM模型来交叉验证FCM自动选取的样本所训练的SVM模型，优化SVM参数；S5、采用优化好的SVM分类器对待处理眼底图像进行分割，即初步提取眼底图像血管；S6、对所述初步提取眼底图像血管进行形态学操作及阈值化处理，实现最终提取眼底血管。2.如权利要求1所述的基于SVM的医疗图像血管识别方法，其特征在于，基于FCM训练样本自动选取以及训练SVM模型包括以下步骤：(1)、灰度化以及初始化参数：像素点的灰度值作为像素点的属性值，初始化参数，设定聚类个数，取c＝2(即分为血管和背景)，以及收敛误差ε；(2)、更新迭代...

【专利技术属性】
技术研发人员：胡启东，李建强，张苓琳，韩赫，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：北京;11