CN104077754B - 基于对称性的视网膜血管滤波增强方法 - Google Patents

基于对称性的视网膜血管滤波增强方法 Download PDF

Info

Publication number
CN104077754B
CN104077754B CN201410315974.9A CN201410315974A CN104077754B CN 104077754 B CN104077754 B CN 104077754B CN 201410315974 A CN201410315974 A CN 201410315974A CN 104077754 B CN104077754 B CN 104077754B
Authority
CN
China
Prior art keywords
image
retinal
scale
retina
filtering
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.)
Active
Application number
CN201410315974.9A
Other languages
English (en)
Other versions
CN104077754A (zh
Inventor
田捷
李国栋
杨彩云
徐敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Automation of Chinese Academy of Science
Original Assignee
Institute of Automation of Chinese Academy of Science
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Institute of Automation of Chinese Academy of Science filed Critical Institute of Automation of Chinese Academy of Science
Priority to CN201410315974.9A priority Critical patent/CN104077754B/zh
Publication of CN104077754A publication Critical patent/CN104077754A/zh
Application granted granted Critical
Publication of CN104077754B publication Critical patent/CN104077754B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)
  • Eye Examination Apparatus (AREA)

Abstract

本发明公开了一种基于对称性的视网膜血管滤波增强方法,包括步骤:对原始视网膜图像进行灰度化处理,再进行二维多尺度高斯滤波;利用滤波后的信息计算梯度合成向量;计算图像每一点在不同尺度下的梯度方向参数、每一点的梯度幅值参数、每一点的灰度参数;计算得到基于对称性的多尺度视网膜血管滤波增强函数,使用该血管函数对视网膜图像进行滤波,达到视网膜血管滤波增强的效果。本发明能对视网膜血管结构进行增强,同时可以有效抑制视盘等非血管结构,对视网膜图像血管分割、视网膜图像分析以及视网膜图像的配准中有重要的辅助处理和实际应用价值。

Description

