一种OCT图像中脉络膜新生血管分割算法制造技术

本发明专利技术公开了一种OCT图像中脉络膜新生血管分割算法，包括以下步骤：S01：对训练图像设计一种结构先验学习方法，构建结构先验矩阵，所述结构先验矩阵用于区分脉络膜新生血管区域和背景区域；S02：基于所述结构先验矩阵将OCT原图像转换为显著性增强图像，用于增强脉络膜新生血管区域的显著性；S03：在所述显著性增强图像上使用多尺度分析，将所述显著性增强图像划分为m个尺度；S04：基于每种尺度训练得到m个训练好的卷积神经网络模型；S05：利用步骤S01、S02、S03对测试图像进行处理，利用步骤S04中训练好的所述卷积神经网络模型进行测试，输出m个分割结果，对m个分割结果进行融合即为最终的分割结果，本发明专利技术可以显著提高OCT图像中脉络膜新生血管的分割精度。

An algorithm for choroidal neovascularization segmentation in OCT images

The invention discloses a choroidal neovascularization in OCT image segmentation algorithm, which comprises the following steps: S01: a structure prior to learning method training image design, construct the structure of the prior matrix, matrix structure prior to distinguish between CNV region and background region for S02; the structure based on the prior matrix OCT of the original image is converted to significantly enhance the image for enhanced CNV region; S03: a significant increase in the use of multi-scale image analysis, the image enhancement was divided into m scale; S04: each scale training m trained convolutional neural network based on the model; S05: dealing with the test image using S01, S02, step S03, tested by the convolution of the trained neural network model in step S04, output m segmentation Results the M segmentation results were fused as the final segmentation results, and the proposed method can significantly improve the segmentation accuracy of choroidal neovascularization in OCT images.

【技术实现步骤摘要】
一种OCT图像中脉络膜新生血管分割算法
本专利技术涉及一种OCT图像中脉络膜新生血管分割算法，属于视网膜图像分割

