眼底图像的动静脉视网膜血管分割方法
【技术领域】
[0001]本发明涉及计算机辅助诊断技术领域，具体涉及一种眼底图像的动静脉视网膜血管分割方法。
【背景技术】
[0002]随着计算机技术中的人工智能领域的快速发展，计算机辅助诊断技术也逐渐发展。计算机辅助诊断技术是指通过影像学、医学图像处理技术以及其他可能的生理、生化手段，结合计算机的分析计算，辅助影像科医师发现病灶，提高诊断的准确率。
[0003]通常医学影像学中计算机辅助诊断分为三步，具体如下:第一步是把病变从正常结构中提取出来；第二步是图像特征的量化；第三步是对数据进行处理并得出结论。
[0004]因为计算机可以全面利用影像信息进行精确的定量计算，去除人的主观性，避免因个人知识和经验的差异而引起的“千差万别”的诊断结果；所以它的结果是不含糊的，是确定的，它使诊断变得更为准确、更为科学。
[0005]随着现代高科技的发展，计算机辅助诊断将与图像处理和PACS系统等技术融合，变得更易于操作、也更趋于准确，其临床应用范围将进一步扩大。
[0006]在医学检测中，眼睛是唯一可无损检测同时信息丰富的器官。研宄指出视网膜血管病变中的血管局限缩窄、弥漫缩窄、动静脉交叉压迫、血管行走改变、铜丝动脉、出血、棉絮斑、硬性渗出以及视网膜神经纤维层缺损与脑卒有显著的相关性。且对于脑卒中的预测，眼底检查仅需40元，而MRI检查则需要上千元，颈动脉超声也需要140元。相比之下眼底检查的性价比最高。眼底图像计算机分析的全自动化的方法，包括可以提供即时的视网膜病变分类，而不需要专家意见，建立以眼底血管视神经预测三高并发症的系统具有其确实的经济意义。因此，视网膜血管的病变检测在对脑卒的辅助检测具有突出作用。其中构建一个动静脉交叉压迫视网膜血管病变的自动检测系统更是其中的关键部分。
[0007]对眼底图像进行血管分割、视盘定位和血管分类(动静脉分裂)是视网膜血管的病变检测的基础，现有的血管分割方法需要人工添加标注信息，自动化程度不高。
【发明内容】
[0008]针对现有技术的不足，本发明提供了一种眼底图像的动静脉视网膜血管分割方法。
[0009]一种眼底图像的动静脉视网膜血管分割方法，包括如下步骤:
[0010](1-1)按照预设的二值化阈值对经过预处理的眼底图像进行二值化处理，并提取二值化处理后的眼底图像中的中心线和边缘，得到血管树；
[0011](1-2)对所述的血管树分叉处做断开处理得到血管段，并对各个血管段进行线分割得到血管，得到原始血管集；
[0012](1-3)确定误分割血管，并从原始血管集中去除得到全局血管集。
[0013]本发明中对眼底图像进行小波变换以对眼底图像进行预处理。
[0014]二值化阈值为二值化处理后为血管的像素点个数占整个眼底图像的像素点比例，通常取值为4?20%。二值化阈值越大，则越宽松。作为优选，所述的二值化阈值为14%。
[0015]所述步骤(1-3)中的误分割血管包括分为两类，第一类误分割血管基于血管两侧的背景差异确定，第二类基于血管形状确定。
[0016]本发明中所述的背景差异值背景颜色的差异，步骤(1-3)中通过如下步骤确定第一类误分割血管:
[0017](al)针对每个血管，提取该血管两侧背景的特征向量；
[0018]任意一侧背景的特征向量根据该侧距离中心线5?10个像素点以内区域中所有像素点(即与中心线的距离小于5?10个像素点)的RGB三个通道上的颜色值得到。
[0019]每侧的特征向量为一个三维向量，分别表示血管两侧背景的在RGB三个通道上的颜色值信息。具体实现时获取该侧距离中心线5?10个像素点的区域中的所有像素点在R、G、B三个通道上的颜色值并分别在每个通道上求平均，进而得到该侧的特征向量。
[0020](a2)根据所述特征向量的采用聚类法将所有血管聚为两类，得到的小类即为视盘周围误分割血管。
[0021]针对每个血管，计算两侧特征向量的欧氏距离，然后对所有血管对应的欧式距离进行聚类，即完成对血管的聚类。
[0022]由于K均值聚类(即K-means算法)不需要调整参数，且运行速度较快。作为优选，所述步骤(a2)中采用K均值聚类法将所有血管聚为两类。
[0023]本发明中所述的血管形状实际上指各个血管的连接关系，步骤(1-3)中通过如下步骤确定第二类误分割血管:
[0024]确定划分出原始血管集的眼底图像中的环状结构，针对各个环状结构，若该环状结构中长度最大的血管的长度小于预设的分割长度阈值，则认为该环状结构中所有血管均为第二类误分割血管，进一步进行如下操作:
[0025]确定该环状结构的中心，并计算该中心到长度大于或等于分割长度阈值的血管的最短距离(即该中心到距离其最近的长度大于或等于分割长度阈值的血管的距离)，认为以该中心为圆心、最短距离为半径的圆形区域内所有血管为第二类误分割血管。
[0026]分割长度阈值的设定根据实际经验值设定，通常为分割长度阈值为α = χ/60?x/45，X为眼底图像的横向大小，即眼底图像沿横轴的长度。
[0027]未作特殊说明，本发明中对长度、距离、图片大小等参数进行衡量时统一以像素点为单位。
[0028]与现有技术相比，本发明在得到原始血管集后进一步还利用血管的背景和形状确定误分割的血管，可以有效去除因为拍照造成的环状反光、视盘周围的非血管的跃阶边缘、斑状病变以及出血病变等原因造成的误分割血管。
【附图说明】
[0029]图1为本实施例的眼底图像；
[0030]图2为本实施例中对眼底图像进行血管分割的流程图；
[0031]图3为血管分割得到的原始血管集的示意图；
[0032]图4为血管分割得到的全局血管集的示意图。
【具体实施方式】
[0033]下面将结合附图和具体实施例对本发明进行详细描述。
[0034]本实施例以图1所示的眼底图像为例来说明眼底图像的动静脉视网膜血管分割方法，该眼底图像的大小为3000 X 3000ο由拍照造成的环状反光、视盘周围的非血管的跃阶边缘、斑状病变以及出血病变等原因，该眼底图像中存在亮环。
[0035]对该眼底图像进行动静脉视网膜血管分割的分类流程如图2所示，包括如下步骤:
[0036](1-1)对眼底图像进行小波变换(IUWT小波)，按照预设的二值化阈值对经过小波变换的眼底图像进行二值化处理，并提取二值化处理后的眼底图像中的中心线和边缘，得到血管树；
[0037](1-2)对血管树分叉处做断开处理得到血管段，并对各个血管段进行线分割得到血管，得到原始血管集。
[0038]对血管树分叉处做断开处理时:当血管树中的血管中心线中多根中心线汇集到一点时，去除中心点(汇集的交叉点)，得到单独的多根血管中心线。
[0039]对各个血管段进行线分割时:以每一根中心线作为一个血管段。血管段为一条曲线，运用图像处理的线分割的传统方法，将曲线用多根直线逼近。得到的多根直线，每根直线即代表一根血管，所有直线的集合即为原始血