针对婴幼儿颅缝早闭的颅骨重塑手术模拟方法

Simulation method of skull remodellingsurgeryfor infant with craniosynostosis

珅宇陆 , ¹ 杨宇罗 , ^1, ^* 文键郑 , ² and 剑宫 ²

珅宇陆

中国科学院自动化研究所（北京 100190）, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, P.R.China

Find articles by 珅宇陆

杨宇罗

中国科学院自动化研究所（北京 100190）, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, P.R.China

Find articles by 杨宇罗

文键郑

中国科学院自动化研究所（北京 100190）, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, P.R.China

Find articles by 文键郑

剑宫

中国科学院自动化研究所（北京 100190）, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, P.R.China 中国科学院自动化研究所（北京 100190）, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, P.R.China 首都医科大学附属北京天坛医院（北京 100070）, Beijing Tiantan Hospital, Capital Medical University, Beijing 100070, P.R.China

Corresponding author.

杨宇罗: [email protected]

罗杨宇，Email： [email protected]

式中 x 为月龄， y 为颅腔容积（单位为 mL）。有限元模拟儿童颅骨生长的方法为：采用标准生长曲线上不同年龄对应的颅内容积（intracranial volume，ICV）作为计算的目标值，对 ICV 进行热膨胀以模拟儿童脑部的生长发育。ICV 与外部颅骨接触后，挤压外部颅骨使颅骨发生形变以模拟在该年龄时颅骨的形状。

本文基于 CT 影像提取了患者的颅骨以及颅脑和硬膜构成的颅腔，并在 Mimics20（Materialise，Leuven，比利时）中对患者进行虚拟手术。术后的颅骨模型经修复和实体化后导入 Hypermesh 13.0（Altair，Troy，美国）中进行有限元模型的前处理，包括对模型划分网格、材料赋值、设置接触方式以及设置约束和加载方式。有限元模拟患儿颅骨生长过程采用的加载方式为颅腔的线性热膨胀，颅腔与颅骨之间的接触对采用罚函数算法。由 ABAQUS 2016（Dassault Systems）对模型进行计算并完成后处理。

2. 实验

本研究通过了首都医科大学北京天坛医院的伦理审查。本文以一例典型的矢状缝早闭儿童为对象，研究针对矢状缝早闭的颅骨重塑手术方案生成方法，并建立有限元模型以模拟术后的生长趋势。该病例为 2017 年 10 月在首都医科大学附属北京天坛医院就诊的一例男性患儿，首次就诊的年龄为 9 个月，图 4 显示的是由患者 CT 影像建立的颅骨模型。患儿矢状缝闭合，头部呈现舟状头的特征。患儿成长到 18 个月时拍摄头部 CT，并于 3 个月后（月龄 21 个月）接受了 π 型颅骨重塑手术，术后一年（月龄 33 个月）再次拍摄了头部 CT。

An external file that holds a picture, illustration, etc. Object name is swyxgcxzz-38-5-932-4.jpg

Open in a separate window

图 4

A three-dimensional model of a 9-month-old child with premature fusion of sagittal suture

9 个月患儿颅骨的三维模型

2.1. 畸形骨瓣区域与理想颅骨模型的配准

通过 Mimics 建立患儿颅骨的点云模型，选取患儿的乳突和顶结节区域作为特征与颅骨数据库内同龄、同性别正常儿童的颅骨进行匹配，得到对应该患儿的理想颅骨模型。图 5 所示为患儿颅骨与正常儿童颅骨数据库的配准结果，图 5 a 下方为该患儿对应正常颅骨模型。

An external file that holds a picture, illustration, etc. Object name is swyxgcxzz-38-5-932-5.jpg

Open in a separate window

图 5

The result of point cloud registration between the child’s skull and the ideal skull

患儿颅骨与理想颅骨配准结果

a.配准区域；b.配准结果

a. registration area; b. registration result

选取患儿的左侧额结节区域点云、颅顶区域点云以及顶结节区域点云与相应理想颅骨的左侧点云进行配准。图 6 为患儿颅骨形态异常区域与理想颅骨模型的匹配结果，患儿额结节点云与正常理想颅骨的顶结节部位形态吻合，患儿颅顶区节点云与正常理想颅骨的额结节部位形态吻合，患儿顶结节点云与正常理想颅骨的颅顶区部位形态吻合。

