CN113643354A - 基于增强分辨率的眼底图像的血管管径的测量装置 - Google Patents
基于增强分辨率的眼底图像的血管管径的测量装置 Download PDFInfo
- Publication number
- CN113643354A CN113643354A CN202110951405.3A CN202110951405A CN113643354A CN 113643354 A CN113643354 A CN 113643354A CN 202110951405 A CN202110951405 A CN 202110951405A CN 113643354 A CN113643354 A CN 113643354A
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- measurement
- image
- pixel point
- vessel
- 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
- 210000004204 blood vessel Anatomy 0.000 title claims abstract description 343
- 238000005259 measurement Methods 0.000 claims abstract description 194
- 230000011218 segmentation Effects 0.000 claims abstract description 121
- 238000000605 extraction Methods 0.000 claims abstract description 16
- 238000012549 training Methods 0.000 claims description 37
- 238000000034 method Methods 0.000 claims description 36
- 238000013528 artificial neural network Methods 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000003708 edge detection Methods 0.000 claims 1
- 238000004422 calculation algorithm Methods 0.000 description 26
- 238000002372 labelling Methods 0.000 description 26
- 238000010586 diagram Methods 0.000 description 19
- 238000000691 measurement method Methods 0.000 description 15
- 230000002792 vascular Effects 0.000 description 9
- 238000012545 processing Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 230000003902 lesion Effects 0.000 description 6
- 210000005252 bulbus oculi Anatomy 0.000 description 5
- 210000001508 eye Anatomy 0.000 description 5
- 230000015654 memory Effects 0.000 description 5
- 238000007781 pre-processing Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 239000012528 membrane Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 210000001525 retina Anatomy 0.000 description 4
- 238000004590 computer program Methods 0.000 description 3
- 206010025421 Macule Diseases 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 210000004220 fundus oculi Anatomy 0.000 description 2
- 210000002189 macula lutea Anatomy 0.000 description 2
- 210000004379 membrane Anatomy 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 210000001927 retinal artery Anatomy 0.000 description 2
- 230000002207 retinal effect Effects 0.000 description 2
- 210000001957 retinal vein Anatomy 0.000 description 2
- 201000004569 Blindness Diseases 0.000 description 1
- 206010061218 Inflammation Diseases 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000004054 inflammatory process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 208000013441 ocular lesion Diseases 0.000 description 1
- 238000011176 pooling Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
- 231100000216 vascular lesion Toxicity 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/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4053—Scaling of whole images or parts thereof, e.g. expanding or contracting based on super-resolution, i.e. the output image resolution being higher than the sensor resolution
-
- 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/12—Edge-based segmentation
-
- 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/155—Segmentation; Edge detection involving morphological operators
-
- 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/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
-
- 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/20—Special algorithmic details
- G06T2207/20036—Morphological image processing
- G06T2207/20044—Skeletonization; Medial axis transform
-
- 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/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- 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/20—Special algorithmic details
- G06T2207/20084—Artificial neural networks [ANN]
-
- 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/30096—Tumor; Lesion
-
- 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
-
- 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/30172—Centreline of tubular or elongated structure
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Geometry (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本公开描述了一种基于增强分辨率的眼底图像的血管管径的测量装置,包括用于获取眼底图像的获取单元;用于生成血管分割图像的分割单元;基于分辨率增强生成增强血管分割图像的增强单元;获取增强血管分割图像对应的连续血管骨架和测量像素点的管径测量方向的提取单元;基于增强血管分割图像、测量像素点、测量像素点的管径测量方向和预设精度,对增强血管分割图像进行插值操作以生成测量像素点对应的血管轮廓的生成单元;并且基于血管轮廓中的血管像素点的个数、预设倍数和预设精度计算血管管径的计算单元。由此,能够对血管管径进行自动超分辨率测量并具有较高的测量准确度。
Description
本申请是申请日为2020年09月04日、申请号为2020109242872、发明名称为眼底图像的血管管径的测量方法及测量装置的专利申请的分案申请。
技术领域
本公开涉及一种基于增强分辨率的眼底图像的血管管径的测量装置。
背景技术
眼底病变通常是指由炎症、肿瘤或血管病变等因素引起的眼部病变。眼底病变不仅会影响患者的视力,严重者还会导致失明,因此应定时进行眼底检查以便能较早地发现眼底病变并及时进行治疗。眼底病变可以通过多种的体征反应。其中眼底中的血管管径的变化是眼底病变的一个重要的体征。因此,在眼底检查中,通常通过对眼底中血管管径进行测量,从而获得血管管径的变化以对眼底病变进行识别。
传统的血管管径的测量方法(例如投影法、检眼镜测像法等)大多依靠专业的眼科医生借助测量工具并根据经验进行测量,虽然测量准确度较高,但耗时较长且无法批量自动处理,不利于眼底病变的识别。因此在眼底图像分析领域,准确度高且可自动测量血管管径的测量方法已经成为一项重要的研究内容。相比于传统的管径测量方法,眼底图像具有采集时间短,获取成本低等优点,因此近年来,越来越多的研究者利用眼底图像进行血管管径的自动测量。特别地,随着计算机图像处理技术的发展,基于眼底图像的血管管径的自动测量已被广泛研究。例如,Gregson等研究者使用血管灰度轮廓曲线矩形拟合的方法对血管管径进行自动测量;Peter等研究者通过搜索血管灰度轮廓曲线的二阶导数零点的方法对血管管径进行自动测量。然而,对血管管径进行自动测量的测量准确度还有待提高。
发明内容
本公开是有鉴于上述的状况而提出的,其目的在于提供一种能够对血管管径进行自动超分辨率测量并具有较高测量准确度的眼底图像的血管管径的测量方法及测量装置。
为此,本公开第一方面提供了一种眼底图像的血管管径的测量方法,其包括:获取包括血管的眼底图像;基于所述眼底图像生成血管分割图像;将所述血管分割图像按照预设倍数进行分辨率增强以生成增强血管分割图像;提取所述增强血管分割图像中的血管骨架并对所述血管骨架进行拟合以获取连续血管骨架和测量像素点的管径测量方向,所述测量像素点为所述连续血管骨架上的多个像素点,所述管径测量方向垂直于所述连续血管骨架在所述测量像素点的切线;基于所述增强血管分割图像、所述测量像素点、所述测量像素点的管径测量方向和预设精度生成所述测量像素点对应的血管轮廓;并且基于所述测量像素点对应的所述血管轮廓中的血管像素点的个数、所述预设倍数和所述预设精度计算所述测量像素点对应的血管管径。在本公开中,通过对血管分割图像进行增强以生成增强血管分割图像,提取增强血管分割图像的血管骨架并进行拟合以获取连续血管骨架和测量像素点的管径测量方向,基于增强血管分割图像、连续血管骨架上的测量像素点、测量像素点的管径测量方向和预设精度生成测量像素点对应的血管轮廓,通过测量像素点对应的血管轮廓中的血管像素点的个数、预设倍数和预设精度计算测量像素点对应的血管管径。在这种情况下,能够增加眼底图像中的血管的分辨率,进而可以利用更多的像素点对血管管径进行测量。由此,能够对血管管径进行自动超分辨率测量并提高测量准确度。
另外,在本公开第一方面所涉及的测量方法中,可选地,利用插值算法并基于所述增强血管分割图像、所述测量像素点、所述测量像素点的管径测量方向和所述预设精度,对所述增强血管分割图像进行插值操作以生成所述测量像素点对应的所述血管轮廓。由此,能够基于预设精度对增强血管图像进行插值操作以获取血管轮廓。
另外,在本公开第一方面所涉及的测量方法中,可选地,各个所述测量像素点对应的所述血管轮廓按照所述测量像素点在所述连续血管骨架上的排列顺序以中心共直线的方式并排以形成拉直血管图像。由此,能够方便地获取各个测量像素点对应的血管轮廓。
另外,在本公开第一方面所涉及的测量方法中,可选地,所述测量像素点对应的所述血管管径l为l=n×s/e,其中,n为所述测量像素点对应的所述血管轮廓中的血管像素点的个数,s为所述预设精度,e为所述预设倍数。由此,能够计算血管管径。
另外,在本公开第一方面所涉及的测量方法中,可选地,所述预设倍数为大于1的整数,所述预设精度为大于0且小于1的小数。由此,能够增加眼底图像血管的分辨率,进而能够提高后续对血管管径进行测量的测量准确度。
另外,在本公开第一方面所涉及的测量方法中,可选地,所述血管轮廓的宽度为最大血管管径的N倍,其中,N为2至5的整数。由此,能够使增强血管分割图像中最宽的血管完整地呈现在血管轮廓中。
另外,在本公开第一方面所涉及的测量方法中,可选地,所述插值算法为三次样条插值算法。由此,能够利用三次样条插值算法对增强血管分割图像进行插值操作。
另外,在本公开第一方面所涉及的测量方法中,可选地,所述测量像素点对应的所述血管轮廓中的血管像素点的个数n通过以下公式计算:n=card({p:p∈P,f(p)>T}),其中,p为所述测量像素点,f(p)为所述测量像素点p对应的像素值,T为预设的阈值参数,P为所述测量像素点p对应的所述血管轮廓中的像素点集合,card表示集合的基数。由此,能够计算测量像素点对应的血管轮廓中的血管像素点的个数。
本公开第二方面提供了一种眼底图像的血管管径的测量装置,其包括:获取单元,其用于获取包括血管的眼底图像;分割单元,其基于所述眼底图像生成血管分割图像;增强单元,其用于将所述血管分割图像按照预设倍数进行分辨率增强以生成增强血管分割图像;提取单元,其用于提取所述增强血管分割图像中的血管骨架并对所述血管骨架进行拟合以获取连续血管骨架和测量像素点的管径测量方向,所述测量像素点为所述连续血管骨架上的多个像素点,所述管径测量方向垂直于所述连续血管骨架在所述测量像素点的切线;生成单元,其基于所述增强血管分割图像、所述测量像素点、所述测量像素点的管径测量方向和预设精度生成所述测量像素点对应的血管轮廓;以及计算单元,其基于所述测量像素点对应的所述血管轮廓中的血管像素点的个数、所述预设倍数和所述预设精度计算所述测量像素点对应的血管管径。在本公开中,增强单元通过对分割单元生成的血管分割图像进行增强以生成增强血管分割图像,提取单元获取连续血管骨架和测量像素点的管径测量方向,生成单元基于增强血管分割图像、连续血管骨架上的测量像素点、测量像素点的管径测量方向和预设精度生成测量像素点对应的血管轮廓,计算单元通过测量像素点对应的血管轮廓中的血管像素点的个数、预设倍数和预设精度计算测量像素点对应的血管管径。在这种情况下,能够增加眼底图像中血管的分辨率,进而可以利用更多的像素点对血管管径进行测量。由此,能够对血管管径进行自动超分辨率测量并提高测量准确度。
另外,在本公开第二方面所涉及的测量装置中,可选地,利用插值算法并基于所述增强血管分割图像、所述测量像素点、所述测量像素点的管径测量方向和预设精度,对所述增强血管分割图像进行插值操作以生成所述测量像素点对应的所述血管轮廓。由此,能够基于预设精度对增强血管图像进行插值操作以获取血管轮廓。
另外,在本公开第二方面所涉及的测量装置中,可选地,所述测量像素点对应的所述血管管径l为l=n×s/e,其中,n为所述测量像素点对应的所述血管轮廓中的血管像素点的个数,s为所述预设精度,e为所述预设倍数。由此,能够计算血管管径。
根据本公开,提供一种能够对血管管径进行自动超分辨率测量并具有较高测量准确度的眼底图像的血管管径的测量方法及测量装置。
附图说明
现在将仅通过参考附图的例子进一步详细地解释本公开,其中:
图1是示出了本公开示例所涉及的眼底图像的血管管径的测量方法的应用场景的示意图。
图2是示出了本公开示例所涉及的眼底图像的血管管径的测量方法的流程图。
图3是示出了本公开示例所涉及的眼底图像的示意图。
图4是示出了本公开示例所涉及的基于人工神经网络的血管分割的训练方法的流程图。
图5(a)是示出了本公开示例所涉及的训练眼底图像的示意图。
图5(b)是示出了本公开示例所涉及的血管标注结果的示意图。
图6是示出了本公开示例所涉及的血管分割图像的示意图。
图7是示出了本公开示例所涉及的局部眼底图像的示意图。
图8是示出了本公开示例所涉及的局部血管分割图像的示意图。
图9是示出了本公开示例所涉及的局部血管骨架的示意图。
图10是示出了本公开示例所涉及的管径测量方向的示意图。
图11是示出了本公开示例所涉及的血管轮廓生成的流程图。
图12是示出了本公开示例所涉及的拉直血管图像的示意图。
图13是示出了本公开示例所涉及的眼底图像的血管管径的测量装置的***框图。
主要标号说明:
1…应用场景,11…操作员,12…终端,13…采集设备,14…人眼,15…服务器,P1…眼底图像,P2…训练眼底图像,P3…血管标注结果,P4…血管分割图像,P5…局部眼底图像,P6…局部血管分割图像,P7…拉直血管图像,L1…血管骨架,L2…连续血管骨架,L3…切线,L4…管径测量方向,2…测量装置,210…获取单元,220…分割单元,230…增强单元,240…提取单元,250…生成单元,260…计算单元,A…血管,B…血管,C…局部血管,C'…拉直局部血管,D1…测量像素点。
具体实施方式
以下,参考附图,详细地说明本公开的优选实施方式。在下面的说明中,对于相同的部件赋予相同的符号,省略重复的说明。另外,附图只是示意性的图,部件相互之间的尺寸的比例或者部件的形状等可以与实际的不同。
需要说明的是,本公开中的术语“包括”和“具有”以及它们的任何变形,例如所包括或所具有的一系列步骤或单元的过程、方法、***、产品或设备不必限于清楚地列出的那些步骤或单元,而是可以包括或具有没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
图1是示出了本公开示例所涉及的眼底图像的血管管径的测量方法的应用场景的示意图。
在一些示例中,本公开涉及的眼底图像的血管管径的测量方法(有时可以简称为测量方法)可以应用于如图1所示的应用场景1中。在应用场景1中,操作员11可以通过控制连接在终端12上的采集设备13,以采集人眼14的眼底图像。在采集设备13完成眼底图像采集后,终端12可以将眼底图像通过计算机网络提交至服务器15。服务器15通过执行本公开示例所涉及的测量方法可以获得眼底图像的血管管径并返回给终端12。在一些示例中,终端12可以显示血管管径。在另一些示例中,血管管径可以作为中间结果存储在终端12或者服务器15的存储器中。
在一些示例中,操作员11可以是具有对眼底图像的血管管径进行测量的专业知识的医师。在另一些示例中,操作员11可以是熟悉如何通过操作终端12自动对血管管径进行测量的普通人员。
在一些示例中,终端12可以包括但不限于是笔记本电脑、平板电脑或台式机等。在一些示例中,终端12可以是包括处理器、存储器、显示屏和采集设备13的对血管管径进行测量的专用设备。
在一些示例中,采集设备13可以包括但不限于是相机等。相机例如可以是彩色眼底相机。采集设备13可以通过串口与终端12相连接或者集成在终端12中。
在一些示例中,人眼14的眼底是指眼球内后部的组织,可以包括眼球的内膜、视网膜、黄斑和血管(视网膜动静脉)。在一些示例中,可以通过监测血管管径的变化对眼底病变进行识别。
在一些示例中,服务器15可以包括一个或多个处理器和一个或多个存储器。其中,处理器可以包括中央处理单元、图形处理单元以及能够处理数据的其它任何电子部件,能够执行计算机程序指令。存储器可以用于存储计算机程序指令。在一些示例中,测量方法可以以计算机程序指令的形式存储在存储器中并由服务器15执行。在一些示例中,服务器15也可以是云端服务器。
如上所述,服务器15可以通过执行本公开示例所涉及的测量方法获得眼底图像的血管管径。图2是示出了本公开示例所涉及的眼底图像的血管管径的测量方法的流程图。图3是示出了本公开示例所涉及的眼底图像的示意图。其中,眼底图像P1为由眼底相机拍摄的眼底图像。
在一些示例中,如图2所示,本公开涉及的眼底图像的血管管径的测量方法可以包括获取眼底图像(步骤S210),进行血管分割(步骤S220),增强分辨率(步骤S230),获取连续血管骨架和管径测量方向(步骤S240),生成血管轮廓(步骤S250)和计算血管管径(步骤S260)。在这种情况下,能够增加眼底图像中的血管的分辨率,进而可以利用更多的像素点对血管管径进行测量。由此,能够对血管管径进行自动超分辨率测量并提高测量准确度。
一般而言,对眼底图像中血管管径测量的宽度(也即像素点个数)为正整数,例如测量宽度可以为1、3、6或8等。而基于分辨率增强和预设精度的超分辨率技术可以进行血管管径的亚像素级别的测量,从而使测量宽度可以为小数,例如测量宽度可以为1.23、3.12、5.63或7.56等。由此,能够更加准确地对血管管径进行测量。
在步骤S210中,可以获取眼底图像。眼底图像可以包括血管。在一些示例中,眼底图像可以为彩色的眼底图像。彩色的眼底图像能够清晰的呈现眼球的内膜、视网膜、黄斑和血管(视网膜动静脉)等眼底信息。在另一些示例中,眼底图像可以为灰度图像。在一些示例中,眼底图像可以为由采集设备13采集的眼底图像。作为眼底图像的示例,例如图3示出了由眼底相机拍摄的眼底图像P1。其中,眼底图像P1包括血管A。在另一些示例中,眼底图像可以为预先存储在终端12或服务器15中的图像。
在一些示例中,在步骤S210中,可以对眼底图像进行预处理。一般而言,由于眼底图像可能存在图像格式、尺寸不同等问题,因此需要对眼底图像进行预处理,使眼底图像转换为固定标准形式的图像。固定标准形式可以指图像的格式相同,尺寸一致。例如,在一些示例中,眼底图像经过预处理后的尺寸可以统一宽度为512或1024像素。
在步骤S220中,可以对步骤S210获取的眼底图像进行血管分割。在一些示例中,基于眼底图像可以生成血管分割图像。如上所述,眼底可以包括眼球的内膜、视网膜、黄斑和血管(视网膜动静脉)。因此,眼底的眼底图像可以包括眼球的内膜、视网膜、黄斑和血管(视网膜动静脉)等图案,对眼底图像进行血管分割可以识别血管区域,从而去除血管外的其他结构的影响。
一般而言,血管分割方法可以分为无监督学习和监督学习方法。无监督学习方法可以先利用图像处理技术对眼底图像进行预处理,然后在经由预处理的眼底图像上采用阈值的方法进行分割。
在一些示例中,监督学习方法可以利用眼底图像及其血管标注结果,通过手动或者自动方式提取眼底图像的特征以对分类器进行训练,从而实现对眼底图像的血管分割。例如,可以通过训练人工神经网络获得的最优模型对眼底图像进行血管分割。
以下结合附图详细说明基于人工神经网络的血管分割的训练方法。图4是示出了本公开示例所涉及的基于人工神经网络的血管分割的训练方法的流程图。图5(a)是示出了本公开示例所涉及的训练眼底图像的示意图。图5(b)是示出了本公开示例所涉及的血管标注结果的示意图。
如图4所示,血管分割的训练方法可以包括获取训练眼底图像和血管标注结果(步骤S221),对训练眼底图像和血管标注结果进行预处理以获取预处理眼底图像和预处理血管标注结果(步骤S222)和基于预处理眼底图像和预处理血管标注结果对人工神经网络进行训练以获取最优模型(步骤S223)。在这种情况下,基于人工神经网络可以自动学习血管特征并输出血管分割结果。
在步骤S221中,可以获取训练眼底图像和血管标注结果。步骤S221中的训练眼底图像可以是拍摄眼底获得的眼底图像。作为训练眼底图像的示例,例如图5(a)示出训练眼底图像P2。其中,训练眼底图像P2包括血管B。在一些示例中,训练眼底图像可以为预先存储的多张眼底图像。在一些示例中,多张训练眼底图像可以组成训练数据集。
一般而言,对于监督学习的人工神经网络,需基于血管标注结果作为真值计算损失函数。因此,需要对训练眼底图像进行标注以获取训练眼底图像的血管标注结果。血管标注结果可以包括血管。在一些示例中,可以由经验丰富的医师利用标注工具对训练眼底图像进行标注以生成血管标注结果。作为血管标注结果的示例,例如图5(b)示出了对训练眼底图像P2进行标注后生成的血管标注结果P3。其中,血管标注结果P3包括血管B。
在步骤S222中,可以对训练眼底图像和血管标注结果进行预处理以获取预处理眼底图像和预处理血管标注结果。在一些示例中,可以对训练眼底图像进行预处理以生成预处理眼底图像。在一些示例中,对训练眼底图像进行预处理可以包括对训练眼底图像进行裁切、降噪、灰度化等。由此,能够突显训练眼底图像中的血管。在一些示例中,可以通过对训练眼底图像进行旋转,缩放、平移等操作。由此,能够增加人工神经网络训练的数据量。在一些示例中,可以对血管标注结果进行预处理以生成预处理血管标注结果。在一些示例中,对血管标注结果进行预处理与对训练眼底图像进行预处理的方式可以一样。
在步骤S223中,可以基于预处理眼底图像和预处理血管标注结果对人工神经网络进行训练以获取最优模型。在一些示例中,可以基于预处理眼底图像和预处理血管标注结果对人工神经网络进行训练以获取最优模型。在一些示例中,预处理血管标注结果可以作为真值计算损失函数,并通过损失函数不断地对人工神经网络进行优化直至损失函数的值(也即损失)收敛达到最优。
在一些示例中,人工神经网络可以是U-Net网络。具体地,U-net网络是一种图像语义分割的人工神经网络,包括特征提取和上采样部分。其中,特征提取部分可以包括多个编码层,第一个编码层的输入为训练眼底图像。在一些示例中,编码层可以包含一系列的卷积层、批归一化层、激活层和池化层。上采样部分可以包括多个解码层。最后一个解码层(即输出层)用于输出训练眼底图像中各个像素点属于血管的概率。在这种情况下,基于训练眼底图像中各个像素点属于血管的概率可以生成血管分割图像。
图6是示出了本公开示例所涉及的血管分割图像的示意图。
如上所述,人工神经网络可以基于训练眼底图像中各个像素点属于血管的概率生成血管分割图像。因此,经由对人工神经网络进行训练获得的最优模型也具有对眼底图像进行血管分割以生成血管分割图像的能力。在一些示例中,可以将步骤S210获取的眼底图像输入最优模型以对眼底图像进行血管分割,进而生成血管分割图像。在一些示例中,血管分割图像可以是二值图像,也即可以仅包括黑和白两种颜色。作为血管分割图像的示例,例如图6示出了对眼底图像P1进行血管分割生成的血管分割图像P4。其中,血管分割图像P4包括血管A。
如上所述,测量方法可以包括步骤S230。在步骤S230中,可以对步骤S220获取的血管分割图像进行分辨率增强。在一些示例,可以按照预设倍数将血管分割图像进行分辨率增强以生成增强血管分割图像。例如,在一些示例中,可以对分辨率大小为140×63的血管分割图像进行10倍分辨率增强,以生成分辨率大小为1400×630的增强血管分割图像。
另外,在一些示例中,可以利用线性插值方法对血管分割图像进行分辨率增强。但本公开的示例不限于此,在另一些示例中,可利用基于深度学习的图像超分辨率方法对血管分割图像进行分辨率增强。
另外,在一些示例中,预设倍数可以为大于1的整数。由此,能够增加眼底图像中血管的分辨率,进而能够提高后续对血管管径进行测量的测量准确度。在一些示例中,预设倍数可以为5至15。例如,预设倍数可以为5、10或15等。
图7示出了本公开示例所涉及的局部眼底图像的示意图。其中,局部眼底图像P5为眼底图像P1的局部图,局部眼底图像P5包含局部血管C。图8是示出了本公开示例所涉及的局部血管分割图像的示意图。其中,局部血管分割图像P6为血管分割图像P4中与局部眼底图像P5对应的局部图,局部血管分割图像P6包含局部血管C。图9是示出了本公开示例所涉及的局部血管骨架的示意图。
在步骤S240中,可以获取连续血管骨架和管径测量方向。
在一些示例中,在获取步骤S230生成的增强血管分割图像后,可以提取增强血管分割图像中的血管骨架。血管骨架可以为血管的中线。作为血管骨架的示例,例如,对图8所示的局部血管分割图像P6进行增强以提取局部血管C的血管骨架。图9示出了增强后的局部血管分割图像P6对应的局部血管C的血管骨架L1,从图9可以看出血管骨架L1是由离散的像素点组成。
在一些示例中,可以采用形态学细化算法对增强血管分割图像进行细化以提取血管骨架。也即将增强血管分割图像中血管宽度向血管中心方向细化成一个像素的宽度,形成血管骨架,并保持增强血管分割图像中的血管形状的基本拓扑结构不变。在一些示例中,形态学细化算法可以包括但不限于Hilditch细化算法、Pavlidis细化算法或Rosenfeld细化算法等。在一些示例中,对增强血管分割图像进行血管骨架提取前可以对增强血管分割图像进行中值滤波运算。由此,能够去除血管骨架末端可能出现的分叉。
如上所述,血管骨架可以由离散的像素点组成。因此,在一些示例中,可以对血管骨架进行拟合以获取连续血管骨架。在一些示例中,可以使用最小二乘三次样条插值算法对血管骨架进行拟合。由此,能够获得连续血管骨架和拟合方程。
图10是示出了本公开示例所涉及的管径测量方向的示意图。
另外,在一些示例中,在步骤S240中,可以获取测量像素点的管径测量方向。测量像素点可以为连续血管骨架上的多个像素点。在一些示例中,管径测量方向可以垂直于连续血管骨架在测量像素点的切线。在一些示例中,可以通过上述拟合方程的一阶导数获取测量像素点的切线。如图10所示,局部血管C的连续血管骨架L2在测量像素点D1的切线可以为切线L3。经过测量像素点D1并垂直于切线L3的直线可以为管径测量方向L4。在步骤S250中,可以生成血管轮廓。在一些示例中,可以基于步骤S230获得的增强血管分割图像、步骤S240获得的测量像素点、步骤S240获得的测量像素点的管径测量方向和预设精度生成测量像素点对应的血管轮廓。
在一些示例中,可以利用插值算法并基于增强血管分割图像、测量像素点、测量像素点的管径测量方向和预设精度,对增强血管分割图像进行插值操作以生成测量像素点对应的血管轮廓。由此,能够基于预设精度对增强血管图像进行插值操作以获取血管轮廓。在一些示例中,插值算法可以为三次样条插值算法。由此,能够利用三次样条插值算法对增强血管分割图像进行插值操作。
另外,在一些示例中,预设精度可以为大于0且小于1的小数。由此,能够增加眼底图像中的血管的分辨率,进而能够提高后续对血管管径进行测量的测量准确度。在一些示例中,预设精度可以为0.01至0.10。例如,预设精度可以为0.01、0.05或0.10等。
以下,结合附图详细说明血管轮廓的生成过程。图11是示出了本公开示例所涉及的血管轮廓生成的流程图。
如图11所示,步骤S250中血管轮廓的生成过程可以包括获取血管轮廓的宽度(步骤S251),获取插值采样间隔(步骤S252),生成插值点(步骤S253),对血管分割图像进行插值操作以确定血管轮廓的像素值(步骤S254)和基于血管轮廓的像素值生成血管轮廓(步骤S255)。
在步骤S251中,可以获取血管轮廓的宽度。在一些示例中,血管轮廓的宽度可以为最大血管管径的N倍。其中,N可以为2到5之间的整数。在一些示例中,血管轮廓的宽度可以为最大血管管径的2倍。在一些示例中,最大血管管径可以是预估的最大血管管径。由此,能够使增强血管分割图像中最宽的血管完整地呈现在血管轮廓中。
在一些示例中,预估的最大血管管径可以为连续血管骨架中的像素点与其最近的非血管像素的距离中的最大距离的2倍。其中距离可以是欧几里德距离。具体地,计算连续血管骨架中的每个像素点与其最近的非血管像素的欧几里德距离,选择最大欧几里德距离的2倍作为预估的最大血管管径。由此,可以粗略估算出最大血管管径。但本公开的示例不限于此,在另一些示例中,可以使用其他方式估算最大血管管径。
在步骤S252中,可以获取插值采样间隔。在一些示例中,基于血管轮廓的宽度、预设精度可以获取插值采样间隔。具体地,假设血管轮廓的宽度可以用wi表示,预设精度可以用s表示。则插值采样间隔的值可以为在-(wi-1)/2至(wi-1)/2之间,以预设精度s为步长递增的值。例如插值采样间隔的值可以为-(wi-1)/2、(wi-1)/2+s、(wi-1)/2+2×s或(wi-1)/2等。
在步骤S253中,可以生成插值点。在一些示例中,可以基于步骤S252获得的插值采样间隔、管径测量方向和连续血管骨架生成插值点。具体地,假设插值采样间隔用inc表示、管径测量方向可以用(dx,dy)表示,连续血管骨架可以用(x,y)表示。则插值点可以表示为(x+dx×inc,y+dy×inc)。在这种情况下,生成的插值点沿着管径测量方向分布。由此,能够增加管径测量方向附近的像素点个数,进而能够提高血管管径的测量准确度。
在步骤S254中,可以对血管分割图像进行插值操作以确定血管轮廓的像素值。在一些示例中,可以使用插值算法,基于步骤S253获得的插值点对血管分割图像进行插值操作以确定血管轮廓的像素值。在一些示例中,插值算法可以为三次样条插值算法。由此,能够利用三次样条插值算法对增强血管分割图像进行插值操作。
图12是示出了本公开示例所涉及的拉直血管图像的示意图。
在步骤S255中,可以基于血管轮廓的像素值生成血管轮廓。在一些示例中,可以基于血管轮廓的像素值输出各个测量像素点对应的血管轮廓。在一些示例中,各个测量像素点对应的血管轮廓可以按照测量像素点在连续血管骨架上的排列顺序以中心共直线的方式并排以形成拉直血管图像。作为拉直血管图像的示例,图12示出了拉直血管图像P7。其中,拉直血管图像P7包括局部血管C对应的拉直局部血管C'。由此,能够方便地获取各个测量像素点对应的血管轮廓。
如上所述,测量方法包括步骤S260。在步骤S260,可以计算血管管径。在一些示例中,可以基于测量像素点对应的血管轮廓中的血管像素点的个数、预设倍数和预设精度计算测量像素点对应的血管管径。其中,测量像素点对应的血管轮廓可以由步骤S250获得。
在一些示例中,测量像素点对应的血管轮廓中的血管像素点的个数n可以通过以下公式计算:
n=card({p:p∈P,f(p)>T}),
其中,p为测量像素点,f(p)为测量像素点p对应的像素值,T为预设的阈值参数,P为测量像素点p对应的血管轮廓中的像素点集合,card表示集合的基数。由此,能够计算测量像素点对应的血管轮廓中的血管像素点的个数。在一些示例中,T可以为0.9。由此,能够区分测量像素点对应的血管轮廓中的血管像素点和非血管像素点。
在一些示例中,测量像素点对应的血管管径l可以为:
l=n×s/e,
其中,n为测量像素点对应的血管轮廓中的血管像素点的个数,s为预设精度,e为预设倍数。由此,能够计算血管管径。
如上所述,预设倍数可以为大于1的整数。预设精度可以为大于0且小于1的小数。在一些示例中,可以利用不同的预设倍数e和预设精度s对同一段血管的血管管径进行测量以获取准确度较高的血管管径的测量结果。
以下以图7所示的局部眼底图像P5中的局部血管C为例,将基于不同预设倍数e和预设精度s的局部血管C的血管管径测量结果与基于人工标注的局部血管C的血管管径测量结果进行对比,对比结果可以参见表1。其中,均值为局部血管C的血管管径的平均值,标准差为局部血管C的血管管径的标准差。
表1血管管径测量结果对比的部分结果
测量参数 | 均值 | 标准差 |
人工标注 | 6.4094 | 0.421 |
e=1,s=1 | 5.9365 | 0.7319 |
e=5,s=0.5 | 6.4457 | 0.4737 |
e=5,s=0.1 | 6.435 | 0.4723 |
e=5,s=0.05 | 6.4341 | 0.4707 |
e=5,s=0.01 | 6.4345 | 0.471 |
e=10,s=0.5 | 6.5142 | 0.4799 |
e=10,s=0.1 | 6.5078 | 0.4791 |
e=10,s=0.05 | 6.5069 | 0.4787 |
e=10,s=0.01 | 6.507 | 0.4788 |
e=15,s=0.5 | 6.5341 | 0.4757 |
e=15,s=0.1 | 6.5293 | 0.4757 |
e=15,s=0.05 | 6.5288 | 0.4754 |
e=15,s=0.01 | 6.5289 | 0.4755 |
从血管管径测量结果对比表1中可以看,针对局部血管C,人工标注的血管管径均值为6.4094,标准差为0.421。在不使用超分辨技术(即e=1,s=1)时,本发明测量出的血管管径的均值为5.9365,标准差为0.7319,与人工标注结果有较大差异。而在使用超分辨技术(即e∈{5,10,15}且s∈{0.01,0.05,0.1,0.5})时,管径测量的均值和标准差都接近于人工标注结果。特别地,在e=5,s=0.05时,血管管径的均值为6.4341,标准差为0.4707,与人工标注的血管管径相当。由此可见,使用超分辨率技术的测量结果与人工标注的结果更接近,说明超分辨率技术有效地提高了血管管径测量的准确度。
在一些示例中,可以将血管管径转成实际尺寸大小。具体地,可以在采集眼底图像时输入一幅包括标尺的图像。标尺上带有刻度。在一些示例中,可以基于图像处理算法对包括标尺的图像进行边缘检测以获得标尺上刻度较清晰的区域,进而可以获取该区域对应的刻度范围,例如刻度范围可以为2um(微米)到100um(微米)。基于该刻度范围和眼底图像中该刻度范围之间的像素点个数可以获得像素点距离。在一些示例中,可以利用该刻度范围获得实际尺寸大小,例如98um(微米)除以眼底图像中该刻度范围之间的像素点个数以获得像素点距离。在这种情况下,通过将血管管径乘以像素点距离可以将血管管径转成实际尺寸大小。
但本公开的示例不限于此,在另一些示例中,基于DICOM(Digital Imaging andCommunications in Medicine,医学数字成像与通讯)格式的眼底图像可以获取像素间隔(眼底图像中两个像素点之间距离代表的实际尺寸大小)。在这种情况下,通过将血管管径乘以像素间隔可以将血管管径转成实际尺寸大小。在另一些示例中,可以不将血管管径转成实际尺寸大小,例如可以将血管管径作为中间结果。
以下,结合图13详细描述本公开所涉及的眼底图像的血管管径的测量装置。本公开中的眼底图像的血管管径的测量装置可以简称为“测量装置”。测量装置用于实施上述测量方法。图13是示出了本公开示例所涉及的眼底图像的血管管径的测量装置的***框图。
在一些示例中,如图13所示,测量装置2可以包括获取单元210、分割单元220、增强单元230、提取单元240、生成单元250和计算单元260。获取单元210可以用于获取眼底图像。分割单元220可以用于对眼底图像进行血管分割。增强单元230可以用于对血管分割图像进行分辨率增强。提取单元240可以用于获取连续血管骨架和管径测量方向。生成单元250可以用于生成血管轮廓。计算单元260可以用于计算血管管径。在这种情况下,能够增加眼底图像中的血管的分辨率,进而可以利用更多的像素点对血管管径进行测量。由此,能够对血管管径进行自动超分辨率测量并提高测量准确度。
在一些示例中,获取单元210可以用于获取眼底图像。眼底图像可以包括血管。在一些示例中,眼底图像可以为彩色的眼底图像。在一些示例中,可以对眼底图像进行预处理。具体内容可以参见步骤S210,此处不再赘述。
在一些示例中,分割单元220可以用于对获取单元210获取的眼底图像进行血管分割。在一些示例中,基于眼底图像可以生成血管分割图像。在一些示例中,可以通过人工神经网络对眼底图像进行血管分割。在这种情况下,基于人工神经网络可以自动学习血管特征并输出血管分割结果。具体内容可以参见步骤S220,此处不再赘述。
在一些示例中,增强单元230可以用于对血管分割图像进行分辨率增强。在一些示例,可以将血管分割图像按照预设倍数进行分辨率增强以生成增强血管分割图像。在一些示例中,可以利用线性插值方法对血管分割图像进行分辨率增强。在一些示例中,预设倍数可以设为大于1的整数。由此,能够增加眼底图像中的血管的分辨率,进而能够提高后续对血管管径进行测量的测量准确度。具体内容可以参见步骤S230,此处不再赘述。
在一些示例中,提取单元240可以用于获取连续血管骨架和管径测量方向。在一些示例中,在增强单元230生成的增强血管分割图像后,可以提取增强血管分割图像中的血管骨架。血管骨架可以为血管的中线。在一些示例中,可以采用形态学细化算法对增强血管分割图像进行细化以提取血管骨架。在一些示例中,对增强血管分割图像进行血管骨架提取前可以对增强血管分割图像进行中值滤波运算。由此,能够去除血管骨架末端可能出现的分叉。在一些示例中,可以对血管骨架进行拟合以获取连续血管骨架。在一些示例中,可以使用最小二乘三次样条插值算法对血管骨架进行拟合。由此,能够获得连续血管骨架和拟合方程。在一些示例中,提取单元240可以用于获取测量像素点的管径测量方向。测量像素点可以为连续血管骨架上的多个像素点。在一些示例中,管径测量方向可以垂直于连续血管骨架在测量像素点的切线。具体内容可以参见步骤S240,此处不再赘述。
在一些示例中,生成单元250可以用于生成血管轮廓。在一些示例中,可以基于增强单元S230获得的增强血管分割图像、提取单元240获得的测量像素点、提取单元240获得的测量像素点的管径测量方向和预设精度生成测量像素点对应的血管轮廓。在一些示例中,可以利用插值算法并基于增强血管分割图像、测量像素点、测量像素点的管径测量方向和预设精度,对增强血管分割图像进行插值操作以生成测量像素点对应的血管轮廓。由此,能够基于预设精度对增强血管分割图像进行插值操作。在一些示例中,插值算法可以为三次样条插值算法。由此,能够利用三次样条插值算法对增强血管分割图像进行插值操作。具体内容可以参见步骤S250,此处不再赘述。
在一些示例中,在生成单元250中,预设精度可以为大于0且小于1的小数。由此,能够增加眼底图像中的血管的分辨率,进而能够提高后续对血管管径进行测量的测量准确度。在一些示例中,血管轮廓的宽度可以为最大血管管径的N倍,其中,N可以为2至5的整数。在一些示例中,各个测量像素点对应的血管轮廓按照测量像素点在连续血管骨架上的排列顺序可以以中心共直线的方式并排以形成拉直血管图像。具体内容可以参见步骤S250,此处不再赘述。
在一些示例中,计算单元260可以用于计算血管管径。在一些示例中,可以基于测量像素点对应的血管轮廓中的血管像素点的个数、预设倍数和预设精度计算测量像素点对应的血管管径。在一些示例中,测量像素点对应的所述血管轮廓中的血管像素点的个数n可以通过以下公式计算:
n=card({p:p∈P,f(p)>T}),
其中,p为测量像素点,f(p)为测量像素点p对应的像素值,T为预设的阈值参数,P为测量像素点p对应的血管轮廓中的像素点集合,card表示集合的基数。由此,能够计算测量像素点对应的血管轮廓中的血管像素点的个数。在一些示例中,T可以为0.9。由此,能够区分测量像素点对应的血管轮廓中的血管像素点和非血管像素点。在一些示例中,测量像素点对应的血管管径l可以为:
l=n×s/e,
其中,n为测量像素点对应的血管轮廓中的血管像素点的个数,s为预设精度,e为预设倍数。由此,能够计算血管管径。具体内容可以参见步骤S260,此处不再赘述。
在本公开中,通过对血管分割图像进行增强以生成增强血管分割图像,提取增强血管分割图像的血管骨架并进行拟合以获取连续血管骨架和管径测量方向,基于增强血管分割图像、连续血管骨架上的测量像素点、测量像素点的管径测量方向和预设精度生成测量像素点对应的血管轮廓,通过测量像素点对应的血管轮廓中的血管像素点的个数、预设倍数和预设精度计算测量像素点对应的血管管径。在这种情况下,能够增加眼底图像中的血管的分辨率,进而可以利用更多的像素点对血管管径进行测量。由此,能够提高血管管径的测量准确度。
虽然以上结合附图和示例对本公开进行了具体说明,但是可以理解,上述说明不以任何形式限制本公开。本领域技术人员在不偏离本公开的实质精神和范围的情况下可以根据需要对本公开进行变形和变化,这些变形和变化均落入本公开的范围内。
Claims (10)
1.一种基于增强分辨率的眼底图像的血管管径的测量装置,其特征在于,包括获取单元、分割单元、增强单元、提取单元、生成单元和计算单元;所述获取单元用于获取包括血管的眼底图像;所述分割单元基于所述眼底图像生成血管分割图像;所述增强单元用于将所述血管分割图像按照预设倍数进行分辨率增强以生成增强血管分割图像;所述提取单元用于获取所述增强血管分割图像对应的连续血管骨架和测量像素点的管径测量方向,所述测量像素点为所述连续血管骨架上的多个像素点,所述管径测量方向垂直于所述连续血管骨架在所述测量像素点的切线;所述生成单元基于所述增强血管分割图像、所述测量像素点、所述测量像素点的管径测量方向和预设精度,对所述增强血管分割图像进行插值操作以生成所述测量像素点对应的血管轮廓,所述预设精度为插值采样的步长;以及所述计算单元基于所述测量像素点对应的血管轮廓中的血管像素点的个数、所述预设倍数和所述预设精度计算所述测量像素点对应的血管管径。
2.根据权利要求1所述的测量装置,其特征在于:
通过将所述血管管径乘以像素点距离将所述血管管径转成实际尺寸大小,其中,所述像素点距离为所述眼底图像中两个像素点之间的距离代表的实际尺寸大小。
3.根据权利要求2所述的测量装置,其特征在于:
在采集所述眼底图像时输入一幅包括带有刻度的标尺的图像,对所述包括带有刻度的标尺的图像进行边缘检测以获得标尺上刻度清晰的区域对应的刻度范围,基于该刻度范围和所述眼底图像中该刻度范围之间的像素点个数获得所述像素点距离。
4.根据权利要求1所述的测量装置,其特征在于:
所述测量像素点对应的所述血管管径l为
l=n×s/e,
其中,n为所述测量像素点对应的血管轮廓中的血管像素点的个数,s为所述预设精度,e为所述预设倍数。
5.根据权利要求4所述的测量装置,其特征在于:
所述测量像素点对应的血管轮廓中的血管像素点的个数n通过以下公式计算:
n=card({p:p∈P,f(p)>T}),
其中,p为所述测量像素点,f(p)为所述测量像素点p对应的像素值,T为预设的阈值参数,P为所述测量像素点p对应的血管轮廓中的像素点集合,card表示集合的基数。
6.根据权利要求5所述的测量装置,其特征在于:
T为0.9。
7.根据权利要求1所述的测量装置,其特征在于:
所述预设倍数为大于1的整数,所述预设精度为大于0且小于1的小数。
8.根据权利要求1所述的测量装置,其特征在于:
所述血管轮廓的宽度为最大血管管径的N倍,其中,N为2至5的整数。
9.根据权利要求1所述的测量装置,其特征在于,
提取所述增强血管分割图像中的血管骨架并对所述血管骨架进行拟合以获取所述连续血管骨架和所述测量像素点的管径测量方向。
10.根据权利要求1所述的测量装置,其特征在于:
训练人工神经网络以获得最优模型并利用所述最优模型对所述眼底图像进行血管分割以生成所述血管分割图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110951405.3A CN113643354B (zh) | 2020-09-04 | 2020-09-04 | 基于增强分辨率的眼底图像的血管管径的测量装置 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010924287.2A CN112288794B (zh) | 2020-09-04 | 2020-09-04 | 眼底图像的血管管径的测量方法及测量装置 |
CN202110951405.3A CN113643354B (zh) | 2020-09-04 | 2020-09-04 | 基于增强分辨率的眼底图像的血管管径的测量装置 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010924287.2A Division CN112288794B (zh) | 2020-09-04 | 2020-09-04 | 眼底图像的血管管径的测量方法及测量装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113643354A true CN113643354A (zh) | 2021-11-12 |
CN113643354B CN113643354B (zh) | 2023-10-13 |
Family
ID=74419886
Family Applications (3)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010924287.2A Active CN112288794B (zh) | 2020-09-04 | 2020-09-04 | 眼底图像的血管管径的测量方法及测量装置 |
CN202110950244.6A Active CN113643353B (zh) | 2020-09-04 | 2020-09-04 | 眼底图像的血管管径的增强分辨率的测量方法 |
CN202110951405.3A Active CN113643354B (zh) | 2020-09-04 | 2020-09-04 | 基于增强分辨率的眼底图像的血管管径的测量装置 |
Family Applications Before (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010924287.2A Active CN112288794B (zh) | 2020-09-04 | 2020-09-04 | 眼底图像的血管管径的测量方法及测量装置 |
CN202110950244.6A Active CN113643353B (zh) | 2020-09-04 | 2020-09-04 | 眼底图像的血管管径的增强分辨率的测量方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US20240005545A1 (zh) |
EP (1) | EP4187489A1 (zh) |
CN (3) | CN112288794B (zh) |
WO (1) | WO2022048171A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115100222A (zh) * | 2022-08-24 | 2022-09-23 | 首都医科大学附属北京朝阳医院 | 分离动静脉血管的图像处理方法、装置、存储介质及终端 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112288794B (zh) * | 2020-09-04 | 2021-09-07 | 深圳硅基智能科技有限公司 | 眼底图像的血管管径的测量方法及测量装置 |
CN117576096A (zh) * | 2024-01-16 | 2024-02-20 | 成都泰盟软件有限公司 | 一种基于图像识别的自动测量血管管径的方法及装置 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05264232A (ja) * | 1991-12-09 | 1993-10-12 | Kokuritsu Jiyunkankibiyou Center Souchiyou | 造影血管高精度管径計測装置 |
CN101984916A (zh) * | 2010-11-17 | 2011-03-16 | 哈尔滨工程大学 | 基于数字图像处理技术的血管管径测量方法 |
JP2014193225A (ja) * | 2013-03-29 | 2014-10-09 | Nidek Co Ltd | 眼底画像処理装置および眼底画像処理プログラム |
US20160287071A1 (en) * | 2015-04-03 | 2016-10-06 | Oregon Health & Science University | Methods of measuring total retinal blood flow using en face doppler oct |
CN106934816A (zh) * | 2017-03-23 | 2017-07-07 | 中南大学 | 一种基于elm的眼底图像视网膜血管分割方法 |
JP2018114031A (ja) * | 2017-01-16 | 2018-07-26 | 大日本印刷株式会社 | 眼底画像処理装置 |
CN108537789A (zh) * | 2018-04-13 | 2018-09-14 | 深圳大学 | 管腔识别方法、装置及计算机存储介质 |
CN109166124A (zh) * | 2018-11-20 | 2019-01-08 | 中南大学 | 一种基于连通区域的视网膜血管形态量化方法 |
CN109829942A (zh) * | 2019-02-21 | 2019-05-31 | 韶关学院 | 一种眼底图像视网膜血管管径自动量化方法 |
CN110197493A (zh) * | 2019-05-24 | 2019-09-03 | 清华大学深圳研究生院 | 眼底图像血管分割方法 |
CN111340789A (zh) * | 2020-02-29 | 2020-06-26 | 平安科技(深圳)有限公司 | 眼底视网膜血管识别及量化方法、装置、设备及存储介质 |
Family Cites Families (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101672904B (zh) * | 2009-09-14 | 2012-09-05 | 东南大学 | 用于便携式核磁共振血管检测装置的探头 |
US8355544B2 (en) * | 2011-02-01 | 2013-01-15 | Universidade Da Coruna-Otri | Method, apparatus, and system for automatic retinal image analysis |
WO2015003225A1 (en) * | 2013-07-10 | 2015-01-15 | Commonwealth Scientific And Industrial Research Organisation | Quantifying a blood vessel reflection parameter of the retina |
US9968250B2 (en) * | 2014-02-24 | 2018-05-15 | Image Technologies Corporation | Autofluorescence imaging of macular pigment: image quality criteria and corrections |
JP2015202236A (ja) * | 2014-04-15 | 2015-11-16 | キヤノン株式会社 | 眼底画像解析装置及び解析方法 |
WO2016093140A1 (ja) * | 2014-12-09 | 2016-06-16 | 雫石 誠 | 撮像装置 |
US20170100029A1 (en) * | 2015-10-13 | 2017-04-13 | The University Of North Carolina At Chapel Hill | Compositions and Methods for Analyzing Collateral Density |
CN107292835B (zh) * | 2017-05-31 | 2020-03-13 | 瑞达昇医疗科技(大连)有限公司 | 一种眼底图像视网膜血管自动矢量化的方法及装置 |
WO2018237011A1 (en) * | 2017-06-20 | 2018-12-27 | University Of Louisville Research Foundation, Inc. | SEGMENTATION OF RETINIAN BLOOD VESSELS IN OPTICAL COHERENCE TOMOGRAPHY ANGIOGRAPHY IMAGES |
CN108122236B (zh) * | 2017-12-18 | 2020-07-31 | 上海交通大学 | 基于距离调制损失的迭代式眼底图像血管分割方法 |
CN110310323A (zh) * | 2018-03-20 | 2019-10-08 | 天津工业大学 | 基于Hessian矩阵和二维高斯拟合的视网膜血管管径测量方法 |
CN108961334B (zh) * | 2018-06-26 | 2020-05-08 | 电子科技大学 | 一种基于图像配准的视网膜血管壁厚度测量方法 |
CN109512450A (zh) * | 2018-10-18 | 2019-03-26 | 深圳市孙逸仙心血管医院(深圳市心血管病研究所) | 测量血管血流速度的方法 |
CN109658393B (zh) * | 2018-12-06 | 2022-11-22 | 代黎明 | 眼底图像拼接方法及*** |
DE102018222606A1 (de) * | 2018-12-20 | 2020-06-25 | Siemens Healthcare Gmbh | Verfahren und Vorrichtung zur Detektion eines anatomischen Merkmals eines Blutgefäßabschnittes |
CN109886938B (zh) * | 2019-01-29 | 2023-07-18 | 深圳市科曼医疗设备有限公司 | 一种超声图像血管直径自动测量方法 |
CN109993731A (zh) * | 2019-03-22 | 2019-07-09 | 依未科技(北京)有限公司 | 一种眼底病变分析方法及装置 |
CN110348541B (zh) * | 2019-05-10 | 2021-12-10 | 腾讯医疗健康(深圳)有限公司 | 眼底血管图像分类方法、装置、设备及存储介质 |
CN110163872A (zh) * | 2019-05-14 | 2019-08-23 | 西南科技大学 | 一种hrmr图像分割与三维重建的方法及电子设备 |
CN110652312B (zh) * | 2019-07-19 | 2023-03-14 | 慧影医疗科技(北京)股份有限公司 | 一种血管cta智能分析***及应用 |
CN111062942B (zh) * | 2020-03-16 | 2020-06-12 | 南京景三医疗科技有限公司 | 血管分叉的检测方法、装置及医疗设备 |
CN112001928B (zh) * | 2020-07-16 | 2023-12-15 | 北京化工大学 | 一种视网膜血管分割方法及*** |
CN112288794B (zh) * | 2020-09-04 | 2021-09-07 | 深圳硅基智能科技有限公司 | 眼底图像的血管管径的测量方法及测量装置 |
CN112716446B (zh) * | 2020-12-28 | 2023-03-24 | 深圳硅基智能科技有限公司 | 高血压性视网膜病变的病变特征的测量方法及测量*** |
-
2020
- 2020-09-04 CN CN202010924287.2A patent/CN112288794B/zh active Active
- 2020-09-04 CN CN202110950244.6A patent/CN113643353B/zh active Active
- 2020-09-04 CN CN202110951405.3A patent/CN113643354B/zh active Active
-
2021
- 2021-04-29 WO PCT/CN2021/091223 patent/WO2022048171A1/zh active Application Filing
- 2021-04-29 US US18/043,072 patent/US20240005545A1/en active Pending
- 2021-04-29 EP EP21863241.2A patent/EP4187489A1/en active Pending
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH05264232A (ja) * | 1991-12-09 | 1993-10-12 | Kokuritsu Jiyunkankibiyou Center Souchiyou | 造影血管高精度管径計測装置 |
CN101984916A (zh) * | 2010-11-17 | 2011-03-16 | 哈尔滨工程大学 | 基于数字图像处理技术的血管管径测量方法 |
JP2014193225A (ja) * | 2013-03-29 | 2014-10-09 | Nidek Co Ltd | 眼底画像処理装置および眼底画像処理プログラム |
US20160287071A1 (en) * | 2015-04-03 | 2016-10-06 | Oregon Health & Science University | Methods of measuring total retinal blood flow using en face doppler oct |
JP2018114031A (ja) * | 2017-01-16 | 2018-07-26 | 大日本印刷株式会社 | 眼底画像処理装置 |
CN106934816A (zh) * | 2017-03-23 | 2017-07-07 | 中南大学 | 一种基于elm的眼底图像视网膜血管分割方法 |
CN108537789A (zh) * | 2018-04-13 | 2018-09-14 | 深圳大学 | 管腔识别方法、装置及计算机存储介质 |
CN109166124A (zh) * | 2018-11-20 | 2019-01-08 | 中南大学 | 一种基于连通区域的视网膜血管形态量化方法 |
CN109829942A (zh) * | 2019-02-21 | 2019-05-31 | 韶关学院 | 一种眼底图像视网膜血管管径自动量化方法 |
CN110197493A (zh) * | 2019-05-24 | 2019-09-03 | 清华大学深圳研究生院 | 眼底图像血管分割方法 |
CN111340789A (zh) * | 2020-02-29 | 2020-06-26 | 平安科技(深圳)有限公司 | 眼底视网膜血管识别及量化方法、装置、设备及存储介质 |
Non-Patent Citations (4)
Title |
---|
LUMINTA MORARU ET AL.: "Assessment of vessel diameters for MR brain angiography processed images", 《AIP CONFERENCE PROCEEDINGS》, pages 1 - 7 * |
MAHDIEH REZAEIAN ET AL.: "A novel method for retinal vessel segmentation and diameter measurement using high speed video", 《2019 41ST ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SICUETY》, pages 2781 - 2784 * |
杨嵩: "彩色眼底图像血管参数测量方法研究", 《中国优秀硕士学位论文全文数据库信息科技辑》, pages 138 - 1181 * |
薛岚燕;余轮;林嘉雯;曹新容;郑绍华;: "视网膜血管直径交互式测量方法", 福州大学学报(自然科学版), no. 03, pages 45 - 53 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115100222A (zh) * | 2022-08-24 | 2022-09-23 | 首都医科大学附属北京朝阳医院 | 分离动静脉血管的图像处理方法、装置、存储介质及终端 |
CN115100222B (zh) * | 2022-08-24 | 2022-12-09 | 首都医科大学附属北京朝阳医院 | 分离动静脉血管的图像处理方法、装置、存储介质及终端 |
Also Published As
Publication number | Publication date |
---|---|
US20240005545A1 (en) | 2024-01-04 |
WO2022048171A1 (zh) | 2022-03-10 |
CN113643353B (zh) | 2024-02-06 |
CN112288794B (zh) | 2021-09-07 |
EP4187489A1 (en) | 2023-05-31 |
CN113643353A (zh) | 2021-11-12 |
CN112288794A (zh) | 2021-01-29 |
CN113643354B (zh) | 2023-10-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112288794B (zh) | 眼底图像的血管管径的测量方法及测量装置 | |
CN112716446B (zh) | 高血压性视网膜病变的病变特征的测量方法及测量*** | |
US8855388B2 (en) | Microcalcification detection classification in radiographic images | |
Chakravarty et al. | Joint optic disc and cup boundary extraction from monocular fundus images | |
Luo et al. | Micro-vessel image segmentation based on the AD-UNet model | |
CN112837805B (zh) | 基于深度学习的眼睑拓扑形态特征的提取方法 | |
CN112164043A (zh) | 一种用于多眼底图像拼接的方法及*** | |
CN112750531A (zh) | 一种中医自动化望诊***、方法、设备和介质 | |
CN112562058B (zh) | 一种基于迁移学习的颅内血管模拟三维模型快速建立方法 | |
Maaliw et al. | An enhanced segmentation and deep learning architecture for early diabetic retinopathy detection | |
Rodrigues et al. | Retinal vessel segmentation using parallel grayscale skeletonization algorithm and mathematical morphology | |
CN108665474B (zh) | 一种基于b-cosfire的眼底图像视网膜血管分割方法 | |
CN115393239A (zh) | 一种多模态眼底图像配准和融合方法及*** | |
JPWO2019073962A1 (ja) | 画像処理装置及びプログラム | |
CN113516644B (zh) | 一种血管内光学相干层析成像分支血管识别方法 | |
Pavlov et al. | Methods and computer tools for identifying diabetes-induced fundus pathology | |
CN111583248B (zh) | 一种基于眼部超声图像的处理方法 | |
Liu et al. | Retinal vessel segmentation using densely connected convolution neural network with colorful fundus images | |
CN109816665B (zh) | 一种光学相干断层扫描图像的快速分割方法及装置 | |
CN115690556B (zh) | 一种基于多模态影像学特征的图像识别方法及*** | |
CN111292285B (zh) | 一种基于朴素贝叶斯与支持向量机的糖网病自动筛查方法 | |
Rashid et al. | Segmenting melanoma lesion using single shot detector (SSD) and level set segmentation technique | |
Xiang et al. | Segmentation of retinal blood vessels based on divergence and bot-hat transform | |
Liang et al. | Location of optic disk in the fundus image based on visual attention | |
Ponmalar et al. | Efficient Diabetic Retinopathy Detection using Machine Learning Techniques |
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 |