北京大学学报(医学版) ›› 2021, Vol. 53 ›› Issue (1): 220-226. doi: 10.19723/j.issn.1671-167X.2021.01.034

• 技术方法 • 上一篇    下一篇

基于赋权形态学分析的三维面部对称参考平面构建方法

朱玉佳1,2,赵一姣1,2,郑盛文3,4,温奥楠1,2,傅湘玲3,4,Δ(),王勇1,2,Δ()   

  1. 1.北京大学口腔医学院·口腔医院,口腔医学数字化研究中心 国家口腔疾病临床医学研究中心 口腔数字化医疗技术和材料国家工程实验室 口腔数字医学北京市重点实验室,北京 100081
    2.北京大学口腔医学院·口腔医院口腔修复教研室,北京 100081
    3.北京邮电大学计算机学院(国家示范性软件学院),北京 100876
    4.北京邮电大学可信分布式计算与服务教育部重点实验室,北京 100876
  • 收稿日期:2020-10-10 出版日期:2021-02-18 发布日期:2021-02-07
  • 通讯作者: 傅湘玲,王勇 E-mail:fuxiangling@bupt.edu.cn;kqcadc@bjmu.edu.cn
  • 基金资助:
    国家自然科学基金(81870815);国家自然科学基金(82071171);北京大学口腔医院实验室开放课题(PKUSS20190501);宁夏回族自治区重点研发计划项目(2018BEG02012)

A method for constructing three-dimensional face symmetry reference plane based on weighted shape analysis algorithm

ZHU Yu-jia1,2,ZHAO Yi-jiao1,2,ZHENG Sheng-wen3,4,WEN Ao-nan1,2,FU Xiang-ling3,4,Δ(),WANG Yong1,2,Δ()   

  1. 1. Center of Digital Dentistry, Peking University School and Hospital of Stomatology & National Clinical Research Center for Oral Diseases & National Engineering Laboratory for Digital and Material Technology of Stomatology & Beijing Key Laboratory of Digital Stomatology, Beijing 100081, China
    2. Department of Prosthodontics, Peking University School and Hospital of Stomatology & National Clinical Research Center for Oral Diseases & National Engineering Laboratory for Digital and Material Technology of Stomatology & Beijing Key Laboratory of Digital Stomatology, Beijing 100081, China
    3. School of Computer Science, Beijing University of Posts and Telecommunications (National Pilot Software Engineering School), Beijing 100876, China
    4. Key Laboratory of Trustworthy Distributed Computing and Service, Ministry of Education, Beijing University of Posts and Telecommunications, Beijing 100876, China
  • Received:2020-10-10 Online:2021-02-18 Published:2021-02-07
  • Contact: Xiang-ling FU,Yong WANG E-mail:fuxiangling@bupt.edu.cn;kqcadc@bjmu.edu.cn
  • Supported by:
    National Natural Science Foundation of China(81870815);National Natural Science Foundation of China(82071171);Open Subject Foundation of Peking University School and Hospital of Stomatology(PKUSS20190501);Key R&D Program of Ningxia Hui Autonomous Region of China(2018BEG02012)

摘要:

目的: 建立一种基于三维形态学分析的赋权普氏分析(weighted Procrustes analysis,WPA)方法,实现对下颌偏斜患者面部解剖标志点不对称程度量化赋权的三维面部对称参考平面自动化构建。方法: 临床收集下颌偏斜患者30例,基于本体-镜像关联法的普氏分析(Procrustes analysis,PA)算法和自主研发的WPA算法,分别构建每名患者三维面部模型的对称参考平面(实验组一和二),以专家定义的区域迭代最近点(iterative closest point,ICP)算法(选取对称性良好的区域)构建的对称参考平面为“真值平面”,作为参考组。以实验组与参考组的平面角度误差作为评价指标,分别比较两个实验组与参考组的差异,评价WPA算法构建对称参考平面的临床适用性。结果: 30例下颌偏斜患者WPA算法构建的对称参考平面与真值平面的角度误差为1.53°±0.84°,PA算法构建的对称参考平面与真值平面的角度误差为2.06°±0.86°。对于下颌偏斜大于12 mm的患者,WPA算法比PA算法的平均平面角度误差小0.86°。结论: 基于三维面部形态学分析的WPA算法对下颌偏斜面部不对称畸形患者构建三维面部对称参考平面的临床应用效果比传统PA算法有显著提高,初步达到了口腔临床专家级诊断策略的效果。

关键词: 面部不对称, 对称参考平面, 普氏分析, 欧几里得距离矩阵分析

Abstract:

Objective: To establish a novel method based on three-dimensional (3D) shape analysis and weighted Procrustes analysis (WPA) algorithm to construct a 3D facial symmetry reference plane (SRP), automatically assigning weight to facial anatomical landmarks. The WPA algorithm suitability for commonly observed clinical cases of mandibular deviation were analysed and evaluated.Methods: Thirty patients with mandibular deviation were recruited for this study. The 3D facial SRPs were extracted independently based on original-mirror alignment method. Thirty-two anatomical landmarks were selected from the overall region by three times to obtain the mean coordinate. The SRP of experimental groups 1 and 2 were using the standard Procrustes analysis (PA) algorithm and WPA algorithm, respectively. A reference plane defined by experts based on regional iterative closest point (ICP) algorithm, served as the ground truth. Three experts manually selecting facial regions with good symmetry for original model, and common region was included in the study. The angle error values between the SRP of WPA algorithm in the experimental group 1 and the truth plane were evaluated in this study, and the SRP of PA algorithm of experimental group 2 was calculated in the same way. Statistics and measurement analysis were used to comprehensively evaluate the clinical suitability of the WPA algorithm to calculate the SRP. A paired t-test analysis (two-tailed) was conducted to compare the angles.Results: The average angle error between the SRP of WPA algorithm and the ground truth was 1.53°±0.84°, which was smaller than that between the SRP of PA and the ground truth (2.06°±0.86°). There were significant differences in the angle errors among the groups (P<0.05). For the patients with severe mandibular deviation that the distance between pogonion and facial midline greater than 12 mm, the average angle error of the WPA algorithm was 0.86° smaller than that of the PA algorithm.Conclusion: The WPA algorithm, based on weighted shape analysis, can provide a more adaptable SRP than the standard PA algorithm when applied to mandibular deviation patients and preliminarily simulate the diagnosis strategies of clinical experts.

Key words: Facial asymmetry, Symmetry reference plane, Procrustes analysis, Euclidean distance matrix analysis

中图分类号: 

  • R782.2

图1

研究采用的面部32个解剖标志点示意图"

图2

标志点间距离线段示意图"

图3

一例下颌偏斜患者由WPA算法、PA算法及专家构建的对称参考平面"

表1

PA算法和WPA算法构建对称参考平面的平面角度误差"

Subject number Angle error/(°)
PA WPA
1 2.93 2.56
2 0.66 0.27
3 2.51 2.2
4 2.06 0.85
5 3.05 1.86
6 1.55 0.93
7 2.09 1.38
8 1.17 0.94
9 1.93 1.44
10 2.76 1.88
11 1.67 1.05
12 0.61 0.31
13 2.1 1.94
14 4.22 1.66
15 3.31 3.09
16 2.43 2.39
17 0.63 0.43
18 0.74 0.23
19 1.61 1.29
20 1.3 0.86
21 2.27 0.83
22 2.42 2.16
23 2.29 2.17
24 1.7 1.09
25 3.16 3.01
26 2.35 2.25
27 2.08 2.1
28 1.13 0.29
29 2.39 2.57
30 2.66 1.89
x-±s 2.06±0.86 1.53±0.84

表2

不同程度下颌偏斜畸形数据PA算法和WPA算法构建对称参考平面的平面角度误差分布"

Mandibular deviation
degrees
Patient number Angle error/(°)
PA WPA
≥1 mm, <8 mm 12 1.65 1.16
≥8 mm, <12 mm 10 2.11 1.82
≥12 mm, <20 mm 8 2.60 1.74
[1] Baudouin JY, Tiberghien G. Symmetry, averageness, and feature size in the facial attractiveness of women[J]. Acta Psychol (Amst), 2004,117(3):313-332.
doi: 10.1016/j.actpsy.2004.07.002
[2] Ferrario VF, Sforza C, Ciusa V, et al. The effect of sex and age on facial asymmetry in healthy subjects: a cross-sectional study from adolescence to mid-adulthood[J]. J Oral Maxillofac Surg, 2001,59(4):382-388.
doi: 10.1053/joms.2001.21872 pmid: 11289167
[3] Shaner DJ, Peterson AE, Beattie OB, et al. Assessment of soft tissue facial asymmetry in medically normal and syndrome-affected individuals by analysis of landmarks and measurements[J]. Am J Med Genet, 2000,93(2):143-154.
doi: 10.1002/1096-8628(20000717)93:2<143::aid-ajmg12>3.0.co;2-q pmid: 10869118
[4] Burke PH, Healy MJ. A serial study of normal facial asymmetry in monozygotic twins[J]. Ann Hum Biol, 1993,20(6):527-534.
doi: 10.1080/03014469300002932 pmid: 8257078
[5] 黄佳梦, 李娟, 林军. 面部不对称畸形评估的研究进展[J]. 中国医疗美容, 2019,9(6):121-128.
[6] Kim MS, Lee EJ, Song IJ, et al. The location of midfacial landmarks according to the method of establishing the midsagittal reference plane in three-dimensional computed tomography analysis of facial asymmetry[J]. Imaging Sci Dent, 2015,45(4):227-232.
doi: 10.5624/isd.2015.45.4.227 pmid: 26730370
[7] Lee MS, Chung DH, Lee JW, et al. Assessing soft-tissue characteristics of facial asymmetry with photographs[J]. Am J Orthod Dentofacial Orthop, 2010,138(1):23-31.
doi: 10.1016/j.ajodo.2008.08.029 pmid: 20620830
[8] Haraguchi S, Iguchi Y, Takada K. Asymmetry of the face in orthodontic patients[J]. Angle Orthod, 2008,78(3):421-426.
doi: 10.2319/022107-85.1 pmid: 18416611
[9] O’Grady KF, Antonyshyn OM. Facial asymmetry: three-dimensional analysis using laser surface scanning[J]. Plast Reconstr Surg, 1999,104(4):928-937.
doi: 10.1097/00006534-199909040-00006 pmid: 10654730
[10] Nur RB, Cakan DG, Arun T. Evaluation of facial hard and soft tissue asymmetry using cone-beam computed tomography[J]. Am J Orthod Dentofacial Orthop, 2016,149(2):225-237.
doi: 10.1016/j.ajodo.2015.07.038 pmid: 26827979
[11] Lee JK, Jung PK, Moon CH. Three-dimensional cone beam computed tomographic image reorientation using soft tissues as reference for facial asymmetry diagnosis[J]. Angle Orthod, 2014,84(1):38-47.
doi: 10.2319/112112-890.1 pmid: 23758600
[12] Klingenberg CP, Barluenga M, Meyer A. Shape analysis of symmetric structures: quantifying variation among individuals and asymmetry[J]. Evolution, 2002,56(10):1909-1920.
doi: 10.1111/j.0014-3820.2002.tb00117.x pmid: 12449478
[13] Du S, Xu Y, Wan T, et al. Robust iterative closest point algorithm based on global reference point for rotation invariant registration[J]. PLoS One, 2017,12(11):e0188039.
doi: 10.1371/journal.pone.0188039 pmid: 29176780
[14] Damstra J, Fourie Z, De Wit M, et al. A three-dimensional comparison of a morphometric and conventional cephalometric mid-sagittal planes for craniofacial asymmetry[J]. Clin Oral Investig, 2012,16(1):285-294.
doi: 10.1007/s00784-011-0512-4 pmid: 21271348
[15] Besl PJ, McKay ND. A method for registration of 3-D shapes[J]. IEEE T Pattern Anal, 1992,14(2):239-256.
doi: 10.1109/34.121791
[16] Verhoeven TJ, Coppen C, Barkhuysen R, et al. Three dimensio-nal evaluation of facial asymmetry after mandibular reconstruction: validation of a new method using stereophotogrammetry[J]. Int J Oral Maxillofac Surg, 2013,42(1):19-25.
doi: 10.1016/j.ijom.2012.05.036 pmid: 22939875
[17] De Momi E, Chapuis J, Pappas I, et al. Automatic extraction of the mid-facial plane for cranio-maxillofacial surgery planning[J]. Int J Oral Maxillofac Surg, 2006,35(7):636-642.
doi: 10.1016/j.ijom.2006.01.028 pmid: 16542822
[18] Xiong Y, Zhao Y, Yang H, et al. Comparison between interactive closest point and procrustes analysis for determining the median sagittal plane of three-dimensional facial data[J]. J Craniofac Surg, 2016,27(2):441-444.
doi: 10.1097/SCS.0000000000002376 pmid: 26825747
[19] Richtsmeier JT, Lele S. A coordinate-free approach to the analysis of growth patterns: models and theoretical considerations[J]. Biol Rev Camb Philos Soc, 1993,68(3):381-411.
doi: 10.1111/j.1469-185x.1993.tb00737.x pmid: 8347767
[20] Lele S, Richtsmeier JT. Euclidean distance matrix analysis: a coordinate-free approach for comparing biological shapes using landmark data[J]. Am J Phys Anthropol, 1991,86(3):415-421.
doi: 10.1002/ajpa.1330860307 pmid: 1746646
[21] Verhoeven TJ, Nolte JW, Maal TJ, et al. Unilateral condylar hyperplasia: a 3-dimensional quantification of asymmetry[J]. PLoS One, 2013,8(3):e59391.
doi: 10.1371/journal.pone.0059391 pmid: 23544063
[22] Gateno J, Jajoo A, Nicol M, et al. The primal sagittal plane of the head: a new concept[J]. Int J Oral Maxillofac Surg, 2016,45(3):399-405.
doi: 10.1016/j.ijom.2015.11.013 pmid: 26708049
[23] Chetverikov D, Stepanov D, Krsek P. Robust Euclidean alignment of 3D point sets: the trimmed iterative closest point algorithm[J]. Image Vision Comput, 2005,23(3):299-309.
doi: 10.1016/j.imavis.2004.05.007
[24] Zelditch ML, Swiderski DL, Sheets DH, et al. Geometric morphometrics for biologists: A primer [M]. San Diego: Elsevier Academic Press, 2004: 293-319.
[25] Walters M, Claes P, Kakulas E, et al. Robust and regional 3D facial asymmetry assessment in hemimandibular hyperplasia and hemimandibular elongation anomalies[J]. Int J Oral Maxillofac Surg, 2013,42(1):36-42.
doi: 10.1016/j.ijom.2012.05.021 pmid: 22749574
[26] Claes P, Walters M, Clement J. Improved facial outcome assessment using a 3D anthropometric mask[J]. Int J Oral Maxillofac Surg, 2012,41(3):324-330.
doi: 10.1016/j.ijom.2011.10.019 pmid: 22103995
[27] Claes P, Walters M, Vandermeulen D, et al. Spatially-dense 3D facial asymmetry assessment in both typical and disordered growth[J]. J Anat, 2011,219(4):444-455.
doi: 10.1111/j.1469-7580.2011.01411.x
[28] Wu J, Heike C, Birgfeld C, et al. Measuring symmetry in children with unrepaired cleft lip: defining a standard for the three-dimensional midfacial reference plane[J]. Cleft Palate Craniofac J, 2016,53(6):695-704.
doi: 10.1597/15-053 pmid: 26752127
[29] Perlyn CA, DeLeon VB, Babbs C, et al. The craniofacial phenotype of the Crouzon mouse: analysis of a model for syndromic craniosynostosis using three-dimensional MicroCT[J]. Cleft Palate Craniofac J, 2006,43(6):740-748.
doi: 10.1597/05-212 pmid: 17105336
[30] Richtsmeier JT, Cole TR, Krovitz G, et al. Preoperative mor-phology and development in sagittal synostosis[J]. J Craniofac Genet Dev Biol, 1998,18(2):64-78.
pmid: 9672839
[31] Richtsmeier JT, Lele S. Analysis of craniofacial growth in Crouzon syndrome using landmark data[J]. J Craniofac Genet Dev Biol, 1990,10(1):39-62.
pmid: 2373755
[32] Dufresne C, Richtsmeier JT. Interaction of craniofacial dysmorphology, growth, and prediction of surgical outcome[J]. J Craniofac Surg, 1995,6(4):270-281.
doi: 10.1097/00001665-199507000-00003 pmid: 9020701
[33] Ferrario VF, Sforza C, Pizzini G, et al. Sexual dimorphism in the human face assessed by Euclidean distance matrix analysis[J]. J Anat, 1993,183(Pt 3):593-600.
[34] Ferrario VF, Sforza C, Miani AJ, et al. Maxillary versus mandi-bular arch form differences in human permanent dentition assessed by Euclidean-distance matrix analysis[J]. Arch Oral Biol, 1994,39(2):135-139.
doi: 10.1016/0003-9969(94)90108-2 pmid: 8185498
[35] Ferrario VF, Sforza C, Miani AJ, et al. Human dental arch shape evaluated by Euclidean-distance matrix analysis[J]. Am J Phys Anthropol, 1993,90(4):445-453.
doi: 10.1002/ajpa.1330900405 pmid: 8476003
[36] Ferrario VF, Sforza C, Miani AJ, et al. Dental arch asymmetry in young healthy human subjects evaluated by Euclidean distance matrix analysis[J]. Arch Oral Biol, 1993,38(3):189-194.
doi: 10.1016/0003-9969(93)90027-j pmid: 8489412
[37] 聂琼, 林久祥. Angle Ⅱ1与正常牙合牙弓形态差异——应用欧几里德距离矩阵分析法[J]. 实用口腔医学杂志, 2005,21(5):659-662.
[1] 邱淑婷,朱玉佳,王时敏,王飞龙,叶红强,赵一姣,刘云松,王勇,周永胜. 姿势微笑位口唇对称参考平面的数字化构建及初步应用验证[J]. 北京大学学报(医学版), 2022, 54(1): 193-199.
[2] 朱玉佳,许晴,赵一姣,张磊,付子旺,温奥楠,高梓翔,张昀,傅湘玲,王勇. 深度学习算法辅助构建三维颜面正中矢状平面[J]. 北京大学学报(医学版), 2022, 54(1): 134-139.
[3] 熊玉雪, 杨慧芳, 赵一姣, 王勇. 两种评价面部三维表面数据不对称度方法的比较[J]. 北京大学学报(医学版), 2015, 47(2): 340-343.
[4] 张晓芸, 许天民. 应用锥束计算机体层摄影术评价颅颌面的不对称[J]. 北京大学学报(医学版), 2013, 45(1): 156-161.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 田增民, 陈涛, Nanbert ZHONG, 李志超, 尹丰, 刘爽. 神经干细胞移植治疗遗传性小脑萎缩的临床研究(英文稿)[J]. 北京大学学报(医学版), 2009, 41(4): 456 -458 .
[2] 郭岩, 谢铮. 用一代人时间弥合差距——健康社会决定因素理论及其国际经验[J]. 北京大学学报(医学版), 2009, 41(2): 125 -128 .
[3] 成刚, 钱振华, 胡军. 艾滋病项目自愿咨询检测的技术效率分析[J]. 北京大学学报(医学版), 2009, 41(2): 135 -140 .
[4] 卢恬, 朱晓辉, 柳世庆, 郑杰, 邱晓彦. 白细胞介素2促进宫颈癌细胞系HeLaS3免疫球蛋白G的表达[J]. 北京大学学报(医学版), 2009, 41(2): 158 -161 .
[5] 袁惠燕, 张苑, 范田园. 离子交换型栓塞微球及其载平阳霉素的制备与性质研究[J]. 北京大学学报(医学版), 2009, 41(2): 217 -220 .
[6] 徐莉, 孟焕新, 张立, 陈智滨, 冯向辉, 释栋. 侵袭性牙周炎患者血清中抗牙龈卟啉单胞菌的IgG抗体水平的研究[J]. 北京大学学报(医学版), 2009, 41(1): 52 -55 .
[7] 董稳, 刘瑞昌, 刘克英, 关明, 杨旭东. 氯诺昔康和舒芬太尼用于颌面外科术后自控静脉镇痛的比较[J]. 北京大学学报(医学版), 2009, 41(1): 109 -111 .
[8] 祁琨, 邓芙蓉, 郭新彪. 纳米二氧化钛颗粒对人肺成纤维细胞缝隙连接通讯的影响[J]. 北京大学学报(医学版), 2009, 41(3): 297 -301 .
[9] Jian-wei GU, Emily YOUNG, Zhi-jun PAN, Kevan B. TUCKER, Megan SHPARAGO, Min HUANG, Amelia Purser BAILEY. SD大鼠长期高盐饮食可导致其高血压并改变肾细胞因子基因表达谱[J]. 北京大学学报(医学版), 2009, 41(5): 505 -515 .
[10] 李宏亮*, 安卫红*, 赵扬玉, 朱曦. 妊娠合并高脂血症性胰腺炎行血液净化治疗1例[J]. 北京大学学报(医学版), 2009, 41(5): 599 -601 .