肺部CT图像分割方法与装置.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN111583277A(43)申请公布日2020.08.25(21)申请号202010376872.3(22)申请日2020.05.07(71)申请人三峡大学地址443002湖北省宜昌市大学路8号(72)发明人夏平师冬霞任强雷帮军唐庭龙邹耀斌(74)专利代理机构北京兴智翔达知识产权代理有限公司11768代理人蒋常雪(51)Int.Cl.G06T7/11(2017.01)权利要求书2页说明书7页附图5页(54)发明名称肺部CT图像分割方法与装置(57)摘要本发明实施例提供一种肺部CT图像分割方法及装置，所述方法包括：获取待分割的肺部CT图像；对待分割的肺部CT图像进行双树复小波变换，确定各尺度上的双树复小波的低频分量；利用聚类算法获取最粗尺度上的初始分割；利用高斯模型进行参数估计；基于估计的参数，采用置信度传播算法通过全局能量最小获取尺度n上优化的分割结果；进行尺度内迭代；进行尺度间迭代；输出最终的分割结果。本发明实施例提供的肺部CT图像分割方法及装置，基于复小波域MRF置信度传播对肺部CT图像进行分割，分割结果中不但可以清晰地分割出感兴趣的区域，而且很好的保留了区域的细节部分以及边界信息。CN111583277ACN111583277A权利要求书1/2页1.一种肺部电子计算机断层扫描CT图像分割方法，其特征在于，包括：S1：获取待分割的肺部CT图像；S2：对待分割的肺部CT图像进行双树复小波变换，确定各尺度上的双树复小波的低频分量；S3：利用聚类算法获取最粗尺度上的初始分割；S4：利用高斯模型进行参数估计；S5：基于估计的参数，采用置信度传播算法通过全局能量最小获取尺度n上优化的分割结果；S6：进行尺度内迭代，重复S4和S5直到满足迭代准则，迭代停止；比较每次迭代后，肺结节图像分割结果的全局能量，输出全局能量最小的一次分割结果；S7：进行尺度间迭代，将尺度n的分割结果直接映射到相邻尺度n-1上，作为相邻尺度n-1的初始分割，重复S6，直到获取最细尺度上的分割结果；S8：输出最终的分割结果。2.根据权利要求1所述的肺部CT图像分割方法，其特征在于，所述基于估计的参数，采用置信度传播算法通过全局能量最小获取尺度n上优化的分割结果，具体包括：利用节点与节点之间相互传递信息而更新当前整个马尔科夫随机场MRF的标记状态；每次迭代后，计算图像新分割结果的全局能量。3.根据权利要求1所述的肺部CT图像分割方法，其特征在于，所述双树复小波变换为Q-shift双树复小波变换。4.根据权利要求1所述的肺部CT图像分割方法，其特征在于，所述聚类算法为K-means聚类算法。5.一种肺部CT图像分割装置，其特征在于，包括：获取模块，用于获取待分割的肺部CT图像；变换模块，用于对待分割的肺部CT图像进行双树复小波变换，确定各尺度上的双树复小波的低频分量；初始分割模块，用于利用聚类算法获取最粗尺度上的初始分割；参数估计模块，用于利用高斯模型进行参数估计；优化模块，用于基于估计的参数，采用置信度传播算法通过全局能量最小获取尺度n上优化的分割结果；尺度内迭代模块，用于进行尺度内迭代，重复S4和S5直到满足迭代准则，迭代停止；比较每次迭代后，肺结节图像分割结果的全局能量，输出全局能量最小的一次分割结果；尺度间迭代模块，用于进行尺度间迭代，将尺度n的分割结果直接映射到相邻尺度n-1上，作为相邻尺度n-1的初始分割，重复S6，直到获取最细尺度上的分割结果；输出模块，用于输出最终的分割结果。6.根据权利要求5所述的肺部CT图像分割装置，其特征在于，所述优化模块包括更新子模块和计算子模块；所述更新子模块用于利用节点与节点之间相互传递信息而更新当前整个MRF的标记状态；所述计算子模块用于每次迭代后，计算图像新分割结果的全局能量。2CN111583277A权利要求书2/2页7.根据权利要求5所述的肺部CT图像分割装置，其特征在于，所述双树复小波变换为Q-shift双树复小波变换。8.根据权利要求5所述的肺部CT图像分割装置，其特征在于，所述聚类算法为K-means聚类算法。9.一种电子设备，包括存储器、处理器，以及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时，实现如权利要求1至4任一项所述肺部CT图像分割方法的步骤。10.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，当所述计算机程序被处理器执行时，实现如权利要求1至4任一项所述肺部CT图像分割方法的步骤。3CN111583277A说明书1/7页肺部CT图像分割方法与装置技术领域[0001]本发明涉及图像处理技术领域，尤其涉及一种肺部CT图像分割方法与装置。背景技术[0002]马