删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

修正与融合的肺部分割算法

本站小编 Free考研考试/2020-03-23

窦圣昶, 赵海, 朱宏博, 王彬
东北大学 计算机科学与工程学院, 辽宁 沈阳 110169
收稿日期:2017-08-10
基金项目:辽宁省科技厅软课题(2015401039);辽宁省教育厅重点实验室基金资助项目(LZ2014015)。
作者简介:窦圣昶(1988-), 男, 山东泰安人, 东北大学博士研究生;
赵海(1959-), 男, 辽宁沈阳人, 东北大学教授, 博士生导师。

摘要:为了能在适当的计算复杂度下获得较好的分割效果, 在总结现有方法优缺点的基础上提出了将简单方法进行修正与融合的肺部分割方法.首先使用聚类的方法将CT图像中的像素点分为亮与暗两类, 获得灰度值不同的小区域, 然后区分出初始肺部的大致区域, 最后根据CT图像的具体特性, 使用识别出的初始肺部区域边界附近像素点的梯度以及灰度信息对边界进行修正.使用上述综合方法对CT图像的处理结果表明, 在计算复杂度不是太大的情况下能获得较好的分割效果.
关键词:CT图像肺结节融合聚类梯度
Modify and Integrate Lung Segmentation Algorithm
DOU Sheng-chang, ZHAO Hai, ZHU Hong-bo, WANG Bin
School of Computer Science & Engineering, Northeastern University, Shenyang 110169, China
Corresponding author: DOU Sheng-chang, E-mail: doushengc@qq.com
Abstract: In order to obtain better segmentation results under acceptable computational complexity, the advantages and disadvantages of the existing methods are summarized, and then the lung segmentation method which modifies and integrates the simple methods is proposed. Firstly, the clustering method is used to classify the pixels in CT images into two categories, and the small regions with different gray values are obtained. And then the initial region of the lung is identified. Finally, according to the specific characteristics of the CT image, the boundary is corrected by using the gradient and the gray information of the pixels near the initial lung boundary. The processing results of CT images using the above synthesis method show that the segmentation results can be better when the computational complexity is not too large.
Key words: CT imagespulmonary nodulesintegrateclusteringgradient
肺癌是现在致死率最高的癌症, 由于肺癌的治疗效果与发现的时期密切相关, 因此, 早发现早治疗成为降低肺癌致死率的重要方法[1].高风险人群的肺癌筛查是实现早发现的重要手段, 基于CT影像的筛查方式以其便捷、有效的特性成为肺癌筛查的主要方式.
随着肺癌筛查人群的增多, 生成的CT图像正处于急速膨胀之中, 传统的筛查方法是靠医生人工读取CT图像识别肺癌及其前期病变, 识别的效果受医生的经验以及疲劳程度等状态影响较大, 大量CT图像的辨识工作加重了医生的工作量, 使得医生无法在有效的时间内完成所有CT图像的诊断, 并保证没有漏诊的病例.帮助医生进行CT图像判读的计算机辅助系统成为很多医学以及计算机研究者们研究的重要领域[2-3].
计算机辅助诊断系统是对肺部CT图像进行辨识, 识别出候选的肺结节, 供医生参考, 系统可以分为3个功能部分, 即肺部区域分割、肺部结节识别、肺结节的性质判断.肺部区域分割功能是后续处理步骤的基础, 其分割的效果直接影响到后续处理部分的可靠性[4-5].现有的CT图像上的肺部分割方法多是直接基于灰度阈值的分割方法[6-7], 其执行效果并不理想, 在左右肺叶接近的CT图像上难以将左右肺区域分割开.近期有学者选择使用神经网络的方法对肺部区域进行分割[8-9], 然而, 其计算复杂度急速增长, 需要过于强大的计算硬件以及很长的计算时间, 同时用于训练的数据量非常大, 难以在实际中使用.
本文提出了基于多方法融合的肺部区域识别算法, 以提高肺部区域识别的有效性与准确性.
1 肺分割建模含有较多空气的肺部区域对X射线的吸收较少, 而包围肺部区域的胸墙对X射线有更强的吸收效果, 人体其他组织与肺部区域在CT图像上表现为不同的灰度值, 肺部以及体外的空气区域表现为亮度较暗的区域, 如图 1所示, 同时, 身体之外的物体以及身体内的不同部分对灰度上的区分有较大的干扰, 因此基于灰度值的肺部分割方法可以作为初始的预处理方法使用[10].本文首先使用了灰度聚类的方法对整个图像进行处理, 以获得初始的处理结果, 然后根据人体组织对X射线的吸收特性从处理结果中分割出肺部区域.
图 1(Fig. 1)
图 1 原始CT图像Fig.1 Original CT image

