[0048] Specificity = TN/(TN+FP)
[0049] Accuracy = (TP+TN)/(TP+FN+TN+FP)
[0050] 式中，TP(True Positive)表示真阳性率，FP(False Positive)表示假阳性率， FN(False Negative)表示假阴性率，TN(True Negative)表示真阴性率。
[0051 ] 为了验证本方法，选取DRIVE数据库中的的20例眼底图像作为实验样本，对20例 眼底数码图像分别利用二维匹配滤波法，基于数学形态学的方法，区域增长法，自适应阈值 法，Snake模型，高斯一阶导数匹配滤波法，传统Gabor滤波法，基于分水岭的方法进行血管 分割，评价结果如表1所示；
[0052] 表1. 20例眼底数码图像不同血管分割结果评价均值
[0054] 由上表实验结果可知，本发明血管自动分割算法的特异性和准确率分别为0. 9777 和0. 939,高于现有的几种血管分割算法。说明本发明在实现血管全自动分割的同时，保证 了较好的分割特异性和准确性。
[0055] 本发明以眼底图像血管的自动分割为研究目的，着重阐述了利用Gabor滤波器和 ICA相结合的方法进行视网膜图像血管分割的算法。本发明实现了滤波器相关参数的自动 选取；在保持算法分割有效性的同时，且准确率最大的情况下，实现了阈值自动选取。
[0056] 以上所述，仅为本发明中的【具体实施方式】，但本发明的保护范围并不局限于此，任 何熟悉该技术的人在本发明所揭露的技术范围内，可理解想到的变换或替换，都应涵盖在 本发明的包含范围之内，因此，本发明的保护范围应该以权利要求书的保护范围为准。
【主权项】
1. 基于Gabor滤波器的眼底图像血管自动分割方法，本方法设计了基于Gabor滤波器 的眼底图像血管分割算法，在保证分割结果有效性的同时，有效减少了滤波器的方向，实现 了滤波器参数以及阈值的自动选取； 首先在预处理阶段选取血管对比度最优的绿色通道图像，利用20个不同角度的Gabor滤波器对眼底图像血管边缘进行粗提取，然后，以滤波器的最大响应作为血管的粗定位位 置，在滤波器参数选择中，引入优化的ICA算法实现了参数的自动选取，在二值化阈值选择 中，在准确率最大原则的基础上实现阈值的自动选择，最后，利用形态学处理后得到的图像 掩膜，去除分割边缘假阳性区域，得到血管的最终分割结果； 其特征在于：本方法的特征在于采取以下步骤， 步骤1，原始图像的预处理； 步骤2,利用20个间隔为Γ的Gabor滤波器去获取血管的边缘信息，以具有最大响应 的图像作为血管的粗分割结果； 步骤3,利用优化的ICA算法实现Gabor滤波器参数的最优化自动选取； 步骤4,利用准确率最大原则的系统研究法决定有效的阈值，对初分割图像进行二值 化； 步骤5,基于形态学的二值化图像后处理，得到眼底图像血管最终分割结果。2. 根据权利要求1所述的基于Gabor滤波器的眼底图像血管自动分割方法，其特征在 于：步骤1，读取原始图像，分别获取原始图像的红色（R)、绿色（G)、蓝色（B)通道图像，选 取三个通道图像中血管对比度最强的绿色通道为后续处理的输入图像； 步骤2,Gabor滤波器是一个用于边缘检测的线性滤波器；在空域，一个二维的Gabor滤 波器具有在空间域和频率域同时取得最优局部化的特性，与人类生物视觉特性很相似，因 此能够很好地描述对应于空间频率、空间位置及方向选择性的局部结构信息；Gabor滤波 器是自相似的，也就是说，所有Gabor滤波器都可以从一个母小波经过膨胀和旋转产生；实 际应用中，Gabor滤波器可以在频域的不同尺度，不同方向上提取相关特征；Gabor滤波器 的核函数如下：式中，g(x，y)为滤波器的幅值响应，〇y是χ和y方向的标准差，f。是曲线的频 率； 余弦曲线有一个宽度τ，其中仁=l/τ，L为宽度垂直平方向上的高度； 传统滤波器中，共有范围中的180个间隔Γ的不同角度的核函数被用于眼底 图像中血管的边缘提取；选取DRIVE数据库中的20幅图像，对不同方向个数的滤波器进行 分割结果灵敏性，特异性和准确率的分析实验，结果发现当角度间隔为Γ的Gabor滤波器 大于20个时，分割性能不再随着滤波器个数的增多而变得更优，因此本方法在保证系统性 能的前提下，将传统Gabor滤波器中的180个角度降低为20个；其中角度旋转过程中的的 坐标变换公式如下： xf =xcosΘ+ysinθ y1 = _xsin Θ +ycos θ 式中，（？，太）是每个旋转后的Θ值所对应的坐标；180个不同角度（Θ)的Gabor滤波器ge(X，y)作用于图像I(X，y)将得到一系列的滤波器响应Ge(X，y)，计算公式如下： Ge(X，y)=I(X，y)*ge(X，y) 式中，*代表卷积运算.为有效检测血管位置，在每一个像素点（x，y)，只有最大的响应 才会被保留；保留最大值后的Gabor滤波器的各个方向响应R(x，y)计算公式如下：随后，利用对数函数抑制全响应中的峰值，表达式如下：式中，为对数函数处理后的滤波器相应图像，Ln为对数运算； 步骤3,Gabor滤波器有很多参数，其中最重要的是空间的宽度τ和高度L，分割结果 的好坏很大程度上取决于这两个参数； 步骤4,将图像/0^；)的灰度值阈值化得到二值图像Ω(x，y)，二值化公式如下：式中，是二值化的阈值，二值化的结果是得到血管的分割结果.在阈值的确定过程 中，首先通过图像统计的方法确定能使每幅图像准确率达到最大化的阈值，然后计算所有 图像阈值的平均值为本方法中的阈值； 步骤5,为提取图像中的有效分割区域?(x，y)，本方法将步骤4中的二值化图像Ω(X，y)与一个形态学腐蚀操作后的二值化图像掩膜Μ(X，y)逐像素相乘来去除分割结果 中有效区域外围的假阳性部分，确保整幅图像中只保留眼底区域为有效的分割区域； Φ(x,y)=Ω (χ, y) [Μ (χ, y) Θ Β] 式中，Θ表示形态学的腐蚀操作，Β表示一个圆形的结构元。3. 根据权利要求2所述的基于Gabor滤波器的眼底图像血管自动分割方法，其特征在 于：步骤3中，首先，获得数据库中的训练集中每幅图像的最优参数，取所有图像最优参数 的平均值作为本方法的参数，其次将以上确定的参数用于图像测试集，本方法中数据库的 训练集和测试集是完全分开的；ICA优化方法分别在[0. 1，20](间隔为0. 1)和[0. 1，15] (间隔为〇. 1)范围对每幅图像依次实验获取其各个最优宽度和高度，而后得到各个最优宽 度平均值和高度平均值作为本方法的滤波器宽度和高度值。4. 根据权利要求2所述的基于Gabor滤波器的眼底图像血管自动分割方法，其 特征在于：为了评价本方法对于眼底图像血管分割结果的有效性，本方法引入敏感性 SE(Sensitivity)、特异性SP(Specificity)、精确度ACC(Accuracy)三个评价指标，具体计 算方法如下： Sensitivity=TP/(TP+FN) Specificity=TN/(TN+FP) Accuracy= (TP+TN)/(TP+FN+TN+FP) 式中，TP(TruePositive)表示真阳性率，FP(FalsePositive)表示假阳性率，FN(FalseNegative)表示假阴性率，TN(TrueNegative)表示真阴性率。
【专利摘要】基于Gabor滤波器的眼底图像血管自动分割方法，在保证分割结果有效性的同时，有效减少了滤波器的方向，实现了滤波器参数以及阈值的自动选取。算法首先在预处理阶段选取血管对比度最优的绿色通道图像，利用20个不同角度的Gabor滤波器对眼底图像血管边缘进行粗提取，然后，以滤波器的最大响应作为血管的粗定位位置，在滤波器参数选择中，本发明引入了优化的ICA算法实现了参数的自动选取，在二值化阈值选择中，在准确率最大原则的基础上实现阈值的自动选择，最后，利用形态学处理后得到的图像掩膜，去除分割边缘假阳性区域，得到血管的最终分割结果。
【IPC分类】G06T7/00
【公开号】CN105261015
【申请号】CN201510638912
【发明人】杨春兰, 法尔纳兹, 段彦华, 吴水才, 刘冰
【申请人】北京工业大学
【公开日】2016年1月20日
【申请日】2015年9月29日