发明内容

本发明的目的在于提供一种人工膝关节精准置换的定位块制造方法及系统，以解决上述问题。

为实现上述目的，本发明采用以下技术方案：

一种人工膝关节精准置换的定位块制造方法，包括以下步骤：

步骤1，获取影像学数据，采集下肢全长的CT数据，通过断层数据处理软件建立下肢骨的三维几何模型，导出STL文件；

步骤2：获取人体运动学数据，通过运动捕捉系统来采集步态运动信息和足底反作用力信息，将光学摄像机和测力板采集的参数以C3D格式文件获得；

步骤3：根据步骤1和步骤2中下肢骨的三维几何模型和人体运动学数据，基于轮廓特征拾取技术、个体化骨骼肌肉缩放技术和逆向运动学分析技术，建立个体化骨肌多体动力学模型；

步骤4：在个体化骨肌多体动力学模型模拟定位截骨安装，在动态步态运动仿真过程中分析人工膝关节的不同定位位置所产生的生物力学和运动学响应，获得假体安装位置的安全范围，确定手术定位块定位孔的设计位置；

步骤5：基于人体下肢骨骼的统计模型库，建立了通用的下肢骨骼几何三维模型，针对通用模型的膝关节设计参数化的定位块，通过参数优化或调整实现定位块的个体化建模；

步骤6：基于步骤1中建立的个体化患者下肢骨几何模型，利用形状统计分析方法缩放通用模型，根据缩放后的通用模型几何形状的相关参数调整步骤5中设计的参数化定位块的定位面以得到与目标膝关节几何形状相匹配的定位块；

步骤7：根据步骤4获得的定位块定位孔的设计位置，进一步对步骤6中获得个体化定位块的定位孔进行参数化调整；

步骤8：个体化定位块的3D打印制造。

进一步的，在步骤1中，CT或MRI数据以DICOM格式文件保存。

进一步的，步骤1中采集对象采用平躺姿态，扫描方向从头到脚，扫描层厚建议0.625mm，同时记录扫描设备型号、扫描电压、扫描电流参数。

进一步的，在步骤2中，采用Vicon三维步态分析系统进行人体行走过程中的运动学数据采集和处理；在人体运动采集过程中，在人体骨骼解剖特征的皮肤表面粘贴反光Marker点，利用10台反光摄像机(采集频率100Hz)；Vicon三维步态分析系统捕捉每个Marker点的三维坐标，测量足底反作用力，并通过MX Giganet数据转换控制台用以实现信号同步和连接到PC电脑。

进一步的，具体的：Vicon三维步态分析系统包括反光Marker点、10台光学反光摄像机、MX Giganet数据转换控制台、电脑、MXT型标定架、三块AMTI Gen5/Optima测力台和Vicon Nexus 1.8.5软件；人体运动采集过程中，在人体骨骼解剖特征的皮肤表面粘贴反光Marker点，利用10台反光摄像机捕捉每个Marker点的三维坐标，三块AMTI Gen5/Optima测力台测量足底反作用力，并通过MX Giganet数据转换控制台用以实现信号同步和连接到PC电脑；利用Vicon Nexus 1.8.5软件可以导出C3D格式的文件进行数据储存。

进一步的，在步骤3中，获取人体下肢骨骼模型之后，采用先进的骨变形缩放技术将通用模型的骨骼几何缩放成个体化下肢骨形状，同时定义在骨骼表面的肌肉附着点被利用一个非线性径向基函数RBF的插值框架缩放；定义一个仿射转换来将最初的统计学骨骼模型缩放成个体化的骨骼；基于上述仿射转换定义一个3角谐波RBF缩放方程；在每个骨骼映射缩放时，在骨骼周围定义了一个边界框；最后，基于先前使用的骨标点定义一个反向刚性转换，将缩放后的骨骼几何从患者数据的参考坐标系下转换到Anybody骨肌系统的参考坐标系下；缩放所用的RBF缩放方程如下：

φ(r)＝r³

其中，cj是RBF函数φ的系数，是根据目标骨和被缩放骨的骨标计算而来的；q是多项式p的阶数；

