用于截骨矫形的CT灰度-材料属性赋值有限元建模方法技术

一种用于截骨矫形的CT灰度‑材料属性赋值有限元建模方法，（1）骨骼三维模型重建，（2）骨骼表面模型的曲面构造与实体化，（3）骨骼实体模型的截骨分割；（4）未发生位置变化的骨骼部件实体模型的网格划分与输出，（5）存在位置变化的骨骼部件实体模型的位置调整与基准参考模型的建立，（6）基于CT图像灰度‑材料属性经验公式的网格模型赋值，（7）骨骼部件网格模型的位置调整与有限元建模前处理，得到最终的骨骼有限元模型。本发明专利技术能够实现骨骼部件模型位置任意调整而同时具备基于灰度赋值的非均匀单元材料属性，具有模型材质高度还原、操作简单和有限元计算结果精确等特点。

【技术实现步骤摘要】


本专利技术涉及骨生物力学有限元分析
，特别是涉及需要对骨骼空间位置或形态学变化的CT灰度赋值有限元建模方法。

技术介绍

准确的材料参数赋值是确保有限元模型计算准确性的关键因素。目前，对于骨骼模型的材料参数赋值主要分为两大类，一类是采用均一赋值方法，另一类是根据骨骼模型的CT图像灰度进行空间配准来完成非均匀赋值。均一赋值的方法是将骨骼组织视为均匀的材料属性，粗略的将具有较高弹性模量的皮质骨和较低弹性模量的松质骨分别赋予一种均一的弹性模量参数，该方法操作简便、快速，模型运算过程中占用的计算机运算资源较少，并且可以灵活进行骨骼模型的各种位置变换操作，如旋转截骨、内外翻截骨以及模拟内固定手术等，完全不受原有二维CT图像的空间限制。然而该方法因简化了材料赋值操作，并不符合正常人体骨骼非均匀材料属性分布的特点，以此所获得的有限元分析结果准确性较低，难以反映骨骼受力后的真实应力应变分布情况。
另一种骨骼有限元建模的材料属性赋值方法为基于CT图像灰度-材料属性关系公式的原位赋值方法。该方法要求骨骼模型的空间坐标位置须和原有CT图像所覆盖的范围保持一致，在对骨骼模型进行网格划分以后，未发生任何位置和形状变化的网格模型重新被导回原有CT图像的范围，如此，现有的三维重建软件可根据CT图像-网格单元进行空间位置的一一配对融合，经过灰度公式转换后的非均匀单元弹性模量数值将被赋予到该空间位置的单元中，使得整个骨骼模型的不同部位的单元根据不同的部位的CT图像灰度获得相应大小的弹性模量，从而使整个骨骼模型的材料属性分布更符合人体骨骼的分布规律。采用该赋值的有限元建模方法所计算得到的结果则具有更高的准确性。然而，这一方法同样存在其局限性，对于需要经过特殊变换操作的骨骼模型如旋转截骨和内外翻截骨，不可避免的使模型本身发生空间位置的变化，骨骼各个部位以新的相对位置关系导入到原有的CT图像范围内的时候，原有的CT图像灰度-骨骼模型网格单元的空间位置关系将不再适用，甚至会出现导入的模型超出原有的CT图像范围的情况。在这种情况下，基于CT灰度-材料属性赋值的有限元方法将不再适用。
为此，国外有学者针对模型位置变化后的参数采用了映射配对的方法，这类方法能有效解决模型变动的问题，能用于内固定手术和截骨手术的模拟。然而这一方法需要自行编写相应的计算机程序和相应的空间映射算法，整个过程需要耗费极大的人力物力，并且大大增加了有限元建模的工作量，而空间映射坐标变化尚不能有效解决网格单元缺失的问题，因此该方法的推广和应用均存在很大的局限性。为此，如何运用现有的商业软件来简单实现既可以进行空间位置变换，又能够基于原有CT图像灰度进行骨骼非均匀材料属性赋值的有限元建模方法，仍是当前骨骼生物力学有限元分析的一个技术难题。
为此，针对这一现状，实现一种可用于截骨矫形的CT灰度-材料属性赋值有限元建模方法以克服现有技术不足甚为必要。

技术实现思路