An external file that holds a picture, illustration, etc. Object name is swyxgcxzz-38-5-932-6.jpg

Open in a separate window

图 6

The matching result of the abnormal area of the child with the normal area of the ideal skull

患儿颅骨形态异常区域与理想颅骨模型的匹配结果

2.2. 骨瓣切割与错位拼接虚拟手术

由 Mimics 导入医学数字成像和通信（Digital Imaging and Communications in Medicine，DICOM）格式的 9 个月患儿的 CT 影像。由于有限元模拟生长的模型是通过颅腔膨胀与颅骨发生接触后使颅骨形变得到的，因此我们在 Mimics 中通过 CT 阈值分割的方法 ^{[

21

]} ，提取出患儿的颅骨以及颅腔。其中颅骨的阈值范围 226～3 071 HU，颅腔的 CT 阈值范围为 − 3～139 HU。在颅骨的区域内进行区域增长计算，剥离出完整的颅骨区域。因为颅腔内的 CT 值与皮肤、器官相近，且颅腔内部结构复杂，因此需手工将颅腔区域划分开，并计算出颅骨和颅腔的三维（three-dimensional，3D）模型。

根据点云配准获得的边界线对额结节、颅顶区域、顶结节进行切割，如图 7 所示。额结节骨瓣的割线上方距离冠状缝约为 21 mm，下方距离眼眶上缘约 9 mm。将两侧的额结节切割并从颅骨分离后得到额结节骨瓣。顶骨在距离前囟点 9.92 mm 处为切割定位点，以额结节骨瓣在颅顶区域的正投影尺寸进行切割，得到颅顶区域骨瓣。顶结节区域按标记尺寸进行切割，得到顶结节骨瓣。

An external file that holds a picture, illustration, etc. Object name is swyxgcxzz-38-5-932-7.jpg

Open in a separate window

图 7

Surgical cutting plan for premature cranial suture fusion

颅缝早闭手术切割方案

骨瓣的拼接方案如图 8 所示，将颅顶区域骨瓣作为新的前额，额结节骨瓣移动至顶结节位置处，顶结节骨瓣移动至颅顶。该拼接方案减小了患者额结节的凸出，并利用额结节的自然形态给颅腔提供了生长的空间。将虚拟手术后的点云文件导入逆向工程软件 Geomagic Wrap 2017（3D Systems，Rock Hill，美国）对模型进行实体化 ^{[

22

]} 。Geomagic Wrap 将切割后的颅骨、骨瓣以及颅腔构造为实体模型。将实体化的模型导入 CATIAV5（Dassault System）软件中用连接片固定骨瓣与颅骨。

An external file that holds a picture, illustration, etc. Object name is swyxgcxzz-38-5-932-8.jpg

Open in a separate window

图 8

Scheme diagram of bone plates repositioning

骨瓣拼接方案图

2.3. 有限元模型

为了模拟患儿颅骨术后 9 个月的生长趋势，我们以颅内容积作为患儿术后年龄的判断依据。7～24 个月男童的 Lichtenberg 标准生长曲线 ^{[

20

]} 如图 9 所示，在标准曲线上年龄为 9 个月和 18 个月儿童的颅内容积分别对应 954.5 mL 与 1 123 mL。我们选取的患儿 9 个月的颅腔容积为 918.65 mL。患儿在未经手术干预的情况下成长到 18 个月，通过 CT 测量患儿的颅腔容积为 1 152.43 mL。患儿初始颅腔与 18 个月大的颅腔容积与 Lichtenberg 提出的标准颅腔生长曲线上 9 个月以及 18 个月大的儿童的颅腔体积大致相符。因此，我们使用热膨胀的方法，将颅腔容积增加到约为 1 123 mL，相当于儿童生长到 18 个月的标准颅腔容积。

An external file that holds a picture, illustration, etc. Object name is swyxgcxzz-38-5-932-9.jpg

Open in a separate window

图 9

Standard ICV growth curve for 7-24 months children

7～24 月龄标准颅腔容积生长曲线

