JP5054669B2 - 膝関節形態の特徴付けのための方法及びシステム - Google Patents

Description

本発明は、一般に、研究目的又は特定の個人の病理診断のためにそれぞれの膝関節に対して１回又は複数回行われる、膝関節の特徴付けに関する。

膝の解剖学的研究は、構造パラメータの定量的特徴付けを必要とする。ある特定の例として、識別された母集団の膝の解剖学的研究は、関節の軟骨の損失により特徴付けられ関節運動の損失及び増大する痛みをもたらす、ゆっくりとした進行性疾患である膝の変形性関節症（ＯＡ）の測定を必要とする。膝の主要な２つのコンパートメント（外側及び内側）のうち、内側コンパートメントにかかる体重を支える負荷の方がより大きいため、ＯＡは主に内側コンパートメントに見られる。個体におけるこの疾患の長期的評価は、臨床的特徴及び放射線写真の特徴、主に痛み、障害及び構造上の変化に依存している。疾患緩和治療法が現在開発中であるが、これらは疾患の進行の正確且つ精密な判定に依存することになる。

現在の診療では、膝のＯＡに関する、母集団を対象とした研究、臨床試験、及び疫学調査における主要な評価項目は、図１の放射線写真に示されているように半屈曲位で撮影された放射線写真から大腿骨顆部及び脛骨プラトーの体重を支える両方の面の間で測定される、放射線写真による最小関節裂隙（ｍｉｎｉｍｕｍ Ｊｏｉｎｔ Ｓｐａｃｅ Ｗｉｄｔｈ（ｍＪＳＷ））の代理測定値である。軟骨の進行性消失は、ｍＪＳＷが狭くなっていくことによって測定される。ＯＡは主に内側コンパートメントにおいて見られるため、ｍＪＳＷは通常、このコンパートメントでのみ測定される。

ｍＪＳＷの測定は通常、訓練を受けた医師によって、ライトボックスの上の放射線写真を精査しながら目盛り付きの手持ち式レンズを使って行われる。この方法を使った場合、観察する人間の主観のため、観察者間及び観察者内の有意な変動を避けることが困難である。さらに、ｍＪＳＷは通常内側コンパートメントでしか測定されないため、外側コンパートメントのＯＡが検出を免れる可能性があり、実際、原発性外側コンパートメント疾患を患う患者もいる。最後に、ｍＪＳＷが病気の進行に関する唯一の測定値であるという点で、唯一の指標に依存することよって関節全体への影響が隠される可能性がある。

人間の主観という重大な問題に対処するために、ｍＪＳＷの測定のためにデジタル化された膝の放射線写真のコンピュータ解析が幾人かの著者によって用いられている。Ｘ線像からｍＪＳＷを導出するコンピュータ解析の１つの方法が、Ｄｕｒｙｅａ他によりＭｅｄ．Ｐｈｙｓ．２７，５８０−９１ページ（２０００年３月）に記載されている。さらに、ｍＪＳＷ以外の膝関節空間の特徴付けに対する代理的結果測定が、例えば、Ｄｕｒｙｅａ他のＯｓｔｅｏａｒｔｈｒｉｔｉｓ ＆ Ｃａｒｔｉｌａｇｅ１１，１０２−１１０ページ（２００３年）で研究され、ｍＪＳＷと比較されている。上述の方法は、特徴に基づいた解析技術の例である。

しかしながら、非モデルベースのアプローチを用いて得られる感度よりもさらに高感度なＯＡの定量的測定に寄与する技術で、特に可能性のある治療法の臨床評価を視野に入れると、感度が高いので、より少ない患者で、より迅速に実行できる技術で、放射線写真のｍＪＳＷ測定を補うことが望ましい。

本発明の好ましい実施形態により、モデルのパラメータ化の観点から膝関節を特徴付ける方法が提供される。その方法は、
ａ．膝関節の特定の特徴の空間的位置の集合を表す抽象ベクトル空間を集合として張るパラメータの集合であって、各ベクトルが母平均からの変分を一意的に記述する、パラメータの集合によって特徴付けられる、変形可能な統計的テンプレートを創出するステップと、
ｂ．実行時に、膝の骨格の２次元画像中に特定の特徴に関連する複数の軌跡をフィッティングして、変形可能な統計的テンプレートに関してのパラメータ化を可能にする、ステップと、
ｃ．所与の被検体における複数の軌跡の位置をパラメータ化して、変形可能な統計的テンプレートのパラメータに対する値を導出する、ステップと、
ｄ．複数の軌跡の部分集合、及び／又は変形可能な統計的テンプレートの導出されたパラメータに基づいて、所与の被検体の膝の骨格形態を特徴付けるステップと
を有する。