本专利技术的目的在于解决传统基于CT图像灰度-材料属性赋值的骨骼有限元建模方法不能实现模型的空间位置变换的问题，而该有限元方法建立的骨骼模型既能满足基于CT图像灰度的材料属性非均匀赋值，又可以允许骨骼模型发生相对原有CT图像范围的空间位置变动，方法具有模型材质高度还原、操作简单和有限元计算结果精确等特点。
本专利技术的上述目的通过以下技术措施实现：
提供一种用于截骨矫形的CT灰度-材料属性赋值有限元建模方法，包括如下步骤：
（1）骨骼三维模型重建
采集目标对象骨骼的计算机断层成像数据，将所获得的Dicom格式的断层成像数据导入医学三维重建软件建立骨骼三维模型；
（2）骨骼表面模型的曲面构造与实体化
针对构建的骨骼三维模型，在逆向重建软件中进一步对骨骼三维模型进行表面平滑及去特征化处理，并在骨骼三维模型表面构建曲面片，进而在曲面片基础上构造格栅，从而通过拟合曲面得到骨骼实体模型；
（3）骨骼实体模型的截骨分割
将建立好的骨骼实体模型文件导入CAD软件对骨骼实体模型进行模拟截骨的分割操作，得到多个分割后的骨骼部件模型，以此作为划分单元网格的对象；
（4）未发生位置变化的骨骼部件实体模型的网格划分与输出
将步骤（3）中的未发生位置变化的骨骼部件模型以CAD文件格式输出，并导入到有限元前处理软件中进行网格划分，对未发生位置变化的多个骨骼部件模型分别单独划分网格后并在有限元前处理软件中以网格模型的文件格式输出；
（5）存在位置变化的骨骼部件实体模型的位置调整与基准参考模型的建立
针对步骤（3）中的发生了位置变化的多个骨骼部件模型，在CAD软件中根据临床截骨的标准方法，把发生了位置变化的多个骨骼部件模型进行相应的平移、旋转以及共轴操作，操作过程中分别建立相应的参考平面、参考轴以及参考点，以此建立相应的平面、直线以及空间点参考模型，保存为CAD文件格式；
（6）基于CT图像灰度-材料属性经验公式的网格模型赋值
将步骤（4）中输出的多个骨骼部件的网格模型导入医学三维重建软件中进行材料参数赋值，在软件中分别定义图像CT值HU与表观密度、表观密度与弹性模量的关系公式后，根据分析需要划分材料分布区间数，完成基于CT图像灰度的经验公式赋值过程，最后以有限元模型格式把各个骨骼部件导出并保存；
（7）骨骼部件网格模型的位置调整与有限元建模前处理
将步骤（6）中的各个骨骼部件模型加载至有限元前处理软件中，同时将步骤（5）中的基准参考模型也导入到有限元前处理软件中，在有限元前处理软件中的装配模式下，根据既定的移动距离、旋转角度参数，将骨骼指定的部件模型依照参考面、参考轴以及参考点进行相应距离和角度的移动操作，使各个骨骼部件模型按照预定参数调整到预定的空间位置，得到调整位置后的整个骨骼模型；最后进行模型的载荷、边界条件和接触属性的设置，得到最终的骨骼有限元模型。
优选的，上述步骤（1）中，目标骨骼断层影像数据可通过CT扫描、MRI扫描或micro-CT扫描成像设备获取，医学三维重建软件为Mimics、Simpleware、或者3D-doctor，步骤（1）建立的骨骼三维模型保存文件为STL、IGES或者IGS。
优选的，上述步骤（2）中，骨骼表面模型的曲面构造是在GeomagicStudio软件中完成，模型实体化的操作在CAD软件中完成。
优选的，上述步骤（6）中，具体通过软件Mimics14.0导入多个骨骼部件的网格模型，在软件Mimics14.0中选用材料属性赋值模块完成材料属性赋值经验公式的设定。
本专利技术的一种用于截骨矫形的CT灰度-材料属性赋值有限元建模方法，采用的建模方法，能够实现骨骼部件模型位置任意调整而同时具备基于灰度赋值的非均匀单元材料属性，具有模型材质高度还原、操作简单和有限元计算结果精确等特点。
附图说明
结合附图对本专利技术作进一步的描述，但附图中的内容不构成对本专利技术的任何限制。
图1是本专利技术在实施例2中构建的胫骨三维模型示意图。
图2是本专利技术在实施例2中构建胫骨实体化模型的过程示意图。
图3是本专利技术在实施例2中运用截骨参考平面的方法对胫骨实体模型实施截骨分割示意图。
图4是本专利技术在实施例2中对胫骨分割模型的下半部分依照参考轴进行内旋10°的空间位置变换操作本文档来自技高网...

