一种新的视网膜眼底图像分割方法技术

本发明专利技术提出了一种新的视网膜血管网络分割算法。该算法的新颖点在于结合了多尺度线性检测与利用图像的灰度-梯度共生矩阵计算最佳熵阈值的方法。首先，提取视网膜眼底图像中包含血管轮廓信息丰富的绿色分量，对其进行阴影修正，降噪，CLAHE等预处理；其次，基于血管的形态结构特征，对视网膜图像血管采取多尺度多方向的线性检测，融合不同尺度下的图像响应，从而获得血管图像的特征；最后，基于图像的灰度—梯度共生矩阵，计算图像的最佳熵阈值进行分割。本算法分割准确度高，能提取更多细小的血管，计算速度快，具有很好的鲁棒性。本发明专利技术适用于正常及病变的视网膜眼底图像的分割。

【技术实现步骤摘要】
【专利摘要】本专利技术提出了一种新的视网膜血管网络分割算法。该算法的新颖点在于结合了多尺度线性检测与利用图像的灰度-梯度共生矩阵计算最佳熵阈值的方法。首先，提取视网膜眼底图像中包含血管轮廓信息丰富的绿色分量，对其进行阴影修正，降噪，CLAHE等预处理；其次，基于血管的形态结构特征，对视网膜图像血管采取多尺度多方向的线性检测，融合不同尺度下的图像响应，从而获得血管图像的特征；最后，基于图像的灰度—梯度共生矩阵，计算图像的最佳熵阈值进行分割。本算法分割准确度高，能提取更多细小的血管，计算速度快，具有很好的鲁棒性。本专利技术适用于正常及病变的视网膜眼底图像的分割。【专利说明】
本专利技术涉及计算机图像处理领域，具体为一种基于结合多尺度线性检测与利用灰度一梯度共生矩阵获取最佳熵阈值，进行视网膜眼底图像分割方法，尤其适合正常的视网膜眼底图像分割。
技术介绍
视网膜血管，是全身微循环系统的重要组成部分，其形态结构的变化与糖尿病，高血压，脑血管硬化和冠状动脉硬化等心血管疾病的病情严重程度密切相关。通过提取视网膜血管，分析其特征，例如，血管管径及弯曲度等，并进行相关参数的测量分析，就能够很大程度上对心血管疾病进行预测，从而实施科学的预防性干预和药物治疗。视网膜眼底图像中血管的分割算法还是不少的。Chaudhuri等人在论文((Detect1n of blood vessels in retinal images using two-dimens1nal matched filters))中提出了一种二维匹配滤波的方法，滤波后的血管结构得到增强，但是易丢失血管分叉点和细小的血管。Thitiporn 等人在论文((An efficient blood vessel detect1n algorithmfor retinal images using local entropy thresholding》中提出 了基于局部熵阈值的血管分割方法，方法简单快速，对正常的视网膜血管图像能实现较好的分割，对于病变图像分割效果差。Nguyen 等人在论文《An effective retinal blood vessel segmentat1n methodusing multi_scale line detect1n))中提出了多尺度线性检测器的方法对血管进行检测。此外，监督学习的方法也用于实现血管分割(You在论文《Segmentat1n of retinal bloodvessels using the radial project1n and sem1-supervised approach))中使用了支持向量机的方法，Soares 在论文《Retinal vessel segmentat1n using the2~dgabor waveletand supervised classificat1n》中使用贝叶斯分类器等)。多数方法采用已分割的血管和背景数据得到训练模型，然后用此模型对新图像的每一个像素点进行分类，判断该像素点是血管还是背景。Zana 等人在论文《Segmentat1n of vessel-like patterns usingmathematica lmorphology and curvature evaluat1n》中提出基于血管分段呈线性和分段连接的特点，使用具有线性结构元素的数字形态操作算子来增强血管和区别背景。上述血管分割算法仍然会存在一些缺陷，包括:血管中央反射区分割效果差，分叉点和交叉点的分割效果差；紧靠的血管容易合并，难分割；细小血管容易丢失；在视盘和病理区存在错误的血管检测。
技术实现思路
本专利技术的目的在于克服上述现有技术的不足之处，提出了，该方法能够正确分割紧邻的血管和存在中央反射区的血管，并且能够分割出更多细小的血管。除此之外，由于该方法有很好的鲁棒性。本专利技术是通过以下技术方案实现的，包括以下步骤:S1:对视网膜眼底图像进行图像预处理。彩色眼底图像为RGB图像，选取包含血管轮廓信息丰富且对比度较强的绿色分量进行血管分割。绿色分量的灰度图像光照不均匀，动态范围小，细节不够清晰，故需要对绿色分量图像进行预处理，具体操作包括:阴影修正，降噪，对比度受限的直方图均衡化处理(CLAHE)。采用OSTU阈值方法求解原始眼底图像的红色分量的阈值，根据阈值将图像进行二值化处理，从而得到二值图像。对二值图像执行形态学的“开”操作，“闭”操作及“腐蚀”操作，便可得到掩膜图像，该图像中，对原始眼底图像感兴趣的部分标记为1，其余背景部分标记为O。利用掩膜图像选取眼底图像感兴趣的部分，最后得到预处理后的图像。S2:利用第一步得到的预处理后的图像进行多尺度多方向线性检测。本文采用多尺度多方向的线性检测器对归一化的图像进行检测。较长的线性检测器用于检测中央反射区，较短的线性检测器可以避免出现紧靠的血管不易分割，因此，融合不同尺度的线性检测器便可以综合上述优点，达到很好的分割效果。S3:基于灰度一梯度共生矩阵，计算第二步得到的图像的最佳分割阈值。图像的灰度是构成一幅图像的基础，梯度是构成图像的边缘轮廓的要素。灰度一梯度共生矩阵反映了图像中灰度和梯度的分布规律，还能够表示各像素点与其邻域像素点的空间关系。S4:对第三步得到图像进行后处理。本专利技术选取八连通域对二值图像进行连通域标记，遵循从左到右，从上到下的规则对分割后的图像进行扫描。每一幅眼底图像都设定一个像素个数值，当标记的连通分量的面积小于该设定值时，删除该区域，否则，就保留该连通分量。【专利附图】【附图说明】为了使本专利技术的目的，技术方案及优点更加清楚，下面将结合附图对本专利技术做进一步的详尽阐述，其中:图1为本专利技术实施的总体流程框图。图2为预处理后的图像。图3为多尺度线性检测器图4为灰度-梯度共生矩阵。图5为提取出的血管。【具体实施方式】下面结合附图对本专利技术的优选实施例做详细说明，本实施例在以本专利技术技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但优选实例仅为了说明本专利技术，而不是为了限制本专利技术的保护范围。本实施采用的图像帧来自标准数据库。图1为本专利技术实施例提供的的系统总体框图，如图所示:该系统框图由4个功能模块组成:(I)图像预处理；(2)多尺度线性检测；(3)基于灰度-梯度共生矩阵，计算最佳熵阈值；(4)后处理图2为经过预处理后的图像。预处理主要包括以下步骤:Sll:提取视网膜眼底图像的绿色分量S12:阴影修正S13:降噪处理S14:对比度受限的直方图均衡化处理图3为检测窗口为15X15，检测线长度为11，方向为12个方向的线性检测器。该线性检测器的核心思想就是选取尺度为L的检测线，每隔15度旋转一次，共有12个方向，改变L的数值就可以构成多尺度线性检测器。由于视网膜血管的宽度集中在7-8个像素，因此设置W为15个像素，线性检测器的尺度范围为1-15，步长设置为2。图像的特征响应为=4ιΛ-(I)其中窗口宽度用W表示，L表示线性滤波器的尺度，η为窗口内的像素个数。【权利要求】1.一种新的视网膜眼底图像分割算法，结合了多尺度线性检测与利用图像的灰度-梯度共生矩阵计算最佳熵阈值的方法，其特征在于包括以下步骤: 51:对视网膜眼底图像进行图像预本文档来自技高网...

【技术保护点】
一种新的视网膜眼底图像分割算法，结合了多尺度线性检测与利用图像的灰度‑梯度共生矩阵计算最佳熵阈值的方法，其特征在于包括以下步骤： S1：对视网膜眼底图像进行图像预处理； S2：利用第一步得到的预处理后的图像进行多尺度多方向线性检测； S3：基于灰度—梯度共生矩阵，计算第二步得到的图像的最佳分割阈值； S4：对第三步得到图像进行后处理。

【技术特征摘要】

【专利技术属性】
技术研发人员：张思杰，翟丽红，邓蕊，戴阳，徐鹏，曾孝平，
申请(专利权)人：重庆大学，
类型：发明
国别省市：重庆;85