(57)【要約】  本発明は、Ｘ線イメージングデータ、特に位相コントラストＸ線イメージングデータの処理に関し、例えば、スペクトル情報を利用することによる格子ベースの微分位相コントラストイメージングにおける振動に関連する問題に対処する。マルチエネルギービンＸ線検出器、振動誘起モアレパターンの位相シフトのエネルギー非依存性、物体誘起モアレパターンの位相シフトのエネルギー依存性、及びエネルギービンの関数としての物体誘起位相シフトの依存性を測定できるように設計されたファントムに基づいて、モアレパターンの振動状態を決定する際にスペクトル情報を利用することができることが見出された。特に、撮像装置の基本的な物理構造を変更する必要がない一方で、適切なデータ処理により、振動誘起シフト及び同様のものが、物体の特性から分離されることができ、画像内の振動アーチファクトが低減されることがわかった。
2020512080.tif 000002
【特許請求の範囲】
【請求項１】
エネルギービニングされる微分位相コントラストＸ線イメージングデータ用の処理装置であって、
第１のＸ線エネルギービンに関する物体誘起位相シフトと第２のＸ線エネルギービンに関する物体誘起位相シフトとの間の比例係数を提供する比例係数提供ユニットと、
前記比例係数提供ユニットによって提供される前記比例係数を使用して、前記イメージングデータに含まれる、前記第１のエネルギービンのピクセルの検出強度値及び前記第２のエネルギービンのピクセルの検出強度値から、暗視野信号及び物体誘起位相シフトの少なくとも一方を決定する暗視野及び位相シフト決定ユニットと、
を有する処理装置。
【請求項２】
前記比例係数提供ユニットは、前記第１のエネルギービンの減衰値及び前記第２のエネルギービンの減衰値に依存して、前記比例係数を提供する、請求項１に記載の処理装置。
【請求項３】
前記比例係数提供ユニットが、ルックアップテーブルを有する、請求項２に記載の処理装置。
【請求項４】
前記暗視野及び位相シフト決定ユニットは、前記測定データを取得する見込みが最大になるように、前記第１のエネルギービンの測定データ及び前記第２のエネルギービンの測定データから、前記第１のエネルギービンの減衰値、前記第２のエネルギービンの減衰値、前記第１のエネルギービンの暗視野信号、前記第２のエネルギービンの暗視野信号、物体誘起位相シフト、及び外乱誘起位相シフト、を決定する、請求項１に記載の処理装置。
【請求項５】
前記暗視野及び位相シフト決定ユニットは、
2020512080.tif 000013
に従う最小二乗法による最小化を用い、ここで、Ｂｉ，１、Ｂｉ，２、Ｖｉ，１、Ｖｉ，２、αｉ，１及びαｉ，２は、それぞれ、ブランクスキャン中に取得されるモアレパターンの強度、可視性、及び位相である、請求項１に記載の処理装置。
【請求項６】
前記暗視野及び位相シフト決定ユニットは、前記決定に、隣接ピクセルの外乱誘起位相シフトの間の小さな差を優先する更なる制限を含めるように構成される、請求項１に記載の処理装置。
【請求項７】
エネルギービニングされる微分位相コントラストＸ線イメージングにおいて、第１のＸ線エネルギービンに関する物体誘起位相シフトと第２のＸ線エネルギービンに関する物体誘起位相シフトとの間の比例係数を決定する比例係数決定装置であって、
予め決められた形状を有するファントムを使用して、第１のエネルギービンを使用する微分位相コントラストＸ線イメージングにより得られる位相画像から前記第１のエネルギービンの第１の位相シフトを計算し、前記ファントムを使用して、第２のエネルギービンを使用する微分位相コントラストＸ線イメージングにより得られる位相画像から前記第２のエネルギービンの第２の位相シフトを計算する計算ユニットと、
前記第１の位相シフト及び前記第２の位相シフトから前記比例係数を得る取得ユニットと、
を有する、比例係数決定装置。
【請求項８】
エネルギービニングされる微分位相コントラストＸ線イメージングにおいて、第１のＸ線エネルギービンに関する物体誘起位相シフトと第２のＸ線エネルギービンに関する物体誘起位相シフトとの間の比例係数を決定する比例係数決定システムであって、
請求項７に記載の比例係数決定装置と、
複数の部分を有するファントムであって、各部分は、Ｘ線のビーム方向を横切る方向において予め決められた態様で変化する高さを有し、前記ファントムの前記部分の高さは、前記ビーム方向における前記部分の延長であり、前記ファントムのそれぞれ異なる部分が、それぞれ異なる平均高さを有する、ファントムと、
を有し、
前記比例係数決定装置が、前記第１のエネルギービンの減衰値及び前記第２のエネルギービンの減衰値に依存して、前記比例係数を決定する、比例係数決定システム。
【請求項９】
微分位相コントラストＸ線イメージングシステムであって、
微分位相コントラストイメージングユニットと、
第１及び第２のエネルギービンのＸ線を提供するＸ線源と、
請求項１に記載の処理装置であって、前記微分位相コントラストイメージングユニットから、エネルギービニングされる微分位相コントラストＸ線イメージングデータを受信するように構成される処理装置と、
を有する微分位相コントラストＸ線イメージングシステム。
【請求項１０】
請求項８に記載の比例係数決定システムを更に有し、前記比例係数決定システムは、前記処理装置の前記比例係数提供ユニットに、前記比例係数に関する情報を提供する、請求項９に記載の微分位相コントラストＸ線イメージングシステム。
【請求項１１】
エネルギービニングされる微分位相コントラストＸ線イメージングデータの処理方法であって、
第１のＸ線エネルギービンに関する物体誘起位相シフトと第２のＸ線エネルギービンに関する物体誘起位相シフトとの間の比例係数を提供する比例係数提供ステップと、
前記比例係数提供ユニットによって提供される比例係数を使用して、前記イメージングデータに含まれる前記第１のエネルギービンのピクセルの検出強度値及び前記第２のエネルギービンのピクセルの検出強度値から、暗視野信号及び物体誘起位相シフトの少なくとも一方を決定する決定ステップと、
を有する方法。
【請求項１２】
前記比例係数提供ステップは、
前記第１のＸ線エネルギービンと予め決められた形状を有するファントムとを使用して、微分位相コントラストＸ線イメージングにより得られる位相画像から、前記第１のＸ線エネルギービンの第１の位相シフトを計算し、前記第２のＸ線エネルギービンと前記ファントムとを使用して、微分位相コントラストＸ線イメージングにより得られる位相画像から、前記第２のＸ線エネルギービンの第２の位相シフトを計算する計算ステップと、
前記第１の位相シフト及び前記第２の位相シフトから前記比例係数を得る取得ステップと、
を有する、請求項１１に記載の方法。
【請求項１３】
エネルギービニングされる微分位相コントラストＸ線イメージングデータの処理方法のためのソフトウェアであって、前記ソフトウェアが前記コンピュータ上で実行されるとき、請求項１１に記載の方法の各ステップをコンピュータに実行させるプログラムコード手段を有するソフトウェア。
【請求項１４】
エネルギービニングされる微分位相コントラストＸ線イメージングにおいて、第１のＸ線エネルギービンに関する物体誘起位相シフトと第２のＸ線エネルギービンに関する物体誘起位相シフトとの間の比例係数を決定するためのソフトウェアであって、前記ソフトウェアは、前記ソフトウェアがコンピュータ上で実行されるとき、請求項１２に記載の方法の各ステップを前記コンピュータに実行させるプログラムコード手段を有するソフトウェア。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、Ｘ線イメージングデータ、特に位相コントラストＸ線イメージングデータの処理に関する。
【背景技術】
【０００２】
回折格子ベースの微分位相コントラストイメージングは、近年勢いを増している。そのような技術は、従来の減衰画像、すなわち屈折及び暗視野信号を補完する追加情報を提供することができる。とりわけ、フルフィールドデジタルマンモグラフィ（ＦＦＤＭ）は、この技術の有望なアプリケーションである。最近、最初の臨床ＦＦＤＭユニットが、Philips MicroDose L30に基づくＸ線回折格子ハードウェアを装備することに成功した(T. Koehler et al., Slit-scanning differential x-ray phase-contrast mammography: Proof-of-concept experimental studies, Medical Physics 42, 1959 (2015)参照）。最初の結果は、スリットスキャンシステムがＦＯＶ、スキャン時間、線量に関する重要な要件を満たすことを示す。
【発明の概要】
【発明が解決しようとする課題】
【０００３】
しかしながら、微分位相コントラストイメージングは、機械的な不安定性及び振動の影響を受けやすい。物体のパラメータ（減衰、屈折、暗視野信号）を抽出するために、記録されたデータは通常、物体なしの参照スキャン（「ブランクスキャン」、上記のT. Koehler他の文献を参照）を考慮することにより処理される。その結果、結果的に得られる画像品質は、反復スキャン中に得られるモアレパターンの繰り返し性によって制限される。従来、例えばブランクスキャンと物体スキャンとの間の振動によって引き起こされる偏差は、誤って物体に帰するものとされて、画像アーチファクトをもたらす。
【０００４】
例えば、独国特許出願公開第102014213817A1公報に記述されるように、イメージング装置の構造的な変更により振動等の発生を低減することを目的とする試みがある。しかしながら、このような構造的な変更は、より複雑な機械、すなわちコストの増加につながり、洗練された試みであっても、望ましくない影響を完全に排除することはできない。
【０００５】
本発明の目的は、特に格子ベースの微分位相コントラストイメージングにおける振動及び同様のものに関連する前述の問題に対処するためのアプローチを提供することである。
【課題を解決するための手段】
【０００６】
本発明の第１の態様において、エネルギービニングされる微分位相コントラストＸ線イメージングデータの処理装置であって、第１のＸ線エネルギービンに関する物体誘起位相シフト（object induced phase shift）と第２のＸ線エネルギービンに関する物体誘起位相シフトとの間の比例係数を提供するように構成される比例係数提供ユニットと、比例係数提供ユニットによって提供される比例係数を使用して、前記イメージングデータに含まれる第１のエネルギービンのピクセルの検出強度値及び第２のエネルギービンのピクセルの検出強度値から、暗視野信号及び物体誘起位相シフトの少なくとも一方を決定するように構成される暗視野及び位相シフト決定ユニットと、を有する処理装置が提示される。請求項に記載の本発明の状況において、「暗視野及び位相シフト決定ユニット」の語は、暗視野信号及び物体誘起位相シフトの両方を決定するように構成される単一のユニット、又は暗視野信号と物体誘起位相シフトとを個別に別々のサブユニットにおいて決定するように構成される決定ユニットを意味することができる。
【０００７】
本発明の第２の態様において、エネルギービニングされる微分位相コントラストＸ線イメージングにおいて、第１のＸ線エネルギービンに関する物体誘起位相シフトと第２のＸ線エネルギービンに関する物体誘起位相シフトとの間の比例係数を決定する比例係数決定装置であって、第１のエネルギービン及び予め決められた形状をもつファントムを使用して、微分位相コントラストＸ線イメージングにより得られる位相画像から第１のエネルギービンに関する第１の位相シフトを計算し、第２のエネルギービン及び前記ファントムを使用して、微分位相コントラストＸ線イメージングにより得られる位相画像から第２のエネルギービンの第２の位相シフトを計算するように構成される計算ユニットと、第１の位相シフト及び第２の位相シフトから比例係数を取得するように構成される取得ユニットと、を有する比例係数決定装置が提示される。
【０００８】
本発明の第３の態様において、エネルギービニングされる微分位相コントラストＸ線メージングデータの処理方法であって、第１のＸ線エネルギービンに関する物体誘起位相シフトと第２のＸ線エネルギービンに関する物体誘起位相シフトとの間の比例係数を提供する比例係数提供ステップと、比例係数提供ステップにより提供される比例係数を使用して、前記イメージングデータに含まれる第１のエネルギービンのピクセルの検出強度値及び第２のエネルギービンのピクセルの検出強度値から暗視野信号及び物体誘起位相シフトの少なくとも一方を決定する決定ステップと、を有する処理方法が提示される。
【０００９】
本発明の第４の態様では、エネルギービニングされる位相差コントラストＸ線イメージングにおいて第１のＸ線エネルギービンの物体誘起位相シフトと第２のＸ線エネルギービンの物体誘起位相シフトとの間の比例係数を決定する比例係数決定方法であって、第１のエネルギービン及び予め決められた形状を有するファントムを使用して、微分位相コントラストＸ線イメージングにより得られる位相画像から第１エネルギービンの第１位相シフトを計算し、第２のエネルギービン及び前記ファントムを使用して、微分位相コントラストＸ線イメージングにより得られる位相画像から第２のエネルギービンの第２の位相シフトを計算する計算ステップと、第１の位相シフト及び第２の位相シフトから比例係数を取得する取得ステップと、を有する比例係数決定方法が提示される。
【００１０】
本発明は、例えば、スペクトル情報を利用することにより、格子ベースの微分位相コントラストイメージングにおける振動に関連する問題に対処する。
【００１１】
本願の発明者は、マルチエネルギービンＸ線検出器、振動誘起モアレパターン位相シフトのエネルギー独立性、物体誘起モアレパターン位相シフトのエネルギー依存性、及びエネルギービンの関数として物体誘起位相シフトの依存性を測定できるように設計されたファントム、に基づいて、モアレパターンの振動状態を決定する際にスペクトル情報を利用することができることを見出した。
【００１２】
特に、イメージング装置の基本的な物理構造を変更する必要がない一方で、適切なデータ処理により、振動誘起シフト及び同様のものが、物体の特性から分離されることができ、ゆえに、画像内の振動アーチファクトが低減されることがわかった。
【００１３】
ここで、本発明は、２より多くのエネルギービンも企図されるので、２つのエネルギービンのみの場合に限定されないことに留意されたい。それぞれの検出された強度値を考慮する場合、振動による位相シフトを決定するために３又はそれより多くのエネルギービンが使用されることができ、更なるエネルギービンの使用は、例えば、本願明細書において２つのエネルギービンについて明示的に説明されるように、第１及び第２のエネルギービンを考慮すること、第１及び（更に）第３のエネルギービンを考慮すること、並びに第２及び第３のエネルギービンを考慮することのそれぞれを実施することによって、個別の（ある程度独立した）決定が可能になるという利点を与える。更に、２つのエネルギービンを使用して、振動による位相シフトを識別することができ、これは、第３のエネルギービンを使用した同時微分位相コントラストイメージングの補正に使用される。
【００１４】
好適な実施形態において、比例係数提供ユニットは、第１のエネルギービンの減衰値及び第２のエネルギービンの減衰値に応じて比例係数を提供するように構成される。
【００１５】
第１及び第２のエネルギービンそれぞれの減衰値に応じて比例係数を提供する場合、精度の向上が実現可能であるが、状況に応じて、（例えば影響がほとんどない場合に）個々の異なる減衰について変化しない係数の近似値を提供することもできる。
【００１６】
上述の実施形態の好適な変形例において、比例係数提供ユニットはルックアップテーブルを有する。
【００１７】
ルックアップテーブルは、減衰値をルックアップテーブルのパラメータとして使用することにより、比例係数に迅速にアクセスすることを可能にする。それでもなお、係数を提供する他の方法も可能であり、例えば計算式、又は計算及びルックアップテーブルの組み合わせを提供することも可能である。
【００１８】
好適な実施形態において、暗視野信号及び位相シフト決定ユニットは、選択された画像値に関して測定データを得る見込みが最大化されるように、第１のエネルギービンの測定されたデータｍｉ，１及び第２のエネルギービンの測定されたデータｍｉ，２から、第１のエネルギービンの減衰値Ａ１、第２のエネルギービンの減衰値Ａ２、第１のエネルギービンの暗視野信号Ｄ１、第２のエネルギービンの暗視野信号Ｄ２、物体誘起位相シフトφ１、及び外乱誘起位相シフトψｉを決定するように構成される。
【００１９】
これは、例えば、測定データと期待データとの間の差にペナルティを課す費用関数が最小化されることを含む。
【００２０】
上記の実施形態の好適な変形例において、暗視野信号及び位相シフト決定ユニットは、
2020512080.tif 000003
に従う最小二乗法による最小化を用いるように構成される。ここで、Ｂｉ，１、Ｂｉ，２、Ｖｉ，１、Ｖｉ，２、αｉ，１及びαｉ，２は、それぞれ、ブランクスキャン中に取得されるモアレパターンの強度、可視性、及び位相である。
【００２１】
しかしながら、最小二乗最小化のようなアプローチを使用することは、１つの選択肢にすぎない。与えられた状況のように、光子のカウントを扱い、負のポアソン対数尤度を最小化することができ、ブランクスキャンの強度が変化する場合は、重み付き最小二乗費用関数を使用することが望ましいことがある。
【００２２】
具体的には、上述の実施形態の好適な変形例としても考えられうる更なる好適な実施形態において、暗視野信号及び位相シフト決定ユニットは、決定（例えば、最小二乗最小化）の際に、隣接ピクセルの外乱誘起位相シフトの間のわずかな差を優先するために更なる制限を含むように構成される。特に、最小二乗最小化の場合、これは、予め決められた値より小さい隣接ピクセルの外乱誘起位相シフトの間の絶対差に関する制限に対応しうる。
【００２３】
隣接ピクセルが同様の外乱又は振動に苦しむことは先験的知識と見なすことができ、従って、そのような情報が、決定に含められることができる。振動は、回折格子の互いに対する相対位置を変化させることがあり、この場合、更に干渉計は、そのような変化に非常に敏感であることに留意されたいす（μｍレンジの変化でさえ目に見える効果をもたらすことがわかった）。
【００２４】
好適な実施形態において、エネルギービニングされる微分位相コントラストＸ線イメージングにおいて、第１のＸ線エネルギービンの物体誘起位相シフトと第２のＸ線エネルギービンの物体誘起位相シフトとの間の比例係数を決定する比例係数決定システムが提供され、比例係数決定システムは、本発明による比例係数決定装置と、複数の部分を有するファントムとを有し、ファントムの各部分は、Ｘ線のビーム方向を横切る方向において、予め決められた態様で変化する高さを有し、ファントムの部分の高さは、ビーム方向における前記部分の延長であり、ファントムのそれぞれ異なる部分は、異なる平均高さを有し、比例係数決定装置は、第１のエネルギービンの減衰値及び第２のエネルギービンの減衰値に応じて比例係数を決定するように構成される。
【００２５】
好適な実施形態において、微分位相コントラストＸ線イメージングシステムが提供され、微分位相コントラストＸ線イメージングシステムは、微分位相コントラストイメージングユニットと、第１及び第２のエネルギービンのＸ線を提供するように構成されるＸ線源と、微分位相コントラストイメージングからエネルギービニングされる微分位相コントラストＸ線イメージングデータを受信するように構成される本発明による処理装置と、を有する。
【００２６】
上述の実施形態の好適な変形例において、イメージングシステムは更に、上述の実施形態による比例係数決定システムを有し、比例係数決定システムは、処理装置の比例係数提供ユニットに比例係数に関する情報を提供するように構成される。
【００２７】
比例係数が一定であり（又は定数の値（又は値の組）によって少なくとも十分に近似され）、ゆえに、値を決定するための特定の構成が必要とされない（予め決められた値又は値のセットを記憶する手段を除く）場合が企図され、他の場合には、比例係数に関する情報を明確に決定するように構成される比例係数決定システムが有利であり、これは、例えばエネルギービンのサイズが非常に大きいため、これらのビンにより測定される有効エネルギーが、物体によるビームハードニングのために大きく変化する場合などである。しかしながら、エネルギービンが狭い場合は、単純に式（５）（以下を参照）を使用して比例係数を計算するだけで十分でありうる。
【００２８】
本発明の更なる態様において、エネルギービニングされる微分位相コントラストＸ線イメージングデータの処理方法のためのソフトウェア製品であって、ソフトウェア製品がコンピュータ上で実行されるときに本発明による方法のステップをコンピュータに実行させるプログラムコード手段を有するソフトウェア製品、及びエネルギービニングされる微分位相コントラストＸ線イメージングにおいて、第１のＸ線エネルギービンの物体誘起位相シフトと第２のＸ線エネルギービンの物体誘起位相シフトとの間の比例係数を決定するソフトウェア製品であって、ソフトウェア製品がコンピュータ上で実行されるときに本発明による方法のステップをコンピュータに実行させるプログラムコード手段を有するソフトウェア製品が提示される。
【００２９】
請求項１に記載の処理装置、請求項７に記載の決定ユニット、請求項１１に記載の処理方法、請求項１２に記載の決定方法、並びに請求項１３及び１４に記載のコンピュータプログラムは、特に従属請求項に規定される同様の及び／又は同一の好適な実施形態を有することを理解されたい。
【００３０】
本発明の好適な実施形態は、従属請求項又は上記の実施形態とそれぞれの独立請求項との任意の組み合わせでもあり得ることを理解されたい。
【００３１】
本発明のこれら及び他の態様は、以下に記述される実施形態から明らかになり、それらを参照して説明される。
【図面の簡単な説明】
【００３２】
【図１】本発明の実施形態による微分位相コントラストＸ線イメージングシステムを示す図。
【図２】本発明の実施形態による、エネルギービニングされる微分位相コントラストＸ線イメージングデータの処理方法を示すフロー図。
【発明を実施するための形態】
【００３３】
図１は、本発明の実施形態による微分位相コントラストＸ線イメージングシステムを示す。
【００３４】
システム１０は、例えば第１のエネルギーレンジの光子及び第２のエネルギーレンジの光子を含む第１及び第２のエネルギービンのＸ線を提供するように構成されるＸ線源１５と、対応する微分位相コントラストイメージングユニット２０と、を有する。
【００３５】
図１では、Ｘ線源１５とイメージングユニット２０との間にファントム２５が設けられているが、対象物を撮像するシステム１０の動作時は、ファントム２５ではなく対象物（図示せず）が、Ｘ線源１５とイメージングユニット２０との間に提供される。
【００３６】
ファントム２５は、それぞれ異なる高さのいくつかの柱体２５ａからなる。これらのピラー２５ａの上部には、ビーム方向において投影される材料高さの一定の勾配をもつ三角形がある。投影される電子密度の一定の勾配（ビーム方向の材料高さによって実現される）は、モアレパターンの一定の位相シフトを生成することに留意されたい。
【００３７】
更に、システム１０において、イメージングユニット２０は、比例係数決定装置３０及び処理装置３５に結合される。
【００３８】
比例係数決定装置３０は、計算ユニット４０及び取得ユニット４５を有し、それらの動作については以下で説明する。
【００３９】
処理ユニット３５は、比例係数提供ユニット５０と、暗視野信号及び位相シフト決定ユニット５５とを有し、その動作も以下の記述から明らかになる。
【００４０】
比例係数提供ユニット５０は、比例係数決定装置３０に結合され、そこから比例係数に関するデータ又は情報を受信する。
【００４１】
図２は、本発明の実施形態による、エネルギービニングされる微分位相コントラストＸ線イメージングデータの処理方法を示すフロー図を示す。
【００４２】
以下で更に説明するように、この実施形態の方法は、計算ステップ６０を有し、図１に示されるようにファントムを微分位相コントラストＸ線イメージングすることにより得られる位相画像から、第１のエネルギービンに関する第１の位相シフトが計算され、かかる位相画像は、第１のエネルギービンを使用し又はそれに対応する。計算ステップ６０は更に、第２のエネルギービンを使用してファントムを微分位相コントラストＸ線イメージングすることにより得られる位相画像から、第２のエネルギービンの第２の位相シフトを計算することを含む。
【００４３】
計算ステップ６０において計算された情報に基づいて、次の取得ステップ６５は、減衰に関する影響を考慮しながら、ファントムに基づいて第１の位相シフト及び第２の位相シフトから比例係数を取得することを含む。
【００４４】
得られる比例係数は、比例係数提供ステップ７０において提供され、第１のＸ線エネルギービンに関する物体誘起位相シフトと第２のＸ線エネルギービンに関する物体誘起位相シフトとの間の比例関係を示す。
【００４５】
最後に、比例係数提供ステップ７０により提供される比例係数を使用して、前記イメージングデータに含まれる第１のエネルギービンのピクセルの検出強度値及び第２のエネルギービンのピクセルの検出強度値から、暗視野信号及び物体誘起位相シフトの少なくとも一方を決定する決定ステップ７５が提供される。
【００４６】
以下において、本発明の実施形態の上述の概要が更に説明される。
【００４７】
一般的な幾何学的光線について、ピクセルｉのタルボット−ロー干渉計の背後で測定される強度について一般に使用されるモデルは、
2020512080.tif 000004
【００４８】
ここで、Ｂｉ、Ｖｉ、及びαｉは、ブランクスキャン中に得られるモアレパターンの強度、可視性、及び位相をそれぞれ示す。物体のパラメータＡ、Ｄ、φは、それぞれ減衰、暗視野信号、位相シフトに対応する。物体パラメータは、測定データｍｉを使用して最小二乗最小化によって取得される。
2020512080.tif 000005
【００４９】
しかしながら、上記のように、機械的振動は、以下に対応するモアレパターンの追加の位相シフトψｉを生成する。
2020512080.tif 000006
【００５０】
非スペクトルデータ蓄積を用いる場合、φ及びψｉが、最小化の最中に分離されることができない。しかしながら、２つのエネルギービンｅ＝［１，２］を有するＸ線検出器を使用する場合、式（３）は、以下となる：
2020512080.tif 000007
【００５１】
振動シフトψｉを除くすべてのパラメータは光子エネルギーに依存することに注意されたい。
【００５２】
単一エネルギーのＸ線光子の場合、位相シフトφのエネルギー依存性は以下のとおりである：
2020512080.tif 000008
【００５３】
多色Ｘ線の場合、比例係数ｃは、
2020512080.tif 000009
により両方のエネルギービンにおいて得られる位相シフトφ１とφ２との間の関係を記述し、比例係数ｃは、両方のエネルギービンに対応する有効エネルギースペクトルｓ１及びｓ２に依存する。それは、図１に示されるファントムを使用して測定されることができる。
【００５４】
両方のエネルギービンについて別個に式（２）を使用することにより、両方のビンの位相画像が得られる。位相画像内の個々の柱体に対応する領域において平均化することにより、式（６）を介して各柱体のｃが得られる。高さが異なるため、柱体は、個々異なる程度のビームハードニングを与え、かかるビームハードニングは、両方のビンＡ１，Ａ２について減衰信号の中で符号化され、減衰画像から抽出されることができる。この方策に従って、ｃ（Ａ１，Ａ２）についてルックアップテーブルを得ることが可能である。
【００５５】
わずかに異なる実施形態（図示しない）において、エネルギービンの平均エネルギーの商に対応する多色Ｘ線の場合も、単一エネルギーＸ線光子の場合と同様に比例係数ｃが等しくなり、それに対応して、ｃ（Ａ１，Ａ２）についてそのようなルックアップテーブルが得られうる。
【００５６】
物体パラメータから振動誘起シフトを分離するために、両方のエネルギービンに関して式（３）を明示的に記述し、式（６）を使用する。
2020512080.tif 000010
【００５７】
両方のビンからの測定データｍｉ，ｅを使用する最小二乗最小化によって、物体の位相φ１は、振動成分ψｉから分離されることができる：
2020512080.tif 000011
【００５８】
上記の説明は、ピクセルごとに物体誘起位相シフトと振動誘起位相シフトとを分離するのに十分な情報をすでに提供していることに注意されたい。しかしながら、更なる事前知識を追加することによって、推定は更にロバストにされることができ、すなわち、データが隣接する検出器ピクセルと同時に取得される場合、振動誘起位相は、隣接ピクセル同士で同じ（又はほぼ同じ）はずである。スキャン方向に垂直な画像内のラインのピクセルにｘのインデックスが付けられる場合、最小化は次のような式で表現できる：
2020512080.tif 000012
ここでβは、振動誘起位相の変動がどれくらい許容されるかを制御する正則化パラメータである。この場合、隣接するピクセルの振動位相もまた考慮される。太字のシンボルは、ベクトルに言及するために使用される。最後の項は、隣接ピクセルの外乱誘起位相シフトの間の小さな差を優先し、これは、これらの差が振動する格子の共通の原因であるという先験知識と一致する。
【００５９】
本願発明者らは、マルチエネルギービンＸ線検出器に基づくスペクトル情報を利用することにより、機械的不安定性及び変動に伴う問題を回避することを提案する。モアレパターンの振動状態は、振動誘起モアレパターン位相シフトのエネルギー非依存性と、物体誘起シフトのエネルギー依存性を利用して決定されることができる。本願発明者らは、物体誘起位相シフトのエネルギー依存性を測定するためのファントムを提案する。適切なデータ処理により、振動誘起シフトが、物体の特性から分離されることにより、画像のアーチファクトを最小限に抑えることができる。
【００６０】
本発明は、図面及び前述の説明において詳しく図示され説明されたが、そのような図示及び説明は、説明的又は例示的であり、限定的なものではないと考えられるべきである。本発明は開示された実施形態に限定されない。
【００６１】
特に、上記の説明は、２つのエネルギービンの使用に言及しているが、本発明は、そのような方法に限定されると理解されるべきではなく、例えば追加のエネルギービンの場合、それに応じて追加の比例係数が提供されることができることを述べておく。
【００６２】
開示される実施形態に対する他の変更形態は、図面、開示、及び添付の特許請求の範囲の検討から、請求項に記載の本発明を実施する際に当業者によって理解され、達成されることができる。
【００６３】
特許請求の範囲において、「含む、有する（comprising）」という語は、他の構成要素又はステップを除外するものではなく、不定冠詞「a」又は「an」は複数性を除外しない。
【００６４】
単一のプロセッサ、デバイス、又は他のユニットが、特許請求の範囲に記載されるいくつかのアイテムの機能を果たすこともできる。特定の手段が相互に異なる従属請求項に記載されているという単なる事実は、これらの手段の組み合わせが有利に使用されることができないことを示すものではない。
【００６５】
提供、決定、計算、取得、及び処理などの操作は、コンピュータプログラムのプログラムコード手段として、及び／又は専用ハードウェアとして実現されることができる。
【００６６】
コンピュータプログラムは、他のハードウェアと一緒に又はその一部として供給される光学記憶媒体又は固体媒体のような適切な媒体に記憶／配布されることができるが、他の形態で、例えばインターネット又は他の有線若しくは無線の電気通信システムを介して、配布されることもできる。
【００６７】
請求項中の如何なる参照符号も請求項の範囲を制限するものとして解釈されるべきでない。
【図１】
2020512080.tif 000014
【図２】
2020512080.tif 000015
【国際調査報告】
2020512080.tif 000016
2020512080.tif 000017