CN113506284A - 一种眼底图像微血管瘤检测装置、方法及存储介质 - Google Patents
一种眼底图像微血管瘤检测装置、方法及存储介质 Download PDFInfo
- Publication number
- CN113506284A CN113506284A CN202110847212.3A CN202110847212A CN113506284A CN 113506284 A CN113506284 A CN 113506284A CN 202110847212 A CN202110847212 A CN 202110847212A CN 113506284 A CN113506284 A CN 113506284A
- Authority
- CN
- China
- Prior art keywords
- image
- microangioma
- candidate
- area
- candidate region
- 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
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000001514 detection method Methods 0.000 title claims description 28
- 210000001525 retina Anatomy 0.000 claims abstract description 7
- 239000013598 vector Substances 0.000 claims description 35
- 238000013459 approach Methods 0.000 claims description 9
- 230000010339 dilation Effects 0.000 claims description 8
- 238000012549 training Methods 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 7
- 238000012216 screening Methods 0.000 claims description 7
- 230000011218 segmentation Effects 0.000 claims description 6
- 239000003550 marker Substances 0.000 claims description 5
- 230000002207 retinal effect Effects 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 4
- 206010012601 diabetes mellitus Diseases 0.000 abstract description 2
- 238000012545 processing Methods 0.000 abstract description 2
- 210000004204 blood vessel Anatomy 0.000 description 9
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 4
- 238000013135 deep learning Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 206010012689 Diabetic retinopathy Diseases 0.000 description 2
- 208000017442 Retinal disease Diseases 0.000 description 2
- 206010038923 Retinopathy Diseases 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 208000035269 cancer or benign tumor Diseases 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 238000013136 deep learning model Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 201000011066 hemangioma Diseases 0.000 description 1
- 231100000915 pathological change Toxicity 0.000 description 1
- 230000036285 pathological change Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000004256 retinal image Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
- G06T5/30—Erosion or dilatation, e.g. thinning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/187—Segmentation; Edge detection involving region growing; involving region merging; involving connected component labelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/194—Segmentation; Edge detection involving foreground-background segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Quality & Reliability (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Image Analysis (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本发明公开了一种眼底图像微血管瘤检测装置、方法及存储介质,涉及医学图像处理及机器视觉应用领域。本发明首先对输入的彩色视网膜图像提取待检测图像,再通过小目标去除、测地膨胀等操作得到微血管瘤候选区模板和正负样本图片;针对微血管瘤候选区提取能量特征,同时针对微血管瘤形态手工设计特征,最后进行将传统特征和手工设计的特征进行级联;将特征和标签送入已训练好的分类器进行分类,从而从糖网图像中检测出微血管瘤所在位置。本方法能够通过糖尿病眼底图像检测出微小目标微血管瘤,具有较高的准确率,能够辅助眼底医生更方便的观察到微血管瘤的存在。
Description
技术领域
本发明涉及眼底图像微血管瘤检测的一种眼底图像微血管瘤检测装置、方法 及存储介质。
背景技术
糖尿病患者晚期会引发视网膜发生病变,眼底图像微血管瘤 (Microaneurysm,MA)是糖尿病视网膜病变(糖网)的初期症状,因此,实现 眼底图像的微血管瘤检测并及时治疗有助于防止视网膜病变进一步加深。人工 对眼底图像的微血管瘤检测主要依赖于眼科医生直接对视网膜图像进行观测, 但由于视网膜结构复杂、微血管瘤面积微小且局部对比度较低,人眼观测微血 管瘤费时费力,任务量巨大,且偏远地区缺乏有经验的眼科医生,通过计算机 视觉技术实现视网膜微血管瘤的自动检测,有助于缓解眼科医生压力,同时有助于医疗资源下沉,具有深远的医学意义。
目前眼底图像的微血管瘤检测方法主要包括基于深度学习的方法和基于分类 器的方法。基于深度学习的方法主要是采用深度学习模型构建端到端的卷积神 经网络,如像素级目标分割的语义分割网络、标记目标所在区域的目标检测网 络。李英采用了SSD目标检测网络实现了微血管瘤的检测,但准确率不高,同 时深度学习网络框架由于参数量大且效果不稳定,难以集成上软件进行实际使 用。基于分类器的方法首先会提取MA候选区,再对候选区进行特征建模与分 类。Orlando等首先采用背景估计法提取MA候选区,然后对候选区提取纹理特 征、灰度特征、体特征和深度特征,然后送入分类器进行分类;Dasht采用LCF 滤波器提取微血管瘤候选区,并将滤波响应值与传统特征融合,共同作为训练 分类器所需要的特征。这些方法并未针对MA本身特性进行分析,因此依然存在 准确率不高、鲁棒性不强的特点。因此,目前基于计算机图像处理技术的微血 管瘤检测方法还存在着鲁棒性不高、检测准确率不高、难以集成使用等问题。
发明内容
本发明所要解决的技术问题就是提供一种眼底图像微血管瘤的检测方法,能 够较为精确的检测出眼底图像上的微血管瘤,排除血管、背景噪声等干扰结 构,有利于医生发现微血管瘤所在位置,进而给出诊断和治疗,防止患者病变 加深。
为了解决上述技术问题,达到上述目的,本发明所采用的技术方案如下所 述。
一种眼底图像微血管瘤检测方法,包括如下步骤:
步骤1:输入眼底图像,提取包含微血管瘤信息的待检测图像,对待检测图 像进行小目标去除,得到模糊糖网图像,再对待处理图像与模糊糖网图像反复 进行测地膨胀,得到糖网背景图像,转入步骤2;
步骤2:将糖网背景图像与待处理图像相减,并通过归一化与特定灰度分 割,得到微血管瘤候选区模板图,转入步骤3;
步骤3:对微血管瘤候选区模板图进行连通域分析,计算各连通域的面积以 及中心坐标,以各中心坐标为图像中心从输入眼底图像的绿色通道提取一定大 小的图像片,并使用对应的连通域面积进行筛选,去除连通域面积较少和较大 对应的图像片,得到微血管瘤候选区图像,转入步骤4;
步骤4:利用步骤3中的微血管瘤候选区图像,设计手工特征提取器,提取 手工特征,得到最终的特征向量,转入步骤5;
步骤5:将步骤4中的各候选区的特征向量和对应的类别标签送入分类器进 行训练,并利用训练好的模型对测试时的微血管瘤候选区域特征向量进行分 类,判断每一个候选区域的类别,最终输出微血管瘤在眼底图像上的中心坐 标。
上述技术方案中,步骤1中具体有以下几个步骤:
步骤1.1:从输入彩色眼底图像提取绿色通道图像,并对其进行反射得到待 检测图像I;
步骤1.2:对待检测图像I采用滤波器进行小目标去除,得到模糊糖网图像 Ivague;
步骤1.3:将模糊糖网图像Ivague作为图像L,将待检测图像I作为图像T,通 过式(1)经过反复测地膨胀后得到视网膜背景图像Ibackground,式(1)如下:
上述,B其中表示大小为3×3值为1的结构元,表示采用结构元B对L 的膨胀操作,∩表示两图像空间相应元素中最小灰度形成的阵列,表示 标记图像L关于模板图像T的一次测地膨胀操作,整个式子迭代运算,将一次测 地膨胀操作的结果作为下次测地膨胀标记图像,并循环往复直到结果不再发生 变换。
上述技术方案中,所述步骤2具体有以下几个步骤:
步骤2.1:将步骤1中的待检测图像I减去步骤1中的视网膜背景图像 Ibackground得到Idif,并对Idif进行归一化处理得到Inormal;
步骤2.2:设定阈值t1,对Inomal进行分割,像素大于t1则置1,否则置0, 最终得到微血管瘤候选区模板图Icandidate。
上述技术方案中,所述步骤3具体有以下几个步骤:
步骤3.1:对步骤2得到的微血管瘤候选区模板图Icandidate进行连通域分 析,计算各连通域的面积以及对应的中心坐标,筛选出连通的域面积大于Smin小 于Smax的中心坐标集合centers=c1,c2,...,cn,其中ci表示第i个连通域的中心坐 标,i∈{1,2,3,...,n},n表示微血管瘤候选区个数;
步骤3.2:通过步骤3.1得到的中心坐标集合centers,以每个坐标作为图 像片中心,从步骤1中的待检测图像I提取大小k×k的图像片,构成微血管瘤候 选区图像Ipatches=p1,p2,...,pn,其中pi表示第i个候选区图像,i∈{1,2,3,...,n}。
上述技术方案中,所述步骤4具体有以下几个步骤:
步骤4.1:对步骤3.2得到的微血管瘤候选区图像提取用于描述灰度信息的 能量特征,主要包括灰度平均值、方差、偏度、对比度、熵等;将能量特征定 义为attrib1:
步骤4.2:针对微血管瘤图像一定程度上具有旋转不变性,将候选区图像pi顺时针旋转90°得到旋转后的候选区图像将pi和通过相同顺序平铺为k2维向量分别得到vi和通过式(2)的结果衡量其旋转不变性,并将该特征定 义为attrib2;式(2)如下:
其中,vi={vi1,vi2,vi3,...,vik 2},vij表示vi中的第j个元素,表示中的第j个元素,k2表示向量维度,数值与单张候选区图像的元素个数相等;例如,类似 微血管瘤形态的候选区可以简单描述为经过旋转后依然为 两者通过相同顺序平铺得到[0,1,0,1,0,1,0,1,0]。类似血管形态的候选 区可以描述为经过旋转后得到两者通过相同顺序得到 [2,0,0,0,2,0,0,0,2]和[0,0,2,0,2,0,2,0,0]。采用式(2)对其计算,若候选区为微血 管瘤则计算结果接近1,若为血管则会远离1;
步骤4.3:针对候选区图像上,微血管瘤区域像素会集中在候选区图像的中 央,且灰度值低于背景区域,而背景噪声较低灰度值形成的区域会随机分布在 候选区图像中的各个区域,设定阈值t2,对候选区图像pi进行分割,像素大于t1则置0,否则置1,得到低灰度像素区域li。
步骤4.4:对步骤4.3得到的低灰度像素区域li进行连通域分析,得到连通 域个数m,并计算得到各连通域像素面积Ai=Ai1,Ai2,Ai3,...,Aim,并通过式(3) 计算各连通域像素面积占总面积之比Pi=Pi1,Pi2,Pi3,...,Pim,式(3)如下式所 示:
其中,Aij表示Ai中的第j个连通域的面积,Pij其中表示Pi中的第j个连通域面 积占总体连通域面积之比;
步骤4.5:通过式(4)计算候选区pi对应的低灰度像素区域的混乱程度Hi, 并将其作为attrib3,式(4)如下所示:
其中,主要对混乱程度进行归一;用于描述混乱程 度,若m为1,表示低灰度区域只有一个连通区域,此时 则表示当前候选区低灰度区域图像单一;当m大于1,Ai中若存在一个连通域远大于其他连通域面积,混乱程度依然趋近于0;若存在多 个面积相差不大的连通域,则的值趋近于log2 m,经过归一后 该值趋近于1;
步骤4.6:将步骤4.5得到的attrib3和步骤4.2得到的attrib2依次接在步 骤4.1得到的attrib1的后面,形成最终的特征向量。
上述技术方案中,所述步骤5具体有以下几个步骤:
步骤5.1:输入多张眼底图像通过步骤1,2,3,4得到大量微血管瘤候选区 的最终的特征向量,并相应打上标签,若该候选区为微血管瘤则标记为1,若 不为微血管瘤则标记为0,将特征向量和标签一同送入分类器进行训练,得到 训练好的分类器,转入步骤5.2;
步骤5.2:输入需要检测的彩色眼底图像,通过步骤1,2,3,4得到该图像 的微血管瘤候选区图像以及对应的最终特征向量,利用步骤5.1训练好的分类 器对其预测,得到各候选区图像的类别,转入步骤5.3。
步骤5.3:对于分类为微血管瘤的候选区图像,同时记录中心坐标集合 centers中对应的坐标,可在输入彩色眼底图像对应的位置依次标出,最终实现 微血管瘤的检测。
本发明还提供了一种眼底图像微血管瘤检测装置,其特征在于,包括如下 步骤:
糖网背景图像模块:输入眼底图像,提取包含微血管瘤信息的待检测图像, 对待检测图像进行小目标去除,得到模糊糖网图像,再对待处理图像与模糊糖 网图像反复进行测地膨胀,得到糖网背景图像;
微血管瘤候选区模板图模块:将糖网背景图像与待处理图像相减,并通过归 一化与特定灰度分割,得到微血管瘤候选区模板图;
微血管瘤候选区图像模块:对微血管瘤候选区模板图进行连通域分析,计算 各连通域的面积以及中心坐标,以各中心坐标为图像中心从输入眼底图像的绿 色通道提取一定大小的图像片,并使用对应的连通域面积进行筛选,去除连通 域面积较少和较大对应的图像片,得到微血管瘤候选区图像;
最终的特征向量模块:利用微血管瘤候选区图像,设计手工特征提取器,提 取手工特征,得到最终的特征向量;
结果输出模块:将各候选区的特征向量和对应的类别标签送入分类器进行训 练,并利用训练好的模型对测试时的微血管瘤候选区域特征向量进行分类,判 断每一个候选区域的类别,最终输出微血管瘤在眼底图像上的中心坐标。
本发明还提供了一种存储介质,所述存储介质上存储有一种眼底图像微血管 瘤检测程序,所述一种眼底图像微血管瘤检测程序被处理器执行时实现一种眼 底图像微血管瘤检测方法的步骤。
因为本发明采用上述技术手段,因此具备以下有益效果:
在本文算法中,针对微血管瘤的微小特性和视网膜背景结构设计了步骤 1,2,3,能够实现微血管瘤的候选区提取,在候选区中即包含微血管瘤,同时也 包含其他结构,相比较其他的候选区提取技术,该技术能够将一些体型较大结 构如软性渗出、神经视盘中部分高频结构、大面积的血管等排除在外,对于形 态类似的硬性渗出,也能排除,最终得到的候选区类别少、结构不复杂,有利 于后续的特征建模与分类;并通过步骤4进行进一步特征建模,本人对大量眼 底图像通过步骤1,2,3得到大量的微血管瘤候选区,通过详细的观察与统计, 我们发现在微血管瘤候选区中,除了正样本微血管瘤以外,主要存在两类负样本,分别为血管和背景噪声,传统算法一般直接采用常规特征对其进行分类。 本文首先采用步骤4.1能量特征来描述微血管瘤的灰度特性,虽然这是常规 的,但也是必要的;然后针对微血管瘤与血管的形态手动设计高区分度的特 征,针对微血管瘤为圆形或椭圆形,旋转之后结构变换不大,而血管旋转之后 除了中心部位以外其他区域的像素几乎发生了翻转,进而根据步骤4.2设计了 旋转不变性特征来区分微血管瘤与血管;然后通过大量实验观察到微血管瘤与 背景噪声在低灰度像素区域的分布显然不同,微血管瘤低灰度像素区域分布单 一,主要集中在中央,而背景噪声在低灰度像素区域的分布杂乱,本文首次提 出了描述图像某灰度区域的结构混乱性特征,如步骤4.3、4.4、4.5所示,该 特征不仅在低灰度区域区分微血管瘤与背景噪声具有巨大的作用,同时能应用 到其他分类场景以及其他灰度区域;最终将所有特征进行简单级联用于模型训 练与分类。常规的分类手段一般会采用灰度特征提取器、纹理特征提取器、体 特征提取器等进行特征提取与融合,该方式存在大量不必要的特征且由于特征 过多造成检测时间常,而本文通过分析正负样本并手动设计了高相关特征,使 得对目标的特征描述更为具体,因此最终模型的检测准确率更高。
附图说明
图1为一种眼底图像微血管瘤检测方法设计流程;
图2为输入眼底图像和待检测图像,其中(a)为眼底图像,(b)为待检测图 像;
图3为候选区域提取示意图,其中(a)为待检测图像,(b)为经过小目标去除后 的图像,(c)为视网膜背景图像,(d)为视网膜背景图像与待检测图之差, (e)为候选区模板图像。
图4为区别微血管瘤与血管所设定的旋转不变性特征示意图;上图为微血管 瘤,下图为血管,可以看到经过旋转后微血管瘤依然能保持一定的不变性,可 采用步骤4.2中的特征进行区分。
图5为17×17的候选区图像;其中(a)为正样本,即微血管瘤,(b)为负样 本;
图6为区别微血管瘤与背景噪声所设定的混乱性特征示意图;(a)为微血 管瘤低灰度区域,(b)为背景噪声低灰度区域,可以采用步骤4.3中的特征进 行区分。
图7为微血管瘤检测标记图。
具体实施方式
下面结合试验例及具体实施方式对本发明作进一步的详细描述。但不应将此 理解为本发明上述主题的范围仅限于以下的实施例,凡基于本发明内容所实现 的技术均属于本发明的范围。
本发明提出了一种眼底图像微血管瘤检测方法,能够检测出眼底图像的微血 管瘤区域,具有较高的特异性和灵敏度,整个算法设计方案流程如图1所示, 包括步骤:
上述技术方案中,所述步骤1中具体有以下几个步骤:
步骤1.1:从输入彩色眼底图像提取绿色通道图像,并对其进行反射得到待 检测图像I。本示例中,输入彩色眼底图像的大小为2544×1696×3。
步骤1.2:对待检测图像I采用滤波器进行小目标去除,得到模糊糖网图像 Ivague;本示例中,采用了15×15的中值滤波器去除了眼底图像小目标。
步骤1.3:将Ivague作为图像L,将I作为图像T,通过式(1)经过迭代测地 膨胀后得到视网膜背景图像Ibackground。式(1)如下:
上述,B其中表示大小为3×3值为1的结构元,表示采用结构元B对L 的膨胀操作,∩表示两图像空间相应元素中最小灰度形成的阵列。表示 标记图像L关于模板图像T的一次测地膨胀操作。整个式子迭代运算,将一次测 地膨胀操作的结果作为下次测地膨胀标记图像,并循环往复直到结果不再发生 变换。
上述技术方案中,所述步骤2具体有以下几个步骤:
步骤2.1:将步骤1中的待检测图像I减去步骤1中的视网膜背景图像 Ibackground得到Idif,并对Idif进行归一化处理得到Inormal;
步骤2.2:设定阈值t1,对Inormal进行分割,像素大于t1则置1,否则置0。 最终得到微血管瘤候选区模板图Icandidate。本示例中,阈值t1的值为0.6。
上述技术方案中,所述步骤3具体有以下几个步骤:
步骤3.1:对步骤2得到的微血管瘤候选区模板图Icandidate进行连通域分 析。计算各连通域的面积以及对应的中心坐标,筛选出连通域的面积大于Smin小 于Smax的中心坐标集合centers=c1,c2,...,cn,其中ci表示第i个连通域的中心坐 标,i∈{1,2,3,...,n},n表示微血管瘤候选区个数。本示例中,Smin=1,Smax=100。
步骤3.2:通过步骤3.1得到的中心坐标集合centers,以每个坐标作为图 像片中心,从步骤1中的待检测图像I提取大小k×k的图像片,构成微血管瘤候 选区图像Ipatches=p1,p2,...,pn,其中pi表示第i个候选区图像。本示例中,k=17 ,即每个候选区图像的大小为17×17。
上述技术方案中,所述步骤4具体有以下几个步骤:
步骤4.1:对步骤3.2得到的微血管瘤候选区图像提取用于描述灰度信息的 能量特征,主要包括灰度平均值、方差、偏度、对比度、熵等;将能量特征定 义为attrib1;
步骤4.2:针对微血管瘤图像一定程度上具有旋转不变性,将候选区图像pi顺时针旋转90°得到将和通过相同顺序平铺为k2维向量分别得到vi和通过式(2)的结果衡量其旋转不变性,并将该特征定义为attrib2;式(2)如
步骤4.3:针对候选区图像上,微血管瘤区域像素会集中在候选区图像的中 央,且灰度值低于背景区域,而背景噪声较低灰度值形成的区域会随机分布在 候选区图像中的各个区域。设定阈值t2,对候选区图像pi进行分割,像素大于t1则置0,否则置1,得到低灰度像素区域li。本示例中,t2=87。
步骤4.4:对步骤4.3得到的低灰度像素区域li进行连通域分析,得到连通 域个数m,并计算得到各连通域像素面积Ai=Ai1,Ai2,Ai3,...,Aim,并通过式(3) 计算各连通域像素面积占总面积之比Pi=Pi1,Pi2,Pi3,...,Pim,式(3)如下式所 示:
步骤4.5:通过式(4)计算候选区pi对应的低灰度像素区域的混乱程度Hi, 并将其作为attrib3,式(4)如下所示:
其中,主要对混乱程度进行归一;用于描述混乱程 度,若m为1,表示低灰度区域只有一个连通区域,此时 则表示当前候选区低灰度区域图像单一;当m大于1,Ai中若存在一个连通域远大于其他连通域面积,混乱程度依然趋近于0;若存在多 个面积相差不大的连通域,则的值趋近于log2 m,经过归一后 该值趋近于1。
步骤4.6:将步骤4.5得到的attrib3和步骤4.2得到的attrib2依次接在步 骤4.1得到的attrib1的后面,形成最终的特征向量。
上述技术方案中,所述步骤5具体有以下几个步骤:
步骤5.1:输入多张眼底图像通过步骤1,2,3,4得到大量微血管瘤候选区 的最终的特征向量,并相应打上标签,若该候选区为微血管瘤则标记为1,若 不为微血管瘤则标记为0,将特征向量和标签一同送入分类器进行训练,得到 训练好的分类器,转入步骤5.2;本示例中,采用lightgbm框架作为模型训练 的分类器,采用gbdt作为拟合方法,针对从4112个微血管瘤候选区提取出来 的特征进行五折交叉验证训练,经过500次迭代后得到模型。
步骤5.2:输入需要检测的彩色眼底图像,通过步骤1,2,3,4得到该图像 的微血管瘤候选区图像以及对应的最终特征向量,利用步骤5.1训练好的分类 器对其预测,得到各候选区图像的类别,转入步骤5.3。
步骤5.3:对于分类为微血管瘤的候选区图像,同时记录中心坐标集合 centers中对应的坐标,可在输入彩色眼底图像对应的位置依次标出,最终实现 微血管瘤的检测。
Claims (8)
1.一种眼底图像微血管瘤检测方法,其特征在于,包括如下步骤:
步骤1:输入眼底图像,提取包含微血管瘤信息的待检测图像,对待检测图像进行小目标去除,得到模糊糖网图像,再对待处理图像与模糊糖网图像反复进行测地膨胀,得到糖网背景图像,转入步骤2;
步骤2:将糖网背景图像与待处理图像相减,并通过归一化与特定灰度分割,得到微血管瘤候选区模板图,转入步骤3;
步骤3:对微血管瘤候选区模板图进行连通域分析,计算各连通域的面积以及中心坐标,以各中心坐标为图像中心从输入眼底图像的绿色通道提取一定大小的图像片,并使用对应的连通域面积进行筛选,去除连通域面积较少和较大对应的图像片,得到微血管瘤候选区图像,转入步骤4;
步骤4:利用步骤3中的微血管瘤候选区图像,设计手工特征提取器,提取手工特征,得到最终的特征向量,转入步骤5;
步骤5:将步骤4中的各候选区的特征向量和对应的类别标签送入分类器进行训练,并利用训练好的模型对测试时的微血管瘤候选区域特征向量进行分类,判断每一个候选区域的类别,最终输出微血管瘤在眼底图像上的中心坐标。
2.根据权利要求1所述的一种眼底图像微血管瘤检测方法,其特征在于,步骤1中具体有以下几个步骤:
步骤1.1:从输入彩色眼底图像提取绿色通道图像,并对其进行反射得到待检测图像I;
步骤1.2:对待检测图像I采用滤波器进行小目标去除,得到模糊糖网图像Ivague;
步骤1.3:将模糊糖网图像Ivague作为图像L,将待检测图像I作为图像T,通过式(1)经过反复测地膨胀后得到视网膜背景图像Ibackground,式(1)如下:
3.根据权利要求1所述的一种眼底图像微血管瘤检测方法,其特征在于,所述步骤2具体有以下几个步骤:
步骤2.1:将步骤1中的待检测图像I减去步骤1中的视网膜背景图像Ibackground得到Idif,并对Idif进行归一化处理得到Inormal;
步骤2.2:设定阈值t1,对Inormal进行分割,像素大于t1则置1,否则置0,最终得到微血管瘤候选区模板图Icandidate。
4.根据权利要求1所述的一种眼底图像微血管瘤检测方法,其特征在于,所述步骤3具体有以下几个步骤:
步骤3.1:对步骤2得到的微血管瘤候选区模板图Icandidate进行连通域分析,计算各连通域的面积以及对应的中心坐标,筛选出连通的域面积大于Smin小于Smax的中心坐标集合centers=c1,c2,...,cn,其中ci表示第i个连通域的中心坐标,i∈{1,2,3,...,n},n表示微血管瘤候选区个数;
步骤3.2:通过步骤3.1得到的中心坐标集合centers,以每个坐标作为图像片中心,从步骤1中的待检测图像I提取大小k×k的图像片,构成微血管瘤候选区图像Ipatches=p1,p2,...,pn,其中pi表示第i个候选区图像,i∈{1,2,3,...,n}。
5.根据权利要求1所述的一种眼底图像微血管瘤检测方法,其特征在于,所述步骤4具体有以下几个步骤:
步骤4.1:对步骤3.2得到的微血管瘤候选区图像提取用于描述灰度信息的能量特征,主要包括灰度平均值、方差、偏度、对比度、熵等;将能量特征定义为attrib1;
步骤4.2:针对微血管瘤图像一定程度上具有旋转不变性,将候选区图像pi顺时针旋转90°得到旋转后的候选区图像将pi和通过相同顺序平铺为k2维向量分别得到vi和通过式(2)的结果衡量其旋转不变性,并将该特征定义为attrib2;式(2)如下:
步骤4.3:针对候选区图像上,微血管瘤区域像素会集中在候选区图像的中央,且灰度值低于背景区域,而背景噪声较低灰度值形成的区域会随机分布在候选区图像中的各个区域,设定阈值t2,对候选区图像pi进行分割,像素大于t1则置0,否则置1,得到低灰度像素区域li;
步骤4.4:对步骤4.3得到的低灰度像素区域li进行连通域分析,得到连通域个数m,并计算得到各连通域像素面积Ai=Ai1,Ai2,Ai3,...,Aim,并通过式(3)计算各连通域像素面积占总面积之比Pi=Pi1,Pi2,Pi3,...,Pim,式(3)如下式所示:
其中,Aij表示Ai中的第j个连通域的面积,Pij其中表示Pi中的第j个连通域面积占总体连通域面积之比;
步骤4.5:通过式(4)计算候选区pi对应的低灰度像素区域的混乱程度Hi,并将其作为attrib3,式(4)如下所示:
其中,主要对混乱程度进行归一;用于描述混乱程度,若m为1,表示低灰度区域只有一个连通区域,此时Pij=1, 则表示当前候选区低灰度区域图像单一;当m大于1,Ai中若存在一个连通域远大于其他连通域面积,混乱程度依然趋近于0;若存在多个面积相差不大的连通域,则的值趋近于log2m,经过归一后该值趋近于1;
步骤4.6:将步骤4.5得到的attrib3和步骤4.2得到的attrib2依次接在步骤4.1得到的attrib1的后面,形成最终的特征向量。
6.根据权利要求1所述的一种眼底图像微血管瘤检测方法,其特征在于,所述步骤5具体有以下几个步骤:
步骤5.1:输入多张眼底图像通过步骤1-4得到大量微血管瘤候选区的最终的特征向量,并相应打上标签,若该候选区为微血管瘤则标记为1,若不为微血管瘤则标记为0,将特征向量和标签一同送入分类器进行训练,得到训练好的分类器,转入步骤5.2;
步骤5.2:输入需要检测的彩色眼底图像,通过步骤1-4得到该图像的微血管瘤候选区图像以及对应的最终特征向量,利用步骤5.1训练好的分类器对其预测,得到各候选区图像的类别,转入步骤5.3;
步骤5.3:对于分类为微血管瘤的候选区图像,同时记录中心坐标集合centers中对应的坐标,可在输入彩色眼底图像对应的位置依次标出,最终实现微血管瘤的检测。
7.一种眼底图像微血管瘤检测装置,其特征在于,包括如下步骤:
糖网背景图像模块:输入眼底图像,提取包含微血管瘤信息的待检测图像,对待检测图像进行小目标去除,得到模糊糖网图像,再对待处理图像与模糊糖网图像反复进行测地膨胀,得到糖网背景图像;
微血管瘤候选区模板图模块:将糖网背景图像与待处理图像相减,并通过归一化与特定灰度分割,得到微血管瘤候选区模板图;
微血管瘤候选区图像模块:对微血管瘤候选区模板图进行连通域分析,计算各连通域的面积以及中心坐标,以各中心坐标为图像中心从输入眼底图像的绿色通道提取一定大小的图像片,并使用对应的连通域面积进行筛选,去除连通域面积较少和较大对应的图像片,得到微血管瘤候选区图像;
最终的特征向量模块:利用微血管瘤候选区图像,设计手工特征提取器,提取手工特征,得到最终的特征向量;
结果输出模块:将各候选区的特征向量和对应的类别标签送入分类器进行训练,并利用训练好的模型对测试时的微血管瘤候选区域特征向量进行分类,判断每一个候选区域的类别,最终输出微血管瘤在眼底图像上的中心坐标。
8.一种存储介质,其特征在于,所述存储介质上存储有一种眼底图像微血管瘤检测程序,所述一种眼底图像微血管瘤检测程序被处理器执行时实现如权利要求1-6所述的一种眼底图像微血管瘤检测方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110847212.3A CN113506284B (zh) | 2021-07-26 | 2021-07-26 | 一种眼底图像微血管瘤检测装置、方法及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110847212.3A CN113506284B (zh) | 2021-07-26 | 2021-07-26 | 一种眼底图像微血管瘤检测装置、方法及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113506284A true CN113506284A (zh) | 2021-10-15 |
CN113506284B CN113506284B (zh) | 2023-05-09 |
Family
ID=78014031
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110847212.3A Active CN113506284B (zh) | 2021-07-26 | 2021-07-26 | 一种眼底图像微血管瘤检测装置、方法及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113506284B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114529724A (zh) * | 2022-02-15 | 2022-05-24 | 推想医疗科技股份有限公司 | 图像目标的识别方法、装置、电子设备及存储介质 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105069803A (zh) * | 2015-08-19 | 2015-11-18 | 西安交通大学 | 一种基于彩色图像的糖尿病病变体微血管瘤的分类装置 |
CN107590941A (zh) * | 2017-09-19 | 2018-01-16 | 重庆英卡电子有限公司 | 拍照式混合火焰探测器及其探测方法 |
CN109977930A (zh) * | 2019-04-29 | 2019-07-05 | 中国电子信息产业集团有限公司第六研究所 | 疲劳驾驶检测方法及装置 |
CN110276356A (zh) * | 2019-06-18 | 2019-09-24 | 南京邮电大学 | 基于r-cnn的眼底图像微动脉瘤识别方法 |
US20200085290A1 (en) * | 2017-05-04 | 2020-03-19 | Shenzhen Sibionics Technology Co., Ltd. | Artificial neural network and system for identifying lesion in retinal fundus image |
US20200160521A1 (en) * | 2017-05-04 | 2020-05-21 | Shenzhen Sibionics Technology Co., Ltd. | Diabetic retinopathy recognition system based on fundus image |
CN111259680A (zh) * | 2020-02-13 | 2020-06-09 | 支付宝(杭州)信息技术有限公司 | 二维码图像二值化处理方法及装置 |
WO2020140198A1 (zh) * | 2019-01-02 | 2020-07-09 | 深圳市邻友通科技发展有限公司 | 指甲图像分割方法、装置、设备和存储介质 |
WO2020199773A1 (zh) * | 2019-04-04 | 2020-10-08 | 京东方科技集团股份有限公司 | 图像检索方法及装置和计算机可读存储介质 |
CN111914874A (zh) * | 2020-06-09 | 2020-11-10 | 上海欣巴自动化科技股份有限公司 | 一种目标检测方法和系统 |
-
2021
- 2021-07-26 CN CN202110847212.3A patent/CN113506284B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105069803A (zh) * | 2015-08-19 | 2015-11-18 | 西安交通大学 | 一种基于彩色图像的糖尿病病变体微血管瘤的分类装置 |
US20200085290A1 (en) * | 2017-05-04 | 2020-03-19 | Shenzhen Sibionics Technology Co., Ltd. | Artificial neural network and system for identifying lesion in retinal fundus image |
US20200160521A1 (en) * | 2017-05-04 | 2020-05-21 | Shenzhen Sibionics Technology Co., Ltd. | Diabetic retinopathy recognition system based on fundus image |
CN107590941A (zh) * | 2017-09-19 | 2018-01-16 | 重庆英卡电子有限公司 | 拍照式混合火焰探测器及其探测方法 |
WO2020140198A1 (zh) * | 2019-01-02 | 2020-07-09 | 深圳市邻友通科技发展有限公司 | 指甲图像分割方法、装置、设备和存储介质 |
WO2020199773A1 (zh) * | 2019-04-04 | 2020-10-08 | 京东方科技集团股份有限公司 | 图像检索方法及装置和计算机可读存储介质 |
CN109977930A (zh) * | 2019-04-29 | 2019-07-05 | 中国电子信息产业集团有限公司第六研究所 | 疲劳驾驶检测方法及装置 |
CN110276356A (zh) * | 2019-06-18 | 2019-09-24 | 南京邮电大学 | 基于r-cnn的眼底图像微动脉瘤识别方法 |
CN111259680A (zh) * | 2020-02-13 | 2020-06-09 | 支付宝(杭州)信息技术有限公司 | 二维码图像二值化处理方法及装置 |
CN111914874A (zh) * | 2020-06-09 | 2020-11-10 | 上海欣巴自动化科技股份有限公司 | 一种目标检测方法和系统 |
Non-Patent Citations (2)
Title |
---|
ORLANDO J I等: ""An ensemble deep learning based approach for red lesion detection in fundus images"" * |
刘尚平等: ""荧光视网膜图像的照度均衡及自适应血管增强算法"" * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114529724A (zh) * | 2022-02-15 | 2022-05-24 | 推想医疗科技股份有限公司 | 图像目标的识别方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113506284B (zh) | 2023-05-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110276356B (zh) | 基于r-cnn的眼底图像微动脉瘤识别方法 | |
Wang et al. | Hard exudate detection based on deep model learned information and multi-feature joint representation for diabetic retinopathy screening | |
CN109472781B (zh) | 一种基于串行结构分割的糖尿病视网膜病变检测系统 | |
Zhao et al. | Retinal vessels segmentation based on level set and region growing | |
Yavuz et al. | Blood vessel extraction in color retinal fundus images with enhancement filtering and unsupervised classification | |
Tian et al. | Multi-path convolutional neural network in fundus segmentation of blood vessels | |
Kande et al. | Segmentation of exudates and optic disk in retinal images | |
CN108986106A (zh) | 面向青光眼临床诊断的视网膜血管自动分割方法 | |
Pathan et al. | Automated segmentation and classifcation of retinal features for glaucoma diagnosis | |
Solís-Pérez et al. | Blood vessel detection based on fractional Hessian matrix with non-singular Mittag–Leffler Gaussian kernel | |
Mahapatra et al. | A novel framework for retinal vessel segmentation using optimal improved frangi filter and adaptive weighted spatial FCM | |
David et al. | A Comprehensive Review on Partition of the Blood Vessel and Optic Disc in RetinalImages | |
Jayanthi et al. | Automatic diagnosis of retinal diseases from color retinal images | |
Sharma et al. | Machine learning approach for detection of diabetic retinopathy with improved pre-processing | |
Senapati | Bright lesion detection in color fundus images based on texture features | |
Lyu et al. | Deep tessellated retinal image detection using Convolutional Neural Networks | |
Kanca et al. | Learning hand-crafted features for k-NN based skin disease classification | |
Athira et al. | Automatic detection of diabetic retinopathy using R-CNN | |
CN113011340B (zh) | 一种基于视网膜图像的心血管手术指标风险分类方法及系统 | |
CN113506284A (zh) | 一种眼底图像微血管瘤检测装置、方法及存储介质 | |
Syed et al. | Detection of tumor in MRI images using artificial neural networks | |
Gou et al. | A novel retinal vessel extraction method based on dynamic scales allocation | |
Jana et al. | A semi-supervised approach for automatic detection and segmentation of optic disc from retinal fundus image | |
Purwanithami et al. | Hemorrhage diabetic retinopathy detection based on fundus image using neural network and FCM segmentation | |
CN113269756A (zh) | 基于多尺度匹配滤波与粒子群优化视网膜血管分割方法、装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |