CN106372593A - 一种基于血管收敛的视盘区定位方法 - Google Patents
一种基于血管收敛的视盘区定位方法 Download PDFInfo
- Publication number
- CN106372593A CN106372593A CN201610776889.1A CN201610776889A CN106372593A CN 106372593 A CN106372593 A CN 106372593A CN 201610776889 A CN201610776889 A CN 201610776889A CN 106372593 A CN106372593 A CN 106372593A
- Authority
- CN
- China
- Prior art keywords
- optic disc
- blood vessel
- disc area
- method based
- center
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/10—Human or animal bodies, e.g. vehicle occupants or pedestrians; Body parts, e.g. hands
- G06V40/18—Eye characteristics, e.g. of the iris
- G06V40/197—Matching; Classification
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/24—Aligning, centring, orientation detection or correction of the image
- G06V10/245—Aligning, centring, orientation detection or correction of the image by locating a pattern; Special marks for positioning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/26—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion
- G06V10/267—Segmentation of patterns in the image field; Cutting or merging of image elements to establish the pattern region, e.g. clustering-based techniques; Detection of occlusion by performing operations on regions, e.g. growing, shrinking or watersheds
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/30—Noise filtering
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Ophthalmology & Optometry (AREA)
- Human Computer Interaction (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本发明涉及一种基于血管收敛的视盘区定位方法,用于对眼底图像的视盘区进行定位,所述方法包括下列步骤:利用最小二乘分类器分割出眼底图像的初始视盘区;根据血管收敛位置计算分割出的初始视盘区的中心;计算初始视盘区的中心的各个方向上梯度值最大的点,构成点集,利用快速椭圆拟合算法处理点集,得到视盘区边界并确定视盘区。与现有技术相比,本发明具有鲁棒性强、能克服各种干扰因素、可以准确定位视盘区中心以及定位清晰明了等优点。
Description
技术领域
本发明涉及图像处理技术领域,尤其是涉及一种基于血管收敛的视盘区定位方法。
背景技术
视盘在眼底图像中具有较明显的特征:(1)通常表现为近似圆形的黄色或白色亮斑;(2)血管从视盘中心延伸出来;(3)延伸出的血管在方向上具有相似性。
1997年,Zhengliu等人在发表的“Automatic image analysis of fundusphotograph”中选取眼底图像中最大亮度值的2%作为阈值,所留下的区域即视盘所在位置;Chrastek等人在发表的“Optic disc segmentation in retinal images”中选取图像中亮度均值最大的点作为视盘中心。这两种方法只是单纯利用视盘区域最亮这一特征,当图像中存在亮度不均、渗出病变或视盘亮度丢失等情况时,定位的准确性受到极大影响。
2002年,Sinthanayothin等人在发表的“Automated location of the opticdisc,fovea,and retinal blood vessels from digital color fundus images”中则利用较亮的视盘中存在较暗的血管这一特征,计算各像素一定区域的亮度方差,将平均亮度方差最大的一点作为视盘的近似中心位置。很显然,这样可以排除部分高亮度病变或光照不均等得影响。Osareh等人在发表的“Comparison of color spaces for optic disclocation in retinal images”中利用视盘亮度和形状都具有相似性这一特点提出了基于模板匹配的方法,模板匹配后相关系数最大的点即视盘中心。
上面四种方法都只是利用视盘区域内的亮度和形状特征,不能很好的应用于视盘存在病变、视盘亮度丢失等情况的眼底图像中。另外一类方法是基于视盘与血管结构关系的。
2003年,Hoover等人在发表的“Locating the optic nerve in a retinal imageusing the fuzzy convergence of the blood vessels”中利用血管汇集于视盘中心这一特征,提出了一种计算像素点的血管汇集程度的算法。这种方法的缺点在于它易在血管分叉的地方给出错误的定位而且计算复杂度高。Forachia等人在发表的“Detection ofoptic disc in retinal images by means of a geometrical model of vesselstructure”中不仅利用血管汇集特征,并且综合考虑了整个血管***的结构方向特征,即主血管从视盘延伸出来的路径近似于两条相切的几何抛物线。最终采用几何模型方法实现视盘定位。由于光照不均、病变、血管干扰等干扰因素,这些方法不能准确定位视盘。
发明内容
本发明的目的是针对上述问题提供一种基于血管收敛的视盘区定位方法。
本发明的目的可以通过以下技术方案来实现:
一种基于血管收敛的视盘区定位方法,用于对眼底图像的视盘区进行定位,所述方法包括下列步骤:
1)利用最小二乘分类器分割出眼底图像的初始视盘区;
2)根据血管收敛位置计算步骤1)中分割出的初始视盘区的中心;
3)计算步骤2)中初始视盘区的中心的各个方向上梯度值最大的点,构成点集,利用快速椭圆拟合算法处理点集,得到视盘区边界并确定视盘区。
所述最小二乘分类器包括CSK***。
所述步骤1)具体为:
11)将眼底图像划分为多个大小为M×N的图像块;
12)对步骤11)中划分得到的每一个图像块中的每一个像素pm,n进行预处理,得到处理后的像素G(m,n),其中1≤m≤M,1≤n≤N;
13)利用最小二乘分类器遍历所有图像块,将图像块与视盘区模板进行匹配检测,并根据步骤12)中得到的G(m,n)计算每个图像块的检测分数;
14)选取最大检测分数对应的图像块,即为分割出的眼底图像的初始视盘区。
所述处理后的像素G(m,n)具体为:
其中,σ为正态分布的标准差。
所述检测分数具体为:
其中,St为每一个图像块的检测分数,为离散傅里叶变换,为离散傅里叶反变换,κ为点积操作,a(m,n)具体为:
其中,p为图像块,λ为正则化参数。
所述步骤2)具体为:
21)对步骤1)中得到的初始视盘区依次进行中值滤波和二值化血管分割,完成血管的提取;
22)利用垂直的矩形窗口从左到右遍历初始视盘区,结合提取的血管确定视盘区中心的横坐标;
23)利用水平的矩形窗口从上到下遍历初始视盘区,结合提取的血管确定视盘区中心的纵坐标。
所述初始视盘区的中心的横坐标为血管聚集程度D(v)最小时对应的垂直窗口v的中心横坐标,所述血管聚集程度D(v)具体为:
其中,n为窗口的数量,Mi为第个血管连接区域的像素个数,M为垂直窗口中像素的总个数。
所述初始视盘区的中心的纵坐标为Gabor滤波器强度值R(h)最大时对应的水平窗口h的中心纵坐标,所述血管聚集程度R(h)具体为:
其中,n为窗口的数量,(x,y)为像素的坐标,l(x,y)为坐标为(x,y)的像素的RGB值,g(x,y)为第h个窗口的Gabor滤波器的中值,Np为第h个水平窗口中的像素个数。
所述步骤3)具体为:
31)以步骤2)中得到的初始视盘区的中心为端点,向四周20个等间距的角度作20条射线;
32)分别求出20条射线上梯度值最大的点,构成点集;
33)利用快速椭圆拟合算法处理点集,得到视盘区边界并确定视盘区。
所述梯度值具体为:
G(i,j)=I(i,j)-I(i,j+1)
其中,G(i,j)为第i个方向上第j个像素点的梯度值,I(i,j)为第i个方向上第j个像素点的亮度值,I(i,j+1)为第i个方向上第j+1个像素点的亮度值。
与现有技术相比,本发明具有以下有益效果:
(1)对视盘区域检测方法是基于最小二乘分类器和模板匹配,鲁棒性强,能克服光照不均、病变、血管干扰等干扰因素。
(2)利用血管从视盘中心延伸出来的形态学特点,通过提取血管并计算血管收敛位置能准确定位视盘区中心。
(3)通过计算各个方向上梯度值最大的点并椭圆模拟出视盘边界,清晰明了。
(4)采用CSK***作为最小二乘分类器,相较于其他最小二乘分类器,CSK***速度最快,提高了视盘区域检测的速度。
(5)在分割视盘区时对每一个图像块中的每一个像素pm,n进行预处理,防止了眼底图像局部信息的扩散。
(6)首先通过模板匹配确定初始视盘区,再通过血管收敛确定视盘区的中心,既克服了模板匹配情况下会受到病变和视盘亮度丢失的情况,又克服了血管收敛时在血管分叉的地方给出错误定位这种情况,因此整体定位精度高效果好。
附图说明
图1为本发明的方法流程图;
图2为生成视盘区模板的实例图,其中(2a)为眼底图像,(2b)为生成的视盘区模板;
图3为视盘区边界椭圆模拟示意图,其中(3a)为点集图,(3b)为利用快速椭圆拟合算法处理点集得到的视盘区边界图;
图4为不同眼底图像检测视盘区的实例演示图,其中(4a)为眼底图像1,(4b)为眼底图像2,(4c)为眼底图像3,(4d)为眼底图像4,(4e)为眼底图像5,(4f)为眼底图像6;
图5为健康的眼底图像的视盘区检测结果和手工检测结果对比图,其中(5a)为眼底图像7,(5b)为眼底图像8,(5c)为眼底图像9,(5d)为眼底图像10;
图6为不健康的眼底图像的视盘区检测结果和手工检测结果对比图,其中(6a)为眼底图像11,(6b)为眼底图像12,(6c)为眼底图像13,(6d)为眼底图像14。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1所示,本实施例提供了一种基于血管收敛的视盘区定位方法,用于对眼底图像的视盘区进行定位,该方法包括下列步骤:
1)利用最小二乘分类器分割出眼底图像的初始视盘区;
2)根据血管收敛位置计算步骤1)中分割出的初始视盘区的中心;
3)计算步骤2)中初始视盘区的中心的各个方向上梯度值最大的点,构成点集,利用快速椭圆拟合算法处理点集,得到视盘区边界并确定视盘区。
上述步骤的具体步骤如下:
1)使用最小二乘分类器分割出视盘区域。
如图2所示,本实施例中最小二乘分类器使用的是Henriques J F等人在“Exploiting the circulant structure of tracking-by-detection with kernels”中提出的CSK tracker,这个分类器是最快的特征分类器。分割视盘区域是基于模板匹配的方法,左眼和右眼的视盘区模板是根据Marc Lalonde等人在“Fast and Robust Optic DiscDetection Using Pyramidal Decomposition and Hausdorff-Based TemplateMatching”提出的方法提取出来的。本发明通过对比每一个图像块和模板来判断每一个图像块是否是视盘区区域,图像块的大小是基于模板大小的(假设图像块大小为Mx N,1≤m≤M,1≤n≤N),在我们方法中,图像块的大小是80x80。为了防止眼底图像局部信息的扩散,对图像块中的每一个像素pm,n使用一个高斯函数,函数表达式如下:
式中G(m,n)表示每一个像素pm,n处理之后的信息,σ代表正态分布的标准差,在我们的方法中σ的值取0.2。
本发明在对比每一个图像块和模板的过程中计算最大检测分数,检测分数最大的图像块即为视盘区域,检测分数计算公式为:
式中St表示每一个图像块的检测分数,“x”操作为点乘,表示离散傅里叶变化,离散傅里叶反变换,κ是点积操作,定义如下:
a(m,n)系数定义如下:
式中p表示图像块,λ是正则化参数,λ>0,其他变量定义如上。
在最小二乘分类器遍历一张眼底图像的每一个图像块时,计算检测分数,检测分数最大的图像块就认为是视盘区。
2)通过计算血管收敛位置确定视盘区中心。
因为血管从视盘中心延伸出来,所以通过提取血管并计算血管收敛位置即可确定视盘区中心。确定视盘中心包括三个步骤:血管提取、确定横坐标和确定纵坐标。
RGB图像的亮度信息和血管信息在R通道和G通道,本发明选取RGB图像的G通道进行血管提取。在提取出来的视盘区域中使用中值滤波器,该中值滤波器使用直方图数据和边缘增强技术。然后再使用二值化技术进行血管分割,二值化的阈值是从其他眼底图像得到的,在本发明中二值化阈值为150。通过上述操作,即可提取出血管。
在提取出血管之后,通过一个垂直的矩形窗口从左到右遍历视盘区域来确定视盘中心的横坐标,窗口的宽度是血管宽度的两倍(血管宽度从模板中得到),窗口的高度和视盘区的高度相同。D(v)用来量化血管的聚集程度,定义为:
式中D(v)表示血管在第v个垂直窗口的聚集程度,n表示窗口的数量,Mi表示第i个血管连接区域的像素个数,M表示垂直窗口v中像素的总个数。当越小,是绝对值越大的负数。当D(v)最小时,视盘中心的横坐标就在第v个垂直窗口。视盘中心的横坐标因此确定。
在提取视盘横坐标之后,通过一个水平的矩形窗口从上到下遍历视盘区域来确定视盘中心的纵坐标,窗口的宽度和视盘区的宽度相同,窗口的高度和血管的宽度相同。首先使用中值滤波器去除不均匀的亮度信息和图像伪影,本发明使用Gabor滤波器,Gabor滤波器的中值为:
式中,g(x,y)表示第h个窗口的Gabor滤波器的中值,δ表示标准差,λ表示Gabor滤波器的波长,表示Gabor滤波器的相位偏移,γ表示椭圆率,x'和y'定义如下:
x'=x×cosθ+y×sinθ
y'=-x×sinθ+y×cosθ
式中(x,y)表示像素坐标,θ是Gabor滤波器的弧度。本方法中θ的值设置为π/2。
R(h)用来衡量水平窗口中Gabor滤波器值的强度,定义为:
式中l(x,y)表示(x,y)处像素的RGB值,Np表示第h个水平窗口中的像素个数。视盘中心的亮度值最大,Gabor滤波器值也最大,R(h)也最大,我们取R(h)最大的水平窗口作为视盘中心所在窗口,视盘中心的纵坐标由此确定。
3)计算各个方向上梯度值最大的点,对这些点使用快速椭圆拟合算法,拟合出视盘边界。
如图3所示,在得到比较清晰连续的眼底局部图像后我们需要将视盘的边界检测出来。首先我们以第二步中检测出的视盘中心为端点,向四周20个等间距的角度作20条射线,检测第i条射线上的第j个点的梯度值。如下:
G(i,j)=I(i,j)-I(i,j+1)
其中G(i,j)是第i个方向上第j个像素点的梯度值,I(i,j)表示第i个方向第j个像素点的亮度值,I(i,j+1)表示第i个方向上j+1个像素点的亮度值。分别求出20个方向上梯度值最大的点,因为这些点基本上是在视盘边界上,故将其作为我们后面拟合工作的候选点集合(用seeds表示)。得到边界点集合seeds之后,我们运用“最小二乘法”的思想用椭圆拟合这些点,即可得到我们需要的确定的视盘边界。
依据上述步骤,对于我们收集的眼底图像进行分析,我们选取Digital RetinalImages for Vessel Extraction(DRIVE)数据集和Non-fluorescein Images for VesselExtraction(NIVE)数据集。所有试验均在PC计算机上实现,该PC计算机的主要参数为:中央处理器Intel(R)Core(TM)2Duo CPU [email protected],内存2GB。
定位结果如图4~图6所示,结果显示,本发明对两个数据库(DRIVE和NIVE)中图片进行视盘区定位的精准度为97.65%。准确度方面明显优于其他算法。与此同时本发明算法还有计算量小的优点,处理每一张眼底图像的响应时间为1.6126s。相对于其他用了复杂的数学模型和数学变换的方法,我们的算法简单高效,效率高,鲁棒性强。综合看来我们的算法是先进的。
Claims (10)
1.一种基于血管收敛的视盘区定位方法,用于对眼底图像的视盘区进行定位,其特征在于,所述方法包括下列步骤:
1)利用最小二乘分类器分割出眼底图像的初始视盘区;
2)根据血管收敛位置计算步骤1)中分割出的初始视盘区的中心;
3)计算步骤2)中初始视盘区的中心的各个方向上梯度值最大的点,构成点集,利用快速椭圆拟合算法处理点集,得到视盘区边界并确定视盘区。
2.根据权利要求1所述的基于血管收敛的视盘区定位方法,其特征在于,所述最小二乘分类器包括CSK***。
3.根据权利要求1所述的基于血管收敛的视盘区定位方法,其特征在于,所述步骤1)具体为:
11)将眼底图像划分为多个大小为M×N的图像块;
12)对步骤11)中划分得到的每一个图像块中的每一个像素pm,n进行预处理,得到处理后的像素G(m,n),其中1≤m≤M,1≤n≤N;
13)利用最小二乘分类器遍历所有图像块,将图像块与视盘区模板进行匹配检测,并根据步骤12)中得到的G(m,n)计算每个图像块的检测分数;
14)选取最大检测分数对应的图像块,即为分割出的眼底图像的初始视盘区。
4.根据权利要求3所述的基于血管收敛的视盘区定位方法,其特征在于,所述处理后的像素G(m,n)具体为:
其中,σ为正态分布的标准差。
5.根据权利要求3所述的基于血管收敛的视盘区定位方法,其特征在于,所述检测分数具体为:
其中,St为每一个图像块的检测分数,为离散傅里叶变换,为离散傅里叶反变换,κ为点积操作,a(m,n)具体为:
其中,p为图像块,λ为正则化参数。
6.根据权利要求1所述的基于血管收敛的视盘区定位方法,其特征在于,所述步骤2)具体为:
21)对步骤1)中得到的初始视盘区依次进行中值滤波和二值化血管分割,完成血管的提取;
22)利用垂直的矩形窗口从左到右遍历初始视盘区,结合提取的血管确定视盘区中心的横坐标;
23)利用水平的矩形窗口从上到下遍历初始视盘区,结合提取的血管确定视盘区中心的纵坐标。
7.根据权利要求6所述的基于血管收敛的视盘区定位方法,其特征在于,所述初始视盘区的中心的横坐标为血管聚集程度D(v)最小时对应的垂直窗口v的中心横坐标,所述血管聚集程度D(v)具体为:
其中,n为窗口的数量,Mi为第个血管连接区域的像素个数,M为垂直窗口中像素的总个数。
8.根据权利要求7所述的基于血管收敛的视盘区定位方法,其特征在于,所述初始视盘区的中心的纵坐标为Gabor滤波器强度值R(h)最大时对应的水平窗口h的中心纵坐标,所述血管聚集程度R(h)具体为:
其中,n为窗口的数量,(x,y)为像素的坐标,l(x,y)为坐标为(x,y)的像素的RGB值,g(x,y)为第h个窗口的Gabor滤波器的中值,Np为第h个水平窗口中的像素个数。
9.根据权利要求1所述的基于血管收敛的视盘区定位方法,其特征在于,所述步骤3)具体为:
31)以步骤2)中得到的初始视盘区的中心为端点,向四周20个等间距的角度作20条射线;
32)分别求出20条射线上梯度值最大的点,构成点集;
33)利用快速椭圆拟合算法处理点集,得到视盘区边界并确定视盘区。
10.根据权利要求9所述的基于血管收敛的视盘区定位方法,其特征在于,所述梯度值具体为:
G(i,j)=I(i,j)-I(i,j+1)
其中,G(i,j)为第i个方向上第j个像素点的梯度值,I(i,j)为第i个方向上第j个像素点的亮度值,I(i,j+1)为第i个方向上第j+1个像素点的亮度值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610776889.1A CN106372593B (zh) | 2016-08-30 | 2016-08-30 | 一种基于血管收敛的视盘区定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610776889.1A CN106372593B (zh) | 2016-08-30 | 2016-08-30 | 一种基于血管收敛的视盘区定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106372593A true CN106372593A (zh) | 2017-02-01 |
CN106372593B CN106372593B (zh) | 2019-12-10 |
Family
ID=57902051
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610776889.1A Active CN106372593B (zh) | 2016-08-30 | 2016-08-30 | 一种基于血管收敛的视盘区定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106372593B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107292835A (zh) * | 2017-05-31 | 2017-10-24 | 瑞达昇科技(大连)有限公司 | 一种眼底图像视网膜血管自动矢量化的方法及装置 |
CN108205395A (zh) * | 2018-01-16 | 2018-06-26 | 安徽慧视金瞳科技有限公司 | 一种标定点中心坐标精准定位方法 |
CN112712521A (zh) * | 2021-01-18 | 2021-04-27 | 佛山科学技术学院 | 一种基于全局梯度搜索的眼底视盘自动定位方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2011022783A1 (en) * | 2009-08-28 | 2011-03-03 | Centre For Eye Research Australia | Feature detection and measurement in retinal images |
CN102112044A (zh) * | 2008-05-14 | 2011-06-29 | 科学、技术与研究机构 | 自动杯盘比测量*** |
CN102843957A (zh) * | 2009-08-24 | 2012-12-26 | 新加坡保健服务集团有限公司 | 一种用于检测视盘出血的方法及*** |
CN102842136A (zh) * | 2012-07-19 | 2012-12-26 | 湘潭大学 | 一种综合血管分布和视盘外观特性的视盘投影定位方法 |
US20130222767A1 (en) * | 2011-10-12 | 2013-08-29 | Jun Cheng | Methods and systems for detecting peripapillary atrophy |
CN104463140A (zh) * | 2014-12-23 | 2015-03-25 | 天津工业大学 | 一种彩色眼底图像视盘自动定位方法 |
CN104545792A (zh) * | 2014-12-31 | 2015-04-29 | 浙江大学 | 眼底图像的动静脉视网膜血管视盘定位方法 |
CN105718919A (zh) * | 2016-02-02 | 2016-06-29 | 吉林大学 | 一种眼底图像视盘定位方法及*** |
EP3048580A1 (en) * | 2015-01-20 | 2016-07-27 | Ulma Innovacion, S.L. | Method of extracting the optic disc of a retinal image |
-
2016
- 2016-08-30 CN CN201610776889.1A patent/CN106372593B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102112044A (zh) * | 2008-05-14 | 2011-06-29 | 科学、技术与研究机构 | 自动杯盘比测量*** |
CN102843957A (zh) * | 2009-08-24 | 2012-12-26 | 新加坡保健服务集团有限公司 | 一种用于检测视盘出血的方法及*** |
WO2011022783A1 (en) * | 2009-08-28 | 2011-03-03 | Centre For Eye Research Australia | Feature detection and measurement in retinal images |
US20130222767A1 (en) * | 2011-10-12 | 2013-08-29 | Jun Cheng | Methods and systems for detecting peripapillary atrophy |
CN102842136A (zh) * | 2012-07-19 | 2012-12-26 | 湘潭大学 | 一种综合血管分布和视盘外观特性的视盘投影定位方法 |
CN104463140A (zh) * | 2014-12-23 | 2015-03-25 | 天津工业大学 | 一种彩色眼底图像视盘自动定位方法 |
CN104545792A (zh) * | 2014-12-31 | 2015-04-29 | 浙江大学 | 眼底图像的动静脉视网膜血管视盘定位方法 |
EP3048580A1 (en) * | 2015-01-20 | 2016-07-27 | Ulma Innovacion, S.L. | Method of extracting the optic disc of a retinal image |
CN105718919A (zh) * | 2016-02-02 | 2016-06-29 | 吉林大学 | 一种眼底图像视盘定位方法及*** |
Non-Patent Citations (6)
Title |
---|
ANA SALAZAR-GONZALEZ 等: "Segmentation of the Blood Vessels and Optic Disk in Retinal Images", 《IEEE JOURNAL OF BIOMEDICAL AND HEALTH INFORMATICS》 * |
FLORIN ROTARU 等: "Optic Disc Identification Methods for Retinal Images", 《COMPUTER SCIENCE JOURNAL OF MOLDOVA》 * |
JOA˜O F. HENRIQUES 等: "Exploiting the Circulant Structure of Tracking-by-detection with Kernels", 《ECCV 2012》 * |
张先杰,张贵英: "眼底图像中视盘的自动定位方法研究", 《电脑知识与技术》 * |
肖志涛 等: "基于眼底结构特征的彩色眼底图像视盘定位方法", 《中国生物医学工程学报》 * |
高玮玮 等: "一种基于视网膜主血管方向的视盘定位及提取方法", 《激光生物学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107292835A (zh) * | 2017-05-31 | 2017-10-24 | 瑞达昇科技(大连)有限公司 | 一种眼底图像视网膜血管自动矢量化的方法及装置 |
CN107292835B (zh) * | 2017-05-31 | 2020-03-13 | 瑞达昇医疗科技(大连)有限公司 | 一种眼底图像视网膜血管自动矢量化的方法及装置 |
CN108205395A (zh) * | 2018-01-16 | 2018-06-26 | 安徽慧视金瞳科技有限公司 | 一种标定点中心坐标精准定位方法 |
CN108205395B (zh) * | 2018-01-16 | 2021-03-23 | 安徽慧视金瞳科技有限公司 | 一种标定点中心坐标精准定位方法 |
CN112712521A (zh) * | 2021-01-18 | 2021-04-27 | 佛山科学技术学院 | 一种基于全局梯度搜索的眼底视盘自动定位方法 |
CN112712521B (zh) * | 2021-01-18 | 2023-12-12 | 佛山科学技术学院 | 基于全局梯度搜索的眼底视盘自动定位方法及其存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN106372593B (zh) | 2019-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US12035971B2 (en) | Automated determination of arteriovenous ratio in images of blood vessels | |
Sadeghi et al. | Detection and analysis of irregular streaks in dermoscopic images of skin lesions | |
Kovács et al. | A self-calibrating approach for the segmentation of retinal vessels by template matching and contour reconstruction | |
CN102525381B (zh) | 图像处理装置、图像处理方法及计算机可读取的记录装置 | |
Liao et al. | Automatic segmentation for cell images based on bottleneck detection and ellipse fitting | |
CN103048329B (zh) | 一种基于主动轮廓模型的路面裂缝检测方法 | |
CN109409190A (zh) | 基于梯度直方图和Canny边缘检测器的行人检测方法 | |
Hatanaka et al. | Vertical cup-to-disc ratio measurement for diagnosis of glaucoma on fundus images | |
CN108615239B (zh) | 基于阈值技术和灰度投影的舌图像分割方法 | |
CN104766046A (zh) | 一种利用交通标识颜色及形状特征进行检测与识别算法 | |
CN106204555A (zh) | 一种结合Gbvs模型和相位一致性的视盘定位方法 | |
Fiorini et al. | Automatic Generation of Synthetic Retinal Fundus Images. | |
CN103914708A (zh) | 基于机器视觉的食品品种检测方法及*** | |
CN104408711A (zh) | 一种基于多尺度区域融合的显著区域检测方法 | |
CN104036521A (zh) | 一种新的视网膜眼底图像分割方法 | |
CN108186051A (zh) | 一种从超声图像中自动测量胎儿双顶径长度的图像处理方法及处理*** | |
CN106408566A (zh) | 一种胎儿超声图像质量控制方法及*** | |
CN106682678A (zh) | 一种基于支撑域的图像角点检测和分类方法 | |
CN110310323A (zh) | 基于Hessian矩阵和二维高斯拟合的视网膜血管管径测量方法 | |
CN104331877A (zh) | 融合色调梯度的彩色图像边缘检测方法 | |
CN106297492A (zh) | 一种教育玩具套件及利用颜色和轮廓识别编程模块的方法 | |
CN106372593A (zh) | 一种基于血管收敛的视盘区定位方法 | |
CN104574312A (zh) | 标靶图像求圆心的方法和装置 | |
CN106355596B (zh) | 一种融合均匀颜色信息和复合感受野模型的边缘检测方法 | |
CN103279960B (zh) | 一种基于x射线背散射图像的人体隐藏物的图像分割方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |