股骨头旋转中心计算方法与系统

技术领域

本发明涉及手术导航系统技术领域，特别是涉及一种股骨头旋转中心计算方法与系统。

背景技术

全膝关节置换术(Total Knee Replacement，TKR)是解决影响患者活动功能的膝关节严重病损的最有效的手术之一，能够有效提高患者的生活质量。但是目前传统的人工膝关节置换手术非常依赖于临床医生的经验，由于假体松动、脱位、断裂和感染噪声的失败率达5％～8％。相比于传统的人工膝关节置换手术，膝关节置换手术导航系统将计算机图像处理与机器人精准规划结合起来，对受力与运动分析进行智能化评估，辅助临床医生完成膝关节置换手术，能够减少手术损伤，减少手术时间，提高手术成功率与质量。

全膝关节置换术的目标是手术后患者的股骨力线与胫骨力线完全重合，所以膝关节置换手术导航系统需要计算得到股骨头中心从而获取股骨力线来评价手术是否成功，因此，研究股骨头旋转中心计算方法与系统意义重大。

发明内容

本发明要解决的技术问题是提供一种股骨头旋转中心计算方法与系统，能够准确计算得到股骨头旋转中心点，简单、方便、鲁棒性好。

为解决上述技术问题，本发明提供了一种股骨头旋转中心计算方法，所述方法包括：记录摇动大腿时股骨头上多个标志物的空间位置数据；使用随机抽样一致算法删除空间位置数据的奇异点；采用高斯消元法拟合得到空间位置数据的初始球心坐标；计算所有标志物的球心坐标的平均值；将球心的平均值作为初始值，以柯西核函数为鲁棒核函数，以删除奇异点后的所有跟踪标志物的信息为基础，建立球形方程，进行非线性迭代优化，得到球心坐标。

在一些实施方式中，采用高斯消元法拟合得到空间位置数据的初始球心坐标，包括：根据构造的函数，利用高斯消元法拟合得到空间位置数据的初始球心坐标。

在一些实施方式中，构造的函数具有如下的形式：

其中，x

在一些实施方式中，根据构造的函数，利用高斯消元法拟合得到空间位置数据的初始球心坐标，包括：通过求解由构造的函数推导得到的方程组，拟合得到空间位置数据的初始球心坐标。

在一些实施方式中，方程组具有如下形式：

其中，(u

此外，本发明还提供了一种股骨头旋转中心计算系统，所述系统包括：一个或多个处理器；存储装置，用于存储一个或多个程序，当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现根据前文所述的股骨头旋转中心计算方法。

采用这样的设计后，本发明至少具有以下优点：

通过高斯消元能够快速获取球心坐标的大概位置，再通过非线性的迭代方法快速、准确计算得到球心坐标，避免非线性迭代进入局部最优的缺陷。

附图说明

上述仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，以下结合附图与具体实施方式对本发明作进一步的详细说明。

图1是股骨头旋转中心计算的流程图。

具体实施方式

以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。

本发明提出一种股骨头中心计算方法与系统，该方法能够准确计算得到股骨头旋转中心点，简单、方便、鲁棒性好。

股骨头中心的确认是股骨力线确认的基础，也是膝关节导航系统术中假体摆放规划的前提。本发明提出一种膝关节导航系统的股骨头中心计算方法与系统。

在股骨上安装跟踪标志物支架，在固定髋关节不动的情况下，完成摇大腿的操作，记录跟踪标志物支架中每一个跟踪标志物的空间位置。

首先，对于每一个标志物的位置，使用随机抽样一致(RANSAC)算法计算标志物不同集合的方差，删除奇异点。

然后对于每一个跟踪标志物的点集，构造一个函数

目标是计算得到(x

采用高斯消元法拟合得到球心坐标。设：

简化得到

解方程得到(u

再计算所有跟踪标志物的球心的平均值。

将球心的平均值作为初始值，以柯西(CauchyLoss)核函数为鲁棒核函数，以删除奇异点后的所有跟踪标志物的信息为基础，建立球形方程，进行非线性迭代优化，得到球心。

建立的计算球心非线性最小二乘方程为：

柯西(CauchyLoss)核函数表达式为：ρ(s)＝log(1+s)

本发明的优势在于，通过高斯消元能够快速获取球心坐标的大概位置，再通过非线性的迭代方法快速、准确计算得到球心坐标，避免非线性迭代进入局部最优的缺陷。

以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，本领域技术人员利用上述揭示的技术内容做出些许简单修改、等同变化或修饰，均落在本发明的保护范围内。