基于对称性的视网膜血管滤波增强方法
技术领域
本发明涉及图像处理技术领域,特别涉及一种基于对称性的视网膜血管滤波增强方法。
背景技术
视网膜是一种重要的人体组织,是人体循环系统的重要组成部分,也是人体少数可以通过非侵入方式获得图像的组织结构。通过眼底照相技术,可以直接获取到视网膜血管的图像,如图2A、图3A、图4A和图5A所示,这些血管信息在全身很多疾病的诊断与治疗中起着至关重要的作用。作为人体循环系统的一个重要组成部分,视网膜血管网络形态结构的变化也可以间接了解到很多全身的其它疾病情况。通过了解视网膜血管的各种情况,例如动静脉管径变化、走行是否弯曲、交叉压迫的程度、反光带的宽窄、有无血管鞘和视网膜出血等,在一定程度上可以了解到脑、心、肾等重要器官的血管状况,为高血压的诊断、分期、治疗及疗效评估提供可靠的依据。因此,视网膜图像中血管的检测,有助于为疾病的早期诊断和治疗提供科学依据,具有很重要的临床使用价值。视网膜血管存在大量的噪声和病变组织,这对于视网膜血管的分割产生很大的影响,因此在对视网膜血管分割之前一个很关键的步骤是对血管图像进行滤波和增强,经过滤波后的图像再进行分割就会得到一个较好的结果,所以对于视网膜血管图像的滤波增强一直是一个重点和难点问题。本发明提出的基于对称性的视网膜血管滤波增强方法能够有效地对视网膜图像的血管结构起到滤波和增强的效果,同时有效的抑制视盘等非血管结构,为视网膜图像的后期处理起到一个关键的作用。
发明内容
由上可知,视网膜图像血管分割是视网膜处理的一个重要的组成部分, 对视网膜血管进行分割的很关键的一个步骤是对血管图像进行滤波增强,同时抑制非血管结构的响应,通过滤波增强后再对视网膜血管图像进行分割可以得到更好地分割效果。本发明的目的就在于提供一种基于对称性的视网膜血管滤波增强方法,使得视网膜图像的血管结构能够得到增强,同时有效抑制非血管结构的响应。
为达到上述目的,本发明提供了一种基于对称性的视网膜血管滤波增强方法,该方法包括以下步骤:
步骤S1:获得视网膜灰度图像I(x,y);
步骤S2:用二维多尺度高斯滤波器G(x,y,s)对于所述视网膜灰度图像I(x,y)进行滤波,得到高斯滤波图像VG(x,y,s),其中,x,y表示图像中某一点的坐标,s为对应的尺度;
步骤S3:基于得到的高斯滤波图像VG(x,y,s)计算旋转后的梯度合成向量T(x,y,s);
步骤S4:计算所述视网膜灰度图像I(x,y)中的每一点在不同尺度下的梯度方向参数D(x,y,s);
步骤S5:计算所述视网膜灰度图像I(x,y)中的每一点在不同尺度下的梯度幅值参数M(x,y,s);
步骤S6:计算所述视网膜灰度图像I(x,y)中的每一点在不同尺度下的灰度参数I(x,y,s);
步骤S7:根据所述梯度合成向量G(x,y,s)、梯度方向参数D(x,y,s)、梯度幅值参数M(x,y,s)、灰度参数I(x,y,s),计算得到基于对称性的多尺度视网膜血管滤波增强函数F(x,y,s),利用所述血管滤波增强函数对于所述视网膜灰度图像进行滤波,得到经过滤波和增强的视网膜图像。
本发明提供的基于对称性的视网膜血管滤波增强方法,是视网膜血管分割的预处理技术,通过对视网膜血管的滤波增强,可以使得视网膜血管的分割更加准确,减少噪声的干扰,对视网膜图像的后期处理,比如视网膜图像配准,视网膜血管跟踪等相关方面有重要的应用价值。
另外,本发明提供的基于对称性的视网膜血管滤波增强方法,能够为视网膜疾病的诊断提供最重要的视网膜血管分析结果,对青光眼、糖尿病、高血压等疾病的筛选和预防具有重要的作用。
综上,本发明旨在对视网膜图像的血管结构进行滤波和增强,以便于视网膜图像的分割和配准等后处理。实验结果显示,本发明方法对于视网膜图像的血管结构有着明显的增强效果,同时有效地抑制了视网膜图像的非血管结构。
附图说明
图1为本发明基于对称性的视网膜血管滤波增强方法的流程图;
图2A为眼底照相机拍摄的视网膜图像,图2B为根据本发明方法对图2A所示视网膜图像进行滤波增强得到的视网膜图像;
图3A为眼底照相机拍摄的另一幅视网膜图像,图3B为根据本发明方法对图3A所示视网膜图像进行滤波增强得到的视网膜图像;
图4A为眼底照相机拍摄的另一幅视网膜图像,图4B为根据本发明方法对图4A所示视网膜图像进行滤波增强得到的视网膜图像;
图5A为眼底照相机拍摄的另一幅视网膜图像,图5B为根据本发明方法对图5A所示视网膜图像进行滤波增强得到的视网膜图像。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
本发明提出一种基于对称性的视网膜血管滤波增强方法,该方法的基本思路为:用二维多尺度高斯滤波器对视网膜灰度图像I(x,y)进行滤波,得到高斯滤波图像VG(x,y,s);利用高斯滤波图像VG(x,y,s)计算旋转后的梯度合成向量G(x,y,s);计算图像I(x,y)每一点在不同尺度下的梯度方向参数D(x,y,s);计算图像I(x,y)每一点在不同尺度下的梯度幅值参数M(x,y,s);计算图像I(x,y)每一点在不同尺度下的灰度参数I(x,y,s);计算得到基于对称性的多尺度视网膜血管滤波增强函数F(x,y,s),使用该血管函数对视网膜图像进行滤波,从而达到视网膜血管滤波增强的效果。
图1为本发明基于对称性的视网膜血管滤波增强方法的流程图,如图1所示,所述方法包括以下步骤:
步骤S1:获得视网膜灰度图像I(x,y);
获得视网膜灰度图像的方法有很多,比如通过眼底照相机拍摄得到视网膜彩色图像,然后对于该视网膜彩色图像进行灰度化(比如取绿色通道)就可以得到视网膜灰度图像。视网膜灰度图像的获取不是本发明的重点,因此,本发明对其不作过多描述。
步骤S2:用二维多尺度高斯滤波器G(x,y,s)对于所述视网膜灰度图像I(x,y)进行滤波,得到高斯滤波图像VG(x,y,s),其中,x,y表示图像中某一点的坐标,s为对应的尺度;
其中,经过所述二维多尺度高斯滤波器G(x,y,s)滤波后,在x方向得到的滤波结果表示为在y方向得到的滤波结果表示为其中,*表示卷积操作,
在本发明一实施例中,所述二维多尺度高斯滤波器G(x,y,s)的大小为6*6,均方误差σ=3,尺度S从0.2到3。
步骤S3:基于得到的高斯滤波图像VG(x,y,s)计算所述视网膜灰度图像的梯度合成向量T(x,y,s);
其中,T(x,y,s)=gx(x,y,s)+gy(x,y,s),在本发明一实施例中,采用平行四边形法则对于x、y方向上得到的滤波向量进行相加操作,得到所述视网膜图像的高斯梯度和幅值信息,并且在相加之前对于所述x、y方向上得到的滤波向量分别进行正向和逆向旋转90°,然后再进行相加处理。
步骤S4:计算所述视网膜灰度图像I(x,y)中的每一点在不同尺度下的梯度方向参数D(x,y,s);
其中,所述梯度方向参数D(x,y,s)表示为:
D(x,y,s)=|cos(∠(gx(x,y,s),gy(x,y,s)))|,
其中,∠表示两个向量的角度。
考虑到尺度因子S的取值范围为从0.2到3,因此,在本发明一实施 例中,按照0.4的阶梯进行跳跃,得到所述视网膜灰度图像I(x,y)中的每一点在多个尺度下的梯度方向因子。
步骤S5:计算所述视网膜灰度图像I(x,y)中的每一点在不同尺度下的梯度幅值参数M(x,y,s);
其中,所述梯度幅值参数M(x,y,s)表示为:
其中,δm为梯度幅值标准差。
考虑到尺度因子S的取值范围为从0.2到3,因此,在本发明一实施例中,按照0.4的阶梯进行跳跃,计算得到所述视网膜灰度图像I(x,y)中的每一点在多个尺度下的梯度幅值因子。
步骤S6:计算所述视网膜灰度图像I(x,y)中的每一点在不同尺度下的灰度参数I(x,y,s);
其中,所述灰度参数I(x,y,s)表示为:
其中,δi为灰度标准差;i1(x,y,s),i2(x,y,s)为对称的两个像素点的灰度值。
考虑到尺度因子S的取值范围为从0.2到3,因此,在本发明一实施例中,按照0.4的阶梯进行跳跃,计算得到所述视网膜灰度图像I(x,y)中的每一点在多个尺度下的灰度因子。
步骤S7:根据所述梯度合成向量G(x,y,s)、梯度方向参数D(x,y,s)、梯度幅值参数M(x,y,s)、灰度参数I(x,y,s),计算得到基于对称性的多尺度视网膜血管滤波增强函数F(x,y,s),利用所述血管滤波增强函数对于所述视网膜灰度图像进行滤波,得到经过滤波和增强的视网膜图像。
其中,多尺度视网膜血管滤波增强函数F(x,y,s)表示为:
其中,s表示尺度,smin,smax分别表示尺度空间中的最小和最大尺度。 为验证本发明的有效性和实用性,接下来在两个国际公认的视网膜图像数据库(STARE和DRIVE)上进行了实验。这两个数据库都分别提供20幅视网膜图像用于算法测试,两个手动分割的血管网络数据集作为参考。选择每个数据库中的2组图片进行展出,如图2A、图3A、图4A、图5A所示。其中,图2A、图3A、图4A、图5A为眼底照相机拍摄的视网膜图像,图2B、图3B、图4B、图5B为经过本发明方法滤波和增强后得到的视网膜图像。由图2B、图3B、图4B、图5B所示的实验结果可知,本发明方法能够有效地对视网膜图像的血管结构进行滤波增强,同时有效的对非血管结构进行抑制。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (9)

