视网膜图像的视盘定位方法
【专利摘要】本发明公开了一种视网膜图像的视盘定位方法，包括以下步骤：a.擦除背景中的明亮斑块，同时保留血管特征；b.基于血管汇聚趋势特征确定视盘中心的候选位置；c.基于血管方向特征从所述视盘中心的候选位置中确定视盘中心的最终位置。该方法能够解决视盘定位易受到背景中各类病损及其它噪声所带来的干扰的问题，提高复杂背景下的视盘定位的精准度和速度。
【专利说明】视网膜图像的视盘定位方法
【技术领域】
[0001]本发明涉及一种视网膜图像的视盘定位方法。
【背景技术】
[0002]基于视网膜图像的生理特征的定位，包括视盘定位、血管定位等，定位的精准度将在很大程度上影响到后续对生理特征位置相应的病损判别过程的效率和精度。目前，视盘定位方法的研究主要集中在寻找一种图像特征，该特征能够在视盘与背景干扰目标上有差异性较大的表现。目前提出的比较有效的特征包括视盘灰度分布特征、视盘几何形状特征、视盘训练样本统计特征等等。基于这些特征分析的算法能够在背景不太复杂的情况下对视盘进行较好的定位，尤其是在无病损存在的视网膜图像中，定位效果很好。然而，在视网膜存在病变的视网膜图像中，表现各异的多种病损分布在背景中，将极大地干扰视盘的正常定位。由于视网膜病变病损表现差异性很大，在多种病损的干扰下，采取单一特征无法很好地实现视盘的精准定位。另外，除了病损给视盘定位带来的干扰，还存在着其它的干扰目标，如形态各异的血管网以及拍摄过程产生的噪声等。因此，这种复杂背景给快速、准确地定位视盘带来了很大的困难。
【发明内容】
[0003]本发明的主要目的在于提供一种视网膜图像的视盘定位方法，解决视盘定位易受到背景中各类病损及其它噪声所带来的干扰的问题，提高复杂背景下的视盘定位的精准度和速度。
[0004]为实现上述目的，本发明采用以下技术方案:
[0005]一种视网膜图像的视盘定位方法，包括以下步骤:
[0006]a.擦除背景中的明亮斑块，同时保留血管特征；
[0007]b.基于血管汇聚趋势特征确定视盘中心的候选位置；
[0008]c.基于血管方向特征从所述视盘中心的候选位置中确定视盘中心的最终位置。
[0009]在优选的技术方案中:
[0010]步骤a中，采用最近邻差值与高斯滤波算法处理所述明亮斑块和所述血管特征。
[0011]步骤b中，采用形态学滤波对血管汇聚趋势特征进行分析以确定所述视盘中心的候选位置。
[0012]步骤c中，在所述视盘中心的候选位置上进行血管方向特征匹配滤波，将取得最小滤波响应的位置确定为视盘中心的最终位置。
[0013]步骤a包括:
[0014](I)基于形态学重构的明亮物体检测，并将明亮物体的灰度值置为零得到填充种子图像Iseed ；
[0015](2)在填充种子图像Iseed灰度值为零的地方使用高斯滤波进行图像平滑；
[0016](3)图像平滑后在填充种子图像Iseed灰度值为零的地方进行最近邻插值得到填充图像Ipad ；
[0017](4)取填充图像Ipad与原始图像Ig的最小值作为填充结果Ipad =min (Ipad, Ig)；
[0018](5)重复(2) - (4)直至填充完毕。
[0019]步骤b包括:
[0020](I)提取视网膜血管并骨骼化，得到血管中心图；
[0021](2)去除所述血管中心图中所有的分支点与端点；
[0022](3)沿血管方向延长所有分支点与端点；
[0023](4)对延长后的血管中心图进行形态学膨胀；
[0024](5)对膨胀后的血管中心图进行二值化并与初始血管中心图相交，得到所述视盘中心的候选位置。
[0025]步骤c包括:
[0026]使用多个不同尺度的血管方向匹配滤波器，以多尺度滤波响应值累计和最小的位置作为所述视盘中心的最终位置。
[0027]步骤c中，使用4个不同尺度的血管方向匹配滤波器。
[0028]步骤c中，所述血管方向匹配滤波器是基于视盘内与周边血管方向的统一特征构造的。
[0029]本发明的有益技术效果:
[0030]本发明针对复杂背景下的视盘定位容易受到背景中各类病损及其它噪声带来的干扰问题，提供了一种复杂背景下的基于多特征融合的视盘分级定位方法，显著提高了视盘定位的效率和精度，能够迅速、准确的定位视盘。本发明基于视盘血管汇聚特征与视盘血管方向特征，有效地结合视盘的多种特征，由粗到精进行视盘分级定位。结合视盘的多种特征构造的特征集在视盘与背景干扰目标上有良好的互斥性。本发明不仅有效的结合了视盘的多个特征，还具有较低的运算复杂度，实现了复杂背景下快速、准确的盘定位。
[0031]本发明的优点主要体现于以下方面:
[0032]1.通过擦除背景中的明亮斑块，同时保留视盘中的血管特征，能够有效地去除视网膜图像中明亮病损对视盘定位的干扰；
[0033]2.通过视盘中心候选位置的快速粗定位，能够快速地得到视盘候选位置，保证了视盘定位的实时性；
[0034]3.在视盘粗定位的基础上进行血管方向匹配滤波，不仅大大提高了定位效率，而且有效地结合血管汇聚趋势与血管方向两个特征，提高了定位准确度，从而实现了视盘的快速、准确的定位。
【专利附图】
【附图说明】
[0035]图1是本发明视盘定位方法的流程图；
[0036]图2是去除明亮斑块同时保留血管特征的结果示意图；
[0037]图3是视盘候选位置粗定位的结果示意图；
[0038]图4是视盘最终位置精定位的结果示意图。【具体实施方式】
[0039]以下结合附图对本发明的实施例作详细说明。应该强调的是，下述说明仅仅是示例性的，而不是为了限制本发明的范围及其应用。
[0040]参阅图1，在一些实施例里，视网膜图像的视盘定位方法主要包括以下步骤。
[0041]1、去除明亮病损干扰
[0042]由于视网膜病变中的明亮病损和视盘均表现为明亮的黄色斑块，视盘地定位容易受到明亮病损的干扰。实施例中利用插值算法，快速有效地擦除背景中的明亮斑块，同时保留视盘中的血管特征，从而去除视网膜图像中明亮病损对视盘定位的干扰，并且保证在后面的步骤中以血管特征定位视盘。
[0043]优选地，插值算法主要包括以下步骤:(I)基于形态学重构的明亮物体检测，并将明亮物体的灰度值置为零得到填充种子图像Iseed ； (2)在Iseed灰度值为零的地方使用高斯滤波进行图像平滑；(3)图像平滑后在Iseed灰度值为零的地方进行最近邻插值得到填充图像Ipad ； (4)取填充图像Ipad与原始图像Ig的最小值作为填充结果Ipad =min(Ipad,Ig) ； (5)重复(2)-(4)直至填充完毕。图2所示为在一个示范例中采用本算法处理后的结果。可以看到，明亮物体被擦除，而视盘血管特征得到很好地保留。结合最近邻差值与高斯滤波对明亮病损进行擦除，可以使得擦除后的区域边缘梯度最小化，从而去除明亮病损对血管分割及血管方向分析的干扰。
[0044]2、视盘中心候选位置的快速粗定位
[0045]在明亮病损被去除的基础上，接下来可先进行视盘位置的粗定位。由于粗定位的精度不需要很高，优选采用基于形态学的视盘快速定位算法，可以保证定位的实时性，快速得到视盘候选位置。
[0046]优选地，视盘中心候选位置的快速粗定位过程主要包括以下步骤:(1)提取视网膜血管并骨骼化，得到血管中心图；(2)去除血管中心图中所有的分支点与端点；(3)沿血管方向延长所有分支点与端点；(4)对延长后的血管中心图进行形态学膨胀；(5)对膨胀后的血管中心图进行二值化并与初始血管中心图相交，得到视盘候选位置。图3所示为在一个示范例中采用快速粗定位处理后的结果。可以看到，视盘候选位置的像素数量很少。采用形态学滤波对血管汇聚趋势进行分析，有利于快速得到视盘中心候选位置。
[0047]3、视盘中心候选位置的精定位
[0048]在视盘候选位置粗定位的基础上，在视盘候选位置上使用血管方向匹配滤波器，实现视盘中心候选位置的精定位。
[0049]优选地，采用血管方向匹配滤波器，基于视盘内与周边血管方向的统一特征构造，该滤波器窗口尺寸与视盘尺寸为同一量级，仅在视盘候选位置上使用该滤波器，不仅大大提高了定位效率，并且有效地结合了血管汇聚趋势与血管方向两个特征，实现了视盘快速、准确的定位。通过在视盘候选位置上使用血管方向匹配滤波，可以将取得最小滤波响应的位置为视盘中心的最终位置。
[0050]在一个实施例中，可使用4个不同尺度的血管方向匹配滤波器，以多尺度滤波响应值累计和最小的位置作为视盘中心的最终定位位置。
[0051]图4所示为在一个示范例中视盘中心的最终定位位置。
[0052]以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属【技术领域】的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。
【权利要求】
1.一种视网膜图像的视盘定位方法，其特征在于，包括以下步骤: a.擦除背景中的明亮斑块，同时保留血管特征； b.基于血管汇聚趋势特征确定视盘中心的候选位置； c.基于血管方向特征从所述视盘中心的候选位置中确定视盘中心的最终位置。
2.如权利要求1所述的视网膜图像的视盘定位方法，其特征在于，步骤a中，采用最近邻差值与高斯滤波算法处理所述明亮斑块和所述血管特征。
3.如权利要求1所述的视网膜图像的视盘定位方法，其特征在于，步骤b中，采用形态学滤波对血管汇聚趋势特征进行分析以确定所述视盘中心的候选位置。
4.如权利要求1所述的视网膜图像的视盘定位方法，其特征在于，步骤c中，在所述视盘中心的候选位置上进行血管方向特征匹配滤波，将取得最小滤波响应的位置确定为视盘中心的最终位置。
5.如权利要求1至4任一项所述的视网膜图像的视盘定位方法，其特征在于，步骤a包括: (1)基于形态学重构的明亮物体检测，并将明亮物体的灰度值置为零得到填充种子图像 Iseed ； (2)在填充种子图像Iseed灰度值为零的地方使用高斯滤波进行图像平滑； (3)图像平滑后在填充种子图像Iseed灰度值为零的地方进行最近邻插值得到填充图像 Ipad ； (4)取填充图像Ipad与原始图像Ig的最小值作为填充结果Ipad= min(Ipad, Ig)； (5)重复(2)-(4)直至填充完毕。
6.如权利要求1至4任一项所述的视网膜图像的视盘定位方法，其特征在于，步骤b包括: (1)提取视网膜血管并骨骼化，得到血管中心图； (2)去除所述血管中心图中所有的分支点与端点； (3)沿血管方向延长所有分支点与端点； (4)对延长后的血管中心图进行形态学膨胀； (5)对膨胀后的血管中心图进行二值化并与初始血管中心图相交，得到所述视盘中心的候选位置。
7.如权利要求1至6任一项所述的视网膜图像的视盘定位方法，其特征在于，步骤c包括: 使用多个不同尺度的血管方向匹配滤波器，以多尺度滤波响应值累计和最小的位置作为所述视盘中心的最终位置。
8.如权利要求7所述的视网膜图像的视盘定位方法，其特征在于，步骤c中，使用4个不同尺度的血管方向匹配滤波器。
9.如权利要求4或7所述的视网膜图像的视盘定位方法，其特征在于，步骤c中，所述血管方向匹配滤波器是基于视盘内与周边血管方向的统一特征构造的。
【文档编号】G06T7/00GK103971369SQ201410204389
【公开日】2014年8月6日 申请日期:2014年5月14日 优先权日:2014年5月14日
【发明者】李钦, 卢瑞祥, 周迎春, 李名兆 申请人:深圳市计量质量检测研究院