技术介绍
现有的脉络膜新生血管自动分割技术大部分基于眼底荧光素血管造影图像。相比较眼底荧光素血管造影，OCT图像具有无创、高速、高分辨率、三维成像等优点，其对于老年性变性黄斑等临床常见眼科疾病的辅助诊断具有更重要的临床意义。目前尚未有基于OCT图像的脉络膜新生血管分割算法。OCT图像中的脉络膜新生血管分割面临诸多挑战：纹理变化较大、存在灰度不同质性、形状和大小不一致、边界模糊、存在大量斑噪声等。这些问题使得传统的方法很难取得较为精准的分割效果。卷积神经网络具有强大的学习能力，其在医学图像分割(例如，MR脑图像的灰白质分割、电子显微镜下细胞膜的分割、乳腺病理图像的有丝分裂检测等)中已取得了巨大的成功。所以考虑将该框架用于OCT图像中的脉络膜新生血管分割任务中。然而，由于OCT图像特点较为复杂，直接使用卷积神经网络对脉络膜新生血管进行分割存在两个不足：(1)传统的方法首先需要将图像分成若干小块(patch)，然后基于patch训练卷积神经网络分割模型。然而，patch与patch是存在一些结构相关性(例如局部相似性)，传统的方法在训练卷积神经网络分割模型时并未考虑到这种有效的结构信息，从而限制了模型的分割性能。(2)传统的模型是建立在单一尺度patch上的。脉络膜新生血管的大小是不固定的，所以单一尺度的patch很难获取有效的上下文信息，从而影响了分割精度。
技术实现思路
本专利技术所要解决的技术问题是，提供一种能够提高OCT图像中脉络膜新生血管的分割精度的基于多尺度结构先验卷积神经网络的OCT图像中脉络膜新生血管分割算法。为解决上述技术问题，本专利技术采用的技术方案为：一种OCT图像中脉络膜新生血管分割算法包括以下步骤：S01：对训练图像设计一种结构先验学习方法，构建结构先验矩阵，所述结构先验矩阵用于区分脉络膜新生血管区域和背景区域，包括以下步骤：a：对训练图像进行超像素分割得到若干超像素；b：提取特征；c：根据所述训练图像的正确标注将每个所述超像素标记为2类，分别为脉络膜新生血管区域和背景区域；d：使用标记好的若干所述超像素构造字典；e：对所有所述超像素进行分类，得到全局结构先验图；f：基于所述全局结构先验图，计算局部相似结构先验，求得所述结构先验矩阵；S02：基于所述结构先验矩阵将OCT原图像转换为显著性增强图像，用于增强脉络膜新生血管区域的显著性；S03：在所述显著性增强图像上使用多尺度分析，将所述显著性增强图像划分为m个尺度；S04：基于每种尺度训练得到m个训练好的卷积神经网络模型；S05：利用步骤S01中步骤a、b、e、f对测试图像进行超像素分割、提取特征、分类以及计算结构先验矩阵，利用步骤S02基于S05中得到的所述结构先验矩阵对测试图像进行图像转换，利用步骤S03将转换后的所述测试图像划分为m个尺度，利用步骤S04中训练好的所述卷积神经网络模型进行测试，输出m个分割结果，对m个分割进行融合即为最终的分割结果。使用SLIC算法进行超像素分割。所述特征包括每个所述超像素的平均灰度值、基于共生矩阵的纹理特征以及局部灰度特征。使用K-means算法构造词典，假设在训练集中有N个病人的数据，每个数据包括2分类，使用K-means聚成K类，则N个病人的数据共可聚成2KN类，获得2KN个聚类中心，所述聚类中心组成字典D，如公式(1)所示：D＝[C1,1,C1,2…C1,K,B1,1,B1,2,..B1,K,…Cn,k..Bn,k…CN,1,..CN,K,BN,1,..BN,K](1)式中，Cn,k表示第n个病人的来自脉络膜新生血管区域的第k个聚类中心；Bn,k表示第n个病人的来自背景区域的第k个聚类中心，n＝1,2,…N,k＝1,2,…K。使用稀疏表示对每个所述超像素进行分类，分类过程形式化如公式(2)所示，式中，x为要求的稀疏系数，y为所述超像素；argminx||x||1subjecttoDx＝y(2)公式(2)使用SLEP工具箱进行求解，获得x的解；使用公式(3)获得所述超像素的分类结果，式中xi表示第i类的稀疏系数，i＝1,2,…2KN；ri(y)＝||y-Dxi||2(3)根据公式(3)计算得到2KN个ri(y)，当r(y)的值最小时，此时的类别就是所述超像素的类别；对每幅图像的所有所述超像素分类，即可得到全局空间结构先验图。基于所述全局空间结构先验图，利用高斯概率密度函数计算局部相似结构先验，如下式所示：M(a,b)＝exp(-(cor(a,b)-c)2/u2)(4)式中，cor(a,b)是所述全局空间结构先验图中每个点的坐标，c是所述全局空间结构先验图中脉络膜新生血管的中心，u作为所述全局空间结构先验图中脉络膜的半径，是根据中心c与脉络膜边界点的距离的平均值求得，求得的矩阵M即为所述结构先验矩阵。使用所述结构先验矩阵M对图像进行转换，转换公式如公式(5)：Is＝MI0(5)式中，I0是原图像，Is是显著性增强图像。步骤S05中所述融合方法采用最大投票准则。本专利技术所达到的有益效果：结构先验矩阵能够获取patch之间的结构相关性信息，利用结构先验学习方法来计算图像中的结构先验矩阵，可用于增强来自脉络膜新生血管区域和背景区域的patch之间的相关性，多尺度分析能够获取不同尺度的有效上下文信息，进一步提高OCT图像中脉络膜新生血管的分割精度。附图说明图1是基于多尺度结构先验卷积神经网络的脉络膜新生血管分割流程图；图2是典型的卷积神经网络框架图；图3是脉络膜新生血管分割示例图。具体实施方式下面结合附图对本专利技术作进一步描述。以下实施例仅用于更加清楚地说明本专利技术的技术方案，而不能以此来限制本专利技术的保护范围。实施例1一种OCT图像中脉络膜新生血管分割算法，如图1所示，本方法主要包括2个阶段：训练阶段和测试阶段，具体步骤如下：(1)训练阶段，主要包括“结构先验学习”和“多尺度结构先验卷积神经网络训练”两部分，具体步骤如下：S01：对训练图像设计一种结构先验学习方法，构建结构先验矩阵，所述结构先验矩阵用于区分脉络膜新生血管区域和背景区域；所述结构先验学习方法包括以下步骤：a：对训练图像进行超像素分割得到若干超像素区域，采用SLIC算法进行分割；b：提取特征，所述特征包括每个所述超像素的平均灰度值、基于共生矩阵的纹理特征以及局部灰度特征；c：根据所述训练图像的正确标注将每个所述超像素标记为2类，分别为脉络膜新生血管区域和背景区域；d：使用标记好的若干所述超像素构造字典，使用K-means算法构造词典，假设在训练集中有N个病人的数据，每个数据包括2分类，使用K-means聚成K类，则N个病人的数据共可聚成2KN类，获得2KN个聚类中心，所述聚类中心组成字典D，如公式(1)所示：D＝[C1,1,C1,2…C1,K,B1,1,B1,2,..B1,K,…Cn,k..Bn,k…CN,1,..CN,K,BN,1,..BN,K](1)式中，Cn,k表示第n个病人的来自脉络膜新生血管区域的第k个聚类中心；Bn,k表示第n个病人的来自背景区域的第k个聚类中心；e：对所有所述超像素进行分类，得到全局结构先验，使本文档来自技高网...

