基于深度神经网络的颅内血管树分割方法和分割装置与流程

1.本技术涉及医学影像分割技术领域，特别是涉及一种基于深度神经网络的颅内血管树分割方法和分割装置。


背景技术：

2.脑血管疾病由于其高发病率和死亡率已成为人类健康的主要威胁。在目前的临床应用中，脑血管疾病的可视化成像诊断技术包括数字减影血管造影（dsa）、磁共振血管造影（mra）、经颅多普勒（tcd）。dsa作为诊断脑血管疾病的金标准，被广泛应用于临床，因为它允许高分辨率的评估，并可以为诊断提供直观的血管结构信息。然而，由于脑血管分支的复杂性，脑血管dsa中不仅有许多微小的毛细血管，而且颈内动脉（ica）也较厚。这就要求在脑血管分割任务中应考虑不同尺度的分割目标。同时，在dsa成像过程中，造影剂的浓度会不断变化，扩散不均匀，可能导致脑血管dsa图像中出现伪影和噪声。不仅如此，颅内全血管树往往有众多的细小血管，用深度学习方法训练网络需要大量的数据，因此标注数据往往成为限制该技术方案的重要问题。


技术实现要素：

3.基于此，有必要基于上述技术问题，提供一种基于深度神经网络的颅内血管树分割方法。
4.本技术基于深度神经网络的颅内血管树分割方法，包括利用训练完成的深度神经网络从输入的dsa影像中分割出颅内血管树，训练过程包括：对含有血管的dsa图像进行采样，从采样图像中筛选获得血管密集图像，将所述血管密集图像划分为小血管样本图和大血管样本图；根据典型角度造影中不同尺寸血管的相对位置关系，对各所述血管密集图像进行铺列组合，用于形成样本图像；在所述样本图像中获得每个像素点的离散度信息，所述离散度信息用于表示预定区域内血管像素值的离散程度；将所述样本图像、以及所述样本图像中每个像素点对应的所述离散度信息，结合制作为输入数据集，所述输入数据集用于训练所述深度神经网络。
5.可选的，所述采样图像的滑动步长，为所述采样图像尺寸长或宽的一半；颅内血管树分割方法包括：对不含血管的dsa图像进行采样，筛选获得含有噪声分布区域的负样本图，所述负样本图用于补充铺列、以形成所述样本图像。
6.可选的，筛选获得含有噪声分布区域的负样本图，具体包括：基于不含血管的dsa图像，对所有像素点的灰度值进行由大到小的排布分析，获得灰度值排名在后的、10%数量的像素点，作为待选位置；在所述待选位置中，选择负样本像素点，进而获得以所述负样本像素点为中心的负样本图。
7.可选的，所述颅内血管树分割方法包括将dsa序列帧图像划分为含有血管的dsa图像、以及不含血管的dsa图像，具体包括：基于dsa序列帧图像中的每帧dsa图像，统计整幅图像中全部像素点的灰度值，对灰度值进行由大到小的排布分析；获得灰度值排名在前的、第一百分比数量的灰度值均值；获得灰度值排名在后的、第二百分比数量的灰度值中位数值，所述第二百分比与第一百分比的加和为一；将所述灰度值均值和所述灰度值中位数值做差，若差值大于给定阈值，则判定该帧dsa图像为含有血管的dsa图像，若差值小于或等于给定阈值，则判定该帧dsa图像为不含血管的dsa图像。
8.可选的，所述采样图像大小为32*32，所述血管密集图像的判断标准为采样图像内血管像素点的占比率大于预期值，所述小血管样本图和大血管样本图依据所述血管密集图像中血管像素点的离散度划分；所述第一百分比为10%，所述第二百分比为90%。
9.可选的，在所述样本图像中获得每个像素点的离散度信息，所述离散度信息用于表示预定区域内血管像素值的离散程度，具体包括：针对每个像素点，获得以该像素点为中心的第一区域；在各所述第一区域中，筛选获得包含血管像素点的第二区域；将所述第二区域中血管像素点的离散度信息，作为该像素点的离散度信息。
10.可选的，在各所述第一区域中，筛选获得包含血管像素点的第二区域，具体包括：基于各所述第一区域内的图像，统计图像中全部像素点的灰度值，对灰度值进行由大到小的排布分析；获得灰度值排名在前的、第一百分比数量的灰度值均值；获得灰度值排名在后的、第二百分比数量的灰度值中位数值，所述第二百分比与第一百分比的加和为一；将所述灰度值均值和所述灰度值中位数值做差，若差值大于给定阈值，则判定该第一区域内的图像含有血管像素点，筛选获得包含血管像素点的第二区域。
11.可选的，所述深度神经网络为生成对抗网络，所述生成对抗网络用于对所述样本图像进行相似性判别和重构，通过相似性判别的样本图像用于制作输入数据集。
12.可选的，所述训练过程包括：获得所述样本图像中每个像素点对应的所述离散度信息，获得统计区间；根据所述离散度信息的分布，将所述统计区间按划分为若干个子区间；针对每个子区间，分别计算深度神经网络的损失函数，以进行所述深度神经网络的训练。
13.本技术还提供一种基于深度神经网络的颅内血管树分割装置，包括存储器、处理器及存储在存储器上的计算机程序，所述处理器执行所述计算机程序以实现本技术所述的基于深度神经网络的颅内血管树分割方法的步骤。
14.本技术还提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现本技术所述的基于深度神经网络的颅内血管树分割方法的步骤。
15.本技术还提供一种计算机程序产品，包括计算机指令，该计算机指令被处理器执行时实现本技术所述的基于深度神经网络的颅内血管树分割方法的步骤。
16.本技术基于深度神经网络的颅内血管树分割方法至少具有以下效果：本技术的样本图像是铺列组合形成的，不同的血管密集图像可通过排列组合形成不同的样本图像。在同样数量的标注数据下，本技术可提供更多的训练样本，可以大大降低标注成本需求数量以及数据采集难度。
17.本技术用于训练深度神经网络的输入数据集，不仅包括样本图像，还包括离散度信息，能够使深度神经网络在降采样过程中额外的关注离散程度高的样本，进而提高深度神经网络的分割精度。
18.本技术通过负样本提高训练完成的深度神经网络的抗噪声和干扰能力，从而极大地减弱无关区域的噪声影响，对主要血管及细小分支都有较好的分割效果。
19.本技术通过对不同离散度信息进行子区间划分，并分别针对不同的的子区间，分区计算损失函数，能够进一步提高抗噪声和干扰能力，对主要血管及细小分支都有较好的分割效果。
附图说明
20.图1为本技术一实施例中基于深度神经网络的颅内血管树分割方法的流程示意图；图2为本技术一实施例中基于深度神经网络的颅内血管树分割方法的流程示意图；图3为本技术一实施例中基于深度神经网络的颅内血管树分割方法的流程示意图；图4为一个实施例中计算机设备的内部结构图。
具体实施方式
21.如何利用深度神经网络尽可能少的利用标注从dsa图像中自动准确地分割脑血管成为一个具有挑战性的问题。本技术旨在克服现有技术存在的问题和缺陷，提出了一种基于深度神经网络的颅内血管树分割方法。
22.为了使本技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处描述的具体实施例仅用以解释本技术，并不用于限定本技术。
23.参见图1，本技术一实施例中提供一种基于深度神经网络的颅内血管树分割方法，包括利用训练完成的深度神经网络从输入的dsa影像中分割出颅内血管树，训练过程包括：对含有血管的dsa图像进行采样，从采样图像中筛选获得血管密集图像，将血管密集图像划分为小血管样本图和大血管样本图；根据典型角度造影中不同尺寸血管的相对位置关系，对各血管密集图像进行铺列组合，用于形成样本图像；在样本图像中获得每个像素点的离散度信息，离散度信息用于表示预定区域内血管像素值的离散程度；
将样本图像、以及样本图像中每个像素点对应的离散度信息，结合制作为输入数据集，输入数据集用于训练深度神经网络。
24.参见图2和图3，在一个实施例中，还提供一种基于深度神经网络的颅内血管树分割方法，用于进一步地解释本技术上一实施例的详细流程。包括（1）灰度信息分析与分patch采样；（2）基于gan的分patch样本重构；（3）包含离散度分析的分支监督。其中，“patch”意为采样窗口。“分patch”意为多个采样窗口分别执行采样过程，例如通过滑窗等形式。“分patch样本”意为采样窗口获得的采样图像。“分patch样本重构”意为利用各个采样窗口所获得的采样图像完成样本重构，以形成样本图像。
25.（1）灰度信息分析与分patch采样，对应步骤s100～步骤s200。
26.步骤s100，将dsa序列帧图像划分为含有血管的dsa图像、以及不含血管的dsa图像。步骤s100具体包括步骤s110～步骤s140。
27.步骤s110，基于dsa序列帧图像中的每帧dsa图像，统计整幅图像中全部像素点的灰度值，对灰度值进行由大到小的排布分析。
28.步骤s120，获得灰度值排名在前的、第一百分比数量的灰度值均值。
29.步骤s130，获得灰度值排名在后的、第二百分比数量的灰度值中位数值，第二百分比与第一百分比的加和为一。例如第一百分比为10%，第二百分比为90%。
30.步骤s140，将灰度值均值和灰度值中位数值做差，若差值大于给定阈值，则判定该帧dsa图像为含有血管的dsa图像，若差值小于或等于给定阈值，则判定该帧dsa图像为不含血管的dsa图像。
31.具体地，对dsa序列帧图像进行灰度信息分析，统计整幅图像中的全部像素点的像素值，找出灰度值最高的前10%，求出其均值；再将后90%灰度值的中位数值与均值比较。若二者的差值小于给定阈值，则认为此图像中不包含任何血管。若差值大于给定值，则认为此图像中包含血管。
32.给定阈值可根据第一百分比和第二百分比的取值，适应性调整。可以理解，对于dsa图像而言，流经造影剂的血管区域灰度值较低，无血管的背景区域灰度值相对较高。相比于不含血管的dsa图像，含有血管的dsa图像其像素点的灰度值信息受到血管区域的较低灰度值影响，中位数也较低。含有血管的dsa图像，在执行步骤s140时差值较大。本实施例利用上述流程，能够dsa图像自动划分为含有血管和不含血管两个类别。
33.步骤s200，对含有血管的dsa图像进行采样，从采样图像中筛选获得血管密集图像，将血管密集图像划分为小血管样本图和大血管样本图。
34.采样图像的滑动步长可以是采样图像尺寸长或宽的一半。例如采样图像大小为32*32，滑动步长为16。血管密集图像的判断标准为采样图像内血管像素点的占比率大于预期值，例如大于50%。
35.具体地，对含有血管的dsa图像进行采样，例如针对步骤s130中排名在后的、第二百分比数量的灰度值所对应的像素点进行采样，即对后90%的像素点集进行采样。采样方式例如从左上角按照从左到右从上到下的方式进行滑窗，滑窗大小为采样patch大小32*32，滑动步长为16，滑动后先判定采样图像内包含的像素点集内像素点的个数占比，如果大于50%，则保留。对于保留的采样图像，再分析采样图像内像素点集内像素点的坐标离散程度，若离散程度相对高，则判定为小血管样本图，若离散程度相对低，则判定为大血管样本图。
小血管样本图和大血管样本图依据血管密集图像中血管像素点的离散度划分。大血管样本图是指采样图像中所包含的血管，是由单根或少数根相对粗壮的血管分支组成；小血管样本图是指采样图像中所包含的血管由多根相对细小的血管分支组成。相对高低的判断标准，可由操作者设置或通过有限次实验测试获得。
36.（2）基于gan的分patch样本重构，对应步骤s310～步骤s320。
37.步骤s310，根据典型角度造影中不同尺寸血管的相对位置关系，对各血管密集图像进行铺列组合，用于形成样本图像。
38.本领域技术人员可知，典型角度造影中，大血管分布在中间和中下部位，中动脉区域也分布有一些大血管，中动脉附近存在大量的小血管，整体形成不同尺寸血管的相对位置关系。因此在形成样本图像时，可在大血管分布的位置相应地铺列大血管样本图，在小血管分布的位置相应地铺列小血管样本图。
39.本实施例中的样本图像是铺列组合形成的，不同的标注数据可通过排列组合形成不同的样本图像，标注数据包括大血管样本图、小血管样本图、以及下文负样本图。在同样数量的标注数据下，本实施例可提供更多的训练样本，可以大大降低标注成本需求数量以及数据采集难度。相比于一般的深度神经网络训练，所需的标注数据量需求比较小，能够利用更少的标注数据更好地提取复杂影像的颅内血管。
40.步骤s320，对不含血管的dsa图像进行采样，筛选获得含有噪声分布区域的负样本图，负样本图用于补充铺列、以形成样本图像。负样本图可以和血管密集图像的大小相同，以便于铺列形成样本图像，例如二者均为32*32像素。
41.在步骤s320中，筛选获得含有噪声分布区域的负样本图，具体包括：步骤s321，基于不含血管的dsa图像，对所有像素点的灰度值进行由大到小的排布分析，获得灰度值排名在后的、10%数量的像素点，作为待选位置；步骤s322，在待选位置中，选择负样本像素点，进而获得以负样本像素点为中心的负样本图。
42.本实施例通过负样本提高训练完成的深度神经网络的抗噪声和干扰能力，从而极大地减弱无关区域的噪声影响，对主要血管及细小分支都有较好的分割效果。
43.本实施例所指的是噪声是指在深度神经网络分割血管区域时，容易对深度神经网络产生干扰的区域。可以理解，对于不含血管的dsa图像，其像素值较低的区域，更加邻近于血管部分的像素值，该区域即作为噪声分布区域的至少一部分。本实施例将噪声分布区域作为负样本，使得深度神经网络在训练过程中能够具有识别噪声分布区域的能力，降低误分割的可能性。
44.具体地，噪声分布区域包括：灰度值最高的前10%，中间区域的10%，最低的后10%共三部分。在三部分的像素点集合中随机选取10个像素点，以像素点为中心进行采样，采样的patch大小为32*32，采出的样本归类为负样本。
45.本实施例中，为了充分的保留血管树的总体信息，利用分patch样本进行随机组合来构建新样本图像。样本图像尺寸例如可以是为512*512，需要由256个32*32的分patch样本组成。
46.首先，按照颅内血管造影的通常情况，以典型角度造影为例，大血管分布在中间和中下部位，对给定区间内从中间向下填充大血管分patch样本（填充大血管样本图），这些分patch样本随机选自于不同的影像序列。对于中动脉区域，也需要向两侧填充一定数量的大
血管分patch样本。其次，中动脉附近存在大量的小血管样本，需要用小血管分patch样本（小血管样本图）来填充。最后，剩余区域用负样本进行填充。
47.进一步地，本实施例中的深度神经网络为生成对抗网络（generative adversarial network，gan），生成对抗网络用于对样本图像进行相似性判别和重构，通过相似性判别的样本图像用于制作输入数据集。生成对抗网络可用于提高样本图像的可靠性。
48.由于上述构造的样本图像会和真实的病例有一定的差距，因此通过gan网络将这些构造样本与真实病例进行相似性判别和重构，最终生成满足需求的类真实样本。
49.（3）包含离散度分析的分支监督，对应步骤s400～步骤s600。
50.步骤s400，在样本图像中获得每个像素点的离散度信息，离散度信息用于表示预定区域内血管像素值的离散程度。
51.步骤s400，具体包括：步骤s410，针对每个像素点，获得以该像素点为中心的第一区域；步骤s420，在各第一区域中，筛选获得包含血管像素点的第二区域；步骤s430，将第二区域中血管像素点的离散度信息，作为该像素点的离散度信息。可以理解，在各第一区域，不含血管像素点的区域，无须计算离散度信息。第一区域的大小例如可以是32*32像素。
52.进一步地，步骤s420，在各第一区域中，筛选获得包含血管像素点的第二区域，具体包括步骤s421～步骤s424。步骤s421～步骤s424的判定原则，与步骤s100及其子步骤的判定原则基本相同。
53.步骤s421，基于各第一区域内的图像，统计图像中全部像素点的灰度值，对灰度值进行由大到小的排布分析；步骤s422，获得灰度值排名在前的、第一百分比数量的灰度值均值；步骤s423，获得灰度值排名在后的、第二百分比数量的灰度值中位数值，第二百分比与第一百分比的加和为一；步骤s424，将灰度值均值和灰度值中位数值做差，若差值大于给定阈值，则判定该第一区域内的图像含有血管像素点，筛选获得包含血管像素点的第二区域。
54.步骤s500，将样本图像、以及样本图像中每个像素点对应的离散度信息，结合制作为输入数据集，输入数据集用于训练深度神经网络。
55.本实施例中用于训练深度神经网络的输入数据集，不仅包括样本图像，还包括离散度信息，能够使深度神经网络在降采样过程中额外的关注离散程度高的样本，进而提高深度神经网络的分割精度。
56.本实施例利用上述的样本图像进行深度学习网络训练，网络结构设计采用常见的unet架构。并对细小分支处理做了额外的设计，在深度学习网络的层中在通道信息中融合了坐标信息的离散度，该数值表示附近第二区域内（32*32大小）范围内有效像素值（血管像素点的像素值）的离散程度，有效像素值的定义与判断每帧dsa图像是否存在血管的判定准则相同，可用于在降采样过程中额外的关注离散程度高的样本。
57.步骤s600，根据样本图像中每个像素点的离散度信息，分区计算损失函数。具体包括：步骤s610，获得样本图像中每个像素点对应的离散度信息，获得统计区间；步骤s620，根据离散度信息的分布，将统计区间按划分为若干个子区间；步骤s630，针对每个子区间，分别计算深度神经网络的损失函数，以进行深度神经网络的训练。
58.本实施例中，样本图像由于每个像素点都具有一个离散度信息，则样本图像的不同像素点整体形成离散度信息的渐变。本实施例通过对不同离散度信息进行子区间划分，
并分别针对不同的的子区间，分区计算损失函数，能够使得深度神经网络能够关注到细小血管。同时，也能够进一步提高抗噪声和干扰能力，对主要血管及细小分支都有较好的分割效果。
59.应该理解的是，虽然图1-3的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图1-3中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
60.在一个实施例中，提供了一种基于深度神经网络的颅内血管树分割装置，可采用计算机设备，该计算机设备可以是终端，其内部结构图可以如图4所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种基于深度神经网络的颅内血管树分割方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
61.在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：基于深度神经网络的颅内血管树分割方法，包括利用训练完成的深度神经网络从输入的dsa影像中分割出颅内血管树，训练过程包括：对含有血管的dsa图像进行采样，从采样图像中筛选获得血管密集图像，将血管密集图像划分为小血管样本图和大血管样本图；根据典型角度造影中不同尺寸血管的相对位置关系，对各血管密集图像进行铺列组合，用于形成样本图像；在样本图像中获得每个像素点的离散度信息，离散度信息用于表示预定区域内血管像素值的离散程度；将样本图像、以及样本图像中每个像素点对应的离散度信息，结合制作为输入数据集，输入数据集用于训练深度神经网络。
62.在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：基于深度神经网络的颅内血管树分割方法，包括利用训练完成的深度神经网络从输入的dsa影像中分割出颅内血管树，训练过程包括：对含有血管的dsa图像进行采样，从采样图像中筛选获得血管密集图像，将血管密集图像划分为小血管样本图和大血管样本图；根据典型角度造影中不同尺寸血管的相对位置关系，对各血管密集图像进行铺列
组合，用于形成样本图像；在样本图像中获得每个像素点的离散度信息，离散度信息用于表示预定区域内血管像素值的离散程度；将样本图像、以及样本图像中每个像素点对应的离散度信息，结合制作为输入数据集，输入数据集用于训练深度神经网络。
63.在一个实施例中，提供了一种计算机程序产品，包括计算机指令，该计算机指令被处理器执行时实现以下步骤：基于深度神经网络的颅内血管树分割方法，包括利用训练完成的深度神经网络从输入的dsa影像中分割出颅内血管树，训练过程包括：对含有血管的dsa图像进行采样，从采样图像中筛选获得血管密集图像，将血管密集图像划分为小血管样本图和大血管样本图；根据典型角度造影中不同尺寸血管的相对位置关系，对各血管密集图像进行铺列组合，用于形成样本图像；在样本图像中获得每个像素点的离散度信息，离散度信息用于表示预定区域内血管像素值的离散程度；将样本图像、以及样本图像中每个像素点对应的离散度信息，结合制作为输入数据集，输入数据集用于训练深度神经网络。
64.本实施例中，计算机程序产品包括程序代码部分，以用于当计算机程序产品由一个或多个计算装置执行时，执行本技术各实施例中基于深度神经网络的颅内血管树分割方法的步骤。计算机程序产品可被存储在计算机可读记录介质上。还可经由数据网络（例如，通过ran、经由因特网和/或通过rbs）提供计算机程序产品以便下载。备选地或附加地，该方法可被编码在现场可编程门阵列（fpga）和/或专用集成电路（asic）中，或者功能性可借助于硬件描述语言被提供以便下载。
65.本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器（rom）、可编程rom（prom）、电可编程rom（eprom）、电可擦除可编程rom（eeprom）或闪存。易失性存储器可包括随机存取存储器（ram）或者外部高速缓冲存储器。作为说明而非局限，ram以多种形式可得，诸如静态ram（sram）、动态ram（dram）、同步dram（sdram）、双数据率sdram（ddrsdram）、增强型sdram（esdram）、同步链路（synchlink） dram（sldram）、存储器总线（rambus）直接ram（rdram）、直接存储器总线动态ram（drdram）、以及存储器总线动态ram（rdram）等。
66.以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。不同实施例中的技术特征体现在同一附图中时，可视为该附图也同时披露了所涉及的各个实施例的组合例。
67.以上所述实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来
说，在不脱离本技术构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。因此，本技术专利的保护范围应以所附权利要求为准。