本発明の説明又は付随のいずれの請求項においても、ある集合を構成する元に関連して「部分集合」という用語が用いられる場合、「部分集合」という用語は、元の集合全体の真部分集合又は完全な部分集合を包含すると理解されるべきである。さらに、「地点（ｌｏｃｕｓ）」という語は、特徴を定義する（ベクトル空間内の）点からある「距離」の内にある複数の点の１つを意味するために用いられる。「距離」は、一般に、ベクトル空間上で定義されるノルムを意味する。

本発明の他の実施形態によれば、２次元画像は、膝関節を貫通する放射線の透過によって、例えば膝関節を通過するＸ線の透過によって得られる放射線写真とすることができる。この方法は、変形可能な統計的テンプレートのパラメータに関しての骨格形態の特徴付けに基づいて、治療様式の臨床転帰を予測するステップをさらに含むことができる。

本発明のさらなる実施形態によれば、変形可能な統計的テンプレートを創出するステップは、膝の２次元画像の集合の統計的解析に基づいてパラメータの集合を定義することを含むことができ、実行時に特定のモデルの特徴に関連する複数の軌跡をフィッティングするステップは、操作者の入力を受け取ることを含むことができる。前述の方法の（ｃ）及び（ｄ）のステップは、経時的な骨格形態の進展を記述するために、継続的な時点において繰り返されてもよい。

本発明のさらに他の実施形態によれば、ステップ（ｃ）の継続的な反復の間に、臨床的介入を行う追加のステップを含むことができる。位置をパラメータ化するステップは、複数の軌跡の位置の逐次近似を含むことができる。膝の骨格形態を特徴付けるステップは、内側コンパートメント又は外側コンパートメントの特定の位置での関節裂隙を推定すること、或いは、内側コンパートメント及び外側コンパートメントの少なくとも１つにおいて最小関節裂隙を推定することを含むことができる。

本発明のさらなる態様は、骨格形態に関して膝関節を特徴付けるための方法を提供する。この方法は、
ａ．実行時に、膝の骨格の２次元画像中に特定の特徴に関連する複数の軌跡をフィッティングして、変形可能な統計的テンプレートのパラメータ化を可能にする、ステップと、
ｂ．所与の被検体における複数の軌跡の位置をパラメータ化して、変形可能な統計的テンプレートのパラメータを導出する、ステップと、
ｃ．複数の軌跡の部分集合、及び／又は変形可能な統計的テンプレートの導出されたパラメータに基づいて、骨格形態を特徴付けるステップと
を有する。

最後に、本発明によれば、被検体の膝関節の骨格形態を特徴付けるためにコンピュータ・システム上で使用するためのコンピュータ・プログラム製品が提供される。このコンピュータ・プログラム製品は、膝関節の２次元画像中の特定の骨格特徴に関する複数の軌跡をコンピュータ・メモリに格納するプログラム・コードと、所与の被検体における複数の軌跡の位置をパラメータ化して、変形可能な統計的テンプレートのパラメータを導出するためのプログラム・コードと、複数の軌跡の部分集合、及び／又は変形可能な統計的テンプレートの導出されたパラメータに基づいて、骨格形態を特徴付けるためのプログラム・コードと、を有する。

本発明は、添付の図面と共に以下の説明を参照することによって、より容易に理解されるであろう。

本発明の好ましい実施形態によれば、デジタル化された膝の放射線写真のコンピュータ解析が、例えば、手書きの注釈のある膝の放射線写真の多数の実例の統計的解析によって生成された変形可能な統計的テンプレートを用いて、自動的に又は半自動的に行われる。