在 Hypermesh 软件中对颅骨模型划分网格，使用四面体网格对颅骨进行划分。连接片为二维（two-dimensional，2D）结构，采用四边形网格对连接片进行划分。骨瓣有限元模型有节点 17 101 个、一阶四面体网格 51 994 个；颅骨有限元模型有节点 56 837 个、四面体网格 193 225 个；颅腔有限元模型有节点 163 100 个、四面体网格 801 679 个。综合多个文献的研究，模型材料的赋值见表 1 ^{[

1

,

9

-

10

]} 。颅腔的材料属性较难获取 ^{[

23

]} ，由于与颅骨接触面主要是硬脑膜，我们参考硬脑膜对颅腔进行赋值。颅腔相对于颅骨材料较软，我们假设弹性模量为 500 MPa。硬脑膜的泊松比参考范围为 0.23～0.45 ^{[

24

]} ，我们假设颅腔的泊松比为 0.24。颅腔外表面主要颅骨和 6 个骨瓣的网格采用颅骨材料进行赋值。连接片与头骨的连接设置为点连接，在软件中创建 3DConnector，用连接片将骨瓣和颅骨连接。

表 1

Model material assignment

模型材料赋值

元件	弹性模量/MPa	泊松比
颅骨	2 500	0.3
连接片	5 000	0.25
颅腔	500	0.24

Open in a separate window

颅骨内部摩擦力较小，假设摩擦系数为 0.2。分别建立颅腔外表面与颅骨内表面、颅腔外表面与 6 个骨瓣内表面的接触对，接触对设置为面面接触。约束条件参考倪健 ^{[

25

]} 的研究，约束枕骨大孔 ^{[

22

]} 的 6 个自由度以及颅腔与颅底接触的平面。在颅腔中心处 31 个节点放置温度载荷，进行线性的热膨胀，颅腔内部容积增长至 1 113.73 mL，对应标准颅内容积生长曲线约为 18 个月。

3. 结果

3.1. 颅骨骨瓣切割与错位拼接

通过点云匹配的算法，我们在 120 例正常儿童的颅骨数据库中匹配到了对应于该患儿的理想模板。通过矢状缝早闭患儿的颅骨外形特征分析，患儿颅骨形态异常区域主要在额结节、颅顶区域以及顶结节。通过与理想颅骨的配准可以得出：患儿的额结节与模板的顶结节相似度最高，患儿的颅顶区域与模板的额结节相似度最高，患儿的顶结节与模板的颅顶区域相似度最高，在配准中心区域的匹配精度为 0.3 mm。由此我们得出该患儿的颅骨错位拼接方案：将额结节、颅顶区、顶结节骨瓣分别错位拼接到顶结节、额结节、颅顶区域。

3.2. 有限元模拟术后颅骨生长

图 10 为有限元计算的应变分布与原始模型的对比图，其中图 10 b 和 10 d 为模拟患儿颅骨生长到约 18 个月的结果。由图 10 a、 10 b 可以看出，颅顶区域的应变较大，给予颅腔沿高度生长的空间，形态上相较于术后更加符合理想颅骨的形状。颅骨经切割的区域形变较大，颅骨与骨瓣之间的过渡更加平滑，并使得颅腔可以沿横向进行生长，缓解了舟状头患者颅骨只沿纵向生长的情况。除了颅顶与顶结节区域外，其余区域形变较小且较为均匀。从图 10 c、 10 d 可以看出，正视图中生长后的颅骨更加光滑，横向尺寸相较于手术后有了明显的变化。颅骨两侧的形变量较为一致，生长后的颅骨对称性较好。

An external file that holds a picture, illustration, etc. Object name is swyxgcxzz-38-5-932-10.jpg

Open in a separate window

图 10

ABAQUS calculation result

ABAQUS 计算结果

a.虚拟术后颅骨的斜侧视图；b.模拟生长后颅骨的斜侧视图；c.虚拟术后颅骨的正视图；d.模拟生长后颅骨的正视图

a. side view of the skull after virtual surgery; b. side view of the skull after simulated growth; c. front view of the skull after virtual surgery; d. front view of the skull after simulated growth