1.1 CT图像的聚类处理CT图像上的聚类处理是将图像中的像素点按照灰度值进行聚类以获得初始像素点分类信息的处理方法及其过程.Senthilkumar等[7]使用最小最大值的方法来确定初始的聚类中心, 该方法是一个递归的过程, 在此问题中初始聚类中心的选择对结果没有影响, 对求解过程的影响同样很小, 因此本文在使用该方法的过程中在保证结果相同的情况下舍弃了不必要的步骤.本文的算法过程:
1) 设置初始聚类中心.设置两个类别的聚类中心c1与c2,使得c1 < c2.
2) 计算图像I的每个点到两个聚类中心的距离矩阵d1和d2:
(1)
(2)
其中:I是图像矩阵;C1与C2分别是由聚类中心c1和c2扩展成的与I等大小的矩阵.
3) 由对应像素点处距离d1与d2的大小关系获得两个聚类结果矩阵K1与K2:
(3)
(4)
4) 由矩阵K1与K2的值更新聚类中心c1与c2.
(5)
(6)
其中, L1与L2分别是K1与K2中不等于零的像素点的个数.
5) 当c1与c1update以及c2与c2update间的差别足够小时结束, 否则重复执行2)~4)步.
聚类算法流程图如图 2所示.
图 2(Fig. 2)
图 2 聚类流程图Fig.2 Clustering flow chart

使用上述算法对图 1进行聚类,结果如图 3所示.在图 3中, 将图 1中的区域聚类成两个区域, 图 3a是聚类结果中的K1矩阵的可视化显示图像, 原来灰度值较小的区域保持原值, 灰度值较大的区域像素点的灰度值被置为零, 图 3b是聚类结果中的K2矩阵的可视化显示图像, 保留原值的是灰度值较大的区域, 其他区域被置为零, 两图之和即是原来的CT图像.
图 3(Fig. 3)
图 3 聚类结果Fig.3 Clustering result (a)—K1; (b)—K2.

1.2 肺部区域的初始分割对CT图像进行聚类处理后获得的是离散点的形式, 由于图像灰度的变化特性, 除个别明显区域外, 还不能形成统一的连通识别结果区域, 同时, 由于肺部高亮度肺结节以及血管的影响, 难以直接由聚类结果识别出肺部区域.本文使用基于胸墙的连通高亮度大区域特性实现胸墙的识别, 首先识别出胸墙区域, 进而实现肺部区域的初始识别.
在CT图像中, 包围肺部区域的胸墙区域由于对肺部的支撑作用, 是由质地较密的组织组成, 在图像中表现为较大区域的连续高亮度区域, 因此可以使用大面积的联通区域特性对其进行识别, 识别出的胸墙区域如图 4所示.
图 4(Fig. 4)
图 4 胸墙区域图Fig.4 Chest wall area figure

图 4中已经没有了影响肺部区域分割的小的区域, 并且将肺部CT图像分割成了明显的3部分区域, 对图 4进行处理能很容易获得初始的肺部区域.
2 肺区域边界修正模型由图 4可以明显地看出, 使用聚类后的结果分割出的肺部区域在左右肺比较近的部分无法区分左右肺叶, 会将其识别为一个区域, 这会对后续的处理造成影响, 因此需要对其进行处理, 本文提出了基于梯度向量场与边界灰度相结合的肺边界修正算法.
在识别出包含肺部的大致区域后, 在识别结果中存在两个问题, 即左右肺叶的分割问题以及肺部区域模糊边缘处会形成错误边界的问题.然而单一的基于梯度的方法无法解决该问题, 本文提出将梯度与灰度相结合的肺边界修正方法, 将像素点的梯度和灰度值结合在一起作为新的特征值用于边界识别.像素点的特征值为梯度与灰度的加权和:
(7)
式中:V为像素的灰度值; G为像素点处的梯度值; 实验中对ab取值为常数1.在初始边界上寻找初始的种子节点, 作为边界扩展的初始点, 在该像素点的梯度向量一侧寻找T值最大的像素点, 以获得下一个边界点.递归执行上述过程, 直到识别出的边界点形成环, 即获得了包含目标区域的边界点序列.算法的伪代码如下,结果如图 5所示.
表 (Table )
算法描述
Begin(算法开始)
输入初始分割肺部区域图像Image
计算图像的梯度Gradient
选取初始边界Edges上的种子节点P
While新边界NewEdges未形成封闭区域
????计算种子节点一侧180°内的像素点的特征值T
????选取T值最大的节点作为下一个种子节点
End While
End(算法结束)



图 5(Fig. 5)
图 5 肺部边界图Fig.5 Lung boundary figure

