JP2016513529A - マーカーまたはストラットの投影を用いて物体の位置を決定する方法 - Google Patents

マーカーまたはストラットの投影を用いて物体の位置を決定する方法 Download PDF

Info

Publication number
JP2016513529A
JP2016513529A JP2015562961A JP2015562961A JP2016513529A JP 2016513529 A JP2016513529 A JP 2016513529A JP 2015562961 A JP2015562961 A JP 2015562961A JP 2015562961 A JP2015562961 A JP 2015562961A JP 2016513529 A JP2016513529 A JP 2016513529A
Authority
JP
Japan
Prior art keywords
markers
ray
coordinates
struts
ray source
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2015562961A
Other languages
English (en)
Other versions
JP6209624B2 (ja
Inventor
ミハイエフ、コンスタンチン・エフゲニー
ヴヴィデンスキー、ペトル・スタニスラヴォヴィチ
モレンコ、アレクサンドル・アレクサンドロヴィチ
シヴァチェフ、ドミトリー・アレクサンドロヴィチ
サムチェコフ、ミハイル
チェルカシン、アレクサンドル
Original Assignee
アメイ テクノロジーズ インコーポレイテッド
アメイ テクノロジーズ インコーポレイテッド
テキサス スコティッシュ ライト ホスピタル フォー チルドレン
テキサス スコティッシュ ライト ホスピタル フォー チルドレン
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by アメイ テクノロジーズ インコーポレイテッド, アメイ テクノロジーズ インコーポレイテッド, テキサス スコティッシュ ライト ホスピタル フォー チルドレン, テキサス スコティッシュ ライト ホスピタル フォー チルドレン filed Critical アメイ テクノロジーズ インコーポレイテッド
Publication of JP2016513529A publication Critical patent/JP2016513529A/ja
Application granted granted Critical
Publication of JP6209624B2 publication Critical patent/JP6209624B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • G06T19/20Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • G06T7/85Stereo camera calibration
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B17/00Surgical instruments, devices or methods, e.g. tourniquets
    • A61B17/56Surgical instruments or methods for treatment of bones or joints; Devices specially adapted therefor
    • A61B17/58Surgical instruments or methods for treatment of bones or joints; Devices specially adapted therefor for osteosynthesis, e.g. bone plates, screws, setting implements or the like
    • A61B17/60Surgical instruments or methods for treatment of bones or joints; Devices specially adapted therefor for osteosynthesis, e.g. bone plates, screws, setting implements or the like for external osteosynthesis, e.g. distractors, contractors
    • A61B17/62Ring frames, i.e. devices extending around the bones to be positioned
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10116X-ray image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/20Indexing scheme for editing of 3D models
    • G06T2219/2004Aligning objects, relative positioning of parts

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Architecture (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computer Graphics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Surgery (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biomedical Technology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

マーカーを備える物体と、固定部材と、ストラットとを用いて、身体部分の3Dモデルが作成される。本方法は、X線源とX線撮像装置との間に配置された身体部分と物体との、第1および第2のレントゲン画像に関連する。本方法は、その後、マーカーまたはストラットの投影の間の距離の第1および第2の組を決定し、マーカー間またはストラット間の事前に決定された距離と、それらの投影の間の距離の第1および第2の組とを用いて、X線撮像装置に対する、X線源および物体の第1および第2の3D位置を決定する。本方法は、その後、2つの向きにおけるX線撮像装置に対するマーカーまたはストラットの3D位置を用いて、3D参照フレームにおいて、第1および第2の3D物体投影を位置合わせする。本方法は、3D物体投影に基づいて、物体の3Dモデルを作成する。【選択図】図2A

Description

[技術分野]
本開示は、概してX線撮像法の分野に関し、より詳細には、空間における物体の3次元(3D)モデルを2次元(2D)レントゲン画像に基づいて作成することに関する。
[背景]
3D空間における物体のモデル化には、いくつかの有用な用途がある。物体の3Dモデルにより、各物体の相互の向きをより容易に可視化し分析することができる。モデル化のこの側面は、整形外科学において(または、とくに骨変形(bone deformities)を分析する際に)、とくに有用である。コンピュータ断層撮影法(CT)は、人体組織の3D表現を生成するために、整形外科分野において用いられてきた従来技法の一つである。別の従来技術は、2Dのレントゲン画像の支援による骨変形の可視化および分析を含む。まず、変形した骨セグメントの放射線像が直交視点(orthogonal views)において取得される。続いて、撮像された骨セグメントの2D線形表現を作成し、そのような線形表現を変形平面に投影することにより、変形を分析することができる。代替として、2Dのレントゲン画像における変形した骨セグメントのアウトラインを手作業で決定し、これを外挿して、変形した骨セグメントの3Dモデルを構築してもよい。
[概要]
本開示は、身体部分の3Dモデルを作成する方法であって、前記身体部分は物体に結合され、前記物体は、前記物体に沿って、事前に決定された距離において、複数のマーカーを備える、方法を提供する。一実施形態では、前記方法は、
1)X線源とX線撮像装置との間に配置された前記身体部分および前記物体の、第1のレントゲン画像を受信することであって、前記第1のレントゲン画像は、
‐前記身体部分と、
‐前記物体と、
‐前記複数のマーカーと
の画像を含む、第1のレントゲン画像を受信することと、
2)前記X線源と前記X線撮像装置との間に配置された前記身体部分および前記物体の、第2のレントゲン画像を受信することであって、前記第2のレントゲン画像は、
‐前記身体部分と、
‐前記物体と、
‐前記複数のマーカーと
の画像を含む、第2のレントゲン画像を受信することと、
3)前記第1のレントゲン画像上の前記複数のマーカーの投影の間の距離の第1の組を決定することと、
4)前記複数のマーカー間の前記事前に決定された距離と、前記第1のレントゲン画像上の前記複数のマーカーの投影の間の距離の前記第1の組とを用いて、前記X線撮像装置に対する、前記X線源の第1の3D位置と、前記物体の第1の3D位置とを決定することと、
4)前記第2のレントゲン画像上の前記複数のマーカーの投影の間の距離の第2の組を決定することと、
5)前記複数のマーカー間の前記事前に決定された距離と、前記第2のレントゲン画像上の前記複数のマーカーの投影の間の距離の前記第2の組とを用いて、前記X線撮像装置に対する、前記X線源の第2の3D位置と、前記物体の第2の3D位置とを決定することと、
6)前記第1および第2の向きにおける前記X線撮像装置に対する、前記複数のマーカーの前記各3D位置を用いて、3D参照フレームにおける前記第1および第2の3D物体投影を位置合わせすることと、
7)前記第1および第2の3D物体投影に基づいて、前記3D参照フレームにおける撮像された物体の3Dモデルを作成することと
を備える。一実施形態では、前記複数のマーカーは、複数のストラットが少なくとも1つのリングに接続される複数のジョイントを含む。一実施形態では、前記物体は整形外科固定器である。一実施形態では、前記複数のマーカーは5つのマーカーを含み、前記X線源の3D位置と物体の3D位置とは数学的関係によって決定される。他の実施形態では、前記複数のマーカーは4つのマーカーを含み、前記X線源の3D位置と物体の3D位置とは数学的関係によって決定される。別の実施形態は、物体の3Dモデルを作成する方法であって、
前記物体は物体に結合され、
前記物体は、事前に決定された長さの複数のストラットを備え、
前記複数のストラットは、それぞれ、事前に決定された諸寸法の少なくとも2つの固定部材に接続され、
前記方法は、
1)前記X線源とX線撮像装置との間に配置された前記物体の、第1のレントゲン画像を受信することであって、前記第1のレントゲン画像は、
‐前記物体と、
‐事前に決定された長さの前記複数のストラットであって、前記複数のストラットは、それぞれ少なくとも2つの固定部材に、2つの接続点で接続され、前記2つの接続点の間の距離は事前に決められている、前記複数のストラットと
の画像を含む、第1のレントゲン画像を受信することと、
2)前記X線源と前記X線撮像装置との間に配置された前記物体の、第2のレントゲン画像を受信することであって、前記第2のレントゲン画像は、
‐前記物体と、
‐事前に決定された長さの前記複数のストラットであって、前記複数のストラットは、それぞれ少なくとも2つの固定部材に、2つの接続点で接続され、前記2つの接続点の間の距離は事前に決められている、前記複数のストラットと
の画像を含む、第2のレントゲン画像を受信することと、
3)前記第1のレントゲン画像上の、前記複数のストラットの長手軸の投影の第1の組を決定することと、
4)前記複数のストラットの前記接続点間の前記事前に決定された距離と、前記第1のレントゲン画像上の前記複数のストラットの長手軸の投影の前記第1の組とを用いて、前記X線撮像装置に対する、前記X線源の第1の3D位置と、前記物体の第1の3D位置とを決定することと、
5)記第2のレントゲン画像上の、前記複数のストラットの長手軸の投影の第2の組を決定することと、
6)前記複数のストラットの前記接続点間の前記事前に決定された距離と、前記第2のレントゲン画像上の前記複数のストラットの長手軸の投影の前記第2の組とを用いて、前記X線撮像装置に対する、前記X線源の第2の3D位置と、前記物体の第2の3D位置とを決定することと
7)前記第1および第2の向きにおける前記X線撮像装置に対する、前記複数のストラットの前記各3D位置を用いて、3D参照フレームにおける前記第1および第2の3D物体投影を位置合わせすることと、
8)前記第1および第2の3D物体投影に基づいて、前記3D参照フレームにおける撮像された物体の3Dモデルを作成することと
を備える、方法を提供する。一実施形態では、前記物体は整形外科固定器である。一実施形態では、前記複数のストラットは5つのストラットを含み、前記X線源の3D位置と物体の3D位置とは数学的関係によって決定される。他の実施形態では、前記複数のストラットは4つのストラットを含み、前記X線源の3D位置と物体の3D位置とは数学的関係によって決定される。
一実施形態では、方法は、
前記第1のレントゲン画像において、撮像された身体部分の、第1身体部分アウトラインを特定することと、
前記第2のレントゲン画像において、撮像された身体部分の、第2身体部分アウトラインを特定することと、
前記第1身体部分アウトラインから、前記X線源の前記第1の3D位置への、第1の3D身体部分投影を準備することと、
前記第2身体部分アウトラインから、前記X線源の前記第2の3D位置への、第2の3D身体部分投影を準備することと、
前記第1の身体部分投影および前記第2の身体部分投影に基づき、前記3D参照フレームにおいて前記撮像された身体部分の3Dモデルを作成することと
を備える。別の実施形態では、方法は、
前記3D参照フレームにおいて傾斜軸を特定することであって、前記傾斜軸は、
‐前記第1の向きにおける前記X線源の前記第1の位置に対応する、前記3D参照フレームにおける第1の3D位置と、
‐前記第2の向きにおける前記X線源の前記第2の位置に対応する、前記3D参照フレームにおける第2の3D位置と、
の間を通過する、傾斜軸を特定することと、
前記傾斜軸を通過し、前記3D参照フレームにおける前記撮像された身体部分の前記第1および第2の3D投影を通過する、1つ以上の交差平面を特定することと、
前記1つ以上の交差平面のそれぞれについて、次のステップa)〜c)を実行すること:
a)前記3D参照フレームにおいて、前記第1の3D身体部分投影および前記第2の3D身体部分投影の間の交点を1つ以上と、前記交差平面とを特定すること、
b)前記交差平面内で各前記交点を接続する1つ以上の多角形を準備すること、
c)前記1つ以上の多角形のそれぞれの内部に1つ以上の閉曲線を準備することであって、前記1つ以上の閉曲線は、前記交差平面における前記撮像された身体部分の断面図に対応する、1つ以上の閉曲線を準備すること、
と、
前記撮像された身体部分の3Dモデルを形成するために、前記3D参照フレームにおいて、前記閉曲線のそれぞれを接続する表面を準備することと
をさらに備える。
本開示およびその利点のより完全な理解のために、添付図面と併せて以下の記載への参照が提供される。添付図面において、類似の参照番号は類似の特徴を示す。
整形外科固定器の一実施形態の斜視図である。 整形外科固定器の別の一実施形態の斜視図である。 本開示による、撮像システムと、第1の向きおよび第2の向きにおける物体との斜視図である。 本開示による、撮像システムと、第1の向きおよび第2の向きにおける物体との、別の斜視図である。 本開示による、撮像システムと、ある向きにおけるマーカーを伴う物体との斜視図である。 本開示による、X線源と撮像装置へのマーカーの投影との座標の概略図である。 本開示による、撮像システムと、ある向きにおけるストラットを伴う物体との斜視図である。 本開示による、X線源と、撮像装置へのストラット接続点の投影との座標の概略図である。 本開示による例示的アルゴリズムを図示するフロー図である。 本開示による、3Dフレームワークにおいて撮像された物体を図示する概略図である。 物体の例示的な3Dモデルを図示する概略図である。 物体の3Dモデルを作成するための例示的システムを図示する概略ブロック図である。 物体の3Dモデルを決定し、外部固定器を調整するための例示的セットアップを図示する概略図である。 本開示による、物体の第1および第2のレントゲン画像を図示する概略図である。 本開示による、光源の向きと、図12Aに示すそれぞれ対応するレントゲン画像とを図示する概略図である。 本開示による、傾斜軸をそれぞれ通過する複数の交差平面と、図12Aに示す第1および第2のレントゲン画像とを図示する概略図である。 本開示による、図12Aに示す第1のレントゲン画像からの3D物体投影と、第2のレントゲン画像からの3D物体投影の間の複数の交差部分を図示する概略図である。 本開示による、図12Dに示す交差平面内の交差部分を接続する複数の多角形を図示する概略図である。 本開示による、図12Eに示す1つ以上の多角形のそれぞれに含まれる複数の1つ以上の閉曲線を図示する概略図である。 本開示による、図12Fに示す各閉曲線を接続する表面を図示する概略図である。 本開示による、物体に取り付けられる複数の物体マーカーを図示する概略図である。 本開示による、様々な光源位置に対する第1および第2のレントゲン画像からの複数の基準(fiducial)の投影を図示する概略図である。 本開示による、マーカーの投影の3Dモデルを図示する概略図である。 本開示による、図13Cに示すマーカーの投影の3Dモデルを図示する、別の概略図である。 本開示による、3D空間内のマーカーの別のモデルを図示する概略図である。
[詳細な説明]
3Dモデルを生成するための従来技術は、多数の欠点を有する。CTスキャンは、人体組織の3D表現を作成するために組み合わせることができる断面図の組を生成する。しかしながら、CTスキャンを整形外科用途に用いることは、いくつかの制限により現実的ではない可能性がある。CTスキャンの間、患者は比較的大量の放射線を受け、CTスキャンを繰り返し用いることは患者を過度の放射線に曝露して健康リスクを与える可能性がある。さらに、CTスキャンは比較的高価であり、望ましくない変形を起こし得る金属を撮像するには適切ではない。また、患者はCTスキャンの間静止することを要求され、患者が幼い子供である場合には麻酔が必要となる可能性がある。しかしながら、麻酔の使用は処置のコストを増加させ、追加の健康リスクを与える可能性がある。
別の従来技術は、変形した骨セグメントのアウトラインを2Dレントゲン画像において手作業で決定することと、2Dのアウトラインを外挿して変形した骨セグメントの3Dモデルを構築することとを含む。しかしながら、そのような技術を用いて作成されたモデルの精度には、様々な要素が悪影響を与え得る。第一に、変形した骨セグメントの線形表現を投影することは、3D空間における骨セグメントの周囲(girth)の詳細を明らかにはせず、内科医(physician)に、骨変形を十分に修正しない処置を指示させるかもしれない。さらに、従来技術によって作成されたモデルは、各レントゲン画像が直交位置で取得されたという想定に基づいており、そうでない場合にはモデルの精度に悪影響を生じる。レントゲン画像を取得するための直交位置を推定するよう技師を訓練することはできるが、小さいヒューマンエラーは不可避であり、したがって従来技術によって生成されたモデルを不正確なものにする。さらに、X線源から撮像装置まで進むX線の拡大効果により、レントゲン画像内の物体は実際のサイズより大きく現れる。この拡大効果に対応するために、既知の寸法の参照マーカー(複数可)が、関心領域に近接する物体上に正確に配置されなければならない(この参照マーカーの既知の寸法が、拡大効果を決定し対応するために用いられる)。ここでも、参照マーカーの配置における、人間による不可避の不正確さが、精度の低下をもたらす可能性がある。
従来技術における上述の誤差に起因して、取得される線形パラメータおよび角度パラメータは、真のパラメータではなく投影となる。投影は物体の真のサイズや形状とは対応せず、物体の真の形状と比べると歪んでいる。そのような技術は、選択された物体の3D空間における各点の座標を正確に決定するためには適切ではなく、また、整形外科用途においては、そのような方法は、組織セグメントの望ましい伸延、圧縮、変位または他の移動を正確に計算するには適切ではない。
本開示は、レントゲン画像を用いて物体の3Dモデルを作成するための技術を提供する。X線の拡大効果に対応するよう正確に配置されたマーカーを用いる必要性を、本開示の技術が回避する可能性があるということを、当業者は本開示から理解するであろう。また、本開示の技術は、直交位置で取得されたレントゲン画像を要しない可能性があり、様々な相対的向きにおいて取得されたレントゲン画像に適している可能性がある。さらに、本開示の技術は、レントゲン画像を取得する時に撮像装置に配置されたマーカーの使用を必要としない可能性がある。そして、本開示の技術は、基準を伴うマーカーを用いる必要性を回避する可能性がある。
本開示の各実施形態により、2Dレントゲン画像に基づく物体の正確な3Dモデル化が可能となる。これらの実施形態は、既知の幾何学的配置を持つ物体を用いて、身体各部(骨等)の位置を決定する可能性がある。物体の一実施形態は、外部固定器(固定部材、ストラットおよび/またはマーカーを含む)であってもよい。図1Aは、2つの固定部材101と、3つのストラット102と、6つのマーカー103とを含む、外部固定器100(環状固定器)の実施形態である。6つのマーカー103は、2つの固定部材101が3つのストラット102に接続する点である。他の実施形態は、外部固定器100の異なる部分に配置されたマーカーを有してもよい。また、外部固定器100は、固定部材101に取り付けられた追加のマーカー104を用いてもよい。外部固定器100は、身体部分105を取り囲み、身体部分105の破損(fracture)が治癒可能となるように身体部分105を固定するために使用されてもよい。固定部材101は、本実施形態ではリングであるが、他の実施形態は、六角形状、矩形状、五角形状または他の適切な形状に成形された固定部材101を有してもよい。さらに、本実施形態は3つのストラット102を用いるが、他の実施形態は任意の適切な数(4つ、5つまたは6つ等)のストラットを用いてもよい。
図1Bは、外部固定器110の別実施形態である。これは、2つの固定部材101と、6つのストラット102と、12個のマーカー113とを含み、12個のマーカーは、2つの固定部材101が6つのストラット102に接続する位置に配置されている。ここでも、他の実施形態は、固定部材101の異なる部分に配置された追加のマーカー104を有してもよい。図1A〜図1Bに図示される各実施形態は、単に例示的なものであり、本明細書に開示される様々な設計要素または既知の技術に従って変更してもよいということが理解される。
本開示により、空間内で、X線源、身体部分、および物体の位置を導出するために、物体(外部固定器等)の既知の幾何学的配置および2Dレントゲン画像へのその投影を含む数学的モデルを用いて、身体部分を3Dモデル化することが可能となる。空間内で身体部分の位置を決定することにより、内科医または他の医療スタッフメンバーは、処置中の骨の最適な固定のために物体(外部固定器等)を調整してもよい。また、彼らは他の医療目的の身体部分の既知の位置を用いてもよい。
いくつかの実施形態は、物体の3Dモデル化を作成するために複数のマーカーを用いる数学的モデルを採用する。一実施形態では、複数のマーカーはさらに基準を含んでもよい。しかしながら、この複数のマーカーは、ストラットが固定部材に接する位置のマーカーまたは基準を伴うマーカーに限らない。さらに、複数のマーカーは、5つのマーカーまたは4つのマーカーを含んでもよい。
他の実施形態は、物体の3Dモデル化を作成するために、複数のマーカーに代えて複数のストラットを用いる数学的モデルを採用してもよい。一実施形態では、複数のストラットは、固定部材に接続されたストラットを含んでもよく、複数のストラットが5つのストラットまたは4つのストラットを含んでもよい。
[マーカーの投影の使用]
本明細書において開示される技術の一実施形態は、X線源と撮像装置との間に配置された身体部分と物体との、第1および第2のレントゲン画像を受信することを含む。身体部分は物体に結合している(coupled)。図2Aは、第1および第2の向き220,230において第1および第2のレントゲン画像をそれぞれ取得するために動作可能な概略図である。第1および第2のレントゲン画像を取得するために、X線源202と撮像装置204との間に、物体210に取り囲まれた身体部分201が配置される。物体210は整形外科固定器であってもよく、より具体的には図2Aに示すようなヘキサポッド(hexapod)であってもよい。物体210は、物体210に沿って事前に決定された距離に複数のマーカー203を備える。第1のレントゲン画像を生成するために、物体210に取り囲まれた身体部分201と、X線源202と、撮像装置204とは、互いに対して第1の向き220にある。第2のレントゲン画像は、物体210に取り囲まれた身体部分201を、X線源202および撮像装置204に対して新たな第2の向き230まで回転させることによって、または、図2Aに示すように、X線源202および撮像装置204を、身体部分201の周りに新たな第2の向き230まで回転させることによって、生成されてもよい。したがって、第1のレントゲン画像は、第1の向き220における、身体部分201と、物体210と、複数のマーカー203との画像を含む。同様に、第2のレントゲン画像は、第2の向き230における、身体部分201と、物体210と、複数のマーカー203との画像を含む。本実施形態では、マーカー203はストラット205が固定部材206に接続する位置にあるが、他の実施形態では、マーカー203は物体の他の部分に配置されてもよい。固定部材206はリングまたは他の適切な形状として構成することができる。
本実施形態は、X線撮像装置204上に描画される、第1のレントゲン画像の複数のマーカー203の投影の第1の組と、X線撮像装置204上に描画される、第2のレントゲン画像の複数のマーカー203の投影の第2の組とを決定する。図2Bは、第1のレントゲン画像220に描画され、また、第2のレントゲン画像230に描画される、2つのマーカーの投影の決定を図示する。本実施形態は、その後、複数のマーカー間の事前に決定された距離と、第1のレントゲン画像220上に描画される複数のマーカーの投影の第1の組とを用いて、X線撮像装置に対する、X線源202の第1の3D位置と、物体210の第1の3D位置とを決定してもよい。同様に、本技術は、複数のマーカー間の事前に決定された距離と、第2のレントゲン画像230上の複数のマーカーの投影の第2の組とを用いて、X線撮像装置に対する、X線源202の第2の3D位置と、物体210の第2の3D位置とを決定してもよい。一実施形態では、複数のマーカーの投影の第1および第2の組は、より詳細には、それぞれ、第1および第2のレントゲン画像上の複数のマーカー間の距離の投影であってもよい。X線源202および物体の、第1および第2の3D位置を決定するために、様々な数学的モデルを採用することができる。これらのモデルは、後により詳しく記述する。
X線源および物体の、第1および第2の3D位置が決定されると、その後、本技術は、第1および第2の向きにおけるX線撮像装置に対し、複数のマーカーの各3D位置を用いて、3D参照フレームにおいて第1および第2の3D物体投影を位置合わせ(align)させてもよい。その後、本実施形態は、第1および第2の3D物体投影に基づいて、撮像された身体部分の3Dモデルを3D参照フレームにおいて作成してもよい。別の実施形態は、物体210の3Dモデルを、単独で、または、撮像された身体部分の3Dモデルに加えて、作成してもよい。
[モデル1]
本明細書に開示される技術の一実施形態では、複数のマーカーは物体に関連付けられた5つのマーカーを含んでもよい。図3は、ある向きの、2つの固定部材302と、少なくとも5つのマーカー303とを備える物体310(この場合には整形外科固定器)を示す。これら少なくとも5つのマーカー303は、ストラット306が固定部材302に接する位置にある。他の実施形態では、これら少なくとも5つのマーカー303は、物体310の異なる部分に配置されてもよい。本実施形態では、X線源301の第1の3D位置と、物体310の第1の3D位置とは、以下に論じる数学的モデルを用いて決定されてもよい。この数学的モデルでは、(x,y,z)はX線源301の座標であり、(x0,y0,z0)〜(x4,y4,z4)は5つのマーカー303の座標であり、(X0,Y0,Z0)〜(X4,Y4,Z4)は第1のレントゲン画像上の5つのマーカー303の投影の第1の組の座標であり、l01,l02,l03,l04,l12,l13,l14,l23,l24,l34は5つのマーカー303間の事前に決定された距離であると想定される。このように18個の未知のパラメータがあり、X線源301および物体310の位置を決定するためには18個の方程式の系が必要である。
3点(x,y,z),(xi,yi,zi),(Xi,Yi,Zi)は、X線源301からマーカー303へ、そして撮像装置304へと通過する同一の直線上に位置している。したがって、この直線の方程式は、
Figure 2016513529
と表せる。ここでiは0から4までの数である。
この方程式は、代替的に、2つの方程式により表現することもできる:
Figure 2016513529
本実施形態では5つのマーカーがあり、5つのマーカーそれぞれについて上記の方程式の対があるので、マーカー間の距離を記述する方程式が10個あることになる。マーカーの位置に関与する追加の8つの方程式を、ピタゴラスの定理を用いて導出してもよい。したがって、結果として次の方程式が得られてもよい:
Figure 2016513529
しかしながら、この方程式の組は、各マーカー相互の相対的な位置を記述する方程式をすべて含むわけではない。したがって、系に含まれない次の方程式に関連して、X線源301および物体310の3D位置に対する方程式の組の解をチェックする必要がある:
Figure 2016513529
X線源301の第2の3D位置と、物体310の第2の3D位置とは、以下に論じる、実質的に同様な数学的モデルを用いて決定されてもよい。この数学的モデルでは、(‘x,‘y,‘z)はX線源301の座標であり、(´X0,´Y0,´Z0)〜(´X4,´Y4,´Z4)は第2のレントゲン画像上の5つのマーカーの投影の第2の組の座標であると想定される。
Figure 2016513529
しかしながら、この方程式の組は、各マーカー相互の相対的な位置を記述する方程式をすべて含むわけではない。したがって、系に含まれない次の方程式に関連して、X線源301および物体310の3D位置に対する方程式の組の解をチェックする必要がある:
Figure 2016513529
これによって、X線源301および物体310の、第1および第2の3D位置を求めてもよい。
特定の状況下では、X線源301のZ座標の変化は、レントゲン画像上の各マーカーの投影の比例的変化(proportional change)を生じる場合がある。これらの状況下では、X線源301のZ座標は、方程式の組を解く際に解く必要のない定数パラメータとして設定することができる。ここで、Z座標に対する定数パラメータは、物体がX線源301とレントゲン画像との間にフィットできるようにする大きな数とすべきである。これにより、以下のモデル2〜4に例示するように、X線源301の3D位置と物体310の3D位置とを決定するために、より少ない数のマーカーを用いることができるようになる。
[モデル2]
本技術の本実施形態では、複数のマーカーは4つのマーカーを含んでもよい。本実施形態では、X線源301の第1の3D位置と、物体310の第1の3D位置とは、以下に論じる数学的モデルを用いて決定されてもよい。この数学的モデルでは、(x,y,z)はX線源301の座標であり、(x0,y0,z0)〜(x3,y3,z3)は4つのマーカー303の座標であり、(X0,Y0,Z0)〜(X3,Y3,Z3)は第1のレントゲン画像上の4つのマーカー303の投影の第1の組の座標であり、l01,l02,l03,l12,l13,l23は4つのマーカー303間の事前に決定された距離であると想定される。これらの変数の間の関係は、さらに図4に示される。
3点(x,y,z),(xi,yi,zi),(Xi,Yi,Zi)は、X線源301からマーカー303へ、そして撮像装置304へと通過する同一の直線上に位置している。したがって、この直線の方程式は、
Figure 2016513529
と表せる。
この方程式は、代替的に、2つの方程式により表現することもできる:
Figure 2016513529
ここでiは0から3までの数である。
本実施形態では4つのマーカーがあり、4つのマーカーそれぞれについて上記の方程式の対があるので、マーカー間の距離を記述する方程式が8つあることになる。マーカーの位置に関与する追加の6つの方程式を、ピタゴラスの定理を用いて導出してもよい。したがって、結果として次の方程式が得られてもよい:
Figure 2016513529
X線源301の第2の3D位置と、物体310の第2の3D位置とは、以下に論じる、実質的に同様な数学的モデルを用いて決定されてもよい。この数学的モデルでは、(‘x,’y,’z)はX線源301の座標であり、(´X0,´Y0,´Z0)〜(´X3,´Y3,´Z3)は第2のレントゲン画像上の4つのマーカーの投影の第2の組の座標であると想定される。
Figure 2016513529
これによって、X線源301および物体310の、第1および第2の3D位置を求めてもよい。
[モデル3]
本技術の別の実施形態でもまた、複数のマーカーは4つのマーカーであってもよい。この実施形態では、X線源301の第1の3D位置と、物体310の第1の3D位置とは、以下に論じる数学的モデルを用いて決定されてもよい。この数学的モデルでは、(x,y,z)はX線源301の座標であり、(x0,y0,z0)〜(x3,y3,z3)は4つのマーカー303の座標であり、(X0,Y0,Z0)〜(X3,Y3,Z3)は第1のレントゲン画像上の4つのマーカー303の投影の第1の組の座標であり、l01,l02,l03,l12,l13,l23は4つのマーカー303間の事前に決定された距離であると想定される。これらの変数の間の関係は、さらに図4に示される。
本実施形態は、次のパラメータ方程式を採用する:
x=x0+α*t; y=y0+β*t; z=z0+γ*t
ここで、α,β,γは方向ベクトルを示し、tはある直線上の点(x,y,z)を別の点(たとえば(x0,y0,z0))に対して特徴づけるパラメータである。3点(x,y,z),(xi,yi,zi),(Xi,Yi,Zi)は同一の直線上に位置している。したがって、この直線のパラメータ方程式は、次のように表せる:
xi=x+α*t; yi=y+β*t; zi=z+γ*t
ただしα=Xi−x; β=Yi−y; γ=Zi−z
ここでiは0から3までの数である。
本実施形態は、物体上のマーカー間の距離についての適切な数の方程式(既知の幾何学的配置の接続点間の6つの距離について6つの方程式)を追加することにより、14個の方程式を採用する。これらをマーカー間の距離についての方程式に組み込むことにより、本実施形態は次の6つの方程式を提供する:
Figure 2016513529
本実施形態は、t0〜t3,x,yを、上記6つの方程式を解くことにより決定してもよい。その後、X線源301の第1の3D位置と物体310の第1の3D位置とを決定してもよい。X線源301の第2の3D位置と物体310の第2の3D位置とは、実質的に同様な数学的モデルを用いて決定されてもよい。
3点(‘x,‘y,‘z),(xi,yi,zi),(‘Xi,‘Yi,‘Zi)は同一の直線上に位置している。したがって、この直線のパラメータ方程式は、次のように表せる:
xi=´x+α*t; yi=´y+β*t; zi=´z+γ*t
ただしα=´Xi−x; β=´Yi−y; γ=´Zi−z
ここでiは0から3までの数である。
本実施形態は、物体上のマーカー間の距離についての適切な数の方程式(既知の幾何学的配置の接続点間の6つの距離について6つの方程式)を追加することにより、14個の方程式を採用する。これらをマーカー間の距離についての方程式に組み込むことにより、本実施形態は次の6つの方程式を提供する:
Figure 2016513529
これにより、X線源301および物体(t0〜t3,x、y)の、第1および第2の3D位置を求めてもよい。
[モデル4]
本技術の別の実施形態では、複数のマーカーは4つのマーカーを含んでもよい。この実施形態では、X線源301の第1の3D位置と、物体310の第1の3D位置とは、以下に論じる数学的モデルを用いて決定されてもよい。この数学的モデルでは、(x,y,z)はX線源301の座標であり、(x0,y0,z0)〜(x3,y3,z3)は4つのマーカー303の座標であり、(X0,Y0,Z0)〜(X3,Y3,Z3)は第1のレントゲン画像上の4つのマーカー303の投影の第1の組の座標であり、l01,l02,l03,l12,l13,l23は4つのマーカー303間の事前に決定された距離であると想定される。これらの変数の間の関係は、図4に示される。さらに、各マーカーは、X線源301と、そのマーカーの投影点と、隣接する(neighboring)マーカーの投影点とを通る平面上に配置される。このように、物体のマーカーそれぞれについて
2平面付属方程式(two planes appurtenant equations)が決定され、以下に提供するように、4つのマーカーを伴う物体310について8つの方程式の系を作成してもよい。
Figure 2016513529
ここで、本実施形態は、マーカー間の距離についての適切な数の方程式(たとえば、既知の幾何学的配置固定器の4つの点間の6つの距離について6つの方程式)を追加することにより、次のような14個の方程式の数学的モデルを提供する:
Figure 2016513529
X線源301の第1の3D位置と、物体310の第1の3D位置とは、上記の方程式を解くことにより決定されてもよい。X線源301の第2の3D位置と、物体310の第2の3D位置とは、以下に論じる、実質的に同様な数学的モデルを用いて決定されてもよい。この数学的モデルでは、(‘x,’y,’z)はX線源301の座標であり、(´X0,´Y0,´Z0)〜(´X3,´Y3,´Z3)は第2のレントゲン画像上の4つのマーカーの投影の第2の組の座標であると想定される。
Figure 2016513529
これにより、X線源301および物体の、第1および第2の3D位置を決定してもよい。
[ストラットの投影の使用]
本明細書に開示される技術の、別の実施形態は、X線源と撮像装置との間に配置された物体の、第1および第2のレントゲン画像を受信することを含む。図2Aは、第1の向き220および第2の向き230において、それぞれ第1のレントゲン画像および第2のレントゲン画像を取得するために動作可能な概略図である。第1および第2のレントゲン画像を取得するために、物体210に取り囲まれた身体部分201が、X線源202と撮像装置204との間に配置される。物体210は整形外科固定器であってもよく、より詳細には図2Aに示すヘキサポッドであってもよい。物体210は、事前に決められた寸法の少なくとも2つの固定部材206にそれぞれ接続された、事前に決められた長さの複数のストラット205を備える。この実施形態では固定部材206はリングであるが、他の実施形態は他の形状に成形された固定部材を採用してもよい。第1のレントゲン画像を生成するために、物体210によって取り囲まれた身体部分201と、X線源202と、撮像装置204とが、互いに対して第1の向き220にある。第2のレントゲン画像は、物体210に取り囲まれた身体部分201を、X線源202および撮像装置204に対して新たな第2の向き230まで回転させることによって、または、図2Aに示すように、X線源202および撮像装置204を、身体部分201の周りに新たな第2の向き230まで回転させることによって、生成されてもよい。したがって、第1のレントゲン画像は、身体部分201と、物体210と、事前に決められた寸法の少なくとも2つの固定部材206にそれぞれ接続された事前に決められた長さの複数のストラット205との第1の画像を含む。同様に、第2のレントゲン画像は、身体部分201と、物体210と、事前に決められた寸法の少なくとも2つの固定部材206にそれぞれ接続された事前に決められた長さの複数のストラット205との第2の画像を含む。
本実施形態は、第1のレントゲン画像上に描画される複数のストラット205の投影の第1の組と、第2のレントゲン画像上に描画される複数のストラット205の投影の第2の組とを決定する。一実施形態では、複数のストラット205の投影の第1および第2の組は、より詳細には、複数のストラット205の長手軸の投影を含んでもよい。別の実施形態では、複数のストラット205の投影の第1および第2の組は、ストラット接続点の投影を含む。ストラット接続点とは、ストラット205が固定部材206に接する位置をいう。別の実施形態では、ストラットの長手軸の投影の各組と、ストラット接続点の投影との双方が採用される。図2Bは、第1のレントゲン画像220上に描画される、また、第2のレントゲン画像230上に描画される、あるストラット205の投影の決定を図示する。その後、本実施形態は、複数のストラット205の事前に決定された長さと、第1のレントゲン画像220上に描画される複数のストラット205の投影の第1の組とを用いて、X線撮像装置に対する、X線源202の第1の3D位置と、物体210の第1の3D位置とを決定してもよい。その後、同様に、本技術は、複数のストラット205の事前に決定された長さと、第2のレントゲン画像230上に描画される複数のストラット205の投影の第2の組とを用いて、X線撮像装置に対する、X線源202の第2の3D位置と、物体210の第2の3D位置とを決定してもよい。X線源202および物体210の、第1および第2の3D位置を決定するために、様々な数学的モデルを採用することができる。これらのモデルは、後により詳しく記述される。
X線源202および物体210の、第1および第2の3D位置が決定されると、その後、本技術は、第1および第2の向き220,230において、X線撮像装置204に対し、複数のストラット205の3D位置を用いて、3D参照フレームにおいて第1および第2の3D物体投影を位置合わせしてもよい。その後、本実施形態は、第1および第2の3D物体投影に基づいて、撮像された身体部分の3Dモデルを3D参照フレームにおいて作成してもよい。別の実施形態は、物体210の3Dモデルを、単独で、または、撮像された身体部分の3Dモデルに加えて、作成してもよい。
[モデル5]
本明細書に開示される技術の一例示的実施形態では、複数のストラットは5つのストラットを含んでもよい。図5は、ある向きにおける、2つの固定部材502と5つのストラット506とを備える物体510を示す。本実施形態では、X線源501の第1の3D位置と、物体510の第1の3D位置とは、以下に論じる数学的モデルを用いて決定されてもよい。この数学的モデルでは、(x,y,z)はX線源501の座標であり、(x0,y0,z0)〜(x9,y9,z9)は、固定部材502への10個のストラット接続点503の座標であり、(X0,Y0,Z0)から(X1,Y1,Z1)まで、…、(X8,Y8,Z8)から(X9,Y9,Z9)までは、第1のレントゲン画像上のストラット506の長手軸の投影の第1の組であり、l01,l02,l03,l04,l06,l07,l08,l09,l12,l13,l14,l15,l16,l17,l18,l19,l23,l24,l25,l26,l27,l28,l29は、10個のストラット接続点503間の事前に決定された距離であり、siは未知の比率であり、v(i−1)ix,v(i−1)iy,v(i−1)izはストラット506の長手軸の投影のベクトルであり、ただしv(i−1)ix=Xi−X(i−1),v(i−1)iy=Yi−Y(i−1),v(i−1)iz=Zi−Z(i−1)であると想定される。ここでiは0から9までの数である。したがって、たとえばv01x=X1−X0,v01y=Y1−Y0,v01z=Z1−Z0である。このように33個の未知のパラメータがあり、X線源501および物体510の位置を決定するためには33個の方程式の系が必要である。
さらに、ストラット接続点の投影の座標は、次の方程式によって表せる:
X座標: X0+v01x*s0
Y座標: Y0+v01y*s0
Z座標: Z0+v01z*s0
他のストラット接続点503の座標は同様に決定され、結果として43個の未知のパラメータとなる。したがって、43個の未知のパラメータを持つ数学的モデルを解くために43個の方程式が必要となる。しかしながら、この数学的モデルは無制限な数の解を生成する可能性がある。本実施形態は、解の数を制限するために、投影上に単一の区別可能な点(distinguishable point)を決定してもよい。ストラット端点の1つの投影が、ユーザにより手作業で決定されてもよい。たとえば、ユーザは、s9=0である場合に点の投影(X9,Y9,Z9)を決定し、これにより42個の方程式の数学的モデルが解けるようにしてもよい。
本実施形態では、X線源501と、ストラット接続点503と、ストラット接続点の投影とは、同一の直線上に位置している。したがって、点(x0,y0,z0)について、次の方程式の対が決定可能である:
Figure 2016513529
本実施形態では5つのストラットがあり、各ストラットについて上記の方程式の対があるので、ストラット接続点を記述する方程式が10個あることになる。ストラット接続点の位置に関与する追加の22個の方程式を、ピタゴラスの定理を用いて導出してもよい。したがって、結果として次の方程式が得られてもよい:
Figure 2016513529
しかしながら、この方程式の組は、各マーカー相互の相対的な位置を記述する方程式をすべて含むわけではない。したがって、系に含まれない方程式に関連して、X線源501および物体510の3D位置に対する方程式の組の解をチェックする必要がある:
Figure 2016513529
X線源501の第2の3D位置と、物体510の第2の3D位置とは、以下に論じる、実質的に同様な数学的モデルを用いて決定されてもよい。この数学的モデルでは、(‘x,’y,’z)はX線源301の座標であり、(´X0,´Y0,´Z0)から(´X1,´Y1,´Z1)まで、…、(´X8,´Y8,´Z8)から(´X9,´Y9,´Z9)までは、第2のレントゲン画像上のストラット506の長手軸の投影の第2の組であり、‘siは未知の比率であり、‘v(i−1)ix,‘v(i−1)iy,‘v(i−1)izはストラット506の長手軸の投影のベクトルであり、ただし‘v(i−1)ix=’Xi−‘X(i−1),‘v(i−1)iy=’Yi−‘Y(i−1),‘v(i−1)iz=’Zi−‘Z(i−1)であると想定される。ここでiは0から9までの数である。したがって、たとえば‘v01x=’X1−‘X0,‘v01y=’Y1−‘Y0,‘v01z=’Z1−‘Z0である。
Figure 2016513529
しかしながら、この方程式の組は、ストラット506の位置およびストラット接続点503の位置を記述する方程式をすべて含むわけではない。したがって、系に含まれない方程式に関連して、X線源501および物体510の3D位置に対する方程式の組の解をチェックする必要がある:
Figure 2016513529
これにより、X線源501および物体510の、第1および第2の3D位置を求めてもよい。
特定の状況下では、X線源501のZ座標の変化は、レントゲン画像上の各マーカーの投影の比例的変化(proportional change)を生じる。これらの状況下では、X線源501のZ座標は、方程式の組を解く際に解く必要のない定数パラメータとして設定することができる。ここで、Z座標に対する定数パラメータは、物体510がX線源501とレントゲン画像との間にフィットできるようにする大きな数とすべきである。これにより、以下のモデル6〜8に例示するように、X線源501の3D位置と物体510の3D位置とを決定するために、より少ない数のマーカーを用いることができるようになる。
[モデル6]
本技術の本実施形態では、複数のストラット506は4つのストラット506を含んでもよい。本実施形態では、X線源501の第1の3D位置と、物体510の第1の3D位置とは、以下に論じる数学的モデルを用いて決定されてもよい。この数学的モデルでは、(x,y,z)はX線源501の座標であり、(x0,y0,z0)〜(x7,y7,z7)は、8つのストラット接続点503の座標であり、(X0,Y0,Z0)から(X1,Y1,Z1)まで、…、(X6,Y6,Z6)から(X7,Y7,Z7)までは、第1のレントゲン画像上のストラットの長手軸の投影の第1の組であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、8個のストラット接続点503間の事前に決定された距離であり、siは未知の比率であり、v(i−1)ix,v(i−1)iy,v(i−1)izはストラットの長手軸の投影のベクトルであり、ただしv(i−1)ix=Xi−X(i−1),v(i−1)iy=Yi−Y(i−1),v(i−1)iz=Zi−Z(i−1)であると想定される。ここでiは0から7までの数である。したがって、たとえばv01x=X1−X0,v01y=Y1−Y0,v01z=Z1−Z0である。変数のこのような関係は、さらに図6に示される。
さらに、(´X0,´Y0,´Z0)〜(´X7,´Y7,´Z7)は第1のレントゲン画像上の8つのストラット接続点503の投影の第1の組であると想定され、これは次のように表せる:
X’0=X0+v01x*s0
Y’0=Y0+v01y*s0
Z’0=Z0+v01z*s0
同様の方程式が、固定部材502に対して他のストラット接続点503について決定される。
3点(x,y,z),(xi,yi,zi),(X’i,Y’i,Z’i)は同一の直線上に位置している。したがって、この直線の方程式は次のように表せる:
Figure 2016513529
この方程式は、代替的に、2つの方程式により表現することもできる:
Figure 2016513529
固定部材502に対して8つのストラット接続点503があり、ストラット接続点503それぞれについて上記の方程式の対があるので、ストラット506の長さを記述する方程式が16個あることになる。ストラット接続の位置に関与する追加の18個の方程式を、ピタゴラスの定理を用いて導出してもよい。
Figure 2016513529
X線源501の第2の3D位置と、物体510の第2の3D位置とは、以下に論じる、実質的に同様な数学的モデルを用いて決定されてもよい。この数学的モデルでは、(‘x,’y,’z)はX線源501の座標であり、(‘X0,’Y0,’Z0)から(‘X1,’Y1,’Z1)まで、…、(‘X6,’Y6,’Z6)から(‘X7,’Y7,’Z7)までは、第2のレントゲン画像上のストラット506の長手軸の投影の第2の組であり、(´X´0,´Y´0,´Z´0)〜(´X´7,´Y´7,´Z´7)は、第2のレントゲン画像上の8つのストラット接続点503の投影の第2の組の座標であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、は、8個のストラット接続点503間の事前に決定された距離であり、siは未知の比率であり、‘v(i−1)ix,‘v(i−1)iy,‘v(i−1)izはストラット506の長手軸の投影のベクトルであり、ただし‘v(i−1)ix=’Xi−‘X(i−1),‘v(i−1)iy=’Yi−‘Y(i−1),‘v(i−1)iz=’Zi−‘Z(i−1)であると想定される。ここでiは0から7までの数である。したがって、たとえばv01x=X1−X0,v01y=Y1−Y0,v01z=Z1−Z0である。
Figure 2016513529
これにより、X線源501および物体510の、第1および第2の3D位置を求めてもよい。
[モデル7]
本技術の別の実施形態では、複数のマーカーは4つのストラットを含んでもよい。この実施形態では、X線源501の第1の3D位置と、物体510の第1の3D位置とは、以下に論じる数学的モデルを用いて決定されてもよい。この数学的モデルでは、(x,y,z)はX線源501の座標であり、(x0,y0,z0)〜(x7,y7,z7)は、固定部材502に対する8つのストラット接続点503の座標であり、(X0,Y0,Z0)から(X1,Y1,Z1)まで、…、(X6,Y6,Z6)から(X7,Y7,Z7)までは、第1のレントゲン画像上のストラット506の長手軸の投影の第1の組であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、8個のストラット接続点503間の事前に決定された距離であり、siは未知の比率であり、v(i−1)ix,v(i−1)iy,v(i−1)izはストラットの長手軸の投影のベクトルであり、ただしv(i−1)ix=Xi−X(i−1),v(i−1)iy=Yi−Y(i−1),v(i−1)iz=Zi−Z(i−1)であると想定される。ここでiは0から7までの数である。したがって、たとえばv01x=X1−X0,v01y=Y1−Y0,v01z=Z1−Z0である。変数のこのような関係は図6に示される。
さらに、(X´0,Y´0,Z´0)〜(X´7,Y´7,Z´7)は第1のレントゲン画像上の8つのストラット接続点503の投影の第1の組であると想定され、これは次のように表せる:
X’0=X0+v01x*s0
Y’0=Y0+v01y*s0
Z’0=Z0+v01z*s0
同様の方程式が、固定部材502に対して他のストラット接続点503について決定される。
本実施形態は、次のパラメータ方程式を採用する:
x=x0+α*t; y=y0+β*t; z=z0+γ*t
ここで、α,β,γは方向ベクトルを示し、tはある直線上の点(x,y,z)を別の点(たとえば(x0,y0,z0))に対して特徴づけるパラメータである。3点(x,y,z),(xi,yi,zi),(X’i,Y’i,Z’i)は、X線源601からマーカー603へ、そして撮像装置604へと通過する同一の直線上に位置している。ここでiは0から7までの数である。したがって、この直線のパラメータ方程式は、次のように表せる:
xi=x+α*t; yi=y+β*t; zi=z+γ*t
ただしα=X´i−x; β=Y´i−y; γ=Z´i−z
本実施形態は、固定部材502上のストラット接続点503間の距離についての適切な数の方程式(既知の幾何学的配置の固定器の4つのストラット506の8つのストラット接続点503間の18個の距離について18個の方程式)を追加することにより、18個の方程式を決定する。結果として得られる18個の方程式は、次のように提供される:
Figure 2016513529
本実施形態は、t0〜t7,s0〜s7およびX線源501の第1の3D位置を決定してもよく、続いてストラット接続点503の座標を計算してもよい。X線源501の第2の3D位置と物体510の第2の3D位置とは、以下に論じる、実質的に同様な数学的モデルを用いて決定されてもよい。この数学的モデルでは、(‘x,’y,’z)はX線源501の座標であり、(‘X0,’Y0,’Z0)から(‘X1,’Y1,’Z1)まで、…、(‘X6,’Y6,’Z6)から(‘X7,’Y7,’Z7)までは、第2のレントゲン画像上のストラット506の長手軸の投影の第2の組であると想定される。
さらに、(´X´0,´Y´0,´Z´0)〜(´X´7,´Y´7,´Z´7)は第2のレントゲン画像上の8つのストラット接続点503の投影の第2の組の座標であると想定され、これは次のように表せる:
‘X’0=‘X0+‘v01x*‘s0
‘Y’0=‘Y0+‘v01y*‘s0
‘Z’0=‘Z0+‘v01z*‘s0
同様の方程式が、固定部材502に対して他のストラット接続点503について決定される。
本実施形態は、次のパラメータ方程式を採用する:
´x=x0+α*t; ´y=y0+β*t; ´z=z0+γ*t
3点(‘x,‘y,‘z),(xi,yi,zi),(‘X’i,‘Y’i,‘Z’i)は、X線源601からマーカー603へ、そして撮像装置604へと通過する同一の直線上に位置している。したがって、この直線のパラメータ方程式は、次のように表せる:
xi=´x+α*t; yi=´y+β*t; zi=´z+γ*t
ただしα=´X´i−´x; β=´Y´i−´y; γ=´Z´i−´z
本実施形態は、リング上のストラット接続点間の距離についての適切な数の方程式を追加することにより、18個の方程式を決定する。結果として得られる18個の方程式は、以下に提供される:
Figure 2016513529
これによって、X線源および物体の、第1および第2の3D位置を求めてもよい。
[モデル8]
本技術の別の実施形態では、複数のストラット506は4つのストラットを含んでもよい。この実施形態では、X線源501の第1の3D位置と、物体510の第1の3D位置とは、以下に論じる数学的モデルを用いて決定されてもよい。この数学的モデルでは、(x,y,z)はX線源501の座標であり、(x0,y0,z0)〜(x7,y7,z7)は、固定部材502に対する8つのストラット接続点503の座標であり、(X0,Y0,Z0)から(X1,Y1,Z1)まで、…、(X6,Y6,Z6)から(X7,Y7,Z7)までは、第1のレントゲン画像上のストラット506の長手軸の投影の第1の組であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、8個のストラット接続点503間の事前に決定された距離であると想定される。変数のこのような関係は図6に示される。
さらに、各ストラット接続点503は、X線源501と、当該ストラットの長手軸の投影とを通る平面上に配置される。このように、ストラット接続点503のそれぞれについて1平面付属方程式(one plane appurtenant equations)が決定され、以下に提供するように、4つのストラット506を持つ固定器について8つの方程式の系を作成してもよい。
Figure 2016513529
ここで、本実施形態は、ストラット接続点間の距離についての適切な数の方程式(たとえば、既知の幾何学的配置固定器の4つのストラットの8つのストラット接続点間の18個の距離について18個の方程式)を追加することにより、次のような26個の方程式の数学的モデルを提供する:
Figure 2016513529
X線源501の第1の3D位置と、物体510の第1の3D位置とは、上記の方程式を解くことにより決定されてもよい。X線源501の第2の3D位置と、物体510の第2の3D位置とは、以下に論じる、実質的に同様な数学的モデルを用いて決定されてもよい。この数学的モデルでは、(‘x,’y,’z)はX線源501の座標であり、(´X0,´Y0,´Z0)から(´X1,´Y1,´Z1)まで、…、(´X6,´Y6,´Z6)から(´X7,´Y7,´Z7)までは第2のレントゲン画像上の4つのストラットの投影の第2の組であると想定される:
Figure 2016513529
これにより、X線源501および物体510の、第1および第2の3D位置を決定してもよい。
[X線源の位置の近似]
場合によっては、第1および第2のレントゲン画像の解像度は、マーカーによって第1および第2のレントゲン画像上に作成された影の位置を正確に特定できるようにするには充分でないということが理解される。図2Bを参照すると、小さな誤差が存在してベクトル/軌跡(たとえば242および244)の位置合わせを誤らせ交差しないようにする可能性がある。そのような場合には、X線源202の位置は近似モデルを用いて決定してもよい。例示的な近似モデルでは、ベクトル/軌跡242および244の間の線分(segment、セグメント)の向きおよび位置を決定してもよく、X線源202の位置を表すために線分上の1点を選択してもよい。一実施形態では、選択される線分は、ベクトル/軌跡242および244の共通垂線であってもよく、この共通垂線の中点がX線源202の位置を表すものとして選択されてもよい。ベクトル/軌跡242および244の双方の共通垂線は、ベクトル/軌跡242および244の間の最短線分であってもよく、X線源202の正確な近似を可能にしてもよいし、近似モデルの望ましい精度に応じて他の線分を選択してもよいということが理解される。
一例示的実施形態では、ベクトル/軌跡242および244の共通垂線は、以下に論じる数学的モデルを用いて決定されてもよい。この数学的モデルでは、(x ,y ,z )はマーカー1影(250)の座標であり、(x ,y ,z )はマーカー1(252)の座標であり、(x ,y ,z )はマーカー2影(254)の座標であり、(x ,y ,z )はマーカー2(256)の座標であると想定される。したがって、第1の直線242についての方程式は、次のように表せる:
Figure 2016513529
また、第2の直線244についての方程式は、次のように表せる:
Figure 2016513529
第1の直線242および第2の直線244から結果として得られるベクトルは、それぞれ次のように表せる:
Figure 2016513529
ただし
Figure 2016513529
である。
ベクトルaおよびbを、次の方程式に従って乗算すると、直線242および244の双方に対して垂直なベクトルcが提供される:
Figure 2016513529
ただしi,j,kは座標軸x,y,zに沿う向きの単位ベクトルである。
Figure 2016513529
一実施形態では、X線源202の位置の近似は、ベクトルc内に存在して直線242および244を接続する線分Sを定義することを含んでもよい。したがって線分Sは直線242および244の共通垂線である。これを行う方法の一つは、マーカー1影(250)と、第1の直線242と、ベクトルcとを含む平面Dを構築することである。そのような平面Dに垂直なベクトルは、ベクトル乗算
Figure 2016513529
の積であり、次のように表せる:
Figure 2016513529
このベクトルは、単位長に対して正規化可能であり、次のように表せる:
Figure 2016513529
したがって、座標
Figure 2016513529
を有するマーカー1影(250)を通過し、垂直ベクトル
Figure 2016513529
を有する平面Dは、次の方程式によって表せる:
Figure 2016513529
線分Sの端点の1つは、平面Dが直線244と交わる交点であってもよい。この交点の位置を決定するために、斜辺Gが直線244に沿って延び、マーカー2影254と直線244が平面Dに交わる交点とを接続するような、直角三角形を描いてもよい。さらに、この直角三角形の、斜辺でない第1の辺Rは、平面Dに垂直でマーカー2影254から平面Dに延びるベクトルrによって定義されてもよい。直角三角形の、斜辺でない第2の辺は、平面D内の斜辺Gの投影によって定義されてもよい。
斜辺でない第1の辺Rの長さ(マーカー2影254と平面Dとの距離)は、平面Dの正規化された垂直ベクトルnと、ベクトルrとのスカラー乗算(scalar multiplication)により決定されてもよい。この場合には、このスカラー乗算の積は、次の方程式に示すように「マーカー2影254」の座標について表せる:
Figure 2016513529
さらに、ベクトルrとベクトルbとがなす角φの余弦は、次のように表せる:
Figure 2016513529
したがって、斜辺Gの長さは、斜辺でない第1の辺Rの長さを、この辺Rと斜辺Gとがなす角の余弦で除算することによって決定することができる:
Figure 2016513529
直線244が平面Dと交わる交点の座標を求めるために、マーカー2影254および長さGから、直線244に沿って延びるベクトル
Figure 2016513529
を定義することができる。
Figure 2016513529
ただし
Figure 2016513529
である。
これらの座標は、線分Sの端点の1つを定義する。線分Sの第2の端点の座標を求めるために、同様の計算を実行してもよい。一実施形態では、直線244に沿って平面を定義してもよく、この平面の直線244上の交点を求めてもよい。一実施形態では、線分の両端を定義した後に、X線源202の位置が線分Sの中点に配置されていると近似してもよく、それらの座標の平均値として計算されてもよい:
Figure 2016513529
他の実施形態では、X線源202が近似される位置は、線分Sの両端間の任意の位置としてよいということが理解される。さらに、上に論じた例示的な数学的モデルは、X線源202の位置を近似する効率的かつ正確な方法を提供するが、X線源の位置を近似するために、本開示の原理による他の適切なモデルを用いてもよいということが理解される。
第1および第2の撮像の向き(220,230)におけるX線源202の3D位置が特定されると、様々な技術を用いて撮像された物体の3Dモデルを作成することができる。一実施形態では、第1の撮像の向き220と、第2の撮像の向き230との間の、撮像軸I周りの角変位量が既知である。この実施形態に対応する例は、図12A〜図12Gに図示される。図12Aは、2つの向き(1220,1230)において撮像される物体1201を図示する。2つの相対的向き(1220,1230)における画像は、撮像される物体1201を撮像軸Iの周りに角変位αだけ回転させることにより、または、X線源および撮像装置を撮像軸Iの周りに角変位αだけ回転させることにより、準備することができる。撮像軸Iは、第1の向き1220のX線撮像装置(図示せず)の平面と、第2の向き1230のX線撮像装置の平面とに平行であることが好ましいが、これは必須ではない。2つの向きにおいてこれらの画像を作成することにより、結果として、向き1220,1230にそれぞれ対応する2つのレントゲン画像(1202,1204)が得られる。また、図12Aには、レントゲン画像(1202,1204)それぞれに対するX線源1212の相対的位置が示されている。これらのX線源1212の3D位置は、上述のようにレントゲン画像(1202,1204)上に参照マーカーにより作成される影に基づいて決定してもよいし、当該技術分野で既知の任意の他の技術(たとえば撮像装置104に対するX線源1212の位置を物理的に測定する等)により決定してもよい。
物体1201の3Dモデルの作成における別のステップは、レントゲン画像における撮像された物体1201のアウトラインを決定することである。このコンセプトは図12Bに図示されており、図12Bでは、第1のレントゲン画像1202における撮像された物体のアウトラインがアウトライン1213として特定されている。同様に、第2のレントゲン画像1204における撮像された物体のアウトラインがアウトライン1215として特定されている。レントゲン画像はコンピュータシステムに記憶されるデジタル画像であるが、この処理は画像処理ソフトウェアを用いることにより自動的に実行可能である。別の実施形態によれば、この処理は、レントゲン画像における撮像された物体のアウトラインを、マウス、スタイラスまたは他の任意のトレース装置でトレースすることにより、手作業で実行可能である。撮像された物体のアウトラインと、X線源1212の3D位置とを決定した後に、撮像された物体のアウトラインの投影を作成することができる。第1の向き120における物体アウトライン1213の投影は、図12Bにおいて、第1のレントゲン画像1202におけるアウトライン1213から第1の向き1220におけるX線源1212の3D位置まで通過する投影線308によって示される。同様に、第2の向き1230における物体アウトライン1215の投影は、図12Bにおいて、第2のレントゲン画像1204におけるアウトライン1215から第2の向き1230におけるX線源1212の3D位置まで通過する投影線1216によって示される。
第1および第2のレントゲン画像(1220,1230)それぞれについて撮像された物体の投影が作成されると、それらの投影が互いにどのように交差するかを、向き(1220,1230)の互いに対する相対的位置を用いて決定することができる。これは様々な方法で行うことができる。一実施形態では、各3D投影を、図12A〜図12Eに図示されるx,y,z参照フレーム1250に対応する単一の3D参照フレーム内に組み合わせてもよい。x,y,z参照フレーム1250の原点は、撮像軸Iに沿って、第1の向き1220におけるX線源1212からのX線が撮像軸と垂直に交差し、第2の向き1230におけるX線源1212からのX線が撮像軸と垂直に交差する点1251に配置されてもよい。上述のように、この参照フレーム1250において、角度αは、x軸または撮像軸Iの周りの、2つの向き(1220,1230)の間の角変位量に対応する。この角度αは様々な方法で決定することができる。X線源1212および撮像装置が静止しており物体が回転する実施形態では、図2Aおよび図2Bに示すように、角度αはx軸または撮像軸Iの周りの物体の回転量に対応する。物体が静止しているがX線源1212および撮像装置が物体の周りに回転する代替的実施形態では、図2Aおよび図2Bに示すように、角度αは、撮像軸Iの周りのX線源1212および撮像装置の回転量に対応する。X線源1212および撮像装置を撮像される物体1201の周りに回転する時には、X線源1212の位置が撮像装置の位置に対して固定されていることが好ましい可能性があるが、これは必須ではない。さらに、レントゲン画像(1202,1204)は、互いに実質的に直交する向き(1220,1230)において取得されることが好ましい可能性があるが、これは必須ではない。上記の角変位αとアウトラインの投影とを用いて、レントゲン画像(1202,1204)の相対的位置と、それらの対応する投影とを互いに位置合わせすることができる。
角度βおよびγは、それぞれ、z軸およびy軸周りの、第2のレントゲン画像1204に対する第1のレントゲン画像1202の角変位に対応する。上述のように、実施形態によっては、第1および第2の相対的向き1220および1230が互いに実質的に直交しており、それらの実施形態では角度βおよびγは実質的に0であってもよい。第1および第2の相対的向き1220および1230が互いに実質的に直交していない実施形態では、第1および第2のレントゲン画像(1202,1204)はさらに、様々な手法(本開示に記載される反復的手法を含む)を用いて、角度βおよびγで位置合わせすることができる。レントゲン画像(1202,1204)を角度βおよびγで位置合わせすることは任意選択的であるが、位置合わせすることによってより正確な物体1201の3Dモデルが可能になる可能性があるということが理解される。
図12Cは、撮像された物体1201のレントゲン画像(1202,1204)と、撮像された物体1201の3D物体投影とが、複数の平面(平面1231,1232,1234を含む)と交差し得るということを図示する。これらの平面は、それぞれ、第1の向き1220におけるX線源1212の位置に対応する3D参照フレームにおける第1の3D位置と、第2の向き1230におけるX線源1212の位置に対応する3D参照フレームにおける第2の3D位置とを通過する。また、図12Cは、3D参照フレームにおいて第1および第2の3D位置を通過する傾斜軸1236を図示する。平面1231,1232,1234は、傾斜軸1236の周りにそれぞれ異なる傾斜を有し、第1および第2のレントゲン画像(1202,1204)において撮像された物体1201のアウトラインと交差する。平面1231と、第1および第2のレントゲン画像(1202,1204)における画像アウトラインとの交点の位置は、点318でマークされている。平面1232と、第1および第2のレントゲン画像(1202,1204)における画像アウトラインとの交点の位置は、点1240でマークされている。平面1234と、第1および第2のレントゲン画像(1202,1204)における画像アウトラインとの交点の位置は、点1244でマークされている。
図12Dでは、交点1238,1240,1242は、それぞれ、対応する向き(120,150)で3D参照フレーム1250においてX線源1212の位置に接続されている。したがって、レントゲン画像1202における交点1238,1240,1242は、直線324によって、3D参照フレーム1250における第1の3D位置に接続されており、これは第1の向き1220におけるX線源1212の位置に対応する。同様に、レントゲン画像1204における交点1238,1240,1242は、直線1256によって、3D参照フレーム1250における第2の3D位置に接続されており、これは第2の向き1230におけるX線源1212の位置に対応する。さらに、第1および第2のレントゲン画像(1202,1204)における点の組1238と交差する4本の直線は、また、3D参照フレーム1250における多角形1258を形成するよう互いに交差する。同様にして、第1および第2のレントゲン画像(1202,1204)における点の組1240と交差する4本の直線は、また、3D参照フレーム1250における多角形1258を形成するよう互いに交差する。さらに、第1および第2のレントゲン画像(1202,1204)における点の組1242と交差する4本の直線は、また、3D参照フレーム1250における多角形1258を形成するよう互いに交差する。
図12Eにおいて、多角形1258を定義する処理は、十分な解に到達するか、または撮像された物体1201との交点がそれ以上発見されなくなるまで、傾斜軸1236と位置合わせされた1つ以上の平面について繰り返されてもよい。これらの多角形1258は、それぞれ、参照フレーム1250における、第1および第2の向き(1202,1204)における3D物体投影の間の交差部分に対応する。
3D投影の交差部分に対応する多角形1258の系列を作成した後に、各多角形1258を、図12Fに図示される撮像された物体1201の断面形状に対応する閉曲線(たとえば楕円)1290へと変換してもよい。多角形1258の系列を閉曲線1290に変換する前に、撮像された物体1201の概略的な形状および向きが既知であることが好ましいが、これは必須ではない。たとえば、撮像された物体(骨等)が、概して楕円の断面形状を有している場合には、各多角形1258は、その多角形内に配置される閉曲線1290(楕円等)によって置き換えることが可能である。一方、撮像された物体が非対称な形状を有している場合には、物体の正確な3Dモデルを作成するために、撮像された物体についての他の情報(たとえば、形状、断面、向き、等)を用いることができる。各多角形1258が対応する形状(たとえば閉曲線または楕円)1290によって置き換えられると、これらの形状を接続する表面を準備することができる。この表面は、図12Gに図示するように、撮像された物体の正確な3Dモデル1292を表す場合がある。いくつかの実施形態では、3Dモデル332の精度は、画像ライブラリに格納された既知の形状に従ってモデル1292を修正することによって改善することができる。
上述のように、第1および第2の相対的な向き1220および1230が実質的に直交しない場合には、角度βおよびγは、本開示の原理に従って反復的手法を用いて決定されてもよい。一例示的実施形態では、レントゲン画像(1202,1204)は、まずレントゲン画像(1202,1204)を既知のαで位置合わせし、その後、レントゲン画像(1202,1204)を様々な角度βおよびγで位置合わせすることによって撮像された物体1201の様々なテスト3Dモデルを作成し、最後に第1および第2のレントゲン画像1202および1204において撮像された物体1201のアウトラインに実質的にマッチする2D投影を作成する3Dモデルを特定することによって、角度βおよびγに向けられてもよい。物体1201のテストモデルは、精度をより良くするために、図12C〜図12Gに関連して上述された手法に従って作成されてもよい。しかしながら、当該技術分野で既知の任意の適切なモデル化技術に従って様々なモデルを生成してもよいということが理解される。
別の実施形態では、物体の3Dモデルは、2つの撮像の向き(1220,1230)の間の角変位αが既知でない時であっても、固定された参照フレームにおいて作成可能である。この実施形態に対応する例は、図13A〜図13Eに図示される。直前に述べた、角変位αが既知である処理と同様に、物体の2つのレントゲン画像が異なる複数の向きにおいて準備される。これらのレントゲン画像は、それぞれ、物体マーカーが取り付けられた物体1301の画像を含む。少なくとも1つの物体マーカー1342は、物体1301に直接的または間接的に取り付けられてもよく、物体マーカー1342の数は、それぞれが含む基準(fiducial)の数に応じて変化してもよい。一例示的実施形態では、少なくとも1つの物体マーカー1342に、全部で少なくとも3つの基準が含まれる。この少なくとも3つの基準は、1つの物体マーカー1342に封入されてもよく(enclosed)、または複数の物体マーカー1342(たとえば2つまたは3つの物体マーカー1342等)にわたって分散してもよい。別の例示的実施形態では、少なくとも1つの物体マーカー1342に、全部で4つ以上の基準が含まれる。この4つ以上の基準は、1つのマーカー1342に封入されてもよく(enclosed)、または、基準の具体的な数に応じて、複数(たとえば2つ、3つ、4つ、またはそれ以上)の物体マーカー1342にわたって分散してもよい。
代表的な物体マーカー1342が取り付けられた物体1301の例が図13Aに図示される。図13Aでは、撮像された物体1301は、それぞれに取り付けられる4つの物体マーカー1342を含み、これらの物体マーカー1342はそれぞれ基準1344を含む。固定された参照フレームにおいて物体1301のモデルを作成するために物体マーカー1342を用いる第1の手法では、少なくとも3つの基準1344が直接的または間接的に少なくとも1つの撮像された物体1301に取り付けられていれば、物体マーカー1342の数および種類は異なってもよい。たとえば、一実施形態では、1つの物体マーカー1342が3つの基準1344を備えてもよい。別の例示的実施形態では、それぞれ2つの基準1344を備えた物体マーカー1342を2つ用いてもよい。さらに別の例示的実施形態では、基準1344を1つずつ備えた物体マーカー1342を3つ用いてもよい。いくつかの実施形態では3つの基準1344を用いてもよいが、以下に論じる理由から、4つ以上の基準1344を用いる実施形態がより望ましいということが理解される。さらに、固定された参照フレームにおいて物体1301のモデルを作成するために物体マーカー1342を用いる第1の手法によれば、基準1344の互いに対する位置は事前に決定されているということが理解される。一例示的実施形態では、基準1344間の線分の長さおよび向きを決定するために、測定値を取得してもよい。別の実施形態では、基準1344の互いに対する位置を事前に決定できるように、物体マーカー1342を事前に決定された向きに配置してもよい。この場合、基準1344の間の線分は数学的に決定可能である。
図13Aに図示する実施形態では、物体マーカー1342の各画像が、それぞれ対応するレントゲン画像1302および1304に示されている。2つのレントゲン画像を受信する際に、X線撮像装置に対するX線源1312の3D位置が、本開示に開示される原理に従って、撮像の向き(1320,1330)のそれぞれについて決定されてもよい。とくに、このような決定は、上記の各実施形態に関連して記載されたものと同一の方法において、参照マーカーおよび基準の使用に基づいてもよい。同様に、第1および第2のレントゲン画像(1302,1304)における、撮像された物体1301のアウトラインと、基準1344の影点1306とは、直前の実施形態において記載されたものと同一の技術を用いて特定されてもよい。ここで、撮像された物体1301の3Dモデルを準備するために、物体マーカー1342および基準1344を用いて、異なるステップを利用してもよい。
一般的に、固定された参照フレームにおける物体1301のモデルを作成するために物体マーカー1342を用いる第1の手法は、図13Bに示すように、レントゲン画像(1302,1304)における影点1306と、それらがそれぞれ対応する撮像の向き(1320,1330)におけるX線源1312の位置とを接続する投影線1340を構成することを含む。各レントゲン画像(1302,1304)に対する基準1344の相対的な3D位置は、投影線410の向きと、基準1344間の事前に決定された線分とに基づいて数学的に決定してもよい。また、2つの撮像の向き220,230間の角変位αは、基準1344の3D位置を固定された参照フレームにおいて位置合わせすることによって決定してもよい。変位角度αが特定されると、撮像された物体の3Dモデルを作成する処理は、図12C〜図12Gに関連して記載したものと同一の方法で進行してもよい。
また、各レントゲン画像(1302,1304)に対する基準1344の3D位置の決定は、様々な数学的手法に従って達成可能であるということが理解される。一例示的数学的手法は、図13Cを参照して探索される。上述のように、マーカー(図示せず)は、3つの基準1344の互いに対する位置が事前に決定されたものとなるように、物体(図示せず)に固定されてもよい。図13Cに例示する実施形態では、レントゲン画像304における基準1344の影点1306を用いて、投影線1340(これらは幾何学的に協働して三角錐(three-sided pyramid)を形成してもよい)を構成してもよい。さらに、対応する基準1344の互いに対する位置は事前に決定されているので、各基準1344の3D位置を接続することによって形成される三角形1308の諸寸法もまた数学的に決定可能である。この場合、図13Cに示すように、以下の幾何学的要素が確立可能である:光源1312の座標(L)、影点406の座標(K,M,N)、三角形1308の辺(a,b,c)の長さ。図13Cの例示的手法は、三角形1308の3D位置および向きを決定するために、三角形1308の諸寸法と錐の外側輪郭とがマッチする位置に到達するまで、三角形1308を錐内で回転させかつ「移動」させることを含んでもよい。既知の三角測量および三角法の技術に基づき、三角形1308の位置は、次の方程式の系の解に対応するものであってもよい:
Figure 2016513529
ここで、角度KLM,MLN,KLNはそれぞれα,β,γに対応し、x,y,zは、光源1312と基準1344との距離に対応する。この方程式の系は、数学的には8つの異なる解を持つが、それらの一部は複素数および負の数を含む可能性があり、したがって除外されてもよい。この場合、基準1344の位置を正しく反映し得る解が2つ残る可能性がある。しかしながら、残った2つの解のうちいずれが正しいかを数学的に決定するのは困難である。一実施形態では、双方の解に基づく画像物体の3Dモデルが人間に提示され、その人間が、撮像された物体の向きにマッチするモデルを視覚的に決定して選択してもよい。整形外科用途では、マッチするモデルを選択する人間は内科医であってもよい。
モデル化処理をより速めるために、上述のように正しいモデルを選択する人間の関与を、図13Dおよび図13Eに関連して開示される手法に従って低減または排除してもよい。一般的に、上述の数学的モデルは、追加の基準の考慮を含むよう変更されてもよい。例示される手法によれば、追加の基準1344が1つ用いられて、3つでなく4つの基準1344が用いられ、結果として、3つの影点1306からなる組が4通りあることに基づいて4つの異なる三角錐が構成される。三角錐のそれぞれについていくつかの解が利用可能な場合がある。一実施形態では、異なる複数の解が互いに比較され、当該技術分野で既知の数学的基準に従って最終的な解が選択されてもよい。たとえば、他の解からの偏差(deviation)が最も小さい解を、レントゲン画像1304に対する基準1344の3D位置を決定するために選択してもよい。別の例では、基準1344の3D位置を決定するために、すべての解の平均値を選択してもよい。
上述の各手法は、他のレントゲン画像1302に対する各基準1344の3D位置を決定するために繰り返されてもよい。そうすることによって、基準1344の3D位置を、2つの異なる座標系に対し、上述の手法に従って決定することができる。さらに、各基準を2つの座標系において位置合わせすることによって、図13Bに図示されるように、第1および第2のレントゲン画像の平行移動および回転の向き(x,y,z,α,β,γ)を、単一の固定された参照フレームにおいて決定することができる。いくつかの実施形態では、固定された参照フレームにおいて基準1344の3D位置を決定することによって、また、物体に対する、3D空間における基準1344の事前に決定された相対的な向きが与えられると、固定された参照フレームにおいて物体の3D位置を決定することができる。
上述の各例示的手法は、第2のレントゲン画像1304に対する第1のレントゲン画像1302の平行移動および回転の向き(x,y,z,α,β,γ)に対応する効率的かつ正確な方法を提供するために、3つまたは4つの基準1344を用いて実施してもよいが、本開示の原理に従う他の手法では、これら以外の数の基準1344を用いてもよいということが理解される。正確さおよび/または精密さをより高くできるようにするために、5つまたはそれ以上の基準を用いてもよい。たとえば、図13Aに示す実施形態では、8つの基準を用いてもよい。そのような場合には、3つの基準からなる組の組み合わせが56個存在してもよい。各組み合わせに対する少なくとも2つの可能な解によって、基準1344の位置について少なくとも112個の異なる可能な解が存在し得る。最終的な解は、基準1344の可能な位置すべての数学的分析に基づく、以下の例示的アルゴリズムに従って選択されてもよい:
1)上述のようにして取得された可能な解すべてに基づき、各基準1344の潜在的3D位置をすべて決定する。
2)各基準1344の平均3D位置を決定する。
3)ステップ2で決定された平均3D位置それぞれからの、各基準1344の潜在的3D位置すべての偏差を決定する。
4)ステップ2で決定された平均3D位置それぞれから最も大きくずれている3D位置に対応する、最も尤度の低い3D位置を特定する。
5)最も尤度の低い3D位置を導いた解を排除する。
6)各基準1344の、残る潜在的3D位置のそれぞれの、対応する平均3D位置からの偏差が、閾値(たとえば2mm、5mm、10mm、等)より小さくなるまで、ステップ1 5を繰り返す。
7)各基準1344の3D位置を、各基準1344の、残る潜在的3D位置のそれぞれの平均として近似する。
上記のアルゴリズムにより、基準1344の位置について正確な近似が可能になるということ、また、これは本明細書で論じた原理および当該技術分野で既知の任意の数学的技術に従って修正してもよいということが理解される。たとえば、一例示的実施形態では、このアルゴリズムは、各基準1344の潜在的位置間の分散を決定することを含み、平均および分散の双方から偏差に基づいて潜在的解を排除するよう、さらに修正されてもよい。
[実用的考慮]
整形外科装置が物体とともにレントゲン画像に撮像された後、手作業で、または適切なグラフィックソフトウェアを用いて、整形外科装置のアウトラインを決定してもよい。たとえば、内科医が手作業で整形外科装置のアウトラインを取り、そのような情報をコンピュータに入力してもよい。別の実施形態では、整形外科装置のアウトラインはパターン認識ソフトウェアによって自動的に生成されてもよい。その場合、整形外科装置のアウトラインは、本開示に従って、物体の3Dモデルを決定するために用いられてもよい。
いくつかの実施形態では、視認可能な影がデジタルのレントゲン画像上で2つ以上の画素にわたる可能性があるということが理解される。したがって、視認可能な影の正確な位置は、近似モデルを用いて近似されてもよい。図7は、一例示的近似モデル700の手法を図示するフローチャートである。近似モデル700は、複数の被制御画素領域(各被制御画素領域は複数の画素を含む)を定義するためのステップ702を含む。各被制御画素領域の複数の画素は、それぞれ、視認可能な影のそれぞれが配置される可能性が最も高い各位置に対応してもよい。たとえば、ある定義された被制御画素領域は、ある視認可能な影の周辺に3×3グリッドの9つの画素を含んでもよい。別の例では、ある定義された被制御画素領域は、ある視認可能な影の周辺に4×4グリッドの16個の画素を含んでもよい。例示的モデル700は、画素の様々な組(画素の各組は、それぞれ、定義された被制御画素領域からのピクセルを1つ含む)に基づいて、可能な影位置の複数の組み合わせを任意に割り当てるステップ704を含んでもよい。例示的モデル700は、割り当てられた影位置の組み合わせのそれぞれについて、望ましい参照点の位置を決定するステップ706を含んでもよい。たとえば、望ましい参照点は、リングの中心であってもよい。一実施形態では、望ましい参照点に対応する位置を決定するために、影位置の組み合わせがすべて割り当てられ、用いられる。別の実施形態では、望ましい参照点に対応する位置を決定するために、影位置の組み合わせのうち選択された組み合わせのみが割り当てられ、用いられる。例示的モデル700は、さらに、リングの中心の近似された位置を決定するために客観的判定基準を用いて望ましい参照点の第1および第2の位置を処理するステップ708を含んでも良い。例示的実施形態では、ステップ708のこの客観的判定基準は、当該技術分野において既知の、1つ以上の数学測定(平均、中央値、分散、標準偏差またはこれらの任意の組み合わせ、等)を含んでもよい。例示的実施形態では、差が0.01mmより大きい位置はすべて除去してもよい。選択された領域の組み合わせのいずれもその精度を持たない場合には、リングの中心の位置決定において差を最も小さくする組み合わせを用いてもよい。
図8は、合成3D座標系850において固定具(fixture)のモデルを示す概略図である。図2A〜図2Bに関連して上述したように、2つのレントゲン画像に基づいて第1および第2の3D座標系が独立に作成され、それぞれ第1および第2の平面801および802を含んでいる。この第1および第2の3D座標系が組み合わされて、合成3D座標系850が作成される。第1および第2の平面801および802は、第1および第2の3D座標系における第1および第2の参照点の座標が一致するような角度に位置合わせされている。
図9は、固定具900の第1リング901に結合された第1の物体セグメント910と、固定具900の第2リング902に結合された第2の物体セグメント920とのモデルである。このモデルは、上述の方法を用いて生成された合成3D座標系850に基づいている。いくつかの実施形態では、図8のモデルにより、第2のセグメント920の向きに対する第1の骨セグメント910の向きを決定することができる。とくに、このモデルにより、第2リング902に対する第1リング901の様々な向きに基づいて、第1および第2の骨セグメント910および920の相対的な向きを数学的に決定することができる。
本開示は、互いに実質的に直交する撮像の向き、または直交しない向きを用いることを説明してきた。これら2通りの実施形態からの選択は、様々な要素(器具の制限、および、特定の向きを撮像する関心またはその欠如を含む)に依存し得る。さらに、本開示の範囲と整合して、3つ以上の向きを利用してもよい。3つ以上の撮像の向きを用いることにより、フレームおよび組織の3Dモデルの精度を向上することができる。
フレームおよび組織セグメントの3Dモデルが作成されると、内科医または外科医は、より容易に、破損の性質と、望ましい結果を達成するために組織セグメントに作用させるべき固定、圧縮または伸延(または他の力)の程度とを理解することができる。ヘキサポッドリング固定器の3Dモデルは、望ましい固定、圧縮または伸延の命令が自動的に実行されるように、自動化されたフレーム制御装置に結合することもできると考えられる。
上述のように、物体の3Dモデルは物体のレントゲン画像から生成されてもよい。図10は、本開示の原理に従って、撮像された物体(図示せず)の3Dモデルをデジタルに生成するよう動作可能なシステム1000の概略図である。システム1000は、撮像された物体のレントゲン画像を受信するよう動作可能なコンピュータワークステーション1002を含んでもよく、コンピュータワークステーション1002は、様々な補助装置と通信する1つ以上のマイクロプロセッサ/コントローラを含んでもよい。一実施形態では、システム1000は、コンピュータワークステーション1002と通信するX線撮像装置1004を含んでもよく、X線撮像装置1004は、X線源1006から出て撮像される物体を通過するX線を受信するよう動作可能であってもよい。X線撮像装置1004は、レントゲン画像を直接的に生成するよう動作可能であってもよいし、画像データをコンピュータワークステーション1002に送信するよう動作可能であってもよい(その場合にはコンピュータワークステーション1002がX線画像を生成する)。別の実施形態では、システム1000は、ワークステーション1002と通信するスキャナ1008を含んでもよく、スキャナ1008は、X線フィルムをデジタル化されたレントゲン画像へとスキャンするよう動作可能であってもよい。いくつかの実施形態では、システム1000はさらに、ワークステーション1002と通信するディスプレイ1010を含んでもよく、ディスプレイ1010は、LCDディスプレイであってもよく、CRTディスプレイであってもよく、当該技術分野において既知の他の任意の表示装置であってもよい。ワークステーション1002は、デジタル化されたレントゲン画像をユーザに対してディスプレイ1010上に表示するよう構成されてもよく、ユーザは、本開示におけるような、表示されたレントゲン画像に関連する様々なデータ(マーカーまたはストラットの位置、マーカーまたはストラットの互いに対する事前に決定された位置、等)を入力してもよい。一例示的実施形態では、システム1000は、ワークステーション1002と通信する1つ以上の入力装置1012(マウス、ライトペンおよび/またはキーボード等)を含み、ユーザはこの入力装置1012を用いてデータを入力してもよい。ワークステーション1002のマイクロプロセッサまたはコントローラは、ユーザ入力データおよび画像データに基づき、本開示に従って、撮像された物体の3Dモデルを生成してもよい。いくつかの実施形態では、システム1000はさらに、ユーザに様々なモデルデータ、計算結果、画像またはグラフィックスを提供するよう動作可能な出力装置1014(プリンター等)を含んでもよい。システム1000はさらに、様々なモデルデータ、計算結果、画像またはグラフィックスを、後に使用するために格納するための記憶モジュール1016を含んでもよい。
図11は、図10に例示されるシステムの応用の概略図である。本実施形態は、本開示の原理に従って、撮像された物体の3Dモデルをデジタル的に生成する。患者がテーブル1101に座るかまたは横たわっていてもよい。他の実施形態では、患者はテーブル1101でなく椅子に座っていてもよい。テーブル1101の表面は、撮像装置と同一の材料から構成されてもよい。代替的に、撮像装置は、テーブル1101の表面の上、患者の下に配置されてもよい。X線およびレントゲン画像によって検査されるべき領域(本実施形態では脚)は、整形外科固定器1109で取り囲まれている。整形外科固定器1109で取り囲まれた脚は、本開示の原理に従って、X線装置1102によってX線照射される。整形外科固定器1109で取り囲まれた脚は、x,y,z方向において回転可能なX線装置1102によって、異なる複数の向きからX線照射されてもよい。
X線照射されたデータは、インターネットまたは他の適切なネットワークを介して、ケーブル(図示せず)を通じてまたは無線で、ユーザのローカルマシン1107に送信される。ユーザのローカルマシン1107は、本実施形態では通常のデスクトップコンピュータであるが、図10のコンピュータワークステーション1002として例示されるような任意の計算装置であってもよい。ユーザのローカルマシン1107は、X線データを受信し、処理し、記憶するためのプロセッサおよびメモリを備えてもよい。ユーザのローカルマシン1107は、X線データを画像として表示するディスプレイ1103に接続されてもよい。ユーザのローカルマシン1107は、マウス1104、キーボード(図示せず)、およびスキャナ/プリンタ1108に接続されてもよい。スキャナ/プリンタ1108は、X線画像をスキャンするか、または送信されたX線データを印刷するよう動作可能である。X線データは、また、CD、ユニバーサルシリアルバス(USB)ドライブ1105または任意の他の記憶装置(フロッピー(登録商標)ディスク等)によってユーザのローカルマシン1107に供給されてもよい。
患者の脚および固定器の3D位置は、開示された方法に従って決定される。その後、ユーザは、送信されたX線データを処理して、整形外科固定器1109に与えられなければならない必要な調整を決定する。ユーザは、これらの決定に基づき、ユーザのローカルマシン1107に接続されたプログラマブルレンチ1106を用いて、接続点を締めるまたは緩めることにより、自動的に整形外科固定器1109を調整してもよい。代替的に、ユーザ(最も可能性が高い者は内科医または医療スタッフメンバー)は、送信されたX線データに基づき、整形外科固定器1109を手作業で調整してもよい。
本明細書に記載される特定の実施形態は、発明の限定ではなく、例示という方法によって示されているということが理解される。本発明の主要な特徴は、本発明の範囲から逸脱することなく、様々な実施形態において採用可能である。当業者は、本明細書に記載された特定の手順のいくつかの均等物を、せいぜい日常的な実験法を用いることにより、認識し、または把握することができる。そのような均等物は、本発明の範囲内であると考えられ、特許請求の範囲によってカバーされる。
本明細書において言及したすべての刊行物および特許出願は、本発明に関連する技術分野における当業者の水準を示す。すべての刊行物および特許出願は、各刊行物または各特許出願が参照により組み込まれるものとして個別にかつ独立して示された場合と同程度に、参照により本明細書に組み込まれる。
本発明の方法およびシステムは、好適な実施形態に関して記述してきたが、これらの方法およびシステムに対して、また、本明細書に記載される方法の各ステップまたはステップの系列において、本発明のコンセプト、精神および範囲から逸脱することなく、変更を加えることができるということが当業者には明白である。当業者に明白な、すべてのそのような同様な置き換えおよび修正は、添付の特許請求の範囲によって定義される本発明の、精神、範囲およびコンセプトの範囲内にあると考えられる。