在模型骨骼体节缩放的同时，附着在骨骼上肌肉的附着点和几何形状也根据一个线性的几何缩放定理被缩放：

s＝Sp+t

其中，s为被缩放后的点，S为缩放矩阵，p为初始点，t是空间偏移量；

长度参数控制骨骼轴向尺寸的缩放，其他方向根据体重和脂肪参数按比例缩放。脂肪比例根据公式计算：

BMI＝BodyMass/(bodyheight^2)

FatPercent＝(-0.09+0.0149*BMI-0.00009*BMI^2)*100。

进一步的，在步骤3、4中，在Anybody软件中根据个体患者的个体化下肢骨的三维几何特征、个体化步态运动学信息和人工膝关节置换的手术规划，建立人工膝关节置换的个体化骨肌多体动力学模型，基于逆向运动学分析技术求解人体关节的运动学，基于逆向动力学分析技术求解人体关节的接触力学和周围肌肉的作用。

进一步的，一种人工膝关节精准置换的定位块制造系统，包括：

影像学数据获取模块，用于获取影像学数据，采集下肢全长的CT数据，通过断层数据处理软件建立下肢骨的三维几何模型，导出STL文件；

人体运动学数据获取模块，用于获取人体运动学数据，通过运动捕捉系统来采集步态运动信息和足底反作用力信息，将光学摄像机和测力板采集的参数以C3D格式文件获得；

个体化骨肌多体动力学模型建立模块，用于根据步骤1和步骤2中下肢骨的三维几何模型和人体运动学数据，基于轮廓特征拾取技术、个体化骨骼肌肉缩放技术和逆向运动学分析技术，建立个体化骨肌多体动力学模型；

假体安装位置获得模块，用于在个体化骨肌多体动力学模型模拟定位截骨安装，在动态步态运动仿真过程中分析人工膝关节的不同定位位置所产生的生物力学和运动学响应，获得假体安装位置的安全范围，确定手术定位块定位孔的设计位置；

参数优化模块，用于基于人体下肢骨骼的统计模型库，建立了通用的下肢骨骼几何三维模型，针对通用模型的膝关节设计参数化的定位块，通过参数优化或调整实现定位块的个体化建模；

定位块获取模块，用于基于步骤1中建立的个体化患者下肢骨几何模型，利用形状统计分析方法缩放通用模型，根据缩放后的通用模型几何形状的相关参数调整步骤5中设计的参数化定位块的定位面以得到与目标膝关节几何形状相匹配的定位块；

参数二次优化模块，根据步骤4获得的定位块定位孔的设计位置，进一步对步骤6中获得个体化定位块的定位孔进行参数化调整；

打印模块，用于个体化定位块的3D打印制造。

与现有技术相比，本发明有以下技术效果：

1.本发明揭示了患者因素与定位块设计之间的内在映射关系，基于患者的个体化影像学数据和步态运动学数据，利用骨肌多体动力学仿真技术，通过人工膝关节的生物力学和运动学响应，明确了个体化定位块定位孔的设计位置，为人工膝关节精准置换提供了途径，为个体化定位块的精准设计提供了手段。基于本发明设计的个体化定位块，针对性强，精度更高，设计依据更明确。

2.本发明基于所建立的参数化定位块模型，针对个体化患者下肢骨几何模型和仿真计算获得的定位块定位孔的设计位置，通过参数调整或优化能够快速实现个体化定位块的设计，大大缩短了定位块设计的周期，提高了定位块设计的时效性。

3.本发明是依据个体化患者的影像学、运动学和体内关节生物力学信息而建立的，个体化定位块的界面匹配性和定位孔的精准性大大降低了定位误差带来的患者人工膝关节失效风险。以此来计算出膝关节假体最合适的安装定位位置，从而设计出相对应的个体化定位块，相对于传统的定位块的固定位置，本发明的个体化定位块更适合病人的生理特征，减少病人在日常运动中产生的不适感和受力不均等现象，避免了术后假体出现偏移需要再次进行手术修复的情况，极大的减少了术后并发症的出现。