视网膜血管弯曲度确定方法与流程

本发明涉及计算机领域，尤其涉及一种视网膜血管弯曲度确定方法。

背景技术：

随着科学技术的发展，目前已可以通过医学影像中关于结构和组织所反映出的各种特征和细节，进行相关的分析。例如眼底影像中包含了丰富的玻璃体、视网膜、血管网和脉络膜等细节。视网膜血管在正常情况下，其形态跟结构一直处于稳定状态。然而高血压、糖尿病及冠状动脉硬化等心脑血管疾病会引起视网膜血管弯曲度等结构变化。如果能够尽早发现视网膜血管弯曲度变化，可以及时采取措施防止这些疾病的危害。

当前的眼底图像识别方式大多通过眼底照相机获取眼底图像后，靠肉眼识别视网膜血管的弯曲度变化，视网膜血管弯曲度识别效率有待提高。

技术实现要素：

本发明实施例提供一种视网膜血管弯曲度确定方法，能够提高网膜血管弯曲度识别效率。

本发明实施例采用如下技术方案：

一种视网膜血管弯曲度确定方法，包括：

构建全卷积深度神经网络，所述全卷积深度神经网络能够从眼底图像中分割出血管和视盘；

获取目标眼底图像，应用所述全卷积深度神经网络从所述目标眼底图像中分割出所述目标眼底图像中的血管和视盘；

根据所述目标眼底图像中的血管和视盘，确定所述目标眼底图像中血管的弯曲度。

本发明实施例公开一种视网膜血管弯曲度确定方法，构建全卷积深度神经网络，应用全卷积深度神经网络可以从获取的目标眼底图像中分割出目标眼底图像中的血管和视盘。根据目标眼底图像中的血管和视盘，确定目标眼底图像中血管的弯曲度，实现基于全卷积深度神经网络自动确定视网膜血管弯曲度，从而能够提高视网膜血管弯曲度识别效率。

应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本公开。

附图说明

此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本发明的实施例，并与说明书一起用于解释本发明的原理。

图1为本发明实施例示出的视网膜血管弯曲度确定方法的流程图。

图2为本发明实施例示出的眼底图像分割的流程示意图。

图3为本发明实施例示出的分割血管和视盘的示意图。

图4为本发明实施例示出的对血管弯曲度计算的示意图。

图5为本发明实施例示出的全卷积深度神经网络的架构示意图。

图6为本发明实施例示出的特征提取部分和上采样部分的参数配置示意图。

图7a为本发明实施例示出的提取血管的灰度示意图。

图7b为本发明实施例示出的提取血管的二值化示意图。

图7c为本发明实施例示出的提取血管的距离变换示意图。

图7d为本发明实施例示出的提取血管的中心线示意图。

图8a为本发明实施例示出的数字化的曲率验证圆弧示意图。

图8b为本发明实施例示出的圆弧上曲线分布对比的示意图。

图9a为本发明实施例示出的曲率计算原理的示意图。

图9b为本发明实施例示出的血管曲率分布示意图。

具体实施方式

这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本发明的一些方面相一致的装置和方法的例子。

本发明实施例提供一种视网膜血管弯曲度确定方法，通过预训练好的深度学习神经网络能够快速、准确的完成眼底图像和视盘的提取。

如图1所示，本发明实施例提供一种视网膜血管弯曲度确定方法，包括：

s1、构建全卷积深度神经网络，所述全卷积深度神经网络能够从眼底图像中分割出血管和视盘。

s2、获取目标眼底图像，应用所述全卷积深度神经网络从所述目标眼底图像中分割出所述目标眼底图像中的血管和视盘。

s3、根据所述目标眼底图像中的血管和视盘，确定所述目标眼底图像中血管的弯曲度。

本发明实施例提供的视网膜血管弯曲度确定方法，构建全卷积深度神经网络，全卷积深度神经网络能够从眼底图像中分割出血管和视盘，获取目标眼底图像，应用全卷积深度神经网络从目标眼底图像中分割出目标眼底图像中的血管和视盘，根据目标眼底图像中的血管和视盘，确定目标眼底图像中血管的弯曲度，从而实现基于全卷积深度神经网络自动确定视网膜血管弯曲度，从而能够提高网膜血管弯曲度识别效率。

