髋关节中心的寻找方法、膝关节植入物与流程

本发明涉及髋关节中心的寻找方法、以及其于该方法的膝关节植入物放置位置设计。

背景技术：

由计算机和机械臂辅助的手术是帮助外科医生进行更精确手术的尖端技术。这些技术进步正在彻底改变骨科领域，如全膝关节置换术或部分膝关节置换术。

基于这种机器人辅助技术的手术系统可以根据每个病人的解剖结构的独特形状和运动来设计膝关节植入物放置位置。

植入物的对齐和固定是根据几个参数进行的，其中包括股骨的机械轴。

股骨的机械轴定义为连接髋关节中心和膝关节股骨中心的线，其中髋关节中心定义为一个圆的中心点，该圆在冠状面、矢状面和横向面上与股骨头相匹配。

准确估计髋关节和膝关节股骨中心是使植入物与患者的解剖结构正确对齐的重要一步。

在其他已经上市的先进的膝关节置换技术系统中，髋关节中心的位置是通过围绕髋关节中心旋转股骨来确定的，期间保持髋部(髋臼)的固定。这种类型的旋转生成一组点，这些点是通过使用双目跟踪仪识别安装在股骨上的参考架来收集的。膝关节置换手术系统使用这个参考架来跟踪股骨的位置。

在一个完美的系统中，这些收集点属于一个球体，其半径对应于参考架在股骨上的相应位置，其中心为髋关节中心。但在实际系统中，这些点会受到测量误差和噪声的影响，这些误差和噪声是由股骨运动期间髋关节中心的微小变化引起的。这会导致所找到的点的位置不准确。

技术实现要素：

鉴于上述问题而完成本发明，其目的在于提供一种髋关节中心的寻找方法，能够有效地提高髋关节中心的近似精度。

根据本发明的一方面，提供一种髋关节中心的寻找方法，其由计算机和机械臂辅助进行，包括如下步骤：在保持检测对象的髋部固定的同时，在使该检测对象的下肢围绕髋关节中心进行第一旋转的不同位置处分别使股骨以圆周方式进行第二旋转；使用导航追踪仪识别固定安装在所述股骨上的参考架，在所述计算机中利用所收集到的与所述第二旋转对应的每组点生成空间曲线图中的各个椭圆；确定组成所述各个椭圆的每组点的质心；确定分别最适合所述各个椭圆的每个平面；确定所述每个平面的穿过质心的垂线；用最小二乘法求取一个表示所有垂线中的最近点的解，用所找到的最近点代表髋关节中心。

优选地，所述不同位置是所述下肢以所述髋关节中心为原点而围绕预定中心线进行所述第一旋转时的各个位置，所述预定中心线是连接所述髋关节中心和预定中心的线。

优选地，以完成所述第一旋转一圈的期间为一个测定周期，基于该周期内在所述不同位置处的所述第二旋转来生成所述空间曲线图中的各个椭圆。

优选地，在3～15个的所述不同位置处分别进行1～10次的所述第二旋转。

优选地，通过适当的拟合来求出最适合所述各个椭圆的平面。

优选地，所述导航追踪仪包括双目光学跟踪仪。

根据本发明的另一方面，提供一种膝关节植入物，其放置位置根据包括股骨的机械轴的参数而确定，所述股骨的机械轴定义为连接髋关节中心和膝关节股骨中心的线，其中，所述髋关节中心通过以上所述的髋关节中心的寻找方法来确定。

根据本发明，通过确定穿过各个椭圆质心的每个平面的垂线，可以更为精确地求取髋关节中心的最佳最小二乘近似值，从而有利于进行植入物的对齐和固定。

附图说明

图1示意性地示出了股骨附近的骨骼及安装在上面的参考架的转动模式。

图2示出了股骨盘旋运动而在计算机系统中产生的椭圆信息。

图3示出了求解所有椭圆的法线中最近点的示意图。

具体实施方式

