膝关节形态的特征描述方法和系统的制作方法

专利名称：膝关节形态的特征描述方法和系统的制作方法
技术领域：
本发明主要涉及膝关节形态的特征描述，每个膝关节一次或多 次，目的是为了进行研究或对具体的个体进行病理诊断。
背景技术：
膝盖解剖学的研究需要结构参数的定量特征描述。如一个具体例
子中，确定人群中的膝盖解剖学研究需要测量膝盖骨关节炎(OA)， 一种慢性发展的疾病以关节处软骨的损耗为特征，将导致关节活动损 耗并增加疼痛。在膝盖主要的两个膜(外侧和内侧)当中，OA主要
出现在内侧膜，这是由于此处承担较高的承重负荷。针对个体疾病的 纵向评价有赖于临床和射线照相的特征，主要是疼痛，伤残和结构变 化。病情改善疗法目前还在研究中且有赖于对该疾病的病情的精确评 估。
在当前实践中，用于人群研究，临床试验和膝盖OA流行病学研 究的主要前沿点是射线照相的最小关节间隙宽度(mJSW)的替代测 量，在股骨髁的承重表面或半弯曲位置的射线照片中的胫骨畔间进行 测量，如图1的射线照片所示。软骨的逐步损耗通过mJSW的收缩来 测量。由于OA主要出现在内侧膜，所以只在此膜中进行mJSW的测 量。
mJSW的测量通常是由受过训练的医师用手持式渐层滤镜察看 灯箱上的射线照片来完成。采用这种方法，难以避免由于观察者的主 观性所导致的观察者之间的显著差异。此外，由于mJSW通常在内侧 膜中进行测量，外侧膜中的OA可能没被检测到，而事实上，有些病 人患有原发性外侧膜疾病。在这种方法中mJSW是疾病病情的单一测 量，整个关节的情况可能因为依赖单一指标而被掩盖。
为了解决人类主观性所致的主要问题，利用计算机分析数字化的
膝盖射线照片来实现mJSW的测量已经被一些作者采用。用计算机分 析从X射线图像推导mJSW的方法之一由Duryea等人描述，见2000 年3月的27 Med.Phvs., dd.580國91。此外，除mJSW外的描述膝关节 间隙特征的替代结果测量已经进行了研究并与mJSW进行对比，例 如，Duryea等人在11 Osteoarthritis & Cartilage, pp.102-110 ( 2003 ) 中所述。上述方法都是基于特征分析技术的例子。
然而，需要采用一项技术对射线照相的mJSW测量进行补充， 该技术与现有的使用非基于模型的方法相比可给OA的定量测量带来 更高的灵敏度，特别是通过潜在治疗临床评价的观察可发现，凭借灵 敏度的提高，该技术可以更快地进行而且病患减少。