CI 可以在一定程度上表征颅骨的形态，为头颅的最宽径与最长径的比值，正常范围在 75%～85% ^{[

26

]} ，矢状缝早闭患者通常 CI < 0.7 ^{[

27

]} 。图 11 为患儿月龄为 9 个月、18 个月、模拟生长至 18 个月以及真实术后一年的 CI 计算图。选择头骨前后径最大值所在的横断面测量患儿的 CI，9 个月时患儿颅骨原始的 CI 值为 65.31%。患儿在成长至 18 个月时，颅骨的 CI 值降为 64.4%，舟状头的特征更为明显。虚拟手术后患儿颅骨的 CI 值为 66.56%，有限元模拟生长到 18 个月的患儿颅骨的 CI 值增长为 71.50%。在颅腔容积大致相同的条件下，经颅骨重塑手术与未经手术干预成长到 18 个月的患儿 CI 值有着较为显著的差异。

An external file that holds a picture, illustration, etc. Object name is swyxgcxzz-38-5-932-11.jpg

Open in a separate window

图 11

Results of cranial index calculation

头颅指数计算结果

测量在计算过程中颅骨的 CI 值，可以得到 CI 随颅腔容积变化曲线，如图 12 所示。使用该切割方案进行手术后的有限元模拟，生长过程中随着颅腔内的容积增加，CI 值也在逐渐增长，患儿颅骨的形态得到了矫正。患儿在真实手术一年后颅腔容积增加至 1 312.11 mL，患儿 CI 增加至 70.97%，与我们提出的颅骨骨瓣切割、重拼接方案模拟术后 9 个月的 CI 相差 0.53%。从本文的颅骨拼接手术方案模拟的生长效果与患儿 9 个月颅骨对比看出，患儿的颅骨在纵向上的尺寸基本没有变化，患儿额结节的凸出得到了改善。

An external file that holds a picture, illustration, etc. Object name is swyxgcxzz-38-5-932-12.jpg

Open in a separate window

图 12

The changing curve of cranial index with cranial volume

头颅指数随颅内容积增长变化曲线

在不同横切面提取术前、虚拟手术后和生长模拟后在同一坐标系下的颅骨的二维轮廓，保持各层之间的比例不变得到图 13 。其中黑色曲线为 9 个月患儿术前颅骨轮廓，黄色曲线为虚拟手术后颅骨轮廓，蓝色曲线为有限元模拟生长的颅骨轮廓。从术后与术前的轮廓对比图中可以看出，手术减小了患儿额结节的凸出量，在横向弥补了头宽的不足。第一层截面可以观察到颅顶区域骨瓣以及顶结节骨瓣，生长模拟后两块区域骨瓣横向生长趋势明显，在形态上与未切割的颅骨融合度较好。第二层显示，顶结节骨瓣和切割的骨缝带动了颅骨侧面的生长，模拟生长后的颅骨两侧形变量相近，额结节的形变量较小，且相较于术前的额结节凸出量有所减少。第三层截面模拟生长后的颅骨二维轮廓在纵向上与术前尺寸相近的情况下，横向的尺寸明显大于术前。第四层截面接近于颅底，截面上的切割区域较少，三者的差异不明显。由此说明了该套截骨方案可以有效地修复患儿颅骨在形态上的缺陷，在扩充颅内容积的同时达到促进美观的效果。

An external file that holds a picture, illustration, etc. Object name is swyxgcxzz-38-5-932-13.jpg

Open in a separate window

图 13

Two-dimensional contours of the skull on different cross sections

在不同横切面上颅骨的二维轮廓

4. 结论

传统颅缝早闭患者的手术治疗方案常依赖于医生的经验，医生对术后的效果以及患者颅骨的生长趋势很难做出判断。针对这一问题，我们提出了以理想颅骨模型为参考的个性化手术方案生成方法，以及基于有限元模型的手术模拟方法。

我们通过患儿头部 CT 与理想颅骨模型匹配后提出一种针对舟状头的切割与拼接方案，即用患儿较为凸出的额结节替换凸出量不足的顶结节，顶结节替换较为扁平的颅顶以提供生长的空间，颅顶的骨瓣替换额结节，从而改善颅骨的形态，扩充颅脑生长空间。本文以颅腔的膨胀与颅骨的相互作用模拟颅骨生长，以标准颅腔容积曲线将有限元计算结果与患儿的年龄相对应，模拟了一例舟状头患儿在上述手术方案后 9 个月的颅骨生长趋势。结果表明，患儿在术后颅顶区域和顶结节区域的应变较大，对应于颅骨在高度和宽度上的尺寸变化，其余区域的应变较为平均，颅骨两侧对称性保持较好。

