[发明专利]眼底图像的动静脉视网膜血管分割方法

说明书

技术领域

本发明涉及计算机辅助诊断技术领域，具体涉及一种眼底图像的动静脉视网膜血管分割方法。

背景技术

随着计算机技术中的人工智能领域的快速发展，计算机辅助诊断技术也逐渐发展。计算机辅助诊断技术是指通过影像学、医学图像处理技术以及其他可能的生理、生化手段，结合计算机的分析计算，辅助影像科医师发现病灶，提高诊断的准确率。

通常医学影像学中计算机辅助诊断分为三步，具体如下：第一步是把病变从正常结构中提取出来；第二步是图像特征的量化；第三步是对数据进行处理并得出结论。

因为计算机可以全面利用影像信息进行精确的定量计算，去除人的主观性，避免因个人知识和经验的差异而引起的“千差万别”的诊断结果；所以它的结果是不含糊的，是确定的，它使诊断变得更为准确、更为科学。

随着现代高科技的发展，计算机辅助诊断将与图像处理和PACS系统等技术融合，变得更易于操作、也更趋于准确，其临床应用范围将进一步扩大。

在医学检测中，眼睛是唯一可无损检测同时信息丰富的器官。研究指出视网膜血管病变中的血管局限缩窄、弥漫缩窄、动静脉交叉压迫、血管行走改变、铜丝动脉、出血、棉絮斑、硬性渗出以及视网膜神经纤维层缺损与脑卒有显著的相关性。且对于脑卒中的预测，眼底检查仅需40元，而MRI检查则需要上千元，颈动脉超声也需要140元。相比之下眼底检查的性价比最高。眼底图像计算机分析的全自动化的方法，包括可以提供即时的视网膜病变分类，而不需要专家意见，建立以眼底血管视神经预测三高并发症的系统具有其确实的经济意义。因此，视网膜血管的病变检测在对脑卒的辅助检测具有突出作用。其中构建一个动静脉交叉压迫视网膜血管病变的自动检测系统更是其中的关键部分。

对眼底图像进行血管分割、视盘定位和血管分类(动静脉分裂)是视网膜血管的病变检测的基础，现有的血管分割方法需要人工添加标注信息，自动化程度不高。

发明内容

针对现有技术的不足，本发明提供了一种眼底图像的动静脉视网膜血管分割方法。

一种眼底图像的动静脉视网膜血管分割方法，包括如下步骤：

(1-1)按照预设的二值化阈值对经过预处理的眼底图像进行二值化处理，并提取二值化处理后的眼底图像中的中心线和边缘，得到血管树；

(1-2)对所述的血管树分叉处做断开处理得到血管段，并对各个血管段进行线分割得到血管，得到原始血管集；

(1-3)确定误分割血管，并从原始血管集中去除得到全局血管集。

本发明中对眼底图像进行小波变换以对眼底图像进行预处理。

二值化阈值为二值化处理后为血管的像素点个数占整个眼底图像的像素点比例，通常取值为4～20％。二值化阈值越大，则越宽松。作为优选，所述的二值化阈值为14％。

所述步骤(1-3)中的误分割血管包括分为两类，第一类误分割血管基于血管两侧的背景差异确定，第二类基于血管形状确定。

本发明中所述的背景差异值背景颜色的差异，步骤(1-3)中通过如下步骤确定第一类误分割血管：

(a1)针对每个血管，提取该血管两侧背景的特征向量；

任意一侧背景的特征向量根据该侧距离中心线5～10个像素点以内区域中所有像素点(即与中心线的距离小于5～10个像素点)的RGB三个通道上的颜色值得到。

每侧的特征向量为一个三维向量，分别表示血管两侧背景的在RGB三个通道上的颜色值信息。具体实现时获取该侧距离中心线5～10个像素点的区域中的所有像素点在R、G、B三个通道上的颜色值并分别在每个通道上求平均，进而得到该侧的特征向量。

(a2)根据所述特征向量的采用聚类法将所有血管聚为两类，得到的小类即为视盘周围误分割血管。

针对每个血管，计算两侧特征向量的欧氏距离，然后对所有血管对应的欧式距离进行聚类，即完成对血管的聚类。

由于K均值聚类(即K-means算法)不需要调整参数，且运行速度较快。作为优选，所述步骤(a2)中采用K均值聚类法将所有血管聚为两类。

本发明中所述的血管形状实际上指各个血管的连接关系，步骤(1-3)中通过如下步骤确定第二类误分割血管：

确定划分出原始血管集的眼底图像中的环状结构，针对各个环状结构，若该环状结构中长度最大的血管的长度小于预设的分割长度阈值，则认为该环状结构中所有血管均为第二类误分割血管，进一步进行如下操作：