亚星游戏官网-www.yaxin868.com

专利名称：用于跟踪股骨参照系的基于mems的方法和系统的制作方法
技术领域：
本申请涉及计算机辅助的利用微电子机械传感器(MEMS)的整形外科手木。
背景技术：
利用MEMS传感器导航骨骼的关键步骤之ー是相对于所述传感器定位骨骼。对于股骨而言，传感器相对于横轴的取向会受到机械限制，例如，使抓手插入在后髁部下方，以 使传感器横轴与骨骼的横轴对齐。然而，传感器相对于股骨机械轴的取向更加复杂，因为定义该轴的解剖学界标之一，股骨头，隐藏在臀部内部。在光学导航系统中，股骨头是通过移动股骨、并假设骨盆是稳定的、找到骨骼绕其转动的固定枢轴点来定位的。这依赖于具有六个自由度(DOF)的光学传感器，所述六个自由度即，位置的3D0F和取向的3D0F。然而，在MEMS系统中，传感器不自动提供6D0F。6D0F能够通过对陀螺仪和加速度计读数求积分来得到——一种被称作“航位推算”的标准技术——但是这种技术对传感器误差非常敏感并且因而不适用于低成本传感器。存在用于基于轴旋转动力学获取股骨的轴的其他基于陀螺仪的方法。然而，这样的方法需要非常特定的并且不自然的腿运动，其可能难以应用并且受限。

发明内容
因此，本公开的目的是提供一种用于从基于MEMS的数据来定义和跟踪股骨坐标系的新型方法和系统。因此，根据本申请的第一实施例，提供了一种用于在计算机辅助外科手术中跟踪股骨參照系的系统，该系统包括传感器単元，其适于被固定到股骨上，所述传感器単元包括加速度计传感器和陀螺仪传感器，二者均产生取向数据；处理单元，其接收陀螺仪数据和加速度计数据并且包括陀螺仪数据计算器，其用于从至少由于股骨的运动而产生的陀螺仪数据来提供计算的加速度数据，加速度计数据计算器，其用于从由于股骨的运动而产生的加速度计数据来计算测量的加速度数据，以及加速度比较器，其用于从计算的加速度数据与测量的加速度数据之间的比较来将传感器単元的取向与股骨相关联以定义股骨參照系；传感器取向接ロ，其通过跟踪所述股骨參照系来提供针对股骨的取向数据。仍根据第一实施例，所述陀螺仪数据计算器利用针对股骨的估计的轨�：桶刖吨道刺峁┘扑愕募铀俣仁�据。仍根据第一实施例，所述陀螺仪数据计算器利用针对股骨的估计的半径值从所述陀螺仪数据推导向心加速度，以及利用陀螺仪传感器的取向从所述加速度计数据推导重力，来提供计算的加速度数据。仍又根据第一实施例，由所述加速度比较器定义的參照系是三轴股骨坐标系。仍又根据第一实施例，所述比较器对所述比较应用数值优化以定义股骨參照系。仍又根据第一实施例，所述加速度比较器计算股骨长度和股骨轨迹中的至少ー项。仍又根据第一实施例，所述系统包括位于传感器単元和处理单元中的任一个中的滤波器，用于减少来自所述陀螺仪数据和所述加速度计数据中的任一个的噪声。根据本申请的第二实施例，提供了一种用于利用包括被固定到股骨上的加速度计传感器和陀螺仪传感器的传感器单元来跟踪股骨參照系的方法，该方法包括获取由于股 骨的运动而产生的陀螺仪读数和加速度计读数；至少从所述陀螺仪读数计算加速度数据；从所述加速度计读数測量加速度数据；比较计算的加速度数据和测量的加速度数据；利用所述比较将传感器単元的取向与股骨相关联以定义股骨參照系；并且从所述陀螺仪读数以及从所述加速度计读数跟踪所述股骨參照系。仍根据第二实施例，计算加速度数据包括利用针对股骨的轨迹和股骨的半径的估计值。仍又根据第二实施例，计算加速度数据包括利用针对股骨的半径的估计值从所述陀螺仪数据推导向心加速度，以及获得所述陀螺仪的取向以从加速度计读数推导重力。仍又根据第二实施例，将取向进行关联包括利用所述比较执行数值优化以获得所述传感器単元相对于股骨的取向。仍又根据第二实施例，所述方法包括，从所述參照系计算股骨的半径和股骨的轨迹中的至少ー项。仍又根据第二实施例，定义所述股骨參照系包括定义三轴坐标系。仍又根据第二实施例，所述方法在骨骼模型或尸体上执行。