发明内容
根据本发明的优选实施方案，提供了一种按照模型参数化描述膝 关节特征的方法。该方法包括如下步骤
a. 创建由一组参数表征的可变形统计模板，这组参数作为一个 集合跨越一个代表膝关节的指定特征的空间位置集合的抽象向量空 间，每个向量唯一地描述其与人群平均值的偏差；
b. 在运行时在膝关节骨骼的二维图像中拟合与指定特征相关联 的多个轨迹，以允许可变形统计模板的参数化；
c. 对于一个给定对象中的多个轨迹的位置进行参数化，以得到 可变形统计模板的参数值；以及
d. 基于所述多个轨迹的子集和/或所述可变形统计模板的得到 的参数，描述给定对象的膝盖的骨骼形态的特征。
可以理解，本发明的说明书或任何所附权利要求中用到的术语 "子集，，与构成一个集合的元素有关，术语"子集，，被理解为包含整个集 合元素的真子集或全集。此外，术语"轨迹(locus)"用来表示在一个 定义某个特征的点(在向量空间中)的某一"距离，，内的多个点之一。 "距离"通常是指在向量空间上定义的一个范数。
根据本发明的其他实施方案，二维图像可以是通过使放射线穿过
膝关节所得到的射线照片，例如传送x射线透过膝关节。此方法还包
括基于根据可变形统计模板的参数对骨骼形态的特征描述来预测治
疗方法的临床结果。
根据本发明更进一步的实施方案，创建可变形统计模板的步骤包 括基于一组膝盖二维图像的统计分析定义所述一组参数，而在运行时 拟合与指定模型特征相关联的多个轨迹的步骤包括接收操作员的输
入。在连续的时间点重复步骤(c)和(d)，以便描述骨骼形态随时 间的演变。
根据本发明更进一步的实施方案， 一个附加步骤可包括在步骤 (c)的连续迭代之间进行临床干预。位置参数化的步骤可包括多个 轨迹位置的逐次逼近。描述膝盖骨骼形态特征的步骤可包括在内侧或 外侧膜中估计关节间隙在指定位置的宽度，或者至少在内侧和外侧膜 之一中估计最小关节间隙。
本发明另一方面还提供了按照骨骼形态描述膝关节特征的方法。 该方法步骤如下
a. 在运行时在膝盖骨骼的二维图像中拟合与指定特征相关联的 多个轨迹，以允许可变形统计模板的参数化；
b. 对于一个给定对象中的多个轨迹的位置进行参数化，以得到 所述可变形统计模板的参数；以及
c. 基于所述多个轨迹的子集和/或所述可变形统计模板的得到的 参数，描述所述骨骼形态的特征。
最后，才艮据本发明，可以提供一种用于计算机系统上来描述一个 对象的膝关节骨骼形态特征的计算机程序产品。该计算机程序产品包 括用于将与膝关节二维图像中的指定骨骼特征相关联的轨迹存储到 计算机存储器中的程序代码；对于一个给定对象中的多个轨迹的位置 进行参数化，以得到可变形统计模板的参数的程序代码；和基于所述 多个轨迹的子集和/或所述可变形统计模板的得到的参数来描述所述 骨骼形态的特征的程序代码。


