[发明专利]一种从胸部CT图像中提取肺血管的方法及装置

权利要求书

1.一种从胸部CT图像中提取肺血管的方法，其特征在于，包括：

接收输入的n层胸部CT图像，其中n为自然数，获取指定图像层，在所述指定图像层上进行灰度映射获得对应的灰度直方图，根据所述灰度直方图获得预设分割阈值；

在所述指定图像层上选取肺部区域的指定像素点获得初始种子点，根据所述预设分割阈值以及初始种子点，在所述n层胸部CT图像进行3D区域增长，获得不带血管的肺组织区域；

在所述的肺组织区域进行形态学运算，获得带血管的封闭肺组织区域，计算血管阈值；

在所述指定图像层上的带血管的封闭肺组织区域查找大于血管阈值的像素点为初始标记点，获得初始扩散面；

在所述初始扩散面和血管阈值，在n层胸部CT图像的带血管的封闭肺组织区域内进行面扩散获得肺血管计算步骤包括：

01)所述n层胸部CT图像的带血管的封闭肺组织区域的像素点将归为四类:第一类为当前扩散点，第二类为边界点，第三类为待处理完点，第四类为已处理完点，初始均设置为第三类点；

02)把在所述指定图像层上带血管的肺组织区域标记为肺血管的像素点标记为第一类点；

03)将与第一类点相邻区域中的第三类像素点标记为第二类点；

04)遍历计算第一类点与其相邻第二类点的到达时间T,并将第一类点标记为第四类点，计算第二类像素点(i,j,k)的到达时间T(i,j,k)的步骤为：其中D_i、D_j、D_k为像素点(i,j,k)与X、Y、Z三个方向相邻像素点的灰度值差绝对值的最大值，具体为：D_i＝max(|V_i,j,k-V_i-1,j,k|,|V_i,j,k-V_i+1,j,k|),D_j＝max(|V_i,j,k-V_i,j-1,k|,|V_i,j,k-V_i,j+1,k|),D_k＝max(|V_i,j,k-V_i,j,k-1|,|V_i,j,k-V_i,j,k+1|)，其中V_i,j,k为像素点(i,j,k)的灰度值；

05)设置到达时间T大于预设阈值的所述第二类点的到达时间为0,并标记为第四类点，设置到达时间T不大于预设阈值的所述第二类点为第一类点；

06)判断与第一类点相邻区域是否存在第三类点，如存在，返回步骤03)，否则将第一类点设置为第四类点，执行步骤07)；

07)获取所述n层胸部CT图像的带血管的封闭肺组织区域的像素点的到达时间T大于0的像素点，即为肺血管。

2.根据权利要求1所述的方法，其特征在于，在所述指定图像层上选取肺部区域的指定像素点获得初始种子点，根据所述预设分割阈值以及初始种子点，在所述n层胸部CT图像进行3D区域增长,获得不带血管的肺组织区域的步骤包括：

01)在所述指定图像层上选取肺部区域的指定像素点获得初始种子点；

02)从所述初始被标记种子点的26邻域像素点中选取一个像素点；

03)判断所选取的像素点是否已经被标记为标记点，若是，则返回02),否则执行步骤04)；

04)判断所选取的像素点的灰度值是否满足预设要求，若是，则把该像素点标记为标记点,加入标记点集，执行步骤06),否则执行步骤05)；

05)停止标记该像素点，执行步骤06)；

06)判断所述26邻域像素点是否全部判断完毕，若是则执行步骤08),否则返回步骤02)；

08)判断标记点集是否为空，如不为空，则从标记点集中取出一个标记点作为初始被标记种子点，返回步骤02),同时将该点从标记点集中去除,否则执行步骤09)；

09)获取已被标记的像素点集,该像素点集即为不带血管的肺组织区域。