CN114782452B - 一种荧光素眼底血管造影图像的处理系统及装置 - Google Patents
一种荧光素眼底血管造影图像的处理系统及装置 Download PDFInfo
- Publication number
- CN114782452B CN114782452B CN202210715436.3A CN202210715436A CN114782452B CN 114782452 B CN114782452 B CN 114782452B CN 202210715436 A CN202210715436 A CN 202210715436A CN 114782452 B CN114782452 B CN 114782452B
- Authority
- CN
- China
- Prior art keywords
- image
- area
- fluorescein fundus
- retinal
- fluorescein
- 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
Links
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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- 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/10—Image acquisition modality
- G06T2207/10064—Fluorescence image
-
- 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/30041—Eye; Retina; Ophthalmic
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Molecular Biology (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Geometry (AREA)
- Computational Linguistics (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本申请属于计算机视觉技术领域,公开了一种荧光素眼底血管造影图像的处理系统及装置,该系统包括:图像预处理模块,用于获取待处理的荧光素眼底血管造影图像,并对荧光素眼底血管造影图像进行预处理;病变区域及病变区域内无灌注区分割模块,用于通过预先训练的语义分割模型对经过预处理的荧光素眼底血管造影图像进行处理,确定其中的病变区域和病变区域内的无灌注区;缺血指数计算模块,用于获取病变区域的面积值和无灌注区的面积值,并计算得到荧光素眼底血管造影图像对应的适用于临床的缺血指数。本申请可以达到分割荧光素眼底血管造影图像中的病变区域及病变区域中的无灌注区、实现病变量化以及适用于多种视网膜病变的效果。
Description
技术领域
本申请涉及计算机视觉技术领域,尤其涉及一种荧光素眼底血管造影图像的处理系统及装置。
背景技术
视网膜病变有较多分类,比如糖尿病视网膜病变(diabetic retinopathy, DR)、视网膜分支静脉阻塞(Branch Retinal Vein Occlusion,BRVO)、视网膜中央静脉阻塞(Central Retinal Vein Occlusion,CRVO)及视网膜血管炎等;荧光素眼底血管造影(Fundus Fluorescein Angiography,FFA)可以反映视网膜血管及灌注情况。然而,现有技术中只有针对糖尿病视网膜病变(DR)中无灌注区的分割和定位,存在不仅未实现病变的量化,而且无法对其他视网膜病变进行处理,适用性差的问题。
发明内容
本申请提供了一种荧光素眼底血管造影图像的处理系统及装置,可以分割荧光素眼底血管造影图像中的病变区域及病变区域中的无灌注区,实现病变量化且适用于多种视网膜病变。
第一方面,本申请实施例提供了一种荧光素眼底血管造影图像的处理系统,该系统包括:
图像预处理模块,用于获取待处理的荧光素眼底血管造影图像,并对荧光素眼底血管造影图像进行预处理,荧光素眼底血管造影图像为预设的多种视网膜病变中任一种视网膜病变的荧光素眼底血管造影图像;
病变区域及病变区域内无灌注区分割模块,用于通过预先训练的语义分割模型对经过预处理的荧光素眼底血管造影图像进行处理,确定其中的病变区域和病变区域内的无灌注区,语义分割模型是基于多种视网膜病变分别对应的荧光素眼底血管造影图像训练而得到的;
缺血指数计算模块,用于获取病变区域的面积值和无灌注区的面积值,并计算得到荧光素眼底血管造影图像对应的适用于临床的缺血指数。
基于上述实施例,在一些实施方式中,预设的多种视网膜病变包括糖尿病视网膜病变、视网膜分支静脉阻塞、视网膜中央静脉阻塞及视网膜血管炎;
在荧光素眼底血管造影图像为糖尿病视网膜病变的荧光素眼底血管造影图像时,病变区域分割模块确定的病变区域为荧光素眼底血管造影图像中的整个视网膜区域;
在荧光素眼底血管造影图像为视网膜分支静脉阻塞、视网膜中央静脉阻塞或视网膜血管炎的荧光素眼底血管造影图像时,病变区域及病变区域内无灌注区分割模块确定的病变区域为荧光素眼底血管造影图像中的视网膜分支静脉阻塞、视网膜中央静脉阻塞或视网膜血管炎的病变区域。
进一步地,该系统还包括模型训练模块,该模型训练模块包括:
原始图像采集单元,用于采集原始图像以创建原始图像集,原始图像为糖尿病视网膜病变或视网膜分支静脉阻塞的荧光素眼底血管造影图像,原始图像集中包括糖尿病视网膜病变和视网膜分支静脉阻塞的荧光素眼底血管造影图像;
原始图像预处理单元,用于对原始图像集中的每张原始图像进行预处理,并从原始图像集中删除不符合预设标准的原始图像,得到已预处理图像集;
图像标注单元,用于对已预处理图像集中的每张已预处理图像进行标注,得到标注图像集;其中,在已预处理图像为糖尿病视网膜病变的荧光素眼底血管造影图像时,对无灌注区进行标注;在已预处理图像为视网膜分支静脉阻塞的荧光素眼底血管造影图像时,分别对无灌注区和视网膜分支静脉阻塞病变的区域进行标注;
图像分配单元,用于按照预设比例将标注图像集中的各张标注图像随机分配到训练集和第一测试集;
模型构建单元,用于构建一个基于卷积神经网络的语义分割模型;
模型训练单元,用于通过训练集中的标注图像对语义分割模型进行训练;
模型验证单元,用于获取第二测试集,并分别使用第一测试集和第二测试集中的图像对训练后的语义分割模型进行验证,以得到预先训练的语义分割模型;第二测试集中包括多种视网膜病变中除糖尿病视网膜病变和视网膜分支静脉阻塞以外的至少一种视网膜病变的荧光素眼底血管造影图像。
进一步地,该系统还包括图像显示模块;
图像显示模块用于显示荧光素眼底血管造影图像,并在荧光素眼底血管造影图像中使用不同的指定颜色显示病变区域和无灌注区;和/或,图像显示模块用于显示病变区域的面积值、无灌注区的面积值及荧光素眼底血管造影图像对应的适用于临床的缺血指数。
进一步地,荧光素眼底血管造影图像对应的适用于临床的缺血指数为荧光素眼底血管造影图像中的无灌注区面积占病变区域面积的比例值。
进一步地,预设的多种视网膜病变包括糖尿病视网膜病变、视网膜分支静脉阻塞、视网膜中央静脉阻塞以及视网膜血管炎中的至少两种。
第二方面,本申请另一实施例提供了一种荧光素眼底血管造影图像的处理装置,该装置包括计算机设备,计算机设备包括处理器和存储器,存储器中存储有计算机程序,处理器调用存储器中的计算机程序时执行如下步骤:
获取待处理的荧光素眼底血管造影图像,并对荧光素眼底血管造影图像进行预处理,荧光素眼底血管造影图像为预设的多种视网膜病变中任一种视网膜病变的荧光素眼底血管造影图像;
通过预先训练的语义分割模型对经过预处理的荧光素眼底血管造影图像进行处理,确定其中的病变区域和病变区域内的无灌注区,语义分割模型是基于多种视网膜病变分别对应的荧光素眼底血管造影图像训练而得到的;
获取病变区域的面积值和无灌注区的面积值,并计算得到荧光素眼底血管造影图像对应的适用于临床的缺血指数。
进一步地,处理器调用存储器中的计算机程序时还执行如下步骤:
采集原始图像以创建原始图像集,原始图像为糖尿病视网膜病变或视网膜分支静脉阻塞的荧光素眼底血管造影图像,原始图像集中包括糖尿病视网膜病变和视网膜分支静脉阻塞的荧光素眼底血管造影图像;
对原始图像集中的每张原始图像进行预处理,并从原始图像集中删除不符合预设标准的原始图像,得到已预处理图像集;
对已预处理图像集中的每张已预处理图像进行标注,得到标注图像集;其中,在已预处理图像为糖尿病视网膜病变的荧光素眼底血管造影图像时,对无灌注区进行标注;在已预处理图像为视网膜分支静脉阻塞的荧光素眼底血管造影图像时,分别对无灌注区和视网膜分支静脉阻塞病变的区域进行标注;
按照预设比例将标注图像集中的各张标注图像随机分配到训练集和第一测试集;
构建一个基于卷积神经网络的语义分割模型;
通过训练集中的标注图像对语义分割模型进行训练;
获取第二测试集,并分别使用第一测试集和第二测试集中的图像对训练后的语义分割模型进行验证,以得到预先训练的语义分割模型;第二测试集中包括多种视网膜病变中除糖尿病视网膜病变和视网膜分支静脉阻塞以外的至少一种视网膜病变的荧光素眼底血管造影图像。
进一步地,处理器调用存储器中的计算机程序时还执行如下步骤:
显示荧光素眼底血管造影图像,并在荧光素眼底血管造影图像中使用不同的指定颜色显示病变区域和无灌注区;和/或,显示病变区域的面积值、无灌注区的面积值及荧光素眼底血管造影图像对应的适用于临床的缺血指数。
第三方面,本申请又一实施例提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现如下步骤:
获取待处理的荧光素眼底血管造影图像,并对荧光素眼底血管造影图像进行预处理,荧光素眼底血管造影图像为预设的多种视网膜病变中任一种视网膜病变的荧光素眼底血管造影图像;
通过预先训练的语义分割模型对经过预处理的荧光素眼底血管造影图像进行处理,确定其中的病变区域和病变区域内的无灌注区,语义分割模型是基于多种视网膜病变分别对应的荧光素眼底血管造影图像训练而得到的;
获取病变区域的面积值和无灌注区的面积值,并计算得到荧光素眼底血管造影图像对应的适用于临床的缺血指数。
综上所述,与现有技术相比,本申请提供的技术方案带来的有益效果至少包括:
本申请提供的一种荧光素眼底血管造影图像的处理系统,可以通过图像预处理模块获取待处理的荧光素眼底血管造影图像,并对荧光素眼底血管造影图像进行预处理;然后通过语义分割模型对经过预处理的荧光素眼底血管造影图像进行处理,确定其中的病变区域和病变区域内的无灌注区,并获取病变区域的面积值和无灌注区的面积值,计算得到荧光素眼底血管造影图像对应的适用于临床的缺血指数,即患者可视范围内的视网膜的缺血程度的值,提供了更直观的缺血指数以供临床参考。其中,语义分割模型是基于多种视网膜病变的荧光素眼底血管造影图像训练而得到的,因此该系统能适用于多种视网膜病变的荧光素眼底血管造影图像处理,具备更好的适用性。
附图说明
图1为本申请一个示例性实施例提供的荧光素眼底血管造影图像的处理系统的结构图。
图2为本申请另一个示例性实施例提供的荧光素眼底血管造影图像的处理系统的结构图。
图3为本申请另一个示例性实施例提供的模型训练模块的结构图。
图4为本申请又一个示例性实施例提供的荧光素眼底血管造影图像的处理系统的结构图。
图5为本申请一个示例性实施例提供的荧光素眼底血管造影图像的处理装置执行步骤的流程图。
图6为本申请一个示例性实施例提供的荧光素眼底血管造影图像的处理装置执行的模型训练步骤的流程图。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。
实施例一:
请参见图1,本申请的一个实施例提供了一种荧光素眼底血管造影图像的处理系统,该系统包括:
图像预处理模块100,用于获取待处理的荧光素眼底血管造影图像,并对荧光素眼底血管造影图像进行预处理,荧光素眼底血管造影图像为预设的多种视网膜病变中任一种视网膜病变的荧光素眼底血管造影图像。
其中,待处理的荧光素眼底血管造影图像可采用55度的荧光素眼底血管造影图像(FFA图像);预设的多种视网膜病变可以包括糖尿病视网膜病变(DR)、视网膜分支静脉阻塞(BRVO)、视网膜中央静脉阻塞(CRVO)以及视网膜血管炎中的至少两种,因此,待处理的荧光素眼底血管造影图像可以是上述任一种视网膜病变的患者的FFA图像。
具体地,图像预处理模块100与病变区域及病变区域内无灌注区分割模块200连接,图像预处理模块100可以是眼底照相机,将拍摄的FFA图像直接发送到病变区域及病变区域内无灌注区分割模块200进行处理;图像预处理模块100也可以是图像输入设备,例如与眼底照相机连接的图像输入数据接口,该图像输入数据接口从眼底照相机处获取到待处理的FFA图像,再将获取到的FFA图像发送给病变区域及病变区域内无灌注区分割模块200。
病变区域及病变区域内无灌注区分割模块200,用于通过预先训练的语义分割模型对经过预处理的荧光素眼底血管造影图像进行处理,确定其中的病变区域和病变区域内的无灌注区,语义分割模型是基于多种视网膜病变对应的荧光素眼底血管造影图像训练而得到的。
其中,视网膜的无灌注区(non-perfusion,NP),指视网膜血管闭塞导致相应灌注的区域视网膜没有视功能,NP的出现反映视网膜处在缺血、缺氧的状态;NP常见于DR、BRVO、CRVO及视网膜血管炎等视网膜血管性疾病,NP位于病变区域内。在临床上,NP可用于衡量视网膜疾病的缺血程度。因此,病变区域及病变区域内无灌注区分割模块200通过预先训练的语义分割模型对病变区域内的NP也进行了定位和分割,病变区域及病变区域内无灌注区分割模块200可以使用计算机设备,计算机设备包括但不限于服务器、台式机及笔记本电脑等。
具体地,病变区域及病变区域内无灌注区分割模块200内部设有预先训练的语义分割模型,该语义分割模型基于多种视网膜病变的FFA图像训练而得到,上述多种视网膜病变包括但不限于DR、BRVO、CRVO及视网膜血管炎;病变区域及病变区域内无灌注区分割模块200通过该语义分割模型对获取到的FFA图像进行处理,处理过程包括对FFA图像中病变区域和NP的识别和分割,以确定FFA图像中的病变区域和NP。
缺血指数计算模块300,用于获取病变区域的面积值和无灌注区的面积值,并计算得到荧光素眼底血管造影图像对应的适用于临床的缺血指数。
其中,荧光素眼底血管造影图像对应的适用于临床的缺血指数为荧光素眼底血管造影图像中的无灌注区面积占病变区域面积的比例值,即适用于临床的缺血指数等于无灌注区面积值/病变区域面积值,临床上可以根据适用于临床的缺血指数辅助评估患者疾病的严重程度。
具体地,缺血指数计算模块300可以是与病变区域及病变区域内无灌注区分割模块200一体的计算机设备,也可以是只与病变区域及病变区域内无灌注区分割模块200通信连接的另一台计算机设备。
本实施例中,预设的多种视网膜病变包括DR、BRVO、CRVO及视网膜血管炎。
在本实施例的一些实施方式中,待处理的荧光素眼底血管造影图像为DR的FFA图像,此时,病变区域及病变区域内无灌注区分割模块200通过语义分割模型确定的病变区域为FFA图像中的整个视网膜区域,则FFA图像对应的适用于临床的缺血指数等于无灌注区面积/整个视网膜区域。
在本实施例的另一些实施方式中,待处理的荧光素眼底血管造影图像为BRVO、CRVO或视网膜血管炎的FFA图像时,此时,病变区域及病变区域内无灌注区分割模块200通过语义分割模型确定的病变区域为FFA图像中的BRVO、CRVO或视网膜血管炎的病变区域。
上述实施例提供的一种荧光素眼底血管造影图像的处理系统,可以通过图像预处理模块获取待处理的荧光素眼底血管造影图像,并对荧光素眼底血管造影图像进行预处理;然后通过预先训练的语义分割模型对经过预处理的荧光素眼底血管造影图像进行处理,确定其中的病变区域和病变区域内的无灌注区,并获取病变区域的面积值和无灌注区的面积值,计算得到荧光素眼底血管造影图像对应的适用于临床的缺血指数,即患者可视范围内的视网膜的缺血程度的值,提供了更直观的适用于临床的缺血指数以供临床参考。其中,语义分割模型是基于多种视网膜病变的荧光素眼底血管造影图像训练而得到的,因此该系统能适用于多种视网膜病变的荧光素眼底血管造影图像处理,具备更好的适用性。
实施例二:
基于上述实施例一,在本实施例中,请参见图2和图3,该系统还包括模型训练模块400,该模型训练模块400包括:
原始图像采集单元402,用于采集原始图像以创建原始图像集。
其中,原始图像可为DR或BRVO的FFA图像,原始图像集中包括DR和BRVO的FFA图像。
具体地,该原始图像采集单元402采集诊断为DR或BRVO的55度的FFA图像作为原始图像,并基于采集到的大量原始图像创建原始图像集。
原始图像预处理单元404,用于对原始图像集中的每张原始图像进行预处理,并从原始图像集中删除不符合预设标准的原始图像,得到已预处理图像集。
其中,预设标准可以为预设的模型输入图像标准,比如要求的输入图像尺寸为512*512,那么将不符合该尺寸要求的原始图像进行删除。
具体地,原始图像预处理单元404通过对原始图像进行预处理,提高样品图像质量,并筛除低质量的图像,得到图像质量更好的已预处理图像集;上述图像预处理可采用任一现有技术中成熟的技术,同时本说明书中未作详细描述的内容均属于本领域技术人员公知的现有技术。
图像标注单元406,用于对已预处理图像集中的每张已预处理图像进行标注,得到标注图像集。
其中,在已预处理图像为DR的FFA图像时,对NP进行标注;在已预处理图像为BRVO的FFA图像时,分别对NP和BRVO病变的区域进行标注。
具体地,图像标注单元406可以采用深度学习图像标注工具Labelme对已预处理图像集中的每张已预处理图像进行标注,标注内容可以包括NP和病变区域,对于诊断为DR的已预处理图像,只标注该已预处理图像中的NP;对于诊断为BRVO的已预处理图像,标注该已预处理图像中的NP和病变区域;已预处理图像的标注结果为JSON格式,存储内容为每个标注区域的标签名及该标注区域边缘点的坐标位置信息。因此,在完成对已预处理图像的标注后,需要通过常用的json2img工具将标注结果转换为图像格式,在转换过程中,可以将背景区域的像素值设置为0,NP的像素值设置为1,BRVO病变区域的像素值设置为2,背景区域指的是既非NP也非BRVO病变区域的视网膜区域。
图像分配单元408,用于按照预设比例将标注图像集中的各张标注图像随机分配到训练集和第一测试集。
具体地,图像分配单元408可以按照80%和20%的占比将已预处理图像集中标注后的各张标注图像随机分配到训练集和第一测试集,除了8:2,上述预设比例还可以是7:3、9:1等任意大于1的比例值。
模型构建单元410,用于构建一个基于卷积神经网络的语义分割模型。
模型训练单元412,用于通过训练集中的标注图像对语义分割模型进行训练。
其中,语义分割模型使用的算法可以为U-Net,以VGG-16为后端,可以定位NP和BRVO病变区域的边界。
在本实施例的一种具体实施方式中,构建的语义分割模型可以包括以数据传递顺序依次连接的输入层、下采样层、上采样层和输出层。
其中,输入层要求输入图像的尺寸为512*512。
下采样层包括13层卷积层,具有特征提取功能;具体地,下采样层可以包括:第1、2层为卷积层,卷积核数量为64、卷积核大小为3×3,步长为1,padding方式为same,且在卷积层后会有一个relu激活层以提高网络模型的非线性拟合能力;第3层为最大值池化层,窗口大小为2×2,步长为2,池化层输出的特征图大小会变小;第4、5层为卷积层,为保证特征数量不减少,卷积核数量会增加变为128;第6层为池化层,参数配置与第3层相同的;第7、8、9层为卷积层,卷积核数量为256;第10层为最大值池化层,参数配置与第3层相同;第11、12、13为卷积层,卷积核数量为512;第14层为池化层,参数配置与第3层相同的;第15、16、17为卷积层,卷积核数量为512;第18层为池化层,参数配置与第3层相同的;第19层为drop层,随机drop率为0.5。
上采样层起到的作用是还原图像细节、提取目标区域;具体地,上采样层可以包括:第20层为转置卷积层,卷积核大小为2×2,该层可将特征图大小放大;第21层为卷积层,随着特征图大小的放大,卷积核数量也会减少,这层的卷积核数量为512,卷积核大小为2×2,该层的输出会与下采样层第14层的输出特征图进行拼接形成最终的输出,通过拼接可保留更多的图像细节,提升病灶区域分割结果的精细度;第22、23层为卷积层,随着特征图大小的放大,卷积核数量也会减少,这两层的卷积核数量为512,卷积核大小为3×3;第24层为转置卷积层,参数与第19层相同;第25层为卷积层,卷积核大小为2×2,卷积核数量为256,该层的输出会与下采样层第10层的输出特征图进行拼接形成最终的输出;第26、27层为卷积层,卷积核数量为256,卷积核大小为3×3;第28层为转置卷积层,参数与第18层相同;第29层为卷积层,卷积核大小为2×2,卷积核数量为128,该层的输出会与下采样层第6层的输出特征图进行拼接生成该层的输出;第30、31层为卷积层,卷积核数量为128,卷积核大小为3×3;第32层为转置卷积层,参数与第18层相同;第33层为卷积层,卷积核大小为2×2,卷积核数量为64,该层的输出会与下采样层第3层的输出特征图进行拼接生成该层的输出;第34、35层为卷积层,卷积核数量为64,卷积核大小为3×3。
输出层用于输出最终分割结果,输出层具体可以包括:第36层为转置卷积层,参数与第18层相同;第37层为卷积层,卷积核数量为类别数目2;38层为softmax层,对输出结果进行归一化,得到造影图像每个像素分别属于每个类别的概率。
上述的语义分割模型各种层的搭建是通过keras框架中的layers模块实现的,其中,卷积层的构建方法为layers.Conv2D,池化层的构建方法为layers.MaxPooling2D,转置卷积层的构建方法为layers.Upsampling2D,高层特征与低层特征的拼接采用layers.Concatenate方法。
模型验证单元414,用于获取第二测试集,并分别使用第一测试集和第二测试集中的图像对训练后的语义分割模型进行验证,以得到预先训练的语义分割模型,第二测试集中包括多种视网膜病变中除糖尿病视网膜病变和视网膜分支静脉阻塞以外的至少一种视网膜病变的荧光素眼底血管造影图像。
其中,模型验证单元414可以利用第一测试集和第二测试集中的图像作为训练数据,对多流语义分割模型进行优化训练,具体包括:在语义分割模型中分批次加载训练数据,得到该模型的输出结果后,通过对模型输出结果进行计算获得每批次数据的准确率与交叉熵损失;在训练过程中不断保存更优效果的网络模型权重文件,通过反向传播算法利用交叉熵损失对网络模型权重进行优化,使网络模型的损失值不断降低;在训练迭代次数达到程序设定值时退出训练循环,得到预先训练的语义分割模型;使用该语义分割模型对荧光素眼底血管造影图像的NP和BRVO病变区域进行处理时,选取像素点的预测值为1的区域作为NP候选区,选取像素点的预测值为2的区域作为BRVO病变区域候选区。
在上述实施例中,语义分割模型是基于多种视网膜病变分别对应的荧光素眼底血管造影图像训练而得到的,因此上述系统能够适用于多种视网膜病变的荧光素眼底血管造影图像处理,适用性好,应用范围广。
实施例三:
基于上述实施例二,在本实施例中,请参见图4,该系统还包括图像显示模块500。
图像显示模块500可以用于显示荧光素眼底血管造影图像,并在荧光素眼底血管造影图像中使用不同的指定颜色显示病变区域和无灌注区,从而实现了无灌注区的直观显示,有助于协助医生确定激光治疗的具体位置,实施精准医疗,也能使激光手术简单化,促进基层医院也具备开展激光治疗的能力;和/或,图像显示模块500用于显示病变区域的面积值、无灌注区的面积值以及荧光素眼底血管造影图像对应的适用于临床的缺血指数,可以自动显示患者视网膜的病变区域面积值、无灌注区面积值以及适用于临床的缺血指数,可帮助医生直观地了解患者视网膜的缺血程度。
实施例四:
请参见图5,本实施例提供了一种荧光素眼底血管造影图像的处理装置,该装置包括计算机设备,该计算机设备包括处理器和存储器,存储器中存储有计算机程序,处理器调用存储器中的计算机程序时执行如下步骤:
步骤S1,获取待处理的荧光素眼底血管造影图像,并对荧光素眼底血管造影图像进行预处理。
其中,荧光素眼底血管造影图像为预设的多种视网膜病变中任一种视网膜病变的荧光素眼底血管造影图像。
步骤S2,通过预先训练的语义分割模型对经过预处理的荧光素眼底血管造影图像进行处理,确定其中的病变区域和病变区域内的无灌注区。
其中,语义分割模型是基于多种视网膜病变分别对应的荧光素眼底血管造影图像训练而得到的。
步骤S3,获取病变区域的面积值和无灌注区的面积值,并计算得到荧光素眼底血管造影图像对应的适用于临床的缺血指数。
具体地,处理器获取所述病变区域的面积值和所述无灌注区的面积值,并根据所述病变区域的面积值和所述无灌注区的面积值计算得到荧光素眼底血管造影图像对应的适用于临床的缺血指数。
对于本实施例中步骤S1~S3的说明,请参见上述各系统实施例中对于图像预处理模块100、病变区域及病变区域内无灌注区分割模块200及缺血指数计算模块300的说明内容,在此不再加以赘述。
在本实施例的一些实施方式中,所述计算机程序被处理器执行时还实现模型训练步骤,请参见图6,模型训练步骤具体包括:
步骤S41,采集原始图像以创建原始图像集。
其中,原始图像为DR或BRVO的FFA图像,原始图像集中包括DR和BRVO的FFA图像。
步骤S42,对原始图像集中的每张原始图像进行预处理,并从原始图像集中删除不符合预设标准的原始图像,得到已预处理图像集。
步骤S43,对已预处理图像集中的每张已预处理图像进行标注,得到标注图像集。
其中,在已预处理图像为DR的FFA图像时,对NP进行标注;在已预处理图像为BRVO的FFA图像时,分别对NP和BRVO的病变区域进行标注。
步骤S44,按照预设比例将标注图像集中的各张标注图像随机分配到训练集和第一测试集。
步骤S45,构建一个基于卷积神经网络的语义分割模型。
步骤S46,通过训练集中的标注图像对语义分割模型进行训练。
步骤S47,获取第二测试集,并分别使用第一测试集和第二测试集中的图像对训练后的语义分割模型进行验证,以得到预先训练的语义分割模型。
其中,第二测试集中包括多种视网膜病变中DR和BRVO以外的至少一种视网膜病变的荧光素眼底血管造影图像。
对于本实施例中步骤S41~S47的说明,请参见上述各系统实施例中对于模型训练模块400的说明内容,在此不再加以赘述。
在本实施例的一些实施方式中,所述计算机程序被处理器执行时还实现如下步骤:显示荧光素眼底血管造影图像,并在荧光素眼底血管造影图像中使用不同的指定颜色显示病变区域和无灌注区;和/或,显示病变区域的面积值、无灌注区的面积值以及荧光素眼底血管造影图像对应的适用于临床的缺血指数。
对于本实施例中显示步骤的说明,请参见上述各系统实施例中对于图像显示模块500的说明内容,在此不再加以赘述。
具体地,该计算机设备可以是终端,处理器用于提供计算和控制能力,存储器包括非易失性存储介质和内存储器,非易失性存储介质存储有操作系统和计算机程序,内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。
实施例五:
本申请一实施例提供了一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现如上述实施例中荧光素眼底血管造影图像的处理装置执行的步骤。
其中,所述计算机可读存储介质是指存储数据的载体,可以但不限于包括软盘、光盘、硬盘、闪存、优盘和/或记忆棒(Memory Stick)等,所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。
本实施例中提供的计算机可读存储介质的工作过程、工作细节和技术效果,可以参见上文中关于荧光素眼底血管造影图像的处理系统的实施例,于此不再赘述。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于装置或装置实施例而言,由于其基本相似于系统实施例,所以描述得比较简单,相关之处参见系统实施例的部分说明即可。以上所描述的系统及系统实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性劳动的情况下,即可以理解并实施。
以上对本申请所提供的技术方案进行了详细介绍,本文中应用了具体个例对本申请的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本申请的方案及其核心思想;同时,对于本领域的一般技术人员,依据本申请的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本申请的限制。
以上所述仅为本申请的较佳实施例,并不用以限制本申请,凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
Claims (8)
1.一种荧光素眼底血管造影图像的处理系统,其特征在于,所述系统包括:
图像预处理模块,用于获取待处理的荧光素眼底血管造影图像,并对所述荧光素眼底血管造影图像进行预处理,所述荧光素眼底血管造影图像为预设的多种视网膜病变中任一种视网膜病变的55度的荧光素眼底血管造影图像;所述预设的多种视网膜病变包括糖尿病视网膜病变、视网膜分支静脉阻塞、视网膜中央静脉阻塞以及视网膜血管炎;
病变区域及病变区域内无灌注区分割模块,用于通过预先训练的语义分割模型对经过预处理的所述荧光素眼底血管造影图像进行处理,确定其中的病变区域和位于所述病变区域内的无灌注区,所述语义分割模型是基于所述多种视网膜病变分别对应的荧光素眼底血管造影图像训练而得到的;
缺血指数计算模块,用于获取所述病变区域的面积值和所述无灌注区的面积值,并计算得到所述荧光素眼底血管造影图像对应的适用于临床的缺血指数;
在所述荧光素眼底血管造影图像为糖尿病视网膜病变的荧光素眼底血管造影图像时,所述病变区域及病变区域内无灌注区分割模块确定的病变区域为所述荧光素眼底血管造影图像中的整个视网膜区域;
在所述荧光素眼底血管造影图像为视网膜分支静脉阻塞、视网膜中央静脉阻塞或视网膜血管炎的荧光素眼底血管造影图像时,所述病变区域及病变区域内无灌注区分割模块确定的病变区域为所述荧光素眼底血管造影图像中的视网膜分支静脉阻塞、视网膜中央静脉阻塞或视网膜血管炎的病变区域。
2.根据权利要求1所述的系统,其特征在于,所述系统还包括模型训练模块,所述模型训练模块包括:
原始图像采集单元,用于采集原始图像以创建原始图像集,所述原始图像为糖尿病视网膜病变或视网膜分支静脉阻塞的荧光素眼底血管造影图像,所述原始图像集中包括糖尿病视网膜病变和视网膜分支静脉阻塞的荧光素眼底血管造影图像;
原始图像预处理单元,用于对所述原始图像集中的每张原始图像进行预处理,并从所述原始图像集中删除不符合预设标准的原始图像,得到已预处理图像集;
图像标注单元,用于对所述已预处理图像集中的每张已预处理图像进行标注,得到标注图像集;其中,在已预处理图像为糖尿病视网膜病变的荧光素眼底血管造影图像时,对无灌注区进行标注;在已预处理图像为视网膜分支静脉阻塞的荧光素眼底血管造影图像时,分别对无灌注区和视网膜分支静脉阻塞病变的区域进行标注;
图像分配单元,用于按照预设比例将所述标注图像集中的各张标注图像随机分配到训练集和第一测试集;
模型构建单元,用于构建一个基于卷积神经网络的语义分割模型;
模型训练单元,用于通过所述训练集中的标注图像对所述语义分割模型进行训练;
模型验证单元,用于获取第二测试集,并分别使用所述第一测试集和所述第二测试集中的图像对训练后的语义分割模型进行验证,以得到所述预先训练的语义分割模型;所述第二测试集中包括所述多种视网膜病变中除糖尿病视网膜病变和视网膜分支静脉阻塞以外的至少一种视网膜病变的荧光素眼底血管造影图像。
3.根据权利要求1或2任一项所述的系统,其特征在于,所述系统还包括图像显示模块;
所述图像显示模块用于显示所述荧光素眼底血管造影图像,并在所述荧光素眼底血管造影图像中使用不同的指定颜色显示所述病变区域和所述无灌注区;和/或,所述图像显示模块用于显示所述病变区域的面积值、所述无灌注区的面积值及所述荧光素眼底血管造影图像对应的适用于临床的缺血指数。
4.根据权利要求1所述的系统,其特征在于,所述荧光素眼底血管造影图像对应的适用于临床的缺血指数为所述荧光素眼底血管造影图像中的无灌注区面积占病变区域面积的比例值。
5.一种荧光素眼底血管造影图像的处理装置,所述装置包括计算机设备,所述计算机设备包括处理器和存储器,所述存储器中存储有计算机程序,其特征在于,所述处理器调用所述存储器中的计算机程序时执行如下步骤:
获取待处理的荧光素眼底血管造影图像,并对所述荧光素眼底血管造影图像进行预处理,所述荧光素眼底血管造影图像为预设的多种视网膜病变中任一种视网膜病变的55度的荧光素眼底血管造影图像;所述预设的多种视网膜病变包括糖尿病视网膜病变、视网膜分支静脉阻塞、视网膜中央静脉阻塞以及视网膜血管炎;
通过预先训练的语义分割模型对经过预处理的所述荧光素眼底血管造影图像进行处理,确定其中的病变区域和位于所述病变区域内的无灌注区,所述语义分割模型是基于所述多种视网膜病变分别对应的荧光素眼底血管造影图像训练而得到的;
获取所述病变区域的面积值和所述无灌注区的面积值,并计算得到所述荧光素眼底血管造影图像对应的适用于临床的缺血指数;
在所述荧光素眼底血管造影图像为糖尿病视网膜病变的荧光素眼底血管造影图像时,所述病变区域及病变区域内无灌注区分割模块确定的病变区域为所述荧光素眼底血管造影图像中的整个视网膜区域;
在所述荧光素眼底血管造影图像为视网膜分支静脉阻塞、视网膜中央静脉阻塞或视网膜血管炎的荧光素眼底血管造影图像时,所述病变区域及病变区域内无灌注区分割模块确定的病变区域为所述荧光素眼底血管造影图像中的视网膜分支静脉阻塞、视网膜中央静脉阻塞或视网膜血管炎的病变区域。
6.根据权利要求5所述的装置,其特征在于,所述处理器调用所述存储器中的计算机程序时还执行如下步骤:
采集原始图像以创建原始图像集,所述原始图像为糖尿病视网膜病变或视网膜分支静脉阻塞的荧光素眼底血管造影图像,所述原始图像集中包括糖尿病视网膜病变和视网膜分支静脉阻塞的荧光素眼底血管造影图像;
对所述原始图像集中的每张原始图像进行预处理,并从所述原始图像集中删除不符合预设标准的原始图像,得到已预处理图像集;
对所述已预处理图像集中的每张已预处理图像进行标注,得到标注图像集;其中,在已预处理图像为糖尿病视网膜病变的荧光素眼底血管造影图像时,对无灌注区进行标注;在已预处理图像为视网膜分支静脉阻塞的荧光素眼底血管造影图像时,分别对无灌注区和视网膜分支静脉阻塞病变的区域进行标注;
按照预设比例将所述标注图像集中的各张标注图像随机分配到训练集和第一测试集;
构建一个基于卷积神经网络的语义分割模型;
通过所述训练集中的标注图像对所述语义分割模型进行训练;
获取第二测试集,并分别使用所述第一测试集和所述第二测试集中的图像对训练后的语义分割模型进行验证,以得到所述预先训练的语义分割模型;所述第二测试集中包括所述多种视网膜病变中除糖尿病视网膜病变和视网膜分支静脉阻塞以外的至少一种视网膜病变的荧光素眼底血管造影图像。
7.根据权利要求5或6所述的装置,其特征在于,所述处理器调用所述存储器中的计算机程序时还执行如下步骤:
显示所述荧光素眼底血管造影图像,并在所述荧光素眼底血管造影图像中使用不同的指定颜色显示所述病变区域和所述无灌注区;和/或,显示所述病变区域的面积值、所述无灌注区的面积值以及所述荧光素眼底血管造影图像对应的适用于临床的缺血指数。
8.一种计算机可读存储介质,其特征在于,其上存储有计算机程序,所述计算机程序被处理器执行时实现如下步骤:
获取待处理的荧光素眼底血管造影图像,并对所述荧光素眼底血管造影图像进行预处理,所述荧光素眼底血管造影图像为预设的多种视网膜病变中任一种视网膜病变的55度的荧光素眼底血管造影图像;所述预设的多种视网膜病变包括糖尿病视网膜病变、视网膜分支静脉阻塞、视网膜中央静脉阻塞以及视网膜血管炎;
通过预先训练的语义分割模型对经过预处理的所述荧光素眼底血管造影图像进行处理,确定其中的病变区域和位于所述病变区域内的无灌注区,所述语义分割模型是基于所述多种视网膜病变分别对应的荧光素眼底血管造影图像训练而得到的;
获取所述病变区域的面积值和所述无灌注区的面积值,并计算得到所述荧光素眼底血管造影图像对应的适用于临床的缺血指数;
在所述荧光素眼底血管造影图像为糖尿病视网膜病变的荧光素眼底血管造影图像时,所述病变区域及病变区域内无灌注区分割模块确定的病变区域为所述荧光素眼底血管造影图像中的整个视网膜区域;
在所述荧光素眼底血管造影图像为视网膜分支静脉阻塞、视网膜中央静脉阻塞或视网膜血管炎的荧光素眼底血管造影图像时,所述病变区域及病变区域内无灌注区分割模块确定的病变区域为所述荧光素眼底血管造影图像中的视网膜分支静脉阻塞、视网膜中央静脉阻塞或视网膜血管炎的病变区域。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210715436.3A CN114782452B (zh) | 2022-06-23 | 2022-06-23 | 一种荧光素眼底血管造影图像的处理系统及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210715436.3A CN114782452B (zh) | 2022-06-23 | 2022-06-23 | 一种荧光素眼底血管造影图像的处理系统及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114782452A CN114782452A (zh) | 2022-07-22 |
CN114782452B true CN114782452B (zh) | 2022-11-01 |
Family
ID=82422307
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210715436.3A Active CN114782452B (zh) | 2022-06-23 | 2022-06-23 | 一种荧光素眼底血管造影图像的处理系统及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114782452B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115330714A (zh) * | 2022-08-10 | 2022-11-11 | 中山大学中山眼科中心 | 一种荧光素眼底血管造影图像的分期和病变诊断系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110120055A (zh) * | 2019-04-12 | 2019-08-13 | 浙江大学 | 基于深度学习的眼底荧光造影图像无灌注区自动分割方法 |
CN112862831A (zh) * | 2021-02-05 | 2021-05-28 | 苏州大学 | 对抗网络生成方法和眼底荧光造影图像无灌注区分割方法 |
CN112957005A (zh) * | 2021-02-01 | 2021-06-15 | 山西省眼科医院(山西省红十字防盲流动眼科医院、山西省眼科研究所) | 眼底造影图像无灌注区自动识别及激光光凝区域推荐算法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107423571B (zh) * | 2017-05-04 | 2018-07-06 | 深圳硅基仿生科技有限公司 | 基于眼底图像的糖尿病视网膜病变识别系统 |
CN109859172A (zh) * | 2019-01-08 | 2019-06-07 | 浙江大学 | 基于眼底造影图像深度学习的糖网病变无灌注区识别方法 |
CN110097502B (zh) * | 2019-04-16 | 2022-03-15 | 深圳市莫廷影像技术有限公司 | 眼底无灌注区的测量方法及装置及图像处理方法 |
-
2022
- 2022-06-23 CN CN202210715436.3A patent/CN114782452B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110120055A (zh) * | 2019-04-12 | 2019-08-13 | 浙江大学 | 基于深度学习的眼底荧光造影图像无灌注区自动分割方法 |
CN112957005A (zh) * | 2021-02-01 | 2021-06-15 | 山西省眼科医院(山西省红十字防盲流动眼科医院、山西省眼科研究所) | 眼底造影图像无灌注区自动识别及激光光凝区域推荐算法 |
CN112862831A (zh) * | 2021-02-05 | 2021-05-28 | 苏州大学 | 对抗网络生成方法和眼底荧光造影图像无灌注区分割方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114782452A (zh) | 2022-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108198184B (zh) | 造影图像中血管分割的方法和系统 | |
CN107680684B (zh) | 用于获取信息的方法及装置 | |
WO2020199593A1 (zh) | 图像分割模型训练方法、图像分割方法、装置、设备及介质 | |
CN110889826B (zh) | 眼部oct图像病灶区域的分割方法、装置及终端设备 | |
CN107729929B (zh) | 用于获取信息的方法及装置 | |
CN111862044B (zh) | 超声图像处理方法、装置、计算机设备和存储介质 | |
CN109829894A (zh) | 分割模型训练方法、oct图像分割方法、装置、设备及介质 | |
CN107665491A (zh) | 病理图像的识别方法及系统 | |
CN108615236A (zh) | 一种图像处理方法及电子设备 | |
WO2022032823A1 (zh) | 图像分割方法、装置、设备及存储介质 | |
CN112017185B (zh) | 病灶分割方法、装置及存储介质 | |
US11954860B2 (en) | Image matching method and device, and storage medium | |
CN110276408B (zh) | 3d图像的分类方法、装置、设备及存储介质 | |
Vij et al. | A systematic review on diabetic retinopathy detection using deep learning techniques | |
CN111325725A (zh) | 视网膜图像识别方法及装置、电子设备、存储介质 | |
CN114782452B (zh) | 一种荧光素眼底血管造影图像的处理系统及装置 | |
CN111986202A (zh) | 青光眼辅助诊断装置、方法及存储介质 | |
US20220181007A1 (en) | Computerized systems for prediction of geographic atrophy progression using deep learning applied to clinical imaging | |
CN113989206A (zh) | 一种基于轻量化模型的骨龄预测方法及装置 | |
CN113256638A (zh) | 基于卷积神经网络模型的眼球角膜神经分割方法及装置 | |
CN111161240A (zh) | 血管分类方法、计算机设备和可读存储介质 | |
CN115330714A (zh) | 一种荧光素眼底血管造影图像的分期和病变诊断系统 | |
CN115829980A (zh) | 一种眼底照片的图像识别方法、装置、设备及存储介质 | |
WO2019171398A1 (en) | A fundus image analysis system | |
CN111462067B (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 |