结合附图参考下面的详细描述将更容易理解本发明，在附图中 图1A是根据本发明的实施方案的二维毛细管柱端成像设备的示
意图1B是图1A所描绘的本发明实施方案的外视图； 图2示出了根据本发明实施方案的发光二极管(LED )阵列的多 方位视图3A示出了一个典型的光强度预算表； 图3B示出了激励放射线和荧光团吸收的典型光傳范围； 图4是根据本发明实施方案的透视照明光学阵列的透视图； 图5示出了制造光学阵列的注射模型；
图6是根据本发明实施方案的一组有关针对毛细管末端设计光 学成像元件的方程式；
图7是碰撞电离CCD阵列的读出操作的示意图8示出了对荧光团浓度的灵敏度的图解；
图9是根据本发明实施方案的平行毛细管阵列电泳系统的示意
图10示出了根据本发明实施方案的緩冲器阵列的剖面图和透视
图11A示出了根据本发明实施方案的微流通道金属板排列及分 离区域连接毛细管阵列末端以分离分子排出物的顶视图11B和11C分别示出了图11A中微流通道金属板排列沿着线 B-B和A-A的剖面图。
具体实施方案
根据本发明的优选实施方案，使用已产生的可变形统计模板自动 地或半自动地实现数字化的膝盖射线照片的计算机分析，例如，该可 变形统计模板可通过大量手工注释的膝盖射线照片样例的统计分析 产生。
此处所述的和任何所附权利要求中要求保护的发明被应用于通 过使用穿透性的电磁辐射，如X射线，实现人或动物的膝关节成像所 得的数据。通常，代表使射线穿过关节的传送的二维射线照片被使用。 可以理解的是，虽然本发明不限于成像的具体形式，也不限于存储并 操作这些所得图像的具体形式，但具体形式有很多优点，如高能射线
(X射线或伽马射线)有利地提供高空间分辨率。
如下所述，膝关节图像(在此不加限制地都被称为"射线照片，，) 的分析在应用于输入数据的'模型，的基础上进行。此处和任何所附权 利要求中所采用的术语"模型"通常是指用来确定对象或其组成部分 的位置和/或动作的参数的任何数学描述。在此所述的方法应用于任何 膝盖图像数据模型都在要求保护的本发明的范围内。当 一个模型在此 被称为"统计的"时，这可被理解为基于全体对象成员中的参数变化的 分析。
更具体而言，本发明将参考一类模型进行描述，其中模型代表了 在膝盖上指定的一组2D点位置的平均相对位置，及这些相对位置在 正常情况下个体之间或具体个体随着时间推移或由于干预情况，如， 病情进展所致的变化方式的数学描述。本发明的实践建立在"看似合 理的"形态的数学模型存在的基础上，其中形态包含形状和形状变化， 并且还可以包含其他方面的外观，例如建模对象的紋理。采用下述方 法将模型应用于从实际膝盖的图像中所得的数据。然而，该方法不特 定于任何具体的点集放置，如图2所示，完全通过举例的方式，但不 限于此， 一组点(或"界标")20被自动或半自动的放置在胫骨棘、关 节的外围边界、股骨骨节的边缘和胫骨坪上。用作已知直径的校准目 标以允许在标准长度单位下表达测量结果的球轴承25也可在图2所 示的图像中找到。
在本发明实施方案的运行时应用中，指定的点通过被训练识别这 些位置的程序单元从膝关节图像中初步识别出来(在这个过程中被称 为分割)。在半自动分析中，应用的用户被请求定义膝盖射线照片上
被自动或半自动识别出的界标(图2中的大圆点26)的某个子集(真 子集或全集)，以创建可变形统计模板。这些位置按照有效"描述"膝 盖射线照片外观的方式定义，无论是否涉及极值特征。
尽管膝关节的精确形态随着对象不同和时间变化而有所变化，这 些界标大部分仍然可识别。模板是统计的，由于它对全体对象的膝盖 射线照片外观(如上所述)的分布建模(利用平均值和"正态"偏差， 或用作代表性特征的其他要素，但不限于此)。因此，该模板允许根 据有限数量的数值对形态进行参数化，本发明不限于任何具体的参数 化方案。
对于本发明的应用，点的看似合理位置的数学模型可被建立，如 下所述。 一组训练2D数据块从全体膝盖射线照片图像中取得。这些 数据在建模之前可通过手动调整来增加。
为建立一个模型，产生的是2D点的相对位置而不是它们的"绝 对，，空间参照位置。因此，在建模时，第一步通常是使每个2D数据帧 与一个公共参照系校准，可釆用多种标准校准技术之一实现，如Horn 的"Procrustes Analysis"， 《C7os^/ For附 5Wm,/ow 爿6s0/w,e OWe"/fl"V " Z7s/"g t/"i7 gwa,e簡》附》，J. Optical Societv,1987年4月 A4巻，pp.629國42。
该模型提供了一个二维点相对位置在训练数据的多帧之间的偏 差的简洁数学描述。只要数据校准，该数学描迷可通过一种或多种统 计建模技术完成，例如，Johnson和Wichern的"Principal Component Analysis",见Applied Multivariate Statistical Analysis, pp.458-513 (2002年，第5版)。
在一组实施方案中，模型可由2D数据的"平均，，形状及一组描述 形状可以如何改变的数学函数組成。通过向数学函数馈送控制数字向 量或"模型参数"，可生成任何看似合理的2D点坐标集。虽然模型参 数可能会跨越模型变化范围的空间，但是要求保护的本发明在实践时 不必如此。此外，该模型可以是线性的，也就是说，其中运动对应于 沿着直线(而不是诸如圆弧或曲线)移动的点的线性组合。然而，本
发明并不限于此类模型的使用。
现在将参照图3的流程图进行讨论，本发明的实施方案采用如上 所述被训练的模型来分割图像，即，确定描述膝盖骨骼形态特征的2D 点坐标。根据本发明的优选实施方案，产生多个假定输入点位置集， 例如随机地产生或在一组预先定义的可能位置上产生。这些位置集中 的每一个均被测试，以在模型认为其位置与相邻点的"距离""合理"的 情况下确定适当的初始输入点位置集。在本发明的范围内，此处用到 的和任何所附权利要求中的术语"距离"是指模型的参数化变量的范 数，并且可以对应于一个欧几里德范数，但并不一定如此。 一旦确定 了初始输入二维点位置，就可以找到最优校准参数集来使这些位置与 模型中的相应点的位置匹配，即，将所有坐标转换为模型的参照系， 这种方式使得模型参数描述实际图像的概率最大化。如下所迷，迭代 法被用于局部拟合。实现这一分割步骤的算法由T. F. Cootes和C. J. Taylor在《5V",/s/Zc"/ AfWe/s 。/々/ efl/"awceAf^Z/cfl/T/Mflge ^w"/戸's flm/CV 附/;w^ F油rt》Proc. SPIE Medical Imaging, (2001)中描述。
初始点集可被描述为向量X，
X= (x^X2... xn，y l5y 2…y n}，
其中(Xi,yi)是下标为i的点的二维坐标。
这些点在使用例如Horn (如前所述)被校准到模型的参照系时 被描述为向量X'，此处X'是对X应用计算出的校准变换T的结果。 X， = T(X)，
此处T是计算出的变换参数的矩阵。
模型是某种函数F，给定一组输入点坐标，这个函数可产生一个 参数向量b， b = F(X')
在一组实施方案中，模型由二维数据的"平均，，形状及一组描述形 状可以如何改变的数学函数组成，b可用下式计算 b=A(X，- X，m)
此处Xm是"平均"形状的二维点坐标的向量，A是在训练阶段学
习得到的矩阵，该训练阶段使用例如Johnson和Wichern在Applied Multivariate Statistical Analysis, pp.458-513 (2002年，第5版)中 所描述的"Principal Component Analysis"。
为了在给定一组模型参数的情况下估计一组点坐标，使用F的 "逆"F'(多数情况都只能是近似的逆)
Xe' = F'(b)
此处Xe是在给定的一组参数b下的二维坐标的估计值。如果利
用"Principal Component Analysis"来建模，则可以写成
Xe = Xm + A，(b)
此处A，是A的伪逆矩阵，在"Principal Component Analysis"的 情况下A'与A的转置矩阵是相同的。
各种模型拟合算法可被用于完成以上步骤。在一个实施方案中， 该模型值的未加权最小平方估计可如下计算
b=A(X，- X'm)
这些值可用于根据下式来预测模型参照系中的全部2D点向量的
值
Xe' = F，(b)
变换矩阵T的逆矩阵T'用于估计原始参照系中的全部2D点向
量
Xe = T'(Xe')。
随后的这组点现在与模型参照系重新校准，并且重复上述过程。 发现最佳的一组校准参数，接着应用拟合算法来获得最佳的一组模型 参数。然后，这组最佳模型参数用于产生另一組点。
这一迭代过程一直重复到出现收敛(在指定标准内)为止或者直 到指定最大次数的迭代已被完成为止。当迭代完成，最后一步的输出 就是由模型参数拟合的2D点全集的解。
最后的分割(即，射线照片中被指定的初始化点的识别)现在可 被用于从射线照片中提取多个测量结果，如下所列，仅为例举而非全 部列出。
1. 外侧或内侧膜中的mJSW。
2. 外侧或内侧膜中任意位置的JSW, JSW测量的位置可能沿着 胫骨棘到关节边缘的指定位置的 一条线被参数化，因此这种测量结果 应在纵向研究的不同时间点进行比较。
3. —个已定义的测量，具有一个区域的尺寸，描述任一膜对向的 特定区域的特征。
由于可变形统计模板已被用于检测和注释关节，关节形状的参数 化可根据此模板的最优拟合来确定。因此，在抽象意义上，关节的"形 状，，可与"普遍，，的"正常，，关节形状进行比较，或者纵向研究中时间点 上形状参数的其他变化可被当作疾病进展的新测量结果。变化可映射 到一段时间上的特定测量中，无论是否存在医疗干预。 一般而言，医 疗干预的情况在此处和任何所附权利要求中是指"治疗形式"，并且包 括药剂的监控，但不限于此。
软件应用的屏幕截图展示了一些已完成的测量结果，如图4所 示。具体地说，依照本发明，以下测量结果均是可能用于描述膝盖形 态特征的测量结果的例子
软骨面积测量
对于原始图像和后续图像并且对于任一侧膜，测量在任一侧膜端 由胫骨坪，股骨骨节和关节间隙宽度(JSWs)限定的面积。
关节间隙宽度(JSW)外形测量
对于原始图像和后续图像并且对于任一侧膜，沿着胫骨坪的整个 长度测量关节间隙宽度的外形。给定点的JSW是胫骨边界和股骨边 界的最小距离。JSW外形被当作从胫骨棘到关节中轴距离的函数来测 量(毫米)。
最小JSW测量
对于原始图像和后续图像，最小JSW可在内、外侧膜中测量。最小JSW是在给定范围内沿着胫骨碎JSW外形测量结果的最小值。 默认情况下，最小值存在于软骨面积测量结果的内部范围和JSW外
形测量结果的外部范围之间。
等价JSW测量
对于后续图像并且对于任一侧膜，JSW可在估计位置(沿着胫 骨坪)被测量，在这个位置mJSW可在原始图像上计算出来。这个测 量需要研究病人的原始图像和后续图像。
交叉标记
交叉标记是内侧膜或外侧膜的一个标记，它指示了胫骨边界是否 与股骨边界交叉。在这些情况下，如果发现交叉，那么可忽略结果。
形状测量
对于给定的一对原始图像和后续图像，系统可计算与原始膝盖和 后续膝盖之间膝盖整体形状发生变化的程度相对应的一个统计量。统 计量的值越高，原始与后续间发生的改变就越多。形状测量的单位是 毫米，数字对应于分割的膝盖边界上的每个点在原始和后续膝盖之间
(在同一参照系中测量的)发生移动的平均量。这个测量需要研究病 人的原始图像和后续图像。这应当允许"感兴趣"图像的快速鉴别，这 可以确保扩展的手工研究。
这个描述膝关节形态特征的公开方法可被实施为与计算机系统 一起使用的计算机程序产品。这样的实施可能包括一 系列计算机指 令，这些指令可固定在有形介质上，如计算机可读介质(如软盘、 CD-ROM、 ROM或硬盘)，或者可经过调制解调器或其他接口设备
(如经介质连接到网络的通信适配器)传送到计算机系统。介质可以 是有形介质(如，光学或模拟通信线路)或由无线技术实现的介质(如， 微波、红外线或其他传输技术)。这一系列计算机指令包含了上述关 于系统的全部或部分功能。本领域的技术人员应该意识到这些计算机
指令可用许多编程语言来编写，以用于多种计算机体系结构或操作系 统。此外，这些指令可被存储于任何存储设备，如半导体，磁性，光 学或其他存储设备，可用任何通信技术传送，如光学、红外线、微波 或其他传输技术。可预料，这样的计算机软件产品可作为移动介质附 加印刷文档或电子文档(例如，压缩打包软件)来发行，预先装载到
计算机系统(如，在系统ROM或硬盘上)，或者通过网络(如，因 特网或万维网)从服务器或电子公告栏来发行。当然，本发明的一些 实施方案可通过软件(如，计算机程序产品)和硬件结合的方式来实 现。本发明还有一些实施方案可完全通过硬件实现，或者完全通过软 件(如，计算机程序产品)。
本发明的上述实施方案仅仅只是示例性的，本领域的技术人员显 然会发现许多变化和改进。所有这些变化和改进都在所附权利要求限 定的本发明的范围内。
权利要求
1.一种按照模型参数化描述膝关节特征的方法，该方法包括a.创建由一组参数表征的可变形统计模板，这组参数作为一个集合跨越一个代表膝关节的指定特征的空间位置的集合的抽象向量空间，每个向量唯一地描述其与人群平均值的偏差；b.在运行时在膝关节骨骼的二维图像中拟合与指定特征相关联的多个轨迹，以允许可变形统计模板的参数化；c.对于一个给定对象中的多个轨迹的位置进行参数化，以得到所述可变形统计模板的参数值；以及d.基于所述多个轨迹的子集和/或所述可变形统计模板的得到的参数，描述所述给定对象的膝盖的骨骼形态的特征。
2. 根据权利要求1中所述的方法，其中所述二维图像是通过使 放射线穿过膝关节的传送过程所得到的射线照片。
3. 根据权利要求1中所述的方法，其中所述二维图像是通过使 X射线穿过膝关节所得到的射线照片。
4. 根据权利要求1中所述的方法，还包括基于根据可变形统 计模板的参数对骨骼形态的特征描述来预测治疗方法的临床结果。
5. 根据权利要求1中所述的方法，其中创建可变形统计模板的 步骤包括基于一组膝盖二维图像的统计分析定义所述一组参数。
6. 根据权利要求1中所述的方法，其中在运行时拟合与指定模 型特征相关联的多个轨迹的步骤包括接收操作员的输入。
7. 根据权利要求1中所述的方法，还包括在连续的时间点重复 步骤(c)和(d)，以便描述骨骼形态随着时间的演变。
8. 根据权利要求7中所述的方法，还包括在步骤(c)的连续迭 代之间进行临床千预的步骤。
9. 根据权利要求1中所述的方法，其中拟合步骤包括所述多个 轨迹的位置的逐次逼近。
10. —种按照骨骼形态描述膝关节特征的方法，该方法包括 a. 在运行时在膝盖骨骼的二维图像中拟合与指定特征相关联的 多个轨迹，以允许可变形统计模板的参数化；b. 对于一个给定对象中的多个轨迹的位置进行参数化，以得到 所述可变形统计模板的参数；以及c. 基于所述多个轨迹的子集和/或所述可变形统计模板的得到的 参数，描述所述骨骼形态的特征。
11. 根据权利要求l中所述的方法，其中描述骨骼形态特征的步 骤包括估计在内侧或外侧膜中的指定位置的关节间隙宽度。
12. 根据权利要求l中所述的方法，其中描述骨骼形态特征的步 骤包括估计在内侧和外侧膜的至少一个中的最小关节间隙。
13. —种用于计算机系统上来描述一个对象的膝关节骨骼形态特 征的计算机程序产品，该程序产品包括a. 用于将与膝关节二维图像中的指定骨骼特征相关联的轨迹存 储到计算机存储器中的程序代码；b. 对于一个给定对象中的多个轨迹的位置进行参数化，以得到 可变形统计模板的参数的程序代码；和c. 基于所述多个轨迹的子集和/或所述可变形统计模板的得到的 参数来描述所述骨骼形态的特征的程序代码。
全文摘要
本发明涉及一种描述膝关节骨骼形态特征的方法和系统。与骨骼结构模型相关的多个轨迹都是按实物尺寸定做的，且被用于确定多元轨迹在给定对象中的位置参数，从而得到可变形统计模板的参数。接着基于所得可变形统计模板的参数来描述骨骼形态特征。
文档编号G06T3/00GK101189638SQ200680009437
公开日2008年5月28日 申请日期2006年3月16日 优先权日2005年3月24日
发明者凯文·沃尔克, 安东尼·霍姆斯, 简·哈斯拉姆 申请人:欧博达西亚医学有限公司