本明細書に記載され、添付の任意の請求項において特許請求されている発明は、例えばＸ線のような貫通性の電磁放射線の使用による、ヒト又は動物の膝関節の画像化によって得られるデータに適用される。典型的には、膝関節を貫通する放射線の透過を示す２次元放射線写真が用いられる。本発明は、特定の画像化方式又は得られた画像を格納し操作する特定の方式に範囲を限定されるものではないことが理解されるべきであるが、高エネルギー（Ｘ線又はガンマ線）放射線によって有利に提供される高空間分解能のような特定の方式から生じる利点がある。

膝関節の画像（本明細書において、限定することなく「放射線写真」と称する）の解析は、入力データに適用される「モデル」に基づいて、以下に説明されるように進められる。本明細書及び付随のいずれの請求項において用いられる場合にも、「モデル」という用語は、一般に、被検体又はその構成部分の位置及び／又は運動のパラメータ化を規定するいかなる数学的記述をも指す。本明細書で説明されている方法を膝の画像データへのいかなるモデルへと適用することも、特許請求される本発明の範囲内である。モデルが本明細書において「統計的」と称されている場合、被検体の母集団の要素間のパラメータの変動の解析に基づくと理解されるべきである。

さらに詳細には、本発明は、あるクラスのモデルに関して説明され、そのモデルは、膝の２Ｄ点の位置の特定の集合の平均相対位置を、これらの相対位置が、正常な環境において個人間で変動する様子、又は時の経過、若しくは、例えば疾患の進行のような介入する環境に起因して特定の個人において変化する様子の数学的記述と共に表す。本発明の実施は、「妥当と思われる（ｐｌａｕｓｉｂｌｅ）」形態の数学的モデルの存在を前提としており、ここで、形態（ｍｏｒｐｈｏｌｏｇｙ）は、形状及び形状の変動を包含し、モデル化された物体の質感のような外観の他の様相をも包含することができる。以下に説明される方法は、モデルを実際の膝の画像から得られたデータに適用するために用いられる。しかしながら、この方法は、点集合のどのような特別な配置にも特定されておらず、図２において、単に例示としてであり、限定するものとしてではなく、頚骨棘（ｔｉｂｉａｌ ｓｐｉｎｅ）、関節の周縁境界部、大腿骨顆部の縁及び脛骨プラトーの上に自動的に又は半自動的に配置された点の集合（又は「ランドマーク」）２０として図示されている。測定値を長さの標準単位で表すことを可能にするための既知の直径の較正目標として用いられるボールベアリング２５も、図２に示される画像中に示されている。

