[发明专利]一种血管分叉级别判断方法、装置及存储介质

权利要求书

1.一种血管分叉级别判断方法，其特征在于，所述方法包括：

对血管原始图像进行预处理，得到预处理后的血管图像；

根据所述预处理后的血管图像的所有体素点获取血管骨架；

将所述血管骨架中的所有骨架点利用图算法生成最小生成树；

以最小生成树的主干根部为起始点，赋予起始点一个标签值，从起始点开始，以深度优先搜索算法遍历最小生成树中的所有节点，当遍历到分叉点时，分叉点后续的节点改变标签值，直到完成对最小生成树所有的节点的标签值赋值；

根据最小生成树中的节点的标签值得到血管骨架的骨架点的标签值，获取与所述血管原始图像中每个体素点在血管骨架中最近的骨架点，将该骨架点的标签值赋予给对应的体素点，完成对血管原始图像的血管分叉分级。

2.根据权利要求1所述的方法，其特征在于，

所述对血管原始图像进行预处理，得到预处理后的血管图像包括：

提取血管原始图像的最大连通域，对最大连通域血管标签图像在各方向上进行等距重采样；

对等距重采样后的血管原始图像进行表面光滑处理，得到预处理后的血管图像。

3.根据权利要求1所述的方法，其特征在于，

所述将所述血管骨架中的所有骨架点利用图算法生成最小生成树，包括：

获取血管骨架所有骨架点，根据任意两个骨架点之间的空间距离生成血管骨架的邻接矩阵；

根据血管骨架的邻接矩阵生成无向加权图；

获取所述无向加权图的最小生成树。

4.根据权利要求3所述的方法，其特征在于，

所述获取骨架所有的骨架点，根据任意两个骨架点之间的空间距离生成血管骨架的邻接矩阵包括：

获取血管骨架中包含的N个骨架点，每个骨架点均包括其空间坐标，生成每一个骨架点之间的距离组成的N阶方阵，将该方阵作为邻接矩阵。

5.根据权利要求3所述的方法，其特征在于，

所述根据血管骨架的邻接矩阵生成无向加权图包括：

将血管骨架中的每一个骨架点均与其他所有的骨架点连接，获得全连接图；

选取全连接图中任意两个连接的骨架点，根据邻接矩阵中该任意两个骨架点的距离对该任意两个连接的骨架点的连接强弱进行加权，得到无向加权图。

6.根据权利要求1所述的方法，其特征在于，

所述以最小生成树的主干根部为起始点，赋予起始点一个标签值，从起始点开始，以深度优先搜索算法遍历最小生成树中的所有节点，当遍历到分叉点时，分叉点后续的节点改变标签值，直到完成对最小生成树所有的节点的标签值赋值包括：

以初始点为起点，获取每一个节点的度，所述节点的度为与一个节点直接连接的节点的个数之和，自初始点后，若遍历到的每一个节点的度均小于或等于2，则表示该节点不是分叉点，对该节点赋予初始点的标签值，直到遇到一个节点，其度大于或等于3，表示该节点为分叉点，对于分叉点，同样赋予初始点的标签值，而分叉点之后的节点则赋予新的标签值，直到遇到新的分叉点，对新的分叉点之后的节点再次赋予新的标签值，直到完成所有节点的赋值。

7.一种血管分叉级别判断装置，其特征在于，所述装置包括：

标签图像获取模块：用于对血管原始图像进行预处理，得到预处理后的血管图像；

骨架获取模块：用于根据所述预处理后的血管图像的所有体素点获取血管骨架；

最小生成树生成模块：用于将所述血管骨架中的所有骨架点利用图算法生成最小生成树；

赋值模块：用于以最小生成树的主干根部为起始点，赋予起始点一个标签值，从起始点开始，以深度优先搜索算法遍历最小生成树中的所有节点，当遍历到分叉点时，分叉点后续的节点改变标签值，直到完成对最小生成树所有的节点的标签值赋值；

血管分叉级别判断模块：用于根据最小生成树中的节点的标签值得到血管骨架的骨架点的标签值，获取与所述血管原始图像中每个体素点在血管骨架中最近的骨架点，将该骨架点的标签值赋予给对应的体素点，完成对血管原始图像的血管分叉分级。

8.一种存储介质，其特征在于，所述存储介质存储有计算机程序，所述计算机程序被主控器执行时，实现如权利要求1-6任一项所述的一种血管分叉级别判断方法中的各个步骤。