一种应用于视网膜眼底图像的配准方法
【技术领域】
[0001] 本发明属于医学图像处理领域，具体涉及一种应用于视网膜眼底图像的配准方 法。
【背景技术】
[0002] 视网膜眼底图像是眼科医生常用的一种标准的客观诊断影像，它是由眼底照相机 获取。眼底是指内眼后部组织结构的统称，它一般包含视网膜、视神经、脉络膜和黄斑等。 研宄表明，青光眼、糖尿病视网膜病变和冠心病等病可导致视力下降，是导致眼盲的主要原 因。在临床诊疗的过程中，眼科医生需要对同一患者不同时期拍摄的视网膜眼底图像进行 对比分析。由于成像角度和成像环境的不同，这时候需要对两幅图像进行高精度的配准。另 外，由于眼底相机的视角问题，每幅拍摄的范围有限，就需要将多幅图像进行融合来观测眼 部是否出现病变区域。因此，需要对这些视网膜眼底图像进行配准，为后续的视网膜眼底图 像融合以及眼底病变检测来做准备。总而言之，视网膜眼底图像配准对视网膜病变的诊断 起到了重要作用，所以视网膜眼底图像的配准研宄对眼底疾病的早期发现、辅助诊断有着 重要的意义。
[0003] 基于SIFT(ScaleInvariantFeatureTransform)方法作为当前在医学图像的配 准中最广受关注的配准方法之一，其提取图像的特征点十分稳定，同时还解决了图像间发 生平移、旋转、仿射变换以及光照变换等情况下的配准问题。但是不可忽视的一点是，该方 法提取的特征不是人们视觉意义上的角点，在某些对比度高的医学图像中会检测出大量的 特征点，这样将会产生冗余信息，致使在特征描述的时候会严重影响效率；而在对比度比较 低的医学图像中，如视网膜眼底图像中，由于SIFT配准方法对特征点候选点的检测条件较 为严格，部分图像检测出来的特征点数量太少，影响配准的精度。
【发明内容】
[0004]针对经典的SIFT方法在部分视网膜眼底图像无法实现的配准缺陷，本发明设计 了一种应用于眼底视网膜图像的配准方法。本发明在配准过程中引入了新的特征点检测 及特征描述方法，提高眼底视网膜图像特征点检测数量和配准精度。为了进一步的提高配 准的精确度，在配准中引入误匹配点剔除策略。本发明的核心技术方案主要有以下两个方 面：
[0005] 1.引入了新的特征点检测方法及特征描述方法
[0006] 如附图1所示，一个特征结构包含主分叉点和它周边的三个邻居节点，其中，与主 分叉点相连的三个点之间的连线的长度分别为，角度为a a2,a3。定义角度描 述特征，a夕a5/a9,k2= a2/a3/a4,k3= a6/a7/a8,k4= a1(|/an/a12。对对h 和h进行归一化处理，
，则特征描述向量为：
[0007]S=[T1；T2, T3.kr!，kr2,kr3,kr4]
[0008] 2.误匹配点剔除策略的引入
[0009] 为了进一步提高图像配准精度，本发明采用参考图像特征点之间的欧氏距离的比 值与待配准图像相对应的特征点之间的欧氏距离比值相等的准则来剔除误配准点。这样可 以进一步提高配准的精度，误配准点剔除策略过程如下所示：
[0010] (1)初始化，特征配准点数量为Rmax，参考图像与待配准图像之间对应的特征点集 分别为M(i)、N(i)，令i= 1 ;
[0011] (2)0&七8(1)=(1(1(1)，(1+1))/(1"(1)，叭1+1))，其中(1代表的是两点之间欧氏 距离；
[0012] (3)统计出Data⑴中出现频率最高的数值h;
[0013] (4)计算Data⑴中等于h特征点集合M'，N'即为已经剔除的误匹配集合。
[0014] 本发明的有益效果在于：
[0015] 1.提高图像特征点检测的数量
[0016] 本发明能够相对较好的对眼底视网膜图像的进行特征点提取。附图2是SIFT方 法、Harris角点方法、SURF(SpeededUpRobustFeatures)方法和本发明的方法在不同眼 底视网膜图像的特征检测实验结果，可以看出本发明在视网膜眼底图像特征点检测数量上 要优于其它的三种方法。表1是四种方法在四幅不同的视网膜眼底图像中提取特征点数量 的对比。
[0017]表1四种方法的提取特征点数量对比[0018]
[0019] 2.提高配准的精确度
[0020] 本发明能够相对准确的实现对医学图像的配准。附图3是SIFT方法、Harris角 点方法、SURF方法和本发明的方法对眼底视网膜图像特征配准的实验结果，可以看出本发 明在视网膜眼底图像特征点配准对数和配准的精确度要优于其它三种方法。表2所示为选 取平移、旋转和光照变换后的眼底视网膜图像中四种方法在配准中误配准点和配准点组数 的对比。
[0021] 表2眼底图像配准中误配准点和配准点组数的客观评价
[0022]
[0023] 综上所述，本发明设计了一种应用于眼底视网膜图像的配准方法。该方法在SIFT 方法中引入了分叉结构特征描述，在配准点数量和配准精度上都有了较大的提升，能够很 好地实现眼底视网膜图像的配准。
【附图说明】
[0024] 图1为分叉结构示意图；
[0025] 图2为SIFT方法、Harris角点方法、SURF方法和本发明的方法在眼底视网膜图像 的特征检测实验结果；
[0026] 图3为SIFT方法、Harris角点方法、SURF方法和本发明的方法特征配准的实验结 果。
【具体实施方式】
[0027] 下面结合附图及具体实施例对本发明作进一步详细的说明。
[0028] 本发明实现医学图像的配准过程依次经过以下步骤：
[0029] 步骤1 :对两幅待配准的图像进行预处理，使其类型一致；
[0030]步骤2:提取图像中的SIFT特征结构点集，并进行特征描述；
[0031] 步骤3:提取图像中的分叉结构点集，并进行特征描述。特征描述的方法如 附图1所示，一个特征结构包含主分叉点和它周边的三个邻居节点，k1=ai/c^/ag， k2=a2/a3/a4,k3=a6/a7/a8,k4=a1(|/an/a12。对 1,ki进行归一化处理，
，则特征描述向量为：
[0032]S=[T1；r2,r3.kr!，kr2,kr3,kr4]
[0033] 步骤4:依据两幅图像得到的特征点集进行相似度度量，确定配准点对；
[0034] 步骤5 :剔除错误的配准点对，实现两幅图像特征点的配准。其中，利用参考图像 特征点之间的欧氏距离的比值与待配准图像相对应的特征点之间的欧氏距离比值相等的 准则来剔除错误配准点对。
[0035] 最后应说明的是，以上实施实例仅用以说明本发明的技术方案而非限制，尽管参 照较佳的实施实例对本发明进行了详细的说明，但本发明的保护范围并不局限于此，任何 熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的修改或等同替换， 而不脱离本发明技术方案的精神和范围，都应涵盖在本发明的保护范围之内。
【主权项】
1. 一种应用于视网膜眼底图像的配准方法，其特征在于包括w下步骤： 步骤1;对两幅待配准的图像进行预处理，使其类型一致； 步骤2 ;提取图像中的SIFT特征结构点集，并进行特征描述； 步骤3 ;提取图像中的分叉结构点集，并进行特征描述； 步骤4 ;依据两幅图像得到的特征点集进行相似度度量，确定配准点对； 步骤5 ;剔除错误的配准点对，实现两幅图像的配准。2. 根据权利要求1所述的一种应用于视网膜眼底图像的配准方法，其特征在于：在所 述步骤3中使用新的特征点W及特征描述方式，如附图1所示，一个特征结构包含主分叉 点和它周边的=个邻居节点，其中，与主分叉点A相连的=个点之间的连线的长度分别为 12, I3,角度为a。a2,a3,则a1+a5+a9=2n。那么k1=ai/as/a9是一个定值，同 理k2=曰2/曰3/曰4,k3=曰e/曰y/曰8,k4=曰1。/曰ii/曰。。对li和ki进行归一化处理，则特征描述向量为： s= [r i，r 2，1' 3，k' i，k' 2，k' 3，k' 4]3. 根据权利要求1所述的一种应用于视网膜眼底图像的配准方法，其特征在于：在所 述步骤5中引入误配准点剔除策略实现配准，利用参考图像特征点之间的欧氏距离的比值 与待配准图像相对应的特征点之间的欧氏距离比值相等准则来剔除误配准点，该样可W进 一步提高配准的精度。
【专利摘要】一种应用于视网膜眼底图像的配准方法，属于医学图像处理领域。在经典的SIFT（Scale Invariant Feature Transform）方法基础上，提出新的特征点检测及特征描述方法，提高眼底视网膜图像特征点检测数量和配准精度，使得图像配准效果更好。为了进一步的提高配准的精确度，在配准中引入误配准点剔除策略。经过大量的配准实验表明，本发明可以实现眼底视网膜图像的准确、稳定的配准，克服了SIFT方法在部分视网膜眼底图像中无法实现配准的缺陷。
【IPC分类】G06T7/00
【公开号】CN104933715
【申请号】CN201510333654
【发明人】董恩清, 孙华魁, 黄振强, 仲伟冲, 李宇森
【申请人】山东大学（威海）
【公开日】2015年9月23日
【申请日】2015年6月16日