[发明专利]一种医学图像中肺分割的方法及装置

权利要求书

1.一种医学图像中肺分割的方法,其特征在于包括以下步骤:

根据训练集中的M张胸片的肺轮廓，得到平均肺模板，其中M为大于或等于2的整数；

对拟被分割的肺部图像预处理，获得预处理后的肺部图像；

提取预处理后的肺部图像的肺边界的二值图像，并根据广义霍夫变换进行初始定位，获得相对应的霍夫定位位置；

将平均肺模板与霍夫定位位置对齐，获得对齐结果；

所述平均肺模板按照以下步骤获得：

a.从训练集中的M张胸片中进行第一次选取，获取第一肺区形状L₀；

b.从训练集中的尚未被选取的胸片中进行第二次选取，获取N张肺区形状L₁、L₂…L_N；

c.使所述N张肺区形状L₁、L₂…L_N经处理后得到相应的N张肺区形状L₁₁、L₂₁…L_N1，分别与第一肺区形状L₀对齐，计算得到肺区平均形状L_ave，所述处理包括使N张肺区形状L₁、L₂…L_N分别相对于第一肺区形状L₀作旋转、缩放、平移变换；

d.若c步骤的肺区平均形状L_ave相较于第一肺区形状L₀的旋转、缩放、平移小于设定阈值，则将所述肺区平均形状L_ave作为平均肺模板；否则，进入e步骤；

e.将第一肺区形状L₀进行处理、所述N张肺区形状L₁₁、L₂₁…L_N1进行处理后，分别与当前肺区平均形状L_ave-i对齐，计算得到本次肺区平均形状L_ave1-i，所述对第一肺区形状处理包括：使第一肺区的形状L₀分别相对于当前肺区平均形状L_ave-i作旋转、缩放、平移变换，得到相应的肺区形状L_0i，以及使所述N张肺区形状L₁₁、L₂₁…L_N1分别相对于当前肺区平均形状L_ave-i作旋转、缩放、平移变换，得到N张肺区形状L_1i、L_2i…L_Ni；

f.若e步骤计算所得本次肺区平均形状L_ave1-i与当前肺区平均形状L_ave-i的旋转、缩放、平移小于设定阈值，则将e步骤计算所得本次肺区平均形状L_ave1-i作为平均肺模板；否则，将当前肺区平均形状L_ave-i用本次肺区平均形状L_ave1-i替代，将第一肺区的形状L₀用肺区形状L_0i替代，N张第二肺区形状L₁₁、L₂₁…L_N1用N张肺区形状L_1i、L_2i…L_Ni替代，且i＝i+1，其中i大于或等于2的整数，“＝”为赋值号，并返回步骤e；

利用动态规划算法进行分割，并将分割结果反变换回原始坐标系中，完成肺区域分割。

2.根据权利要求1所述的医学图像中肺分割的方法，其特征在于，在M张胸片中所对应的每张胸片的肺轮廓上标记相同数量的特征点，且不同胸片上的特征点的分布一致。

3.根据权利要求2所述的医学图像中肺分割的方法，其特征在于，所述在每张胸片的肺轮廓的拐点处、最高点处及最低点处标记特征点。

4.根据权利要求3所述的医学图像中肺分割的方法，其特征在于，所述在每张胸片的肺轮廓上除拐点处、最高点处及最低点处之外的其他位置等距离标记特征点。