Claims (32)

  1. 身体部分の3Dモデルを作成する方法であって、
    前記身体部分は物体に結合され、
    前記物体は、前記物体に沿って、事前に決定された距離において、複数のマーカーを備え、
    前記方法は、
    第1の向きにおいて、X線源とX線撮像装置との間に配置された前記身体部分および前記物体の、第1のレントゲン画像を受信することであって、前記第1のレントゲン画像は、
    ‐前記身体部分と、
    ‐前記物体と、
    ‐前記複数のマーカーと
    の画像を含む、第1のレントゲン画像を受信することと、
    第2の向きにおいて、前記X線源と前記X線撮像装置との間に配置された前記身体部分および前記物体の、第2のレントゲン画像を受信することであって、前記第2のレントゲン画像は、
    ‐前記身体部分と、
    ‐前記物体と、
    ‐前記複数のマーカーと
    の画像を含む、第2のレントゲン画像を受信することと、
    前記第1のレントゲン画像上の前記複数のマーカーの投影の第1の組を決定することと、
    前記複数のマーカー間の前記事前に決定された距離と、前記第1のレントゲン画像上の前記複数のマーカーの投影の前記第1の組とを用いて、前記第1の向きにおける前記X線撮像装置に対する、前記X線源の第1の3D位置と、前記物体の第1の3D位置とを決定することと、
    前記第2のレントゲン画像上の前記複数のマーカーの投影の第2の組を決定することと、
    前記複数のマーカー間の前記事前に決定された距離と、前記第2のレントゲン画像上の前記複数のマーカーの投影の前記第2の組とを用いて、前記第2の向きにおける前記X線撮像装置に対する、前記X線源の第2の3D位置と、前記物体の第2の3D位置とを決定することと、
    前記第1および第2の向きにおける前記X線撮像装置に対する、前記複数のマーカーの前記各3D位置を用いて、3D参照フレームにおける前記第1および第2の3D物体投影を位置合わせすることと、
    前記第1および第2の3D物体投影に基づいて、前記3D参照フレームにおける撮像された身体部分の3Dモデルを作成することと
    を備える、方法。
  2. 前記複数のマーカーは、複数のストラットが少なくとも1つの固定部材に接続される複数のジョイントを含む、請求項1に記載の方法。
  3. 前記物体は整形外科固定器である、請求項1に記載の方法。
  4. 前記複数のマーカーは5つのマーカーを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    さらに、次の関係のうち1つを満たし、
    Figure 2016513529
    ただし(x,y,z)は前記X線源の座標であり、(x0,y0,z0)〜(x4,y4,z4)は前記5つのマーカーの座標であり、(X0,Y0,Z0)〜(X4,Y4,Z4)は前記第1のレントゲン画像上の前記5つのマーカーの投影の前記第1の組の座標であり、l01,l02,l03,l04,l12,l13,l14,l23,l24,l34は、前記5つのマーカー間の前記事前に決定された距離であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    さらに、次の関係のうち1つを満たし、
    Figure 2016513529
    ただし(‘x,’y,’z)は前記X線源の座標であり、(x0,y0,z0)〜(x4,y4,z4)は前記5つのマーカーの座標であり、(‘X0,’Y0,’Z0)〜(‘X4,’Y4,’Z4)は前記第2のレントゲン画像上の前記5つのマーカーの投影の前記第2の組の座標であり、l01,l02,l03,l04,l12,l13,l14,l23,l24,l34は、前記5つのマーカー間の前記事前に決定された距離である、
    請求項1に記載の方法。
  5. 前記複数のマーカーは4つのマーカーを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x,y,z)は前記X線源の座標であり、(x0,y0,z0)〜(x3,y3,z3)は前記4つのマーカーの座標であり、(X0,Y0,Z0)〜(X3,Y3,Z3)は前記第1のレントゲン画像上の前記4つのマーカーの投影の前記第1の組の座標であり、l01,l02,l03,l04,l12,l13,l14,l23,l24,l34は、前記4つのマーカー間の前記事前に決定された距離であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(‘x,’y,’z)は前記X線源の座標であり、(x0,y0,z0)〜(x3,y3,z3)は前記4つのマーカーの座標であり、(‘X0,’Y0,’Z0)〜(‘X3,’Y3,’Z3)は前記第2のレントゲン画像上の前記4つのマーカーの投影の前記第2の組の座標であり、l01,l02,l03,l04,l12,l13,l14,l23,l24,l34は、前記4つのマーカー間の前記事前に決定された距離である、
    請求項1に記載の方法。
  6. 前記複数のマーカーは4つのマーカーを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x0,y0,z0)〜(x3,y3,z3)は前記4つのマーカーの座標であり、l01,l02,l03,l12,l13,l23は、前記4つのマーカー間の前記事前に決定された距離であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x0,y0,z0)〜(x3,y3,z3)は前記4つのマーカーの座標であり、l01,l02,l03,l12,l13,l23は、前記4つのマーカー間の前記事前に決定された距離である、
    請求項1に記載の方法。
  7. 前記複数のマーカーは4つのマーカーを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x,y,z)は前記X線源の座標であり、(x0,y0,z0)〜(x3,y3,z3)は前記4つのマーカーの座標であり、(X0,Y0,Z0)〜(X3,Y3,Z3)は前記第1のレントゲン画像上の前記4つのマーカーの投影の前記第1の組の座標であり、l01,l02,l03,l12,l13,l23は、前記4つのマーカー間の前記事前に決定された距離であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(‘x,’y,’z)は前記X線源の座標であり、(x0,y0,z0)〜(x3,y3,z3)は前記4つのマーカーの座標であり、(‘X0,’Y0,’Z0)〜(‘X3,’Y3,’Z3)は前記第2のレントゲン画像上の前記4つのマーカーの投影の前記第2の組の座標であり、l01,l02,l03,l12,l13,l23は、前記4つのマーカー間の前記事前に決定された距離である、
    請求項1に記載の方法。
  8. 前記第1のレントゲン画像において、撮像された身体部分の、第1身体部分アウトラインを特定することと、
    前記第2のレントゲン画像において、撮像された身体部分の、第2身体部分アウトラインを特定することと、
    前記第1身体部分アウトラインから、前記X線源の前記第1の3D位置への、第1の3D身体部分投影を準備することと、
    前記第2身体部分アウトラインから、前記X線源の前記第2の3D位置への、第2の3D身体部分投影を準備することと、
    前記第1の身体部分投影および第2の身体部分投影に基づき、前記3D参照フレームにおいて前記撮像された身体部分の3Dモデルを作成することと
    を備える、請求項1に記載の方法。
  9. 前記方法は、
    前記3D参照フレームにおいて傾斜軸を特定することであって、前記傾斜軸は、
    ‐前記第1の向きにおける前記X線源の前記第1の位置に対応する、前記3D参照フレームにおける第1の3D位置と、
    ‐前記第2の向きにおける前記X線源の前記第2の位置に対応する、前記3D参照フレームにおける第2の3D位置と、
    の間を通過する、傾斜軸を特定することと、
    前記傾斜軸を通過し、前記3D参照フレームにおける前記撮像された身体部分の前記第1および第2の3D投影を通過する、1つ以上の交差平面を特定することと、
    前記1つ以上の交差平面のそれぞれについて、次のステップa)〜c)を実行すること:
    a)前記3D参照フレームにおいて、前記第1の3D身体部分投影および前記第2の3D身体部分投影の間の交点を1つ以上と、前記交差平面とを特定すること、
    b)前記交差平面内で各前記交点を接続する1つ以上の多角形を準備すること、
    c)前記1つ以上の多角形のそれぞれの内部に1つ以上の閉曲線を準備することであって、前記1つ以上の閉曲線は、前記交差平面における前記撮像された身体部分の断面図に対応する、1つ以上の閉曲線を準備すること、
    と、
    前記撮像された身体部分の3Dモデルを形成するために、前記3D参照フレームにおいて、前記閉曲線のそれぞれを接続する表面を準備することと
    をさらに備える、請求項8に記載の方法。
  10. 身体部分の3Dモデルを作成する方法であって、
    前記身体部分は物体に結合され、
    前記物体は、事前に決定された長さの複数のストラットを備え、
    前記複数のストラットは、それぞれ、事前に決定された諸寸法の少なくとも2つの固定部材に接続され、
    前記方法は、
    第1の向きにおいて、X線源とX線撮像装置との間に配置された前記身体部分および前記物体の、第1のレントゲン画像を受信することであって、前記第1のレントゲン画像は、
    ‐前記身体部分と、
    ‐前記物体と、
    ‐事前に決定された長さの前記複数のストラットであって、前記複数のストラットは、それぞれ少なくとも2つの固定部材に、2つの接続点で接続され、前記2つの接続点の間の距離は事前に決められている、前記複数のストラットと
    の画像を含む、第1のレントゲン画像を受信することと、
    第2の向きにおいて、前記X線源と前記X線撮像装置との間に配置された前記身体部分および前記物体の、第2のレントゲン画像を受信することであって、前記第2のレントゲン画像は、
    ‐前記身体部分と、
    ‐前記物体と、
    ‐事前に決定された長さの前記複数のストラットであって、前記複数のストラットは、それぞれ少なくとも2つの固定部材に、2つの接続点で接続され、前記2つの接続点の間の距離は事前に決められている、前記複数のストラットと
    の画像を含む、第2のレントゲン画像を受信することと、
    前記第1のレントゲン画像上の、前記複数のストラットの長手軸の投影の第1の組を決定することと、
    前記複数のストラットの前記ストラット接続点間の前記事前に決定された距離と、前記第1のレントゲン画像上の前記複数のストラットの長手軸の投影の前記第1の組とを用いて、前記X線撮像装置に対する、前記X線源の第1の3D位置と、前記物体の第1の3D位置とを決定することと、
    前記第2のレントゲン画像上の、前記複数のストラットの長手軸の投影の第2の組を決定することと、
    前記複数のストラットの前記ストラット接続点間の前記事前に決定された距離と、前記第2のレントゲン画像上の前記複数のストラットの長手軸の投影の前記第2の組とを用いて、前記X線撮像装置に対する、前記X線源の第2の3D位置と、前記物体の第2の3D位置とを決定することと
    前記第1および第2の向きにおける前記X線撮像装置に対する、前記複数のストラットの前記各3D位置を用いて、3D参照フレームにおける前記第1および第2の3D物体投影を位置合わせすることと、
    前記第1および第2の3D物体投影に基づいて、前記3D参照フレームにおける撮像された身体部分の3Dモデルを作成することと
    を備える、方法。
  11. 前記物体は整形外科固定器である、請求項10に記載の方法。
  12. 前記複数のストラットは5つのストラットを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    さらに、次の関係のうち1つを満たし、
    Figure 2016513529
    ただし(x,y,z)は前記X線源の座標であり、(x0,y0,z0)〜(x9,y9,z9)は前記固定部材に対する前記10個のストラット接続点の座標であり、(X0,Y0,Z0)から(X1,Y1,Z1)まで、…、(X8,Y8,Z8)から(X9,Y9,Z9)までは、前記第1のレントゲン画像上の各前記ストラットの長手軸の投影の前記第1の組の座標であり、l01,l02,l03,l04,l06,l07,l08,l09,l12,l13,l14,l15,l16,l17,l18,l19,l23,l24,l25,l26,l27,l28,l29は、前記5つのストラットの前記ストラット接続点間の事前に決定された距離であり、siは未知の比率であり、v(i−1)ix,v(i−1)iy,v(i−1)izは各前記ストラットの長手軸の投影のベクトルであり、ただしv(i−1)ix=Xi−X(i−1),v(i−1)iy=Yi−Y(i−1),v0iz=Zi−Z(i−1)であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    さらに、次の関係のうち1つを満たし、
    Figure 2016513529
    ただし(‘x,’y,’z)は前記X線源の座標であり、(x0,y0,z0)〜(x9,y9,z9)は前記固定部材に対する前記10個のストラット接続点の座標であり、(‘X0,’Y0,’Z0)から(‘X1,’Y1,’Z1)まで、…、(‘X8,’Y8,’Z8)から(‘X9,’Y9,’Z9)までは、前記第2のレントゲン画像上の各前記ストラットの長手軸の投影の前記第2の組の座標であり、l01,l02,l03,l04,l06,l07,l08,l09,l12,l13,l14,l15,l16,l17,l18,l19,l23,l24,l25,l26,l27,l28,l29は、前記5つのストラットの前記ストラット接続点間の事前に決定された距離であり、siは未知の比率であり、‘v(i−1)ix,‘v(i−1)iy,‘v(i−1)izは各前記ストラットの長手軸の投影のベクトルであり、ただし‘v(i−1)ix=’Xi−‘X(i−1),‘v(i−1)iy=’Yi−‘Y(i−1),‘v0iz=’Zi−‘Z(i−1)である、
    請求項10に記載の方法。
  13. 前記複数のストラットは4つのストラットを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x,y,z)は前記X線源の座標であり、(x0,y0,z0)〜(x7,y7,z7)は前記固定部材に対する前記8つのストラット接続点の座標であり、(X0,Y0,Z0)から(X1,Y1,Z1)まで、…、(X6,Y6,Z6)から(X7,Y7,Z7)までは、前記第1のレントゲン画像上の各前記ストラットの長手軸の投影の前記第1の組の座標であり、(X’0,Y’0,Z’0)〜(X’7,Y’7,Z’7)は、前記第1のレントゲン画像上の前記8つのストラット接続点の投影の前記第1の組の座標であり、、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、前記4つのストラットの前記ストラット接続点間の事前に決定された距離であり、siは未知の比率であり、v(i−1)ix,v(i−1)iy,v(i−1)izは各前記ストラットの長手軸の投影のベクトルであり、ただしv(i−1)ix=Xi−X(i−1),v(i−1)iy=Yi−Y(i−1),v0iz=Zi−Z(i−1)であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(‘x,’y,’z)は前記X線源の座標であり、(x0,y0,z0)〜(x7,y7,z7)は前記固定部材に対する前記8つのストラット接続点の座標であり、(‘X0,’Y0,’Z0)から(‘X1,’Y1,’Z1)まで、…、(‘X6,’Y6,’Z6)から(‘X7,’Y7,’Z7)までは、前記第2のレントゲン画像上の各前記ストラットの長手軸の投影の前記第2の組の座標であり、(‘X’0,’Y’0,’Z’0)〜(‘X’7,’Y’7,’Z’7)は、前記第2のレントゲン画像上の前記8つのストラット接続点の投影の前記第2の組の座標であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、前記4つのストラットの前記ストラット接続点間の事前に決定された距離であり、siは未知の比率であり、‘v(i−1)ix,‘v(i−1)iy,‘v(i−1)izは各前記ストラットの長手軸の投影のベクトルであり、ただし‘v(i−1)ix=’Xi−‘X(i−1),‘v(i−1)iy=’Yi−‘Y(i−1),‘v0iz=’Zi−‘Z(i−1)である、
    請求項10に記載の方法。
  14. 前記複数のストラットは4つのストラットを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x0,y0,z0)〜(x7,y7,z7)は前記固定部材に対する前記8つのストラット接続点の座標であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、前記4つのストラットの前記ストラット接続点間の事前に決定された距離であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x0,y0,z0)〜(x7,y7,z7)は前記固定部材に対する前記8つのストラット接続点の座標であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、前記4つのストラットの前記ストラット接続点間の事前に決定された距離である、
    請求項10に記載の方法。
  15. 前記複数のストラットは4つのストラットを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x,y,z)は前記X線源の座標であり、(x0,y0,z0)〜(x7,y7,z7)は前記固定部材に対する前記8つのストラット接続点の座標であり、(X0,Y0,Z0)から(X1,Y1,Z1)まで、…、(X6,Y6,Z6)から(X7,Y7,Z7)までは、前記第1のレントゲン画像上の各前記ストラットの長手軸の投影の前記第1の組の座標であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、前記4つのストラットの前記ストラット接続点間の事前に決定された距離であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(‘x,’y,’z)は前記X線源の座標であり、(x0,y0,z0)〜(x7,y7,z7)は前記固定部材に対する前記8つのストラット接続点の座標であり、(‘X0,’Y0,Z0)から(‘X1,Y1,Z1)まで、…、(‘X6,’Y6,’Z6)から(‘X7,’Y7,’Z7)までは、前記第2のレントゲン画像上の各前記ストラットの長手軸の投影の前記第2の組の座標であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、前記4つのストラットの前記ストラット接続点間の事前に決定された距離である、
    請求項10に記載の方法。
  16. 前記第1のレントゲン画像において、撮像された身体部分の、第1身体部分アウトラインを特定することと、
    前記第2のレントゲン画像において、撮像された身体部分の、第2身体部分アウトラインを特定することと、
    前記第1身体部分アウトラインから、前記X線源の前記第1の3D位置への、第1の3D身体部分投影を準備することと、
    前記第2身体部分アウトラインから、前記X線源の前記第2の3D位置への、第2の3D身体部分投影を準備することと、
    前記第1の身体部分投影および第2の身体部分投影に基づき、前記3D参照フレームにおいて前記撮像された身体部分の3Dモデルを作成することと
    をさらに備える、請求項10に記載の方法。
  17. 前記方法は、
    前記3D参照フレームにおいて傾斜軸を特定することであって、前記傾斜軸は、
    ‐前記第1の向きにおける前記X線源の前記第1の位置に対応する、前記3D参照フレームにおける第1の3D位置と、
    ‐前記第2の向きにおける前記X線源の前記第2の位置に対応する、前記3D参照フレームにおける第2の3D位置と、
    の間を通過する、傾斜軸を特定することと、
    前記傾斜軸を通過し、前記3D参照フレームにおける前記撮像された身体部分の前記第1および第2の3D投影を通過する、1つ以上の交差平面を特定することと、
    前記1つ以上の交差平面のそれぞれについて、次のステップa)〜c)を実行すること:
    a)前記3D参照フレームにおいて、前記第1の3D身体部分投影および前記第2の3D身体部分投影の間の交点を1つ以上と、前記交差平面とを特定すること、
    b)前記交差平面内で各前記交点を接続する1つ以上の多角形を準備すること、
    c)前記1つ以上の多角形のそれぞれの内部に1つ以上の閉曲線を準備することであって、前記1つ以上の閉曲線は、前記交差平面における前記撮像された身体部分の断面図に対応する、1つ以上の閉曲線を準備すること、
    と、
    前記撮像された身体部分の3Dモデルを形成するために、前記3D参照フレームにおいて、前記閉曲線のそれぞれを接続する表面を準備することと
    をさらに備える、請求項16に記載の方法。
  18. 物体の3Dモデルを作成する方法であって、
    前記物体は、前記物体に沿って、事前に決定された距離において、複数のマーカーを備え、
    前記方法は、
    第1の向きにおいて、X線源とX線撮像装置との間に配置された前記物体の、第1のレントゲン画像を受信することであって、前記第1のレントゲン画像は、
    ‐前記物体と、
    ‐前記複数のマーカーと
    の画像を含む、第1のレントゲン画像を受信することと、
    第2の向きにおいて、前記X線源と前記X線撮像装置との間に配置された前記物体の、第2のレントゲン画像を受信することであって、前記第2のレントゲン画像は、
    ‐前記物体と、
    ‐前記複数のマーカーと
    の画像を含む、第2のレントゲン画像を受信することと、
    前記第1のレントゲン画像上の前記複数のマーカーの投影の第1の組を決定することと、
    前記複数のマーカー間の前記事前に決定された距離と、前記第1のレントゲン画像上の前記複数のマーカーの投影の前記第1の組とを用いて、前記第1の向きにおける前記X線撮像装置に対する、前記X線源の第1の3D位置と、前記物体の第1の3D位置とを決定することと、
    前記第2のレントゲン画像上の前記複数のマーカーの投影の第2の組を決定することと、
    前記複数のマーカー間の前記事前に決定された距離と、前記第2のレントゲン画像上の前記複数のマーカーの投影の前記第2の組とを用いて、前記第2の向きにおける前記X線撮像装置に対する、前記X線源の第2の3D位置と、前記物体の第2の3D位置とを決定することと、
    前記第1および第2の向きにおける前記X線撮像装置に対する、前記複数のマーカーの前記各3D位置を用いて、3D参照フレームにおける前記第1および第2の3D物体投影を位置合わせすることと、
    前記第1および第2の3D物体投影に基づいて、前記3D参照フレームにおける撮像された物体の3Dモデルを作成することと
    を備える、方法。
  19. 前記複数のマーカーは、複数のストラットが少なくとも1つの固定部材に接続される複数のジョイントを含む、請求項18に記載の方法。
  20. 前記物体は整形外科固定器である、請求項18に記載の方法。
  21. 前記複数のマーカーは5つのマーカーを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    さらに、次の関係のうち1つを満たし、
    Figure 2016513529
    ただし(x,y,z)は前記X線源の座標であり、(x0,y0,z0)〜(x4,y4,z4)は前記5つのマーカーの座標であり、(X0,Y0,Z0)〜(X4,Y4,Z4)は前記第1のレントゲン画像上の前記5つのマーカーの投影の前記第1の組の座標であり、l01,l02,l03,l04,l12,l13,l14,l23,l24,l34は、前記5つのマーカー間の前記事前に決定された距離であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    さらに、次の関係のうち1つを満たし、
    Figure 2016513529
    ただし(‘x,’y,’z)は前記X線源の座標であり、(x0,y0,z0)〜(x4,y4,z4)は前記5つのマーカーの座標であり、(‘X0,’Y0,’Z0)〜(‘X4,’Y4,’Z4)は前記第2のレントゲン画像上の前記5つのマーカーの投影の前記第2の組の座標であり、l01,l02,l03,l04,l12,l13,l14,l23,l24,l34は、前記5つのマーカー間の前記事前に決定された距離である、
    請求項18に記載の方法。
  22. 前記複数のマーカーは4つのマーカーを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x,y,z)は前記X線源の座標であり、(x0,y0,z0)〜(x3,y3,z3)は前記4つのマーカーの座標であり、(X0,Y0,Z0)〜(X3,Y3,Z3)は前記第1のレントゲン画像上の前記4つのマーカーの投影の前記第1の組の座標であり、l01,l02,l03,l04,l12,l13,l14,l23,l24,l34は、前記4つのマーカー間の前記事前に決定された距離であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(‘x,’y,’z)は前記X線源の座標であり、(x0,y0,z0)〜(x3,y3,z3)は前記4つのマーカーの座標であり、(‘X0,’Y0,’Z0)〜(‘X3,’Y3,’Z3)は前記第2のレントゲン画像上の前記4つのマーカーの投影の前記第2の組の座標であり、l01,l02,l03,l04,l12,l13,l14,l23,l24,l34は、前記4つのマーカー間の前記事前に決定された距離である、
    請求項18に記載の方法。
  23. 前記複数のマーカーは4つのマーカーを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x0,y0,z0)〜(x3,y3,z3)は前記4つのマーカーの座標であり、l01,l02,l03,l12,l13,l23は、前記4つのマーカー間の前記事前に決定された距離であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x0,y0,z0)〜(x3,y3,z3)は前記4つのマーカーの座標であり、l01,l02,l03,l12,l13,l23は、前記4つのマーカー間の前記事前に決定された距離である、
    請求項18に記載の方法。
  24. 前記複数のマーカーは4つのマーカーを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x,y,z)は前記X線源の座標であり、(x0,y0,z0)〜(x3,y3,z3)は前記4つのマーカーの座標であり、(X0,Y0,Z0)〜(X3,Y3,Z3)は前記第1のレントゲン画像上の前記4つのマーカーの投影の前記第1の組の座標であり、l01,l02,l03,l12,l13,l23は、前記4つのマーカー間の前記事前に決定された距離であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(‘x,’y,’z)は前記X線源の座標であり、(x0,y0,z0)〜(x3,y3,z3)は前記4つのマーカーの座標であり、(‘X0,’Y0,’Z0)〜(‘X3,’Y3,’Z3)は前記第2のレントゲン画像上の前記4つのマーカーの投影の前記第2の組の座標であり、l01,l02,l03,l12,l13,l23は、前記4つのマーカー間の前記事前に決定された距離である、
    請求項18に記載の方法。
  25. 物体の3Dモデルを作成する方法であって、
    前記物体は、事前に決定された長さの複数のストラットを備え、
    前記複数のストラットは、それぞれ、事前に決定された諸寸法の少なくとも2つの固定部材に接続され、
    前記方法は、
    第1の向きにおいて、X線源とX線撮像装置との間に配置された前記物体の、第1のレントゲン画像を受信することであって、前記第1のレントゲン画像は、
    ‐前記物体と、
    ‐事前に決定された長さの前記複数のストラットであって、前記複数のストラットは、それぞれ少なくとも2つの固定部材に、2つの接続点で接続され、前記2つの接続点の間の距離は事前に決められている、前記複数のストラットと
    の画像を含む、第1のレントゲン画像を受信することと、
    第2の向きにおいて、前記X線源と前記X線撮像装置との間に配置された前記物体の、第2のレントゲン画像を受信することであって、前記第2のレントゲン画像は、
    ‐前記物体と、
    ‐事前に決定された長さの前記複数のストラットであって、前記複数のストラットは、それぞれ少なくとも2つの固定部材に、2つの接続点で接続され、前記2つの接続点の間の距離は事前に決められている、前記複数のストラットと
    の画像を含む、第2のレントゲン画像を受信することと、
    前記第1のレントゲン画像上の、前記複数のストラットの長手軸の投影の第1の組を決定することと、
    前記複数のストラットの前記ストラット接続点間の前記事前に決定された距離と、前記第1のレントゲン画像上の前記複数のストラットの長手軸の投影の前記第1の組とを用いて、前記X線撮像装置に対する、前記X線源の第1の3D位置と、前記物体の第1の3D位置とを決定することと、
    前記第2のレントゲン画像上の、前記複数のストラットの長手軸の投影の第2の組を決定することと、
    前記複数のストラットの前記ストラット接続点間の前記事前に決定された距離と、前記第2のレントゲン画像上の前記複数のストラットの長手軸の投影の前記第2の組とを用いて、前記X線撮像装置に対する、前記X線源の第2の3D位置と、前記物体の第2の3D位置とを決定することと
    前記第1および第2の向きにおける前記X線撮像装置に対する、前記複数のストラットの前記各3D位置を用いて、3D参照フレームにおける前記第1および第2の3D物体投影を位置合わせすることと、
    前記第1および第2の3D物体投影に基づいて、前記3D参照フレームにおける撮像された物体の3Dモデルを作成することと
    を備える、方法。
  26. 前記物体は整形外科固定器である、請求項25に記載の方法。
  27. 前記複数のストラットは5つのストラットを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    さらに、次の関係のうち1つを満たし、
    Figure 2016513529
    ただし(x,y,z)は前記X線源の座標であり、(x0,y0,z0)〜(x9,y9,z9)は前記固定部材に対する前記10個のストラット接続点の座標であり、(X0,Y0,Z0)から(X1,Y1,Z1)まで、…、(X8,Y8,Z8)から(X9,Y9,Z9)までは、前記第1のレントゲン画像上の各前記ストラットの長手軸の投影の前記第1の組の座標であり、l01,l02,l03,l04,l06,l07,l08,l09,l12,l13,l14,l15,l16,l17,l18,l19,l23,l24,l25,l26,l27,l28,l29は、前記5つのストラットの前記ストラット接続点間の事前に決定された距離であり、siは未知の比率であり、v(i−1)ix,v(i−1)iy,v(i−1)izは各前記ストラットの長手軸の投影のベクトルであり、ただしv(i−1)ix=Xi−X(i−1),v(i−1)iy=Yi−Y(i−1),v0iz=Zi−Z(i−1)であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    さらに、次の関係のうち1つを満たし、
    Figure 2016513529
    ただし(‘x,’y,’z)は前記X線源の座標であり、(x0,y0,z0)〜(x9,y9,z9)は前記固定部材に対する前記10個のストラット接続点の座標であり、(‘X0,’Y0,’Z0)から(‘X1,’Y1,’Z1)まで、…、(‘X8,’Y8,’Z8)から(‘X9,’Y9,’Z9)までは、前記第2のレントゲン画像上の各前記ストラットの長手軸の投影の前記第2の組の座標であり、l01,l02,l03,l04,l06,l07,l08,l09,l12,l13,l14,l15,l16,l17,l18,l19,l23,l24,l25,l26,l27,l28,l29は、前記5つのストラットの前記ストラット接続点間の事前に決定された距離であり、siは未知の比率であり、‘v(i−1)ix,‘v(i−1)iy,‘v(i−1)izは各前記ストラットの長手軸の投影のベクトルであり、ただし‘v(i−1)ix=’Xi−‘X(i−1),‘v(i−1)iy=’Yi−‘Y(i−1),‘v0iz=’Zi−‘Z(i−1)である、
    請求項25に記載の方法。
  28. 前記複数のストラットは4つのストラットを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x,y,z)は前記X線源の座標であり、(x0,y0,z0)〜(x7,y7,z7)は前記固定部材に対する前記8つのストラット接続点の座標であり、(X0,Y0,Z0)から(X1,Y1,Z1)まで、…、(X6,Y6,Z6)から(X7,Y7,Z7)までは、前記第1のレントゲン画像上の前記複数のストラットの長手軸の投影の前記第1の組の座標であり、(X’0,Y’0,Z’0)〜(X’7,Y’7,Z’7)は、前記第1のレントゲン画像上の前記8つのストラット接続点の投影の前記第1の組の座標であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、前記4つのストラットの前記ストラット接続点間の事前に決定された距離であり、siは未知の比率であり、v(i−1)ix,v(i−1)iy,v(i−1)izは各前記ストラットの長手軸の投影のベクトルであり、ただしv(i−1)ix=Xi−X(i−1),v(i−1)iy=Yi−Y(i−1),v0iz=Zi−Z(i−1)であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(‘x,’y,’z)は前記X線源の座標であり、(x0,y0,z0)〜(x7,y7,z7)は前記固定部材に対する前記8つのストラット接続点の座標であり、(‘X0,’Y0,’Z0)から(‘X1,’Y1,’Z1)まで、…、(‘X6,’Y6,’Z6)から(‘X7,’Y7,’Z7)までは、前記第2のレントゲン画像上の各前記ストラットの長手軸の投影の前記第2の組の座標であり、(‘X’0,’Y’0,’Z’0)〜(‘X’7,’Y’7,’Z’7)は、前記第2のレントゲン画像上の前記8つのストラット接続点の投影の前記第2の組の座標であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、前記4つのストラットの前記ストラット接続点間の事前に決定された距離であり、siは未知の比率であり、‘v(i−1)ix,‘v(i−1)iy,‘v(i−1)izは各前記ストラットの長手軸の投影のベクトルであり、ただし‘v(i−1)ix=’Xi−‘X(i−1),‘v(i−1)iy=’Yi−‘Y(i−1),‘v0iz=’Zi−‘Z(i−1)である、
    請求項25に記載の方法。
  29. 前記複数のストラットは4つのストラットを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x0,y0,z0)〜(x7,y7,z7)は前記固定部材に対する前記8つのストラット接続点の座標であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、前記4つのストラットの前記ストラット接続点間の事前に決定された距離であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x0,y0,z0)〜(x7,y7,z7)は前記固定部材に対する前記8つのストラット接続点の座標であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、前記4つのストラットの前記ストラット接続点間の事前に決定された距離である、
    請求項25に記載の方法。
  30. 前記複数のストラットは4つのストラットを含み、
    さらに、前記X線源の前記第1の3D位置および前記物体の前記第1の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(x,y,z)は前記X線源の座標であり、(x0,y0,z0)〜(x7,y7,z7)は前記固定部材に対する前記8つのストラット接続点の座標であり、(X0,Y0,Z0)から(X1,Y1,Z1)まで、…、(X6,Y6,Z6)から(X7,Y7,Z7)までは、前記第1のレントゲン画像上の各前記ストラットの長手軸の投影の前記第1の組の座標であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、前記4つのストラットの前記ストラット接続点間の事前に決定された距離であり、
    さらに、前記X線源の前記第2の3D位置および前記物体の前記第2の3D位置は、次の関係により決定され、
    Figure 2016513529
    ただし(‘x,’y,’z)は前記X線源の座標であり、(x0,y0,z0)〜(x7,y7,z7)は前記固定部材に対する前記8つのストラット接続点の座標であり、(‘X0,’Y0,Z0)から(‘X1,Y1,Z1)まで、…、(‘X6,’Y6,’Z6)から(‘X7,’Y7,’Z7)までは、前記第2のレントゲン画像上の各前記ストラットの長手軸の投影の前記第2の組の座標であり、l01,l02,l03,l04,l06,l07,l12,l13,l14,l15,l16,l17,l23,l24,l25,l26,l27は、前記4つのストラットの前記ストラット接続点間の事前に決定された距離である、
    請求項25に記載の方法。
  31. 命令が符号化されたコンピュータメモリ装置であって、
    前記命令は、コンピュータに、以下のステップ:
    第1の向きにおいて、X線源とX線撮像装置との間に配置された前記身体部分および前記物体の、第1のレントゲン画像を受信することであって、前記第1のレントゲン画像は、
    ‐前記身体部分と、
    ‐前記物体と、
    ‐前記複数のマーカーと
    の画像を含む、第1のレントゲン画像を受信することと、
    第2の向きにおいて、前記X線源と前記X線撮像装置との間に配置された前記身体部分および前記物体の、第2のレントゲン画像を受信することであって、前記第2のレントゲン画像は、
    ‐前記身体部分と、
    ‐前記物体と、
    ‐前記複数のマーカーと
    の画像を含む、第2のレントゲン画像を受信することと、
    前記第1のレントゲン画像上の前記複数のマーカーの投影の第1の組を決定することと、
    前記複数のマーカー間の前記事前に決定された距離と、前記第1のレントゲン画像上の前記複数のマーカーの投影の前記第1の組とを用いて、前記第1の向きにおける前記X線撮像装置に対する、前記X線源の第1の3D位置と、前記物体の第1の3D位置とを決定することと、
    前記第2のレントゲン画像上の前記複数のマーカーの投影の第2の組を決定することと、
    前記複数のマーカー間の前記事前に決定された距離と、前記第2のレントゲン画像上の前記複数のマーカーの投影の前記第2の組とを用いて、前記第2の向きにおける前記X線撮像装置に対する、前記X線源の第2の3D位置と、前記物体の第2の3D位置とを決定することと、
    前記第1および第2の向きにおける前記X線撮像装置に対する、前記複数のマーカーの前記各3D位置を用いて、3D参照フレームにおける前記第1および第2の3D物体投影を位置合わせすることと、
    前記第1および第2の3D物体投影に基づいて、前記3D参照フレームにおける撮像された身体部分の3Dモデルを作成することと
    を実行させる、命令が符号化されたコンピュータメモリ装置。
  32. コンピュータコードセグメントであって、
    前記コンピュータコードセグメントは、コンピュータに、以下のステップ:
    第1の向きにおいて、X線源とX線撮像装置との間に配置された前記身体部分および前記物体の、第1のレントゲン画像を受信することであって、前記第1のレントゲン画像は、
    ‐前記身体部分と、
    ‐前記物体と、
    ‐前記複数のマーカーと
    の画像を含む、第1のレントゲン画像を受信することと、
    第2の向きにおいて、前記X線源と前記X線撮像装置との間に配置された前記身体部分および前記物体の、第2のレントゲン画像を受信することであって、前記第2のレントゲン画像は、
    ‐前記身体部分と、
    ‐前記物体と、
    ‐前記複数のマーカーと
    の画像を含む、第2のレントゲン画像を受信することと、
    前記第1のレントゲン画像上の前記複数のマーカーの投影の第1の組を決定することと、
    前記複数のマーカー間の前記事前に決定された距離と、前記第1のレントゲン画像上の前記複数のマーカーの投影の前記第1の組とを用いて、前記第1の向きにおける前記X線撮像装置に対する、前記X線源の第1の3D位置と、前記物体の第1の3D位置とを決定することと、
    前記第2のレントゲン画像上の前記複数のマーカーの投影の第2の組を決定することと、
    前記複数のマーカー間の前記事前に決定された距離と、前記第2のレントゲン画像上の前記複数のマーカーの投影の前記第2の組とを用いて、前記第2の向きにおける前記X線撮像装置に対する、前記X線源の第2の3D位置と、前記物体の第2の3D位置とを決定することと、
    前記第1および第2の向きにおける前記X線撮像装置に対する、前記複数のマーカーの前記各3D位置を用いて、3D参照フレームにおける前記第1および第2の3D物体投影を位置合わせすることと、
    前記第1および第2の3D物体投影に基づいて、前記3D参照フレームにおける撮像された身体部分の3Dモデルを作成することと
    を実行させるよう動作可能である、コンピュータコードセグメント。