下面结合附图详细描述本发明的示例性实施例。下文描述的和附图示出的示例性实施例旨在教导本发明的原理，使本领域技术人员能够在若干不同环境中和对于若干不同应用实施和使用本发明。因此，本发明的保护范围由所附的权利要求来限定，示例性实施例并不意在、并且不应该被认为是对本发明保护的范围的限制性描述。

为了便于说明，如图1所示，将抱持着检测对象(尤其是患者)的下肢使其(尤其是股骨1)以髋关节中心o1为原点而使其机械轴(连接髋关节中心o1和膝关节股骨中心o2的线)围绕预定中心线(连接髋关节中心o1和某一预定中心oi的线)的转动称为公转re(对应于“第一旋转”)。，而将在某一公转位置bn(例如b，n为任意适当的自然数)处使股骨1以髋关节中心o1为原点而而使其机械轴围绕连接髋关节中心o1和该公转位置bn(例如b)或其附近的位置点的线(下称自转轴)的转动称为自转ron(例如ro)。本领域技术人员可以理解，在此自转并不是使股骨1绕自身的机械轴旋转，而是为了便于区分理解而作的定义(对应于“第二旋转”)。

简而言之，在使股骨1大幅度地公转re一圈的过程中，还分别在任一位置b(或bn)附近使股骨1小幅度地自转一圈。进行自转时的各公转位置bn可以为离散的，例如彼此大致均匀地间隔开，间隔距离越小则近似精度越高。另外，公转re的一圈并不一定是一个严格意义上的圆或椭圆或类似螺旋的包络线(因此在上下文中皆统称为椭圆，但并不限于此)，也不一定要位于横向于预定中心线的同一平面上。自转时的膝关节股骨中心o2的轨迹中心也不一定要严格位于任一位置b(或bn)上。公转、自转的方向也不受限定。

在此，设经过多处(优选3～15个，或5～10个)的自转ron而完成一圈公转re的期间，为一个测定周期。

由此，在根据本发明的系统中，也能在保持髋部2固定的期间通过围绕髋关节中心o1旋转股骨1来确定髋关节中心o1的位置。

与以往技术的不同处包括：在各公转位置处还进行自转，由此取代以往仅对应于各公转位置而生成的点(即所谓收集点)数据，而生成对应于各公转位置bn的小椭圆数据，从而显著提高髋关节中心o1的定位精度。

<导航追踪系统>

例如使患者以仰卧位姿躺在手术台上，在其股骨的解剖结构上进行参考架3的固定安装，例如侵入式。

即，在参考架3上，例如安装有3个以上的反光球4(或靶标)用于追踪定位，还可包括例如三根针脚用以放置在股骨面上，利用钛钉等螺纹件施力以将针脚扎入骨面，将参考架3的例如底座牢牢固定在股骨1上。

通过适当调整参考架3的方位，确定原始配准位置后即可以利用导航追踪仪(例如双目光学跟踪仪)(未示出)进行测定。将导航追踪仪所获得的信息输入计算机，经由数据存储、处理等，而得到对应于参考架3的各移动轨迹位置的空间曲线图信息。

在使用双目跟踪仪识别安装在股骨上的参考架来收集各组点数据的时候，对于股骨与摄像机(未示出)坐标系之间的、参考架与摄像机坐标系之间的、股骨与参考架坐标系之间的基于坐标变换原理的变换关系，以及对于当前位置的股骨坐标系的重构等骨科导航设备，不是本发明的重点，在此不再赘述。

<测定方法>

如图1所示，本发明人提出了一种新的寻找髋关节中心o1的方法，即，在关于参考架3的配准完成之后，在例如用专用夹具固定髋部2(例如使骨盆、髋臼保持不动)的同时，在上述公转re过程中的各位置bn，以大致圆周方式旋转股骨1(自转)几次，从而在计算机软件系统的空间曲线图中生成多个由若干收集点(对应于自转时的各收集位置)围成的小椭圆e，如图2、图3中所示。

