【誤訳訂正書】
【提出日】平成20年5月12日(2008.5.12)
【誤訳訂正１】
【訂正対象書類名】明細書
【訂正対象項目名】全文
【訂正方法】変更
【訂正の内容】
【書類名】      明細書
【発明の名称】  画像データセットにおいてボリュームの密度を決める方法及びシステム
【特許請求の範囲】
【請求項１】
検査すべき目的物の画像データセットにおいてボリュームの密度を決める方法であって、前記ボリュームは、当該ボリュームが固有の特徴パラメタを持つ材料を含むことにより周囲物と区別される、方法であって、
それぞれ異なる発生位置において、前記ボリューム及び前記周囲物照射用のＸ線を発生させるステップであって、前記各発生位置において、第１のエネルギー及び第２のエネルギーのＸ線を発生させる、ステップ、
二次元のＸ線検出器において、前記第１のエネルギー及び前記第２のエネルギーのＸ線に対応する前記目的物の透過画像の第１の組を得るステップ、
二重エネルギー法に従って第１の較正を実行し、第１のパラメタの組と第２のパラメタの組とを得るステップ、
二重エネルギー法に従って第２の較正を実行し、前記周囲物に特徴的な第３のパラメタと前記材料に特徴的な第４のパラメタとを得るステップ、
前記第１のパラメタの組と前記第２のパラメタの組とを用いて、前記透過画像の第１の組に対し第１の演算を行い、密度画像の組を得るステップ、
前記密度画像の組に前記第３のパラメタを適用し、前記密度画像の組から前記周囲物を削除するステップ、
前記密度画像の組の情報を用いて３Ｄ再構成アルゴリズムを実行し、前記ボリュームの画像を得るステップ、及び
前記ボリュームの画像に前記第４のパラメタを適用し、前記ボリュームにおける前記材料の密度を計算するステップ
を含む方法。
【請求項２】
請求項１に記載の方法において、透過画像の第２の組を得るために、前記密度画像に対し演算を行い、フェルドカンプアルゴリズムが前記３Ｄ再構成アルゴリズムとして用いられる方法。
【請求項３】
請求項１に記載の方法において、前記第３のパラメタは、軟組織の削除角度であり、前記第４のパラメタは、骨密度因子である方法。
【請求項４】
請求項３に記載の方法において、使用すべき前記骨密度因子の値は、ルックアップテーブルで調べられる方法。
【請求項５】
請求項３に記載の方法において、関心のあるボリュームは、前記ボリュームの画像において規定され、前記関心のあるボリュームは、海綿骨のボリュームの少なくとも一部に対応する方法。
【請求項６】
請求項１に記載の方法において、前記密度画像の組のピクセル値は、ダイナミックレンジを高めるために、因子Ｍで乗算される方法。
【請求項７】
請求項１に記載の方法において、第１及び第２の較正材料は、前記二重エネルギー法に従って前記第１の較正を実行するために使用され、前記第１の較正材料の実効原子番号は、１０よりも低いのに対し、前記第２の較正材料の実効原子番号は、１０よりも大きい方法。
【請求項８】
請求項１に記載の方法において、決められるべき前記密度は、ミネラル密度であり、前記第１のエネルギーのＸ線を生じさせるためにＸ線管に必要な電圧は、５０〜８０ｋＶの範囲にあるのに対し、前記第２のエネルギーのＸ線を生じさせるためにＸ線管に必要な電圧は、１００から１５０ｋＶの範囲にある方法。
【請求項９】
請求項８に記載の方法において、前記第１のエネルギー及び前記第２のエネルギーに対応するＸ線を生じさせる前記Ｘ線管の電圧差は、４０ｋＶから５０ｋＶの間にある方法。
【請求項１０】
請求項１乃至９の何れか一項に記載の方法において、Ｘ線源は前記目的物の周りを約１８０°の弧に沿って回転され、各エネルギーに対し、約５０の透過画像が前記弧において異なる角度で取得される方法。
【請求項１１】
請求項１乃至１０の何れか一項に記載の方法において、低いエネルギー照射に対する放射Ｘ線量は、対応する高いエネルギー照射に対する放射線量と同じ大きさの次数である方法。
【請求項１２】
検査すべき目的物の画像データセットにおいてボリュームの密度を決めるシステムであって、前記ボリュームは、当該ボリュームが固有の特徴パラメタを持つ材料を含むことにより周囲物と区別される、システムであって、
それぞれ異なる発生位置において、前記ボリューム及び前記周囲物照射用のＸ線を発生させる手段であって、前記各発生位置において、第１のエネルギー及び第２のエネルギーのＸ線を発生させる手段と、
前記第１のエネルギー及び前記第２のエネルギーのＸ線に対応する前記目的物の透過画像の第１の組を得るための二次元のＸ線検出器と、
二重エネルギー法に従って第１の較正を実行し、第１のパラメタの組と第２のパラメタの組とを得る手段と、
二重エネルギー法に従って第２の較正を実行し、前記周囲物に特徴的な第３のパラメタと前記材料に特徴的な第４のパラメタとを得る手段と、
前記第１のパラメタの組と前記第２のパラメタの組とを用いて前記透過画像の第１の組に対し第１の演算を行い、密度画像の組を得る手段と、
前記密度画像の組に前記第３のパラメタを適用し、前記密度画像の組から前記周囲物を削除する手段と、
前記密度画像の組の情報を用いて３Ｄ再構成アルゴリズムを実行し、前記ボリュームの画像を得る手段と、
前記ボリュームの画像に前記第４のパラメタを適用し、前記ボリュームにおける前記材料の密度を計算する手段と、
を有するシステム。
【請求項１３】
請求項１２に記載のシステムにおいて、
前記Ｘ線を発生させるＸ線源及び前記Ｘ線のエネルギーを変化させる手段を備えるＸ線装置と、
前記第１の較正及び前記第２の較正を行う較正手段と、
を有するシステム。
【請求項１４】
請求項１３に記載のシステムにおいて、前記較正手段は、適当な第１の較正ファントムと適当な第２の較正ファントムとを含み、さらに前記パラメタを計算する他の算術手段を含むシステム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、検査すべき目的物の画像データセットにおいて、ボリュームの密度を決める方法に関する。前記ボリュームは、このボリュームが固有の特徴パラメタ（μ）を持つ材料を含むことにより周囲物と区別することができ、前記方法は、異なる照射位置において、前記ボリューム及び前記周囲物をＸ線に曝すステップを含み、各照射位置において、Ｘ線の第１のエネルギー及び第２のエネルギーで照射が行われ、Ｘ線検出器においてこれらのエネルギーに対応する前記目的物の透過画像の第１の組を得る。
【０００２】
この種類の方法は、海綿骨(trabecular bone)の骨密度が決められるべき骨密度計測の分野に用いられる。この場合、撮影すべきボリュームは、それの特徴的な線形吸収係数（μ）に基づいて前記周囲物（軟組織、皮質骨）と区別され得る海綿骨を含んでいる。この海綿骨の骨密度の値は、例えば骨粗しょう症を診察するような患者の他の診察にも用いることが可能である。従って、骨密度の値は、撮影すべき対象区域に存在する他の組織による誤りが無く、定量的に決められ得ることが重要である。この種類の組織は、例えば撮影すべき対象区域に存在する皮質骨及び血管の起こり得る石灰化である。
【０００３】
【従来の技術】
骨密度計測を実行する方法は、米国特許番号第US5,778,045号から既知である。この既知の方法に従って、検査すべき患者の透過画像は、Ｘ線を生じさせるＸ線源を含むＸ線装置によって形成される。このＸ線源は、目的物の周りを回転可能であり、生じたＸ線は、細いビームを形成するように平行にされる。この既知の装置は、異なる２つのエネルギーでＸ線を生じさせる手段、すなわち、いわゆる二重エネルギー法を実行する手段を備える。この場合、各照射位置において２つのエネルギーで照射が行われる。既知の方法に従い、次に、両方のエネルギーに対する照射ボリューム内の全部の骨の吸収度に関連する分量を得るように演算が行われる。最終画像に対する軟組織の寄与を制限するために、この計算された分量は、対象区域内の画像平面における距離の関数として図示される。この軟組織は、しきい値を計算された分量に適用することによりグラフィック的に削除される。
【０００４】
【発明が解決しようとする課題】
既知の方法の欠点は、軟組織に対応するボリュームが正確に削除されないことである。確かに、計算された分量は、軟組織の吸収に対応する特徴的な作用、すなわち低く比較的に平坦なラインと、骨の吸収に対応する僅かに高い区域とを示す。脊椎に対しては、例えば前記計算された分量の作用の結果としてのグラフィック的表示は、脊椎に対応するほぼ高い値の区域を示し、この脊椎の左側及び右側に対しては、軟組織に対応する２つの低い値の区域を示す。前記計算された分量の大きさの絶対値は、脊椎の両側にある軟組織に対応する区域において同じ値のものであることはありそうにもない。しきい値に対し、１つの値しか選択されないとき、軟組織が画像から完全には削除されないか、又は骨のボリュームの一部が失われるかが必ず起こる。既知の方法の他の欠点は、患者の断面にのみ照射が行われることである。ボリューム情報を取得するためには、Ｘ線源の移動を繰り返すことが必要であり、この手続きを時間のかかる操作にする。既知の方法の他の欠点は、骨密度に関するデータが二次元でしか示されないという事実にある。
【０００５】
【課題を解決するための手段】
本発明の目的は、分析すべき画像において、軟組織の寄与が正確に且つほぼ完全に削除されるような方法を提供することである。本発明の他の目的は、海綿骨に対応するボリュームの画像を提供することであり、この海綿骨の骨密度の絶対値は、前記画像内に示される。この目的のために、本発明は、前記方法が
−Ｘ線検出器として二次元の検出器を用いるステップ、
−二重エネルギー法に従って第１の較正を実行し、パラメタ(xpl)の第１の組とパラメタ(xal)の第２の組とを得るステップ、
−二重エネルギー法に従って第２の較正を実行し、周囲物に特徴的な第３のパラメタ（φ）と材料に特徴的な第４のパラメタ（ｆ）とを得るステップ、
−前記パラメタ(xpl)の第１の組と前記パラメタ(xal)の第２の組とを用いながら、透過画像の第１の組に対し第１の演算を行い、密度画像の組を得るステップ、
−前記密度画像の組に前記第３のパラメタ（φ）を適用し、前記密度画像の組から前記周囲物を削除するステップ、
−前記密度画像の組の情報を用いて３Ｄ再構成アルゴリズムを実行し、前記ボリュームの画像を得るステップ、及び
−前記ボリュームの画像に第４のパラメタ（ｆ）を適用し、前記ボリューム内の前記材料密度を計算するステップ、
も含むことを特徴とする。
【０００６】
二次元のＸ線検出器を利用しているので、実際には、Ｘ線源が目的物の周りを一回移動するだけで十分であり、この移動中、Ｘ線の第１及び第２のエネルギーを使用しながら、弧に沿った多数の位置において照射が行われる。しかしながら、既知の方法は、細いビーム形状と対応するアレイ検出器とを利用している。このような断層撮影画像を作成すること自体は既知であり、当業者には明白なことである。最初にＸ線源が患者の周りを移動し、その間、第１のエネルギーだけで照射を行うことも可能である。この場合、補足的な第２の断層撮影の照射が第２のＸ線エネルギーで行わなければならない。この場合、連続して患者に照射する、すなわち弧上の対応する位置において第１及び第２のエネルギーで照射を行うことも可能である。二次元のＸ線検出器を使用する他の利点は、再構成すべきボリューム画像データにおいて、立方体のボクセルからなる結果としての等方性分布となることである。これは、３Ｄ再構成の正確さを高め、調査すべき目的物に適用されなければならない最小の照射線量も減少させる。適切な検出器の特徴的なフォーマットは、大体１５×１５ｃｍ2である。本発明による方法は、二重エネルギー法に従って較正を実行することも含んでいる。この種類の較正自体は、H. Neale Cardinal 及び A. Fenster著 "An accurate method for direct dual-energy calibration and decomposition", Med. Phys. 17(3), 1990, p327から既知である。この方法は、材料のＸ線吸収が２つの基本材料の線形結合により表されるという洞察に基づいている。本発明に従って、二重エネルギー法を用いた較正パラメタ(xpl)の第１の組と較正パラメタ(xal)の第２の組との取得は、この目的のために特に設計され、２つの基本材料から成る較正ファントムを利用しながら行われる。本発明の方法に従う第２の較正は、検査すべき材料の周囲物の構成を実質的に表す第２の較正ファントムを用いて行われる。二重エネルギー法に従い、第２の較正ファントムの透過画像を形成した後、これら画像はさらに、前記周囲物を２つの基本材料の線形結合に分解することを達成するように処理される。このような分解演算の状況において、前記透過画像は、密度画像に変換される。前記周囲物は、２つの基本材料の直交空間におけるベクトルとして表される。このベクトルはこのとき、水平軸（例えば第１の基本材料）に対して与えられる角度（α）を囲む。このベクトルの絶対長は、前記周囲物の材料の厚さに及び密度に比例する。例えば、第１の基本材料から成る軸へのこのベクトルの投影は、二重エネルギー分解に従う周囲物における第１の基本材料の寄与に対応する。この周囲物を削除するために、このベクトルの垂直ライン上へ投影することが単に必要である。この角度は、削除角φ＝９０＋αと規定され、本発明に従う方法における第３のパラメタとしても用いられる。これは、目的物の二重エネルギー透過画像から得られるベクトルが２つの基本材料の直交空間内に構築されなければならず、前記目的物の画像から前記周囲物の寄与を削除するために、続いて対応する削除角度に投影されなければならないことを意図している。結果としてのベクトルは、材料の寄与を排他的に含んでおり、すなわち前記周囲物を含まない。この投影が行われるとき、材料ベクトルの長さの絶���値は、もはやこの材料の密度を表さない。前記材料の密度の絶対値を決めることを可能にするため、既知の密度からなる材料のファントムを使用しながら、本発明による方法を用いて他の較正が実行される。このファントムの透過画像も二重エネルギー法により取得され、その後、このファントムに特徴的なベクトルが２つの基本材料の空間内に構成され、前記削除角度に投影される。前記ファントムにおける材料の既知の密度と、結果としてのベクトルの長さの絶対値との比は、後段において目的物データに適用されるべき第４の因子（ｆ）を生じさせる。先に記載された二重エネルギー法に従う目的物の密度画像から前記周囲物を削除した後、３Ｄ再構成が結果としての情報に基づいて実行され、前記ボリュームの画像を得る。前記ボリュームの骨密度の絶対値を計算するために、本発明の方法に従って、第４の因子（ｆ）が前記ボリュームの画像に適用される。初めに、関心のある撮影されたボリュームにおける周囲物の寄与は、二重エネルギー法に従う正確な較正の適用により削除され、次に、ユーザは、材料の密度分布とみなす三次元情報を与えられる。本発明による方法を用いて確かな結果を達成するために、前記較正及び前記目的物の透過画像の形成に対しても適切に安定したＸ線エネルギーを用いることが必要であることに注意されたい。Ｘ線エネルギーの安定性の正確な制御自体は既知である。較正自体は、目的物の透過画像を形成する前、形成中又は形成した後に行うことが可能である。前記較正が送信照射中に実行される場合、較正ファントムは、撮影すべきボリュームに隣接して位置決めされるように、Ｘ線装置のテーブルトップに導入される。２つの異なるＸ線のエネルギーを持つＸ線源は、それ自体が既知である方式、例えばＸ線源と目的物との間にフィルタを用いることで実現可能である。
【０００７】
本発明による方法のバージョンは、密度画像に対し演算を行い、透過画像の第２の組を得て、フェルドカンプアルゴリズム(Feldkamp algorithm)が３Ｄ再構成アルゴリズムとして用いられることを特徴とする。このフェルドカンプアルゴリズムが前記ボリュームの画像の３Ｄ再構成を得るのに適用される場合、入力データとして透過情報を供給することが必要である。透過画像の第１の組の二重エネルギー分解が密度画像の組を形成することになるので、透過画像の第２の組を計算するために単一の演算が行われることが望ましい。つまり、このセットはさらに、フェルドカンプアルゴリズム用の入力データとして用いられる。密度情報及び透過情報は、互いに指数的に関係している。
【０００８】
本発明による方法の他のバージョンは、第３のパラメタ（φ）が軟組織の削除角度であり、第４のパラメタ（ｆ）は、骨密度因子であることを特徴とする。本発明による方法が骨密度の分野に使用される場合、軟組織ファントムは、第２の寄与し得るファントムとして使用され、このファントムは、軟組織の構成を実質的に表している。軟組織の削除角度は、第３のパラメタ（φ）として使用される。この軟組織の削除角を用いた結果として、この軟組織の寄与が患者の透過画像から削除される。既に説明されたように、結果としての因子の長さは、前記材料の１つに対する前記透過画像からの削除計算が実行されるとき、残存する材料の絶対密度をもはや表さない。残存する骨の密度を計算するために、既知のＣａ含有量のファントムは、第４の因子（ｆ）を計算するための較正ファントムとして用いられる。この場合、第４の因子（ｆ）は、骨密度因子に等しい。
【０００９】
本発明による方法の他のバージョンは、使用すべき骨密度因子の値がルックアップテーブルで調べられることを特徴とする。骨密度因子が先験的決められるとき、この因子の値は、コンピュータのメモリに記憶されたルックアップテーブルに記憶される。このテーブルは、使用されるファントムの既知の全Ｃａ含有量のＣａに関する対応する線形吸収係数との比を含む。これら分量の間に線形レベル以上のもの(supralinearity)が起こる場合、Ｃａファントムに関し２つ以上の較正測定を実行することが得策である。
【００１０】
本発明による方法の他のバージョンは、関心のあるボリュームは、このボリュームの画像において規定され、この関心のあるボリュームは、海綿骨のボリュームの少なくとも一部に対応していることを特徴とする。既に説明されたように、正確な較正及び分解が二重エネルギー法を介して実行されるとき、ボリュームの結果としての画像は、石灰化物質だけを含んでいる。例えば、近くにある血管が石灰化した狭窄を持つ区域を含むとき、このような区域は、結果としての画像にも撮影されるだろう。これらの区域が正確な軟組織の削除のために骨のボリュームには関連しないので、これら区域は、簡単に無視されるか、再構成アルゴリズムのアプリケーションによってボリュームデータから切り取られる。前記骨に対応するボリュームは、皮質骨及び海綿骨を一緒に有する。骨粗しょう症の場合、海綿骨にのみ脱石灰化が起こるので、トラッキング分析が海綿骨に対応するボリュームのみを選択することが有利である。このサブ選択、すなわち関心のあるボリュームの選択は、手動又は既知のアルゴリズムで実行される。
【００１１】
本発明による方法の他のバージョンは、密度画像の組のピクセル値がダイナミックレンジを高めるために、因子（Ｍ）で乗算されることを特徴とする。密度画像の組のピクセル値に対する実験上に決められたスケール因子の応用は、高められたダイナミック画像を生じさせることが可能であることが分かる。例えば、０．７５のスケール因子の応用は、２００mg/mlのＣａ含有量を持つファントムにとって魅力的な結果を生じさせる。
【００１２】
本発明による方法の他のバージョンは、第１及び第２の較正材料が、二重エネルギー法に従って第１の較正を実行するために使用され、第１の較正材料の実効原子番号は、１０よりも低いのに対し、第２の較正材料の実効原子番号は、１０よりも高いことを特徴とする。本発明による方法が骨密度の分野に使用される場合、第１の較正材料の原子番号は１０よりも低いのに対し、第２の較正材料の原子番号は１０よりも高い２つの較正材料を使用することは、魅力的な分解結果を提供する。例えばポリメタクリル酸メチル(polymethyl methacrylate)のような合成材料は、第１の較正材料の実施例であり、アルミニウムは第２の較正材料の実施例である。アルミニウムは、骨のＸ線吸収特性に十分対応するＸ線吸収特性を持つので使用され、合成材料は、軟組織を表す。明らかに、他の較正材料を使用することも可能である。その上、較正目的にとって最適な第１及び第２の構成材料の実効原子番号の差は、約６に達することが分かっている。
【００１３】
本発明の他のバージョンは、決められるべき密度が骨密度であり、第１のエネルギーから成るＸ線を生じさせるためにＸ線管に必要とされる電圧は、５０〜８０ｋＶの範囲にあるのに対し、第２のエネルギーから成るＸ線を生じさせるためにＸ線管に必要とされる電圧は、１００〜１５０ｋＶの範囲にあることを特徴とする。骨密度計測の分野における応用にとって、Ｘ線管の陽極に用いられる最適な電圧は、前記第１のエネルギーに対しては５０〜８０ｋＶの範囲、前記第２のエネルギーに対しては１００〜１５０ｋＶの範囲にあるべきである。このキロボルト（ｋＶ）値は、発生する制動輻射(bremsstrahlung)の連続における最大エネルギーを決め、最小値は、Ｘ線が曝される場所をフィルタリングすることにより決められる。
【００１４】
本発明による方法の他のバージョンは、第１及び第２のエネルギーに対応するＸ線を生じさせるＸ線管電圧間の差が４０ｋＶと５０ｋＶとの間となるように選択されることを特徴とする。このステップは、第１及び第２のエネルギーの上記差が特に骨密度計測に適用されるような二重エネルギー法に従う材料の分解に適しているという事実の認識に基づいている。
【００１５】
本発明による方法の他のバージョンは、Ｘ線源が前記目的物の周りの約１８０°の弧に沿って回転可能であり、各エネルギーに対し、約５０の透過画像が上記回転中に前記弧上の異なる角度で取得されることを特徴とする。このステップは、この透過画像の数は、前記ボリュームの画像を得るために正確な３Ｄ再構成を実行するのに十分な情報を提供するという事実の認識に基づいている。
【００１６】
本発明による方法の他のバージョンは、低いエネルギー照射に対し放射される放射線量は、対応する高いエネルギー照射に対し放射される放射線量と同じ大きさの次数であることを特徴とする。このステップは、低いＸ線エネルギーと高いＸ線エネルギーとを備える透過画像を達成するために、概ね同じ放射線量がほぼ等しい信号レベル、例えば１ｍＲに対し出力されなければならないという事実の認識に基づいている。このステップは、高いエネルギー照射と低いエネルギー照射とから成る放射線量の比が全ての対応画像に対しほぼ一定のままであることを保証する。絶対放射線量は、Ｘ線検出器の最適なダイナミックレンジを達成するために、Ｘ線源の軌道に沿って変化され、患者の特徴的な断面が異なる方向において異なる厚さを持つ事実を考慮すれば、異なる方向からのノイズ寄与を平衡にすることに注意されたい。各Ｘ線エネルギーに対し同じ放射線量を送出するために、Ｘ線管に必要な電流の先験的に計算された値の組を使用し、この値は、ルックアップテーブルに記憶される。
【００１７】
本発明は、Ｘ線を生じさせるＸ線源と透過画像を検出する対応するＸ線検出器とを具備するＸ線装置を含むシステムにも関する。このＸ線検出器は、二次元の検出器であり、前記Ｘ線装置は、Ｘ線のエネルギーを変化させる第１の手段、二重エネルギー法に従って較正を行う較正手段及び透過画像に演算を行い、３Ｄ再構成アルゴリズムも実行する算術手段も含んでいる。
【００１８】
【発明の実施の形態】
本発明のこれら及び他の特徴は、対応する数字が対応する要素を意味する図を参照して以下に詳細に説明する。
【００１９】
図１は、本発明による方法を説明するフローチャートを図示している。第１の較正７の間に計算されるパラメタ(xpl),(xal)の組は、二重エネルギー法に従って、低いＸ線エネルギー及び高いＸ線エネルギーで取得される透過画像１、３の組に適用される。これは、２つの密度画像１１、１３を得るための非線形変換５となる。この第１の較正７は、図２ａ、２ｂ、２ｃを参照して以下に詳細に説明される。前記画像１１は、二重エネルギー法による第１の基本材料への分解画像である。画像１３は、二重エネルギー法による第２の基本材料への分解画像である。第２の較正９の間に決められ、図２ｃを参照して説明される軟組織削除角度φは、この軟組織を削除するために、前記密度画像１１、１３の結合体に適用される。得られた密度画像１５の組は、骨だけを含んでいる。骨３０の３Ｄ画像は、例えばフェルドカンプアルゴリズムのような３Ｄ再構成アルゴリズム２０を用いることにより形成される。骨のＣａ含有量の値は、画像３０を既知のＣａ含有量を持つ第３の較正ファントムにより決められる較正因子ｆに当てられることにより決められる。この対応する３Ｄ画像は、機能ブロック５０で表される。
【００２０】
図２ａは、二重エネルギー法に従って前記分解を実行するのに適当な第１のファントムを表した図である。この実施例において、この第１のファントムは、アルミニウムのウェッジ(wedge) （ｃ）から成り、これは、合成材料のウェッジ（ａ，ｂ）上に置かれる。この第１のファントムは、アルミニウムが骨のＸ線吸収特性にほぼ対応するＸ線吸収特性を持つのに対し、合成材料は軟組織の吸収作用を示すので、骨密度測定に用いられるとき、特に材料の分解を行うのに適している。較正パラメタ(xpl, xal)の組の取得自体は、既知であり、H. Neale Cardinal及びA. Fenster著、"An accurate method for direct dual-energy calibration and decomposition" Mad. Phys, 17(3), 1990 p.327に記載されている。この方法に従って計算される較正パラメタ(xpl, xal)の第１及び第２の組は、他の較正ファントム及び患者の透過画像に対しても実行されるような、第１の基本材料（合成材料）の寄与に対応する第１画像及び第２の基本材料（アルミニウム）の寄与に対応する第２画像を得るように、全ての後続する分解演算中にも用いられる。
【００２１】
図２ｂは、二重エネルギー法に従って軟組織削除角度を決める実施例を表した図である。初めに、例えば軟組織のファントムから成る第２のファントムの二重エネルギー画像が形成される。第１の較正を介して決められるパラメタxpl, xalの組を使用する場合、分解演算は、前記軟組織から合成材料及びアルミニウム７１、７２へ行われる。その上、合成材料及びアルミニウム７１、７２によるこれらの寄与に基づいて、全体ベクトル７３が合成材料／アルミニウム空間内に作られ、このベクトルは、軟組織を表す。前記ベクトル７３と水平軸との間の角度αは、この場合、軟組織に対し特徴的角度である特徴的な材料の角度を決める。軟組織による寄与を削除するために、ベクトル７３は、垂直線７４に投影される。前記水平軸とライン７４とにより囲まれる角度は、軟組織削除角度、すなわちφ＝α＋９０を構成する。この角度φは、本発明の方法に従って第３のパラメタとしてさらに使用される。
【００２２】
図２ｃは、患者データの組から軟組織を削除するのを図で説明したものである。二重エネルギー法に従って患者の透過画像を取得した後、パラメタ(xpl, xal)の第１及び第２の組は、材料の分解を決めるために適用される。次いで、全体ベクトル７８は、合成材料７５及びアルミニウム７６による関連する寄与に基づいて分解される。このベクトル７８は、軟組織７８’及び骨７８"による前記寄与の和も形成する。前記軟組織の寄与は、ベクトル７８をライン７４上に投影することにより削除される。この結果としてのベクトル７９は、骨による前記寄与を示す。
【００２３】
図３は、本発明による方法を実行するシステム１００を図示している。このシステムは、Ｘ線源１０１と２つのエネルギーからなるＸ線を生じさせる手段１０２を具備するＸ線装置１１０を含む。この生じたＸ線は、広域な空間角で放射される。手段１０２は、Ｘ線管に適用すべき、必要とされるキロボルト（ｋＶ）値を計算する手段１０２ａと、前記生じたＸ線に必要なフィルタリング演算を行う手段１０２ｂとを含む。前記Ｘ線装置１１０は、照射されるべき目的物１１４の周りを弧に沿ってＸ線源を移動させるための駆動手段１０３も含んでいる。生成される透過画像は、コンピュータ（図示せず）のメモリ１１５に記憶されるように二次元の検出器１１２によって検出される。本システム１００は、較正ファントム（図示せず）とパラメタxpl, xal, φ, ｆの計算を行うための他の算術手段１２２とを含む較正手段１２０も含んでいる。コンピュータ（図示せず）内に記憶されたコンピュータプログラムは、前記算術手段１２２の実施例を構成している。本システム１００は、骨の３Ｄ画像を形成するために３Ｄ再構成演算を行う算術手段１３０も含んでいる。この算術手段１３０の実施例は、例えば、フェルドカンプアルゴリズムに基づき、コンピュータ（図示せず）内に記憶されるコンピュータプログラムによって形成される。
【図面の簡単な説明】
【図１】  本発明による方法を説明するフローチャートの概略図である。
【図２ａ】  二重エネルギー法に従う分解を行う適切なファントムの概略表示である。
【図２ｂ】  二重エネルギー法に従う削除角の決定を実施例を概略的に説明する。
【図２ｃ】  患者のデータの組から軟組織の削除を概略的に説明する。
【図３】  本発明による方法を実行するシステムを概略的に示す。