
为了理解人体运动,生物力学研究者致力于以完美的精度追踪骨骼,将每块骨骼都视为受运动定律支配的刚体。然而,一个巨大的障碍阻碍了这一理想的实现:我们无法直接看到骨骼。相反,我们追踪放置在皮肤上的标记点,但不幸的是,皮肤的运动与骨骼并非完全同步。这种差异导致了一个普遍存在的问题,即软组织伪影(STA),它是运动分析这台“机器中的幽灵”,系统性地破坏我们的数据,并挑战我们结论的有效性。这种误差不是简单的噪声,而是一种复杂的、结构化的信号,可能导致对关节功能和肌肉力量的错误解读。
本文全面概述了软组织伪影,引导读者从其基本性质到为克服它而开发的复杂方法。在“原理与机制”部分,我们将剖析STA是什么,其欺骗性如何导致运动学和动力学计算中的系统性误差,以及测量这种难以捉摸的伪影所涉及的挑战。随后,“应用与跨学科联系”部分将展示STA在临床诊断和体育表现等关键领域的实际影响,并揭示同样的基本问题如何在看似无关的学科中出现和解决,从而展示精确测量的普适性原理。
为了理解我们如何运动,生物力学研究者梦想着一个可以完美精确地追踪身体中每一块骨骼的世界。毕竟,骨骼是一件宏伟的机械工程杰作——一个由杠杆和支柱组成的系统。分析这个系统的基本出发点是刚体假设:即每块骨骼都是一个不可改变、不可变形的物体。如果我们能够追踪这些刚性骨骼,牛顿发现的运动定律将使我们能够清晰地计算出肌肉产生的力和力矩。
但是,正如科学中常有的情况一样,混乱的现实阻碍了我们优雅的梦想。我们无法直接看到骨骼。相反,我们把标记点放在皮肤上,用高速摄像机或惯性传感器来追踪它们。问题就在于此:皮肤不是骨骼。
想象一下,试图通过追踪粘在保龄球表面的果冻块来了解球的精确轨迹。当球沿着球道滚动时,果冻会摇晃、拉伸和抖动。你记录下的运动将是球的真实路径和果冻自身“舞蹈”的结合。果冻相对于球的这种不必要的额外运动,本质上就是软组织伪影(STA)。
在生物力学中,我们的肌肉、脂肪和皮肤就是果冻,而骨骼则是保龄球。STA的正式定义是皮肤标记点(或传感器)相对于下方骨骼的非刚性运动。当你踢球时,股四头肌隆起,会把你大腿上的标记点向前推。当你的脚着地时,冲击波会在你的肉体中传播,导致皮肤振动。这些都不是骨骼本身的运动,但对我们的测量系统来说,它们看起来都像是运动。刚体假设被打破了,一个“幽灵”进入了我们的机器。
这个“幽灵”仅仅是随机噪声,就像老式电视上的静电吗?它远比这更微妙和具有欺骗性。要理解这一点,我们必须首先明白,并非所有的不确定性都是生而平等的。在科学中,我们区分两种基本类型的不确定性:偶然不确定性(aleatory)和认知不确定性(epistemic)。
软组织伪影是认知和偶然成分的复杂混合体。认知成分是系统性的、可重复的:每次你踢球时,你的股四头肌都会以可预测的方式隆起,从而产生一种与运动本身相关的、一致的误差模式。偶然成分是固有的随机性:由于振动和肌肉激活的微小变化,每次迈步时皮肤的抖动可能会略有不同。这个随机部分可以通过对多次重复的运动进行平均来减少,因为波动倾向于相互抵消。然而,系统性部分带来了更大的挑战。就像相机校准中的误差(这纯粹是认知性的)一样,它会引入一种存在于每一次测量中的偏差。平均对此成分毫无帮助。
但故事还不止于此。STA不仅仅是任何随机噪声;它不是“白噪声”,白噪声具有平坦的功率谱,意味着其能量均匀分布在所有频率上。相反,STA是有色噪声。它的功率集中在特定的频率上,通常是与运动本身相关的频率。股四头肌的隆起与踢球动作同步发生。脚部撞击产生的振动发生在撞击的频率上。这意味着伪影的“指纹”与我们希望测量的运动本身交织在一起。
这种结构为我们如何识别它提供了线索。与真实白噪声(其在不同时刻之间不相关)不同,STA具有时间相关性。它的信号也倾向于在空间上是相干的——同一块皮肤上的标记点会一起摆动。这与相机传感器的随机、独立的噪声有着根本的不同。
皮肤上几毫米的微小摆动似乎微不足道。但通过严谨的数学视角,这些小误差会级联成重大且常常违反直觉的问题。
当我们将摆动的皮肤标记点视为刚性附着在骨骼上时,我们制造了一个令人不安的幻觉:骨骼本身似乎在拉伸、收缩和弯曲。考虑一个简单的任务:使用两个标记点测量大腿骨的长度。假设真实长度为,但STA在标记点之间测量的距离中引入了一个小的振荡误差,使得测量长度为。假设这个误差在一个步态周期内的平均值为零。
你可能会认为,如果误差平均为零,那么平均测量长度将是真实长度。但事实并非如此。如果我们使用一种常见且稳健的度量,如均方根(RMS)值来计算长度,一个系统性偏差就会凭空出现。测量的平方长度是。当我们对此进行时间平均时,中间项变为零,但最后一项总是正的。它的平均值,即误差的方差,大于零。结果是,平方长度的平均值大于。因此,估计的RMS长度系统性地大于真实长度。仅仅是摆动行为,即使是对称的,也使骨骼看起来更长!
这会带来严重的后果,因为许多关键的生物力学参数,如转动惯量(物体对旋转加速度的阻力),通常是根据节段长度估算的。由于转动惯量与长度的平方成正比(),这个长度上的微小正偏差在我们的动力学计算中会变成一个大得多的正偏差。
运动捕捉的主要目标通常是计算关节角度。在这里,STA同样造成了严重破坏。像大腿这样的节段的朝向是根据多个标记点的相对位置计算的。一个标记点位置的微小误差可能会被放大为计算出的角度的巨大误差。对于放置在关节旋转中心附近的标记点,这种敏感性特别高。
更糟糕的是,STA可以通过一种称为运动学串扰的现象产生虚假的运动。关节旋转通常在三个维度上描述:屈伸(前后)、外展-内收(侧向)和内-外旋(扭转)。如果STA导致我们计算出的旋转轴与真实的解剖轴略有错位,一个平面上的大的、真实的运动可能会“泄漏”或“串扰”到另一个平面。例如,一个大的髋关节屈曲运动可能被错误地解释为屈曲和一个小的、人为的内旋的组合,仅仅因为大腿上的测量轴因皮肤摆动而倾斜了。
误差并不仅限于几何学。为了理解运动的原因——肌肉产生的力和力矩——生物力学研究者会进行逆向动力学分析。这需要计算加速度,即对我们的位置数据求二阶时间导数。
微分是一种高通滤波器:它会放大高频内容。来自STA的任何微小、快速的摆动,在位置数据中可能几乎不被察觉,但在加速度数据中会变成巨大的尖峰。这对我们的力计算造成的后果是灾难性的。如在一个简化的膝关节模型中所示,计算出的关节力矩误差不仅与角度误差成正比,而且与其二阶导数成正比。对于一个以频率振荡的伪影,力矩中的误差会被乘以一个因子。这意味着对于更快的运动,我们计算出的力的误差不仅仅是增长——它是爆炸性增长。
面对这些严峻的挑战,我们该如何前进?我们如何能信任任何从皮肤上进行的测量?
最直接、尽管也最极端的方法是完全绕过皮肤。在所谓的骨针研究中,研究人员通过外科手术将无菌钢针直接插入受试者的骨骼,并将标记点附着在钢针上。这是我们的“金标准”——它提供了对骨骼真实运动的直接、明确的测量,并使我们能够同时测量邻近皮肤标记点的STA。虽然这对研究来说是无价的,但这种侵入性程序显然不适用于临床或大规模研究。
对于像惯性测量单元(IMU)这样的可穿戴传感器,问题更加困难。绑在大腿上的IMU感受到的是自身的加速度,这是骨骼加速度和伪影加速度的总和。在没有外部参考的情况下,传感器面临着一个深刻的模糊性:单个传感器从根本上无法区分真实的骨骼运动和传感器相对于骨骼的运动。
这导致了我们解决问题方法的范式转变。我们不再忽略这个幽灵,而是试图对其进行建模。先进的方法不再假设一簇标记点是刚性的。相反,它们采用复杂的优化技术,同时估计骨骼最可能的刚性运动和皮肤标记点最合理的非刚性变形。这些算法通过最小化一个平衡多个标准的成本函数来工作:
这种方法将问题从一个简单的测量转变为一个复杂的侦探故事。通过将有缺陷的证据(标记点数据)与可靠的物理原理(平滑性、近刚性)相结合,我们可以重建一个关于骨骼层面真实情况的更可信的描述。我们可以开始,小心翼翼地,驯服机器中的幽灵。
既然我们已经探讨了软组织伪影的原理和机制,我们可以退后一步,问一个关键问题:这个现象在哪些领域真正重要?它仅仅是专家们的技术烦恼,还是在更广泛的科学技术世界中具有深远的影响?正如我们将看到的,理解和应对软组织伪影不仅仅是一个边缘挑战;它是一个核心主题,推动了从手术室到牙科椅等一系列令人惊讶的学科的创新。它迫使我们更聪明,更深入地思考测量生命世界的意义,并在此过程中,揭示了测量、计算和物理定律之间美妙的相互作用。
也许最直接、最引人注目的应用在于对人体本身的研究。在临床步态分析中,我们试图诊断运动障碍。在运动生物力学中,我们旨在提高表现和预防伤害。在这两个领域,我们都依赖于精确测量骨骼的运动学——即运动。但我们的工具,通常是追踪皮肤上标记点的摄像机,并不能看到骨骼。它们看到的是皮肤,而我们知道,皮肤是一件滑溜、移动的外衣。
考虑评估前交叉韧带(ACL)功能的挑战,它是膝关节的关键稳定器。ACL的长度及其在运动过程中的变化,是其健康状况和关节稳定性的关键指标。临床医生和研究人员通过追踪股骨和胫骨来估计这个长度。一个计算模拟可以鲜明地揭示忽略STA的危险。当我们对膝关节的运动学进行建模,然后在我们的虚拟皮肤标记点上添加仅几毫米的、与运动相关的真实皮肤运动时,所导致的ACL长度估计误差可能达到几毫米——这是韧带总长度变化的很大一部分。这样的误差可能导致错误的诊断或误导性的手术计划。
这个问题是如此根本,以至于它迫使我们重新思考我们甚至如何定义关节的旋转轴。使用皮肤上的解剖标志点的方法,注定会受到STA的破坏。一个远为优雅和稳健的解决方案是功能性方法。我们不观察静止的身体,而是在运动中观察它。对于一个像铰链一样工作的关节,比如膝关节在简单摆动时,空间中必然存在一个唯一的轴,胫骨相对于股骨围绕该轴旋转。尽管每个单独的皮肤标记点都在以其自身的伪影运动抖动和滑动,但它们的集体运动必须仍然遵循骨骼的潜在刚性旋转。
通过使用一种数学程序,寻找能够最好地解释腿上所有标记点随时间变化的单一旋转轴,我们可以有效地“平均掉”混乱的皮肤运动,并找到骨骼旋转的真实、稳定的轴 [@problem__id:4197892]。这种方法是在混乱中寻找秩序的优美典范。它认识到,虽然任何单个标记点的位置都不可靠,但所有标记点运动之间的关系掌握着通往潜在真相的关键。
我们如何确定软组织伪影到底是个多大的问题?要量化一个误差,我们需要一个“金标准”——一个我们知道是真实的测量。对于人体运动,金标准是完全绕过皮肤,直接观察骨骼。这就是双平面透视(biplane fluoroscopy)的领域,这项技术本质上是一部三维X射线电影。它使我们能够看到真实的骨骼运动学,不受软组织这层面纱的干扰。
当我们比较用光学运动捕捉(OMC)——追踪皮肤标记点——和双平面透视(BF)同时测量的运动学时,结果令人警醒。在一次典型的膝关节运动分析中,从皮肤标记点估计的关节瞬时旋转中心(ICR)的位置可能偏离几厘米,而用BF追踪骨骼的误差则在毫米级别。用OMC测量的关节轴方向可能偏离超过十度,而用BF则不到一度。这种巨大的差异并非由摄像机引起;它几乎完全是由软组织伪影造成的。这一比较提供了无可辩驳的证据,证明STA通常不是一个小的二阶效应,而是运动捕捉研究中最大的单一误差来源。
通过一个柔顺、移动的层来测量一个刚性结构的斗争并非生物力学所独有。同样的基本原理和类似的解决方案在截然不同的领域中出现。
想象一位牙医使用口内扫描仪为患者创建全牙弓数字模型。扫描仪是一个小型摄像机,在牙齿和牙龈上扫过,捕捉数千张小的、重叠的3D图像,然后将它们“拼接”在一起。牙齿提供了稳定、刚性的基础。然而,牙龈和脸颊是软组织。当牙医扫描时,这些组织可能会移动和变形。如果牙医采用长的、连续的“之”字形路径,微小的配准误差和由于组织漂移引起的不匹配会累积起来,就像在运动捕捉中一样。当扫描从一侧的后磨牙到另一侧的后磨牙时,累积的误差可能达到数百微米,使得数字模型在临床上毫无用处。
解决方案是什么?“分段”扫描路径。牙医扫描一小组刚性牙齿,然后提起扫描仪,在另一组稳定的牙齿上“重新锚定”,然后再扫描下一段。这打破了误差累积的长链。每次短扫描的漂移最小,并且通过重新锚定到稳定的“骨性”标志点,全局精度得以保持。这完美地类比了用于缓解人体运动中STA的策略。
这种刚性耦合的主题在冲击生物力学这一关键领域再次出现。为了预测脑损伤的风险,我们需要测量冲击过程中头部的剧烈旋转。将传感器放在头皮上充满了困难;头皮会滑动,导致典型的STA问题。一个好得多的解决方案是使用配备仪器的护齿器。牙齿与颅骨刚性耦合。嵌入在定制护齿器中的传感器实际上与骨骼有直接、牢固的连接,绕过了有问题的软组织层。因此,其对头部加速度的测量远比皮肤传感器的测量更准确可靠。这个教训是普遍的:要测量一个刚体的运动,就要将你的传感器附着在那个刚体上。
面对这个普遍的挑战,科学家和工程师们已经发展出了一套优雅的解决方案层级体系,从简单的几何原理到复杂的、基于物理的算法。
数量上的优势:标记点簇 最简单的防御方法是不依赖于单个标记点。通过将一个由多个标记点组成的刚性簇附着到一个肢体节段上,我们可以利用平均的力量。节段的姿态是通过找到能够最好地同时对齐簇中所有标记点的旋转和平移来估计的。每个标记点下皮肤运动的独立、随机分量倾向于相互抵消,从而得到对下方骨骼朝向更稳定和准确的估计。正如人们可能直观地猜测的那样,朝向估计的不确定性,通过使用更多的标记点()和将它们散布得更远(增加它们的空间离散度)而减少。一个正式的分析证实了这种直觉,表明在第一近似下,朝向方差与这两者都成反比: 这个简单的原理是现代基于标记点的运动捕捉的基础。
智能算法:稳健估计器 标准的最小二乘拟合,即上述简单平均所做的事情,将所有数据点视为同等有效。它很容易被单个大的、虚假的运动——产生显著、非随机伪影的皮肤抽搐或肌肉收缩——所欺骗。一种更智能的方法是使用*稳健估计器*。像迭代重加权最小二乘法(IRLS)这样的算法就像一位有洞察力的法官。它首先计算一个最佳拟合姿态,然后检查每个标记点的残差。如果一个标记点是一个显著的离群值——如果它测量的位置与模型预测的位置相差很远——算法就会产生怀疑。在下一次迭代中,它会给予那个离群值更小的权重,实际上是说:“你的证词不那么可信。” 这个过程重复进行,自适应地降低离群值的影响,直到它收敛到一个反映可靠标记点的真实共识,而不是来自不可靠标记点的干扰噪声的姿态估计。
大一统:基于模型的传感器融合 最终的解决方案是认识到我们拥有的信息不仅仅是标记点的位置。我们还有其他传感器,比如测量角速度和线性加速度的惯性测量单元(IMU)。更重要的是,我们有物理定律。一个肢体节段不是一个神奇的物体;它是一个必须遵守牛顿运动定律的物理实体。一个由测力台测量的外力,必须产生相应的加速度()。
最先进的方法,通常在称为卡尔曼平滑器的框架中实现,是融合所有这些信息源。它试图找到一个单一的、动力学一致的运动“故事”,这个故事能最好地解释我们测量到的一切。它同时估计节段的位置、速度和加速度,并使其与光学标记点数据、IMU数据、测力台数据以及控制性的牛顿-欧拉运动方程最大程度地保持一致。它甚至可以共同估计并校正不同测量系统之间的时间同步误差等干扰因素。这种整体方法将问题从简单的几何拟合转变为一个约束物理优化问题,产生一个不仅噪声更小,而且在根本上更忠实于运动底层物理的运动估计。
从一个简单的麻烦到一个推动深刻科学和工程创新的驱动力,软组织伪影教会了我们宝贵的一课。它提醒我们,生命世界是复杂和柔顺的,要透过噪声看到其内在的美妙力学,需要测量、计算和物理定律的不断深化融合。