1.一种基于对称性的视网膜血管滤波增强方法,其特征在于,该方法包括以下步骤:
步骤S1:获得视网膜灰度图像I(x,y);
步骤S2:用二维多尺度高斯滤波器G(x,y,s)对于所述视网膜灰度图像I(x,y)进行滤波,得到高斯滤波图像VG(x,y,s),其中,x,y表示图像中某一点的坐标,s为对应的尺度;
步骤S3:基于得到的高斯滤波图像VG(x,y,s)计算所述视网膜灰度图像的梯度合成向量T(x,y,s);
步骤S4:计算所述视网膜灰度图像I(x,y)中的每一点在不同尺度下的梯度方向参数D(x,y,s);
步骤S5:计算所述视网膜灰度图像I(x,y)中的每一点在不同尺度下的梯度幅值参数M(x,y,s);
步骤S6:计算所述视网膜灰度图像I(x,y)中的每一点在不同尺度下的灰度参数I(x,y,s);
步骤S7:根据所述梯度合成向量G(x,y,s)、梯度方向参数D(x,y,s)、梯度幅值参数M(x,y,s)、灰度参数I(x,y,s),计算得到基于对称性的多尺度视网膜血管滤波增强函数F(x,y,s),利用所述血管滤波增强函数对于所述视网膜灰度图像进行滤波,得到经过滤波和增强的视网膜图像;
其中,多尺度视网膜血管滤波增强函数F(x,y,s)表示为:
F ( x , y , s ) = Σ s = s m i n s m a x | G ( x , y , s ) | · D ( x , y , s ) · M ( x , y , s ) · I ( x , y , s ) ,
其中,s表示尺度,smin,smax分别表示尺度空间中的最大和最小尺度。
2.根据权利要求1所述的方法,其特征在于,用二维多尺度高斯滤波器G(x,y,s)对于所述视网膜灰度图像I(x,y)进行滤波后,x方向上得到的滤波结果为y方向得到的滤波结果为其中,*表示卷积操作。
3.根据权利要求1所述的方法,其特征在于,所述二维多尺度高斯滤波器G(x,y,s)的大小为6*6,均方误差σ=3,尺度S从0.2到3。
4.根据权利要求1所述的方法,其特征在于,所述步骤S3中,T(x,y,s)=gx(x,y,s)+gy(x,y,s)。
5.根据权利要求1所述的方法,其特征在于,所述步骤S3中,采用平行四边形法则对于x、y方向上得到的滤波向量进行相加操作。
6.根据权利要求1或5所述的方法,其特征在于,所述步骤S3中,首先对于x、y方向上得到的滤波向量分别正向和逆向旋转90°,然后再进行相加处理。
7.根据权利要求1所述的方法,其特征在于,所述步骤S4中,所述梯度方向参数D(x,y,s)表示为:
D(x,y,s)=|cos(∠(gx(x,y,s),gy(x,y,s)))|,
其中,∠表示两个向量的角度。
8.根据权利要求1所述的方法,其特征在于,所述步骤S5中,所述梯度幅值参数M(x,y,s)表示为:
M ( x , y , s ) = exp ( - ( | | g x ( x , y , s ) | | - | | g y ( x , y , s ) | | ) 2 2 δ m 2 ) ,
其中,δm为梯度幅值标准差。
9.根据权利要求1所述的方法,其特征在于,所述步骤S6中,所述灰度参数I(x,y,s)表示为:
I ( x , y , s ) = exp ( - ( | | i 1 ( x , y , s ) | | - | | i 2 ( x , y , s ) | | ) 2 2 δ i 2 )
其中,δi为灰度标准差;i1(x,y,s),i2(x,y,s)为对称的两个像素点的灰度值。
CN201410315974.9A 2014-07-03 2014-07-03 基于对称性的视网膜血管滤波增强方法 Active CN104077754B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410315974.9A CN104077754B (zh) 2014-07-03 2014-07-03 基于对称性的视网膜血管滤波增强方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410315974.9A CN104077754B (zh) 2014-07-03 2014-07-03 基于对称性的视网膜血管滤波增强方法

