[发明专利]一种在类球体区域约束下进行气管壁缺口封堵的肺气道树分割方法

权利要求书

1.一种在类球体区域约束下进行气管壁缺口封堵的肺气道树分割方法，其特征在于包括以下步骤：

步骤1：读入胸部CT图像，在主气管中设定初始种子点；

步骤2：检测连通域并标记，将种子点留存矩阵中的单连通种子点集合依次压入堆栈，该堆栈用于明确各种子点集合独立分割的顺序；

步骤3：依次提取单连通种子点集合，建立类球体区域约束模型；

步骤4：依次对建立的类球体区域进行分割，进行泄露检测并调整参数；

步骤5：依次定位各类球体区域内的气管壁低亮度缺口，提高原图像该位置体素点的CT值以封堵气管壁缺口，重新分割该区域内图像；

步骤6：依次保存分割结果，依次保存区域边缘分割结果至种子点留存矩阵P；

步骤7：若种子点留存矩阵P为空，则分割结束，否则返回步骤2；

所述步骤2中连通域检测和标记的方法如下：

P_Label(P)＝[P_area(k),k],k＝1,2,3,...,n.

其中，P为种子点留存矩阵，k为堆栈中的排序值和标记值，n为该堆栈中k的最大值，P_area(k)为标记k的单连通种子点集合，满足如下性质：

(1)

(2)

(3) max[P_area(k)]＝1；

(4)

所述步骤3中，建立类球体区域约束模型的过程如下：

(1)计算单连通种子点集合体素点总数，并计算区域半径R：

其中，f_area(x,y,z,k)为种子点集合的二值图像函数，ceil为向上取整操作；

(2)计算初始阈值：

分割中的阈值用Threshold(k)表示，初始阈值为种子点集合对应原图像位置体素点的CT值变化范围；

(3)数据切割以优化算法：

首先寻找种子点曲面边界，然后将边界扩张R(k)个体素点距离，且不超出CT图像边界；

(4)实现约束区域及区域边界的生成：

以R(k)为半径的三维球形算子Sphere，对种子点集合曲面进行单次膨胀操作，可得近似球体的约束区域T_region(k)，区域外体素的CT值取最低值：-2048HU，

用Imdilate表示形态学膨胀操作，则约束区域边缘为：

T_edge(k)＝Imdilate[R(k)]-Imdilate[R(k)-1]。

2.根据权利要求1所述的肺气道树分割方法，其特征在于，所述步骤4中，区域内分割采用波面传播法提取待分割点，并使用阈值法对待分割点进行判定。

3.根据权利要求2所述的肺气道树分割方法，其特征在于，分割过程即迭代过程，迭代次数和每一次迭代的结果均可参与泄漏判定，其中，迭代次数受区域半径R(k)限制，每一次迭代结果的体素点总量受种子点集合内元素总量限制，区域半径R(k)限制和种子点集合内元素总量限制的参数设置保留弹性；若判定泄漏发生，则该小段分割返回初始状态并分割阈值衰减，为保持分割精细程度，设定Threshold(k)-1；当迭代结果不再增加时，区域内分割结束，分割结果表示为X(k)；此时提取区域边缘上分割结果，计算其面积以判定是否泄漏，此外该小段分割结果的体积也可作为泄漏判定的条件，区域边缘面积限定法和区域内分割结果体积限定法同样参考种子点集合内元素总量和区域半径R(k)。

4.根据权利要求1所述的肺气道树分割方法，其特征在于，所述步骤5中，对步骤4的分割结果X(k)进行形态学闭操作，其结果减去X(k)可得存在低亮度气管壁缺口的气管壁图像，再对该图像进行形态学闭操作，其结果减去该图像，可得准确的气管壁缺口位置。

5.根据权利要求4所述的肺气道树分割方法，其特征在于，设定缺口封堵开关，重复气管壁缺口封堵操作，在同一小段气道上实现多个气管壁缺口的封堵，当缺口封堵检测时区域内分割阈值不再增加或超出细小气道阈值范围Threshold(k)＝200时，关闭缺口封堵开关，跳出循环，封堵结束。