使用这些小椭圆e，可以得到一个很好的近似的髋关节中心o。

更具体地，例如，在每一个公转位置，以使股骨自转的方式进行摇动，使用通过双目跟踪仪识别参考架所收集到的点，生成空间曲线图中的由若干小椭圆e(对应于各自转)围成的一个大椭圆f(对应于公转)。

如果是在一个完美的理想系统中，髋关节中心在每次运动中保持固定，没有数据测量误差，构成椭圆e的每个点p都属于一个球体，其半径对应于股骨上的相应位置，中心在髋关节(即，球心对应于髋关节中心)。同样，穿过髋关节中心和椭圆中心的每一条垂直线都将垂直于该点p所属的椭圆e。

然而，在一个真实的系统中，由于患者下肢尤其是骨盆或髋臼会存在晃动、以及多个反光球的加工累计误差等，从椭圆的每个点p到髋关节中心o1的距离可能不是恒定的，穿过椭圆中心的该椭圆的垂直线不一定会穿过髋关节中心。

本发明人研究后认为，在这种情况下，为了提高查找髋关节中心点(hipjointcenter)的精度，通过使用到每个椭圆中心m的垂线l(normalline)，尽管无法确定它们是否穿过髋关节中心，但如果能确保测量误差较小，则这些垂线l将至少非常接近所期望的点(髋关节中心o1)地通过逼近髋关节中心o1的最近点o。该最近点o也可看作是这些法线中的各条法线彼此之间在空间上最邻近的位置。

而且，当这些垂线l分别垂直于每个椭圆e时，可以使用最小二乘法在这一组法线中找到三维空间中的最近点o。

遵循这一逻辑思想，本发明人找到一种最可行的方法。

该方法所基于的算法有五个基本步骤：

1.在保持髋部2固定的同时，在不同的位置以圆周方式旋转股骨1几次。

2.确定组成椭圆e的每组点的质心m。

3.确定最适合每个椭圆e的平面。

4.确定穿过椭圆质心m的每个平面的垂线l。

5.在最小二乘法中找到一个表示所有垂线l中最近点o的解。

所找到的最近点o将代表对患者髋关节中心o1的最佳最小二乘近似值。

在此，该“最近点”o包括在理想的情况下刚好通过髋关节中心o1的点，也视为与下肢的股骨头表面相匹配的球体的中心。

<效果>

为了验证上述算法的方法的有效性，使用软件算法生成了100个不同的数据集。每个数据集由一个代表髋关节中心的点和一组椭圆组成，这些椭圆代表以圆周方式转动股骨时的不同位置。随后，在数据中加入随机噪声，验证了该算法确定髋关节中心位置的有效性，其中，获得了距离髋关节实际中心0.2mm的平均偏差。

由此可知，上述算法可用于膝关节置换术中髋关节中心的估计。

关于以圆周方式旋转股骨几次时的“不同的位置”，可以为在股骨1围绕髋关节中心o1旋转(公转)时机械轴与预定中心线的夹角a为恒定值时的每一个位置，也可以为夹角不是恒定值时的每一个位置。优选地，获得包络成大致环状的多个小椭圆e以提高寻找精度。

根据本发明的上述算法，即使难以实现该夹角a为恒定值的旋转，也能得到上述最佳最小二乘近似值，从而有效地确定最近点。

同样地，对应于公转、自转的采集数据，在空间曲线图中也将生成大致椭圆或大致圆环形(因此在上下文中皆统称为椭圆，但并不限于此)。

每一个椭圆也并不严格在一个平面上，尽管如此，也可通过常规的算法确定组成各个椭圆e的每组点的质心(m)，并通过适当的拟合方法等来求出适合每个椭圆e的平面。

此外，上述100个不同的数据集，对应于100个测定周期，即经过多处的自转而完成100圈公转时的信息。其中，待一个周期的数据检测结束，调用数据处理程序，获取一个周期内的髋关节中心的最近点，存储在计算机中，并显示在显示界面上。依次循环处理、存储和显示。