患儿 9 个月时的 CI 指标为 65.31%，虚拟术后与真实术后的 CI 指标分别为 71.50% 和 70.97%，证明了本文所提手术方案的有效性与模拟方法的可靠性，在满足临床术前规划要求、降低对医生水平和经验依赖的同时，提高了手术方案的质量与标准化程度。本文的后续工作将围绕颅缝早闭手术的验证性和标准化展开，将针对多位舟状头患儿病例开展方法的验证性研究，以评估模拟方法的可靠性以及有效性。手术过程的标准化将以点云配准的颅骨骨瓣切割方案为指导，通过颅骨重塑手术机器人辅助医生定位骨瓣的切割边界，从而缩短手术时间，提高手术的效率。

利益冲突声明：本文全体作者均声明不存在利益冲突。

Funding Statement

国家重点研发计划（2017YFE0121200）

Funding Statement

National Key R&D Program of China

References

1. Malde O, Libby J, Moazen M An overview of modelling craniosynostosis using the finite element method. Mol Syndromol. 2019; 10 (1-2):74–82. doi: 10.1159/000490833. [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]

2. Zhang G, Tan H, Qian X, et al A systematic approach to predicting spring force for sagittal craniosynostosis surgery. J Craniofac Surg. 2016; 27 (3):636–643. doi: 10.1097/SCS.0000000000002590. [ PubMed ] [ CrossRef ] [ Google Scholar ]

3. Hukki J, Saarinen P, Kangasniemi M Single suture craniosynostosis: diagnosis and imaging. Front Oral Biol. 2008; 12 :79–90. [ PubMed ] [ Google Scholar ]

4. 沈卫民, 王刚, 吴玉新, 等先天性颅缝早闭的颅骨成形—附 6 例报道实用美容整形外科杂志 2002;(4):15–18. [ Google Scholar ]

5. 冯胜之, 张涤生, 穆雄铮, 等先天性颅缝早闭症的治疗中华整形烧伤外科杂志 1995;(6):406–411. [ Google Scholar ]

6. 贺佳宾. 基于模板匹配的颅缝早闭手术切割方案规划方法研究. 北京: 中国科学院大学, 2020.

7. Borghi A, Rodriguez-Florez N, Rodgers W, et al Spring assisted cranioplasty: A patient specific computational model. Med Eng Phys. 2018; 53 :58–65. doi: 10.1016/j.medengphy.2018.01.001. [ PubMed ] [ CrossRef ] [ Google Scholar ]

8. Wolanski W, Larysz D, Gzik M, et al Modeling and biomechanical analysis of craniosynostosis correction with the use of finite element method. Int J Numer Method Biomed Eng. 2014; 29 (9):916–925. [ PubMed ] [ Google Scholar ]

9. Jiang X, Jia Y, Ning W, et al. Skull mechanics study of PI procedure plan for craniosynostosis correction based on finite element method// 4th International Conference on Bioinformatics and Biomedical Engineering. Chengdu: IEEE, 2010: 1-4.

10. Malde O, Cross C, Lim L, et al Predicting calvarial morphology in sagittal craniosynostosis. Sci Rep. 2020; 10 (1):3. doi: 10.1038/s41598-019-55224-5. [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]

11. Johnson D, Wilkie A Craniosynostosis. Eur J Hum Genet. 2011; 19 (4):369–376. doi: 10.1038/ejhg.2010.235. [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]

12. Heuzé Y, Boyadjiev S A, Marsh J L, et al New insights into the relationship between suture closure and craniofacial dysmorphology in sagittal nonsyndromic craniosynostosis. J Anat. 2010; 217 (2):85–96. [ PMC free article ] [ PubMed ] [ Google Scholar ]

13. Kapp-Simon K A, Speltz M L, Cunningham M L, et al Neurodevelopment of children with single suture craniosynostosis: a review. Childs Nerv Syst. 2007; 23 (3):269–281. doi: 10.1007/s00381-006-0251-z. [ PubMed ] [ CrossRef ] [ Google Scholar ]