【技术保护点】
用于截骨矫形的CT灰度‑材料属性赋值有限元建模方法，其特征在于，包括如下步骤：（1）骨骼三维模型重建采集目标对象骨骼的计算机断层成像数据，将所获得的Dicom格式的断层成像数据导入医学三维重建软件建立骨骼三维模型；（2）骨骼表面模型的曲面构造与实体化针对构建的骨骼三维模型，在逆向重建软件中进一步对骨骼三维模型进行表面平滑及去特征化处理，并在骨骼三维模型表面构建曲面片，进而在曲面片基础上构造格栅，从而通过拟合曲面得到骨骼实体模型；（3）骨骼实体模型的截骨分割将建立好的骨骼实体模型文件导入CAD软件对骨骼实体模型进行模拟截骨的分割操作，得到多个分割后的骨骼部件模型，以此作为划分单元网格的对象；（4）未发生位置变化的骨骼部件实体模型的网格划分与输出将步骤（3）中的未发生位置变化的骨骼部件模型以CAD文件格式输出，并导入到有限元前处理软件中进行网格划分，对未发生位置变化的多个骨骼部件模型分别单独划分网格后并在有限元前处理软件中以网格模型的文件格式输出；（5）存在位置变化的骨骼部件实体模型的位置调整与基准参考模型的建立针对步骤（3）中的发生了位置变化的多个骨骼部件模型，在CAD软件中根据临床截骨的标准方法，把发生了位置变化的多个骨骼部件模型进行相应的平移、旋转以及共轴操作，操作过程中分别建立相应的参考平面、参考轴以及参考点，以此建立相应的平面、直线以及空间点参考模型，保存为CAD文件格式；（6）基于CT图像灰度‑材料属性经验公式的网格模型赋值将步骤（4）中输出的多个骨骼部件的网格模型导入医学三维重建软件中进行材料参数赋值，在软件中分别定义图像CT值HU与表观密度、表观密度与弹性模量的关系公式后，根据分析需要划分材料分布区间数，完成基于CT图像灰度的经验公式赋值过程，最后以有限元模型格式把各个骨骼部件导出并保存；（7）骨骼部件网格模型的位置调整与有限元建模前处理将步骤（6）中的各个骨骼部件模型加载至有限元前处理软件中，同时将步骤（5）中的基准参考模型也导入到有限元前处理软件中，在有限元前处理软件中的装配模式下，根据既定的移动距离、旋转角度参数，将骨骼指定的部件模型依照参考面、参考轴以及参考点进行相应距离和角度的移动操作，使各个骨骼部件模型按照预定参数调整到预定的空间位置，得到调整位置后的整个骨骼模型；最后进行模型的载荷、边界条件和接触属性的设置，得到最终的骨骼有限元模型。...

【技术特征摘要】
1.用于截骨矫形的CT灰度-材料属性赋值有限元建模方法，其特征在于，包括如下步骤：（1）骨骼三维模型重建采集目标对象骨骼的计算机断层成像数据，将所获得的Dicom格式的断层成像数据导入医学三维重建软件建立骨骼三维模型；（2）骨骼表面模型的曲面构造与实体化针对构建的骨骼三维模型，在逆向重建软件中进一步对骨骼三维模型进行表面平滑及去特征化处理，并在骨骼三维模型表面构建曲面片，进而在曲面片基础上构造格栅，从而通过拟合曲面得到骨骼实体模型；（3）骨骼实体模型的截骨分割将建立好的骨骼实体模型文件导入CAD软件对骨骼实体模型进行模拟截骨的分割操作，得到多个分割后的骨骼部件模型，以此作为划分单元网格的对象；（4）未发生位置变化的骨骼部件实体模型的网格划分与输出将步骤（3）中的未发生位置变化的骨骼部件模型以CAD文件格式输出，并导入到有限元前处理软件中进行网格划分，对未发生位置变化的多个骨骼部件模型分别单独划分网格后并在有限元前处理软件中以网格模型的文件格式输出；（5）存在位置变化的骨骼部件实体模型的位置调整与基准参考模型的建立针对步骤（3）中的发生了位置变化的多个骨骼部件模型，在CAD软件中根据临床截骨的标准方法，把发生了位置变化的多个骨骼部件模型进行相应的平移、旋转以及共轴操作，操作过程中分别建立相应的参考平面、参考轴以及参考点，以此建立相应的平面、直线以及空间点参考模型，保存为CAD文件格式；（6）基于CT图像灰度-材料属性经验公式的网格模型赋值将步骤（4）中输出的多个骨骼部件的网格模型导入医学三维重建软件中进行材料参数赋值，在软件中分别定义图像CT值HU与表观密度、表观密度与弹性模量的关系公式后，...

【专利技术属性】
技术研发人员：欧阳汉斌，黄文华，谢普生，钱堃，杨洋，黄华军，
申请(专利权)人：南方医科大学，
类型：发明
国别省市：广东;44