基于深度学习的动脉瘤检测方法

1.本发明涉及图像处理技术领域，具体涉及一种基于深度学习的动脉瘤检测方法。


背景技术：

2.动脉瘤是由于动脉壁的病变或者损伤形成动脉壁局限性或者弥漫性的扩张而形成，根据动脉瘤出现部位不同可分为周围动脉瘤、腹主动脉瘤、胸腹主动脉瘤、主动脉夹层动脉瘤、内脏动脉瘤等，主要表现为体表搏动性肿块、动脉瘤压迫周围神经或破裂时出现剧烈疼痛，瘤腔内血栓或斑块脱落致远端动脉栓塞产生肢体、器官缺血或坏死等。动脉瘤由于早期症状不明显，因此目前动脉瘤在我国的筛查率较低，一方面是由于患者在病变早期无明显症状对病情不够重视，另一方面从客观上则是因为医疗条件的不足。
3.临床上对动脉瘤病变的检测主要依赖于人工检测，首先利用医疗设备采集患者动脉图像，然后需要经验丰富的医生对动脉图像进行评估，在人工检测过程中既耗费医疗设备等物理资源，同时对临床医生的专业性有较高的要求。现有技术中借助深度学习网络模型只是简单的识别出含有动脉瘤的区域，但是并未能够直接反映出动脉瘤病变形状，难以快速的判断出动脉瘤的危害程度。


技术实现要素：