14. Salti S, Tombari F, Stefano L D SHOT: Unique signatures of histograms for surface and texture description. Comput Vis Image Underst. 2014; 125 :251–264. doi: 10.1016/j.cviu.2014.04.011. [ CrossRef ] [ Google Scholar ]

15. Rusu R B, Blodow N, Beetz M. Fast point feature histograms (FPFH) for 3D registration// IEEE International Conference on Robotics and Automation. Kobe: IEEE, 2009: 3212-3217.

16. van Veelen M L C, Rooda O H J E, de Jong T, et al Results of early surgery for sagittal suture synostosis: long-term follow-up and the occurrence of raised intracranial pressure. Childs Nerv Syst. 2013; 29 (6):997–1005. doi: 10.1007/s00381-013-2024-9. [ PubMed ] [ CrossRef ] [ Google Scholar ]

17. Vollmer D G, Jane J A, Park T S, et al Variants of sagittal synostosis: strategies for surgical correction. J Neurosurg. 1984; 61 (3):557–562. doi: 10.3171/jns.1984.61.3.0557. [ PubMed ] [ CrossRef ] [ Google Scholar ]

18. Massimi L, Caldarelli M, Tamburrini G, et al Isolated sagittal craniosynostosis: definition, classification, and surgical indications. Childs Nerv Syst. 2012; 28 (9):1311–1317. doi: 10.1007/s00381-012-1834-5. [ PubMed ] [ CrossRef ] [ Google Scholar ]

19. Libby J, Marghoub A, Johnson D, et al Modelling human skull growth: a validated computational model. J Roy Soc Interf. 2017; 14 (130):20170202. doi: 10.1098/rsif.2017.0202. [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]

20. Abbott A H, Netherway D J, Niemann D B, et al CT-determined intracranial volume for a normal population. J Craniof Surg. 2000; 11 (3):211–223. doi: 10.1097/00001665-200011030-00002. [ PubMed ] [ CrossRef ] [ Google Scholar ]

21. 赵波涛. 车辆碰撞事故中颅骨有限元精细建模及碰撞伤害研究. 长春: 吉林大学, 2015.

22. 陈丽丹, 杨斌, 倪健, 等单侧冠状缝早闭症整复手术的三维有限元生物力学分析中华整形外科杂志 2019; 35 (2):117–123. [ Google Scholar ]

23. Swider P, Delanoé F, Jalbert F, et al Mechanical properties of fused sagittal sutures in scaphocephaly. Clin Biomechan. 2021; 86 :105369. doi: 10.1016/j.clinbiomech.2021.105369. [ PubMed ] [ CrossRef ] [ Google Scholar ]

24. 何黎民. 头颈三维有限元模型构建及颅脑直接冲击伤致伤机制研究. 上海: 第二军医大学, 2004.

25. 倪健. 基于三维数字化技术的先天性颅缝早闭症整复外科治疗. 北京: 北京协和医学院, 2017.

26. Li X, Zhu W, He J, et al Application of computer assisted three-dimensional simulation operation and biomechanics analysis in the treatment of sagittal craniosynostosis. J Clin Neurosci. 2017; 44 :323–329. doi: 10.1016/j.jocn.2017.06.032. [ PubMed ] [ CrossRef ] [ Google Scholar ]

27. 刘天甲. 基于 MIMICS 软件的婴幼儿颅缝早闭个性化手术设计研究. 上海: 上海交通大学, 2019.

Articles from Sheng Wu Yi Xue Gong Cheng Xue Za Zhi = Journal of Biomedical Engineering are provided here courtesy of West China Hospital of Sichuan University

针对婴幼儿颅缝早闭的颅骨重塑手术模拟方法

Simulation method of skull remodellingsurgeryfor infant with craniosynostosis

珅宇 陆

杨宇 罗

文键 郑

剑 宫

2. 实验

2.1. 畸形骨瓣区域与理想颅骨模型的配准

2.2. 骨瓣切割与错位拼接虚拟手术

2.3. 有限元模型

表 1

3. 结果

3.1. 颅骨骨瓣切割与错位拼接

3.2. 有限元模拟术后颅骨生长

4. 结论

Funding Statement

Funding Statement

References

珅宇陆

杨宇罗

文键郑

剑宫