基于水平集的牙齿CT图像分割技术作者：汪葛王远军来源：《计算机应用》2016年第03期摘要：牙齿的计算机断层扫描（CT）图像中存在边界模糊、相邻牙齿粘连等情况，且拓扑结构较为复杂，要实现准确的牙齿分割非常困难。
对传统的牙齿CT图像分割方法，特别是近年来用于牙齿分割的水平集方法进行介绍，对其水平集函数中各能量项进行研究，并通过对比实验体现水平集方法的优越性。
基于水平集的牙齿CT图像分割方法中水平集函数的能量项主要包括：竞争能量项、梯度能量项、形状约束能量项、全局先验灰度能量项、局部灰度能量项。
实验结果表明基于混合模型的水平集方法分割效果最佳，切牙与磨牙分割准确率分别为88.92%和92.34%，相比自适应阈值和传统水平集方法，分割准确率总体提升10%以上。
在综合利用图像信息和先验知识的基础上，通过对水平集函数中能量项进行优化和创新，有望进一步提高分割的准确率。
关键词：牙齿锥形束计算机断层扫描（CBCT）图像；图像分割；水平集；能量项；混合模型中图分类号： TP391.413 文献标志码：A0引言近年来，人体组织器官的可视化技术已经成为计算机辅助诊断的重要工具。
由于牙齿的计算机断层扫描（Computer Tomography， CT）图像同时包含了牙冠和牙根的解剖信息，为重建完整的牙齿模型提供了可靠数据。
而牙齿的分割作为牙齿模型重建工作中的一个重要步骤，对牙齿分割方法的研究具有重要意义。
牙齿的形状和牙根的具体位置等信息对牙齿的正畸手术、种植手术、根管治疗等临床操作非常重要。
通常手术前需要对这些信息进行手动测量和获取，这个过程往往非常耗时，而且其准确性也不能达到非常高的要求。
因此，通过获得牙齿的三维数字模型可以为口腔疾病的诊断和手术治疗方案的制定等提供完整的解剖信息，极大地提高诊断的准确性和手术的成功率。
为了获得精确的三维信息，则必须要求对牙齿CT图像进行准确分割。
目前基于水平集的图像分割算法已广泛应用于医学图像分割领域，这一类算法可以很好地解决外形复杂和拓扑结构变化剧烈的情况。
本文主要针对水平集方法（Level Set Method， LSM）在牙齿CT图像分割应用中的研究进展进行综述和讨论。
3结语本文针对基于水平集的牙齿CT图像分割技术的研究进展进行探讨。
首先对传统的分割方法进行概述，然后在简要介绍水平集理论之后对应用于牙齿CT图像分割的水平集方法进行详细介绍，并对其水平集函数的能量项进行深入研究，最后，通过分割实验结果详细比较了当前典型的三种方法的精确性。
通过对水平集函数中能量项的改进与创新，基于混合模型的水平集方法在分割CT牙齿图像上已取得了较为满意的结果。
但对于初始层的分割，需要操作人员手动选择初始层的位置，如何减少主观影响、提高算法的自动化程度是未来改进的方向；且目前此类方法的研究针对的都是无金属伪影和阻生牙干扰的情况，为满足临床实际的需要，未来还需深入研究各能量项的混合模型，以实现对牙冠、牙根的精确分割。
参考文献：[1]MEMON A， GODWARD S， WILLIAMS D， et al. Dental xrays and the risk of thyroid cancer： a casecontrol study [J]. Acta Oncologica， 2010， 49（4）： 447-453.[2]LONGSTRETH W T， Jr， PHILLIPS L E， RANGSHOLT M， et al. Dental Xrays and the risk of intracranial meningioma： a populationbased casecontrol study [J]. Cancer， 2004， 100（5）： 1026-1034.[3]GAO H， CHAE O. Individual tooth segmentation from CT images using level set method with shape and intensity prior [J]. Pattern Recognition， 2010， 43（7）： 2407-2417.[4]RAHIMI A， BENDELS L K， KLEIN R， et al. 3D reconstruction of dental specimens from 2D histological images and μCTScans [J]. Computer Metho ds in Biomechanics and Biomedical Engineering， 2005， 8（3）： 167-176.[5]HEO H， CHAE OS. Segmentation of tooth in CT images for the 3D reconstruction of teeth [C]// Proceedings of the Imaging Processing： Algorithms and Systems Ⅲ， SPIE 5298. Bellingham： SPIE， 2004： 455-466.[6]WU X， GAO H， HEO H， et al. Improved BSpline contour ftting using genetic algorithm for the segmentation of dental computerized tomography image sequences [J]. Journal of Imaging Science and Technology， 2007，51（4）：328-336.[7]CHEN T， METAXAS D. A hybrid framework for 3D medical image segmentation [J]. Medical Image Analysis， 2005， 9（6）： 547-565.[8]KEYHANINEJAD S， ZOROOFI R A， SETAREHDAN S K， et al. Automated segmentation of teeth in multislice CT images [C]// VIE 2006： Proceedings of the 2006 IET International Conference on Visual Information Engineering. Piscataway， NJ： IEEE， 2006：339-344.[9]HOSNTALAB M， ZOROOFI R A， TEHRANIFARD A A， et al. Segmentation of teeth in CT volumetric dataset by panoramic projection and variational level set [J]. International Journal of Computer Assisted Radiology and Surgery， 2008， 3（3/4）： 257-265.[10]KEUSTERMANS J， VANDERMEULEN D， SUETENS P. Integrating statistical shape models into a graph cut framework for tooth segmentation [C]// MLMI 2012： Proceedings of the Third International Workshop on Machine Learning in Medical Imaging， LNCS 7588. Berlin：Springer， 2012： 242-249.[11]HIEW L T， ONG S H， FOONG K W C. Tooth segmentation from conebeam CT using graph cut [C]// APSIPA ASC 2010： Proceedings of the 2010 2nd AsiaPacific Signal and Information Processing Association Annual Summit and Conference. Biopolis， Singapore： AsiaPacific Signal and Information Processing Association， 2010： 272-275.[12]MARKEL D， ZAIDI H， EI NAQA I. Novel multimodality segmentation using level sets and JensenRényi divergence [J]. Medical Physics， 2013， 40（12）： 121908.【1934-1936】[13]ABDOLI M， DIERCKX R A， ZAIDI H. Contourletbased active contour model for PET image segmentation [J]. Medical Physics， 2013，40（12）：1934-1936.[14]QIU W， YUAN J， UKWATTA E， et al. Threedimensional prostate segmentation using level set with shape constraint based on rotational slices for 3D endfiring TRUS guided biopsy [J]. Medical Physics， 2013， 40（7）：1303-1305.[15]FRITSCHER K D， PERONI M， ZAFFINO P， et al. Automatic segmentation of head and neck CT images for radiotherapy treatment planning using multiple atlases， statistical appearance models， and geodesic active contours [J]. Medical Physics， 2014， 41（5）： 051910.[16]OSHER S， SETHIAN J A. Fronts propagating with curvaturedependent speed： algorithms based on HamiltonJacobi formulations [J]. Journal of Computational Physics， 1988， 79（1）：12-49.[17]YAU HT， YANG TJ， CHEN YC. Tooth model reconstruction based upon data fusion for orthodontic treatment simulation [J]. Computers in Biology and Medicine， 2014， 48： 8-16.[18]JI D X， ONGA S H， FOONG K W C. A levelset based approach for anterior teeth segmentation in cone beam computed tomography images [J]. Computers in Biology and Medicine，2014， 50： 116-128.[19]GAN Y， XIA Z， XIONG J， et al. Toward accurate tooth segmentation from computed tomography images using a hybrid level set model [J]. Medical Physics， 2015， 42（1）： 14-27.[20]GAO H， CHAE O. Touching tooth segmentation from CT image sequences using coupled level set method [C]// VIE 2008： Proceedings of the 2008 5th International Conference on Visual Information Engineering. Piscataway， NJ： IEEE， 2008： 382-387.[21]DUFOUR A， SHININ V， TAJBAKHSH S， et al. Segmenting and tracking fluorescent cells in dynamic 3D microscopy with coupled active surfaces [J]. IEEE Transactions on Image Processing， 2005， 14（9）： 1396-1410.[22]CHEN H， JAIN A K. Tooth contour extraction for matching dental radiographs [C]// ICPR 04： Proceedings of the 2004 17th International Conference on Pattern Recognition. Washington，DC： IEEE Computer Society， 2004， 3： 522-525.[23]PARAGIOS N， DERICHE R. Geodesic active regions： a new framework to deal with frame partition problems in computer vision [J]. Journal of Visual Communication and Image Representation， 2002， 13（1/2）： 249-268.[24]CHAN T， ZHU W. Level set based shape prior segmentation [C]// CVPR 2005： Proceedings of the 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition. Washington， DC： IEEE Computer Society， 2005： 1164-1170.[25]PLUEMPITIWIRIYAWEJ C， MOURA J M F， WU YJ L， et al. STACS： new active contour scheme for cardiac MR image segmentation [J]. IEEE Transactions on Medical Imaging，2005， 24（5）： 593-603.[26]LI C， KAO CY， GORE J C， et al. Implicit active contours driven by local binary fitting energy [C]// CVPR 2007： Proceedings of the 2007 IEEE Computer Society Conference on Computer Vision and Pattern Recognition. Washington， DC： IEEE Computer Society， 2007： 1-7.[27]WANG L， HE L， MISHRA A， et al. Active contours driven by local Gaussian distribution fitting energy [J]. Signal Processing， 2009， 89（12）： 2435-2447.[28]LI C， XU C， GUI C， et al. Level set evolution without reinitialization： a new variational formulation [C]// CVPR 2005： Proceedings of the 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition. Washington， DC： IEEE Computer Society， 2005：430-436.BackgroundThis work is partially supported by the National Natural Science Foundation of China （61201067） and Scientific Research Innovation Project of Shanghai Municipal Education Committee （13YZ069）.WANG Ge， born in 1992， M. S. candidate. His research interests include medical image segmentation.WANG Yuanjun， born in 1980， Ph. D.， associate professor. His research interests include biomedical engineering， medical image processing and analysis.龙源期刊网 。