<标志点>

作为将参考架固定安装在股骨上的方式，包括以下直接、间接的方式。

以上对通过直接、侵入式固定在股骨上的参考架形成标志点进行定位的方式进行了说明，其中参考架的个数可以为一个，或两个以上。

另一方面，作为间接、非侵入式，例如，也可以在下肢皮肤上贴有一定数量(超过三个)的标志点，用来定位下肢的位姿，此时期望这些标志点在下肢上构成一个刚体。如此，也视为固定安装在股骨上。此时，也可以通过平移关节表面标记点的位置来进行测定。由于标志点可以固定在关节附近的任意位置，计算简单并且需要最少数目的标志点。只是，由于皮肤表面曲率受三维人体模型姿态的影响较大，可能会因由该方法计算得到的位置不够稳定而出现误差。其中，例如，在下肢转动时，皮肤不可避免地会发生微小位移，致使下肢上标志点的相对位置发生改变，最终也会导致下肢位姿数据中含有噪声。但最终都可通过包含自转的算法而有效提高了寻找精度。

<优点>

由此，通过根据本发明的上述方法，与仅使下肢并从而仅使股骨进行公转、而在各公转位置处不进行自转的方法相比，可以在空间曲线图中生成沿一个椭圆(对应公转的大椭圆)的周向排列的多个小椭圆e(对应自转的小椭圆)，而不是仅得到沿该大椭圆的周向排列的多个点。

而且，通过确定穿过椭圆质心的每个平面的垂线，可以更为精确地求取髋关节中心的最佳最小二乘近似值。

关于上述公转过程中的自转位置，其个数优选使得所获取的小椭圆彼此相邻接，或更优选地，至少部分彼此重叠。

而且，关于每个小椭圆的生成，也可以重复自转1次之上，将所获得平均的结果作为所需要的小椭圆。

<植入物>

本文所描述的植入物可包括用于膝关节置换的植入部和假体。该植入部可用于在内部固定骨折或受损的破损部分。该假体可用于替代骨头的整体或一部分。

基于根据本发明所确定的髋关节中心，可用于生成包括股骨的机械轴的股骨几何参数，从而可用于设计适于该患者或该患者所属群体的膝关节置换植入物的放置位置，有利于进行植入物的对齐和固定。

如此，通过计算机辅助的骨科手术系统，可以精确地将髋关节中心的最近点显示在屏幕上，从而可以根据每个患者的解剖结构的独特形状来对齐放置膝关节植入物，提高手术精度，减少术中创伤。

根据本发明的上述方法，适用于计算机和机械臂辅助的全膝关节置换术或部分膝关节置换术(统称膝关节置换术)，某些场景中，还可以省略使用ct等医疗设备进行关节置换前检查。对于典型的假体植入手术，可以精确设计术中假体的安放位置，降低对医生经验控制的要求，减少人为因素造成的误差来源。有利于假体位置的准确选择，手术精度的提高，以及术后的恢复顺利。

在本发明中，用于寻找髋关节中心的系统的各单元可由相应的硬件或软件单元实现，各单元可以为独立的软、硬件单元，也可以集成为一个软、硬件单元，在此不用以限制本发明。各单元的具体实施方式可基于本领域技术人员的知识而实现，在此不再赘述。

在上述说明中，将股骨的机械轴定义为连接髋关节中心o1和膝关节股骨中心o2的线。这并不影响关于对下肢机械轴的定义(通过髋关节中心、膝关节股骨中心和踝关节中心的轴线)，反而有利于下肢机械轴的确定。

而且，本发明的用于寻找髋关节中心的方法，也能够适用于寻找膝关节股骨中心。

应当指出，上面的说明仅是示例性的，本领域技术人员可以根据上述说明对本发明实施例做出各种修改和变型，这些修改和变型均在本发明的保护范围之内。