JP2004518501A - 画像データの組におけるボリュームの密度を決める方法及びシステム - Google Patents

Abstract

二重エネルギー法を利用して骨密度計測を実行する方法であり、当該ボリューム内にある周囲の軟組織が較正手順を用いて透過画像から正確に削除される。この方法は、１つの実行中にデータ取得を実行するように、二次元のＸ線検出器を具備するＸ線装置を利用する。

Description

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

Claims (13)

1. 検査すべき目的物の画像データの組においてボリュームの密度を決める方法であって、前記ボリュームは、当該ボリュームが寄与し得る特徴パラメタを持つ材料を含むことにより環境と区別され、
   −異なる照射位置において前記環境と共に前記ボリュームをＸ線に照射させ、各照射位置において、前記Ｘ線の第１のエネルギー及び第２のエネルギーに対応する前記目的物の透過画像の第１の組をＸ線検出器で得るために、前記Ｘ線の前記第１のエネルギー及び前記第２のエネルギーで照射を行うステップを含む方法において、
   −二次元の検出器を前記Ｘ線検出器として使用するステップ、
   −二重エネルギー法に従って第１の較正を実行し、パラメタの第１の組とパラメタの第２の組とを得るステップ、
   −二重エネルギー法に従って第３の較正を実行し、前記環境を特徴とする第３のパラメタと前記材料を特徴とする第４のパラメタとを得るステップ、
   −前記パラメタの第１の組と前記パラメタの第２の組とを用いながら、透過画像の前記第１の組に対し第１の演算を行い、密度画像の組を得るステップ、
   −前記密度画像の組に前記第３のパラメタを適用し、前記密度画像の組から前記環境を削除するステップ、
   −前記密度画像の組の情報を用いて３Ｄ再構成アルゴリズムを実行し、前記ボリュームの画像を得るステップ、及び
   −前記ボリュームの画像に前記第４のパラメタを適用し、前記ボリュームにおける前記材料の密度を計算するステップ
   も含む方法。
2. 請求項１に記載の方法において、透過画像の第２の組を得るために、前記密度画像に対し演算を行い、フェルドカンプアルゴリズムが３Ｄ再構成アルゴリズムとして用いられる方法。
3. 請求項１に記載の方法において、前記第３のパラメタは、軟組織の削除角度であり、前記第４のパラメタは、骨材料密度因子である方法。
4. 請求項３に記載の方法において、使用すべき前記骨密度因子の値は、ルックアップテーブルで調べられる方法。
5. 請求項３に記載の方法において、該当するボリュームは、前記ボリュームの画像において規定され、前記該当するボリュームは、柱状骨のボリュームの少なくとも一部に対応する方法。
6. 請求項１に記載の方法において、前記密度画像の組のピクセル値は、ダイナミックレンジを高めるために、因子Ｍで乗算される方法。
7. 請求項１に記載の方法において、第１及び第２の較正材料は、前記二重エネルギー法に従って前記第１の較正を実行するために使用され、前記第１の較正材料の実効原子番号は、１０よりも低いのに対し、前記第２の較正材料の実効原子番号は、１０よりも大きい方法。
8. 請求項１に記載の方法において、決められるべき前記密度は、骨密度であり、前記第１のエネルギーのＸ線を生じさせるためにＸ線管に必要な電圧は、５０〜８０ｋＶの範囲にあるのに対し、前記第２のエネルギーのＸ線を生じさせるためにＸ線管に必要な電圧は、１００から１５０ｋＶの範囲にある方法。
9. 請求項８に記載の方法において、前記第１のエネルギー及び前記第２のエネルギーに対応するＸ線を生じさせるＸ線管の電圧差は、４０ｋＶと５０ｋＶとの間にある方法。
10. 請求項１乃至９の何れか一項に記載の方法において、Ｘ線源は前記目的物の周りを約１８０°の弧に沿って回転され、各エネルギーに対し、約５０の透過画像が前記弧において異なる角度で取得される方法。
11. 請求項１乃至１０の何れか一項に記載の方法において、低いエネルギー照射に対する放射Ｘ線量は、対応する高いエネルギー照射に対する放射線量と同じ大きさの次数である方法。
12. 請求項１に記載の方法を実行するシステムにおいて、前記システムは、Ｘ線を生じさせるＸ線源を備えるＸ線装置と透過画像を検出する対応するＸ線検出器とを有し、前記Ｘ線検出器は二次元の検出器であり、前記Ｘ線のエネルギーを変化させる第１の手段、前記二重エネルギー法に従って較正を実行する較正手段、及び前記透過画像に演算を行い、３Ｄ再構成アルゴリズムも実行する算術手段を備えるシステム。
13. 請求項１２に記載のシステムにおいて、前記較正手段は、適当な第１の較正ファントムと適当な第２の較正ファントムとを含み、さらに前記パラメタを計算する他の算術手段を含むシステム。