(57)【特許請求の範囲】
【請求項１】
Ｘ線源及びピクセル化Ｘ線検出器を含み、被写体である物質にＸ線を照射してＸ線画像を生成するＸ線装置であって、該装置は：
物質の種類及び物質の厚さを表す値に関する第１データベース、
物質の種類及び物質の厚さに関連付けられる散乱放射線値に関する第２データベース、
ｉ）前記Ｘ線検出器の出力信号と、前記第１データベースの値とを比較して、前記第１データベースから最も近い物質の種類及び物質の厚さを出力し；
ｉｉ）前記第２データベースから、ステップｉにおける出力、又は、アルゴリズムにおける前の繰返しで識別された物質の種類及び物質の厚さ、と関連付けられる散乱放射線を導出し；
ｉｉｉ）前記散乱放射線を、前記Ｘ線検出器の出力信号から除去し；
ｉｖ）ステップｉｉｉの出力に対して、物質の種類及び物質の厚さ識別ステップを実行し；
ｖ）アルゴリズムのステップｉｉ、ｉｉｉ、及び、ｉｖの少なくとも１回の更なる繰返しにおいて、識別ステップで識別された物質の種類及び物質の厚さを利用する、
アルゴリズムを実行するよう構成されたプロセッサ
を更に含み、
修正されたＸ線出力信号が最適化されていると判定されたとき、前記繰返しは停止され、前記修正されたＸ線出力信号が最適化されているか否かを判定するステップは、
現繰返しにおいて識別された物質の種類及び物質の厚さを示す値を、前の繰返しにおいて識別された物質の種類及び物質の厚さを示す値、又は複数の前の繰返しにおいて識別された物質の種類及び物質の厚さを示す値の平均と比較すること；及び
前記現繰返しにおいて識別された物質の種類及び物質の厚さを示す値が、前記前の繰返しにおいて識別された物質の種類及び物質の厚さを示す値、又は複数の前の繰返しにおいて識別された物質の種類及び物質の厚さを示す値の平均の閾値範囲内であるかを判定すること；を含み、
前記比較された値が前記閾値内にある場合、前記アルゴリズムのステップ（ｉｉ）〜（ｉｖ）の再繰返しは、停止される、
Ｘ線装置。
【請求項２】
前記アルゴリズムは、除去された散乱放射線のＸ線光子が散乱しなければ相互作用したであろう前記検出器の空間的位置の出力信号に、前記除去された散乱放射線を加えることによって、前記Ｘ線検出器の前記出力信号を修正する更なるステップを含む、請求項１に記載のＸ線装置。
【請求項３】
前記Ｘ線出力信号は、前記Ｘ線検出器の各画素に関するＸ線出力信号である、請求項１又は２に記載のＸ線装置。
【請求項４】
前記アルゴリズムは、前記各画素に関するＸ線出力信号に対して実行される、請求項３に記載のＸ線装置。
【請求項５】
物質の種類及び物質の厚さに関連付けられる散乱放射線値の前記第２データベースの散乱値は、散乱カーネルである、請求項１乃至４の何れかに記載のＸ線装置。
【請求項６】
前記データベースは、対象となる物質の範囲の一部に関する素性及び厚さを表す散乱カーネルを追加される、請求項５に記載のＸ線装置。
【請求項７】
前記Ｘ線検出器の出力信号から散乱放射線を除去するステップは、前記検出器の各画素に関して、ステップｉｉで前記データベースから選択された散乱カーネルから散乱カーネルを補間するステップを含む、請求項６に記載のＸ線装置。
【請求項８】
各画素に関連付けられる補間された散乱カーネルからＸ線検出器出力信号全体に対する散乱推定値を生成する更なるステップを含む、請求項７に記載のＸ線装置。
【請求項９】
被写体である物質にＸ線を照射して生成されたＸ線画像におけるコントラスト対ノイズ比を改善する方法であって、該方法は：
ｉ）Ｘ線検出器の出力信号を、物質の種類及び物質の厚さを表す第１データベース値にある値と比較し、前記第１データベースから最も近い物質の種類及び物質の厚さを出力するステップ；
ｉｉ）物質の種類及び物質の厚さに関連付けられる散乱放射線値に関する第２データベースから、ステップｉにおける出力、又は、前の繰返しで識別された物質の種類及び物質の厚さと関連付けられる散乱放射線を導出するステップ；
ｉｉｉ）前記散乱放射線を、前記Ｘ線検出器の前記出力信号から除去するステップ；
ｉｖ）ステップｉｉｉの出力に対して、物質の種類及び物質の厚さ識別ステップを実行するステップ；
ｖ）アルゴリズムのステップｉｉ、ｉｉｉ、及び、ｉｖの少なくとも１回の更なる繰返しにおいて、識別ステップで識別された物質の種類及び物質の厚さを利用するステップ、
を含み、
修正されたＸ線出力信号が最適化されていると判定されたとき、前記繰返しは停止され、前記修正されたＸ線出力信号が最適化されているか否かを判定するステップは、
現繰返しにおいて識別された物質の種類及び物質の厚さを示す値を、前の繰返しにおいて識別された物質の種類及び物質の厚さを示す値、又は複数の前の繰返しにおいて識別された物質の種類及び物質の厚さを示す値の平均と比較するステップ；及び
前記現繰返しにおいて識別された物質の種類及び物質の厚さを示す値が、前記前の繰返しにおいて識別された物質の種類及び物質の厚さを示す値、又は複数の前の繰返しにおいて識別された物質の種類及び物質の厚さを示す値の平均の閾値範囲内であるかを判定するステップ；を含み、
前記比較された値が前記閾値内にある場合、前記アルゴリズムのステップ（ｉｉ）〜（ｉｖ）の再繰返しは、停止される、
方法。
【請求項１０】
前記除去された散乱放射線のＸ線光子が散乱しなければ相互作用したであろう前記検出器の空間的位置での出力信号に、前記除去された散乱放射線を加えることによって、前記Ｘ線検出器の前記出力信号を修正する更なるステップを含む、請求項９に記載の方法。
【請求項１１】
前記Ｘ線検出器の出力信号は、前記Ｘ線検出器の画素からの出力信号であり、請求項９又は１０の方法ステップは、各画素の出力に対して実行される、請求項９又は１０に記載の方法。
【請求項１２】
物質の種類及び物質の厚さに関連付けられる散乱放射線値の前記第２データベースの散乱値は、散乱カーネルである、請求項９乃至１１の何れかに記載の方法。
【請求項１３】
前記データベースは、対象となる物質の範囲の一部に関する素性及び厚さを表す散乱カーネルを追加される、請求項１２に記載の方法。
【請求項１４】
前記Ｘ線検出器の出力信号から散乱放射線を除去するステップは、前記検出器の各画素に関して、ステップｉｉで前記データベースから選択された散乱カーネルから散乱カーネルを補間するステップを含む、請求項１３に記載の方法。
【請求項１５】
各画素に関連付けられる補間された散乱カーネルからＸ線検出器出力信号全体に対する散乱推定値を生成する更なるステップを含む、請求項１４に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、散乱線除去グリッドを通常利用するＸ線機器、特に、散乱Ｘ線をＸ線画像分析に利用するために提供する装置に関する。
【背景技術】
【０００２】
対象物がＸ線照射を受ける際に、Ｘ線光子の一部は吸収され、また一部は、対象物を通過し、散乱されず、Ｘ線検出器に衝突する。これは「直接放射線」と呼ばれる。一部のＸ線は吸収され、残りは散乱される。発生した散乱の強度は、検出器によって検出される直接放射線量を超えることがある。散乱放射線の結果、コントラストが低下し、ノイズが増大するために、画質が悪くなる。
【０００３】
吸収されたＸ線は、Ｘ線画像でコントラストを提供する。散乱Ｘ線光子が、検出器に当たると、散乱Ｘ線光子がどこから来たかを識別できないため、画像におけるランダムノイズが増大される。
【０００４】
散乱の問題を解決するのに用いられる最も広く採用されている技術は、Ｘ線検出器と被験対象物との間に散乱線除去グリッドを設置することである。散乱線除去グリッドは、Ｘ線吸収物質から形成される一連の離間した平行な薄層（ｌａｍｅｌｌａ）を含む。散乱Ｘ線の大部分は、薄層の１つと係合し、吸収される。従って、散乱線除去グリッドが存在する場合にＸ線検出器によって検出されるのは、主に直接放射線となる。
【０００５】
最初の散乱線除去グリッドは、米国特許第１１６４９８７号（Ｂｕｃｋｙ）（特許文献１）に記載されている。
【０００６】
散乱線除去グリッドに関する問題の１つは、検出画像における散乱の影響を軽減するだけでなく、Ｘ線吸収薄層が、直接放射線の一部、即ち薄層の経路で伝搬する光子を吸収してしまうことである。
【０００７】
散乱線除去グリッドにおいて失われた光子を補完する、つまり低下した画質を補うために、Ｘ線束を増加することが慣用技術になっている。しかしながら、これは、Ｘ線画像処理がＸ線感光物質のものである場合、欠点となる。これは、散乱線除去グリッドが存在することで補完するために患者へのＸ線放射線量が増大されなければならない場合、医療画像処理において、最大の懸念事項となる。
【０００８】
Ｘ線画像処理で使用されるＸ線パワーを減少するために、幾つかの試みが行われた。
【０００９】
米国特許第７５５１７１６号（特許文献２）では、散乱線除去グリッドを使用する代わりに、散乱Ｘ線光子を大凡判定するのに数学的方法を使用している。散乱Ｘ線光子を大凡判定するのに数学的方法を利用することによって、散乱線除去グリッドを使用するＸ線装置と比べると、Ｘ線量は減少されることができる、又は信号対ノイズ比は、増大される。
【００１０】
多吸収板を含むＸ線装置が、本出願人の英国特許出願第２４９８６１５号公報（特許文献３）に記載されている。このＸ線装置では、Ｘ線エネルギスペクトルは、多様に摂動される。この装置及び方法は、物質特性を識別可能にする情報を提供する。特許文献３は、参照することにより本明細書に組込まれる。
【００１１】
先行技術は、散乱放射線を大幅に減少するために散乱線除去グリッドを使用することによって、又はＸ線画像から散乱放射線がどのようになるかの推定値を除去するために数学的補正を適用することによって、散乱放射線によって引き起される問題を解決する。
【００１２】
本発明は、直接放射線と散乱放射線の両方を測定することを伴う新たな方法で散乱の問題を解決する。これは、直接放射線と散乱放射線の両方に作用するよう装置に変更を加えることによって、達成される。
【００１３】
散乱放射線は、これまで、Ｘ線画像に悪影響を与えるものと見なされてきた。散乱に関する問題を解決する際の以前の試みは、画像が、Ｘ線検出器に衝突する直接Ｘ線放射線によって出来るだけ形成されるように、散乱放射線の除去を最適化することに重点が置かれていた。
【００１４】
散乱放射線を測定するために提供するだけでなく、本発明は、散乱放射線を、物質及び物質の厚さを識別する際に、及びコントラスト対ノイズ比を改善する際に使用するためにも提供する。
【００１５】
Ｘ線検出器からの測定出力信号を、散乱放射線の独立した測定値とを比較することによって、直接放射線と散乱放射線を含む画像中で散乱放射線を識別でき、且つ画像から散乱放射線を除去できる。それにより、コントラスト対ノイズ比を高められる。医療用途、及びＸ線が、Ｘ線によって損なわれる可能性がある物質を分析するのに使用される他の用途に関連して、本発明は、線量を、同様の水準の画像を作成するのに低減可能にする、又は同じＸ線量を使用して、より良好な画像が生成できる。
【先行技術文献】
【特許文献】
【００１６】
【特許文献１】米国特許第１１６４９８７号
【特許文献２】米国特許第７５５１７１６号
【特許文献３】英国特許出願第２４９８６１５号
【発明の概要】
【課題を解決するための手段】
【００１７】
本発明の第１態様によると、Ｘ線源及び多画素Ｘ線検出器、及びＸ線源と検出器間に、入射するＸ線光子を摂動するよう構成される部材を含むＸ線装置であり、部材は、該部材に入射するＸ線光子を、該Ｘ線光子の入射角度に応じて、異なって摂動するように適合され、部材は、少なくとも部分的にＸ線光子に対して半透明を示す物質から形成された複数の要素を含むＸ線装置であって、該装置は：
物質の種類及び／又は物質の厚さを表す値に関する第１データベース、
物質の種類及び／又は物質の厚さを示す散乱放射線値に関する第２データベース、
ｉ）Ｘ線検出器の各画素の出力信号と第１データベースの値とを比較して、第１データベースから最も近い物質及び／又は厚さを出力し；
ｉｉ）第２データベースから、ステップ（ｉ）における物質の種類及び／又は厚さと関連付けられる散乱放射線を選択し；
ｉｉｉ）散乱放射線を、Ｘ線検出器の出力信号から除去する
アルゴリズムを実行するよう構成されたプロセッサ
を更に含むＸ線装置が提供される。
【００１８】
有利には、アルゴリズムは、除去された散乱放射線のＸ線光子が散乱しなければ相互作用したであろう検出器の空間的位置での出力信号に、除去された散乱放射線を加えることによって、Ｘ線検出器の出力信号を修正するステップを含む。
【００１９】
物質の種類及び／又は厚さを表す値に関する第１データベースは、例えば、散乱線除去グリッドを利用することによって、散乱を除外して、実質的に作成されてもよい。
【００２０】
アルゴリズムは、各画素のＸ線出力信号に対して、実行されてもよい。
【００２１】
好適には、アルゴリズムは、アルゴリズムのステップｉｉ及びｉｉｉの少なくとも１回の更なる繰返しで識別された物質及び／又は厚さを利用して、ステップｉｉｉの出力に対して、物質及び／又は厚さ識別ステップを実行するステップを含む。
【００２２】
有利には、アルゴリズムは、修正されたＸ線出力信号が最適化されているか否かを判定する更なるステップを含む。
【００２３】
修正されたＸ線出力信号が最適化されているか否かを判定するステップは、現繰返しの画像における物質の種類及び厚さを示す値を、前の繰返しの画像における物質の種類及び／又は厚さを示す値、又は多数の前の繰返しにおける物質の種類及び／又は厚さを示す値の平均と比較すること、及び現繰返しの画像における物質の種類及び／又は厚さを示す値が、前の繰返しの物質の種類及び／又は厚さを示す値、又は多数の前の繰返しにおける物質の種類及び／又は厚さを示す値の平均の閾値範囲内であるかを判定することを含んでもよい。比較された値が閾値内にある場合、アルゴリズムのステップ（ｉ）〜ステップ（ｉｉｉ）の再繰返しは、停止されてもよい。
【００２４】
アルゴリズムは、除去された散乱放射線のＸ線光子が散乱しなければ相互作用したであろう検出器の空間的位置での出力信号に、除去された散乱放射線を加えることによって、Ｘ線検出器の出力信号を修正する更なるステップを含んでもよい。
【００２５】
好適には、物質の種類及び／又は厚さを示す散乱放射線値に関する第２データベースの散乱値は、散乱カーネルである。
【００２６】
散乱カーネルは、物質を通るＸ線光子の単一軌道からの散乱を表わすものとして定義されてもよい。
【００２７】
データベースは、対象となる物質の範囲の一部に関する素性及び／又は厚さを表す散乱カーネルを追加されてもよい。
【００２８】
Ｘ線検出器の出力信号から散乱放射線を除去するステップは、検出器の各画素に対して、ステップｉｉでデータベースから選択された散乱カーネルから散乱カーネルを補間するステップを含んでもよく、各画素に関連付けられる補間された散乱カーネルからＸ線検出器の出力信号全体に関する散乱推定値を生成する更なるステップも含んでもよい。
【００２９】
有利には、部材は、多吸収板（ＭＡＰ：ｍｕｌｔｉ−ａｂｓｏｒｐｔｉｏｎ  ｐｌａｔｅ）であり、該ＭＡＰは、複数の異なる領域を含む板であり、異なる領域は、入射Ｘ線エネルギに対して異なる透明性を有する。複数の異なる領域は、隣接する領域が入射Ｘ線エネルギに対して異なる透明性を有する、繰返しパターンで形成されてもよい。
【００３０】
部材は、画像キャプチャ中、静止位置に保持されるのが好ましい。
【００３１】
本発明の第２態様によると、Ｘ線源及び多画素Ｘ線検出器、及びＸ線源と検出器間に、入射するＸ線光子を摂動するよう構成される部材を含むＸ線装置であり、部材は、装置の少なくとも２つの異なる構成を提供するように適合され、部材は、異なる構成で、Ｘ線光子を異なって摂動するＸ線装置を提供する。
【００３２】
好適には、部材は、上記装置の少なくとも２つの異なる構成を提供するように移動可能である。
【００３３】
一実施形態では、部材は、多吸収板（ＭＡＰ）であり、該ＭＡＰは、複数の異なる領域を含む板であり、各領域は、入射Ｘ線エネルギに対して異なる透明性を有する。かかるＭＡＰは、１軸で、２軸又は３軸で移動可能としてもよく、これらの軸は、互いに直交してもよい。
隣接する領域が入射Ｘ線エネルギに対して異なる透明性を有する、繰返しパターンで形成されてもよい。
【００３４】
部材を移動することによって、どれか一つのＸ線光子が衝突する可能性がある部材の部分は、構成別に異なり、従って、装置の構成別に散乱パターンは異なる。
【００３５】
直接放射線の場合、既知の異なる構成間で移動する効果は、常に同じである。即ち、構成（ａ）から構成（ｂ）に部材を移動すると、直接放射線に対して常に同じ効果を有する。構成（ｂ）から構成（ｃ）に部材を移動しても、直接放射線に対して常に同じ効果を有するが、その効果は、構成（ａ）から構成（ｂ）に移動する効果とは異なるかも知れない。この知識を用いて、直接放射線を識別できる、従って、散乱放射線も、直接放射線でないものという理由で、識別できる。
【００３６】
別の実施形態では、部材は、個々の薄層がフレーム内で回転可能に取付けられることを除いて、散乱線除去グリッドに酷似する。フレーム内で薄層を回転させると、装置の構成は変化するが、これは、薄層がＸ線源からＸ線検出器に延伸する軸と平行になる散乱パターンが、薄層が上記軸と平行から例えば５度傾く散乱パターンと異なるためである。
【００３７】
別の実施形態では、ＭＡＰを移動する、又は散乱線除去グリッドの薄層を、回転するように構成する代わりに、部材は、電流等の信号を印加すると平坦構成と屈曲構成との間で変化する形状記憶合金から形成されてもよい。
【００３８】
別の実施形態では、部材は、電流を印加すると大きさを変える圧電物質から形成される。
【００３９】
部材は、フレーム内で回転可能に取付けられる個々の薄層を含んでもよい。
【００４０】
部材は、信号の印加時に形状を変えることによって構成間で変化する形状記憶合金から形成されてもよい。部材は、信号の受信時に平坦構成と屈曲構成との間で変化する形状記憶合金から形成されてもよい。
【００４１】
部材は、電流の印加時に大きさ及び／又は形状を変化させる圧電物質から形成されてもよい。
【００４２】
本発明の別の態様によると、Ｘ線源及び多画素Ｘ線検出器、及びＸ線源と検出器間に、入射するＸ線光子を摂動するよう構成される部材を含むＸ線装置であり、部材は、該部材に入射するＸ線光子を、該Ｘ線光子の入射角度に応じて、異なって摂動するように適合され、部材は、少なくとも部分的にＸ線光子に対して半透明を示す物質から形成された複数の要素を含むＸ線装置が提供される。
【００４３】
好適には、上記要素は、Ｘ線源からＸ線検出器に延伸する軸の方向に、該軸と直角を成す方向よりは、更に延伸する。
【００４４】
散乱Ｘ線光子が、要素の１つに衝突すると、散乱Ｘ線光子が、該要素を通過するか、該要素によって吸収されるか、及び散乱Ｘ線光子が該要素を通過する場合、該光子のエネルギは、Ｘ線光子の該要素に対する入射角に応じて決まる。これは、Ｘ線光子の入射角によって、入射したＸ線光子が通過しなければならない物質の厚さが決まるからである（勿論、該要素の角に最も近くに入射するＸ線光子は、同じ入射角で、該要素の角から最も遠くに入射した場合よりも、物質の厚さが薄くなるかも知れない）。
【００４５】
本発明の他の態様によれば、Ｘ線画像から散乱放射線を除去する方法を提供し、該方法は、本発明の第１態様の装置において、被験物質を、Ｘ線光子に当て、試験中に、部材の構成を変化させるステップ；Ｘ線画像を分析し、Ｘ線画像における直接放射線と散乱放射線を識別し、Ｘ線画像から識別された散乱放射線を除去するステップを含む。
【００４６】
Ｘ線画像から散乱放射線を除去する方法は、物質の種類及び／又は物質の厚さを識別する更なるステップを含んでもよい。
【００４７】
物質の種類及び／又は厚さは、除去された散乱放射線から識別されてもよい。Ｘ線画像は、識別された物質の種類及び／又は厚さに従い、修正されてもよい。
【００４８】
本発明の他の態様によれば、散乱Ｘ線放射線と直接Ｘ線放射線の両方を含む原Ｘ線画像において散乱放射線から物質及び／又は物質の厚さを識別する方法が提供される。
【００４９】
未知の物質及び／又は厚さが、本発明の装置を使用して試験される場合、結果として得られた散乱パターンは、最も近い物質及び／又は厚さが識別され得るデータベース中の散乱パターンと比較され、それにより物質及び／又は厚さを識別できる。
【００５０】
物質及び／又は厚さの識別を実行するために、既知の厚さの既知の物質に、部材の構成が試験中に変化される本発明の第２態様の装置で、Ｘ線光子を当ててもよい。これにより、部材の異なる構成と一致する異なる散乱パターンを、既知の物質及び／又は厚さに対して生成できる。これは、データベースに保存されることができる。これは、画素毎に行われてもよい、即ち、多画素検出器の各画素の出力に対して行われてもよい。
【００５１】
また、既知の物質の種類及び／又は厚さに対する散乱パターンも、データベースに記録され、保存されてもよい。
【００５２】
図面は、本発明の好適実施形態を説明しており、例として挙げたものである。
【図面の簡単な説明】
【００５３】
【図１】本発明による装置の略図である。
【図２】本発明の部材の実施例を示している。
【図３】本発明の別の部材の実施例を示している。
【図４】本発明の更に別の部材の実施例を示している。
【図５】本発明の更に別の部材の実施例を示している。
【図６】本発明の更に別の部材の実施例を示している。
【図７】本発明の更に別の部材の実施例を示している。
【図８ａ】本発明の別の実施形態を示している。
【図８ｂ】本発明の別の実施形態を示している。
【図９】図８ａ及び図８ｂに示された本発明の実施形態の動作に関する方法を示すフローチャートである。
【図１０】図８ａ及び図８ｂに示された本発明の実施形態の動作に関する方法を示す別のフローチャートである。
【図１０ａ】図８ａ及び図８ｂに示された本発明の実施形態の動作に関する別の方法を表すフローチャートである。
【図１１ａ】図８ａ〜図１０に示された本発明の実施形態による処理前に検出された試料の画像を示している。
【図１１ｂ】図８ａ〜図１０に示された本発明の実施形態による処理後に検出された試料の画像を示している。
【図１２ａ】図１１ａに示された画像に関するカウント数対強度のヒストグラムを示している。
【図１２ｂ】図１１ｂに示された画像に関するカウント数対強度のヒストグラムを示している。
【図１２ｃ】散乱再割当前後の多数の対象領域に関するコントラスト対ノイズ比を示すグラフである。
【図１３】板を図８ａ及び図８ｂのように指向させた状態で、上から見た多吸収板の図である。
【発明を実施するための形態】
【００５４】
図１は、本発明による装置１の一般配置図を示している。Ｘ線源２は、検出器３と、この場合多吸収板又は散乱線除去グリッドである部材４が、検出器３の線源側で、検出器３に隣接して位置する状態で、位置合せされる。被験試料５は、線源２と多吸収板（ＭＡＰ）又は散乱線除去グリッド４との間に位置する。
【００５５】
図２〜図７では、部材に関する多数の異なる選択し得る仕組み（ａｒｒａｎｇｅｍｅｎｔ）を示している。
【００５６】
図２では、部材４は、一連の左右対称の段４ａ、４ｂ、４ｃを含むＭＡＰであり、それによりＭＡＰの中心部分４ｃは、ＭＡＰの最外側部分と比べると、物質の厚さが３倍になっている。ＭＡＰは、図２の左側に示されているように、横に移動されてもよい、又は図２の右側に示されているように、上方及び下方に移動されてもよい。次に、左側の図を見ると、本装置の６つの異なる構成を示している。最上より下の各図は、立体遮蔽（ｓｏｌｉｄ  ｓｈａｄｉｎｇ）状態にあるＭＡＰ４の新たな位置を示し、最上図では、陰影状態にあるＭＡＰの位置を示している。勿論、横に移動するＭＡＰ４よりむしろ、それに加えて又はその代わりに、ＭＡＰは、前後に移動してもよい。
【００５７】
右側の図を見ると、ＭＡＰの５つの異なる構成を示している。中央の図は、移動する前のＭＡＰを示している。中央図より上の図は、垂直方向に２段上方へ移動したところを示し、中央図より下の図は、垂直方向に２段下方に移動したところを示している。
【００５８】
図２の図面から分かるように、ＭＡＰが水平方向又は垂直方向に移動されると、特定の経路で伝搬するＸ線光子が通過しなければならない物質の厚さは、ＭＡＰの位置が変化するのに応じて、変化する。
【００５９】
図３では、部材４は、互いに重なる複数の板６ａ、６ｂから成る。各板６ａ、６ｂは、グリッド状の物質で作製される。この実施形態では、上板６ａは、下板６ｂより薄くなっている。最上図では、グリッド６ａ、６ｂは整列されている。最上図のすぐ下の図では、上グリッド６ａは、左へ移動されており、上グリッド６ａの物質が、下グリッドの間隙部に部分的に重なっている。上グリッドが更に移動されると、上グリッド６ａは、下グリッド６ｂの間隙部に重なるように移動し、更に移動すると、最も下の２図に示されるように、上グリッド６ａは、上グリッド６ａの物質の３部分が下グリッドの物質と、下グリッドの最右部分が上グリッドによって被覆されない状態で、整列されるように移動する。勿論、この実施形態で重要なことは、板６ａと板６ｂ間の相対移動である。かかる移動は、上及び／又は下板６ａ、６ｂに関するものである。    
【００６０】
ここでもまた、当業者は、被験物質によって散乱されるＸ線光子が、上下グリッドが互いに対してどのような位置であるかに応じて異なる影響を受けることが分かるであろう。
【００６１】
図４では、部材４は、フレーム（図示せず）内に、該フレームに対する相対的回転が可能であるように、取付けるよう構成される薄層７ａ、８ａの上下組７、８を含む。薄層７ａ、８ａが回転すると、特定の経路で伝搬するＸ線光子（散乱Ｘ線光子又は直接Ｘ線光子とされる）が通過しなければならない物質の厚さは、変化する。図４の４図から分かるように、薄層７ａ、８ａは、薄層７ａ、８ａがＸ線源から検出器に延伸する軸と略同じ平面に存在する位置から、薄層７ａ、８ａが上記軸に対して垂直になる位置に移動できる。    
【００６２】
薄層７ａ、８ａを形成する物質は、Ｘ線光子に対して不透明としてもよい、又はＸ線光子に対して半透明としてもよい。薄層を半透明にする効果については、図５を参照して、以下で更に詳細に説明される。
【００６３】
図５は、直交する２方向に延伸する薄層９ａ、１０ａを有する散乱線除去グリッドの形をした部材４を示している。左から右に延伸する薄層９ａは、Ｘ線光子に対して半透明であるが、薄層９ａに対して垂直に延伸する薄層１０ａは、Ｘ線光子に対して不透明である。線源からＸ線検出器に延伸する軸に対して微小角度の半透明な薄層に衝突する散乱Ｘ線光子は、かなり大きな角度で、例えば、上記軸に対して垂直に傾斜する半透明な薄層に衝突する散乱Ｘ線光子よりも、遥かに厚い物質を通過する必要があるだろう。かかるＸ線光子は、殆ど減衰せずに、薄層を通過する可能性がある。従って、散乱Ｘ線光子は、半透明薄膜に対する入射角に応じて減衰されると理解される。
【００６４】
図５に示される部材４は、静的又は動的配置で、使用されてもよい。
【００６５】
動的配置では、部材４は、移動されてもよいように取付けられ、この移動は、図２及び図３を参照して説明されたのと同様な方法で行われてもよい。これらの図面を参照して説明されたように、図５を参照して説明された部材４を特定の方法で移動する直接Ｘ線光子への効果は、常に同じであり、この予想される変化は、直接Ｘ線光子、従って散乱Ｘ線光子を識別するのに使用されることができる。
【００６６】
静的配置では、部材は、異なる構成を有さないため、該部材の直接Ｘ線光子に対する効果は、可変ではない。しかしながら、散乱Ｘ線光子を識別可能である。
【００６７】
異なる物質の種類及び厚さに関して、散乱Ｘ線光子のパターンは、物質及び厚さに特有であることが知られている。異なる既知の厚さの異なる既知の物質を試験することによって、Ｘ線検出器出力に関するデータベースが、構築できる。かかるＸ線検出器出力は、各画素で記録される信号を含む。未知の物質の試料が装置１で試験される場合、試料５の物質種類及び／又は厚さの識別は、隣接する画素で受信された信号、又は隣接する画素で受信された信号間の差異の形をしたデータを、信号又は異なる既知の物質種類及び／又は厚さと関連付けられた信号間の差異に関するデータベースと比較し、検出された信号又は信号間の差異と最も合致するデータベースからの情報に基づいて、試料の物質種類及び／又は厚さに関する判定を行うことによって、達成されてもよい。
【００６８】
図６と図７は、部材１の構成が、該部材を形成する物質に対して物理的変化を与えることによって変化される点で同様である。図６は、４薄層１１ａから形成された部材１を示しており、各層は、付与される電荷に応じて形状を変化する圧電物質製である。一構成（最上図）では、薄層１１ａは、深くて狭く、平行な薄層間の間隙が、個々の薄層１１ａの幅より大幅に広くなっている。電荷が付与され、圧電物質が形状を変化すると、各薄層は、より浅く、幅広になる。付与される電荷が更に変化されると、物質は、図６の最下図で示されたように、薄層１１ａが纏って連続した平面となるまで、形状を変化する。
【００６９】
図７の実施形態では、部材１は、形状記憶合金から形成され、最上図で示された平面構成と、湾曲領域１２ａ、１２ｂを提供する波形構成との間で動く。印加される電流に応じて、部材１は、図示された２構成間で形状を変化する。
【００７０】
図２、図３、図６及び図７では、同じ矢印ａ、ｂ、ｃ及びｄが示されている。これらは、Ｘ線光子の経路、及び該光子が部材の構成に応じて衝突する部材１の物質を表している。当業者は、部材の構成を変えることで、Ｘ線光子に対する部材の効果が大きくなり得ることを容易に理解するであろう。
【００７１】
図８ａは、本発明の別の実施形態による装置１を示している。Ｘ線源２は、多吸収板（ＭＡＰ）４が検出器３の線源側で検出器に隣接して位置した状態で、検出器３と位置合せされている。ＭＡＰ４は、図２を参照して説明した種類のもの、又は図１３を参照して説明したようなＭＡＰ１６としてもよい。図２で示された本発明の実施形態と、図８ａ及び図８ｂで示された実施形態との主な相違点は、ＭＡＰ４の動きである。図２の実施形態では、ＭＡＰ４は、線源２及び検出器３に対して移動するのに対して、図８ａ及び図８ｂの実施形態では、ＭＡＰ４は静止している。
【００７２】
図１３で示されたＭＡＰ１６は、異なる領域１６ａ〜１６ｉを含む。板の異なる領域は、異なるＸ線吸収能力を有する。これらの領域は、異なる厚さを有してもよい、又は、例えば、異なる物質から形成されてもよい。
【００７３】
図８ａは、被験試料５に入射するＸ線光子に起きることを示している。上記で説明されたように、Ｘ線源２から放出されたＸ線光子ｍの中、一部のＸ線光子ｍ’は、被験試料５を真直ぐ通過して、検出器３に衝突しており、これらのＸ線光子は、「直接放射線」を表し、一部は吸収され、一部のＸ線光子ｎは散乱される。
【００７４】
散乱Ｘ線光子ｎは、望ましくない。本発明のこの実施形態は、図８ａで示された散乱Ｘ線光子ｎを画像から除去可能にし、疑似直接放射線ｎ’として、散乱Ｘ線光子ｎが関連付けられた直接放射線ｍ’に再割当可能にする。その結果、図８ｂにおいて、Ｘ線光子が衝突して生じる検出器３の任意の一画素に関する出力は、直接放射線ｍ’及び再割当された疑似直接放射線ｎ’を含むことになる。これは、コントラストを高め、従って、検出器３によって生成された画像のコントラスト対ノイズ比を高めると共に、より望ましい画像を提供するが、というのも、散乱Ｘ線光子ｎが、散乱線除去グリッドを用いた場合でのようには、除去されず、Ｘ線光子が散乱しなければＸ線光子が相互作用したであろう検出器の空間的位置での出力信号に追加されるからである。そのため、結果的に、画像コントラストを改善できる。
【００７５】
検出された画像から散乱光子ｎを除去し、該散乱光子ｎを疑似直接放射線ｎ’として再割当するために、検出器３によって検出された信号は、以下のように処理される。
【００７６】
異なる物質及び該物質の異なる厚さの散乱パターンに関するデータベースが、作成される。これは、図１０で示されたシミュレーションモデルを使用して行われ、該シミュレーションモデルでは、異なる物質及び該物質の異なる厚さが、Ｘ線光子のペンシルビームが該物質に入射する際に、どのように散乱するかをシミュレートする。シミュレーションモデルを使用する代わりに、データベース用データは、試料を異なる物質の種類及び／又は厚さのものとして、実際の装置及び試料を使用して収集されることもできる。
【００７７】
多画素検出器の１画素と位置合せされるペンシルビームを使用することによって、線源が位置合せされる画素に達するＸ線光子が、「直接放射線」の大部分、又は直接放射線の経路から微小角度を通り散乱された散乱放射線を表すこと、及び他の画素によって検出された全てのＸ線光子が、入射したペンシルビームから来た散乱Ｘ線光子であることを確信できる。従って、特定の厚さのあらゆる物質の特徴的な散乱パターンを確立できる。この処理を、同じ物質の異なる厚さに対して繰返すことによって、同じ物質の異なる厚さに対する特徴的なパターンは、確立されることができる。同様に、異なる物質及び／又は該物質の異なる厚さに対する処理を繰返して、特徴的な散乱パターンに関するデータベースは、構築されることができる。
【００７８】
Ｘ線システム及びＸ線物理学のモンテカルロモデルは、粒子が被験物を通過するのをシミュレートするＧＥＡＮＴ４と呼ばれるソフトウェアパッケージを使用して、作成された。このモデルは、異なる物質及び該物質の厚さに入射するＸ線光子のペンシルビームの結果をシミュレートするために使用された。異なる物質及び該物質の異なる厚さに対する結果は、記録される。
【００７９】
各散乱パターンは、図１０でグラフ１に示されたものと同様の形を有する。散乱と関連して、これは、散乱カーネルとしても知られている。いかなる散乱カーネルに関しても、ペンシルビームと位置合せされない検出器３の画素に入射した全てのＸ線光子によって表される散乱光子は、識別され、画像から除去されることができる。更に画像を改善するために、除去された散乱放射線は、直接放射線としてペンシルビームと直接位置合せされた画素（複数可）で記録された光子強度に再び加えられてもよい。従って、画素（複数可）で結果として得られる出力強度は、直接放射線の出力信号と、散乱されなければ除去された散乱放射線のＸ線光子が相互作用したであろう検出器の空間的位置での出力信号との合計となる。
【００８０】
しかしながら、対象となる全ての物質及び該物質の厚さに特徴的な散乱カーネルに関するデータベースを作成して、使用するには、膨大なデータ保存容量と極めて強力な処理能力を必要とするであろう、或いは、かかるデータベースに保存されたデータを使用するには、極めて低速になるであろう。
【００８１】
そのため、本発明のこの実施形態のモデルは、それら自体が特にシミュレートされていない物質及び該物質の厚さの組合せに対する散乱カーネルを導出するために、補間法を使用する。補間の技法、及びこの実施形態においてエミュレーションは、それ自体新しいものではない。
【００８２】
図９は、本発明の処理が、どのようにしてＸ線検出器３によって記録されたデータから散乱Ｘ線光子を除去し、疑似直接放射線として散乱光子を再割当するかを示すフローチャートである。
【００８３】
図１０で説明されたフローチャートは、図９で示された処理について、より詳細に説明しており、特に、散乱のない画像を齎すサブステップが説明されている。    
【００８４】
ステップ２０では、エネルギ強度の形のデータが、検出器３の各画素に関して、検出器３によって記録される。
【００８５】
ステップ２１では、物質の種類及び物質の厚さを識別するサブプロセスが、例えば、特許文献３に記載された技術を使用して実行されるが、該技術では、ピクセル化検出器の出力強度が、該出力強度を物質の素性（ｉｄｅｎｔｉｔｙ）及び厚さと相関させる値に関するデータベース３０と比較される。従って、このステップの出力は、検出器３の各画素と関連付けられる物質識別データベースからの物質及び厚さである。用語「物質」は、物質の組合せを含む。例えば、このステップの出力は、検出器３の特定の画素で識別された物質が、筋肉と骨というようにしてもよい。事実上、このステップ２１の出力は、物質の素性及び／又は厚さに関する第１推定値である。
【００８６】
ステップ２２では、識別された物質の散乱カーネルが、ステップ２２ａでエミュレートされる、即ち、散乱カーネルに関するデータベース３１から、散乱カーネルが、検出器３の各画素に関して、ステップ２１で識別された物質に対してエミュレート（導出）される。個々の散乱カーネルは、ステップ２２ｂで纏めてコンボリューションされる。この散乱カーネルのコンボリューションは、画像の散乱推定値として知られるものを形成しており、図１０のグラフ２に示されている。グラフ２から分かるように、散乱推定値には、検出器の各画素と関連付けられる散乱値が存在する。ステップ２２ｃでは、画像を表す検出されたＸ線強度に関する散乱は、検出器の各画素に関する検出器の出力値から散乱強度値を減算することによって、除去される。ステップ２２ｃの出力は、画像としてもよい、又は画像として表され得る強度値としてもよい散乱のない画像２２ｄとして示される。
【００８７】
Ｘ線画像における改善は、処理に関してこの時点で止めることができる。しかしながら、物質の素性及び厚さの第１推定値が正しい可能性が低く、そのため、物質の識別及び／又は厚さの第１推定値に基づいて画像から散乱放射線を除去だけする以上に、画像を改善することが好ましい。
【００８８】
ステップ２２では、初期校正が、物質及び物質の厚さに関するデータベースへの検出器の出力を校正するために実行される。
【００８９】
ステップ２３では、物質識別ステップが、ステップ２２の結果得られた検出器の各画素に関する出力強度に対して実行される（これらの出力強度は、画像を形成する）。この物質識別ステップでは、ステップ２２からの検出器の各画素に関する出力強度は、事前に識別された物質及び／又は厚さについて確立された検出器の出力強度に関するデータベースと比較される。
【００９０】
ステップ２４では、ステップ２３における物質及び／又は厚さの識別が最適化されているか否かについて判定が行われる。これは、まず、ステップ２３における物質及び／又は厚さを表す値（又は、該値と相関するパラメータ）を、ステップ２１における物質及び／又は厚さの識別を表わす値と比較し、ステップ２３の結果がステップ２１における値の閾値内にあるか否かを確立することによって、行われる。閾値内にない可能性が高い、答えがＮＯの場合、ステップ２２及び２３は、現在繰返し中のステップ２２及び２３において物質及び／又は厚さを表す値が、収束するまで、即ち、前の繰返しにおける、又は前の繰返し数の平均における物質及び／又は厚さを表す値の閾値内になるまで、繰返される。答えがＹＥＳの場合、処理はステップ２５に進み、ステップ２５では、散乱は、再割当された散乱を有する画像を作成するために、再割当される、又は再割当されず、散乱のない画像が作成される。
【００９１】
各繰返しに関して、散乱除去ステップ２２は、ステップ２０からの散乱を有するＸ線画像について実行されるが、各繰返しでは、物質及び／又は厚さの識別は、最初の繰返し後の全繰返しに関して、前の繰返しのステップ２３からとなる。散乱除去ステップ後の、ステップ２３からの物質及び／又は厚さの識別を利用することによって、次の散乱除去ステップ２２では、ステップ２３の出力は、繰返しの数が増えるほど繰返し当たりの改善は減少しつつ、通常繰返しの数だけ改善するであろう。第２及びそれ以降の繰返しにおける物質の識別は、ステップ２１における物質及び／又は厚さの識別に使用される第１データベースで、或いは物質の種類及び／又は厚さを表す値に関する更なるデータベースで、実行されてもよい。かかる更なるデータベース又は第１データベース自体は、散乱を除外して、実質的に収集されたデータを使用して、追加されてもよい。ＭＡＰは、ステップ２１における物質の素性及び／又は厚さの第１推定値に必要なだけである。
【００９２】
物質の識別及び／又は厚さを表す値は、コントラスト、散乱カーネル、散乱推定値を含んでもよい。
【００９３】
散乱の影響を軽減しながら、ステップ２３における物質の識別は、収束後の散乱のない画像２４’によって表される、より正確な結果を出すことができる。
【００９４】
ステップ２１〜２４は、物質の種類及び厚さの識別が収束するまで、繰返される。勿論、この比較は、以前の物質値の数の平均と行われることもある。しかしながら、動作原理は、明確である。ステップ２１〜２４は、識別された物質の種類及び厚さ（又は、それらと相関されるパラメータ）に関する大きな変化が全くなくなるまで、再び繰返される。
【００９５】
識別された物質及び／又は厚さに関する大きな変化が全く無くなると、処理はステップ２５に進み、ステップ２５では、出力は、除去された散乱放射線のＸ線光子が、散乱しなければ相互作用したであろう検出器の空間的位置での出力信号に、除去された散乱放射線を追加することによって、散乱が再割当された状態の画像となる。再割当ステップ２５の結果は、ステップ２５’における再割当された画像によって表される。
【００９６】
図１０ａは、別の処理について説明しており、より少ないステップを使用するが、より大きな処理能力と、より大きな散乱カーネルに関するデータベース３１’を必要とする。幾つかの物質及び／又は厚さに関する散乱カーネルデータベースを作成し、使用する代わりに、散乱カーネルデータベース３１’は、データベース３１’が対象となる全ての物質及び／又は厚さの散乱カーネルを含むと考えられるような、多数の物質及び／又は厚さが追加される。従って、本実施形態におけるエミュレーションによる補間ステップ２２ａは、必要なく、散乱の推定値は、ステップ２２ｂで、データベース３１’の散乱カーネル値から直接作成される。
【００９７】
ステップ２２ｃ及びステップ２２ｄは、図１０を参照して説明された通りであり、処理の出力が、ステップ２２ｄからの散乱のない画像となっている。
【００９８】
しかしながら、より少ないステップを必要とするものの、必要とされるデータ採取、データ保存、及び処理能力は、今のところ、図９及び図１０に示された処理が好ましいことを意味している。
【００９９】
図１１ａ〜図１２ｂは、図８ａ〜図１０で説明された本発明の処理を実行する効果について、説明している。
【０１００】
ＧＥＡＮＴ４でモンテカルロＸ線シミュレーションを使用するシミュレーション実験では、かなりの散乱放射線を発生するであろう重複するアルミニウムとプラスチック（ＰＭＭＡ）片のシミュレーション試料が、本技法の有効性を説明するのに使用された。同シミュレーションでは、画像の３領域−Ｒ１、Ｒ２及びＲ３を画定する３試料を使用した。Ｒ１は、ＰＭＭＡとして知られる物質であり、厚さは、１８．５ｍｍであった。Ｒ２は、１８．５ｍｍのＰＭＭＡ及び１２．７ｍｍのアルミニウムを使用した。Ｒ３は、１８．５ｍｍのＰＭＭＡ及び６．３５ｍｍのアルミニウムを使用した。
【０１０１】
散乱除去前に得られる画像は、隣接する領域間、Ｒ１とＲ２間でコントラスト対ノイズ比（ＣＮＲ：ｃｏｎｔｒａｓｔ  ｔｏ  ｎｏｉｓｅ）が１６．１という、低コントラストを示している。
【０１０２】
図９及び図１０で説明された処理は、Ｘ線光子が散乱しなければ相互作用したであろう画素位置に散乱放射線を再割当するために使用された。図１１ａに示された画像と図１１ｂに示された画像とのコントラストに関する視覚的な改善は、特にＲ１とＲ２との間で、明らかである。信号対ノイズ比は、５８．７（非修正の画像より４倍近く良好）であった。また、図１２ｂのヒストグラムは、改善を示している。Ｒ１、Ｒ２及びＲ３に関するピークは、ここでは、よく離間している。
【０１０３】
また、図１２ｃも、本発明によるＸ線信号処理の結果、コントラスト対ノイズ比が改善したことを示している。グラフでは、Ｒ１とＲ２の間及びＲ１とＲ３の間のコントラスト対ノイズ比を示している。これは、標準測度であり、２領域の標準偏差間の差を求め、これを１領域の標準偏差で割って、算出される。
【０１０４】
本発明は、散乱放射線を、より良好に識別でき、従って、直接放射線と散乱放射線を含む画像から一層容易に除去できる。
【０１０５】
更に、本発明は、コントラスト対ノイズ比を、散乱放射線を除去することによって、又は散乱放射線を除去し、除去された散乱放射線のＸ線光子が散乱されなければ相互作用したであろう検出器の空間的位置での出力信号に、除去された散乱放射線を再び加えることによって、改善できる。
【０１０６】
コントラスト対ノイズ比を改善することによって、Ｘ線投与量は、減少されるかも知れない、又は試料又は患者へのＸ線投与量を増大させずにコントラスト対ノイズ比を高めるかも知れない。医療においては、投与量を減少することは、患者に害を与える可能性を低くし、そのため極めて望ましいことは、よく理解されている。
【０１０７】
従って、本発明の装置及び方法は、医療画像処理に使用されるＸ線量を減少できる。
【０１０８】
異なる実施形態の特徴は、該特徴が記載された実施形態だけに限らず、本明細書に記載された他の実施形態で使用されてもよい。
【図１】
0006598087.tif 000002
【図２】
0006598087.tif 000003
【図３】
0006598087.tif 000004
【図４】
0006598087.tif 000005
【図５】
0006598087.tif 000006
【図６】
0006598087.tif 000007
【図７】
0006598087.tif 000008
【図８ａ】
0006598087.tif 000009
【図８ｂ】
0006598087.tif 000010
【図９】
0006598087.tif 000011
【図１０】
0006598087.tif 000012
【図１０ａ】
0006598087.tif 000013
【図１１ａ】
0006598087.tif 000014
【図１１ｂ】
0006598087.tif 000015
【図１２ａ】
0006598087.tif 000016
【図１２ｂ】
0006598087.tif 000017
【図１２ｃ】
0006598087.tif 000018
【図１３】
0006598087.tif 000019