【技术保护点】
一种OCT图像中脉络膜新生血管分割算法，其特征是，包括以下步骤：S01：对训练图像设计一种结构先验学习方法，构建结构先验矩阵，所述结构先验矩阵用于区分脉络膜新生血管区域和背景区域，包括以下步骤：a：对训练图像进行超像素分割得到若干超像素；b：提取特征；c：根据所述训练图像的正确标注将每个所述超像素标记为2类，分别为脉络膜新生血管区域和背景区域；d：使用标记好的若干所述超像素构造字典；e：对所有所述超像素进行分类，得到全局结构先验图；f：基于所述全局结构先验图，计算局部相似结构先验，求得所述结构先验矩阵；S02：基于所述结构先验矩阵将OCT原图像转换为显著性增强图像，用于增强脉络膜新生血管区域的显著性；S03：在所述显著性增强图像上使用多尺度分析，将所述显著性增强图像划分为m个尺度；S04：基于每种尺度训练得到m个训练好的卷积神经网络模型；S05：利用S01中步骤a、b、e、f对测试图像进行超像素分割、提取特征、分类以及计算结构先验矩阵，利用S02基于S05中得到的所述结构先验矩阵对测试图像进行图像转换，利用S03将转换后的所述测试图像划分为m个尺度，利用S04中训练好的所述卷积神经网络模型进行测试，输出m个分割结果，对m个分割结果进行融合即为最终的分割结果。...

【技术特征摘要】
1.一种OCT图像中脉络膜新生血管分割算法，其特征是，包括以下步骤：S01：对训练图像设计一种结构先验学习方法，构建结构先验矩阵，所述结构先验矩阵用于区分脉络膜新生血管区域和背景区域，包括以下步骤：a：对训练图像进行超像素分割得到若干超像素；b：提取特征；c：根据所述训练图像的正确标注将每个所述超像素标记为2类，分别为脉络膜新生血管区域和背景区域；d：使用标记好的若干所述超像素构造字典；e：对所有所述超像素进行分类，得到全局结构先验图；f：基于所述全局结构先验图，计算局部相似结构先验，求得所述结构先验矩阵；S02：基于所述结构先验矩阵将OCT原图像转换为显著性增强图像，用于增强脉络膜新生血管区域的显著性；S03：在所述显著性增强图像上使用多尺度分析，将所述显著性增强图像划分为m个尺度；S04：基于每种尺度训练得到m个训练好的卷积神经网络模型；S05：利用S01中步骤a、b、e、f对测试图像进行超像素分割、提取特征、分类以及计算结构先验矩阵，利用S02基于S05中得到的所述结构先验矩阵对测试图像进行图像转换，利用S03将转换后的所述测试图像划分为m个尺度，利用S04中训练好的所述卷积神经网络模型进行测试，输出m个分割结果，对m个分割结果进行融合即为最终的分割结果。2.根据权利要求1所述的一种OCT图像中脉络膜新生血管分割算法，其特征是，使用SLIC算法进行超像素分割。3.根据权利要求2所述的一种OCT图像中脉络膜新生血管分割算法，其特征是，所述特征包括每个所述超像素的平均灰度值、基于共生矩阵的纹理特征以及局部灰度特征。4.根据权利要求1所述的一种OCT图像中脉络膜新生血管分割算法，其特征是，使用K-means算法构造词典，假设在训练集中有N个病人的数据，每个数据包括2分类，使用K-means聚成K类，则N个病人的数据共可聚成2KN类，获得2KN个聚类中心，所述聚类中心组成字典D，如公式(1)所示：D＝[C1,1...

【专利技术属性】
技术研发人员：陈新建，袭肖明，
申请(专利权)人：苏州比格威医疗科技有限公司，
类型：发明
国别省市：江苏,32