本発明の実施形態の実行時の適用の際には、この特定の位置は、これらの位置を識別するように訓練されたプログラム要素によって、膝関節の画像において事前に識別される（本明細書でセグメント化と呼ばれる工程）。半自動的な分析では、このアプリケーションのユーザは、変形可能な統計的テンプレートの創出において、自動的又は半自動的に識別された膝の放射線写真上のランドマークの部分集合(真部分集合又は完全部分集合）（図２の大き目の点２６）を定義するよう求められる。これらの位置は、膝の放射線写真の外観を、極値的特徴に関連して、或いは他の方法で、効果的に「記述」できるような方式で定義される。

膝関節の正確な形態は、被検体間で異なり、且つ経時的に変化するが、これらのランドマークはほとんどの部分で識別可能なままである。このテンプレートは、被検体の集団全体にわたる（上述のような）膝の放射線写真の外観の（平均の「正規分布的（ｎｏｍａｌ）」変動、又は限定ではないが、代表的な特徴として役立つ他のモーメントを有する）分布をモデル化しているという点で統計的である。このように、このテンプレートは、有限個の値に関して形態のパラメータ化を可能にし、本発明は、いずれの特定のパラメータ化のスキームとも無関係である。

ここで説明されるように、本発明における適用のために、点の妥当と思われる位置の数学的モデルを構築することができる。訓練用（ｔｒａｉｎｉｎｇ）２Ｄデータブロックの集合が、膝の放射線画像の集団から取り出される。これらのデータは、モデル構築のプロセスの前に、手動の調整により増強されてもよい。

モデルを構築するためには、空間を参照した２Ｄ点の「絶対」位置よりも、それらの相対位置の方が重大である。それゆえ、モデルの構築における最初のステップは、典型的には、２Ｄデータの各座標系（ｆｒａｍｅ）を共通の基準座標系に位置合わせすることであり、種々の位置合わせ技術のうちの１つ、例えば、Ｈｏｒｎの、Ｊ．Ｏｐｔｉｃａｌ Ｓｏｃｉｅｔｙ Ａ４巻，６２９〜４２ページ（１９８７年４月）のＣｌｏｓｅｄ Ｆｏｒｍ Ｓｏｌｕｔｉｏｎ ｏｆ Ａｂｓｏｌｕｔｅ Ｏｒｉｅｎｔａｔｉｏｎ Ｕｓｉｎｇ Ｕｎｉｔ Ｑｕａｎｔｅｒｎｉｏｎｓに記載されている「プロクラステス（Ｐｒｏｃｒｕｓｔｅｓ）分析」などを用いて達成することができる。

モデルは、訓練用データの座標系間での２Ｄ点の相対位置の変動について、コンパクトな数学的記述を与える。データが位置合わせされると、これは、例えば、Ｊｏｈｎｓｏｎ及びＷｉｃｈｅｎにより、Ａｐｐｌｉｅｄ Ｍｕｌｔｉｖａｒｉａｔｅ Ｓｔａｔｉｓｔｉｃａｌ Ａｎａｌｙｓｉｓ(第５版、２００２年)、４５８〜５１３ページに記載されているような「主成分分析（Ｐｒｉｎｃｉｐａｌ Ｃｏｍｐｏｎｅｎｔ Ａｎａｌｙｓｉｓ）」をはじめとする数種の統計的モデル化技術の１つ又はそれ以上によって行われることができる。

１つの実施形態の組では、モデルは、２Ｄデータについての「平均」形状と共に、その形状がどのように変化し得るかを記述する１組の数学的関数から構成されることができる。制御数（ｃｏｎｔｒｏｌ ｎｕｍｂｅｒ）又は「モデル・パラメータ」のベクトルをこの数学的関数に入力することにより、どのような妥当と思われる２Ｄ点の座標の集合でも生成することができる。モデル・パラメータは、モデル形状の可動域（ｅｘｃｕｒｓｉｏｎ）の空間を張ることができるが、それは特許請求される本発明の実施には必ずしも当てはまらない。さらに、運動が（例えば、弧又は曲線よりも）直線に沿って動く点の線形結合に対応するという意味で、モデルは線形となることもある。しかしながら、本発明はそのようなモデルへの適用可能性に限定されるものではない。

ここで、図３のフロー図を参照して論じられるように、本発明の実施形態は、一旦、前述のように訓練されると、そのモデルを用いて画像をセグメント化する。すなわち、膝の骨格形態を特徴付ける２Ｄ点の座標の位置を確認する。本発明の好ましい実施形態によれば、例えば、無作為に、又は予め定義された可能な位置の集合にわたって、多数の推定入力点位置の集合が生成される。適切な初期入力点位置の集合を識別するために、モデルが「合理的である」と認識する近隣の点からの「距離」内での位置決めに基づいて、これらの集合の各々が検定される。本明細書及び付随の任意の請求項において用いられる「距離」という用語は、モデルのパラメータ化された変数に関するノルムを指し、ユークリッド・ノルムに対応する場合もあるが、本発明の範囲において必須ではない。一旦、初期入力２Ｄ点位置が識別されたら、これらの位置をモデル内の対応する点の位置に適合させるための、すなわち、モデルの基準座標系に対する全ての座標をモデル・パラメータが実際の画像を記述する確率を最大化するような方式で変換するための、最良の位置合わせパラメータの集合が見出される。その後、局所化されたフィッティングのために、以下に述べられるように、反復が用いられる。そのようなセグメント化ステップを実行するためのアルゴリズムは、Ｔ．Ｆ．Ｃｏｏｔｅｓ及びＣ．Ｊ．Ｔａｙｌｏｒの、Ｐｒｏｃ．ＳＰＩＥ Ｍｅｄｉｃａｌ Ｉｍａｇｉｎｇ（２００１年）のＳｔａｔｉｓｔｉｃａｌ Ｍｏｄｅｌｓ ｏｆ Ａｐｐｅａｒａｎｃｅ ｆｏｒ Ｍｅｄｉｃａｌ Ｉｍａｇｅ Ａｎａｌｙｓｉｓ ａｎｄ Ｃｏｍｐｕｔｅｒ Ｖｉｓｉｏｎに記載されている。

初期の点の集合は、ベクトルＸとして表すことができ、
Ｘ＝｛Ｘ₁，Ｘ₂．．．Ｘ_n，ｙ₁，ｙ₂．．．ｙ_n｝
であり、ここで（Ｘ_i，ｙ_i）は、添え字ｉの点の２Ｄ座標である。

点は、例えばＨｏｒｎ（上記）の方法を用いてモデルの基準座標系に位置合わせされた場合、ベクトルＸ’と表され、ここでＸ’は算出された位置合わせ変換ＴをＸに適用した結果であり、
Ｘ’＝Ｔ（Ｘ）
であり、ここでＴは算出された変換パラメータの行列である。

モデルはある関数Ｆであり、これは、入力点座標の集合が与えられたならば、パラメータのベクトルｂを生成する。
ｂ＝Ｆ（Ｘ’）

１つの実施形態の組において、モデルは、２Ｄデータについての「平均」形状と共に、その形状がどのように変化し得るかを記述する１組の数学的関数から構成され、ｂは以下の式を用いて計算され、
ｂ＝Ａ（Ｘ’− Ｘ’_m）
ここで、Ｘ_mは「平均的」形状についての２Ｄ点座標のベクトルであり、Ａは、例えば、Ｊｏｈｎｓｏｎ及びＷｉｃｈｅｎによって、Ａｐｐｌｉｅｄ Ｍｕｌｔｉｖａｒｉａｔｅ Ｓｔａｔｉｓｔｉｃａｌ Ａｎａｌｙｓｉｓ(第５版、２００２年)、４５８〜５１３ページに記載された「Ｐｒｉｎｃｉｐａｌ Ｃｏｍｐｏｎｅｎｔ Ａｎａｌｙｓｉｓ」の中で説明されているような主成分分析を用いて、訓練段階中に学習された行列である。

モデルパラメータの集合が与えられたときに、点座標の集合を推定するために、Ｆの「逆関数」（ほとんどの場合、単に近似的逆関数でよい）であるＦ’が用いられる。
Ｘ_e’＝Ｆ’（ｂ）
ここで、Ｘ_eは、所与のパラメータの集合ｂについての２Ｄ座標の推定値である。モデルが主成分分析を用いて構築されている場合、これは以下のように表される。
Ｘ_e＝Ｘ_m＋Ａ’（ｂ）
ここで、Ａ’はＡの擬逆関数であり、これは、主成分分析の場合、Ａの転置に等しい。

前述のステップを達成するために、種々のモデル・フィッティング・アルゴリズム用いることができる。一実施形態において、モデル値の単純な重みづけされていない最小二乗推定値が計算される。
ｂ＝Ａ（Ｘ’− Ｘ_m’）
これらの値は、次式
Ｘ_e’＝Ｆ’（ｂ）
に従って、モデルの基準座標系における２Ｄ点ベクトル全体の値を予測するために用いられる。

変換行列Ｔの逆関数Ｔ’は、元の基準座標系の２Ｄ点ベクトル全体を推定するために用いられる。
Ｘ_e＝Ｔ’（ Ｘ_e’）

この、次の点の集合をここでモデルの基準座標系で再度位置決めし、説明してきたプロセスが繰り返される。最良の位置決めパラメータの集合が見出され、次に、フィッティング・アルゴリズムが適用されて、最良のモデル・パラメータの集合を導出する。次に、この最良のモデル・パラメータの集合を用いて、別の点の集合を生成する。

この反復プロセスは（指定された基準内に）収束するまで、そうでなければ指定された最大数の反復が遂行されるまで、繰り返される。反復が終了したとき、最終ステップの出力が、このモデル・パラメータによってフィッティングされた２Ｄ点の完全な集合に対する解である。

最終的なセグメント化（すなわち、放射線写真における特定の初期化点の識別）をここで用いて、以下に概要を述べるように、放射線写真から多数の測定値を抽出することができるが、これらは例として提示されているだけであり、包括的な集合として提示されているものではない。
１．外側コンパートメント又は内側コンパートメントのいずれかのｍＪＳＷ。
２．外側コンパートメント又は内側コンパートメントの任意の位置でのＪＳＷ。ＪＳＷ測定の位置は、頚骨棘から関節の縁部の特定の位置までの線に沿ってパラメータ化されることができ、それゆえ、この測定値は長期的研究における種々の時点で比較されることができる。
３．面積の次元を持ち、いずれかのコンパートメントで境界を定められた特定領域を特徴づける、定義された測定値。

変形可能な統計的テンプレートを用いて関節を検出して注釈をつけているため、関節形状のパラメータ化は、このテンプレートの最適化されたフィッティングから決定されることができる。したがって、抽象的な意味での関節の「形状」は「正常な」関節の「領域（ｕｎｉｖｅｒｓｅ）」の形状と比較されてもよく、あるいは長期的研究の時点における形状パラメータの変化を疾患の進行の新規な測定値として用いることができる。ある期間の経過に亘って、医療的介入の存在下、又はそれ以外で、特定の測定値における変化をマッピングすることができる。医療的介入の場合は、最も一般的な意味で、本明細書及び付随する任意の請求項において「治療様式」として言及され、医薬の投与を含むが、それに限定されるものではない。

実行されているこれらの測定のいくつかを示すソフトウェア・アプリケーションの画面例が図４に示される。特に、以下の測定が、本発明によって膝の形態の特徴付けのために用いることができる測定の例である。

軟骨面積の測定
基準画像及び追跡画像の両方に対して、且つ、どちらのコンパートメントに対しても、脛骨プラトー、大腿骨顆部及びコンパートメントの両端の関節裂隙（ＪＳＷ）によって定められる面積が測定される。

関節裂隙（ＪＳＷ）プロファイルの測定
基準画像及び追跡画像の両方に対して、且つ、どちらのコンパートメントに対しても、脛骨プラトーの全長に沿って、関節裂隙（ＪＳＷ）の輪郭を測定することができる。所与の点でのＪＳＷは、脛骨の境界と大腿骨の境界との間の最小距離である。ＪＳＷプロファイルは、関節の中央軸（ｍｅｄｉａｌ ａｘｉｓ）に沿った脛骨棘からの距離の関数として（ミリメートル（ｍｍ）で）測定される。

最小ＪＳＷの測定
基準画像及び追跡画像の両方に対して、最小ＪＳＷが、内側コンパートメント及び外側のコンパートメントの両方で測定される。最小ＪＳＷは、脛骨プラトーに沿った所与の範囲内のＪＳＷプロファイル測定値の最小値である。デフォルトでは、最小値は、軟骨面積測定値の内部範囲（ｉｎｎｅｒ ｅｘｔｅｎｔ）とＪＳＷ輪郭測定値の外部範囲（ｏｕｔｅｒ ｅｘｔｅｎｔ）との間で求められる。

等価ＪＳＷの測定
追跡画像に対して、且つ、どちらのコンパートメントに対しても、ＪＳＷは、基準画像においてｍＪＳＷが算定された（脛骨プラトーに沿った）推定位置で測定することができる。この測定は、ひとりの患者についての基準画像及び追跡画像の両方が研究に付されていることを必要とする。

クロスオーバー・フラグ
クロスオーバー・フラグは、内側コンパートメント又は外側コンパートメントのいずれかで、脛骨の境界が大腿骨の境界と重なっているかどうかを示すフラグである。クロスオーバーが検出された場合、その結果は無視することが望ましい。

形状の測定
所与の一対の基本画像及び追跡画像に対して、システムは、膝の全体的な形状が基本画像の膝と追跡画像の膝との間でどれくらい変化したかに対応する統計量を計算することができる。この統計量の値が高ければ高いほど、基本と追跡との間でより大きな変化があることになる。形状測定の単位はミリメートル（ｍｍ）であり、数字は、同じ基準座標系で測定された、セグメント化された膝の境界上の各点が基本画像の膝と追跡画像の膝との間で移動した量の平均値に相当する。この測定は、ひとりの患者についての基本画像及び追跡画像の両方が研究に付されていることを必要とする。これは、拡張された手作業の調査の正当な理由となり得る「興味深い」画像の迅速な識別を可能にするはずである。

膝関節の形態を特徴付けるための開示された方法は、コンピュータ・システムと共に使用するためのコンピュータ・プログラム製品として実施されることができる。そのような実施は、コンピュータ可読媒体（例えば、ディスケット、ＣＤ−ＲＯＭ、ＲＯＭ、又は固定ディスク）のような有形媒体上に固定されるか、又は媒体を超えたネットワークに接続された通信アダプタのようなモデム又は他のインターフェース機器を介してコンピュータ・システムに伝達可能な、一連のコンピュータ命令を含むことができる。この媒体は、有形媒体（例えば、光学又はアナログの通信線）又は無線技術で実施される媒体（例えば、マイクロ波、赤外線、又は他の伝送技術）とすることができる。一連のコンピュータ命令は、システムに関して本明細書で前述された機能の全て又は一部を具現化する。当業者は、そのようなコンピュータ命令は、多くのコンピュータ・アーキテクチャ又はオペレーティング・システムとともに使用するための多くのプログラミング言語で記述することができるということを理解すべきである。さらに、そのような命令は、半導体、磁気、光学、又は他の記憶装置といった任意の記憶装置に格納されることができ、光学、赤外線、マイクロ波、又は他の伝送技術といった任意の通信技術を用いて伝送されることができる。そのようなコンピュータ・プログラム製品は、添付の印刷された文書又は電子文書と共に取り外し可能媒体として配布され（例えば、シュリンク包装されたソフトウエア）、コンピュータ・システムに予め組み込まれ（例えば、システムＲＯＭは固定ディスク上に）、或いはネットワーク（例えば、インターネット又はワールド・ワイド・ウェブ）上のサーバ又は電子掲示板から配布されることができることが予期される。もちろん、本発明の幾つかの実施形態は、ソフトウエア（例えば、コンピュータ・プログラム製品）とハードウエアとの組み合わせとして実施されることができる。さらに他の本発明の実施形態は、完全にハードウエアとして、又は完全にソフトウエア（例えば、コンピュータ・プログラム製品）として実施される。

説明された本発明の実施形態は、単に例示として意図され、当業者には多くの変更及び修正が意図される。そのような全ての変更及び改変は、付随の特許請求の範囲で定義される本発明の範囲内であると意図される。

従来の標準的な、膝の変形性関節症の代理測定値としての、内側コンパートメントのｍＪＳＷ測定法である。 本発明の適用のために規定した、典型的な解剖学的ランドマーク群を示す。 本発明の好ましい態様の、膝の画像データへの適用を示すフローチャートである。 本発明の好ましい態様による、膝関節の形態の具体的な測定のソフトウェア分析を示す。

Claims (14)

1. 膝の骨格形態の変化を測定する方法であって、
   ａ．１組の入力点座標が与えられたときにパラメータのベクトルを生成する関数であるモデルを、トレーニング組の膝画像から生成するステップと、
   ｂ．実行時に、膝の骨格の２次元画像中に特定の特徴に関連する複数の地点をフィッティングして、変形可能な統計的テンプレートのパラメータ化を可能にすることにより、適切な組の初期入力点位置を識別するステップと、
   ｃ．最良の組のモデルパラメータを反復して計算して、前記フィッティングした地点に関する前記モデルのノルムを最小化することにより、所与の被検体における前記複数の地点の位置をパラメータ化して、前記変形可能な統計的テンプレートの前記パラメータを導出する、ステップと、
   ｄ．前記変形可能な統計的テンプレートの前記導出された最良の組のモデルパラメータに基づいて前記膝の骨格形態を測定するステップであって、（i)前記導出された最良の組のモデルパラメータに基づいて１組の点を生成する段階、及び、（ii）該１組の点に基づいて、関節裂隙の測定値を、頚骨棘から前記関節の縁部における特定の位置までの線に沿ってパラメータ化された位置において、抽出する段階、により、前記膝の骨格形態を測定するステップと、
   ｅ．第１の時点及び第２の時点において前記ステップｂから前記ステップｄ（ii）までを実行するステップと、
   ｆ．前記第１の時点と前記第２の時点との間における前記位置における前記関節裂隙の測定値の変化を決定するステップと、
   を含むことを特徴とする方法。
