CN102908120A - 眼底图像的配准方法、眼底图像中视盘神经、血管的测量方法及眼底图像的拼接方法 - Google Patents
眼底图像的配准方法、眼底图像中视盘神经、血管的测量方法及眼底图像的拼接方法 Download PDFInfo
- Publication number
- CN102908120A CN102908120A CN2012103808406A CN201210380840A CN102908120A CN 102908120 A CN102908120 A CN 102908120A CN 2012103808406 A CN2012103808406 A CN 2012103808406A CN 201210380840 A CN201210380840 A CN 201210380840A CN 102908120 A CN102908120 A CN 102908120A
- Authority
- CN
- China
- Prior art keywords
- fundus image
- eye fundus
- registration
- subject
- marginal point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Eye Examination Apparatus (AREA)
Abstract
本发明公开了一种眼底图像的配准方法、眼底图像中视盘神经、血管的测量方法及眼底图像的拼接方法。该配准方法包括:确定模板眼底图像和待配准眼底图像;在模板眼底图像上选定感兴趣区域;遍历构成待配准眼底图像内血管的所有边缘点,并同时计算sg值;通过判断sg值,确定配准的成功与失败。眼底图像中视盘神经、血管的测量方法及眼底图像的拼接方法均是基于配准方法实现的。本发明配准方法配准速度快、准度与精度高,基于配准方法实现的本发明视盘神经和血管的测量方法得到的测量结果精确度高,基于配准方法实现的本发明拼接方法可对多眼底图像进行精确拼接。
Description
技术领域
本发明涉及一种眼底图像的配准方法、眼底图像中视盘神经和血管的测量方法及眼底图像的拼接方法,属于眼底图像处理技术领域。
背景技术
随着科技的飞速发展,电子、印刷、医药、食品等领域中越来越广泛地应用到了机器视觉技术。需要提及的是,尽管机器视觉技术应用于医学图像分析的时间已很久,并有许多成功的应用范例,但是,机器视觉技术应用于眼科图像的研究到目前为止还是处于起步阶段,相当不完善。
目前,眼科图像研究主要是针对眼底图像的配准、眼底图像中的视盘神经和血管的测量以及眼底图像的拼接。
已有的配准方法主要有基于灰度的配准方法、基于边缘的配准方法、基于轮廓的配准方法。从实践中可以发现,基于灰度的配准方法只能容忍拍摄眼底时的线性光照变化以及刚体变换,而基于边缘或轮廓的配准方法存在一个主要缺陷,即需要提取图像中的眼底边缘,而这就导致了目标识别只能适用于非常小范围的光照变化,如果降低眼底图像的对比度,提取的边缘点就会越来越少,图像中的遮掩就会越来越多,因而在低对比度的眼底图像中便可能会找不到目标物体,而由于眼底图像的特性,特别是存有白内障的眼底图像,其对比度一般会比其他图像的对比度更低。总之,当面对眼底图像中的遮掩、几何畸变、模糊、严重非线性光照变化等干扰时,这些已有的配准方法便束手无策,或者,配准结果很不理想。
而眼底图像中的视盘神经和血管的测量以及眼底图像的拼接是基于眼底图像配准来实现的,因此,目前也没有出现可精确测量眼底图像中的视盘神经和血管的方法以及进行眼底图像之间的准确拼接的方法。
发明内容
本发明的目的在于提供一种眼底图像的配准方法、眼底图像中视盘神经和血管的测量方法及眼底图像的拼接方法,该配准方法配准速度快、准度与精度高,基于该配准方法实现的视盘神经和血管的测量方法得到的测量结果精确度高,基于该配准方法实现的拼接方法可对多眼底图像进行精确拼接。
为了实现上述目的,本发明采用了以下技术方案:
一种眼底图像的配准方法,其特征在于,它包括步骤:
步骤1:确定模板眼底图像和需要配准的待配准眼底图像;
步骤2:在该模板眼底图像上选定一个感兴趣区域;
步骤3:对构成该待配准眼底图像内血管的所有边缘点进行遍历,在遍历过程中,结合该感兴趣区域,对遍历到的各个该边缘点计算评价相似度的全局显性度量值sg;
步骤4:对于该待配准眼底图像中的血管的该边缘点,若存在全局显性度量值sg大于等于设定配准阈值的边缘点,则表明该待配准眼底图像与该模板眼底图像配准成功,若不存在全局显性度量值sg大于等于设定配准阈值的边缘点,则表明该待配准眼底图像与该模板眼底图像配准失败。
在所述步骤3中,所述待配准眼底图像中血管的每个所述边缘点对应的全局显性度量值sg的求取步骤为:
设定选定的所述感兴趣区域内的血管所具有边缘点的个数为n,所述感兴趣区域内的每个该边缘点的坐标为(ri,ci),i=1,2,...,n,所述感兴趣区域内的每个该边缘点的邻域均具有k个邻域点;
设定所述待配准眼底图像中血管的任一所述边缘点坐标为(r,c);
通过将所述待配准眼底图像中血管的各个所述边缘点坐标代入下式1)和2)中,计算出所述待配准眼底图像中血管的每个所述边缘点对应的全局显性度量值sg:
在式1)中,为所述模板眼底图像中的所述感兴趣区域内血管的第i个边缘点的第j个邻域点的横纵坐标基于广义霍夫变换构建的R-Table空间下归一化的梯度方向,j=1,2,...,k,为基于所述模板眼底图像中的所述感兴趣区域内血管的第i个边缘点进行的欧几里德变换,所述待配准眼底图像中血管的边缘点的点方向向量,TC为作用于模板眼底图像中感兴趣区域内血管的边缘点邻域的欧几里德变换。
一种基于所述的眼底图像的配准方法实现的眼底图像中视盘神经的测量方法,其特征在于,它包括步骤:
通过所述的眼底图像的配准方法找到与模板眼底图像配准成功的待配准眼底图像;
对该待配准眼底图像中的视盘神经进行测量,包括如下步骤:
将该待配准眼底图像等分成N个扇区;
在每个扇区中,利用二维度量模型,沿视盘方向寻找视盘神经的外轮廓和内轮廓的边缘,然后,在每个边缘位置用最小二乘法进行位姿的亚像素测量,从而完成对视盘神经的测量。
一种基于所述的眼底图像的配准方法实现的眼底图像中血管的测量方法,其特征在于,它包括步骤:
通过所述的眼底图像的配准方法找到与模板眼底图像配准成功的待配准眼底图像;
对该待配准眼底图像中的血管进行测量,包括如下步骤:
在该待配准眼底图像中,通过动脉血管、静脉血管的颜色与背景颜色的差别,以及利用直线的同极性特征,结合血管边缘梯度的变化,通过高斯模型与血管剖面的匹配,利用最小二乘法对血管边缘进行亚像素逼近,从而测量出血管的直径。
一种基于所述的眼底图像的配准方法实现的眼底图像的拼接方法,其特征在于,它包括步骤:
通过所述的眼底图像的配准方法找到与模板眼底图像配准成功的多个待配准眼底图像,该多个待配准眼底图像分别为从不同方位对一个人的眼球进行拍摄后得到的眼底图像;
将该多个待配准眼底图像映射到同一个空间中,并进行几何校正;
通过极性不同性特征,对每个该待配准眼底图像中的血管分支点进行定位;
通过定位好的血管分支点对该多个待配准眼底图像进行拼接,形成眼底整体图像。
本发明的优点是:
本发明眼底图像的配准方法利用广义hough(霍夫)变换和多尺度空间的不变特征实现了对模板眼底图像与待配准眼底图像间的配准,并且,当模板眼底图像与待配准眼底图像之间存在分辨率不同、尺寸不同、几何畸变等干扰,待配准眼底图像中存在遮挡、模糊等干扰时,本发明眼底图像的配准方法也能很好地完成配准工作,且配准速度快、准度与精度高。另外,K-NN最近邻分类法与随机抽样一致性可实现配准样本的筛选,提高了配准准度。
基于本发明眼底图像的配准方法实现的眼底图像中视盘神经的测量方法、眼底图像中血管的测量方法可以分别很好地实现对视盘神经、血管亚像素测量,测量结果精确度高。
基于本发明眼底图像的配准方法实现的眼底图像的拼接方法可精确地将一个人的各个不同方位角度的眼底图像进行拼接,且该拼接为一种带几何校正的拼接,是一种无缝拼接,拼接效果好。
附图说明
图1是本发明眼底图像的配准方法的实现流程图;
图2是本发明配准方法的第一实施例中模板眼底图像的示意图;
图3是本发明配准方法的第一实施例中待配准眼底图像的示意图;
图4是本发明配准方法的第一实施例中模板眼底图像中选定感兴趣区域的示意图;
图5是本发明配准方法的第一实施例中在待配准眼底图像上示出图4中感兴趣区域内血管的分布示意图;
图6是本发明配准方法的第二实施例中模板眼底图像的示意图;
图7是本发明配准方法的第二实施例中待配准眼底图像的示意图;
图8是本发明配准方法的第二实施例中模板眼底图像中选定感兴趣区域的示意图;
图9是本发明配准方法的第二实施例中在待配准眼底图像上示出图8中感兴趣区域内血管的分布示意图;
图10是对某个人的多个待配准眼底图像拼接后形成的一张眼底整体图像示意图;
图11是对某个人的多个待配准眼底图像拼接后形成的眼底整体图像进行颜色融合后得到的图像示意图。
具体实施方式
如图1所示,本发明眼底图像的配准方法包括如下步骤:
步骤1:确定模板眼底图像和需要配准的待配准眼底图像;
步骤2:在该模板眼底图像上选定一个感兴趣区域;
步骤3:对构成该待配准眼底图像内血管的所有边缘点进行遍历(对于该待配准眼底图像内血管边缘点的确定与选取为公知技术,可根据实际要求来确定和选取,不在这里详述),在遍历过程中,结合该感兴趣区域,对遍历到的各个该边缘点计算评价相似度的全局显性度量值sg;
步骤4:对于该待配准眼底图像中的血管的该边缘点,若存在全局显性度量值sg大于等于设定配准阈值的边缘点,则表明该待配准眼底图像与该模板眼底图像配准成功,即该待配准眼底图像与该模板眼底图像为同一个人的眼底图像,若不存在全局显性度量值sg大于等于设定配准阈值的边缘点,则表明该待配准眼底图像与该模板眼底图像配准失败,即该待配准眼底图像与该模板眼底图像不是同一个人的眼底图像。
在本发明配准方法中,需要说明的是,例如,模板眼底图像可以是某个人多年前拍摄的一张眼底图像,而待配准眼底图像是该人现在拍摄的一张眼底图像,两张眼底图像可能不是对眼球同一方位进行的拍摄,而且相对于多年前拍摄的那张眼底图像而言,现在拍摄的这张眼底图像可能因某些外在或内在的原因发生了刚体变换、扭曲、遮挡,或者变得模糊、不清晰,或者图片分辨率发生了变化、图片尺寸发生了变化等。
在步骤2中,该感兴趣区域为对配准有帮助的区域,一般来讲,只要不是仅包含少许血管的区域即可作为感兴趣区域,对该感兴趣区域的大小、形状是没有限制的。在实际中,可将整个模板眼底图像作为感兴趣区域。较佳地,该感兴趣区域中的血管应具有明显分布特点,以提高配准效率。在实际实施中,对于本发明配准方法,只要按上述要求选定一个感兴趣区域,即可得到准度与精度高的配准结果。当然,选定多个感兴趣区域也是可以的。
为了能够很好地理解求取待配准眼底图像中血管的每个边缘点对应的全局显性度量值sg的计算过程,下面先介绍一下存在刚体变换的两个眼底图像(模板图像与配准图像)间的配准。
首先,对于模板图像和配准图像,利用边缘检测结果来表示各自图像上具有的血管,图像上的血管目标只包括轮廓,这样,只要配准图像中存在任何灰度变化,就能将其表示为一个目标。
然后,在模板图像上选定一个感兴趣区域。
设定选定的该感兴趣区域内的血管具有n个边缘点,n为正整数,每个该边缘点的坐标为(ri,ci),i=1,2,...,n。设定配准图像中血管的任一边缘点的坐标为(r,c)。
于是,便可通过下式3)来进行配准:
公式3)为配准图像中血管的边缘点((r,c)为边缘点的坐标)对应的度量函数,在公式3)中,(r,c)为变量,根据公式3)即可求出配准图像中血管的各个边缘点的S值。
在式3)中,为模板图像中的感兴趣区域内血管的第i个边缘点的横纵坐标基于广义霍夫变换(GHT)构建的R-Table空间下归一化的梯度方向,为基于模板图像中的感兴趣区域内血管的第i个边缘点,配准图像中血管的边缘点的点方向向量,为与的点积,为与的欧几里德范数。
对于配准图像中血管的边缘点,若存在S值大于等于设定值的边缘点,则表明配准成功,反之,则配准失败。S值是一个介于0与1之间的数值,其中的设定值可根据配准要求而定,例如为0.6。
上述式3)示出的度量函数对非线性光照变化、局部遮挡和混乱具有很好的兼容性与不变性,理由如下:
第一,模板图像中的边缘点集与配准图像中的稠密梯度向量场相比较,即使有比较大的非线性光照变化,其梯度的幅值与方向向量的误差是一样的,而且一个滞后的阈值或非极值化抑制会完全消除任何光照变化带来的负面影响。
第二,配准图像中随机梯度方向的局部遮挡、噪声、混乱等虽会对度量函数的最大值产生轻微影响,但并不会使其位置发生偏移,因此,度量函数取值的意义就是模板图像中边缘点集配准成功的一个比率。
但是,在实际应用中,除了刚体变换,还存在很多其他非线性变换等。例如,扭曲形变可以看作是一个特殊、复杂的非线性变换。刚体变换只是扭曲形变的特例,是一个全局性的变换。上述度量函数无法适用于可能发生各种各样变化的配准图像与模板图像之间的配准工作。
因此,本发明配准方法采用了如下度量函数,具体为:
在步骤3中,待配准眼底图像中血管的每个边缘点对应的全局显性度量值sg的求取步骤为:
设定选定的该感兴趣区域内的血管所具有边缘点的个数为n(对于该感兴趣区域内血管边缘点的确定与选取为公知技术,不在这里详述),n为正整数,该感兴趣区域内的每个该边缘点的坐标为(ri,ci),i=1,2,...,n,该感兴趣区域内的每个该边缘点的邻域均具有k个邻域点,k为正整数;
设定该待配准眼底图像中血管的任一边缘点坐标为(r,c);
由此,通过将该待配准眼底图像中血管的各个边缘点坐标代入下式1)和2)中,计算出该待配准眼底图像中血管的每个边缘点对应的全局显性度量值sg:
对于公式1),sl(r,c,i)为针对模板眼底图像中感兴趣区域内血管的第i个边缘点,待配准眼底图像中血管的边缘点((r,c)为边缘点的坐标)对应的允许形变的局部隐性度量函数。对于公式2),sg(r,c)为待配准眼底图像中血管的边缘点((r,c)为边缘点的坐标)对应的全局显性度量函数。在公式1)和2)中,r、c为变量,根据公式1)和2)即可求出待配准眼底图像中血管的各个边缘点的sg值。sg为模板眼底图像上所有点局部贡献之和。sg值是一个介于0与1之间的数值。
在式1)中,为该模板眼底图像中的该感兴趣区域内血管的第i个边缘点的第j个邻域点的横纵坐标基于广义霍夫变换(GHT)构建的R-Table空间下归一化的梯度方向,j=1,2,...,k,为基于该模板眼底图像中的该感兴趣区域内血管的第i个边缘点进行的欧几里德变换,该待配准眼底图像中血管的边缘点的点方向向量,TC为作用于模板眼底图像中感兴趣区域内血管的边缘点邻域的欧几里德变换。该欧几里德变换包括各种形式的变化,如刚体变换、扭曲、遮挡、模糊、不清晰、分辨率变化、图片尺寸变化等,其属于公知技术。在公式1)中,为与的点积,为与的欧几里德范数。
对于公式1),其允许每一点的独立作用,即使整体发生了形变,模板眼底图像上的每一个点以及邻域并不会跟着发生相同的非线性形变,而且每一个点的变化可以通过一个局部欧几里德变换来近似。
在实际实施中,该设定配准阈值是根据配准要求而设定的,例如可为0.5、0.7等。
在实际实施中,在步骤4中:
若针对各个边缘点对应求出的全局显性度量值sg中的最大值大于等于0且小于设定配准阈值,则说明待配准眼底图像与模板眼底图像配准失败,不匹配;
若针对各个边缘点对应求出的全局显性度量值sg中的最大值大于等于设定配准阈值且小于1,则说明待配准眼底图像与模板眼底图像配准成功,但不完全匹配,待进一步确定;
若针对各个边缘点对应求出的全局显性度量值sg中的最大值等于1,则说明待配准眼底图像与模板眼底图像配准成功,且完全匹配。
在实际实施中,在所有的边缘点对中,或多或少会存在误配准的情况,因此,在步骤4中,对于待配准眼底图像中的血管的边缘点,若没有全局显性度量值sg等于1的边缘点且全局显性度量值sg大于等于设定配准阈值的边缘点的个数大于等于1,则通过K-NN最近邻分类法(公知技术)与随机抽样一致性(公知技术)对sg值大于等于设定配准阈值的边缘点进行筛选,筛除掉sg值计算错误的边缘点,然后再进行如下判断:
对于待配准眼底图像中的血管的边缘点,若筛选后存在全局显性度量值sg大于等于设定配准阈值的边缘点,则表明待配准眼底图像与模板眼底图像配准成功,即该待配准眼底图像与该模板眼底图像为同一个人的眼底图像,若筛选后不存在全局显性度量值sg大于等于设定配准阈值的边缘点,则表明该待配准眼底图像与该模板眼底图像配准失败,即该待配准眼底图像与该模板眼底图像不是同一个人的眼底图像。这种通过K-NN最近邻分类法与随机抽样一致性进行的筛选做法进一步保证了配准的稳定性、准度与精度。
同时,为了提高速度及保证边缘检测算法的稳定性,本发明配准方法还可引入多尺度不变空间,实行金字塔的搜索模式,这样可大大降低运算率,且由于多尺度不变空间的引入,可对不同强度的边缘都有一个很好的兼容。
在本发明配准方法中,需要提及的是,当待配准眼底图像中存在遮挡时,模板眼底图像中选定的感兴趣区域的面积至少为待配准眼底图像中的遮挡部分面积的2倍,这样才可保证配准结果的准度与精度。
若待配准眼底图像与模板眼底图像配准成功,且当待配准眼底图像中存在遮挡时,便可找到待配准眼底图像中的遮挡部分在模板眼底图像中的对应部分,从而对待配准眼底图像中的遮挡部分中的血管进行补全。
举例一:
图2为一模板眼底图像,图3为待配准眼底图像,图3中的眼底图像的拍摄角度与图2不同,相对图2而言,图3发生了透视变换,并且,图3中还存在模糊。根据配准要求,将设定配准阈值定为0.9。
通过本发明配准方法,选定感兴趣区域,如图4,白色长方形框示出的即为感兴趣区域,感兴趣区域内的血管边缘用白色线条表示。
对于图3,求出其存在sg值为0.98的边缘点,0.98大于设定配准阈值0.9,因此,图2示出的模板眼底图像与图3示出的待配准眼底图像为同一个人的眼底图像,如图5,图4中感兴趣区域中的血管在图5中的分布被清晰示出。
举例二:
图6为一模板眼底图像,图7为待配准眼底图像,相对图6而言,图7中的眼底图像发生了扭曲且有黑色遮挡。根据配准要求,将设定配准阈值定为0.8。
通过本发明配准方法,选定感兴趣区域,如图8,白色长方形框示出的即为感兴趣区域,感兴趣区域内的血管边缘用白色线条表示。
对于图7,求出其存在sg值为0.8的边缘点,0.8等于设定配准阈值0.8,因此,图6示出的模板眼底图像与图7示出的待配准眼底图像为同一个人的眼底图像。由于图7示出的待配准眼底图像中存在黑色遮挡,因此,通过本发明配准方法找到图7中遮挡部分在图6中的对应部分,对图7中的遮挡部分内的血管进行补全,遮挡部分内补全后的血管分布如图9所示。
由此可见,对于发生了刚体变换、扭曲、遮挡,或者变得模糊、不清晰,或者图片分辨率发生了变化、图片尺寸发生了变化等的待配准图像而言,通过本发明配准方法均能很好地进行配准工作。
本发明还提出了一种基于上述眼底图像的配准方法实现的眼底图像中视盘神经的测量方法,它包括如下步骤:
通过上述的眼底图像的配准方法找到与模板眼底图像配准成功的待配准眼底图像;
对该待配准眼底图像中的视盘神经进行测量,包括如下步骤:
将该待配准眼底图像等分成N个扇区(N为大于1的正整数,一般N为6);
在每个扇区中,利用二维度量模型(公知技术),沿视盘方向寻找视盘神经的外轮廓和内轮廓的边缘,然后,在每个边缘位置用最小二乘法进行位姿的亚像素测量,从而完成对视盘神经的测量。
本发明还提出了一种基于上述的眼底图像的配准方法实现的眼底图像中血管的测量方法,它包括如下步骤:
通过上述的眼底图像的配准方法找到与模板眼底图像配准成功的待配准眼底图像;
对该待配准眼底图像中的血管进行测量,包括如下步骤:
在该待配准眼底图像中,通过动脉血管、静脉血管的颜色与背景颜色的差别,以及利用直线的同极性特征,结合血管边缘梯度的变化,通过高斯模型(公知技术)与血管剖面的匹配,利用最小二乘法对血管边缘进行亚像素逼近,从而测量出血管的直径。
在实际实施上述两个本发明测量方法中,由于测量中不可避免地存在误差,因此可对某些测量进行修正,实现测量的准确性与稳定性。
另外,本发明还提出了一种基于上述的眼底图像的配准方法实现的眼底图像的拼接方法,它包括如下步骤:
通过上述的眼底图像的配准方法找到与模板眼底图像配准成功的多个待配准眼底图像,该多个待配准眼底图像分别为从不同方位对一个人的眼球进行拍摄后得到的眼底图像;
将该多个待配准眼底图像映射到同一个空间中,并进行几何校正,该几何校正可消除因角度带来的畸变;
通过极性不同性特征,对每个该待配准眼底图像中的血管分支点进行定位;
通过定位好的血管分支点对该多个待配准眼底图像进行拼接,形成眼底整体图像。
在实际实施中,对拼接后得到的眼底整体图像可进行颜色融合(公知技术),以消除拼接前各个待配准眼底图像的边界带来的阴影所造成的影响。
图10示出的是对某个人的多个待配准眼底图像拼接后形成的一张眼底整体图像,而图11示出的是对某个人的多个待配准眼底图像拼接后形成的一张眼底整体图像进行颜色融合后得到的图像。
在本发明中,配准方法、测量方法、拼接方法不仅可以用于人的眼底图像,还可用于动物的眼底图像。
本发明的优点是:
本发明眼底图像的配准方法利用广义hough(霍夫)变换和多尺度空间的不变特征实现了对模板眼底图像与待配准眼底图像间的配准,并且,当模板眼底图像与待配准眼底图像之间存在分辨率不同、尺寸不同、几何畸变等干扰,待配准眼底图像中存在遮挡、模糊等干扰时,本发明眼底图像的配准方法也能很好地完成配准工作,且配准速度快、准度与精度高。另外,K-NN最近邻分类法与随机抽样一致性可实现配准样本的筛选,提高了配准准度。
基于本发明眼底图像的配准方法实现的眼底图像中视盘神经的测量方法、眼底图像中血管的测量方法可以分别很好地实现对视盘神经、血管亚像素测量,测量结果精确度高。
基于本发明眼底图像的配准方法实现的眼底图像的拼接方法可精确地将一个人的各个不同方位角度的眼底图像进行拼接,且该拼接为一种带几何校正的拼接,是一种无缝拼接,拼接效果好。
上述是本发明的较佳实施例及其所运用的技术原理,对于本领域的技术人员来说,在不背离本发明的精神和范围的情况下,任何基于本发明技术方案基础上的等效变换、简单替换等显而易见的改变,均属于本发明保护范围之内。
Claims (9)
1.一种眼底图像的配准方法,其特征在于,它包括步骤:
步骤1:确定模板眼底图像和需要配准的待配准眼底图像;
步骤2:在该模板眼底图像上选定一个感兴趣区域;
步骤3:对构成该待配准眼底图像内血管的所有边缘点进行遍历,在遍历过程中,结合该感兴趣区域,对遍历到的各个该边缘点计算评价相似度的全局显性度量值sg;
步骤4:对于该待配准眼底图像中的血管的该边缘点,若存在全局显性度量值sg大于等于设定配准阈值的边缘点,则表明该待配准眼底图像与该模板眼底图像配准成功,若不存在全局显性度量值sg大于等于设定配准阈值的边缘点,则表明该待配准眼底图像与该模板眼底图像配准失败。
2.如权利要求1所述的眼底图像的配准方法,其特征在于:
在所述步骤3中,所述待配准眼底图像中血管的每个所述边缘点对应的全局显性度量值sg的求取步骤为:
设定选定的所述感兴趣区域内的血管所具有边缘点的个数为n,所述感兴趣区域内的每个该边缘点的坐标为(ri,ci),i=1,2,...,n,所述感兴趣区域内的每个该边缘点的邻域均具有k个邻域点;
设定所述待配准眼底图像中血管的任一所述边缘点坐标为(r,c);
通过将所述待配准眼底图像中血管的各个所述边缘点坐标代入下式1)和2)中,计算出所述待配准眼底图像中血管的每个所述边缘点对应的全局显性度量值sg:
3.如权利要求1所述的眼底图像的配准方法,其特征在于:
在所述步骤4中:
若针对各个所述边缘点对应求出的所述全局显性度量值sg中的最大值大于等于0且小于所述设定配准阈值,则说明所述待配准眼底图像与所述模板眼底图像配准失败,不匹配;
若针对各个所述边缘点对应求出的所述全局显性度量值sg中的最大值大于等于所述设定配准阈值且小于1,则说明所述待配准眼底图像与所述模板眼底图像配准成功,但不完全匹配;
若针对各个所述边缘点对应求出的所述全局显性度量值sg中的最大值等于1,则说明所述待配准眼底图像与所述模板眼底图像配准成功,且完全匹配。
4.如权利要求1所述的眼底图像的配准方法,其特征在于:
在所述步骤4中,对于所述待配准眼底图像中的血管的所述边缘点,若没有所述全局显性度量值sg等于1的所述边缘点且所述全局显性度量值sg大于等于所述设定配准阈值的所述边缘点的个数大于等于1,则通过K-NN最近邻分类法与随机抽样一致性对sg值大于等于所述设定配准阈值的所述边缘点进行筛选,筛除掉所述sg值计算错误的所述边缘点,然后进行如下判断:
对于所述待配准眼底图像中的血管的所述边缘点,若筛选后存在所述全局显性度量值sg大于等于所述设定配准阈值的边缘点,则表明所述待配准眼底图像与所述模板眼底图像配准成功,若筛选后不存在所述全局显性度量值sg大于等于所述设定配准阈值的边缘点,则表明该待配准眼底图像与该模板眼底图像配准失败。
5.如权利要求1或2或3或4所述的眼底图像的配准方法,其特征在于:
当所述待配准眼底图像中存在遮挡时,所述模板眼底图像中选定的所述感兴趣区域的面积至少为所述待配准眼底图像中的遮挡部分面积的2倍。
6.一种基于权利要求1所述的眼底图像的配准方法实现的眼底图像中视盘神经的测量方法,其特征在于,它包括步骤:
通过权利要求1所述的眼底图像的配准方法找到与模板眼底图像配准成功的待配准眼底图像;
对该待配准眼底图像中的视盘神经进行测量,包括如下步骤:
将该待配准眼底图像等分成N个扇区;
在每个扇区中,利用二维度量模型,沿视盘方向寻找视盘神经的外轮廓和内轮廓的边缘,然后,在每个边缘位置用最小二乘法进行位姿的亚像素测量,从而完成对视盘神经的测量。
7.一种基于权利要求1所述的眼底图像的配准方法实现的眼底图像中血管的测量方法,其特征在于,它包括步骤:
通过权利要求1所述的眼底图像的配准方法找到与模板眼底图像配准成功的待配准眼底图像;
对该待配准眼底图像中的血管进行测量,包括如下步骤:
在该待配准眼底图像中,通过动脉血管、静脉血管的颜色与背景颜色的差别,以及利用直线的同极性特征,结合血管边缘梯度的变化,通过高斯模型与血管剖面的匹配,利用最小二乘法对血管边缘进行亚像素逼近,从而测量出血管的直径。
8.一种基于权利要求1所述的眼底图像的配准方法实现的眼底图像的拼接方法,其特征在于,它包括步骤:
通过权利要求1所述的眼底图像的配准方法找到与模板眼底图像配准成功的多个待配准眼底图像,该多个待配准眼底图像分别为从不同方位对一个人的眼球进行拍摄后得到的眼底图像;
将该多个待配准眼底图像映射到同一个空间中,并进行几何校正;
通过极性不同性特征,对每个该待配准眼底图像中的血管分支点进行定位;
通过定位好的血管分支点对该多个待配准眼底图像进行拼接,形成眼底整体图像。
9.如权利要求8所述的拼接方法,其特征在于:
对拼接后得到的所述眼底整体图像进行颜色融合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210380840.6A CN102908120B (zh) | 2012-10-09 | 2012-10-09 | 眼底图像的配准方法、眼底图像中视盘神经、血管的测量方法及眼底图像的拼接方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210380840.6A CN102908120B (zh) | 2012-10-09 | 2012-10-09 | 眼底图像的配准方法、眼底图像中视盘神经、血管的测量方法及眼底图像的拼接方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102908120A true CN102908120A (zh) | 2013-02-06 |
CN102908120B CN102908120B (zh) | 2014-09-17 |
Family
ID=47606875
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210380840.6A Active CN102908120B (zh) | 2012-10-09 | 2012-10-09 | 眼底图像的配准方法、眼底图像中视盘神经、血管的测量方法及眼底图像的拼接方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102908120B (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104123711A (zh) * | 2013-04-24 | 2014-10-29 | 上海联影医疗科技有限公司 | 一种三维图像中多器官的定位方法 |
CN104933715A (zh) * | 2015-06-16 | 2015-09-23 | 山东大学(威海) | 一种应用于视网膜眼底图像的配准方法 |
CN105005998A (zh) * | 2015-08-05 | 2015-10-28 | 大连理工大学 | 一种基于多角度序列化图像空间特征点集的脑血管图像分割方法 |
CN105241894A (zh) * | 2015-08-28 | 2016-01-13 | 北京大恒图像视觉有限公司 | 一种用于多工序柔性印刷品检测的模板配准方法 |
CN105934193A (zh) * | 2013-12-23 | 2016-09-07 | Rsbv有限责任公司 | 宽场视网膜图像采集***和方法 |
CN106651827A (zh) * | 2016-09-09 | 2017-05-10 | 浙江大学 | 一种基于sift特征的眼底图像配准方法 |
CN106846293A (zh) * | 2016-12-14 | 2017-06-13 | 海纳医信(北京)软件科技有限责任公司 | 图像处理方法和装置 |
WO2018055545A1 (en) * | 2016-09-23 | 2018-03-29 | International Business Machines Corporation | Prediction of age related macular degeneration by image reconstruction |
CN108629769A (zh) * | 2018-05-02 | 2018-10-09 | 山东师范大学 | 基于最佳兄弟相似度的眼底图像视盘定位方法及*** |
CN108961334A (zh) * | 2018-06-26 | 2018-12-07 | 电子科技大学 | 一种基于图像配准的视网膜血管壁厚度测量方法 |
CN109447948A (zh) * | 2018-09-28 | 2019-03-08 | 上海理工大学 | 一种基于病灶彩色视网膜眼底图像的视盘分割方法 |
CN109767459A (zh) * | 2019-01-17 | 2019-05-17 | 中南大学 | 新型眼底图配准方法 |
CN110037651A (zh) * | 2018-01-15 | 2019-07-23 | 江威 | 眼底图像的质量控制方法及装置 |
CN110544274A (zh) * | 2019-07-18 | 2019-12-06 | 山东师范大学 | 一种基于多光谱的眼底图像配准方法及*** |
CN110689511A (zh) * | 2019-09-19 | 2020-01-14 | 华东师范大学 | 一种大视野图像拼接路径的优化方法 |
CN112070811A (zh) * | 2020-08-04 | 2020-12-11 | 广东工业大学 | 一种基于连续域蚁群算法改进的图像配准方法 |
CN112712521A (zh) * | 2021-01-18 | 2021-04-27 | 佛山科学技术学院 | 一种基于全局梯度搜索的眼底视盘自动定位方法 |
CN115409690A (zh) * | 2021-05-28 | 2022-11-29 | 南京博视医疗科技有限公司 | 一种眼底图像实时映射方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101593351A (zh) * | 2008-05-28 | 2009-12-02 | 中国科学院自动化研究所 | 基于距离变换和刚性变换参数估计的眼底图像配准方法 |
US20090323022A1 (en) * | 2008-06-26 | 2009-12-31 | Canon Kabushiki Kaisha | Medical imaging apparatus |
CN101984916A (zh) * | 2010-11-17 | 2011-03-16 | 哈尔滨工程大学 | 基于数字图像处理技术的血管管径测量方法 |
CN102122359A (zh) * | 2011-03-03 | 2011-07-13 | 北京航空航天大学 | 一种图像配准方法及装置 |
WO2011088578A1 (en) * | 2010-01-22 | 2011-07-28 | Annidis Health Systems Corp. | Registration method for multispectral retinal images |
-
2012
- 2012-10-09 CN CN201210380840.6A patent/CN102908120B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101593351A (zh) * | 2008-05-28 | 2009-12-02 | 中国科学院自动化研究所 | 基于距离变换和刚性变换参数估计的眼底图像配准方法 |
US20090323022A1 (en) * | 2008-06-26 | 2009-12-31 | Canon Kabushiki Kaisha | Medical imaging apparatus |
WO2011088578A1 (en) * | 2010-01-22 | 2011-07-28 | Annidis Health Systems Corp. | Registration method for multispectral retinal images |
CN101984916A (zh) * | 2010-11-17 | 2011-03-16 | 哈尔滨工程大学 | 基于数字图像处理技术的血管管径测量方法 |
CN102122359A (zh) * | 2011-03-03 | 2011-07-13 | 北京航空航天大学 | 一种图像配准方法及装置 |
Non-Patent Citations (1)
Title |
---|
黄琳琳,魏丽芳,潘林,余轮: "新的眼底图像配准方法的研究", 《计算机工程与应用》 * |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104123711B (zh) * | 2013-04-24 | 2016-03-30 | 上海联影医疗科技有限公司 | 一种三维图像中多器官的定位方法 |
CN104123711A (zh) * | 2013-04-24 | 2014-10-29 | 上海联影医疗科技有限公司 | 一种三维图像中多器官的定位方法 |
CN105934193A (zh) * | 2013-12-23 | 2016-09-07 | Rsbv有限责任公司 | 宽场视网膜图像采集***和方法 |
CN104933715A (zh) * | 2015-06-16 | 2015-09-23 | 山东大学(威海) | 一种应用于视网膜眼底图像的配准方法 |
CN105005998A (zh) * | 2015-08-05 | 2015-10-28 | 大连理工大学 | 一种基于多角度序列化图像空间特征点集的脑血管图像分割方法 |
CN105005998B (zh) * | 2015-08-05 | 2018-01-09 | 大连理工大学 | 一种基于多角度序列化图像空间特征点集的脑血管图像分割方法 |
CN105241894A (zh) * | 2015-08-28 | 2016-01-13 | 北京大恒图像视觉有限公司 | 一种用于多工序柔性印刷品检测的模板配准方法 |
CN105241894B (zh) * | 2015-08-28 | 2018-01-19 | 北京大恒图像视觉有限公司 | 一种用于多工序柔性印刷品检测的模板配准方法 |
CN106651827B (zh) * | 2016-09-09 | 2019-05-07 | 浙江大学 | 一种基于sift特征的眼底图像配准方法 |
CN106651827A (zh) * | 2016-09-09 | 2017-05-10 | 浙江大学 | 一种基于sift特征的眼底图像配准方法 |
GB2568623A (en) * | 2016-09-23 | 2019-05-22 | Ibm | Prediction of age related macular degeneration by image reconstruction |
WO2018055545A1 (en) * | 2016-09-23 | 2018-03-29 | International Business Machines Corporation | Prediction of age related macular degeneration by image reconstruction |
GB2568623B (en) * | 2016-09-23 | 2019-09-18 | Ibm | Prediction of age related macular degeneration by image reconstruction |
CN106846293A (zh) * | 2016-12-14 | 2017-06-13 | 海纳医信(北京)软件科技有限责任公司 | 图像处理方法和装置 |
CN110037651A (zh) * | 2018-01-15 | 2019-07-23 | 江威 | 眼底图像的质量控制方法及装置 |
CN108629769B (zh) * | 2018-05-02 | 2020-09-29 | 山东师范大学 | 基于最佳兄弟相似度的眼底图像视盘定位方法及*** |
CN108629769A (zh) * | 2018-05-02 | 2018-10-09 | 山东师范大学 | 基于最佳兄弟相似度的眼底图像视盘定位方法及*** |
CN108961334B (zh) * | 2018-06-26 | 2020-05-08 | 电子科技大学 | 一种基于图像配准的视网膜血管壁厚度测量方法 |
CN108961334A (zh) * | 2018-06-26 | 2018-12-07 | 电子科技大学 | 一种基于图像配准的视网膜血管壁厚度测量方法 |
CN109447948A (zh) * | 2018-09-28 | 2019-03-08 | 上海理工大学 | 一种基于病灶彩色视网膜眼底图像的视盘分割方法 |
CN109767459B (zh) * | 2019-01-17 | 2022-12-27 | 中南大学 | 新型眼底图配准方法 |
CN109767459A (zh) * | 2019-01-17 | 2019-05-17 | 中南大学 | 新型眼底图配准方法 |
CN110544274A (zh) * | 2019-07-18 | 2019-12-06 | 山东师范大学 | 一种基于多光谱的眼底图像配准方法及*** |
CN110544274B (zh) * | 2019-07-18 | 2022-03-29 | 山东师范大学 | 一种基于多光谱的眼底图像配准方法及*** |
CN110689511A (zh) * | 2019-09-19 | 2020-01-14 | 华东师范大学 | 一种大视野图像拼接路径的优化方法 |
CN110689511B (zh) * | 2019-09-19 | 2023-02-03 | 华东师范大学 | 一种大视野图像拼接路径的优化方法 |
CN112070811A (zh) * | 2020-08-04 | 2020-12-11 | 广东工业大学 | 一种基于连续域蚁群算法改进的图像配准方法 |
CN112070811B (zh) * | 2020-08-04 | 2023-07-28 | 广东工业大学 | 一种基于连续域蚁群算法改进的图像配准方法 |
CN112712521A (zh) * | 2021-01-18 | 2021-04-27 | 佛山科学技术学院 | 一种基于全局梯度搜索的眼底视盘自动定位方法 |
CN112712521B (zh) * | 2021-01-18 | 2023-12-12 | 佛山科学技术学院 | 基于全局梯度搜索的眼底视盘自动定位方法及其存储介质 |
CN115409690A (zh) * | 2021-05-28 | 2022-11-29 | 南京博视医疗科技有限公司 | 一种眼底图像实时映射方法及装置 |
CN115409690B (zh) * | 2021-05-28 | 2023-09-29 | 南京博视医疗科技有限公司 | 一种眼底图像实时映射方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN102908120B (zh) | 2014-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102908120B (zh) | 眼底图像的配准方法、眼底图像中视盘神经、血管的测量方法及眼底图像的拼接方法 | |
Kovács et al. | A self-calibrating approach for the segmentation of retinal vessels by template matching and contour reconstruction | |
US8103055B2 (en) | Detection of blobs in images | |
CN108022228A (zh) | 基于SIFT变换和Otsu匹配的彩色眼底图像拼接方法 | |
CN105913093A (zh) | 一种用于文字识别处理的模板匹配方法 | |
CN113826143A (zh) | 特征点检测 | |
CN107292869B (zh) | 基于各向异性高斯核和梯度搜索的图像斑点检测方法 | |
CN107066961B (zh) | 指纹配准方法及装置 | |
CN112164043A (zh) | 一种用于多眼底图像拼接的方法及*** | |
Lebeda et al. | Texture-independent long-term tracking using virtual corners | |
CN108447016A (zh) | 基于直线交点的光学图像和sar图像的匹配方法 | |
CN106682678A (zh) | 一种基于支撑域的图像角点检测和分类方法 | |
Martinikorena et al. | Fast and robust ellipse detection algorithm for head-mounted eye tracking systems | |
CN109657717A (zh) | 一种基于多尺度密集结构特征提取的异源图像匹配方法 | |
CN103632137A (zh) | 一种人眼虹膜图像分割方法 | |
Liu et al. | Retinal image registration via feature-guided Gaussian mixture model | |
CN110009670A (zh) | 基于fast特征提取和piifd特征描述的异源图像配准方法 | |
CN111127532B (zh) | 基于深度学习特征光流的医学图像形变配准方法及*** | |
CN109993090B (zh) | 基于级联回归森林和图像灰度特征的虹膜中心定位方法 | |
CN108320799A (zh) | 一种用于侧向流纸条疾病诊断的图像分析与识别方法 | |
Tang et al. | Retinal image registration based on robust non-rigid point matching method | |
Almasi et al. | Registration of fluorescein angiography and optical coherence tomography images of curved retina via scanning laser ophthalmoscopy photographs | |
AU2020103713A4 (en) | Digital imaging methods and system for processing agar plate images for automated diagnostics | |
Chen et al. | Eyes localization algorithm based on prior MTCNN face detection | |
CN110414521A (zh) | 一种变电站内变压器油位计示数识别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C41 | Transfer of patent application or patent right or utility model | ||
TR01 | Transfer of patent right |
Effective date of registration: 20160422 Address after: 100080 Suzhou Street Haidian District, No. 9, layer, 3 901-01 Patentee after: Beijing Daheng Puxin Medical Technology Co Ltd Address before: 100085 Beijing city Haidian District on the seven Street No. 1 Building 1 room No. 309 Huizhong Patentee before: Beijing Daheng Image Vision Co., Ltd. |