分析实际的CT图像可知, 在分割出的肺部区域的左右两边, 边界坐标值点具有清晰的边界以及区域特性, 同时初始寻找简单, 适于作为种子点, 实验中使用了该点作为种子点.
该方法的优势是, 识别出的边界一定是连续的封闭边界, 不需要后续的边界连续化操作.同时由于该方法的基础是前期获得初始预处理边界, 因此能获得较好的边界轮廓, 识别的肺部区域更精确.通过识别出的边界可以分割出准确的肺部区域范围.
3 实验过程与结果3.1 实验过程文中使用从盛京医院获得的肺部CT图像作为处理对象.CT图像经过如前所述的聚类、肺部识别以及边界修正后获得识别出的肺部区域, 整体的框架如图 6所示.
图 6(Fig. 6)
图 6 算法整体框图Fig.6 Block diagram of the algorithm

3.2 实验结果实验中对多种CT图像使用本文的融合方法进行肺部区域的分割识别, 对几种代表性的肺部CT图像的识别结果如图 7所示.
图 7(Fig. 7)
图 7 多实例肺部分割结果Fig.7 Multi-instance lung segmentation results (a)—原始CT图像;(b)—聚类法获得的初始肺部区域;(c)—修正后的边界;(d)—最后识别结果.

由图中可以看出, 本文将多种方法进行融合的肺部分割方法能很好地分割出CT图像中的肺部区域, 为后续的精确处理提供基础.
3.3 性能分析在聚类以及肺部区域初始识别的部分, 是对图像的像素点进行遍历操作, 其计算复杂度为O(n), 最后的边界修正部分同样是在像素上的操作, 遍历初始肺部区域的边界点, 因此其做大遍历的像素点数量不会超过图像中的像素点总数, 计算复杂度同样是O(n), 因此算法的整体时间复杂度为O(n).
4 结语在肺部分割的问题中, 单一的分割算法无法快速、精确地将肺部区域分割出来, 严重影响后续处理的精确性, 针对这一问题, 本文提出了基于多方法融合的肺部分割方法, 将灰度、梯度等信息进行综合作为识别特征.实验结果表明, 本文提出的方法能获得很好的肺部分割结果.
参考文献
[1]郭薇.面向多种医学图像的肺癌计算机辅助诊断关键技术研究[D].沈阳: 东北大学, 2011.
( Guo Wei.Study on the key techniques of computer-aided diagnosis for lung cancer in medical images[D].Shengyang: Northeastern University, 2011.http://cdmd.cnki.com.cn/Article/CDMD-10145-1015562144.htm)
[2] Jacobs C.Automatic detection and characterization of pulmonary nodules in thoracic CT scans[D].Nijmegen: Radboud University, 2015.
[3]Shen S, Bui A A T, Cong J, et al. An automated lung segmentation approach using bidirectional chain codes to improve nodule detection accuracy[J].Computers in Biology and Medicine, 2015, 57: 139–149.DOI:10.1016/j.compbiomed.2014.12.008
[4]Armato S G, Sensakovic W F. Automated lung segmentation for thoracic CT:impact on computer-aided diagnosis1[J].Academic Radiology, 2004, 11(9): 1011–1021.DOI:10.1016/j.acra.2004.06.005
[5]Dai S, Lu K, Dong J, et al. A novel approach of lung segmentation on chest CT images using graph cuts[J].Neurocomputing, 2015, 168: 799–807.DOI:10.1016/j.neucom.2015.05.044
[6]Do?anay E, Kara S, ?z?elik H K, et al. A hybrid lung segmentation algorithm based on histogram-based fuzzy C-means clustering[J].Computer Methods in Biomechanics and Biomedical Engineering:Imaging & Visualization, 2017(sup 1): 1–11.
[7]Senthilkumar T K, Ganesh E N, Umamaheswari R. Lung nodule volume growth analysis and visualization through auto-cluster k-means segmentation and centroid/shape variance based false nodule elimination[J].Biomedical Research, 2017, 28(5): 1–6.
[8] Harrison A P, Xu Z, George K, et al.Progressive and multi-path holistically nested neural networks for pathological lung segmentation from CT images[C]// International Conference on Medical Image Computing and Computer-Assisted Intervention.Quebec City: Springer, 2017: 621-629.
[9] Arbabshirani M R, Dallal A H, Agarwal C, et al.Accurate segmentation of lung fields on chest radiographs using deep convolutional networks[C]//SPIE Medical Imaging.Qrlando, 2017: 1013305-1013305-6.
[10]Van Rikxoort E M, de Hoop B, Viergever M A, et al. Automatic lung segmentation from thoracic computed tomography scans using a hybrid approach with error detection[J].Medical Physics, 2009, 36(7): 2934–2947.DOI:10.1118/1.3147146

相关话题/肺部 算法

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于改进粒子群算法的取件机械手轨迹综合优化设计
    黄裘俊1,张凯2,宋锦春1,于忠亮11.东北大学机械工程与自动化学院,辽宁沈阳110819;2.沈阳化工大学装备可靠性研究所,辽宁沈阳110142收稿日期:2017-08-16基金项目:国家自然科学基金青年基金资助项目(51605085)。作者简介:黄裘俊(1989-),男,福建三明人,东北大学博士 ...
    本站小编 Free考研考试 2020-03-23
  • 遗传算法优化神经网络的热轧带钢弯辊力预报模型
    王振华,龚殿尧,李广焘,张殿华东北大学轧制技术及连轧自动化国家重点实验室,辽宁沈阳110819收稿日期:2017-08-21基金项目:国家重点研发计划项目(2017YFB0304100);国家自然科学基金资助项目(51634002,51774084,51704067)。作者简介:王振华(1990-) ...
    本站小编 Free考研考试 2020-03-23
  • 一种改进的声发射P波到时自动提取算法
    刘飞跃,杨天鸿,邓文学,张鹏海东北大学资源与土木工程学院,辽宁沈阳110819收稿日期:2017-08-23基金项目:国家自然科学基金资助项目(51604062);国家重点基础研究发展计划项目(2016YFC0801602)。作者简介:刘飞跃(1996-),男,安徽阜阳人,东北大学博士研究生;杨天鸿 ...
    本站小编 Free考研考试 2020-03-23
  • 鱼群算法优化组合核函数GPR的油井动液面预测
    李翔宇1,高宪文1,李琨2,侯延彬11.东北大学信息科学与工程学院,辽宁沈阳110819;2.渤海大学工学院,辽宁锦州121013收稿日期:2015-12-22基金项目:国家自然科学基金资助项目(61573088,61433004,61403040)。作者简介:李翔宇(1982-),男,辽宁沈阳人, ...
    本站小编 Free考研考试 2020-03-23
  • 检测区域自适应调整的TLD多目标跟踪算法
    孟煜,张斌东北大学计算机科学与工程学院,辽宁沈阳110819收稿日期:2015-09-25基金项目:国家科技支撑计划项目(2014BAI17B00);国家关键科技研发基金资助项目(2015BAH09F02,2015BAH47F03);国家自然科学基金资助项目(61572116,61572117,61 ...
    本站小编 Free考研考试 2020-03-23
  • 基于在线判别分布域特征选择的鲁棒跟踪算法
    郭强1,2,吴成东1,赵迎春21.东北大学信息科学与工程学院,辽宁沈阳110819;2.中国刑事警察学院图书馆,辽宁沈阳110854收稿日期:2015-06-05基金项目:国家自然科学基金资助项目(61273078);辽宁省教育厅科学研究一般项目(L2015558);中央高校基本科研业务费专项资金资 ...
    本站小编 Free考研考试 2020-03-23
  • 基于阈值的车联网信标传输功率控制算法
    莫元富1,2,于德新1,包赛男2,高书涛11.吉林大学交通学院,吉林长春130022;2.中国汽车工程研究院,重庆401122收稿日期:2015-10-19基金项目:国家高技术研究发展计划项目(2014BAG03B03)。作者简介:莫元富(1984-),男,辽宁大连人,吉林大学博士研究生;于德新(1 ...
    本站小编 Free考研考试 2020-03-23
  • 基于隐马尔可夫模型和遗传算法的地图匹配算法
    吴刚,邱煜晶,王国仁东北大学信息科学与工程学院,辽宁沈阳110819收稿日期:2015-12-07基金项目:国家自然科学基金资助项目(61332006,61370154);中央高校基本科研业务费专项资金资助项目(N140404009)。作者简介:吴刚(1978-),男,辽宁沈阳人,东北大学副教授;王 ...
    本站小编 Free考研考试 2020-03-23
  • 基于多准则决策方法的WSNs不等簇数据收集算法
    宋晓莹,温涛,孙伟,张启龙东北大学软件中心,辽宁沈阳110819收稿日期:2015-11-05基金项目:国家自然科学基金资助项目(61170169,61170168)。作者简介:宋晓莹(1984-),女,辽宁朝阳人,东北大学博士研究生;温涛(1962-),男,辽宁大连人,东北大学教授,博士生导师。摘 ...
    本站小编 Free考研考试 2020-03-23
  • 基于块分割的新型压缩感知算法
    张娜1,曹琨2,刘亚轩11.中国海洋大学信息科学与工程学院,山东青岛266100;2.河南牧业经济学院信息与电子工程系,河南郑州450044收稿日期:2015-10-27基金项目:山东省自然科学基金资助项目(ZR2014FQ027);中央高校基本科研业务费专项资金资助项目(201513015)。作者 ...
    本站小编 Free考研考试 2020-03-23