在一个实施例中，s1构建的全卷积深度神经网络可以为经过处理及训练的后的全卷积深度神经网络。如图2所示，应用训练图像(眼底图像)得到全卷积网络后进行系列处理后得到全卷积深度神经网络(s1中的全卷积深度神经网络)，眼底图像在输入后通过多任务学习分别进行血管的分割和视盘的分割。具体的，图2中对事先准备好的训练图像进行预处理，预处理过后进行数据扩增，将扩增后的训练作为训练数据，输入到全卷积网络中，根据损失函数计算出二值损失，使用梯度下降算法更新网络权值，根据其在验证集上的最优结果保存模型(s1中的全卷积深度神经网络)。对测试图像进行与训练时相同的预处理，然后将数据输入最优的保存模型，自动分割出视盘和血管。图3所示为眼底图像通过训练出的全卷积深度神经网络分割出的血管和视盘。如图4所示，采用通过多任务学习进行血管和视盘的分割通过数字图像处理完成血管弯曲度计算。具体的，对于一张视网膜眼底图像，通过训练好的全卷积深度神经网络对输入的视网膜眼底照片进行分割，对分割出的视盘根据以视盘为起点中心，向外延伸4个视盘直径距离圆形区域，确定血管曲率计算范围，同时从视网膜眼底照片分割出的血管提取出血管中心线，通过血管树去除分支交叉点后，获得最终需要计算曲折度的roi二值图像，根据roi二值图像获取主血管的分支曲线，提取一张血管树中n段平滑后的分支血管曲线，计算每条曲线上每个像素点上的曲率值。

基于fcn的神经网络模型。该模型的结构分为特征提取部分和上采样两部分。其中特征提取部分，每经过一个池化层就一个尺度，包括原图一共5个尺度。上采样部分的操作为每上采样一次，就和特征提取部分对应的通道数相同尺度融合。基于多任务学习进行眼底视网膜照片的血管和视盘分割。血管分割和视盘分割共享所述全卷积深度神经网络参数。最后输出分割后的血管和视盘照片。训练所述神经网络模型时，对所述神经网络中的每一个卷积在relu前使用batchnormalization函数计算，让网络中间数据的分布尽量规范化在传至下一层。训练所述神经网络模型时，对所述神经网络中的每一个卷积使用leakyrelu方程运算后再传至下一层。训练所述神经网络模型时，使用binarycrossentropyloss作为损失函数。使用带nesterovmomentum的随机梯度下降法作为所述神经网络的损失函数的优化算法。训练所述神经网络模型时，使用l2weightdecay对所述神经网络中的参数进行正则化。血管中心线的提取。通过所述神经网络获取的视网膜血管灰度图二值化处理得到血管中心线图像。视盘中心的确定和视盘直径的测量。通过所述神经获取得到视网膜灰度图像确定视盘中心和视盘直径的测量。

具体的：血管中心线的提取时，通过所述神经网络获取的视网膜血管灰度图进行形态学结合距离变换处理得到血管中心线图像。视盘中心的确定和视盘直径的测量时，通过所述神经获取得到视网膜灰度图像确定视盘中心和视盘直径的测量。曲率测量范围的确定时，以视盘为起点中心，向外延伸4个视盘直径距离圆形区域，确定血管曲率计算范围。血管分支交叉点的去除。通过血管树去除分支交叉点后，获得最终需要计算曲折度的roi二值图像。主血管分支的提取时，通过对所述roi图像进行血管分支提取，获取主血管的分支曲线。曲率计算。提取一张血管树中n段平滑后的分支血管曲线，计算每条曲线上每个像素点上的曲率值。

在一个实施例中，所述构建全卷积深度神经网络包括：

构建所述全卷积深度神经网络的特征提取部分和上采样部分，全卷积深度神经网络的架构如图5所示，其中，特征提取部分和上采样部分的参数配置如图6所示。所述特征提取部分，每经过一个池化层就一个尺度，包括原图一共有5个尺度；所述上采样部分，每上采样一次，就和特征提取部分对应的通道数相同尺度融合，同时该架构在每层的卷积之后增加batchnormalization，保证神经网络中间数据尽量规范化。

在一个实施例中，采用以下方式训练所述全卷积深度神经网络，得到针对医疗眼底图像的血管和视盘分割的全卷积深度神经网络：

使用binarycrossentropyloss作为损失函数；

使用带nesterovmomentum的随机梯度下降法作为优化学习算法；

对所述全卷积深度神经网络中的各参数使用l2weightdecay正则化，防止参数过大导致的过拟合。

通过以上训练过程，得到针对医疗眼底图像的血管和视盘分割的全卷积深度神经网络。

在一个实施例中，采用以下方式训练所述全卷积深度神经网络：