技术特征：
1.基于深度神经网络的颅内血管树分割方法，其特征在于，包括利用训练完成的深度神经网络从输入的dsa影像中分割出颅内血管树，训练过程包括：对含有血管的dsa图像进行采样，从采样图像中筛选获得血管密集图像，将所述血管密集图像划分为小血管样本图和大血管样本图；根据典型角度造影中不同尺寸血管的相对位置关系，对各所述血管密集图像进行铺列组合，用于形成样本图像；在所述样本图像中获得每个像素点的离散度信息，所述离散度信息用于表示预定区域内血管像素值的离散程度；将所述样本图像、以及所述样本图像中每个像素点对应的所述离散度信息，结合制作为输入数据集，所述输入数据集用于训练所述深度神经网络。2.如权利要求1所述的颅内血管树分割方法，其特征在于，所述采样图像的滑动步长，为所述采样图像尺寸长或宽的一半；颅内血管树分割方法包括：对不含血管的dsa图像进行采样，筛选获得含有噪声分布区域的负样本图，所述负样本图用于补充铺列、以形成所述样本图像。3.如权利要求2所述的颅内血管树分割方法，其特征在于，筛选获得含有噪声分布区域的负样本图，具体包括：基于不含血管的dsa图像，对所有像素点的灰度值进行由大到小的排布分析，获得灰度值排名在后的、10%数量的像素点，作为待选位置；在所述待选位置中，选择负样本像素点，进而获得以所述负样本像素点为中心的负样本图。4.如权利要求1所述的颅内血管树分割方法，其特征在于，所述颅内血管树分割方法包括将dsa序列帧图像划分为含有血管的dsa图像、以及不含血管的dsa图像，具体包括：基于dsa序列帧图像中的每帧dsa图像，统计整幅图像中全部像素点的灰度值，对灰度值进行由大到小的排布分析；获得灰度值排名在前的、第一百分比数量的灰度值均值；获得灰度值排名在后的、第二百分比数量的灰度值中位数值，所述第二百分比与第一百分比的加和为一；将所述灰度值均值和所述灰度值中位数值做差，若差值大于给定阈值，则判定该帧dsa图像为含有血管的dsa图像，若差值小于或等于给定阈值，则判定该帧dsa图像为不含血管的dsa图像。5.如权利要求4所述的颅内血管树分割方法，其特征在于，所述采样图像大小为32*32，所述血管密集图像的判断标准为采样图像内血管像素点的占比率大于预期值，所述小血管样本图和大血管样本图依据所述血管密集图像中血管像素点的离散度划分；所述第一百分比为10%，所述第二百分比为90%。6.如权利要求1所述的颅内血管树分割方法，其特征在于，在所述样本图像中获得每个像素点的离散度信息，所述离散度信息用于表示预定区域内血管像素值的离散程度，具体包括：针对每个像素点，获得以该像素点为中心的第一区域；在各所述第一区域中，筛选获得包含血管像素点的第二区域；
将所述第二区域中血管像素点的离散度信息，作为该像素点的离散度信息。7.如权利要求6所述的颅内血管树分割方法，其特征在于，在各所述第一区域中，筛选获得包含血管像素点的第二区域，具体包括：基于各所述第一区域内的图像，统计图像中全部像素点的灰度值，对灰度值进行由大到小的排布分析；获得灰度值排名在前的、第一百分比数量的灰度值均值；获得灰度值排名在后的、第二百分比数量的灰度值中位数值，所述第二百分比与第一百分比的加和为一；将所述灰度值均值和所述灰度值中位数值做差，若差值大于给定阈值，则判定该第一区域内的图像含有血管像素点，筛选获得包含血管像素点的第二区域。8.如权利要求1所述的颅内血管树分割方法，其特征在于，所述深度神经网络为生成对抗网络，所述生成对抗网络用于对所述样本图像进行相似性判别和重构，通过相似性判别的样本图像用于制作输入数据集。9.如权利要求1所述的颅内血管树分割方法，其特征在于，所述训练过程包括：获得所述样本图像中每个像素点对应的所述离散度信息，获得统计区间；根据所述离散度信息的分布，将所述统计区间按划分为若干个子区间；针对每个子区间，分别计算深度神经网络的损失函数，以进行所述深度神经网络的训练。10.基于深度神经网络的颅内血管树分割装置，包括存储器、处理器及存储在存储器上的计算机程序，所述处理器执行所述计算机程序以实现如权利要求1～9任一项基于深度神经网络的颅内血管树分割方法的步骤。

技术总结
本申请公开了一种基于深度神经网络的颅内血管树分割方法和分割装置，颅内血管树分割方法，包括利用训练完成的深度神经网络从输入的DSA影像中分割出颅内血管树，训练过程包括：对含有血管的DSA图像进行采样，从采样图像中筛选获得血管密集图像，划分为小血管样本图和大血管样本图；根据典型角度造影中不同尺寸血管的相对位置关系，对各所述血管密集图像进行铺列组合，用于形成样本图像；在所述样本图像中获得每个像素点的离散度信息，所述离散度信息用于表示预定区域内血管像素值的离散程度；将所述样本图像、以及所述样本图像中每个像素点对应的所述离散度信息，结合制作为输入数据集，所述输入数据集用于训练所述深度神经网络。络。络。


技术研发人员：冷晓畅 陆徐洲 鲁伟 何京松 向建平
受保护的技术使用者：杭州脉流科技有限公司
技术研发日：2023.08.03
技术公布日：2023/9/9