4.本发明提供一种基于深度学习的动脉瘤检测方法，该方法通过在动脉内注射造影剂并使用x射线对其进行成像获得动脉血管图像，而后对图像进行动脉瘤的检测，所以通过深度学习的方法来实现对动脉瘤的检测并获取该动脉瘤的特征标注于图像的动脉瘤旁边，为医生提供参考，从而实现对动脉血管图像的快速判断。
5.本发明的目的是提供一种基于深度学习的动脉瘤检测方法，包括以下步骤：获取血管造影图像；对造影图像通过卷积神经网络获取含动脉瘤的血管区域；获取血管区域中动脉瘤的瘤颈宽度；并获取血管区域中动脉瘤边缘上任意两个像素点的距离，从任意两个像素点的距离中获取最大距离；根据动脉瘤的瘤颈宽度及最大距离获取血管区域中动脉瘤的脆弱程度；根据与血管区域中动脉瘤相对一侧的第一血管边缘上的像素点通过模拟获取动脉瘤所在侧的第二血管边缘；根据动脉瘤边缘上相邻像素点至第二血管边缘的距离差值获取动脉瘤边缘上每个像素点的起伏值；根据动脉瘤边缘上每个像素点的起伏值获取血管区域中动脉瘤的起伏程度；获取动脉瘤边缘两端处对应血管的直径；根据动脉瘤边缘两端处对应血管的直径，动脉瘤的脆弱程度和起伏程度，以及动脉瘤边缘上每个像素点至第二血管边缘的距离获取动脉瘤的危害指数。
6.在一实施例中，所述动脉瘤所在侧的第二血管边缘是按照以下步骤获取：根据动脉瘤边缘两端处的像素点均与其相邻的血管边缘像素点获取动脉瘤边缘
两端处像素点均与其相邻血管边缘像素点连线的第一斜率；根据动脉瘤边缘两端处像素点均与其相邻的血管边缘像素点连线的第一斜率获取垂直于血管并分别经动脉瘤边缘两端处像素点的第一直线斜率；根据分别经动脉瘤边缘两端处像素点的第一直线斜率获取垂直于血管并分别经动脉瘤边缘两端处像素点的两条直线；获取两条直线与血管区域中动脉瘤相对一侧的第一血管边缘各有一个交点；根据该两个交点间与血管区域中动脉瘤相对一侧的第一血管边缘上的像素点通过模拟获取动脉瘤所在侧的第二血管边缘。
7.在一实施例中，所述动脉瘤所在侧的第二血管边缘的获取是通过第一血管边缘上的相邻像素点的斜率模拟而获取的。
8.在一实施例中，根据两条直线与动脉瘤边缘两端处的血管两侧相交的线段长度获取动脉瘤边缘两端处对应血管的直径。
9.在一实施例中，所述血管区域中动脉瘤的起伏程度是按照以下步骤获取：根据动脉瘤边缘上每个像素点的起伏值，通过统计将相同起伏值作为一个起伏级数；根据每个起伏级数中起伏值的个数占起伏值总个数的比例获取血管区域中动脉瘤的起伏程度。
10.在一实施例中，所述血管区域中动脉瘤的起伏程度的计算公式如下：式中，表示动脉瘤的起伏程度；表示第个起伏级数中起伏值的个数占起伏值总个数的比例；表示起伏级数的总数量。
11.在一实施例中，所述动脉瘤的危害指数计算公式如下：式中，表示动脉瘤的危害指数；表示动脉瘤的起伏程度；表示动脉瘤的脆弱程度；、分别表示动脉瘤边缘两端处对应血管的直径；表示动脉瘤边缘上每个像素点至第二血管边缘的距离的平均距离。
12.在一实施例中，所述动脉瘤边缘上每个像素点至第二血管边缘的距离的平均距离计算公式如下：式中，表示动脉瘤边缘上第个像素点至第二血管边缘的距离；表示动脉瘤边
缘上的总像素点个数。
13.在一实施例中，所述动脉瘤边缘上每个像素点至第二血管边缘的距离是按照以下步骤获取：获取动脉瘤边缘两端处像素点连线的第二斜率；根据第二斜率获取垂直于动脉瘤边缘两端处像素点连线的第二直线斜率；根据第二直线斜率获取经动脉瘤边缘上每个像素点，并垂直于动脉瘤边缘两端处像素点连线的直线；根据该直线与第二血管边缘及动脉瘤边缘上的交点获取动脉瘤边缘上每个像素点至第二血管边缘的距离。
14.在一实施例中，根据动脉瘤边缘上每个像素点至第二血管边缘的距离获取动脉瘤边缘上相邻像素点至第二血管边缘的距离差值。
15.本发明的有益效果是：本发明提供的一种基于深度学习的动脉瘤检测方法，该方法通过构建深度网络模型，利用学习算法让模型自动学习出好的特征表示，避免了医生手动提取特征的繁琐流程。基于卷积神经网络获取湖造影图像中的含有动脉瘤的血管区域，在通过计算动脉瘤的脆弱程度以及凸起程度来获得动脉瘤的特征，基于这些特征进而分析该动脉瘤的危害程度，再将分析获取的危害指数将其数值标注于造影图像中对应的动脉瘤旁边为医生提供参考，从而实现对动脉血管造影图像的快速判断。
附图说明
16.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
17.图1为本发明的一种基于深度学习的动脉瘤检测方法的实施例总体步骤的流程示意图。
18.图2为含动脉瘤的血管区域的边缘图。
19.图3为模拟出原始的动脉血管壁的第二血管边缘图。
具体实施方式
20.下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
21.本发明在动脉瘤检测过程中，通过在动脉内注射造影剂并使用x射线对其进行成像获得动脉血管图像，而后对图像进行动脉瘤的检测，所以通过深度学习的方法来实现对动脉瘤的检测并获取该动脉瘤的特征标注于图像的动脉瘤旁边为医生提供参考从而实现对动脉血管图像的快速判断。
22.本发明提供的一种基于深度学习的动脉瘤检测方法，参见图1~3所示，包括以下步骤：
s1、获取血管造影图像；在本实施例中，需要实现对动脉瘤的检测与特征的获取，所以需要获取患者的动脉血管造影图像，对患者动脉血管的获取则需要对其进行动脉血管造影检查，通过对患者动脉注入造影剂，然后通过血管造影机对其扫描成像获得患者的动脉血管造影图像。
23.s2、对造影图像通过卷积神经网络获取含动脉瘤的血管区域；基于深度学习的动脉瘤检测体系基础理论：通过深度学习检测血管造影图像中的动脉瘤，通过对动脉瘤的特征进行获取并将特征值标注于图像的该动脉瘤旁边为医生提供参考，从而实现快速地对患者病情的诊断。
24.在本实施例中，由于血管造影图像中的一些毛细血管图像存在分辨率较低、噪声较多且背景复杂等问题，一定程度上降低了动脉瘤的检测精度。首先对获取的血管造影图像进行预处理；中值滤波对噪声有良好的滤除作用，特别是在滤除噪声的同时，能够保护信号的边缘，使之不被模糊。所以选取中值滤波对血管造影图像进行滤除噪声处理，使用的中值滤波如下：式中，分别为原始图像和处理后图像。为二维模板，为模板内的点，其中为二维模板，通常为3
×
3，5
×
5区域，也可以是不同的形状，如线状，圆形，十字形，圆环形等。
25.再通过直方图均衡化处理使得动脉血管在血管造影图像中更为清晰。
26.对于图像中的动脉瘤，通过深度学习的方法将其检测出来，深度学习的思想和人的思维较为类似，通过各个层次对图像的处理和特征提取再传往下一层，最后完成分类或识别的任务。
27.获取训练好的卷积神经网络来识别血管造影图像中含动脉瘤的血管区域；卷积神经网络训练过程是，采用医院已经人为标记过的血管造影图像作为训练的数据集，将数据集中70%用于训练，20%用于测试，10%用于验证，由此训练卷积神经网络来识别血管造影图像中的动脉瘤，若测试集与验证集的准确率都达到99%以上则训练完成，此时将血管造影图像输入，则会输出标记的含动脉瘤的血管区域，训练使用的损失函数为交叉熵损失函数。
28.s3、获取血管区域中动脉瘤的瘤颈宽度；并获取血管区域中动脉瘤边缘上任意两个像素点的距离，从任意两个像素点的距离中获取最大距离；根据动脉瘤的瘤颈宽度及最大距离获取血管区域中动脉瘤的脆弱程度；需要说明的是，在本实施例中，通过canny边缘检测获得含动脉瘤的血管区域的边缘图像；参见图2所示，表示含动脉瘤的血管区域的边缘图像；图2中，动脉血管100上突出有动脉瘤200，其中a和b表示动脉瘤边缘两端处像素点；其动脉血管壁因发生病变形状酷似肿瘤因而被称为动脉瘤，通过动脉瘤的脆弱程度以及凸起程度等来获得动脉瘤的特征，根据这些特征进而分析该动脉瘤的危害程度。
29.首先获取该动脉瘤的瘤颈宽度，其瘤颈宽度为动脉瘤占据动脉血管起始位置的距
离，也就是计算a和b之间的距离，其瘤颈宽度的计算公式如下：式中，表示瘤颈宽度，和为动脉瘤占据血管的起始位置的坐标，也就是a处像素点和b处像素点的坐标。
30.在本实施例中，通过计算该动脉瘤的最大径来判断其脆弱程度，早期的动脉瘤仅是轻微的突起，若其瘤壁韧性降低时，经过血流的冲击其形状改变导致最大径变大，此时就易发生破裂而危害患者的生命。早期的动脉瘤，经过逐渐的病变导致瘤壁韧性降低；早期的动脉瘤最大径为瘤颈宽度，韧性降低之后其最大径可能会发生改变从而大于瘤颈宽度，通过动脉瘤的最大径的变化来判断其韧性程度。由此，通过计算动脉瘤边缘两两像素点的距离获取血管区域中动脉瘤边缘上任意两个像素点的距离，从而获得一系列动脉瘤边缘点之间的距离：式中，为动脉瘤的边缘点距离的序列；从任意两个像素点的距离中获取最大距离，将该最大距离为最大径如下：式中，为序列中的最大值，从而获得此时动脉瘤的脆弱程度：式中，表示最大径；表示瘤颈宽度；表示动脉瘤的脆弱程度。
31.s4、根据与血管区域中动脉瘤相对一侧的第一血管边缘上的像素点通过模拟获取动脉瘤所在侧的第二血管边缘；根据动脉瘤边缘上相邻像素点至动脉瘤边缘两端处像素点连线的距离差值获取动脉瘤边缘上每个像素点的起伏值；根据动脉瘤边缘上每个像素点的起伏值获取血管区域中动脉瘤的起伏程度；并获取动脉瘤边缘两端处对应血管的直径；动脉瘤所在侧的第二血管边缘是按照以下步骤获取：根据动脉瘤边缘两端处的像素点均与其相邻的血管边缘像素点获取动脉瘤边缘两端处像素点均与其相邻血管边缘像素点连线的第一斜率；根据动脉瘤边缘两端处像素点均与其相邻的血管边缘像素点连线的第一斜率获取垂直于血管并分别经动脉瘤边缘两端处像素点的第一直线斜率；根据分别经动脉瘤边缘两端处像素点的第一直线斜率获取垂直于血管并分别经动脉瘤边缘两端处像素点的两条直线；获取两条直线与血管区域中动脉瘤相对一侧的第一血管边缘各有一个交点；根据该两个交点间与血管区域中动脉瘤相对一侧的第一血管边缘上的像素点通过模拟获取动脉瘤所在侧的第二血管边缘。动脉瘤所在侧的第二血管边缘的获取是通过第一血管边缘上的相邻像素点的斜率模拟而获取的。
32.在本实施例中，动脉瘤的起伏程度也对患者有着重要的影响，若动脉瘤的瘤壁起伏较大，说明抗血流冲击的能力越低，越容易破裂，其突起高度为动脉瘤边缘上的点到理想情况下无动脉瘤时对应血管壁上的点的距离，因此要获得突起程度需要先获取理想状态下的动脉血管壁的情况；计算动脉瘤边缘两端处像素点均与其相邻血管边缘像素点连线的第一斜率；见图2所示，也就是分别计算a点与其相邻血管边缘像素点连线的斜率，设与a点相邻血管边缘像素点的坐标为；和b点与其相邻血管边缘像素点连线的斜率，设与b点相邻血管边缘像素点的坐标为；计算公式如下：式中，表示a点与其相邻血管边缘像素点连线的斜率；表示与a点相邻血管边缘像素点的坐标；表示a点的坐标；式中，表示b点与其相邻血管边缘像素点连线的斜率；表示与b点相邻血管边缘像素点的坐标；表示b点的坐标；参见图3所示，表示模拟出原始的动脉血管壁的第二血管边缘图；计算获取垂直于血管并经动脉瘤边缘a点处的斜率为；垂直于血管并经动脉瘤边缘b点处的斜率为；根据斜率为和获取两直线p1和p2，这两条直线p1和p2与血管区域中动脉瘤相对一侧的第一血管边缘u1各有一个交点，这两交点之间的像素点则为血管壁上的像素点，同时该两条直线与血管两侧相交的线段长度为血管直径，分别记为、。再计算这两交点之间的像素点斜率模拟出动脉瘤所在侧的第二血管边缘u2，其计算公式如下：式中，表示第一血管边缘u1上第个像素点的斜率，表示第个像素点的坐标，表示第个像素点的坐标，通过每一个像素点的斜率在动脉瘤两端之间（a与b点之间）模拟出原始的动脉血管壁的第二血管边缘u2。
33.另外，动脉瘤边缘上每个像素点至动脉瘤边缘两端处像素点连线的距离是按照以下步骤获取：获取动脉瘤边缘两端处像素点连线的第二斜率，也就是瘤颈宽度所在直线的斜
率，计算公式如下：式中，为瘤颈宽度所在直线的斜率；，分别表示动脉瘤边缘两端处像素点的坐标，也就是a和b点处的坐标；根据第二斜率获取垂直于动脉瘤边缘两端处像素点连线的第二直线斜率；具体是过动脉瘤边缘上的每一点第二直线斜率为；根据第二直线斜率获取经动脉瘤边缘上每个像素点，并垂直于动脉瘤边缘两端处像素点连线的直线；根据该直线与第二血管边缘及动脉瘤边缘上的交点获取动脉瘤边缘上每个像素点至动脉瘤边缘两端处像素点连线的距离；也就是表现为动脉瘤上与其对应血管壁上的交点的距离则为动脉瘤的凸起高度；根据动脉瘤边缘上每个像素点至第二血管边缘的距离获取动脉瘤边缘上相邻像素点至第二血管边缘的距离差值。
34.具体的，动脉瘤的起伏值计算公式为：式中，表示动脉瘤边缘上的第个像素点的起伏值，表示动脉瘤边缘上的第个像素点至动脉瘤边缘两端处像素点连线的距离；表示动脉瘤边缘上的第个像素点至动脉瘤边缘两端处像素点连线的距离；血管区域中动脉瘤的起伏程度是按照以下步骤获取：根据动脉瘤边缘上每个像素点的起伏值，通过统计将相同起伏值作为一个起伏级数；根据每个起伏级数中起伏值的个数占起伏值总个数的比例获取血管区域中动脉瘤的起伏程度。
35.血管区域中动脉瘤的起伏程度的计算公式如下：式中，表示动脉瘤的起伏程度；表示第个起伏级数中起伏值的个数占起伏值总个数的比例；表示起伏级数的总数量。
36.s5、根据动脉瘤边缘两端处对应血管的直径，动脉瘤的脆弱程度和起伏程度，以及动脉瘤边缘上每个像素点至动脉瘤边缘两端处像素点连线的距离获取动脉瘤的危害指数。
37.动脉瘤的危害指数计算公式如下：
式中，表示动脉瘤的危害指数；表示动脉瘤的起伏程度；表示动脉瘤的脆弱程度；、分别表示动脉瘤边缘两端处对应血管的直径；表示动脉瘤边缘上每个像素点至第二血管边缘的距离的平均距离。其中，的值越大说明动脉瘤抗血流冲击的能力越低，越容易破裂，所表示的动脉瘤起伏程度，早期动脉瘤其表面较为光滑，起伏程度值较低，而起伏程度越大说明动脉瘤表面起伏程度越复杂，越可能是出于中后期的动脉瘤，其破裂的可能性也越大，也就反映的值也越大，从而能够综合反映动脉瘤抗血流冲击的能力；对于、为动脉瘤边缘两端处对应血管的直径，为动脉瘤边缘上每个像素点至第二血管边缘的距离的平均距离，也就是反映动脉瘤的平均高度，则表示动脉瘤的平均相对高度，值越大说明动脉瘤高度越高，血管承受的压力就越大；基于综合反映动脉瘤抗血流冲击的能力及动脉瘤的平均相对高度对动脉瘤的影响因素综合考虑获得动脉瘤的危害指数。
38.此外，动脉瘤边缘上每个像素点至第二血管边缘的距离的平均距离计算公式如下：式中，表示动脉瘤边缘上第个像素点至动脉瘤边缘两端处像素点连线的距离；表示动脉瘤边缘上的总像素点个数。
39.从而获得了动脉瘤的危害指数。
40.通过上述步骤获得的危害指数将其数值标注于造影图像中对应的动脉瘤旁边为医生提供参考，医生可以结合图中动脉瘤旁标注的危害指数便于对患者病情更快做出判断。
41.综上，本发明提供的一种基于深度学习的动脉瘤检测方法，该方法通过构建深度网络模型，利用学习算法让模型自动学习出好的特征表示，避免了医生手动提取特征的繁琐流程。基于卷积神经网络获取湖造影图像中的含有动脉瘤的血管区域，在通过计算动脉瘤的脆弱程度以及凸起程度来获得动脉瘤的特征，基于这些特征进而分析该动脉瘤的危害程度，再将分析获取的危害指数将其数值标注于造影图像中对应的动脉瘤旁边为医生提供参考，从而实现对动脉血管造影图像的快速判断。
42.以上仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。