获取训练数据集和测试数据集，训练数据集和测试数据集中分别包含20张输入训练图像的文件夹和20张标记好的眼底血管标签数据、20张标记好的眼底视盘标签数据文件夹；

对数据集中的眼底图像进行预处理，针对全卷积深度神经网络所进行的训练，遍历训练数据集中的20张输入图像文件夹的所有图像，用训练样本训练全卷积深度神经网络。

在一个实施例中，参考图4所示，所述根据所述目标眼底图像中的血管和视盘，确定所述目标眼底图像中血管的弯曲度包括：

根据视盘中心、视盘直径确定曲率测量范围；

提取血管中心线；

根据所述血管中心线及所述曲率测量范围，提取血管分支；

对所述血管分支进行曲率计算。

在一个实施例中，所述提取血管中心线包括：

为提取血管中心线，首先对通过深度学习获得的视网膜血管灰度图像进行二值化处理，如图7a所示。

使用自适应大津阈值算法对图像进行全局阈值处理，并采用形态学处理去除孤立或细小的噪声点，如图7b所示。

使用形态学处理得到骨骼化图像，如图7c所示，并通过欧氏距离变换计算血管上每个像素到离它最近的背景像素的距离获得血管距离图像，得到保持距离参数的血管中心线图像，如图7d所示。

在一个实施例中，所述根据视盘中心、视盘直径确定曲率测量范围包括：

对深度学习获取得到视网膜视盘灰度图像进行二值化处理，利用主活动轮廓cv(chan-vese)模型将视盘灰度图像分割成前景和背景形成模板，对模板中近似圆形的视盘区域计算面积估计视盘中心和等效圆直径；

视盘中心为起点中心，向外延伸4个视盘直径距离圆形区域，即为血管曲率的计算范围。

在一个实施例中，所述根据所述血管中心线及所述曲率测量范围，提取血管分支包括：

去除血管分支交叉点；

血管树去除分支交叉点后，获得最终需要计算曲折度的roi图像，将所有分支打断并分别统计每条分支曲线上所有像素的曲率值；

使用canny算子提取血管边缘线并进行角点检测，保留分支曲线上像素个数总数大于12的分支，对每段分支血管曲线进行样条函数拟合，得到平滑后的血管分支曲线。图8a所示为一段数字图像形式的圆弧，图中点表示该弧线上的像素点位置。图8b所示实线表示圆弧上的真实曲线分布，虚线表示拟合平滑后的曲线分布。

在一个实施例中，所述血管分支交叉点的去除包括：

将在视网膜象限区域并二值化处理后的血管中心线图作为感兴趣区域roi图像；

对roi图像进行形态学运算找到血管树中所有的分支交叉点，并构建一个指定半径2的平面圆盘形的结构元素；

对所有的分支交叉点进行膨胀操作并反相，获得一个交叉点模板图，roi图像与交叉点模板图相与得到去除分支点的血管树图。

由于血管树纵交复杂且分支多，在曲率计算之前，本实施例对血管之间的交叉点进行处理。

在一个实施例中，所述对所述血管分支进行曲率计算包括：

如图9a所示，提取一张血管树中n段平滑后的分支血管曲线，一条血管被定义成曲线，对所述每条曲线进行曲线拟合后，近似圆上的某段圆弧，采用k对，k>7(一般取8～12拟合效果较好)像素点的大小作为链码；确定好计算区域后，采用最小二乘法对血管曲线段进行圆拟合，获得对应像素点上的等效圆半径r，曲率值即为1/r。如图9b所示为计算后的血管曲率分布。

本发明实施例通过预训练好的深度学习神经网络能够快速、准确的完成眼底图像和视盘的提取。

本发明实施例提供的视网膜血管弯曲度确定方法，构建全卷积深度神经网络，全卷积深度神经网络能够从眼底图像中分割出血管和视盘，获取目标眼底图像，应用全卷积深度神经网络从目标眼底图像中分割出目标眼底图像中的血管和视盘，根据目标眼底图像中的血管和视盘，确定目标眼底图像中血管的弯曲度，从而实现基于全卷积深度神经网络自动确定视网膜血管弯曲度，从而能够提高网膜血管弯曲度识别效率。

以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。

本领域技术人员在考虑说明书及实践这里公开的公开后，将容易想到本公开的其它实施方案。本申请旨在涵盖本公开的任何变型、用途或者适应性变化，这些变型、用途或者适应性变化遵循本公开的一般性原理并包括本公开未公开的本技术领域中的公知常识或惯用技术手段。