基于前向波气管分割的气管树提取方法
【专利摘要】本发明公开了基于前向波气管分割的气管树提取方法，包括：步骤A：获取胸部CT数据；步骤B：创建气管树结构图像和临时数据池；步骤C：将预设的根结点初始化为第一结点集合；步骤D：遍历第一结点集合中每个结点的像素值的26?邻域，将26?邻域中CT值与像素值之差小于第一阈值的结点，添加到第二结点集合；步骤E：当第二结点集合的数量等于零，或者第二结点集合的数量与第一结点集合的数量之商大于等于第二阈值时，执行步骤F；步骤F：检查第一结点集合是否存在父结点；若不存在，则输出临时数据池中所有结点的像素值和CT值，获得提取后的气管树结构图像。采用本发明实施例能快速准确地提取气管树结构，分割结果精确且耗时少。
【专利说明】
基于前向波气管分割的气管树提取方法
技术领域
[0001]本发明涉及计算机图像处理技术领域，尤其涉及一种基于前向波气管分割的气管树提取方法。
【背景技术】
[0002]随着现代社会工业化的发展，空气污染加重，人体肺部的发病率越来越高，且类型繁多，如肺癌、气管炎、肺气肿、哮喘等。而随着多排螺旋CT的发展，影像科医生可以通过一次扫描获得患者的影像信息，并结合影像处理技术对疾病进行分析诊断，CT成像被认为是诊断肺部疾病的金标准。但是由于气管树的结构复杂，对比度差，且容易受到噪声、容积效应的影响，如在肺部气管内部充满空气时，通过人体组织CT值可知，肺部气管在CT图像中CT值约为-1000HU，但受噪声或容积效应影响时，肺部区域的CT值会在一个范围内波动，且肺部在吸气和呼气时候的CT值也不尽相同，那么可假设气管的灰度在-1200HU至-800HU之间，或者与种子点灰度相差200HU以内的即可视为是气管上的像素。
[0003]为了排除上述影响，在CT图像的影像处理时常使用气管分割法进行处理。目前最为常见的气管分割方法皆采用区域增长的方式，该方法中的阈值选择对分割结果影响较大，当阈值选择不恰当时，不仅会造成分割时间过长，而且容易导致分割失败，图像处理效果不甚理想。
【发明内容】
[0004]本发明实施例提出一种基于前向波气管分割的气管树提取方法，能快速准确地提取气管树结构，分割结果精确且耗时少。
[0005]本发明实施例提供一种基于前向波气管分割的气管树提取方法，包括:
[0006]步骤A:获取胸部CT数据；
[0007]步骤B:根据所述胸部CT数据，创建气管树结构图像和临时数据池;其中，所述气管树结构图像包括:若干个结点以及所述若干个结点所对应的像素值和CT值;所述临时数据池用于存储第一结点合集中所有结点的像素值和CT值；
[0008]步骤C:将预设的根结点初始化为所述第一结点集合；
[0009]步骤D:遍历所述第一结点集合中每个结点的像素值的26-邻域，将所述26-邻域中CT值与像素值之差小于第一阈值的结点，添加到第二结点集合；
[0010]步骤E:当所述第二结点集合的数量等于零，或者所述第二结点集合的数量与所述第一结点集合的数量之商大于等于第二阈值时，执行步骤F;
[0011 ]步骤F:检查所述第一结点集合是否存在父结点；若不存在，则输出所述临时数据池中所有结点的像素值和CT值，获得提取后的气管树结构图像。
[0012]进一步的，所述步骤E还包括:当所述第二结点集合与所述第一结点集合的商小于第二阈值时，执行步骤G;
[0013]步骤G:检查所述第二结点集合中各结点的像素是否具有区域连通性;若是，则将所述第二结点集合中所有结点添加到所述第一结点集合中，并清空所述第二结点集合，再返回步骤D;否则，将非连通区域中结点属性为左分支的结点添加到所述第一结点集合中，将非连通区域中结点属性为右分支的结点添加到第三结点集合中，并清空所述第二结点集合，再执行步骤H;
[0014]步骤H:检查所述第一结点集合中结点属性为左分支的结点是否均被处理;若是，则将所述第三结点集合添加到所述第一结点集合中，再执行步骤I;否则，返回步骤D;
[0015]步骤1:检查所述第一结点集合中结点属性为右分支的结点是否均被处理;若是，则返回步骤F，否则返回步骤D。
[0016]进一步的，所述步骤F还包括:若所述第一结点集合存在父结点，则将所述第一结点集合中除所述父结点外所有结点添加到第四结点集合，并将所述父结点初始化为所述第一结点集合，再执行步骤H;
[0017]其中，所述第四结点集合中所有结点的像素值和CT值存储在所述临时数据池中。
[0018]进一步的，在所述步骤A之后，步骤B之前还包括:
[0019]使用阈值法去除所述胸部CT数据的背景，仅保留胸腔及胸腔内的器官数据。
[0020]实施本发明实施例，具有如下有益效果:
[0021 ]本发明实施例提供的基于前向波气管分割的气管树提取方法，在获取胸部CT数据后，从预设的根结点开始，依次查找上一次迭代过的所有结点的像素的26-邻域中符合条件的像素，并一直朝着肺门方向推进，直到气管末梢再退回到上一级分支处向另一支气管末梢推进，同一条分支上只限单方向增长，直到找到所有的支气管，获得提取后的气管树结构图像。相比于现有技术使用区域增长的方法提取气管树，本发明无需使用阈值来区分气管像素，而是通过考虑待判定像素与相邻的气管像素的相似性来判定该像素是否属于气管，使得分割结果更具有连续性及鲁棒性。
【附图说明】
[0022]图1是本发明提供的基于前向波气管分割的气管树提取方法的一种实施例的流程示意图；
[0023]图2是本发明提供的基于前向波气管分割的气管树提取方法的另一种实施例的流程不意图。
【具体实施方式】
[0024]下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0025]参见图1，是本发明提供的基于前向波气管分割的气管树提取方法的一种实施例的流程示意图，该方法包括以下步骤A至步骤F，具体如下:
[0026]步骤A:获取胸部CT数据。
[0027]在本实施例中，在所述步骤A之后，步骤B之前还包括:使用阈值法去除所述胸部CT数据的背景，仅保留胸腔及胸腔内的器官数据。
[0028]步骤B:根据胸部CT数据，创建气管树结构图像和临时数据池;其中，气管树结构图像包括:若干个结点以及若干个结点所对应的像素值和CT值；临时数据池用于存储第一结点合集中所有结点的像素值和CT值。
[0029]在本实施例中，临时数据池用于存储每次迭代过程中找到的结点的像素值和CT值，或者只存储该找到的结点，其像素值和CT值可通过关联关系查询相应的CT数据而获得。
[0030]步骤C:将预设的根结点初始化为第一结点集合。
[0031]在本实施例中，该预设的根结点可以是横断面最顶层的气管组织上的种子结点。[0032 ]步骤D:遍历第一结点集合中每个结点的像素值的26-邻域，将26-邻域中CT值与像素值之差小于第一阈值的结点，添加到第二结点集合。
[0033]在本实施例中，将26-邻域中CT值与像素值之差小于第一阈值的结点，添加到第二结点集合，具体为将灰度相近的结点添加到第二结点集合。第一阈值的优选值为200HU。
[0034]步骤E:当第二结点集合的数量等于零，或者第二结点集合的数量与第一结点集合的数量之商大于等于第二阈值时，执行步骤F。
[0035]在本实施例中，第二阈值的优选值为1.4。
[0036]步骤F:检查第一结点集合是否存在父结点；若不存在，则输出临时数据池中所有结点的像素值和CT值，获得提取后的气管树结构图像。
[0037]作为本实施例的一种举例，参见图2，图2是本发明提供的基于前向波气管分割的气管树提取方法的另一种实施例的流程示意图。如图2所示，步骤E为:判断第二结点集合的数量是否等于零，或者第二结点集合的数量与第一结点集合的数量之商是否大于等于第二阈值，若是，执行步骤F;当第二结点集合与第一结点集合的商小于第二阈值时，执行步骤G。
[0038]步骤G:检查第二结点集合中各结点的像素是否具有区域连通性;若是，则将第二结点集合中所有结点添加到第一结点集合中，并清空第二结点集合，再返回步骤D;否则，将非连通区域中结点属性为左分支的结点添加到第一结点集合中，将非连通区域中结点属性为右分支的结点添加到第三结点集合中，并清空所述第二结点集合，再执行步骤H;
[0039]步骤H:检查第一结点集合中结点属性为左分支的结点是否均被处理;若是，则将第三结点集合添加到第一结点集合中，执行步骤I;否则，返回步骤D;
[0040]步骤1:检查第一结点集合中结点属性为右分支的结点是否均被处理;若是，则返回步骤F，否则返回步骤D。
[0041]在本举例中，步骤F还包括:若第一结点集合存在父结点，则将第一结点集合中除父结点外所有结点添加到第四结点集合，并将父结点初始化为第一结点集合，再执行步骤H;其中，第四结点集合中所有结点的像素值和CT值存储在临时数据池中。
[0042]在本实施例中，采用树形结构记录气管分支形状，主要是利用气管的生物学形态结构和树的数据结构完全一致，因此采用二叉树结构来记录气管的分支情况具有天然的优势。另外采用深度优先方式进行遍历，可在运算过程中实时跟踪某一分支的处理情况，不会发生重复运算，有助于运算速度的提升。
[0043]在本实施例中，本发明利用气管进入肺门的解剖学信息作为提取的停止条件。相比于现有的区域生长分割方法因为气管与肺实质的灰度几乎一致，使得提取的结果包含整个胸部体数据中的呼吸系统。而本发明考虑到气管进入肺实质处(即肺门)时截面积呈迅速扩张趋势，以及气管上可能存在肿瘤、局部扩张等因素，若相邻截面面积之超过1.4，即可认为气管进入肺实质，并使用这一特点判断气管末梢的位置，与气管的实际解剖学结构相符入口 O
[0044]综上所述，本发明实施例提供的基于前向波气管分割的气管树提取方法，在获取胸部CT数据后，从预设的根结点开始，依次查找上一次迭代过的所有结点的像素的26-邻域中符合条件的像素，并一直朝着肺门方向推进，直到气管末梢再退回到上一级分支处向另一支气管末梢推进，同一条分支上只限单方向增长，直到找到所有的支气管，获得提取后的气管树结构图像。相比于现有技术使用区域增长的方法提取气管树，本发明无需使用阈值来区分气管像素，而是通过考虑待判定像素与相邻的气管像素的相似性来判定该像素是否属于气管，使得分割结果更具有连续性及鲁棒性。
[0045]本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存储记忆体(Random AccessMemory，RAM)等。
[0046]以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。
【主权项】
1.一种基于前向波气管分割的气管树提取方法，其特征在于，包括: 步骤A:获取胸部CT数据； 步骤B:根据所述胸部CT数据，创建气管树结构图像和临时数据池;其中，所述气管树结构图像包括:若干个结点以及所述若干个结点所对应的像素值和CT值;所述临时数据池用于存储第一结点合集中所有结点的像素值和CT值； 步骤C:将预设的根结点初始化为所述第一结点集合； 步骤D:遍历所述第一结点集合中每个结点的像素值的26-邻域，将所述26-邻域中CT值与像素值之差小于第一阈值的结点，添加到第二结点集合； 步骤E:当所述第二结点集合的数量等于零，或者所述第二结点集合的数量与所述第一结点集合的数量之商大于等于第二阈值时，执行步骤F; 步骤F:检查所述第一结点集合是否存在父结点；若不存在，则输出所述临时数据池中所有结点的像素值和CT值，获得提取后的气管树结构图像。2.根据权利要求1所述的基于前向波气管分割的气管树提取方法，其特征在于，所述步骤E还包括:当所述第二结点集合与所述第一结点集合的商小于第二阈值时，执行步骤G; 步骤G:检查所述第二结点集合中各结点的像素是否具有区域连通性;若是，则将所述第二结点集合中所有结点添加到所述第一结点集合中，并清空所述第二结点集合，再返回步骤D;否则，将非连通区域中结点属性为左分支的结点添加到所述第一结点集合中，将非连通区域中结点属性为右分支的结点添加到第三结点集合中，并清空所述第二结点集合，再执行步骤H; 步骤H:检查所述第一结点集合中结点属性为左分支的结点是否均被处理;若是，则将所述第三结点集合添加到所述第一结点集合中，再执行步骤I;否则，返回步骤D; 步骤1:检查所述第一结点集合中结点属性为右分支的结点是否均被处理;若是，则返回步骤F，否则返回步骤D。3.根据权利要求2所述的基于前向波气管分割的气管树提取方法，其特征在于，所述步骤F还包括:若所述第一结点集合存在父结点，则将所述第一结点集合中除所述父结点外所有结点添加到第四结点集合，并将所述父结点初始化为所述第一结点集合，再执行步骤H; 其中，所述第四结点集合中所有结点的像素值和CT值存储在所述临时数据池中。4.根据权利要求1至3任一所述的基于前向波气管分割的气管树提取方法，其特征在于，在所述步骤A之后，步骤B之前还包括: 使用阈值法去除所述胸部CT数据的背景，仅保留胸腔及胸腔内的器官数据。
【文档编号】G06T7/00GK105913422SQ201610214676
【公开日】2016年8月31日
【申请日】2016年4月7日
【发明人】唐宁, 龙杰, 陈凯
【申请人】广州医科大学附属肿瘤医院