2. 前記２次元画像が前記膝関節を貫通する放射線の透過によって得られる放射線写真である、請求項１に記載の方法。
3. 前記２次元画像が前記膝関節を通過するＸ線の透過によって得られる放射線写真である、請求項１に記載の方法。
4. 前記モデルが、前記トレーニング組の膝画像の統計的解析に基づいて生成される、請求項１に記載の方法。
5. 実行時に特定のモデルの特徴に関連する複数の地点をフィッティングする前記ステップが、操作者の入力を受け取ることを含む、請求項１に記載の方法。
6. 前記フィッティングするステップが、前記複数の地点の位置の逐次近似を含む、請求項１に記載の方法。
7. 膝の骨格形態の変化を測定する方法であって、
   ａ．実行時に、膝の骨格の２次元画像中に特定の特徴に関連する複数の地点をフィッティングして、変形可能な統計的テンプレートのパラメータ化を可能にすることにより、適切な組の初期入力点位置を識別するステップと、
   ｂ．最良の組のモデルパラメータを反復して計算して、前記フィッティングした地点に関するモデルのノルムを最小化することにより、所与の被検体における前記複数の地点の位置をパラメータ化して、前記変形可能な統計的テンプレートの前記パラメータを導出する、ステップと、
   ｃ．前記変形可能な統計的テンプレートの前記導出された最良の組のモデルパラメータに基づいて、前記膝の骨格形態を測定するステップであって、（i)前記導出された最良の組のモデルパラメータに基づいて１組の点を生成する段階、及び、（ii）該１組の点に基づいて、関節裂隙の測定値を、頚骨棘から前記関節の縁部における特定の位置までの線に沿ってパラメータ化された位置において、抽出する段階、により、前記膝の骨格形態を測定するステップと、と
   ｄ．第１の時点及び第２の時点において前記ステップａから前記ステップｃ（ii）までを実行するステップと、
   ｅ．前記第１の時点と前記第２の時点との間における前記位置における前記関節裂隙の測定値の変化を決定するステップと、
   を含むことを特徴とする方法。
