CN104573712B - 基于眼底图像的动静脉视网膜血管分类方法 - Google Patents
基于眼底图像的动静脉视网膜血管分类方法 Download PDFInfo
- Publication number
- CN104573712B CN104573712B CN201410850207.8A CN201410850207A CN104573712B CN 104573712 B CN104573712 B CN 104573712B CN 201410850207 A CN201410850207 A CN 201410850207A CN 104573712 B CN104573712 B CN 104573712B
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- main
- fundus image
- eye fundus
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 210000001210 retinal vessel Anatomy 0.000 title claims abstract description 17
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 241
- 210000003733 optic disk Anatomy 0.000 claims abstract description 40
- 230000002792 vascular Effects 0.000 claims description 43
- 230000011218 segmentation Effects 0.000 claims description 41
- 239000013598 vector Substances 0.000 claims description 27
- 210000003462 vein Anatomy 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 238000004040 coloring Methods 0.000 claims description 7
- 239000012141 concentrate Substances 0.000 claims description 6
- 239000000284 extract Substances 0.000 claims description 5
- 210000001367 artery Anatomy 0.000 claims description 4
- 239000012634 fragment Substances 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 2
- 238000009792 diffusion process Methods 0.000 abstract description 2
- 239000008280 blood Substances 0.000 description 11
- 210000004369 blood Anatomy 0.000 description 11
- 125000004122 cyclic group Chemical group 0.000 description 11
- 230000003902 lesion Effects 0.000 description 9
- 208000032843 Hemorrhage Diseases 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000000740 bleeding effect Effects 0.000 description 4
- 238000004195 computer-aided diagnosis Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 230000002207 retinal effect Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000002577 ophthalmoscopy Methods 0.000 description 2
- 206010000234 Abortion spontaneous Diseases 0.000 description 1
- 206010008190 Cerebrovascular accident Diseases 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 229920000742 Cotton Polymers 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000002490 cerebral effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 210000000416 exudates and transudate Anatomy 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 208000015994 miscarriage Diseases 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 210000001328 optic nerve Anatomy 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 208000000995 spontaneous abortion Diseases 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/12—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
- G06V10/457—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components by analysing connectivity, e.g. edge linking, connected component analysis or slices
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Heart & Thoracic Surgery (AREA)
- Probability & Statistics with Applications (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- Ophthalmology & Optometry (AREA)
- Biophysics (AREA)
- Multimedia (AREA)
- Eye Examination Apparatus (AREA)
Abstract
本发明公开了一种基于眼底图像的动静脉视网膜血管分类方法，包括：(1)获取眼底图像的全局血管集和视盘定位信息，所述的全局血管集为所述眼底图像中所有血管的集合，所述的视盘定位信息包括所述眼底图像的视盘中心；(2)根据所述的全局血管集和视盘定位信息确定主血管，并对主血管进行分类得到主血管分类信息；(3)利用所述的主血管分类信息采用基于SAT的广度搜索算法对所述的全局血管集中的血管进行分类得到全局分类信息。本发明首先获取视盘周围的主血管的分类信息，并从主血管开始基于SAT的广度搜索算法外扩扩散得到所有血管，实现了一个完整的自动血管分类方法，无需人工干预，且分类精度高。
Description
技术领域
本发明涉及计算机辅助诊断技术领域，具体涉及一种基于眼底图像的动静脉视网膜血管分类方法。
背景技术
随着计算机技术中的人工智能领域的快速发展，计算机辅助诊断技术也逐渐发展。计算机辅助诊断技术是指通过影像学、医学图像处理技术以及其他可能的生理、生化手段，结合计算机的分析计算，辅助影像科医师发现病灶，提高诊断的准确率。
通常医学影像学中计算机辅助诊断分为三步，具体如下：第一步是把病变从正常结构中提取出来；第二步是图像特征的量化；第三步是对数据进行处理并得出结论。
因为计算机可以全面利用影像信息进行精确的定量计算，去除人的主观性，避免因个人知识和经验的差异而引起的“千差万别”的诊断结果；所以它的结果是不含糊的，是确定的，它使诊断变得更为准确、更为科学。
随着现代高科技的发展，计算机辅助诊断将与图像处理和PACS***等技术融合，变得更易于操作、也更趋于准确，其临床应用范围将进一步扩大。
在医学检测中，眼睛是唯一可无损检测同时信息丰富的器官。研究指出视网膜血管病变中的血管局限缩窄、弥漫缩窄、动静脉交叉压迫、血管行走改变、铜丝动脉、出血、棉絮斑、硬性渗出以及视网膜神经纤维层缺损与脑卒有显著的相关性。且对于脑卒中的预测，眼底检查仅需40元，而MRI检查则需要上千元，颈动脉超声也需要140元。相比之下眼底检查的性价比最高。眼底图像计算机分析的全自动化的方法，包括可以提供即时的视网膜病变分类，而不需要专家意见，建立以眼底血管视神经预测三高并发症的***具有其确实的经济意义。因此，视网膜血管的病变检测在对脑卒的辅助检测具有突出作用。其中构建一个动静脉交叉压迫视网膜血管病变的自动检测***更是其中的关键部分。
对眼底图像进行血管分割、视盘定位和血管分类(动静脉***)是视网膜血管的病变检测的基础，现有的血管分割方法需要人工添加标注信息，自动化程度不高。
发明内容
针对现有技术的不足，本发明提供了一种基于眼底图像的动静脉视网膜血管分类方法。
一种基于眼底图像的动静脉视网膜血管分类方法，包括如下步骤：
(1)针对所述眼底图像，基于血管分割法获取若干个不同的全局血管集；
(2)针对每个全局血管集，基于模糊收敛算法获取对应的最终收敛区域，并根据各个最终收敛区域的重叠情况确定所述眼底图像的视盘中心；
(3)选择一个全局血管集，根据该全局血管集和视盘定位信息确定主血管，并对主血管进行分类得到主血管分类信息；
(4)利用所述的主血管分类信息采用基于SAT的广度搜索算法对所述的全局血管集中的血管进行分类得到全局分类信息。
本发明中步骤(3)中选择对应为血管的像素点数量最多的全局血管集，相应的步骤(4)中即对该选择的全局血管集进行分类。
本发明获取若干个不同的全局血管集的目的在于通过多个全局血管集确定收敛区域，根据各个全局血管集对应的收敛区域之间的重叠情况获取视盘中心，有利于保证得到的视盘中心的准确性。
所述步骤(1)基于血管分割法通过如下步骤获取若干个不同的全局血管集时，首先设定若干个不同的二值化阈值，然后针对每个二值化阈值进行如下操作：
(1-1)按照当前二值化阈值对所述的眼底图像进行二值化处理，并提取二值化处理后的眼底图像中的中心线和边缘，得到血管树；
(1-2)对所述的血管树分叉处做断开处理得到血管段，并对各个血管段进行线分割得到血管，得到原始血管集；
(1-3)确定误分割血管，并从原始血管集中去除即得到全局血管集。
本发明中二值化阈值为二值化处理后为血管的像素点个数占整个眼底图像的像素点比例，通常取值为4～20％。二值化阈值越大，则越宽松。作为优选，所述的二值化阈值为14％。
所述步骤(1-3)中的误分割血管包括分为两类，第一类误分割血管基于血管两侧的背景差异确定，第二类基于血管形状确定。
本发明中所述的背景差异值背景颜色的差异，步骤(1-3)中通过如下步骤确定第一类误分割血管：
(a1)针对每个血管，提取该血管两侧背景的特征向量；
任意一侧背景的特征向量根据该侧距离中心线5～8个像素点以内区域中的所有像素点(即与中心线的距离小于5～10个像素点)的RGB三个通道上的颜色值得到。
每侧的特征向量为一个三维向量，分别表示血管两侧背景的在RGB三个通道上的颜色值信息。具体实现时获取该侧距离中心线小于5～10个像素点的区域中的所有像素点在R、G、B三个通道上的颜色值并分别在每个通道上求平均，进而得到该侧的特征向量。
(a2)根据所述特征向量采用聚类法将所有血管聚为两类，得到的小类即为视盘周围误分割血管。
针对每个血管，计算两侧特征向量的欧氏距离，然后对所有血管对应的欧式距离进行聚类，即完成对血管的聚类。
由于K均值聚类(即K-means算法)不需要调整参数，且运行速度较快。作为优选，所述步骤(a2)中采用K均值聚类法将所有血管聚为两类。
本发明中所述的血管形状实际上指各个血管的连接关系，步骤(1-3)中通过如下步骤确定第二类误分割血管：
确定划分出原始血管集的眼底图像中的环状结构，针对各个环状结构，若该环状结构中长度最大的血管的长度小于预设的分割长度阈值，则认为该环状结构中所有血管均为第二类误分割血管，进一步进行如下操作：
确定该环状结构的中心，并计算该中心到长度大于或等于分割长度阈值的血管的最短距离(即该中心到距离其最近的长度大于或等于分割长度阈值的血管的距离)，认为以该中心为圆心、最短距离为半径的圆形区域内所有血管为第二类误分割血管。
分割长度阈值的设定根据实际经验值设定，本发明中该环状分割长度阈值α＝x/60～x/45。
本发明的血管分割方法，得到原始血管后进一步利用血管的背景和形状确定误分割的血管，可以有效去除因为拍照造成的环状反光、视盘周围的非血管的跃阶边缘、斑状病变以及出血病变等原因造成的误分割血管，大大提高了血管分割结果(即得到的全局血管集)的准确性。
所述步骤(2)中获取当前全局血管集对应的收敛区域时进行如下操作：
(2-1)对当前全局血管集中每一个血管，使用模糊收敛算法获取该血管的收敛区域；
(2-2)统计眼底图像每个像素点所对应的收敛区域的个数作为该像素点的投票值；
(2-3)选取投票值大的前n个像素点，对选取的n个像素点使用基于八连接的区域连通算法得到若干个连通区域，并以面积最大的连通区域作为最终收敛区域。
n的大小根据眼底图像的大小设定，作为优选，所述n的取值为1000～3000。
所述步骤(2)通过如下方法根据各个最终收敛区域的重叠情况确定所述眼底图像的视盘中心：
判断是否存在至少l个最终收敛区域的重叠区域，其中l＝k/3～k/2，K为二值化阈值的个数：
若存在，则以面积最大的重叠区域的中心坐标作为视盘定位信息；
否则，以采用特定模板匹配法得到视盘定位信息。
所述步骤(3)通过如下方法确定主血管：
以距离视盘中心为若干个像素点以内的区域作为视盘邻近区域，以所述的视盘邻近区域内长度大于预设的分类长度阈值的血管作为主血管。
本发明中以距离视盘中心为R个像素点以内的区域，即以视盘中心为圆心，以R为半径的区域作为视盘邻近区域，以确定的视盘邻近区域内长度大于预设的分类长度阈值的血管作为主血管。
其中，半径R和分类长度阈值的大小根据眼底图片的大小和实际情况决定。作为优选，所述R的取值为100～150，所述的分类长度阈值为50～65。
所述步骤(3)通过如下步骤对主血管进行分类得到主血管分类信息：
(3-1)获取各个主血管的平均管径，指定平均管径最大的主血管为静脉血管；
(3-2)将各个主血管切割为若干片段，得到相应的主血管切片；
(3-3)提取各个主血管切片的特征向量，并基于所述的特征向量采用聚类法将所述的主血管切片聚为两类，并以将静脉血管对应的主血管切片所在的类作为静脉血管，另一类作为动脉血管；
(3-4)针对每个主血管，以较多主血管切片所在的类作为该主血管的分类结果。
所述步骤(3-3)中通过如下方法提取各个主血管切片的特征向量：
获取距离主血管切片的血管中心若干个像素点以内的区域中所有像素点的颜色信息，并以该区域内所有像素点的颜色信息的均值作为该主血管切片的特征向量。
作为优选，提取各个主血管切片的特征向量时获取距离主血管切片的血管中心5～8个像素点以内的区域中所有像素点的颜色信息。
对全局血管集进行广度搜索，基于SAT的广度搜索算法在搜索过程中使用三条约束条件进行血管类别的传递：约束条件1：十字交叉的两根血管分别标记为两类血管；约束条件2：三岔结构中的三根血管标记为；约束条件3：三岔结构中的三根血管其中一根血管如果和剩余的两根血管夹角之和小于或等于270度时，判定为约束条件2的三岔结构，即三根血管为同一类血管，否则不做判定。
利用上述三个约束条件可以较好的区分由于血管分割时存在的漏分割而造成的交叉误判断成三叉的情况，提高分类精度。另一方面，基于该约束条件可信度高的血管先得到分类结果，可信度低的血管对于可信度高的血管传递不到的区域做补正的方式可以使得全局的血管标注的可信度提高，从而提高分类效果。
本发明现行指定较粗的主血管为静脉血管，因此在整个SAT的广度搜索算法中血管越粗，得到的该血管的分类结果的可信度越高。
与现有技术相比，本发明首先获取视盘周围的主血管的分类信息，并从主血管开始基于SAT的广度搜索算法外扩扩散得到所有血管，实现了一个完整的自动血管分类方法，无需人工干预，且分类精度高。
附图说明
图1为本实施例的眼底图像；
图2为基于广度搜索算法的眼底图像的动静脉视网膜血管分类的流程图；
图3为本实施例中对眼底图像进行血管分割的流程图；
图4为血管分割得到的原始血管集的示意图；
图5为血管分割得到的全局血管集的示意图；
图6本实施例中对眼底图像进行视盘定位的流程图；
图7为本实施例眼底图像的动静脉视网膜血管分类的全局分类信息的示意图。
具体实施方式
下面将结合附图和具体实施例对本发明进行详细描述。
本实施例以图1所示的眼底图像为例来说明基于眼底图像的动静脉视网膜血管分类方法，该眼底图像的大小为3000×3000。由拍照造成的环状反光、视盘周围的非血管的跃阶边缘、斑状病变以及出血病变等原因，该眼底图像中存在亮环。
对该眼底图像进行动静脉视网膜血管分类，分类流程如图2所示，包括如下步骤：
(1)获取眼底图像的全局血管集(即最终血管集)和视盘定位信息，全局血管集为眼底图像中所有血管的集合，视盘定位信息包括眼底图像的视盘中心；
本实施例中通过对眼底图像进行血管分割获取眼底图像的全局血管集，具体流程如图3所示，包括如下步骤：
(1-1)对眼底图像进行小波变换(IUWT小波)，按照预设的二值化阈值对经过小波变换的眼底图像进行二值化处理，并提取二值化处理后的眼底图像中的中心线和边缘，得到血管树；
(1-2)对血管树分叉处做断开处理得到血管段，并对各个血管段进行线分割得到血管，得到原始血管集。
对血管树分叉处做断开处理时：当血管树中的血管中心线中多根中心线汇集到一点时，去除中心点(汇集的交叉点)，得到单独的多根血管中心线。
对各个血管段进行线分割时：以每一根中心线作为一个血管段。血管段为一条曲线，运用图像处理的线分割的传统方法，将曲线用多根直线逼近。得到的多根直线，每根直线即代表一根血管，所有直线的集合即为原始血管集。
(1-3)确定误分割血管，本实施例中误分割血管得到第一类误分割血管和第二类误分割血管，从原始血管集合中删除第一类误分割血管和第二类误分割血管，则得到全局血管集(即最终血管集)。
对于环状反光造成的误分割，其分割出的血管相对与正常血管具有是由小段的血管组成的环的结构特点。
对于视盘周围的跃阶边缘造成的误分割，其分割出的血管在RGB色彩空间和结构上并没有特别的特点。其误分割血管为视盘周围的背景组成，因为其靠近视盘，而视盘周围的背景颜色相对远离视盘周围的背景来说和普通血管颜色具有相识性；从结构上来说由于其是孤立存在，和视盘周围血管混杂在一起也很难从结构上区分出，如果强行从结构上做判定容易造成大量的误判。但是其血管段两侧的背景在RGB色彩空间上来说具有较大的色差，这是因为其两侧背景一边由视盘而另外一边由普通背景组成。相比一般的血管，其两侧背景都是由普通背景或者都是由视盘组成。
对于斑状病变以及出血病变造成的误分割，其分割出的血管在颜色上是由普通背景组成，不具有特殊特点。但是其结构相对正常血管来说显得特别杂乱，不具有较长的血管形成的树状结构，多为多个小的环状结构和一些细碎的小血管组合而成。
基于以上分析，本实施例中基于血管两侧的背景差异确定第一类误分割血管：
(a1)针对每个血管，提取该血管两侧背景的特征向量；
获取该侧距离中心线10个像素点之类的区域中的所有像素点在R、G、B三个通道上的颜色值并分别在每个通道上求平均，进而得到该侧的特征向量。
每侧的特征向量实际上为一个三维向量，分别表示血管两侧背景的在RGB三个通道上的颜色值信息。
(a2)采用K均值聚类法将特征向量聚为两类，根据特征向量与血管的对应关系将所有血管分为两类，由于误判概率通常不会过高，因此得到的小类(即血管含量较少的血管)即为第一类误分割血管。
本实施例中通过基于血管形状确定第二类误分割血管：
(b1)确定划分出原始血管集的眼底图像中的环状结构。
具体实现时可以构建无向图G＝(V，E)，V为所有血管中心线的两个端点的集合，E为所有血管的中心线的集合，利用该无向图G＝(V，E)确定环状结构。
(b2)针对各个环状结构，若该环状结构中长度最大的血管的长度小于预设的分割长度阈值α，其中α＝x/60～x/45，(本实施例中分割长度阈值α＝x/50，x为眼底图像的横向大小，即x＝3000)，则该环状结构中所有的血管为第二类误分割血管，进一步进行如下操作：
确定该环状结构的中心，并计算该中心到长度大于或等于α的血管的最短距离(即该中心到距离其最近的长度大于或等于α的血管的距离)，以该中心为圆心、最短距离为半径的圆形区域内所有血管为第二类误分割血管。
本实施例中二值化阈值为二值化处理后为血管的像素点个数占整个眼底图像的像素点比例，通常取值为4～20％。二值化阈值越大，则越宽松。
本实施例中使用六个不同的二值化阈值，分别为4％、6％、8％、10％、12％和14％。针对每个二值化阈值均进行步骤(1-1)～(1-3)，分别对应6个全局血管集合。
本实施例中二值化阈值为14％时得到的原始血管集如图4所示，对应得到的全局血管集的示意图如图5所示。可以看出，通过去除误分割血管可以有效消除由拍照造成的环状反光、视盘周围的非血管的跃阶边缘、斑状病变以及出血病变等原因造成的干扰，提高血管分割的精确度。
本实施例中通过对眼底图像进行视盘定位得到视盘定位信息，具体流程如图6所示，针对每个全局血管集进行如下操作：
(1-2)针对当前全局血管集中每一个血管，使用模糊收敛算法获取该血管的收敛区域；
(1-3)统计眼底图像每个像素点所属于的收敛区域的个数作为该像素点的投票值，并根据各个像素点的投票值构建一个投票矩阵，对投票矩阵进行均值滤波，均值滤波时采用的均值滤波器的大小为6×6。
本实施例中构建的投票矩阵中的各个元素与眼底图像中的像素点一一对应，为对应的像素点的投票值。
(1-4)根据滤波后投票矩阵选取投票值大的前n个像素点(本实施例中n＝3000)，对选取的n个像素点使用基于八连接的区域连通算法得到若干个连通区域，以各个全局血管集对应的面积最大的连通区域作该全局血管集的最终收敛区域，判断是否存在至少l个最终收敛区域的重叠区域，其中l＝k/2，k为预设的二值化阈值的个数，即l＝3：
若存在，则以面积最大的重叠区域的中心坐标作为视盘定位信息；
否则，以采用特定模板匹配法得到视盘定位信息。
(2)根据二值化阈值最大(即最宽松)的全局血管集和视盘定位信息确定主血管，并对主血管进行分类得到主血管分类信息；
本实施例中通过如下方法确定主血管：
以距离视盘中心100个像素点以内的区域作为视盘邻近区域，以确定的视盘邻近区域内长度大于预设的分类长度阈值(本实施例中预设的分类长度阈值为60)的血管作为主血管。
本实施例中如下步骤对主血管进行分类得到主血管分类信息：
(2-1)获取各个主血管的平均管径，指定平均管径最大的主血管为静脉血管；
(2-2)将各个主血管切割为若干片段，得到相应的主血管切片；
本实施例中分割时沿着主血管的中心线，每一个像素点即为一个切片，进而将各个主血管切割(即血管切片)为若干片段。
(2-3)沿该主血管切片的血管中心向两侧分别获取5个像素点的颜色信息，并以所有像素点的颜色信息的均值作为该主血管切片的特征向量；
本实施例中该像素点的RGB值和HSL值，即得到的特征向量为6维向量，各维分别对应该像素点在R、G、B以及H、S、L通道上的颜色值。
然后，基于特征向量采用K均值聚类法将所有主血管切片聚为两类，并以将静脉血管对应的主血管切片所在的类作为静脉血管类，另一类作为动脉血管类。
(2-4)针对每个主血管，以较多主血管切片所在的类作为该主血管的分类结果。
例如对于任意一个主血管，其对应的主血管切片中有A％在动脉血管类中，B％在静脉血管类中，若A大于B，则认为该血管为动脉血管，若A小于B则认为该血管为静脉血管，否则，任意指定。
(3)利用主血管分类信息采用基于SAT的广度搜索算法对全局血管集中的血管进行分类得到全局分类信息。
对全局血管集进行广度搜索，基于SAT的广度搜索算法在搜索过程中使用三条约束条件进行血管类别的传递：约束条件1：十字交叉的两根血管分别标记为两类血管；约束条件2：三岔结构中的三根血管标记为；约束条件3：三岔结构中的三根血管其中一根血管如果和剩余的两根血管夹角之和小于或等于270度时，判定为约束条件2的三岔结构，即三根血管为同一类血管，否则不做判定。
图7为本实施例中采用基于SAT的广度搜索算法得到的全局分类信息，可以未分类(即分类后仍然没有分类信息)的血管数量较少，大大较小了因血管分割引起的三岔结构中的血管分类无法进行的情况。
未作特殊说明，本实施例中所有流程图中圆角框表示得到的结果，方角矩形表示操作。
以上所述的具体实施方式对本发明的技术方案和有益效果进行了详细说明，应理解的是以上所述仅为本发明的最优选实施例，并不用于限制本发明，凡在本发明的原则范围内所做的任何修改、补充和等同替换等，均应包含在本发明的保护范围之内。
Claims (6)
1.一种眼底图像的动静脉视网膜血管分类方法，其特征在于，包括如下步骤：
(1)针对所述眼底图像，基于血管分割法获取若干个不同的全局血管集；
(2)针对每个全局血管集，基于模糊收敛算法获取对应的最终收敛区域，并根据各个最终收敛区域的重叠情况确定所述眼底图像的视盘中心；
(3)选择一个全局血管集，根据该全局血管集和视盘定位信息确定主血管，并对主血管进行分类得到主血管分类信息；
(4)利用所述的主血管分类信息采用基于SAT的广度搜索算法对所述的全局血管集中的血管进行分类得到全局分类信息；
所述步骤(1)基于血管分割法通过如下步骤获取若干个不同的全局血管集时，首先设定若干个不同的二值化阈值，然后针对每个二值化阈值进行如下操作：
(1-1)按照当前二值化阈值对所述的眼底图像进行二值化处理，并提取二值化处理后的眼底图像中的中心线和边缘，得到血管树；
(1-2)对所述的血管树分叉处做断开处理得到血管段，并对各个血管段进行线分割得到血管，得到原始血管集；
(1-3)确定误分割血管，并从原始血管集中去除即得到全局血管集；
所述步骤(1-3)中的误分割血管包括分为两类，第一类误分割血管基于血管两侧的背景差异确定，第二类基于血管形状确定；
所述步骤(2)中获取当前全局血管集对应的收敛区域时进行如下操作：
(2-1)对当前全局血管集中每一个血管，使用模糊收敛算法获取该血管的收敛区域；
(2-2)统计眼底图像每个像素点所对应的收敛区域的个数作为该像素点的投票值；
(2-3)选取投票值大的前若干个像素点，对选取的若干个像素点使用基于八连接的区域连通算法得到若干个连通区域，并以面积最大的连通区域作为最终收敛区域；
所述步骤(2)通过如下方法根据各个最终收敛区域的重叠情况确定所述眼底图像的视盘中心：
判断是否存在至少l个最终收敛区域的重叠区域，其中l＝k/3～k/2，K为二值化阈值的个数：
若存在，则以面积最大的重叠区域的中心坐标作为视盘定位信息；
否则，以采用特定模板匹配法得到视盘定位信息；
基于SAT的广度搜索算法在搜索过程中使用三条约束条件进行血管类别的传递：约束条件1：十字交叉的两根血管分别标记为两类血管；约束条件2：三岔结构中的三根血管标记为；约束条件3：三岔结构中的三根血管其中一根血管如果和剩余的两根血管夹角之和小于或等于270度时，判定为约束条件2的三岔结构，即三根血管为同一类血管，否则不做判定。
2.如权利要求1所述的眼底图像的动静脉视网膜血管分类方法，其特征在于，所述步骤(3)通过如下方法确定主血管：
以距离视盘中心为若干个像素点以内的区域作为视盘邻近区域，以所述的视盘邻近区域内长度大于预设的分类长度阈值的血管作为主血管。
3.如权利要求2所述的眼底图像的动静脉视网膜血管分类方法，其特征在于，所述的分类长度阈值为50～65。
4.如权利要求3所述的眼底图像的动静脉视网膜血管分类方法，其特征在于，所述步骤(3)通过如下步骤对主血管进行分类得到主血管分类信息：
(3-1)获取各个主血管的平均管径，指定平均管径最大的主血管为静脉血管；
(3-2)将各个主血管切割为若干片段，得到相应的主血管切片；
(3-3)提取各个主血管切片的特征向量，并基于所述的特征向量采用聚类法将所述的主血管切片聚为两类，并以将静脉血管对应的主血管切片所在的类作为静脉血管，另一类作为动脉血管；
(3-4)针对每个主血管，以较多主血管切片所在的类作为该主血管的分类结果。
5.如权利要求4所述的眼底图像的动静脉视网膜血管分类方法，其特征在于，所述步骤(3-3)中通过如下方法提取各个主血管切片的特征向量：
获取距离主血管切片的血管中心若干个像素点以内的区域中所有像素点的颜色信息，并以该区域内所有像素点的颜色信息的均值作为该主血管切片的特征向量。
6.如权利要求5所述的眼底图像的动静脉视网膜血管分类方法，其特征在于，提取各个主血管切片的特征向量时获取距离主血管切片的血管中心5～8个像素点以内的区域中所有像素点的颜色信息。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410850207.8A CN104573712B (zh) | 2014-12-31 | 2014-12-31 | 基于眼底图像的动静脉视网膜血管分类方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410850207.8A CN104573712B (zh) | 2014-12-31 | 2014-12-31 | 基于眼底图像的动静脉视网膜血管分类方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104573712A CN104573712A (zh) | 2015-04-29 |
| CN104573712B true CN104573712B (zh) | 2018-01-16 |
Family
ID=53089732
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410850207.8A Active CN104573712B (zh) | 2014-12-31 | 2014-12-31 | 基于眼底图像的动静脉视网膜血管分类方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104573712B (zh) |
Families Citing this family (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104809480B (zh) * | 2015-05-21 | 2018-06-19 | 中南大学 | 一种基于分类回归树和AdaBoost的眼底图像视网膜血管分割方法 |
| CN106204555B (zh) * | 2016-06-30 | 2019-08-16 | 天津工业大学 | 一种结合Gbvs模型和相位一致性的视盘定位方法 |
| CN106529420B (zh) * | 2016-10-20 | 2019-07-19 | 天津大学 | 综合眼底图像边缘信息和亮度信息的视盘中心定位方法 |
| CN106846301B (zh) * | 2016-12-29 | 2020-06-23 | 北京理工大学 | 视网膜图像分类方法及装置 |
| CN106991718B (zh) * | 2017-03-31 | 2022-02-15 | 上海健康医学院 | 一种基于明暗度恢复重建眼底三维结构的方法 |
| CN107203758A (zh) * | 2017-06-06 | 2017-09-26 | 哈尔滨理工大学 | 糖尿病人视网膜血管图像分割方法 |
| CN107229937A (zh) * | 2017-06-13 | 2017-10-03 | 瑞达昇科技(大连)有限公司 | 一种视网膜血管分类方法及装置 |
| CN108182680B (zh) * | 2017-12-28 | 2021-12-28 | 中科微光医疗研究中心(西安)有限公司 | 一种基于ivoct图像的分叉血管的角度自动识别方法 |
| CN108073918B (zh) * | 2018-01-26 | 2022-04-29 | 浙江大学 | 眼底视网膜的血管动静脉交叉压迫特征提取方法 |
| CN108230322B (zh) * | 2018-01-28 | 2021-11-09 | 浙江大学 | 一种基于弱样本标记的眼底特征检测装置 |
| CN110276763B (zh) * | 2018-03-15 | 2021-05-11 | 中南大学 | 一种基于可信度和深度学习的视网膜血管分割图生成方法 |
| CN108764286B (zh) * | 2018-04-24 | 2022-04-19 | 电子科技大学 | 一种基于迁移学习的血管图像中特征点的分类识别方法 |
| CN108803994B (zh) * | 2018-06-14 | 2022-10-14 | 四川和生视界医药技术开发有限公司 | 视网膜血管的管理方法及视网膜血管的管理装置 |
| CN109635862B (zh) * | 2018-12-05 | 2021-08-24 | 合肥奥比斯科技有限公司 | 早产儿视网膜病plus病变分类方法 |
| CN111696089B (zh) * | 2020-06-05 | 2023-06-16 | 上海联影医疗科技股份有限公司 | 一种动静脉确定方法、装置、设备和存储介质 |
| CN111932554B (zh) * | 2020-07-31 | 2024-03-22 | 青岛海信医疗设备股份有限公司 | 一种肺部血管分割方法、设备及存储介质 |
| CN112734785B (zh) * | 2021-01-28 | 2024-06-07 | 依未科技(北京)有限公司 | 亚像素级眼底血管边界确定的方法、装置、介质和设备 |
| CN112734828B (zh) * | 2021-01-28 | 2023-02-24 | 依未科技(北京)有限公司 | 一种眼底血管中心线确定的方法、装置、介质和设备 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101393644A (zh) * | 2008-08-15 | 2009-03-25 | 华中科技大学 | 一种肝门静脉血管树建模方法及其*** |
| CN102014731A (zh) * | 2008-04-08 | 2011-04-13 | 新加坡国立大学 | 视网膜图像分析***和方法 |
| CN102346911A (zh) * | 2010-07-28 | 2012-02-08 | 北京集翔多维信息技术有限公司 | 在数字减影血管造影图像序列中分割血管的方法 |
-
2014
- 2014-12-31 CN CN201410850207.8A patent/CN104573712B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102014731A (zh) * | 2008-04-08 | 2011-04-13 | 新加坡国立大学 | 视网膜图像分析***和方法 |
| CN101393644A (zh) * | 2008-08-15 | 2009-03-25 | 华中科技大学 | 一种肝门静脉血管树建模方法及其*** |
| CN102346911A (zh) * | 2010-07-28 | 2012-02-08 | 北京集翔多维信息技术有限公司 | 在数字减影血管造影图像序列中分割血管的方法 |
Non-Patent Citations (2)
| Title |
|---|
| Blood vessel classification into arteries and veins in retinal images;Claudia Nieuwenhuis等;《Proceedings of SPIE - The International Society for Optical Engineering》;20070331;第1-8页，附图5和附图8 * |
| Locating the Optic Nerve in a Retinal Image Using the Fuzzy Convergence of the Blood Vessels;Adam Hoover等;《IEEE TRANSACTIONS ON MEDICAL IMAGING》;20030831;第22卷(第8期);第952-956页，附图6、9 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104573712A (zh) | 2015-04-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104573712B (zh) | 基于眼底图像的动静脉视网膜血管分类方法 | |
| CN104537669B (zh) | 眼底图像的动静脉视网膜血管分割方法 | |
| CN104545792B (zh) | 眼底图像的动静脉视网膜血管视盘定位方法 | |
| CN104573716A (zh) | 基于广度搜索算法的眼底图像的动静脉视网膜血管分类方法 | |
| AU2019205013B2 (en) | Method and system for identification of cerebrovascular abnormalities | |
| CN111046835A (zh) | 一种基于区域特征集合神经网络的眼底照多病种检测*** | |
| Seoud et al. | Automatic grading of diabetic retinopathy on a public database | |
| Muramatsu et al. | Automated detection and classification of major retinal vessels for determination of diameter ratio of arteries and veins | |
| CN111222361A (zh) | 一种高血压视网膜血管改变特征数据分析的方法和*** | |
| CN108764342B (zh) | 一种对于眼底图中视盘和视杯的语义分割方法 | |
| CN114926477A (zh) | 一种基于深度学习的脑肿瘤多模态mri影像分割方法 | |
| CN109636805A (zh) | 一种基于分类先验的宫颈图像病变区域分割装置及方法 | |
| SG172159A1 (en) | A method and system for determining the position of an optic cup boundary | |
| CN109919915A (zh) | 基于深度学习的视网膜眼底图像异常区域检测方法及设备 | |
| CN111797900B (zh) | 一种oct-a图像的动静脉分类方法及装置 | |
| Li et al. | Vessel recognition of retinal fundus images based on fully convolutional network | |
| CN107203758A (zh) | 糖尿病人视网膜血管图像分割方法 | |
| CN110136151A (zh) | 基于ct图像获取肝脏体积的***及方法 | |
| Joshi et al. | Identification and reconnection of interrupted vessels in retinal vessel segmentation | |
| CN110246126A (zh) | 一种从肺部ct图像中提取末端支气管树的方法 | |
| Zhou et al. | Computer aided diagnosis for diabetic retinopathy based on fundus image | |
| CN116152492A (zh) | 基于多重注意力融合的医学图像分割方法 | |
| Yang et al. | Detection of microaneurysms and hemorrhages based on improved Hessian matrix | |
| CN111210436B (zh) | 晶状体分割方法、装置及存储介质 | |
| Jung et al. | Evaluating a deep-learning system for automatically calculating the stroke ASPECT score |
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 | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20191008 Address after: 310000 room 410-2, floor 4, building 3, No. 9, Jiuhuan Road, Jianggan District, Hangzhou City, Zhejiang Province Patentee after: Hangzhou Qiushi innovative health technology Co., Ltd. Address before: 310027 Hangzhou, Zhejiang Province, Xihu District, Zhejiang Road, No. 38, No. Patentee before: Zhejiang University |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: Classification of arteriovenous retinal vessels based on fundus images Effective date of registration: 20210423 Granted publication date: 20180116 Pledgee: Xixi sub branch of Bank of Hangzhou Co.,Ltd. Pledgor: Hangzhou Qiushi innovative health technology Co.,Ltd. Registration number: Y2021330000329 |
|
| PC01 | Cancellation of the registration of the contract for pledge of patent right | ||
| PC01 | Cancellation of the registration of the contract for pledge of patent right |
Date of cancellation: 20211216 Granted publication date: 20180116 Pledgee: Xixi sub branch of Bank of Hangzhou Co.,Ltd. Pledgor: Hangzhou Qiushi innovative health technology Co.,Ltd. Registration number: Y2021330000329 |