图I是图示说明用于关于基于MEMS的传感器単元定义股骨坐标系的根据本申请的实施例的方法的流程图；以及图2是图示说明用于关于基于MEMS的传感器単元定义股骨坐标系的根据本申请的另ー实施例的系统的框图。
具体实施例方式本文使用MEMS指代微电子机械传感器，诸如加速度计和陀螺仪。本文使用CAS指代计算机辅助外科手木。參考附图并具体參考图1，在I中图示了用于关于基于MEMS的传感器単元定义股骨坐标系的方法，用于连续跟踪股骨坐标系，或其他股骨參照系(例如，股骨的机械轴)。作为非限制性范例，利用下文所述的方法定义的坐标系可以用于执行涉及股骨的连续操作，例如，如在美国专利申请No. 12/410884和相应的PCT申请公布No. WO 2009/117833中所描述的，在此通过弓I入将其并入本文。利用方法I获得的股骨坐标系可以被用在计算机辅助外科手术(CAS)中，用在执行对股骨、对骨盆和/或对胫骨的改变中，以及用在为操作者提供诸如肋骨长度数据的外科手术參数中。因此，可以使用方法I来找到股骨的长度，即，股骨头的中心与机械轴进入点之间的距离，以及从基于MEMS的传感器単元10的取向来跟踪所述股骨坐标系(图2)。基于MEMS的传感器単元具有陀螺仪传感器12和加速度传感器14两者(图2)，并被连接到CAS系统以提供来自这两种类型的传感器的读数，根据所述读数，股骨坐标系被定义并被连续跟踪。陀螺仪传感器12和加速度计传感器14中的每个都至少提供沿三个自由度的取向数据。根据步骤2，传感器单元10被固定到股骨。在实施例中，传感器单元10被固定在膝盖处的机械进入点处。当被固定在股骨的膝盖端时，假设传感器単元10的运动处在以股骨头的中心为中心、半径等于股骨的长度的球上。然后如下所述地找到匹配陀螺仪读数和加速度计读数的最佳球。要获得的两个要素是传感器単元10相对于股骨坐标系的取向和球的半径。如果股骨与传感器坐标系之间的变换由——分别围绕x、y和z轴的三个相继的旋转一入i、入2 和X3来表示，那么需要优化四个參数。然而，当把传感器単元10固定到股骨后，传感器单元10可以与股骨后髁对齐以限制旋转，使得传感器単元10的X轴当被投影到股骨横平面上时与股骨的X轴匹配一 X3恒为零并且仅需要优化三个參数(半径、入i和入2)。也考虑了传感器単元10与股骨之间的其他关系(例如，使用其他解剖学界标)，只要已知传感器单元10相对于股骨的机械进入点的位置和取向。根据图I中的步骤3，执行股骨的特定运动。基于MEMS的传感器単元10的陀螺仪传感器12以角速度的形式给出关于取向变化的信息。所述特定的运动是股骨旋转的停止序列，以消除重力对测量的加速度的影响，如下文中所描述的。为了计算在轨迹上方世界中的绝对传感器取向，必须已知初始取向。如果传感器単元10是稳定的，能够从加速度计传感器14的读数计算炉、(和0。￠9、(和0分别是X-轴、Y-轴和Z-轴与地面的夹角。前面所提及的旋转的停止序列允许对陀螺仪的f、^和0角进行重新初始化。如果假设初始偏航(Vtl)为0，那么可以计算传感器単元10在世界中的初始取向传■器。根据步骤4，从陀螺仪传感器读数以及从传感器単元10的初始取向(例如，从稳定时的加速度传感器14获得的)计算股骨的加速度。根据陀螺仪传感器12的特性，沿三个传感器轴的角偏移u*0对应于关于单位矢量u的0弧度的持续旋转。这也被应用于角速度。可以利用如下方程在间隔期间提取角速度以及绕其旋转的轴 t= I (gx, gy, gz) t
「 .(Sx , g' ' SzVt =--。
Cl),如果假设在时间间隔t内为常数，那么在该时间间隔期间的角偏移为Ot = + 1 At O如果将以轴和角度的形式表示的这ー旋转被转换为旋转矩阵Rt，那么利用如下方程计算传感器单元10的新的取向イ_器纖,中，=式 传感器卿冲ノ通过转换，假设质点(股骨上的參考点)处在机械轴进入点的位置处并且与股骨坐标系对齐。传感器单兀10在质点坐标系(传感器在帛ま巾)中的取向直接对应于由入i和入2定义的一机械轴方向并且对于所有轨迹是恒定的。其可以由如下矩阵表示
cos(A|) sin(sin(ノぃ)cos(/L) * sin(A,)
传感器顏点中= 0COs(A1)— Sin(A2) 。
-sin(,しi) cos(A ) * sin(/l1 j cos(a ) * cos(a:)如果通过将传感器単元10与后髁部对齐来限制旋转，那么应用这种变换，而如果·其通过例如将传感器単元10与前后轴对齐来限制，那么应当对其进行修正。在已知传感器单元10与质点之间的关系的情况下，可以利用下式从传感器単元
的取向直接计算质点在世界坐标系中的取向
质リ、_料=传感器錄界巾,*质ム丨、細*器中，
质点爾^,=传感器,,書*质点/細削ぐ1。按照定义，股骨机械轴对应于股骨坐标系的z轴，SP :
V
MecMx、师,丨,，=质点餓别ゝ* 0，
I其中，MeekAxismm^3xI 的向量。如果假设股骨头的中心为世界的中心，S卩(0，0，0)，那么质点的坐标可以由下式给出
Pt = —4^ ■在坦界中。利用运动的方程计算质点的计算的加速度
r-2 *( , - /).丨)Vt = ~^~- V1^，
V, - V,,at =~
At其中，at是切向加速度(或线性加速度[即，假设向心加速度不是线性加速度的一部分])。为了计算这些值，提供针对股骨半径和针对初始传感器轨迹的估计值。假设初始的速度(Vtl)和加速度( )为O。根据步骤5，将从陀螺仪传感器2测量的质点加速度与加速度计传感器14的读数进行比较，以确定误差(即，代价(cost))。候选解决方案的代价是基于通过加速度计观测到的加速度(a-)和从质点运动计算的那些加速度(计算的加速度)之间的误差的。由传感器单元10感测到的加速度包括质点的加速度和重力加速度,从而利用股骨执行运动的停止序列来将重力加速度从其他加速度中区分出来。利用传感器単元10的计算的取向，通过如下地将这两个加速度投影到传感器轴上来获得预期的加速度(aexp)
ちexp, = (0,0, G) 传感器韻则ゝ+Oノ传感器.界,丨
flvexp, =(OAG) 传感器,,鎌种，+ a,,
fifZexpr =の,0，G) 传感器2鍾界中,+a #传感器z在肚界中，
·
其中，G是重力常数。代价被定义为计算的加速度与测量的加速度之间的差的最小平方误差，如以下方程所描述的
n-At --
代价=(A Obs,- a^xplT + obS, —ay exp, T +(az obSra- exp,)"，其中，n是轨迹的持续时间，并且A t是两次采集之间的时间间隔。因为ー些数据，或者数据的一部分，更倾向于是准确的或者有区分性的，所以对代价应用权重因子。因为取向的误差(飘移)随时间増加，对处在轨迹的开始时的误差给予更
多的权重。代价方程因此变为
代价=Zi* ホaX obSt ~ aX exp, )' + (ay ObSt — aV exp, ) + ( - obh ~ °z exp,)"。改进代价函数的另ー种方式是不对应现实的或可能的设置，或者不匹配补偿信息。在该方法中，对不像是或者不可能的半径或机械轴取向进行惩罚。如果半径低于或高于定义的阈值(例如，低于40cm或高于60cm)，施加与差成比例的惩罚。将惩罚添加到总的代价并且如下表述如果半ィ5>半ィ5 max—代价半径=(半ィ5 max-半ィ5) *k半径，如果半径く半径min —代价半彳5 =半径-半径J其中，k￥e是半径惩罚常数。如果X1或X2高于定义的阈值，那么施加同样的惩罚。作为用于计算切向加速度(步骤4)并将这一备选切向加速度与加速度计传感器14的读数进行比较的备选方法，可以在不使用如上所述的运动方程的情况下计算质点的切向加速度at，并且因而不依赖于采集时间。这ー备选计算方法和代价函数源于这样的原理，即，当质点绕固定的旋转中心旋转时，感测到的加速度来源于三种不同的源。ー个第一源是向心加速度，^。向心加速度是当前角速度和半径的函数，也就是将质点保持在其球面轨迹上所需的加速度。第二源是切向加速度，^。这是随时间改变质点的角速度所需的加速度。第三源是重力加速度,^，尽管不是影响质点运动的实际的力(如果是这样，其将不会保持在球形轨迹上)，但是其还是被传感器14感测到。这取决于G和当前设备倾斜的角度。根据备选步骤4和5，需要得到X1和X2的值，如果其被应用到传感器単元10，那么将会对传感器单元10进行取向以使Z轴与机械轴对齐。因此
权利要求
1.一种用于在计算机辅助外科手术中跟踪股骨參照系的系统，包括 传感器単元，其适于被固定到股骨上，所述传感器単元包括加速度计传感器和陀螺仪传感器，二者均产生取向数据； 处理单元，其接收陀螺仪数据和加速度计数据，并且包括 陀螺仪数据计算器，其用于至少从由于所述股骨的运动而产生的所述陀螺仪数据来提供计算的加速度数据； 加速度计数据计算器，其用于从由于所述股骨的运动而产生的所述加速度计数据来计算测量的加速度数据；以及 加速度比较器，其用于从所述计算的加速度数据与所述测量的加速度数据之间的所述比较将所述传感器単元的取向与所述股骨相关联，以定义股骨參照系； 传感器取向接ロ，其用于从对所述股骨參照系的跟踪来提供针对所述股骨的取向数据。
2.根据权利要求I所述的系统，其中，所述陀螺仪数据计算器利用针对所述股骨的估计的轨�：桶刖吨道刺峁┧�述计算的加速度数据。
3.根据权利要求I所述的系统，其中，所述陀螺仪数据计算器通过利用针对所述股骨的估计半径值从所述陀螺仪数据推导向心加速度，以及利用所述陀螺仪传感器的取向从所述加速度计数据推导重力，来提供所述计算的加速度数据。
4.根据权利要求I至3中的任一项所述的系统，其中，由所述加速度比较器定义的所述參照系是三轴股骨坐标系。
5.根据权利要求I至4中的任一项所述的系统，其中，所述加速度比较器对所述比较应用数值优化，以定义所述股骨參照系。
6.根据权利要求I至5所述的系统，其中，所述加速度比较器计算股骨长度和股骨轨迹中的至少ー项。
7.根据权利要求I至6所述的系统，还包括位于所述传感器単元和所述处理単元中的任ー个中的滤波器，用于减少来自所述陀螺仪数据和所述加速度计数据中的任一个的噪声。
8.一种用于利用包括加速度计传感器和陀螺仪传感器的被固定到股骨上的传感器单元来跟踪股骨參照系的方法，包括 获得由于所述股骨的运动而产生的陀螺仪读数和加速度计读数； 至少从所述陀螺仪读数计算加速度数据； 从所述加速度计读数測量加速度数据； 比较计算的加速度数据和测量的加速度数据； 利用所述比较将所述传感器単元的取向与所述股骨相关联以定义股骨參照系；以及 从所述陀螺仪读数以及从所述加速度计读数跟踪所述股骨參照系。
9.根据权利要求8所述的方法，其中，计算加速度数据包括利用针对所述股骨的轨�：退�述股骨的半径的估计值。
10.根据权利要求8所述的方法，其中，计算加速度数据包括利用针对所述股骨的所述半径的估计值从所述陀螺仪读数推导向心加速度，以及获得所述陀螺仪的取向以从所述加速度计读数推导重力。
11.根据权利要求8至10中的任一项所述的方法，其中，对取向进行关联包括利用所述比较执行数值优化以获得所述传感器単元相对于所述股骨的取向。
12.根据权利要求8至11中的任一项所述的方法，还包括，从所述參照系计算所述股骨的半径和所述股骨的轨迹中的至少ー项。
13.根据权利要求8至12中的任一项所述的方法，其中，定义所述股骨參照系包括定义三轴坐标系。
14.根据权利要求8至13中的任一项所述的方法，其中，所述方法是在骨骼模型或尸体上执行的。
全文摘要
一种用于在计算机辅助外科手术中跟踪股骨参照系的系统，其包括传感器单元。所述传感器单元适于被固定到股骨上。所述传感器单元包括产生取向数据的加速度计传感器和陀螺仪传感器。处理单元接收陀螺仪数据和加速度计数据。所述处理单元包括陀螺仪数据计算器，其用于提供由于股骨的运动而产生的计算的加速度数据；加速度计数据计算器，其用于计算由于股骨的运动而产生的测量的加速度数据；以及加速度比较器，其用于将传感器单元的取向与股骨相关联以定义股骨参照系。所述股骨参照系是从计算的加速度数据与测量的加速度数据之间的比较定义的。传感器取向接口通过从对所述股骨参照系的跟踪提供针对股骨的取向数据。还提供了一种用于跟踪股骨参照系的方法。
文档编号G01P15/18GK102843988SQ201080064991
公开日2012年12月26日 申请日期2010年7月27日 优先权日2010年3月2日
发明者M·瓦林, C·普罗克斯, L-P·阿米奥, B·法拉迪约, J·德朗尚 申请人:奥尔索夫特公司