8. 前記骨格形態を測定する前記ステップが、内側コンパートメント及び外側コンパートメントのうちの少なくとも一方において関節裂隙を推定することを含む、請求項１に記載の方法。
9. 前記骨格形態を測定する前記ステップが、内側コンパートメント及び外側コンパートメントのうちの少なくとも一方において最小関節裂隙を推定することを含む、請求項１に記載の方法。
10. 被検体の膝関節の骨格形態の変化を測定するためにコンピュータ・システム上で使用するためのコンピュータ・プログラムであって、
    ａ．膝関節の２次元画像中の特定の骨格特徴に関する複数の地点をコンピュータ・メモリに格納するプログラム・コードと、
    ｂ．最良の組のモデルパラメータを反復して計算して、前記地点に関するモデルのノルムを最小化することにより、所与の被検体における前記複数の地点の位置をパラメータ化して、変形可能な統計的テンプレートのパラメータを導出するためのプログラム・コードと、
    ｃ．前記変形可能な統計的テンプレートの前記導出された最良の組のモデルパラメータに基づいて、前記膝の前記骨格形態を測定するためのプログラム・コードであって、（i）前記導出された最良の組のモデルパラメータに基づいて１組の点を生成する段階、及び、（ii）該１組の点に基づいて、関節裂隙の測定値を、頚骨棘から前記関節の縁部における特定の位置までの線に沿ってパラメータ化された位置において、抽出する段階、により、前記膝の骨格形態を測定するプログラム・コードと、
    ｄ．前記プログラム・コードａから前記プログラム・コードｃまでを第１の時点及び第２の時点において実行させるプログラム・コードと、
    ｅ．前記第１の時点と前記第２の時点との間における前記位置における前記関節裂隙の測定値の変化を決定するプログラム・コードと、
    を含むことを特徴とするコンピュータ・プログラム。
11. 前記骨格形態を測定するプログラム・コードが、内側コンパートメント及び外側コンパートメントのうちの少なくとも一方において関節裂隙を推定するプログラム・コードを含む、請求項１０に記載のコンピュータ・プログラム。
12. 前記骨格形態を測定するプログラム・コードが、内側コンパートメント及び外側コンパートメントのうちの少なくとも一方において最小関節裂隙を推定するプログラム・コードを含む、請求項１０に記載のコンピュータ・プログラム。
13. 前記骨格形態を測定するステップが、内側コンパートメント及び外側コンパートメントのうちの少なくとも一方において関節裂隙を推定することを含む、請求項７に記載の方法。
14. 前記骨格形態を測定するステップが、内側コンパートメント及び外側コンパートメントのうちの少なくとも一方において最小関節裂隙を推定することを含む、請求項７に記載の方法。

Images