JP2015562961A 2013-03-15 2013-03-15 マーカーまたはストラットの投影を用いて物体の位置を決定する方法 Active JP6209624B2 (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/RU2013/000203 WO2014142703A1 (en) 2013-03-15 2013-03-15 Method of determining the position of an object using projections of markers or struts

Publications (2)

Publication Number Publication Date
JP2016513529A true JP2016513529A (ja) 2016-05-16
JP6209624B2 JP6209624B2 (ja) 2017-10-04

Family

ID=49304291

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015562961A Active JP6209624B2 (ja) 2013-03-15 2013-03-15 マーカーまたはストラットの投影を用いて物体の位置を決定する方法

Country Status (8)

Country Link
US (1) US9959683B2 (ja)
EP (1) EP2973410B1 (ja)
JP (1) JP6209624B2 (ja)
CN (1) CN105051786B (ja)
AU (1) AU2013382253B2 (ja)
BR (2) BR112015022014A2 (ja)
IL (1) IL241476B (ja)
WO (1) WO2014142703A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2019509790A (ja) * 2016-07-14 2019-04-11 エーエムディーティー ホールディングス インコーポレイテッドAmdt Holdings, Inc. 外部骨固定システム

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9204937B2 (en) * 2013-02-19 2015-12-08 Stryker Trauma Gmbh Software for use with deformity correction
US10251705B2 (en) 2016-06-02 2019-04-09 Stryker European Holdings I, Llc Software for use with deformity correction
JP6765092B2 (ja) * 2016-06-20 2020-10-07 国立大学法人千葉大学 位置校正用プレート及びこれを用いた撮像装置並びに位置校正方法
US10835318B2 (en) 2016-08-25 2020-11-17 DePuy Synthes Products, Inc. Orthopedic fixation control and manipulation
US11638623B2 (en) 2016-09-26 2023-05-02 Texas Scottish Rite Hospital For Children Radiography aid for an external fixator
US10874433B2 (en) 2017-01-30 2020-12-29 Stryker European Holdings I, Llc Strut attachments for external fixation frame
WO2019040829A1 (en) * 2017-08-24 2019-02-28 Amdt Holdings, Inc. METHODS AND SYSTEMS FOR DETERMINING ADJUSTMENT REQUIREMENTS OF EXTERNAL FIXING DEVICES
US11457953B2 (en) 2017-08-31 2022-10-04 J & A Medical Llc External fixation alignment gauge
EP3457353B1 (en) * 2017-09-18 2020-11-25 Siemens Healthcare GmbH Method and system for obtaining a true shape of objects in a medical image
USD888947S1 (en) 2018-03-14 2020-06-30 J&A Medical LLC External fixation alignment gauge
CN108937987B (zh) * 2018-05-22 2021-07-02 上海联影医疗科技股份有限公司 一种确定模体中标记物位置的方法和系统
US10743918B2 (en) * 2018-06-27 2020-08-18 Texas Scottish Rite Hospital For Children External fixation connection rod with female attachment
AU2019309393A1 (en) * 2018-07-24 2021-03-11 Amdt Holdings, Inc. Methods and systems of registering a radiographic image and a 3D model of an external fixation device
CN109009376B (zh) * 2018-08-10 2019-12-17 天津大学 骨折并联外固定支架空间位姿的自动识别方法
US11334997B2 (en) 2020-04-03 2022-05-17 Synthes Gmbh Hinge detection for orthopedic fixation
CN112120772B (zh) * 2020-08-24 2023-11-24 北京工业大学 一种用于胫骨畸形矫正的定角度旋转中心并联外固定器
CN112129789A (zh) * 2020-09-22 2020-12-25 云南电网有限责任公司电力科学研究院 一种用于导线缺陷定位的方法及装置
EP4163821A1 (en) * 2021-10-08 2023-04-12 Vito NV A method and system for performing characterization of one or more materials

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003530177A (ja) * 2000-04-07 2003-10-14 カーネギー メロン ユニヴァーシティ コンピューター支援による骨延長
WO2011146703A1 (en) * 2010-05-19 2011-11-24 Synthes Usa, Llc Orthopedic fixation with imagery analysis

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
AU9088801A (en) * 2000-09-14 2002-03-26 Univ Leland Stanford Junior Assessing the condition of a joint and devising treatment
BR0116855B1 (pt) * 2001-02-07 2012-06-12 processo para estabelecer uma representaÇço tridimensional virtual de um osso ou fragmento de osso a partir de imagens de raios-x.
AU2003297280B2 (en) * 2002-11-14 2010-04-08 Extraortho, Inc Method for using a fixator device
EP1639548B1 (en) * 2003-06-24 2008-10-08 Philips Intellectual Property & Standards GmbH Device to generate a three-dimensional image of a moved object
US8265770B2 (en) * 2003-10-02 2012-09-11 Medtronic, Inc. Driver circuitry switchable between energy transfer and telemetry for an implantable medical device
FR2882245B1 (fr) * 2005-02-21 2007-05-18 Gen Electric Procede de determination du deplacement 3d d'un patient positionne sur une table d'un dispositif d'imagerie
US7756244B2 (en) * 2006-06-22 2010-07-13 Varian Medical Systems, Inc. Systems and methods for determining object position
EP2245986B1 (de) * 2008-08-22 2013-10-16 BrainLAB AG Pyramidenförmige Röntgenmarkervorrichtung
WO2010104567A1 (en) * 2009-03-10 2010-09-16 Stryker Trauma Sa External fixation system
US8503745B2 (en) * 2009-05-13 2013-08-06 Medtronic Navigation, Inc. System and method for automatic registration between an image and a subject
DE102009049818A1 (de) * 2009-10-19 2011-04-21 Siemens Aktiengesellschaft Verfahren zur Ermittlung der Projektionsgeometrie einer Röntgenanlage

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003530177A (ja) * 2000-04-07 2003-10-14 カーネギー メロン ユニヴァーシティ コンピューター支援による骨延長
WO2011146703A1 (en) * 2010-05-19 2011-11-24 Synthes Usa, Llc Orthopedic fixation with imagery analysis
JP2013526377A (ja) * 2010-05-19 2013-06-24 ジンテス ゲゼルシャフト ミット ベシュレンクテル ハフツング 画像解析を伴う整形外科用の固定

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2019509790A (ja) * 2016-07-14 2019-04-11 エーエムディーティー ホールディングス インコーポレイテッドAmdt Holdings, Inc. 外部骨固定システム
JP2022031641A (ja) * 2016-07-14 2022-02-22 エーエムディーティー ホールディングス インコーポレイテッド 外部骨固定システム
JP7390347B2 (ja) 2016-07-14 2023-12-01 エーエムディーティー ホールディングス インコーポレイテッド 外部骨固定システム

Also Published As

Publication number Publication date
IL241476A0 (en) 2015-11-30
AU2013382253B2 (en) 2019-11-21
BR112015023127B1 (pt) 2022-02-08
EP2973410A1 (en) 2016-01-20
WO2014142703A1 (en) 2014-09-18
BR112015022014A2 (pt) 2017-08-22
CN105051786B (zh) 2020-03-17
AU2013382253A1 (en) 2015-11-05
BR112015023127A2 (ja) 2017-08-22
US9959683B2 (en) 2018-05-01
EP2973410B1 (en) 2024-04-24
JP6209624B2 (ja) 2017-10-04
US20160042571A1 (en) 2016-02-11
IL241476B (en) 2020-07-30
CN105051786A (zh) 2015-11-11

Similar Documents

Publication Publication Date Title
JP6209624B2 (ja) マーカーまたはストラットの投影を用いて物体の位置を決定する方法
EP2606467B1 (en) Method and system for roentgenography-based modeling
JP7269222B2 (ja) 外部固定装置の調整処方を決定するための方法とシステム
Schulze et al. Precision and accuracy of measurements in digital panoramic radiography
Xiaojun et al. Image guided oral implantology and its application in the placement of zygoma implants
US10705677B2 (en) Method and system for planning implant component position
Brown et al. Computer-assisted location of reference points in three dimensions for radiographic cephalometry
US20210007806A1 (en) A method for obtaining 3-d deformity correction for bones
KR20160004862A (ko) 환자맞춤형 치아교정 모의시술과 이를 통한 시뮬레이션 및 치아 교정장치 또는 치아교정 시술유도장치 제작방법
CN110891488A (zh) 矢状旋转确定
CN110462681A (zh) 用于物理到图像/图像到物理配准和图像验证的多个表面
Andre et al. Effect of radiographic landmark identification errors on the accuracy of three-dimensional reconstruction of the human spine
Novosad et al. Three-dimensional (3-D) reconstruction of the spine from a single X-ray image and prior vertebra models
Lemieux et al. Precision and accuracy of suggested maxillary and mandibular landmarks with cone-beam computed tomography for regional superimpositions: An in vitro study
Lagacé et al. 3D reconstruction of the scapula from biplanar radiographs
Enciso et al. Precision, repeatability and validation of indirect 3D anthropometric measurements with light-based imaging techniques
Elad et al. Three-dimensional measurement of biological surfaces
Zeng et al. Low‐dose three‐dimensional reconstruction of the femur with unit free‐form deformation
Dong et al. Theories and algorithms for 3-D root canal model construction
TWI718578B (zh) 植入物表面映射與映射展開方法
JP3667727B2 (ja) 物体の位置測定方法及びその装置
Tanbour et al. A four-wall virtual reality visualization of patient-specific anatomy: Creating full user immersive experience from computed tomography scans
CN115222885A (zh) 双平面下流程自适应的骨骼结构空间重建方法及装置
Aoki et al. 3D head model construction of individuals utilizing standard model and photogrammetry
Maruyama et al. CT image-based training system for positioning in radiography

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20160205

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20160302

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20161227

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20170110

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20170331

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20170815

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20170911

R150 Certificate of patent or registration of utility model

Ref document number: 6209624

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313117

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250