发明内容
本申请的主要目的在于提供一种气管树提取方法及装置,以解决相关技术中在基于肺部CT图像提取气管树的精度较低的问题。
根据本申请的一个方面,提供了一种气管树提取方法,包括:从肺部计算机体层摄影CT图像中识别支气管区域;基于所述支气管区域获取第一阈值;在所述支气管区域根据种子点以及第一阈值进行区域增长;判断是否出现泄露;在出现泄露的情况下,对出现泄露的部分进行填充;基于填充后的所述支气管区域获取所述第一阈值,并使用获取到的所述第一阈值重新在填充后的所述支气管区域进行区域增长,直至不再出现泄露或进行区域增长的次数达到预设次数时,获得气管树。
可选地,所述从肺部CT图像中识别支气管区域,包括:对所述肺部CT图像进行检测,得到空气组织区域以及人体组织区域;对所述人体组织区域进行检测,得到肺部区域;根据所述肺部区域与所述空气组织区域得到支气管区域。
可选地,所述基于所述支气管区域获取第一阈值,或基于填充后的所述支气管区域获取所述第一阈值,包括:计算所述肺部区域的像素中大于初始阈值的像素的第一平均值以及小于所述初始阈值的像素的第二平均值;计算所述第一平均值与所述第二平均值之间的差值的绝对值;在所述绝对值小于第三阈值的情况下,确定所述第一平均值与所述第二平均值之间的差值为进行区域增长的阈值。
可选地,所述方法还包括:在对出现泄露的部分进行填充之后,以填充后的支气管区域的边缘作为初始轮廓线进行活动轮廓处理。
可选地,所述判断是否出现泄露,包括:在每次进行区域增长后,计算得到的气管树的体积,如果相邻两次计算得到的体积的差值大于预设数值,则确定出现泄露。
根据本申请的另一个方面,提供了一种气管树提取装置,包括:识别模块,用于从肺部计算机体层摄影CT图像中识别支气管区域;获取模块,用于基于所述支气管区域获取第一阈值;增长模块,用于在所述支气管区域根据种子点以及第一阈值进行区域增长;判断模块,用于判断是否出现泄露;填充模块,用于在出现泄露的情况下,对出现泄露的部分进行填充;获得模块,用于基于填充后的所述支气管区域获取所述第一阈值,并使用获取到的所述第一阈值重新在填充后的所述支气管区域进行区域增长,直至不再出现泄露时或进行区域增长的次数达到预设次数,获得气管树。
可选地,所述识别模块,包括:第一检测单元,用于对所述肺部CT图像进行检测,得到空气组织区域以及人体组织区域;第二检测单元,用于对所述人体组织区域进行检测,得到肺部区域;获得单元,用于根据所述肺部区域与所述空气组织区域得到支气管区域。
可选地,所述获取模块,以及所述获得模块均包括如下单元:第一计算单元,用于计算所述肺部区域的像素中大于初始阈值的像素的第一平均值以及小于所述初始阈值的像素的第二平均值;第二计算单元,用于计算所述第一平均值与所述第二平均值之间的差值的绝对值;确定单元,用于在所述绝对值小于第三阈值的情况下,确定所述第一平均值与所述第二平均值的之间的差值作为进行区域增长的阈值。
可选地,所述装置还包括:处理模块,用于在对出现泄露的部分进行填充之后,以填充后的支气管区域的边缘作为初始轮廓线进行活动轮廓处理。
可选地,所述判断模块包括:判断单元,用于在每次进行区域增长后,计算得到的气管树的体积,如果相邻两次计算得到的体积的差值大于预设数值,则确定出现泄露。
本申请的实施例提供的方案,基于第一阈值在支气管区域进行区域增长,找到泄露部分,在对泄露部分进行填充后,继续以第一阈值在支气管区域进行区域增长,直至不再出现泄露时,获得气管树,该方案通过对泄露进行填充,以及在填充后继续进行区域增长,可以尽可能的获取较多的气管树,从而实现在未对肺部CT图像进行预处理的基础上也可以保证气管树分割的精确程度。
具体实施方式
下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例,都属于本申请保护的范围。
本申请的实施例提供了一种气管树提取方法,图1是该方法的流程图,如图1所示,该方法包括如下处理:
步骤101:从肺部计算机体层摄影CT图像中识别支气管区域;
在本实施例中,从肺部CT图像中识别支气管区域的操作可以包括:对肺部CT图像进行检测,得到空气组织区域以及人体组织区域;对得到的人体组织区域进行检测,得到肺部区域;根据肺部区域与空气组织区域得到支气管区域。
步骤102:基于支气管区域获取第一阈值;
在本实施例中,第一阈值即为用于进行区域增长的阈值,图2示出了在本实施例中获取第一阈值的流程,如图2所示,基于支气管区域获取第一阈值,或基于填充后的支气管区域获取第一阈值的操作可以包括:
设置初始阈值T,示例的,该阈值可以设置为-1000,计算肺部区域的像素中大于初始阈值的像素的第一平均值T1以及小于初始阈值的像素的第二平均值T2;计算第一平均值T1与第二平均值T2的之间的差值的绝对值,即abs(T1+T2)/2;在计算出的绝对值小于第三阈值eps(例如10∧(-7))的情况下,确定第一平均值与第二平均值的差值作为进行区域增长的阈值。
步骤103:在支气管区域根据种子点以及第一阈值进行区域增长;
其中,种子点可以预先设置,也可以由计算机自动选择。
步骤104:判断是否出现泄露;
在本实施例中,判断是否出现泄露的操作具体可以包括:在每次进行区域增长后,计算得到的气管树的体积,如果相邻两次计算得到的体积的差值大于预设数值,则确定出现泄露,示例的,该预设数值可以是200,单位为立方厘米,即如果相邻两次计算得到的体积的差值大于200立方厘米,则确定出现泄露。
步骤105:在出现泄露的情况下,对出现泄露的部分进行填充;
对出现泄露的部分进行填充可以使得下一次进行区域增长时,该泄露部分不会被识别出来,。
步骤106:基于填充后的支气管区域获取第一阈值,并使用获取到的第一阈值重新在填充后的支气管区域进行区域增长,直至不再出现泄露时或进行区域增长的次数达到预设次数,获得气管树。
可选地,本实施例的气管树提取方法还包括:在对出现泄露的部分进行填充之后,以填充后的支气管区域的边缘作为初始轮廓线进行活动轮廓处理,以对进行区域增长分割得到的区域进行改进。
图3示出了本申请的实施例的另一种气管树提取方法,如图3所示,该方法包括如下处理:
获取肺部CT数据;
对肺部CT数据进行初始分割,得到空气组织区域以及人体组织区域;
对得到的人体组织区域进行检测,得到肺部区域;
对得到的空气组织区域进行区域提取,得到肺部区域组织中的空气区域,即支气管区域;
对支气管区域进行支气管提取,得到支气管树,以及初始肺区域;
对支气管树进行检测,得到肺叶和分割后的支气管树;
对初始肺区域进行检侧,得到左肺区域以及右肺区域。
图4示出了本实施例在进行气管树提取过程中的支气管提取算法的流程,如图4所示,该流程包括:
基于肺区域中的空气区域,即支气管区域进行种子点检测,得到种子点;
使用阈值th对支气管区域进行三维区域增长;
在进行三维区域增长的过程中检测是否出现泄露;
在出现泄露的情况下,检测泄露区域;
对泄露区域进行标记,以使该泄露区域不参与后续进行的区域增长过程;
获取新的阈值th’=th+1,基于标记后的支气管区域重新进行区域增长,直至进行区域增长的次数达到预设次数N,得到分割后的气管树。
本申请的实施例还提供了一种气管树提取装置,该装置用于实现本申请的气管数据提取方法,图5是该装置的框图,如图5所示,该装置50包括如下组成部分:
识别模块51,用于从肺部CT图像中识别支气管区域;
获取模块52,用于基于支气管区域获取第一阈值;
增长模块53用于在支气管区域根据种子点以及第一阈值进行区域增长;
判断模块54,用于判断是否出现泄露;
填充模块55,用于在出现泄露的情况下,对出现泄露的部分进行填充;
获得模块56,用于基于填充后的支气管区域获取第一阈值,并使用获取到的第一阈值重新在填充后的支气管区域进行区域增长,直至不再出现泄露时或进行区域增长的次数达到预设次数,获得气管树。
可选地,上述识别模块51进一步可以包括:第一检测单元,用于对肺部CT图像进行检测,得到空气组织区域以及人体组织区域;第二检测单元,用于对人体组织区域进行检测,得到肺部区域;获得单元,用于根据肺部区域与空气组织区域得到支气管区域。
可选地,上述获取模块52,以及获得模块56均包括如下单元:第一计算单元,用于计算肺部区域的像素中大于初始阈值的像素的第一平均值以及小于初始阈值的像素的第二平均值;第二计算单元,用于计算第一平均值与第二平均值之间的差值的绝对值;确定单元,用于在绝对值小于第三阈值的情况下,确定第一平均值与第二平均值的之间的差值作为进行区域增长的阈值。
可选地,本申请的气管树提取装置还可以进一步包括:处理模块,用于在对出现泄露的部分进行填充之后,以填充后的支气管区域的边缘作为初始轮廓线进行活动轮廓处理。
可选地,上述判断模块54具体可以包括:判断单元,该单元用于在每次进行区域增长后,计算得到的气管树的体积,如果相邻两次计算得到的体积的差值大于预设数值,则确定出现泄露。
本申请的实施例提供的方案,基于第一阈值在支气管区域进行区域增长,找到泄露部分,在对泄露部分进行填充后,重新以第一阈值在支气管区域进行区域增长,可以尽可能的获取较多的气管树,从而实现在未对肺部CT图像进行预处理的基础上也可以保证气管树分割的精确程度。此外,通过活动轮廓法对进行区域增长得到区域进行改进,也进一步提高了气管树分割的效果。
以上所述仅为本申请的实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。