Publications (2)

Publication Number Publication Date
CN104077754A CN104077754A (zh) 2014-10-01
CN104077754B true CN104077754B (zh) 2017-03-01

Family

ID=51598996

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410315974.9A Active CN104077754B (zh) 2014-07-03 2014-07-03 基于对称性的视网膜血管滤波增强方法

Country Status (1)

Country Link
CN (1) CN104077754B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104899862A (zh) * 2015-04-01 2015-09-09 武汉工程大学 基于全局与局部阈值的视网膜血管分割算法
CN104732499A (zh) * 2015-04-01 2015-06-24 武汉工程大学 基于多尺度多方向的视网膜图像增强算法
CN105654439B (zh) * 2015-12-31 2018-09-11 中国科学院自动化研究所 基于梯度调整的血管图像增强方法
CN106558031B (zh) * 2016-12-02 2018-12-28 北京理工大学 一种基于成像模型的彩色眼底图的图像增强方法
CN118014881B (zh) * 2024-04-08 2024-06-28 山东拓新电气有限公司 基于计算机视觉的盾构机导向系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
《A New Approach to Automated Retinal Vessel Segmentation Using Multiscale Analysis》;Qin Li等;《18th International Conference on Pattern Recognition (ICPR"06)》;20061231(第4期);77-80 *
《Retinal Image Graph-Cut Segmentation Algorithm Using Multiscale Hessian-Enhancement-Based Nonlocal Mean Filter》;Jian Zheng等;《Computational & Mathematical Methods in Medicine》;20130411(第1期);84-104 *
《视网膜血管图像增强与检测》;黎飞;《中国优秀硕士学位论文全文数据库 信息科技辑》;20130215(第02期);I138-1614 *
《视网膜血管增强方法研究》;廖苗;《中国优秀硕士学位论文全文数据库 信息科技辑》;20140515(第05期);I138-2893 *

Also Published As

Publication number Publication date
CN104077754A (zh) 2014-10-01

Similar Documents

Publication Publication Date Title
JP2022065067A (ja) 灌流デジタルサブトラクション血管造影
US10769477B2 (en) Method, apparatus, device and storage medium for extracting a cardiovisceral vessel from a CTA image
CN104077754B (zh) 基于对称性的视网膜血管滤波增强方法
CN107292835B (zh) 一种眼底图像视网膜血管自动矢量化的方法及装置
CN105160660A (zh) 基于多特征高斯拟合的活动轮廓血管提取方法及系统
Su et al. A new method for linear feature and junction enhancement in 2D images based on morphological operation, oriented anisotropic Gaussian function and Hessian information
CN109671049B (zh) 一种医学图像处理方法、系统、设备、存储介质
WO2018049598A1 (zh) 一种眼底图像增强方法及系统
CN109754388B (zh) 一种颈动脉狭窄程度计算方法、装置及存储介质
CN102393956B (zh) 基于径向对称性变换的视网膜图像滤波方法
CN106504199B (zh) 一种眼底图像增强方法及系统
CN103295195B (zh) 低对比度图像的脉管增强的方法及其系统
Karasulu Automatic extraction of retinal blood vessels: a software implementation
JP6603814B2 (ja) 網膜画像処理
Aruchamy et al. Automated glaucoma screening in retinal fundus images
Ben Gharsallah et al. Phase congruency-based filtering approach combined with a convolutional network for lung CT image analysis
Samagaio et al. Optical coherence tomography denoising by means of a fourier butterworth Filter-Based approach
Abdallah et al. Restoration of retinal images using anisotropic diffusion like algorithms
Sabaz et al. ROI detection and vessel segmentation in retinal image
Abirami et al. Automated brain tumor detection and identification using medical imaging
CN103020910B (zh) 基于各向异性和各向同性的视网膜图像特征点增强方法
CN113538419A (zh) 一种图像处理方法和系统
Tache Vessels enhancement in x-ray angiograms
CN112017132A (zh) 基于最大曲率法和多尺度Hessian矩阵的静脉图像增强方法
Huang et al. Non-local Hessian-based weighted filter for 2D noisy angiogram image enhancement with cardiovascular tree structure preservation

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant