(57)【要約】  本発明は、画像をフィルタリングするための画像処理装置に関する。当該装置は、同じオブジェクトの第１及び第２画像を取得する画像入力を有し、第１及び第２画像は複数のボクセルを有し、ノイズ共分散によって相関し、各ボクセルは信号値とノイズ値とを含むボクセル値を有する。第１画像をフィルタリングするジョイントバイラテラルフィルタが備えられ、第１画像の第１ボクセルは相対的なボクセルに固有のウェイトを含むフィルタ関数によりフィルタリングされ、第１画像の第１ボクセルの信号値が第１画像の第２ボクセルの信号値と同一であり、第２画像の第１ボクセルの信号値が第２画像の第２ボクセルの信号値と同一であると仮定すると、ウェイトは、第１画像における第１ボクセル及び第２ボクセルのボクセル値と、第２画像における第１ボクセル及び第２ボクセルのボクセル値とを取得する尤度を含む。フィルタリングされた画像は画像出力において提供される。
2016506267.tif 000002
【特許請求の範囲】
【請求項１】
画像をフィルタリングする画像処理装置であって、
同じオブジェクトの第１及び第２画像を取得する画像入力であって、前記第１及び第２画像は複数のボクセルを有し、ノイズ共分散によって相関し、各ボクセルは信号値とノイズ値とを含むボクセル値を有する、画像入力と、
前記第１画像をフィルタリングするジョイントバイラテラルフィルタであって、前記第１画像の第１ボクセルは相対的なボクセルに固有のウェイトを含むフィルタ関数によりフィルタリングされ、前記第１画像の第１ボクセルの信号値が前記第１画像の第２ボクセルの信号値と同一であり、前記第２画像の第１ボクセルの信号値が前記第２画像の第２ボクセルの信号値と同一であると仮定すると、前記ウェイトは、前記第１画像における前記第１ボクセル及び第２ボクセルのボクセル値と、前記第２画像における第１ボクセル及び第２ボクセルのボクセル値とを取得する尤度を含む、ジョイントバイラテラルフィルタと、
前記フィルタリングされた画像を提供する画像出力と、
を有する画像処理装置。
【請求項２】
前記フィルタは、前記ノイズ共分散を記述するノイズ共分散行列を含む尤度を利用するよう構成される、請求項１記載の画像処理装置。
【請求項３】
前記フィルタは、
【数１】
2016506267.tif 000057
又は
【数２】
2016506267.tif 000058
の形式の共分散を利用するよう構成される、請求項２記載の画像処理装置。
【請求項４】
前記フィルタは、
【数３】
2016506267.tif 000059
の形式の尤度を利用するよう構成され、Ｓは前記共分散行列の逆行列であり、ａ１及びａ２は前記第１画像の第１及び第２ボクセル値であり、ｂ１及びｂ２は前記第２画像の第１及び第２ボクセル値であり、Δａ＝ａ１−ａ２，Δｂ＝ｂ１−ｂ２である、請求項２記載の画像処理装置。
【請求項５】
前記フィルタは、
【数４】
2016506267.tif 000060
の形式の尤度を利用するよう構成され、Ｓは前記共分散行列の逆行列であり、ａ１及びａ２は前記第１画像の第１及び第２ボクセル値であり、ｂ１及びｂ２は前記第２画像の第１及び第２ボクセル値であり、Δａ＝ａ１−ａ２，Δｂ＝ｂ１−ｂ２である、請求項２記載の画像処理装置。
【請求項６】
前記第１画像のフィルタリングされた第１ボクセル値は、
【数５】
2016506267.tif 000061
として取得され、ｗｄはバイラテラルフィルタリングにおける距離重み付けである、請求項２記載の画像処理装置。
【請求項７】
前記画像入力は、前記オブジェクトの異なる性質を示すように、前記第１画像と前記第２画像とを提供するよう構成される、請求項１記載の画像処理装置。
【請求項８】
前記画像入力は、スペクトルＣＴ画像を提供するよう構成され、前記第１画像は前記オブジェクトの第１の物理的性質を示し、前記第２画像は前記オブジェクトの第２の物理的性質を示し、前記第１の物理的性質と前記第２の物理的性質とは互いに異なる、請求項７記載の画像処理装置。
【請求項９】
前記第１画像と前記第２画像との一方はコンプトン画像であり、前記第１画像と前記第２画像との他方は光電画像である、請求項７記載の画像処理装置。
【請求項１０】
前記第１画像と前記第２画像との一方はコンプトン画像又は光電画像であり、前記第１画像と前記第２画像との他方はコンプトン画像と光電画像との組み合わせである、請求項７記載の画像処理装置。
【請求項１１】
前記画像入力は、微分位相コントラストＣＴ画像を提供するよう構成され、前記第１画像と前記第２画像との一方は前記オブジェクトの吸収を示し、前記第１画像と前記第２画像との他方は前記オブジェクトにより誘導される前記オブジェクトをトラバースする放射線の位相シフトを示す、請求項１記載の画像処理装置。
【請求項１２】
画像をフィルタリングする画像処理方法であって、
同じオブジェクトの第１及び第２画像を取得するステップであって、前記第１及び第２画像は複数のボクセルを有し、ノイズ共分散によって相関し、各ボクセルは信号値とノイズ値とを含むボクセル値を有する、取得するステップと、
前記第１画像をジョイントバイラテラルフィルタリングするステップであって、前記第１画像の第１ボクセルは相対的なボクセルに固有のウェイトを含むフィルタ関数によりフィルタリングされ、前記第１画像の第１ボクセルの信号値が前記第１画像の第２ボクセルの信号値と同一であり、前記第２画像の第１ボクセルの信号値が前記第２画像の第２ボクセルの信号値と同一であると仮定すると、前記ウェイトは、前記第１画像における前記第１ボクセル及び第２ボクセルのボクセル値と、前記第２画像における第１ボクセル及び第２ボクセルのボクセル値とを取得する尤度を含む、ジョイントバイラテラルフィルタリングするステップと、
前記フィルタリングされた画像を提供するステップと、
を有する画像処理方法。
【請求項１３】
コンピュータ上で実行されると、前記コンピュータに請求項１２記載の方法のステップを実行させるためのプログラムコード手段を有するコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理装置及び画像をフィルタリングするための対応する方法に関する。
【背景技術】
【０００２】
患者又はオブジェクトの単一の画像のみでなく、正確に同じジオメトリにおけるそれらの２以上を提供する新規なイメージングモダリティが出現している。これの２つの具体例は、複素数の屈折率の実部及び虚部が測定される微分位相コントラストイメージング（ＤＰＣＩ；ＵＳ２０１２／００９９７０２Ａ１などに記載される）と、光電効果、コンプトン散乱及びおそらく造影剤によるＸ線の減衰が復元されるスペクトルＣＴ（ＵＳ２００８／０２５３５０４Ａなどに記載される）とである。
【０００３】
ＷＯ２０１１／１４５０４０Ａ１では、ジョイントバイラテラルフィルタリングが、画像ベースノイズ除去のための周知のバイラテラルフィルタの一般化として提案された。特に、オブジェクトの第１画像及び同一のオブジェクトの第２画像を提供する画像提供ユニットと、第１画像の画像値の間の第１類似度と、第２画像の画像値の間の第２類似度とに依存して第１画像をフィルタリングするフィルタリングユニットとを有する画像処理装置が説明されている。これは、オブジェクトが人間又は動物である場合、ノイズのために第１画像及び第２画像の一方の画像値が乱れたとしても、画像値がオブジェクトの同一部分、例えば、オブジェクトの組織又は骨物質に属する確率に依存して第１画像をフィルタリングすることを可能にし、これにより、例えば、ＣＴ画像のエッジ保存画像ベースノイズ除去などを向上させる。
【０００４】
しかしながら、既知の方法及び装置は、各種画像におけるノイズが相関しうることを認めていない。特に、ＤＰＣＩについては、位相コントラスト画像内のトランスアキシャル画像において強い空間的な相関があり（低周波数においてピークを有する周知のノイズパワースペクトルに等しい）、スペクトルＣＴについては、フォト画像とコンプトン画像とにおける同じ空間位置においてノイズの強力な反相関がある。
【０００５】
同じオブジェクトの２つの画像が利用可能であるが、異なる空間解像度においてノイズにおける異なる空間的相関によって反映される他の状況がまたあってもよい。具体例は、ある状況では、解剖学的（Ｔ２加重化など）及びパラメータ的画像（Ｔ１）が取得される磁気共鳴イメージング（ＭＲＩ）であってもよい。
【発明の概要】
【０００６】
画像をフィルタリングするための改良された画像処理装置及び方法と、特にジョイント画像ベースノイズ除去のための方法とを提供することが、本発明の課題である。
【０００７】
本発明の第１の態様では、画像をフィルタリングする画像処理装置であって、
同じオブジェクトの第１及び第２画像を取得する画像入力であって、前記第１及び第２画像は複数のボクセルを有し、ノイズ共分散によって相関し、各ボクセルは信号値とノイズ値とを含むボクセル値を有する、画像入力と、
前記第１画像をフィルタリングするジョイントバイラテラルフィルタであって、前記第１画像の第１ボクセルは相対的なボクセルに固有のウェイトを含むフィルタ関数によりフィルタリングされ、前記第１画像の第１ボクセルの信号値が前記第１画像の第２ボクセルの信号値と同一であり、前記第２画像の第１ボクセルの信号値が前記第２画像の第２ボクセルの信号値と同一であると仮定すると、前記ウェイトは、前記第１画像における前記第１ボクセル及び第２ボクセルのボクセル値と、前記第２画像における第１ボクセル及び第２ボクセルのボクセル値とを取得する尤度を含む、ジョイントバイラテラルフィルタと、
前記フィルタリングされた画像を提供する画像出力と、
を有する画像処理装置が提供される。
【０００８】
本発明の更なる態様では、対応する画像処理方法が提供される。
【０００９】
本発明のまた更なる態様では、コンピュータ上で実行されると、コンピュータに処理方法のステップを実行させるためのプログラムコード手段を有するコンピュータプログラムと共に、プロセッサにより実行されると、ここに開示される方法を実行させるコンピュータプログラムプロダクトを格納する非一時的なコンピュータ可読記憶媒体が提供される。
【００１０】
本発明の好適な実施例は、従属形式の請求項に規定される。請求される方法、コンピュータプログラム及び媒体は、従属形式の請求項に規定され、請求される装置と類似及び／又は同一の好適な実施例を有すると理解される。
【００１１】
本発明によると、改良されたジョイントバイラテラルフィルタリングのためのアプローチが提案される。提案されたジョイントバイラテラルフィルタアプローチは、基本的には２つの項、距離重み付け及びレンジ重み付けを含む。この第１のものは、純粋に幾何学的に動機付けされ、変更されない。レンジフィルタは、本発明により変更されることが提案される画像位置に基づく重み付けを導入する。特に、画像間及び異なる画像位置間のノイズ相関が、改良されたジョイント画像ベースノイズ除去を提供するため、本発明により考慮される。
【００１２】
このコンテクストにおいて、“相関（ｉｎｔｅｒｒｅｌａｔｅｄ）”は、当該相関が同じ画像のボクセル間の相関に限定可能であるか、又は異なる画像における同じボクセル位置の間の相関に限定されてもよいように理解されるであろう。
【００１３】
好適な実施例では、前記フィルタは、前記ノイズ共分散を記述するノイズ共分散行列を含む尤度を利用するよう構成される。共分散を利用することによって、ノイズがオン方向に高く相関するが（ＤＰＣＩのためのイン・プレーン（ｉｎ−ｐｌａｎｅ）など）、ｚ方向には大きくは相関しないケースなどを解決することが可能になる。このような状況では、レンジフィルタの定数はイン・プレーン又はアウト・オブ・プレーン（ｏｕｔ−ｏｆ−ｐｌａｎｅ）ノイズ特性のみに調整可能であるため、従来の（ジョイント）バイラテラルフィルタリングは良好に実行されない。
【００１４】
好ましくは、前記フィルタは、
【００１５】
【数１】
2016506267.tif 000003
又は
【００１６】
【数２】
2016506267.tif 000004
の形式の共分散を利用するよう構成される。これらの共分散行列は効果的であると示した。第２共分散行列は扱うのがより容易である。
【００１７】
以下の特別なケースが特に興味がある。微分位相コントラストイメージングについて、再構成された屈折率の虚部（第１画像）と実部（第２画像）との間にはほとんどノイズ相関がなく、
【００１８】
【数３】
2016506267.tif 000005
を導き（用語はＣの一般的なバージョンに関し、すなわち、上述された２つのバージョンの上方の共分散行列）、一方、屈折率の実部においてトランスアキシャル画像内に強力なノイズ相関があり、すなわち、ｃ３４＞０である。共分散行列は、もちろん、ボクセルポジションに依存し、例えば、第２画像の第２ボクセルが異なるトランスアキシャル画像にある場合、ノイズは無相関であり、すなわち、
【００１９】
【数４】
2016506267.tif 000006
である。
【００２０】
スペクトルイメージングについて、各画像内にノイズの空間的な相関はほとんどなく（すなわち、コンプトン及び光電画像）、
【００２１】
【数５】
2016506267.tif 000007
を導き、また、異なる画像における異なる位置の間には相関がなく、すなわち、
【００２２】
【数６】
2016506267.tif 000008
であるが、同じ位置において２つの画像における画像値の間には強力な負の相関があり、ｃ１３＜０，ｃ２４＜０である。
【００２３】
さらに、実施例では、前記フィルタは、
【００２４】
【数７】
2016506267.tif 000009
の形式の尤度を利用するよう構成され、Ｓは前記共分散行列の逆行列であり、ａ１及びａ２は前記第１画像の第１及び第２ボクセル値であり、ｂ１及びｂ２は前記第２画像の第１及び第２ボクセル値であり、Δａ＝ａ１−ａ２，Δｂ＝ｂ１−ｂ２である。
【００２５】
さらに、第１画像のフィルタリングされた第１ボクセル値は、
【００２６】
【数８】
2016506267.tif 000010
として取得され、ｗｄは標準的な（ジョイント）バイラテラルフィルタリングにおける（通常の）距離ウェイティングである。距離｜ｉ−ｊ｜の増加により単調減少する何れかの関数、例えば、
【００２７】
【数９】
2016506267.tif 000011
が利用されてもよい。
【００２８】
実際的な状況では、第１及び第２画像の少なくとも１つは非白色ノイズを有し、第１画像と第２画像との間にはノイズ相関が存在し、及び／又は第１画像及び／又は第２画像の第１及び第２ボクセル値におけるノイズ共分散は同じである。
【００２９】
好ましくは、同じフィルタが第２画像に適用される。さらに、提案されたアイデアは、３つ（又はより多くの）画像による状況に直接的な方法で拡張可能である。
【００３０】
他の実施例によると、前記画像入力は、前記オブジェクトの異なる性質を示すように、前記第１画像と前記第２画像とを提供するよう構成される。好ましくは、前記画像入力は、スペクトルＣＴ画像を提供するよう構成され、前記第１画像は前記オブジェクトの第１の物理的性質を示し、前記第２画像は前記オブジェクトの第２の物理的性質を示し、前記第１の物理的性質と前記第２の物理的性質とは互いに異なる。
【００３１】
例えば、前記第１画像と前記第２画像との一方はコンプトン画像であり、前記第１画像と前記第２画像との他方は光電画像である。他の実施例では、前記第１画像と前記第２画像との一方はコンプトン画像又は光電画像であり、前記第１画像と前記第２画像との他方はコンプトン画像と光電画像との組み合わせである。
【００３２】
またさらに、前記画像入力は、微分位相コントラストＣＴ画像を提供するよう構成され、前記第１画像と前記第２画像との一方は前記オブジェクトの吸収を示し、前記第１画像と前記第２画像との他方は前記オブジェクトにより誘導される前記オブジェクトをトラバースする放射線の位相シフトを示す。
【図面の簡単な説明】
【００３３】
本発明の上記及び他の態様は、後述される実施例を参照して明らかになるであろう。
【図１】図１は、本発明による画像処理装置の実施例を概略的及び例示的に示す。
【図２】図２は、画像入力ユニットの実施例を概略的及び例示的に示す。
【図３】図３は、多色放射ソースの発光スペクトルを概略的及び例示的に示す。
【図４】図４は、表示されたエネルギー範囲内のｋエッジによる物質の光電効果、コンプトン効果及びトータルの減衰の断面図のエネルギー依存を概略的及び例示的に示す。
【図５】図５は、微分位相コントラストＣＴシステムの実施例の要素を概略的及び例示的に示す。
【図６】図６は、画像を処理する画像処理方法の実施例を示すフローチャートを例示的に示す。
【発明を実施するための形態】
【００３４】
図１は、本発明による画像処理装置１の実施例を概略的に示す。それは、同一のオブジェクトの第１及び第２画像を取得する画像ユニット２を有し、第１及び第２画像は複数のボクセルを有し、ノイズ共分散により相関し、各ボクセルは信号値及びノイズ値を有するボクセル値を有する。装置１は更に、第１画像をフィルタリングするバイラテラルフィルタ３を有し、第１画像の第１ボクセルは相対的なボクセルに固有のウェイトを含むフィルタ関数によりフィルタリングされる。第１画像の第１ボクセルの信号値が第２画像の第１ボクセルの信号値と同じであり、第１画像の第２ボクセルの信号値が第２画像の第２ボクセルの信号値と同じであると仮定すると、当該ウェイトは、当該第１ボクセルと第２画像における第２ボクセルとのボクセル値を取得し、第２画像における第１ボクセルと第２ボクセルとを取得する尤度（ｌｉｋｅｌｉｈｏｏｄ）を含む。装置１は更に、フィルタリングされた画像を表示するためのディスプレイなどであってもよい、フィルタリングされた画像を提供する画像出力４を有する。
【００３５】
第１画像と第２画像とは、例えば、オブジェクトのスライスのみなど、オブジェクト全体又はオブジェクトの一部のみとすることが可能な同一のオブジェクトを示す２次元、３次元又は４次元画像とすることができる。オブジェクトが人間の心臓などの挙動するオブジェクトである場合、第１画像と第２画像とは、好ましくは、同じ挙動状態におけるオブジェクトを示す。本実施例の第１画像及び第２画像は、プロジェクションデータに基づき、特に、ＣＴプロジェクションデータに基づき再構成され、ここで、第１画像及び第２画像は同一のプロジェクションデータから再構成される。
【００３６】
画像入力ユニット２は、それらがオブジェクトの異なる性質、特に、異なる物理的性質を示すように、第１画像と第２画像とを提供するよう構成される。実施例では、画像入力ユニット２は、スペクトルＣＴ（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）画像を提供するよう構成され、第１画像はオブジェクトの第１の物理的性質を示し、第２画像はオブジェクトの第２の物理的性質を示し、第１の物理的性質と第２の物理的性質とは互いに異なる。例えば、第１画像と第２画像との一方は光電画像とすることができ、あるいは、第１画像と第２画像との一方はコンプトン画像又は光電画像とすることができ、第１画像と第２画像との他方はコンプトン画像と光電画像との組み合わせとすることができる。
【００３７】
画像入力ユニット２は格納ユニットとすることが可能であり、そこでは、第１画像と第２画像とが格納されるか、あるいは、画像入力ユニット２はまた、データ接続を介し第１画像及び／又は第２画像を受信することを可能にする画像受信ユニットとすることができる。画像入力ユニット２はまた、第１画像と第２画像とを生成するスペクトルＣＴシステムなどのイメージングモダリティとすることができる。このようなスペクトルＣＴシステムは、以下において図２を参照して概略的及び例示的に説明される。
【００３８】
スペクトルＣＴシステムは、ｚ方向にパラレルに延びる回転軸Ｒの周りに回転可能なガントリ１４を含む。Ｘ線チューブなどの放射ソース１５は、ガントリ１４に搭載される。本実施例では、放射ソース１５は多色放射線を発する。放射ソース１５には、放射ソース１５により発せられた放射から円錐放射ビーム１７を形成するコリメータ装置１６が備えられる。他の実施例では、コリメータ装置１６は、扇形形状などを有する他の形状を有する放射ビームを形成するよう構成可能である。
【００３９】
放射線は、円筒形の検査ゾーン５における関心領域の患者又は技術的オブジェクトなどのオブジェクト（図２に図示せず）をトラバースする。関心領域をトラバースした後、放射ビーム１７は、本実施例では、２次元検出面を有する検出装置６に入射し、ここで、検出装置６はガントリ１４に搭載される。他の実施例では、検出装置６は１次元検出面を有することが可能である。
【００４０】
検出装置６は、好ましくは、エネルギー分解検出要素を有する。エネルギー分解検出要素は、好ましくは、例えば、入射光子をカウントする原理などにより動作し、異なるエネルギーウィンドウにおける光子数を示す信号を出力する光子計数検出要素である。このようなエネルギー分解検出装置は、例えば、Ｌｌｏｐａｒｔ  Ｘ．，ｅｔ  ａｌ．，“Ｆｉｒｓｔ  ｔｅｓｔ  ｍｅａｓｕｒｅｍｅｎｔｓ  ｏｆ  ａ  ６４ｋ  ｐｉｘｅｌ  ｒｅａｄｏｕｔ  ｃｈｉｐ  ｗｏｒｋｉｎｇ  ｉｎ  ａ  ｓｉｎｇｌｅ  ｐｈｏｔｏｎ  ｃｏｕｎｔｉｎｇ  ｍｏｄｅ”，Ｎｕｃｌ．Ｉｎｓｔ．ａｎｄ  Ｍｅｔｈ．Ａ，５０９（１−３）：１５７−１６３，２００３及びＬｌｏｐａｒｔ，Ｘ．，ｅｔ  ａｌ．，“Ｍｅｄｉｐｉｘ２：Ａ  ６４−ｋ  ｐｉｘｅｌ  ｒｅａｄｏｕｔ  ｃｈｉｐ  ｗｉｔｈ  ５５  μｍ  ｓｑｕａｒｅ  ｅｌｅｍｅｎｔｓ  ｗｏｒｋｉｎｇ  ｉｎ  ａ  ｓｉｎｇｌｅ  ｃｏｕｎｔｉｎｇ  ｍｏｄｅ”，ＩＥＥＥ  Ｔｒａｎｓ．Ｎｕｃｌ．Ｓｃｉ．４９（５）：２２７９−２２８３，２００２に記載されている。好ましくは、各検出要素が少なくとも２つの異なるエネルギーウィンドウのための少なくとも２つのエネルギー分解された検出信号を提供するように、エネルギー分解検出要素が構成される。しかしながら、イメージングシステムの感度及びノイズロウバストネスを向上させるため、さらに高いエネルギー分解能を有することが好ましい。
【００４１】
ガントリ１４は、好ましくは、モータ７によって一定であるが、調整可能な角速度により駆動される。更なるモータ８は、検査ゾーン５の患者テーブルに配置された患者などのオブジェクトを、回転軸方向又はｚ軸にパラレルにずらすため設けられる。放射ソース１５と検査ゾーン５、特に関心領域とが螺旋の軌跡に沿って互いに対して移動するように、これらのモータ７，８は制御ユニット９などにより制御される。また、オブジェクト又は検査ゾーン５、特に関心領域は移動せず、放射ソース１５が回転し、すなわち、放射ソース１５が関心領域に対して円の軌跡に沿って移動することが可能である。検出装置６により取得されるデータは、本実施例では、提供されたエネルギーに依存した検出データから関心領域の第１及び第２画像を少なくとも再構成する再構成ユニット１０に提供されるエネルギーに依存した検出データである。また、再構成ユニット１０は、好ましくは制御ユニット９により制御される。再構成された画像は、少なくとも第１画像をフィルタリングするフィルタリングユニット３に提供される。
【００４２】
再構成ユニット１０は、検出データから少なくとも２つの減衰コンポーネントを決定する計算ユニット１２を有する。検出データの減衰コンポーネントは、１つのみ又は複数の物理的効果によって引き起こされ、例えば、ある減衰コンポーネントはコンプトン効果により引き起こされるコンポーネントでありうるものであり、他の減衰コンポーネントは光電効果により引き起こされるコンポーネントでありうる。更なる減衰コンポーネントは関心領域に存在するｋエッジにより引き起こされるコンポーネントでありうる。あるいは、又はさらに、減衰コンポーネントは、例えば、関心領域内のある物質の吸収などによって引き起こされるコンポーネントでありうる。例えば、減衰コンポーネントはある物質の吸収により引き起こされるコンポーネントでありうるものであり、他の減衰コンポーネントは他の物質の吸収により引き起こされるコンポーネントでありうる。
【００４３】
再構成ユニット１０は更に、対応する減衰コンポーネント画像が再構成されるように、検出データの計算された減衰コンポーネントをバックプロジェクションするバックプロジェクションユニット１３を有する。例えば、コンプトン画像は、第１画像と第２画像との一方として再構成可能であり、光電画像は、第１画像と第２画像との他方として再構成可能である。
【００４４】
計算ユニット１２への入力は、複数のエネルギービン、本実施例では、最小で３つのエネルギービンのエネルギー分解検出データｄｉである。これらの検出データｄｉは、第ｉエネルギービンｅｉのスペクトル感度Ｄｉ（Ｅ）を示す。さらに、多色放射ソース１５の発光スペクトルＴ（Ｅ）は、一般に既知であるか、又は測定可能である。多色放射ソースのこのような発光スペクトルＴ（Ｅ）の具体例は、図３において概略的及び例示的に示される。計算ユニット１２では、検出データｄｉの生成は、スペクトルによる光電効果Ｐ（Ｅ）、スペクトルによるコンプトン効果Ｂ（Ｅ）及びスペクトルによるｋエッジＫ（Ｅ）の線形結合としてモデル化される。
【００４５】
スペクトルＰ（Ｅ）、Ｂ（Ｅ）及びＫ（Ｅ）は、図４において例示的及び概略的に示される。
【００４６】
検出データの生成は、以下の連立方程式によりモデル化できる。
【００４７】
【数１０】
2016506267.tif 000012
ここで、ρｐｈｏｔｏ，ρＣｏｍｐｔｏｎ及びρｋ−ｅｄｇｅはそれぞれ、光電コンポーネント、コンプトンコンポーネント及びｋエッジコンポーネントの密度長積（ｄｅｎｓｉｔｙ  ｌｅｎｇｔｈ  ｐｒｏｄｕｃｔ）である。
【００４８】
少なくとも３つの検出信号ｄ１，ｄ２，ｄ３が少なくとも３つのエネルギービンｅ１，ｅ２，ｅ３について利用可能であるため、３つの密度長積である３つの未知を有する少なくとも３つの方程式の連立方程式が構成され、従って、計算ユニット１２において既知の数値方法により解くことができる。３より多くのエネルギービンが利用可能である場合、測定結果のノイズ統計を考慮する最大尤度アプローチを利用することが好ましい。一般に、３つのエネルギービンで十分である。しかしながら、感度及びノイズロウバストネスを向上させるため、より多くのエネルギービンについてより多くの検出データを有することが好ましい。各減衰コンポーネントについて、すなわち、各密度長積について、オブジェクトの減衰コンポーネント画像は、各密度長積をバックプロジェクションすることによって再構成可能である。あるいは、又はさらに、異なる密度長積の画像、又はバックプロジェクション前の異なる密度長積の画像が、合成された減衰コンポーネントにより生成されるオブジェクトの画像である合成画像を生成するため合成可能である。
【００４９】
他の実施例では、画像入力ユニット２は、微分位相コントラストＣＴ画像を提供するよう構成可能であり、第１画像及び第２画像の一方はオブジェクトの吸収を示し、第１画像及び第２画像の他方はオブジェクトにより誘導される、オブジェクトをトラバースする放射位相シフトを示す。特に、第１画像は、好ましくは、オブジェクトのトータルの減衰係数を示し、第２画像は、好ましくは、オブジェクトの屈折率の実部を示す。また、これらの第１及び第２画像は、画像入力ユニット２であってもよい格納ユニットに格納可能であるか、又は画像入力ユニット２であってもよい画像受信ユニットを介し受信可能である。しかしながら、画像入力ユニット２はまた、第１画像及び第２画像を生成し、図５を参照して以下で例示的及び概略的に説明される微分位相コントラストＣＴシステムとすることができる。
【００５０】
図５は、微分位相コントラストＣＴシステムの一例となる実施例の３次元表現を示す。やや大きなＸ線ソース３２はソースグレイティング３４に隣接して配置される。Ｘ線ソース３２は発せられる放射の波長に関するそれのサイズのため非干渉であると考えられるため、ソースグレイティング３４は、複数の単独のコヒーラントなＸ線ソースを提供するのに利用される。
【００５１】
Ｘ線放射３５は、Ｘ線ソース３２から光軸３７の方向に発せられ、おそらくＸ線の扇形ビーム又は円錐型ビームを構成する。図５において、Ｘ線ビームの各形状は示されていない。Ｘ線放射３５はオブジェクト３６に到来し、オブジェクト３６を貫通し、その後にビームスプリッタグレイティング３８に到来する。ビームスプリッタグレイティング３８のトレンチ又はギャップは、ビームスプリッタグレイティングの固体エリア、すなわち、ブロッキング領域に関して電磁放射を通過させる位相を変更する。従って、φによる、特にπによる位相シフトが実行される。
【００５２】
解析グレイティング４０は、ビームスプリッタグレイティング３８とＸ線検出装置４２との間に配置される。Ｘ線検出装置の方向にビームスプリッタグレイティング３８から発せられる複数の波は、解析グレイティング４０に到来し、その後、Ｘ線検出装置４２の表面に強度変調パターンを生成する。
【００５３】
解析グレイティング４０に対してビームスプリッタグレイティング３８をシフトすることによって、互いに対して、特にビームスプリッタグレイティング３８又は解析グレイティング４０のグレイティング期間の一部によってグレイティングをずらすことは、画像検出装置４２により取得可能である。従って、送信画像の取得中、いわゆる、位相ステッピングが実行され、ここでは、グレイティング３８，４０の相対位置が、１つの期間において４〜１０のステップなどにおいて変更され、各相対位置について送信画像が測定される。この送信画像系列から、従来の吸収画像が、各検出ピクセルについて測定された強度を単純平均することによって取得される。従来のＸ線送信イメージングと同様に、この測定はＢｅｅｒの法則によってＸ線パスに沿ってトータルの減衰係数に対する線積分に関連する。屈折率の実部に関する情報は、位相ステッピングサイクルにおける測定された強度のダイナミクスの解析により取得される。当該強度は、好ましくは、解析グレイティング４０のグレイティング期間により定期的に変動する。典型的には、正弦波変動が想定できる。サイナスの位相はビームスプリッタグレイティング３８にヒットする位相フロントの勾配に関連する。物理学によって、位相フロントの当該勾配は、オブジェクトの屈折率の実部のプロジェクションの第１導関数に線形に相関する。従って、オブジェクトの屈折率の実部のプロジェクションは、サイナスの位相に基づき決定可能である。オブジェクトの屈折率の実部の決定されたプロジェクションをバックプロジェクションすることによって、オブジェクトの屈折率の実部を示す画像が再構成可能である。微分位相コントラストＣＴシステムは、例えば、オブジェクトのトータルの減衰係数を示す第１画像とオブジェクトの屈折率の実部を示す第２画像とを再構成可能である。既知の気分いそうコントラストＣＴ技術のより詳細な説明は、例えば、参照することによりここに援用される、Ｔｉｍｍ  Ｗｅｉｔｋａｍｐ  ｅｔ  ａｌ．による論文“Ｘ−ｒａｙ  ｐｈａｓｅ  ｉｍａｇｉｎｇ  ｗｉｔｈ  ａ  ｇｒａｔｉｎｇ  ｉｎｔｅｒｆｅｒｏｍｅｔｅｒ”，ＯＰＴＩＣＳ  ＥＸＰＲＥＳＳ，ｖｏｌ．１３，ｎｏ．１６，ｐａｇｅｓ  ６２９６−６３０４（Ａｕｇｕｓｔ  ２００５）などに開示される。
【００５４】
ＷＯ２０１１／１４５０４０Ａ１には、ジョイントバイラテラルフィルタリングが、画像ベースノイズ除去のための周知なバイラテラルフィルタの一般化として提案された。このアプローチでは、フィルタリングユニットは、好ましくは、第１画像が以下のフィルタリング方程式を利用することによってノイズ除去されるように、フィルタリングを実行するよう構成される。
【００５５】
【数１１】
2016506267.tif 000013
ここで、
【００５６】
【数１２】
2016506267.tif 000014
は第１画像のポジション
【００５７】
【数１３】
2016506267.tif 000015
における第１画像要素の平均画像値を示し、Δｉはポジション
【００５８】
【数１４】
2016506267.tif 000016
における第１画像要素を有する第１領域を示し、
【００５９】
【数１５】
2016506267.tif 000017
は第１領域Δｉ内の第１領域画像要素の位置を示し、
【００６０】
【数１６】
2016506267.tif 000018
は、ａ）第１領域画像要素
【００６１】
【数１７】
2016506267.tif 000019
の画像値と、任意的にはポジション
【００６２】
【数１８】
2016506267.tif 000020
における各第１領域要素を包囲する画像要素の画像値と、ｂ）ポジション
【００６３】
【数１９】
2016506267.tif 000021
における第１画像要素の画像値と、任意的には第１画像の第１画像要素を包囲する画像要素の画像値との間の第１類似度を示す。
【００６４】
関数
【００６５】
【数２０】
2016506267.tif 000022
は、第２画像の第１領域に対応する第２領域内のポジション
【００６６】
【数２１】
2016506267.tif 000023
における第２領域画像要素の画像値と、任意的にはポジション
【００６７】
【数２２】
2016506267.tif 000024
における第２領域要素を包囲する画像要素の画像値と、ｂ）第２画像の第１画像における同一ポジションにおける第１画像要素に対応するポジション
【００６８】
【数２３】
2016506267.tif 000025
における第２画像要素の画像値と、任意的には第２画像要素を包囲する画像要素の画像値との間の第２類似度を示す。関数
【００６９】
【数２４】
2016506267.tif 000026
は、ポジション
【００７０】
【数２５】
2016506267.tif 000027
における画像要素の間の距離に依存する距離関数であり、
【００７１】
【数２６】
2016506267.tif 000028
は第１画像におけるポジション
【００７２】
【数２７】
2016506267.tif 000029
における非フィルタリング画像値を示す。
【００７３】
式（２）に示される分母により除算される第１類似度Ｌ１、第２類似度Ｌ２及び距離関数ｗｄの合成は、式（２）により加重平均を実行するためのウェイトとしてみなすことができる。従って、式（２）に従って、フィルタリングユニット３は、好ましくは、第１画像の画像要素の画像値を加重平均するよう構成され、ここで、第１画像の第１画像要素の平均画像値を決定するため、第１画像要素を含む第１領域内の第１領域画像要素の画像値が重み付け及び平均化され、第１領域画像要素の画像値は、ａ）第１領域画像要素の画像値と、任意的には各第１画像要素を包囲する画像要素の画像値と、ｂ）第１画像要素の画像値と、第１画像の第１画像要素を包囲する画像要素の画像値との間の類似度である第１類似度に依存して、また、ａ）第２画像の第１領域に対応する第２領域内の第２領域画像要素の画像値と、任意的には各第２領域画像要素を包囲する画像要素の画像値と、ｂ）第２画像の第１画像要素に対応する第２画像要素の画像値と、任意的には第２画像要素を包囲する画像要素の画像値との間の類似度である第２類似度に依存して、ウェイトにより重み付けされる。
【００７４】
ベクトル
【００７５】
【数２８】
2016506267.tif 000030
は、
【００７６】
【数２９】
2016506267.tif 000031
を中心とする画像パッチとしてみなすことができる。例えば、各第１画像要素を包囲する画像要素の任意的な画像値と、各第２画像要素を包囲する画像要素の任意的な画像値とが考慮される場合、２次元のケースでは、画像パッチはポジション
【００７７】
【数３０】
2016506267.tif 000032
におけるピクセルを中心とする９×９のサブ画像であってもよい。当該パッチはまた他のサイズを有することも可能である。例えば、パッチは５×５又は７×７のサイズを有することが可能である。さらに、３次元パッチがまた利用可能である。任意的な画像値が考慮されない場合、すなわち、ポジション
【００７８】
【数３１】
2016506267.tif 000033
における単一の画像値の間の差分しか考慮されない場合、各ベクトルはそれぞれ、ポジション
【００７９】
【数３２】
2016506267.tif 000034
における単一の各画像値しか有さない。このとき、パッチサイズは１×１である。
【００８０】
フィルタリングが２より多くの画像に基づき実行される場合、フィルタリングユニット３は、好ましくは、以下のフィルタリング方程式に従って第１画像をフィルタリングするよう構成される。
【００８１】
【数３３】
2016506267.tif 000035
ここで、Ｎは第１画像をフィルタリングするため考慮される画像数を示し、インデックスｋは、各要素がｋ番目の画像に属することを示す。
【００８２】
第ｋ類似度は、以下の方程式により定義できる。
【００８３】
【数３４】
2016506267.tif 000036
ここで、
【００８４】
【数３５】
2016506267.tif 000037
はそれぞれ、
【００８５】
【数３６】
2016506267.tif 000038
におけるローカルノイズ推定値を示す。
【００８６】
ノイズ推定値は、実際に同質な領域であると考えられる領域内の標準偏差の決定によって、手動により取得可能である。また、既知の初期的な生データにおけるノイズから始めて良好な近似までの古典的なノイズ伝搬によって空間的に分解されたノイズ推定値を取得することが可能である。既知の古典的なノイズ伝搬は、例えば、参照することによりここに援用される、Ａｄａｍ  Ｗｕｎｄｅｒｌｉｃｈ  ｅｔ  ａｌ．，による論文“Ｉｍａｇｅ  ｃｏｖａｒｉａｎｃｅ  ａｎｄ  ｌｅｓｉｏｎ  ｄｅｔｅｃｔａｂｉｌｉｔｙ  ｉｎ  ｄｉｒｅｃｔ  ｆａｎ−ｂｅａｍ  ｘ−ｒａｙ  ｃｏｍｐｕｔｅｄ  ｔｏｍｏｇｒａｐｈｙ”，Ｉｎｓｔｉｔｕｔｅ  ｏｆ  Ｐｈｙｓｉｃｓ  ａｎｄ  Ｅｎｇｉｎｅｅｒｉｎｇ  ｉｎ  Ｍｅｄｉｃｉｎｅ，ｐａｇｅｓ  ２４７１  ｔｏ  ２４９３（２００８）などに開示される。同一のコンセプトは、空間ＣＴだけでなく、微分位相コントラストＣＴに適用可能である。
【００８７】
ローカルノイズ推定値の代わりに、グローバルノイズ推定値がまた、特にノイズが画像全体で大きくは変動しない場合、利用可能である。
【００８８】
上述の方程式は、ａ）第１領域画像要素と第１画像要素との間の第１類似度と、ｂ）第２領域画像要素と第２画像要素との間の第２類似度との一方が減少する場合、かつ、ａ）第１領域画像要素と第１画像要素との間の第１類似度と、ｂ）第２領域画像要素と第２画像要素との間の第２類似度との他方が変化しないか、若しくは減少する場合、第１領域画像要素の画像値のウェイトは減少することを保証する。
【００８９】
これらの類似度は、特に、各画像がノイズ除去される場合、２つの画像値又は２つの画像パッチがどの程度同じであるかの推定値を提供する。従って、式（４）において定義される類似度はまた、画像値が類似する確率の一緒としてみなすことができる。ノイズよりはるかに小さな相違について、尤度又は確率は１に近く、類似度はノイズレベルより大きな相違についてゼロに低下する。特に、１×１より大きなサイズを有する画像パッチが考慮される場合、類似度のガウスプロトタイプのため、ポジション
【００９０】
【数３７】
2016506267.tif 000039
における画像値の結果としてのウェイトは、
【００９１】
【数３８】
2016506267.tif 000040
の周りの画像パッチがＮ個の画像の何れかにおいて統計的に異なる場合、小さくなるであろう。従って、第１画像における統計的に有意なエッジにわたるスムージングは実行されないことが保証できる。従って、フィルタリングはエッジ保存的である。
【００９２】
空間加重関数としてもみなすことができる距離関数は、以下の方程式により定義できる。
【００９３】
【数３９】
2016506267.tif 000041
ここで、パラメータσｗｄは、平均化中に考慮される近傍の有効サイズを制御する。
【００９４】
以下において、ここに提案されるようなジョイントバイラテラルフィルタリングのための更に向上されたアプローチが概略される。バイラテラルフィルタリングアプローチは、基本的には２つの項、距離重み付け（上記の式（５）及び（６）における“ｗｄ”として参照される）と、レンジ重み付け（上記の式（４）における“Ｌ”として参照される）とを含む。距離重み付けは、純粋に幾何学的に動機付けされ、上述した既知のアプローチと比較して変更されない。レンジフィルタは、上述された既知のアプローチと比較して本発明により変更されることが提案される画像値に基づく重み付けを導入する。
【００９５】
具体例は、２つの画像について説明される。バイラテラルフィルタリングの基本タスクは、２つの異なるピクセルが同じ組織クラスに属する確率、すなわち、基礎となる真のピクセル値が同じである尤度を推定することである。この推定は、２つの画像における同じ２つの異なる位置におけるサンプルである４つ（又はそれ以上）の画像サンプルの評価に基づく。これらのサンプルは、ａ１，ａ２，ｂ１，ｂ２として示され（ａｉ及びｂｉは上述した
【００９６】
【数４０】
2016506267.tif 000042
に対応し、これらは、よりコンパクトな表記を実現するため以下において利用される）、ベクトル（ａ１，ａ２，ｂ１，ｂ２）を構成する。これらのサンプルに対する２次統計量は、（新たに導入された）共分散行列Ｃ
【００９７】
【数４１】
2016506267.tif 000043
によって、これらの測定結果におけるノイズを記述し、ここで、対角要素はサンプルのノイズ共分散を記述し、非対角要素は画像（ｃ１２，ｃ３４）内又は画像（ｃ１３，ｃ１４，ｃ２３，ｃ２４）間のノイズの相関を記述する。行列Ｃは対称的であり、正定値である。この共分散行列によって、測定されたベクトル（ａ１，ａ２，ｂ１，ｂ２）においてノイズ実現（ｎ１，ｎ２，ｎ３，ｎ４）を取得する尤度を記述できる。
【００９８】
【数４２】
2016506267.tif 000044
この測定ではノイズと信号コンテンツとを区別することはできない。従って、ａ１＝ａ２及びｂ１＝ｂ２の仮定の下、どの程度の確率で実際の測定結果（ａ１，ａ２，ｂ１，ｂ２）を観察することになるかが調べられる。この仮定の下では、測定における信号コンテンツは差分Δａ＝ａ１−ａ２，Δｂ＝ｂ１−ｂ２を計算し、これらの差分を生じさせる（ｎ１，ｎ２，ｎ３，ｎ４）の全ての組み合わせに対して積分することによって除外できるため、これは所与の統計量に対して容易に回答できる。表記のコンパクト化のため、Ｓ＝Ｃ−１と設定される。所望の確率は、
【００９９】
【数４３】
2016506267.tif 000045
に従って計算できる。この周辺確率は、ジョイントバイラテラルフィルタ方程式におけるジョイント尤度項として挿入される必要があり、
【０１００】
【数４４】
2016506267.tif 000046
を生じさせる。
【０１０１】
最も一般的な共分散行列について、積分の明示的な計算はかなり複雑になる。しかしながら、共分散行列Ｃが数値的に既知になると、実際の測定結果の確率と逆行列Ｓとを計算することは容易な作業である。
【０１０２】
２つのシンプルなケースが考えられる。４つのサンプルのノイズが無相関である場合、共分散行列とそれの逆行列との双方は対角となり、所望される尤度はガウス関数の積となる。他の興味深いケースは、
【０１０３】
【数４５】
2016506267.tif 000047
の形式の共分散行列についてである。幾何学的には、この共分散行列は、２つの画像位置における画像共分散が同じであり（すなわち、ｃ１１＝ｃ２２，ｃ３３＝ｃ４４）、画像間の共分散が２つの画像位置（ｃ１３＝ｃ２４，ｃ１４＝ｃ２３）に関して対称的であるケースをモデル化する。典型的な状況では、これらの関係は、再構成されたＣＴ画像における分散は空間においてゆっくりとしか変化しないため、ほぼ充足される。
【０１０４】
この増加した対称性はまた、逆行列Ｓにおいて反映され、それは、
【０１０５】
【数４６】
2016506267.tif 000048
を読む。
【０１０６】
明らかに、共分散行列の対称性がまたそれの逆行列に存在し、従って、
【０１０７】
【数４７】
2016506267.tif 000049
として記述できる。この仮定の下、ベクトル（ａ１，ａ２，ｂ１，ｂ２）を測定するための所望の尤度は、
【０１０８】
【数４８】
2016506267.tif 000050
である。
【０１０９】
フィルタの式（１６）に挿入されるような式（１４）に与えられる結合尤度は、従来技術において行われるような２つの画像に対する導出された独立した類似度の積の形態を有さないが、追加的な項、すなわち、差分の積ΔａΔｂがあることに留意されたい。
【０１１０】
結論として、２つの画像ａ及びｂに対する提案されたフィルタは、ボクセルに固有の正規化ファクタである
【０１１１】
【数４９】
2016506267.tif 000051
によって
【０１１２】
【数５０】
2016506267.tif 000052
として記述でき、ｐは式（１４）の尤度であり、ｗｄ（ｉ，ｊ）はバイラテラルフィルタにおける通常の距離重み付け（式（５）及び（６）において
【０１１３】
【数５１】
2016506267.tif 000053
として参照される）である。２つの和は、ボクセルｉの近傍における全てのボクセルに対して実行される。
【０１１４】
要約すると、画像の既知のバイラテラルフィルタリングは、ローカルな平均化のためのボクセルに固有のウェイトを利用する。
【０１１５】
【数５２】
2016506267.tif 000054
ここで、ｗｄは距離に依存する重み付け（ガウスなど）であり、σは画像ａにおける（ローカルな）ノイズレベルである。
【０１１６】
画像ａの既知のジョイントバイラテラルフィルタリングは、ローカルな平均化のためのボクセルに固有のウェイトを利用する。
【０１１７】
【数５３】
2016506267.tif 000055
ここで、本発明との相違を明確に示すため、式（４）においてすでに説明された類似度の指標Ｌ１及びＬ２としてノイズレベルσａ及びσｂによるシンプルなガウス類似指標が慎重に利用される。
【０１１８】
他方、画像ａの提案されたバイラテラルフィルタリングは、全ての画像からのエッジに関する情報とフル２次ノイズ情報とを利用して、
【０１１９】
【数５４】
2016506267.tif 000056
を取得する。ここで、Δａ及びΔｂはそれぞれａｉ−ａｊ及びｂｉ−ｂｊの短縮形であり、Ｓｉｊは画像サンプル（ａｉ，ａｊ，ｂｉ，ｂｊ）のノイズ共分散行列の逆行列である。
【０１２０】
本発明が図面及び上記説明において詳細に図示及び説明されたが、このような図示及び説明は、例示的又は一例であり、限定的でないとみなされるべきであり、本発明は開示された実施例に限定されるものでない。開示された実施例に対する他の変形は、図面、開示及び添付した請求項の考察から請求された発明を当業者が実施することによって理解及び実現可能である。
【０１２１】
請求項において、“有する”という単語は他の要素又はステップを排除せず、“ある”という不定冠詞は複数を排除するものでない。単一の要素又は他のユニットは、請求項に記載された複数のアイテムの機能を実現してもよい。特定の手段が互いに異なる従属形式の請求項に記載されているという事実は、これらの手段の組み合わせが効果的には利用可能でないことを示すものでない。
【０１２２】
コンピュータプログラムは、他のハードウェアと一緒に又は一部として供給される光記憶媒体又はソリッドステート媒体などの適切な非一時的な媒体に格納／配布されてもよいが、また、インターネットや他の有線若しくは無線通信システムなどを介し他の形態により配布されてもよい。
【０１２３】
請求項における参照符号は、範囲を限定するものとして解釈されるものでない。
【図１】
2016506267.tif 000062
【図２】
2016506267.tif 000063
【図３】
2016506267.tif 000064
【図４】
2016506267.tif 000065
【図５】
2016506267.tif 000066
【図６】
2016506267.tif 000067
【手続補正書】
【提出日】平成27年6月23日(2015.6.23)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
画像をフィルタリングする画像処理装置であって、
同じオブジェクトの第１及び第２画像を取得する画像入力であって、前記第１及び第２画像は複数のボクセルを有し、ノイズ共分散によって相関し、各ボクセルは信号値とノイズ値とを含むボクセル値を有する、画像入力と、
前記第１画像をフィルタリングするジョイントバイラテラルフィルタであって、前記第１画像の第１ボクセルは相対的なボクセルに固有のウェイトを含むフィルタ関数によりフィルタリングされ、前記第１画像の第１ボクセルの信号値が前記第１画像の第２ボクセルの信号値と同一であり、前記第２画像の第１ボクセルの信号値が前記第２画像の第２ボクセルの信号値と同一であると仮定すると、前記ウェイトは、前記第１画像における前記第１ボクセル及び第２ボクセルのボクセル値と、前記第２画像における第１ボクセル及び第２ボクセルのボクセル値とを取得する尤度を含み、前記第１画像のフィルタリングされた第１ボクセル値は、
【数１】
2016506267.tif 000068
として取得され、Ｓは共分散行列の逆行列であり、ａ１及びａ２は前記第１画像の第１及び第２ボクセル値であり、ｂ１及びｂ２は前記第２画像の第１及び第２ボクセル値であり、Δａ＝ａ１−ａ２，Δｂ＝ｂ１−ｂ２であり、ｗｄはバイラテラルフィルタリングにおける距離重み付けであり、前記共分散行列は、
【数２】
2016506267.tif 000069
又は
【数３】
2016506267.tif 000070
の形式を有し、対角要素は前記ボクセル値のノイズ分散を記述し、非対角要素は前記画像内又は前記画像間のノイズの相関を記述する、ジョイントバイラテラルフィルタと、
前記フィルタリングされた画像を提供する画像出力と、
を有する画像処理装置。
【請求項２】
前記フィルタは、前記ノイズ共分散を記述するノイズ共分散行列を含む尤度を利用するよう構成される、請求項１記載の画像処理装置。
【請求項３】
前記フィルタは、
【数４】
2016506267.tif 000071
の形式の尤度を利用するよう構成される、請求項２記載の画像処理装置。
【請求項４】
前記画像入力は、前記オブジェクトの異なる性質を示すように、前記第１画像と前記第２画像とを提供するよう構成される、請求項１記載の画像処理装置。
【請求項５】
前記画像入力は、スペクトルＣＴ画像を提供するよう構成され、前記第１画像は前記オブジェクトの第１の物理的性質を示し、前記第２画像は前記オブジェクトの第２の物理的性質を示し、前記第１の物理的性質と前記第２の物理的性質とは互いに異なる、請求項４記載の画像処理装置。
【請求項６】
前記第１画像と前記第２画像との一方はコンプトン画像であり、前記第１画像と前記第２画像との他方は光電画像である、請求項４記載の画像処理装置。
【請求項７】
前記第１画像と前記第２画像との一方はコンプトン画像又は光電画像であり、前記第１画像と前記第２画像との他方はコンプトン画像と光電画像との組み合わせである、請求項６記載の画像処理装置。
【請求項８】
前記画像入力は、微分位相コントラストＣＴ画像を提供するよう構成され、前記第１画像と前記第２画像との一方は前記オブジェクトの吸収を示し、前記第１画像と前記第２画像との他方は前記オブジェクトにより誘導される前記オブジェクトをトラバースする放射線の位相シフトを示す、請求項１記載の画像処理装置。
【請求項９】
画像をフィルタリングする画像処理方法であって、
同じオブジェクトの第１及び第２画像を取得するステップであって、前記第１及び第２画像は複数のボクセルを有し、ノイズ共分散によって相関し、各ボクセルは信号値とノイズ値とを含むボクセル値を有する、取得するステップと、
前記第１画像をジョイントバイラテラルフィルタリングするステップであって、前記第１画像の第１ボクセルは相対的なボクセルに固有のウェイトを含むフィルタ関数によりフィルタリングされ、前記第１画像の第１ボクセルの信号値が前記第１画像の第２ボクセルの信号値と同一であり、前記第２画像の第１ボクセルの信号値が前記第２画像の第２ボクセルの信号値と同一であると仮定すると、前記ウェイトは、前記第１画像における前記第１ボクセル及び第２ボクセルのボクセル値と、前記第２画像における第１ボクセル及び第２ボクセルのボクセル値とを取得する尤度を含み、前記第１画像のフィルタリングされた第１ボクセル値は、
【数５】
2016506267.tif 000072
として取得され、Ｓは共分散行列の逆行列であり、ａ１及びａ２は前記第１画像の第１及び第２ボクセル値であり、ｂ１及びｂ２は前記第２画像の第１及び第２ボクセル値であり、Δａ＝ａ１−ａ２，Δｂ＝ｂ１−ｂ２であり、ｗｄはバイラテラルフィルタリングにおける距離重み付けであり、前記共分散行列は、
【数６】
2016506267.tif 000073
又は
【数７】
2016506267.tif 000074
の形式を有し、対角要素は前記ボクセル値のノイズ分散を記述し、非対角要素は前記画像内又は前記画像間のノイズの相関を記述する、ジョイントバイラテラルフィルタリングするステップと、
前記フィルタリングされた画像を提供するステップと、
を有する画像処理方法。
【請求項１０】
コンピュータ上で実行されると、前記コンピュータに請求項９記載の方法のステップを実行させるためのプログラムコード手段を有するコンピュータプログラム。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理装置及び画像をフィルタリングするための対応する方法に関する。
【背景技術】
【０００２】
患者又はオブジェクトの単一の画像のみでなく、正確に同じジオメトリにおけるそれらの２以上を提供する新規なイメージングモダリティが出現している。これの２つの具体例は、複素数の屈折率の実部及び虚部が測定される微分位相コントラストイメージング（ＤＰＣＩ；ＵＳ２０１２／００９９７０２Ａ１などに記載される）と、光電効果、コンプトン散乱及びおそらく造影剤によるＸ線の減衰が復元されるスペクトルＣＴ（ＵＳ２００８／０２５３５０４Ａなどに記載される）とである。
【０００３】
ＷＯ２０１１／１４５０４０Ａ１では、ジョイントバイラテラルフィルタリングが、画像ベースノイズ除去のための周知のバイラテラルフィルタの一般化として提案された。特に、オブジェクトの第１画像及び同一のオブジェクトの第２画像を提供する画像提供ユニットと、第１画像の画像値の間の第１類似度と、第２画像の画像値の間の第２類似度とに依存して第１画像をフィルタリングするフィルタリングユニットとを有する画像処理装置が説明されている。これは、オブジェクトが人間又は動物である場合、ノイズのために第１画像及び第２画像の一方の画像値が乱れたとしても、画像値がオブジェクトの同一部分、例えば、オブジェクトの組織又は骨物質に属する確率に依存して第１画像をフィルタリングすることを可能にし、これにより、例えば、ＣＴ画像のエッジ保存画像ベースノイズ除去などを向上させる。
【０００４】
ＵＳ２０１１／０６４２９２Ａ１は、バーチャル非コントラスト画像をエンハンスするための方法を開示し、デュアルスキャンＣＴ画像のペアを受信し、既知の組織減衰係数を利用してＣＴ画像のペアからバーチャル非コントラスト画像を計算することを含む。条件付き確率分布が、同一のタイプであるとして、ＣＴ画像のペアとバーチャル非コントラスト画像とのそれぞれにおける第１及び第２ポイントの組織について推定される。組織のための条件付き確率分布は、異なるタイプを有するとして、ＣＴ画像のペアとバーチャル非コントラスト画像とのそれぞれにおける第１及び第２ポイントにおいて推定される。同一のタイプであるとして第１及び第２ポイントの組織の事後確率は、条件付き確率分布から計算され、エンハンスされたバーチャル非コントラスト画像が、同一のタイプであるとして第１及び第２ポイントにおける組織の事後確率を利用して計算される。
【０００５】
しかしながら、既知の方法及び装置は、各種画像におけるノイズが相関しうることを認めていない。特に、ＤＰＣＩについては、位相コントラスト画像内のトランスアキシャル画像において強い空間的な相関があり（低周波数においてピークを有する周知のノイズパワースペクトルに等しい）、スペクトルＣＴについては、フォト画像とコンプトン画像とにおける同じ空間位置においてノイズの強力な反相関がある。
【０００６】
同じオブジェクトの２つの画像が利用可能であるが、異なる空間解像度においてノイズにおける異なる空間的相関によって反映される他の状況がまたあってもよい。具体例は、ある状況では、解剖学的（Ｔ２加重化など）及びパラメータ的画像（Ｔ１）が取得される磁気共鳴イメージング（ＭＲＩ）であってもよい。
【発明の概要】
【０００７】
画像をフィルタリングするための改良された画像処理装置及び方法と、特にジョイント画像ベースノイズ除去のための方法とを提供することが、本発明の課題である。
【０００８】
本発明の第１の態様では、画像をフィルタリングする画像処理装置であって、
同じオブジェクトの第１及び第２画像を取得する画像入力であって、前記第１及び第２画像は複数のボクセルを有し、ノイズ共分散によって相関し、各ボクセルは信号値とノイズ値とを含むボクセル値を有する、画像入力と、
前記第１画像をフィルタリングするジョイントバイラテラルフィルタであって、前記第１画像の第１ボクセルは相対的なボクセルに固有のウェイトを含むフィルタ関数によりフィルタリングされ、前記第１画像の第１ボクセルの信号値が前記第１画像の第２ボクセルの信号値と同一であり、前記第２画像の第１ボクセルの信号値が前記第２画像の第２ボクセルの信号値と同一であると仮定すると、前記ウェイトは、前記第１画像における前記第１ボクセル及び第２ボクセルのボクセル値と、前記第２画像における第１ボクセル及び第２ボクセルのボクセル値とを取得する尤度を含み、前記第１画像のフィルタリングされた第１ボクセル値は、
【０００９】
【数１】
2016506267.tif 000075
として取得され、Ｓは共分散行列の逆行列であり、ａ１及びａ２は前記第１画像の第１及び第２ボクセル値であり、ｂ１及びｂ２は前記第２画像の第１及び第２ボクセル値であり、Δａ＝ａ１−ａ２，Δｂ＝ｂ１−ｂ２であり、ｗｄはバイラテラルフィルタリングにおける距離重み付けであり、前記共分散行列は、
【００１０】
【数２】
2016506267.tif 000076
又は
【００１１】
【数３】
2016506267.tif 000077
の形式を有し、対角要素は前記ボクセル値のノイズ分散を記述し、非対角要素は前記画像内又は前記画像間のノイズの相関を記述する、ジョイントバイラテラルフィルタと、
前記フィルタリングされた画像を提供する画像出力と、
を有する画像処理装置が提供される。
【００１２】
本発明の更なる態様では、対応する画像処理方法が提供される。
【００１３】
本発明のまた更なる態様では、コンピュータ上で実行されると、コンピュータに処理方法のステップを実行させるためのプログラムコード手段を有するコンピュータプログラムと共に、プロセッサにより実行されると、ここに開示される方法を実行させるコンピュータプログラムプロダクトを格納する非一時的なコンピュータ可読記憶媒体が提供される。
【００１４】
本発明の好適な実施例は、従属形式の請求項に規定される。請求される方法、コンピュータプログラム及び媒体は、従属形式の請求項に規定され、請求される装置と類似及び／又は同一の好適な実施例を有すると理解される。
【００１５】
本発明によると、改良されたジョイントバイラテラルフィルタリングのためのアプローチが提案される。提案されたジョイントバイラテラルフィルタアプローチは、基本的には２つの項、距離重み付け及びレンジ重み付けを含む。この第１のものは、純粋に幾何学的に動機付けされ、変更されない。レンジフィルタは、本発明により変更されることが提案される画像位置に基づく重み付けを導入する。特に、画像間及び異なる画像位置間のノイズ相関が、改良されたジョイント画像ベースノイズ除去を提供するため、本発明により考慮される。
【００１６】
このコンテクストにおいて、“相関（ｉｎｔｅｒｒｅｌａｔｅｄ）”は、当該相関が同じ画像のボクセル間の相関に限定可能であるか、又は異なる画像における同じボクセル位置の間の相関に限定されてもよいように理解されるであろう。
【００１７】
好適な実施例では、前記フィルタは、前記ノイズ共分散を記述するノイズ共分散行列を含む尤度を利用するよう構成される。共分散を利用することによって、ノイズがオン方向に高く相関するが（ＤＰＣＩのためのイン・プレーン（ｉｎ−ｐｌａｎｅ）など）、ｚ方向には大きくは相関しないケースなどを解決することが可能になる。このような状況では、レンジフィルタの定数はイン・プレーン又はアウト・オブ・プレーン（ｏｕｔ−ｏｆ−ｐｌａｎｅ）ノイズ特性のみに調整可能であるため、従来の（ジョイント）バイラテラルフィルタリングは良好に実行されない。
【００１８】
本発明によると、前記フィルタは、
【００１９】
【数４】
2016506267.tif 000078
又は
【００２０】
【数５】
2016506267.tif 000079
の形式の共分散を利用するよう構成される。これらの共分散行列は効果的であると示した。第２共分散行列は扱うのがより容易である。
【００２１】
以下の特別なケースが特に興味がある。微分位相コントラストイメージングについて、再構成された屈折率の虚部（第１画像）と実部（第２画像）との間にはほとんどノイズ相関がなく、
【００２２】
【数６】
2016506267.tif 000080
を導き（用語はＣの一般的なバージョンに関し、すなわち、上述された２つのバージョンの上方の共分散行列）、一方、屈折率の実部においてトランスアキシャル画像内に強力なノイズ相関があり、すなわち、ｃ３４＞０である。共分散行列は、もちろん、ボクセルポジションに依存し、例えば、第２画像の第２ボクセルが異なるトランスアキシャル画像にある場合、ノイズは無相関であり、すなわち、
【００２３】
【数７】
2016506267.tif 000081
である。
【００２４】
スペクトルイメージングについて、各画像内にノイズの空間的な相関はほとんどなく（すなわち、コンプトン及び光電画像）、
【００２５】
【数８】
2016506267.tif 000082
を導き、また、異なる画像における異なる位置の間には相関がなく、すなわち、
【００２６】
【数９】
2016506267.tif 000083
であるが、同じ位置において２つの画像における画像値の間には強力な負の相関があり、ｃ１３＜０，ｃ２４＜０である。
【００２７】
さらに、実施例では、前記フィルタは、
【００２８】
【数１０】
2016506267.tif 000084
の形式の尤度を利用するよう構成され、Ｓは前記共分散行列の逆行列であり、ａ１及びａ２は前記第１画像の第１及び第２ボクセル値であり、ｂ１及びｂ２は前記第２画像の第１及び第２ボクセル値であり、Δａ＝ａ１−ａ２，Δｂ＝ｂ１−ｂ２である。
【００２９】
さらに、第１画像のフィルタリングされた第１ボクセル値は、
【００３０】
【数１１】
2016506267.tif 000085
として取得され、ｗｄは標準的な（ジョイント）バイラテラルフィルタリングにおける（通常の）距離ウェイティングである。距離｜ｉ−ｊ｜の増加により単調減少する何れかの関数、例えば、
【００３１】
【数１２】
2016506267.tif 000086
が利用されてもよい。
【００３２】
実際的な状況では、第１及び第２画像の少なくとも１つは非白色ノイズを有し、第１画像と第２画像との間にはノイズ相関が存在し、及び／又は第１画像及び／又は第２画像の第１及び第２ボクセル値におけるノイズ共分散は同じである。
【００３３】
好ましくは、同じフィルタが第２画像に適用される。さらに、提案されたアイデアは、３つ（又はより多くの）画像による状況に直接的な方法で拡張可能である。
【００３４】
他の実施例によると、前記画像入力は、前記オブジェクトの異なる性質を示すように、前記第１画像と前記第２画像とを提供するよう構成される。好ましくは、前記画像入力は、スペクトルＣＴ画像を提供するよう構成され、前記第１画像は前記オブジェクトの第１の物理的性質を示し、前記第２画像は前記オブジェクトの第２の物理的性質を示し、前記第１の物理的性質と前記第２の物理的性質とは互いに異なる。
【００３５】
例えば、前記第１画像と前記第２画像との一方はコンプトン画像であり、前記第１画像と前記第２画像との他方は光電画像である。他の実施例では、前記第１画像と前記第２画像との一方はコンプトン画像又は光電画像であり、前記第１画像と前記第２画像との他方はコンプトン画像と光電画像との組み合わせである。
【００３６】
またさらに、前記画像入力は、微分位相コントラストＣＴ画像を提供するよう構成され、前記第１画像と前記第２画像との一方は前記オブジェクトの吸収を示し、前記第１画像と前記第２画像との他方は前記オブジェクトにより誘導される前記オブジェクトをトラバースする放射線の位相シフトを示す。
【図面の簡単な説明】
【００３７】
本発明の上記及び他の態様は、後述される実施例を参照して明らかになるであろう。
【図１】図１は、本発明による画像処理装置の実施例を概略的及び例示的に示す。
【図２】図２は、画像入力ユニットの実施例を概略的及び例示的に示す。
【図３】図３は、多色放射ソースの発光スペクトルを概略的及び例示的に示す。
【図４】図４は、表示されたエネルギー範囲内のｋエッジによる物質の光電効果、コンプトン効果及びトータルの減衰の断面図のエネルギー依存を概略的及び例示的に示す。
【図５】図５は、微分位相コントラストＣＴシステムの実施例の要素を概略的及び例示的に示す。
【図６】図６は、画像を処理する画像処理方法の実施例を示すフローチャートを例示的に示す。
【発明を実施するための形態】
【００３８】
図１は、本発明による画像処理装置１の実施例を概略的に示す。それは、同一のオブジェクトの第１及び第２画像を取得する画像ユニット２を有し、第１及び第２画像は複数のボクセルを有し、ノイズ共分散により相関し、各ボクセルは信号値及びノイズ値を有するボクセル値を有する。装置１は更に、第１画像をフィルタリングするバイラテラルフィルタ３を有し、第１画像の第１ボクセルは相対的なボクセルに固有のウェイトを含むフィルタ関数によりフィルタリングされる。第１画像の第１ボクセルの信号値が第２画像の第１ボクセルの信号値と同じであり、第１画像の第２ボクセルの信号値が第２画像の第２ボクセルの信号値と同じであると仮定すると、当該ウェイトは、当該第１ボクセルと第２画像における第２ボクセルとのボクセル値を取得し、第２画像における第１ボクセルと第２ボクセルとを取得する尤度（ｌｉｋｅｌｉｈｏｏｄ）を含む。装置１は更に、フィルタリングされた画像を表示するためのディスプレイなどであってもよい、フィルタリングされた画像を提供する画像出力４を有する。
【００３９】
第１画像と第２画像とは、例えば、オブジェクトのスライスのみなど、オブジェクト全体又はオブジェクトの一部のみとすることが可能な同一のオブジェクトを示す２次元、３次元又は４次元画像とすることができる。オブジェクトが人間の心臓などの挙動するオブジェクトである場合、第１画像と第２画像とは、好ましくは、同じ挙動状態におけるオブジェクトを示す。本実施例の第１画像及び第２画像は、プロジェクションデータに基づき、特に、ＣＴプロジェクションデータに基づき再構成され、ここで、第１画像及び第２画像は同一のプロジェクションデータから再構成される。
【００４０】
画像入力ユニット２は、それらがオブジェクトの異なる性質、特に、異なる物理的性質を示すように、第１画像と第２画像とを提供するよう構成される。実施例では、画像入力ユニット２は、スペクトルＣＴ（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）画像を提供するよう構成され、第１画像はオブジェクトの第１の物理的性質を示し、第２画像はオブジェクトの第２の物理的性質を示し、第１の物理的性質と第２の物理的性質とは互いに異なる。例えば、第１画像と第２画像との一方はコンプトン画像とすることができ、第１画像と第２画像との他方は光電画像とすることができ、あるいは、第１画像と第２画像との一方はコンプトン画像又は光電画像とすることができ、第１画像と第２画像との他方はコンプトン画像と光電画像との組み合わせとすることができる。
【００４１】
画像入力ユニット２は格納ユニットとすることが可能であり、そこでは、第１画像と第２画像とが格納されるか、あるいは、画像入力ユニット２はまた、データ接続を介し第１画像及び／又は第２画像を受信することを可能にする画像受信ユニットとすることができる。画像入力ユニット２はまた、第１画像と第２画像とを生成するスペクトルＣＴシステムなどのイメージングモダリティとすることができる。このようなスペクトルＣＴシステムは、以下において図２を参照して概略的及び例示的に説明される。
【００４２】
スペクトルＣＴシステムは、ｚ方向にパラレルに延びる回転軸Ｒの周りに回転可能なガントリ１４を含む。Ｘ線チューブなどの放射ソース１５は、ガントリ１４に搭載される。本実施例では、放射ソース１５は多色放射線を発する。放射ソース１５には、放射ソース１５により発せられた放射から円錐放射ビーム１７を形成するコリメータ装置１６が備えられる。他の実施例では、コリメータ装置１６は、扇形形状などを有する他の形状を有する放射ビームを形成するよう構成可能である。
【００４３】
放射線は、円筒形の検査ゾーン５における関心領域の患者又は技術的オブジェクトなどのオブジェクト（図２に図示せず）をトラバースする。関心領域をトラバースした後、放射ビーム１７は、本実施例では、２次元検出面を有する検出装置６に入射し、ここで、検出装置６はガントリ１４に搭載される。他の実施例では、検出装置６は１次元検出面を有することが可能である。
【００４４】
検出装置６は、好ましくは、エネルギー分解検出要素を有する。エネルギー分解検出要素は、好ましくは、例えば、入射光子をカウントする原理などにより動作し、異なるエネルギーウィンドウにおける光子数を示す信号を出力する光子計数検出要素である。このようなエネルギー分解検出装置は、例えば、Ｌｌｏｐａｒｔ  Ｘ．，ｅｔ  ａｌ．，“Ｆｉｒｓｔ  ｔｅｓｔ  ｍｅａｓｕｒｅｍｅｎｔｓ  ｏｆ  ａ  ６４ｋ  ｐｉｘｅｌ  ｒｅａｄｏｕｔ  ｃｈｉｐ  ｗｏｒｋｉｎｇ  ｉｎ  ａ  ｓｉｎｇｌｅ  ｐｈｏｔｏｎ  ｃｏｕｎｔｉｎｇ  ｍｏｄｅ”，Ｎｕｃｌ．Ｉｎｓｔ．ａｎｄ  Ｍｅｔｈ．Ａ，５０９（１−３）：１５７−１６３，２００３及びＬｌｏｐａｒｔ，Ｘ．，ｅｔ  ａｌ．，“Ｍｅｄｉｐｉｘ２：Ａ  ６４−ｋ  ｐｉｘｅｌ  ｒｅａｄｏｕｔ  ｃｈｉｐ  ｗｉｔｈ  ５５  μｍ  ｓｑｕａｒｅ  ｅｌｅｍｅｎｔｓ  ｗｏｒｋｉｎｇ  ｉｎ  ａ  ｓｉｎｇｌｅ  ｃｏｕｎｔｉｎｇ  ｍｏｄｅ”，ＩＥＥＥ  Ｔｒａｎｓ．Ｎｕｃｌ．Ｓｃｉ．４９（５）：２２７９−２２８３，２００２に記載されている。好ましくは、各検出要素が少なくとも２つの異なるエネルギーウィンドウのための少なくとも２つのエネルギー分解された検出信号を提供するように、エネルギー分解検出要素が構成される。しかしながら、イメージングシステムの感度及びノイズロウバストネスを向上させるため、さらに高いエネルギー分解能を有することが好ましい。
【００４５】
ガントリ１４は、好ましくは、モータ７によって一定であるが、調整可能な角速度により駆動される。更なるモータ８は、検査ゾーン５の患者テーブルに配置された患者などのオブジェクトを、回転軸方向又はｚ軸にパラレルにずらすため設けられる。放射ソース１５と検査ゾーン５、特に関心領域とが螺旋の軌跡に沿って互いに対して移動するように、これらのモータ７，８は制御ユニット９などにより制御される。また、オブジェクト又は検査ゾーン５、特に関心領域は移動せず、放射ソース１５が回転し、すなわち、放射ソース１５が関心領域に対して円の軌跡に沿って移動することが可能である。検出装置６により取得されるデータは、本実施例では、提供されたエネルギーに依存した検出データから関心領域の第１及び第２画像を少なくとも再構成する再構成ユニット１０に提供されるエネルギーに依存した検出データである。また、再構成ユニット１０は、好ましくは制御ユニット９により制御される。再構成された画像は、少なくとも第１画像をフィルタリングするフィルタリングユニット３に提供される。
【００４６】
再構成ユニット１０は、検出データから少なくとも２つの減衰コンポーネントを決定する計算ユニット１２を有する。検出データの減衰コンポーネントは、１つのみ又は複数の物理的効果によって引き起こされ、例えば、ある減衰コンポーネントはコンプトン効果により引き起こされるコンポーネントでありうるものであり、他の減衰コンポーネントは光電効果により引き起こされるコンポーネントでありうる。更なる減衰コンポーネントは関心領域に存在するｋエッジにより引き起こされるコンポーネントでありうる。あるいは、又はさらに、減衰コンポーネントは、例えば、関心領域内のある物質の吸収などによって引き起こされるコンポーネントでありうる。例えば、減衰コンポーネントはある物質の吸収により引き起こされるコンポーネントでありうるものであり、他の減衰コンポーネントは他の物質の吸収により引き起こされるコンポーネントでありうる。
【００４７】
再構成ユニット１０は更に、対応する減衰コンポーネント画像が再構成されるように、検出データの計算された減衰コンポーネントをバックプロジェクションするバックプロジェクションユニット１３を有する。例えば、コンプトン画像は、第１画像と第２画像との一方として再構成可能であり、光電画像は、第１画像と第２画像との他方として再構成可能である。
【００４８】
計算ユニット１２への入力は、複数のエネルギービン、本実施例では、最小で３つのエネルギービンのエネルギー分解検出データｄｉである。これらの検出データｄｉは、第ｉエネルギービンｅｉのスペクトル感度Ｄｉ（Ｅ）を示す。さらに、多色放射ソース１５の発光スペクトルＴ（Ｅ）は、一般に既知であるか、又は測定可能である。多色放射ソースのこのような発光スペクトルＴ（Ｅ）の具体例は、図３において概略的及び例示的に示される。計算ユニット１２では、検出データｄｉの生成は、スペクトルによる光電効果Ｐ（Ｅ）、スペクトルによるコンプトン効果Ｂ（Ｅ）及びスペクトルによるｋエッジＫ（Ｅ）の線形結合としてモデル化される。
【００４９】
スペクトルＰ（Ｅ）、Ｂ（Ｅ）及びＫ（Ｅ）は、図４において例示的及び概略的に示される。
【００５０】
検出データの生成は、以下の連立方程式によりモデル化できる。
【００５１】
【数１３】
2016506267.tif 000087
ここで、ρｐｈｏｔｏ，ρＣｏｍｐｔｏｎ及びρｋ−ｅｄｇｅはそれぞれ、光電コンポーネント、コンプトンコンポーネント及びｋエッジコンポーネントの密度長積（ｄｅｎｓｉｔｙ  ｌｅｎｇｔｈ  ｐｒｏｄｕｃｔ）である。
【００５２】
少なくとも３つの検出信号ｄ１，ｄ２，ｄ３が少なくとも３つのエネルギービンｅ１，ｅ２，ｅ３について利用可能であるため、３つの密度長積である３つの未知を有する少なくとも３つの方程式の連立方程式が構成され、従って、計算ユニット１２において既知の数値方法により解くことができる。３より多くのエネルギービンが利用可能である場合、測定結果のノイズ統計を考慮する最大尤度アプローチを利用することが好ましい。一般に、３つのエネルギービンで十分である。しかしながら、感度及びノイズロウバストネスを向上させるため、より多くのエネルギービンについてより多くの検出データを有することが好ましい。各減衰コンポーネントについて、すなわち、各密度長積について、オブジェクトの減衰コンポーネント画像は、各密度長積をバックプロジェクションすることによって再構成可能である。あるいは、又はさらに、異なる密度長積の画像、又はバックプロジェクション前の異なる密度長積の画像が、合成された減衰コンポーネントにより生成されるオブジェクトの画像である合成画像を生成するため合成可能である。
【００５３】
他の実施例では、画像入力ユニット２は、微分位相コントラストＣＴ画像を提供するよう構成可能であり、第１画像及び第２画像の一方はオブジェクトの吸収を示し、第１画像及び第２画像の他方はオブジェクトにより誘導される、オブジェクトをトラバースする放射位相シフトを示す。特に、第１画像は、好ましくは、オブジェクトのトータルの減衰係数を示し、第２画像は、好ましくは、オブジェクトの屈折率の実部を示す。また、これらの第１及び第２画像は、画像入力ユニット２であってもよい格納ユニットに格納可能であるか、又は画像入力ユニット２であってもよい画像受信ユニットを介し受信可能である。しかしながら、画像入力ユニット２はまた、第１画像及び第２画像を生成し、図５を参照して以下で例示的及び概略的に説明される微分位相コントラストＣＴシステムとすることができる。
【００５４】
図５は、微分位相コントラストＣＴシステムの一例となる実施例の３次元表現を示す。やや大きなＸ線ソース３２はソースグレイティング３４に隣接して配置される。Ｘ線ソース３２は発せられる放射の波長に関するそれのサイズのため非干渉であると考えられるため、ソースグレイティング３４は、複数の単独のコヒーラントなＸ線ソースを提供するのに利用される。
【００５５】
Ｘ線放射３５は、Ｘ線ソース３２から光軸３７の方向に発せられ、おそらくＸ線の扇形ビーム又は円錐型ビームを構成する。図５において、Ｘ線ビームの各形状は示されていない。Ｘ線放射３５はオブジェクト３６に到来し、オブジェクト３６を貫通し、その後にビームスプリッタグレイティング３８に到来する。ビームスプリッタグレイティング３８のトレンチ又はギャップは、ビームスプリッタグレイティングの固体エリア、すなわち、ブロッキング領域に関して電磁放射を通過させる位相を変更する。従って、φによる、特にπによる位相シフトが実行される。
【００５６】
解析グレイティング４０は、ビームスプリッタグレイティング３８とＸ線検出装置４２との間に配置される。Ｘ線検出装置の方向にビームスプリッタグレイティング３８から発せられる複数の波は、解析グレイティング４０に到来し、その後、Ｘ線検出装置４２の表面に強度変調パターンを生成する。
【００５７】
解析グレイティング４０に対してビームスプリッタグレイティング３８をシフトすることによって、互いに対して、特にビームスプリッタグレイティング３８又は解析グレイティング４０のグレイティング期間の一部によってグレイティングをずらすことは、画像検出装置４２により取得可能である。従って、送信画像の取得中、いわゆる、位相ステッピングが実行され、ここでは、グレイティング３８，４０の相対位置が、１つの期間において４〜１０のステップなどにおいて変更され、各相対位置について送信画像が測定される。この送信画像系列から、従来の吸収画像が、各検出ピクセルについて測定された強度を単純平均することによって取得される。従来のＸ線送信イメージングと同様に、この測定はＢｅｅｒの法則によってＸ線パスに沿ってトータルの減衰係数に対する線積分に関連する。屈折率の実部に関する情報は、位相ステッピングサイクルにおける測定された強度のダイナミクスの解析により取得される。当該強度は、好ましくは、解析グレイティング４０のグレイティング期間により定期的に変動する。典型的には、正弦波変動が想定できる。サイナスの位相はビームスプリッタグレイティング３８にヒットする位相フロントの勾配に関連する。物理学によって、位相フロントの当該勾配は、オブジェクトの屈折率の実部のプロジェクションの第１導関数に線形に相関する。従って、オブジェクトの屈折率の実部のプロジェクションは、サイナスの位相に基づき決定可能である。オブジェクトの屈折率の実部の決定されたプロジェクションをバックプロジェクションすることによって、オブジェクトの屈折率の実部を示す画像が再構成可能である。微分位相コントラストＣＴシステムは、例えば、オブジェクトのトータルの減衰係数を示す第１画像とオブジェクトの屈折率の実部を示す第２画像とを再構成可能である。既知の気分いそうコントラストＣＴ技術のより詳細な説明は、例えば、参照することによりここに援用される、Ｔｉｍｍ  Ｗｅｉｔｋａｍｐ  ｅｔ  ａｌ．による論文“Ｘ−ｒａｙ  ｐｈａｓｅ  ｉｍａｇｉｎｇ  ｗｉｔｈ  ａ  ｇｒａｔｉｎｇ  ｉｎｔｅｒｆｅｒｏｍｅｔｅｒ”，ＯＰＴＩＣＳ  ＥＸＰＲＥＳＳ，ｖｏｌ．１３，ｎｏ．１６，ｐａｇｅｓ  ６２９６−６３０４（Ａｕｇｕｓｔ  ２００５）などに開示される。
【００５８】
ＷＯ２０１１／１４５０４０Ａ１には、ジョイントバイラテラルフィルタリングが、画像ベースノイズ除去のための周知なバイラテラルフィルタの一般化として提案された。このアプローチでは、フィルタリングユニットは、好ましくは、第１画像が以下のフィルタリング方程式を利用することによってノイズ除去されるように、フィルタリングを実行するよう構成される。
【００５９】
【数１４】
2016506267.tif 000088
ここで、
【００６０】
【数１５】
2016506267.tif 000089
は第１画像のポジション
【００６１】
【数１６】
2016506267.tif 000090
における第１画像要素の平均画像値を示し、Δｉはポジション
【００６２】
【数１７】
2016506267.tif 000091
における第１画像要素を有する第１領域を示し、
【００６３】
【数１８】
2016506267.tif 000092
は第１領域Δｉ内の第１領域画像要素の位置を示し、
【００６４】
【数１９】
2016506267.tif 000093
は、ａ）第１領域画像要素
【００６５】
【数２０】
2016506267.tif 000094
の画像値と、任意的にはポジション
【００６６】
【数２１】
2016506267.tif 000095
における各第１領域要素を包囲する画像要素の画像値と、ｂ）ポジション
【００６７】
【数２２】
2016506267.tif 000096
における第１画像要素の画像値と、任意的には第１画像の第１画像要素を包囲する画像要素の画像値との間の第１類似度を示す。
【００６８】
関数
【００６９】
【数２３】
2016506267.tif 000097
は、第２画像の第１領域に対応する第２領域内のポジション
【００７０】
【数２４】
2016506267.tif 000098
における第２領域画像要素の画像値と、任意的にはポジション
【００７１】
【数２５】
2016506267.tif 000099
における第２領域要素を包囲する画像要素の画像値と、ｂ）第２画像の第１画像における同一ポジションにおける第１画像要素に対応するポジション
【００７２】
【数２６】
2016506267.tif 000100
における第２画像要素の画像値と、任意的には第２画像要素を包囲する画像要素の画像値との間の第２類似度を示す。関数
【００７３】
【数２７】
2016506267.tif 000101
は、ポジション
【００７４】
【数２８】
2016506267.tif 000102
における画像要素の間の距離に依存する距離関数であり、
【００７５】
【数２９】
2016506267.tif 000103
は第１画像におけるポジション
【００７６】
【数３０】
2016506267.tif 000104
における非フィルタリング画像値を示す。
【００７７】
式（２）に示される分母により除算される第１類似度Ｌ１、第２類似度Ｌ２及び距離関数ｗｄの合成は、式（２）により加重平均を実行するためのウェイトとしてみなすことができる。従って、式（２）に従って、フィルタリングユニット３は、好ましくは、第１画像の画像要素の画像値を加重平均するよう構成され、ここで、第１画像の第１画像要素の平均画像値を決定するため、第１画像要素を含む第１領域内の第１領域画像要素の画像値が重み付け及び平均化され、第１領域画像要素の画像値は、ａ）第１領域画像要素の画像値と、任意的には各第１画像要素を包囲する画像要素の画像値と、ｂ）第１画像要素の画像値と、第１画像の第１画像要素を包囲する画像要素の画像値との間の類似度である第１類似度に依存して、また、ａ）第２画像の第１領域に対応する第２領域内の第２領域画像要素の画像値と、任意的には各第２領域画像要素を包囲する画像要素の画像値と、ｂ）第２画像の第１画像要素に対応する第２画像要素の画像値と、任意的には第２画像要素を包囲する画像要素の画像値との間の類似度である第２類似度に依存して、ウェイトにより重み付けされる。
【００７８】
ベクトル
【００７９】
【数３１】
2016506267.tif 000105
は、
【００８０】
【数３２】
2016506267.tif 000106
を中心とする画像パッチとしてみなすことができる。例えば、各第１画像要素を包囲する画像要素の任意的な画像値と、各第２画像要素を包囲する画像要素の任意的な画像値とが考慮される場合、２次元のケースでは、画像パッチはポジション
【００８１】
【数３３】
2016506267.tif 000107
におけるピクセルを中心とする９×９のサブ画像であってもよい。当該パッチはまた他のサイズを有することも可能である。例えば、パッチは５×５又は７×７のサイズを有することが可能である。さらに、３次元パッチがまた利用可能である。任意的な画像値が考慮されない場合、すなわち、ポジション
【００８２】
【数３４】
2016506267.tif 000108
における単一の画像値の間の差分しか考慮されない場合、各ベクトルはそれぞれ、ポジション
【００８３】
【数３５】
2016506267.tif 000109
における単一の各画像値しか有さない。このとき、パッチサイズは１×１である。
【００８４】
フィルタリングが２より多くの画像に基づき実行される場合、フィルタリングユニット３は、好ましくは、以下のフィルタリング方程式に従って第１画像をフィルタリングするよう構成される。
【００８５】
【数３６】
2016506267.tif 000110
ここで、Ｎは第１画像をフィルタリングするため考慮される画像数を示し、インデックスｋは、各要素がｋ番目の画像に属することを示す。
【００８６】
第ｋ類似度は、以下の方程式により定義できる。
【００８７】
【数３７】
2016506267.tif 000111
ここで、
【００８８】
【数３８】
2016506267.tif 000112
はそれぞれ、
【００８９】
【数３９】
2016506267.tif 000113
におけるローカルノイズ推定値を示す。
【００９０】
ノイズ推定値は、実際に同質な領域であると考えられる領域内の標準偏差の決定によって、手動により取得可能である。また、既知の初期的な生データにおけるノイズから始めて良好な近似までの古典的なノイズ伝搬によって空間的に分解されたノイズ推定値を取得することが可能である。既知の古典的なノイズ伝搬は、例えば、参照することによりここに援用される、Ａｄａｍ  Ｗｕｎｄｅｒｌｉｃｈ  ｅｔ  ａｌ．，による論文“Ｉｍａｇｅ  ｃｏｖａｒｉａｎｃｅ  ａｎｄ  ｌｅｓｉｏｎ  ｄｅｔｅｃｔａｂｉｌｉｔｙ  ｉｎ  ｄｉｒｅｃｔ  ｆａｎ−ｂｅａｍ  ｘ−ｒａｙ  ｃｏｍｐｕｔｅｄ  ｔｏｍｏｇｒａｐｈｙ”，Ｉｎｓｔｉｔｕｔｅ  ｏｆ  Ｐｈｙｓｉｃｓ  ａｎｄ  Ｅｎｇｉｎｅｅｒｉｎｇ  ｉｎ  Ｍｅｄｉｃｉｎｅ，ｐａｇｅｓ  ２４７１  ｔｏ  ２４９３（２００８）などに開示される。同一のコンセプトは、空間ＣＴだけでなく、微分位相コントラストＣＴに適用可能である。
【００９１】
ローカルノイズ推定値の代わりに、グローバルノイズ推定値がまた、特にノイズが画像全体で大きくは変動しない場合、利用可能である。
【００９２】
上述の方程式は、ａ）第１領域画像要素と第１画像要素との間の第１類似度と、ｂ）第２領域画像要素と第２画像要素との間の第２類似度との一方が減少する場合、かつ、ａ）第１領域画像要素と第１画像要素との間の第１類似度と、ｂ）第２領域画像要素と第２画像要素との間の第２類似度との他方が変化しないか、若しくは減少する場合、第１領域画像要素の画像値のウェイトは減少することを保証する。
【００９３】
これらの類似度は、特に、各画像がノイズ除去される場合、２つの画像値又は２つの画像パッチがどの程度同じであるかの推定値を提供する。従って、式（４）において定義される類似度はまた、画像値が類似する確率の一緒としてみなすことができる。ノイズよりはるかに小さな相違について、尤度又は確率は１に近く、類似度はノイズレベルより大きな相違についてゼロに低下する。特に、１×１より大きなサイズを有する画像パッチが考慮される場合、類似度のガウスプロトタイプのため、ポジション
【００９４】
【数４０】
2016506267.tif 000114
における画像値の結果としてのウェイトは、
【００９５】
【数４１】
2016506267.tif 000115
の周りの画像パッチがＮ個の画像の何れかにおいて統計的に異なる場合、小さくなるであろう。従って、第１画像における統計的に有意なエッジにわたるスムージングは実行されないことが保証できる。従って、フィルタリングはエッジ保存的である。
【００９６】
空間加重関数としてもみなすことができる距離関数は、以下の方程式により定義できる。
【００９７】
【数４２】
2016506267.tif 000116
ここで、パラメータσｗｄは、平均化中に考慮される近傍の有効サイズを制御する。
【００９８】
以下において、ここに提案されるようなジョイントバイラテラルフィルタリングのための更に向上されたアプローチが概略される。バイラテラルフィルタリングアプローチは、基本的には２つの項、距離重み付け（上記の式（５）及び（６）における“ｗｄ”として参照される）と、レンジ重み付け（上記の式（４）における“Ｌ”として参照される）とを含む。距離重み付けは、純粋に幾何学的に動機付けされ、上述した既知のアプローチと比較して変更されない。レンジフィルタは、上述された既知のアプローチと比較して本発明により変更されることが提案される画像値に基づく重み付けを導入する。
【００９９】
具体例は、２つの画像について説明される。バイラテラルフィルタリングの基本タスクは、２つの異なるピクセルが同じ組織クラスに属する確率、すなわち、基礎となる真のピクセル値が同じである尤度を推定することである。この推定は、２つの画像における同じ２つの異なる位置におけるサンプルである４つ（又はそれ以上）の画像サンプルの評価に基づく。これらのサンプルは、ａ１，ａ２，ｂ１，ｂ２として示され（ａｉ及びｂｉは上述した
【０１００】
【数４３】
2016506267.tif 000117
に対応し、これらは、よりコンパクトな表記を実現するため以下において利用される）、ベクトル（ａ１，ａ２，ｂ１，ｂ２）を構成する。これらのサンプルに対する２次統計量は、（新たに導入された）共分散行列Ｃ
【０１０１】
【数４４】
2016506267.tif 000118
によって、これらの測定結果におけるノイズを記述し、ここで、対角要素はサンプルのノイズ共分散を記述し、非対角要素は画像（ｃ１２，ｃ３４）内又は画像（ｃ１３，ｃ１４，ｃ２３，ｃ２４）間のノイズの相関を記述する。行列Ｃは対称的であり、正定値である。この共分散行列によって、測定されたベクトル（ａ１，ａ２，ｂ１，ｂ２）においてノイズ実現（ｎ１，ｎ２，ｎ３，ｎ４）を取得する尤度を記述できる。
【０１０２】
【数４５】
2016506267.tif 000119
この測定ではノイズと信号コンテンツとを区別することはできない。従って、ａ１＝ａ２及びｂ１＝ｂ２の仮定の下、どの程度の確率で実際の測定結果（ａ１，ａ２，ｂ１，ｂ２）を観察することになるかが調べられる。この仮定の下では、測定における信号コンテンツは差分Δａ＝ａ１−ａ２，Δｂ＝ｂ１−ｂ２を計算し、これらの差分を生じさせる（ｎ１，ｎ２，ｎ３，ｎ４）の全ての組み合わせに対して積分することによって除外できるため、これは所与の統計量に対して容易に回答できる。表記のコンパクト化のため、Ｓ＝Ｃ−１と設定される。所望の確率は、
【０１０３】
【数４６】
2016506267.tif 000120
に従って計算できる。この周辺確率は、ジョイントバイラテラルフィルタ方程式におけるジョイント尤度項として挿入される必要があり、
【０１０４】
【数４７】
2016506267.tif 000121
を生じさせる。
【０１０５】
最も一般的な共分散行列について、積分の明示的な計算はかなり複雑になる。しかしながら、共分散行列Ｃが数値的に既知になると、実際の測定結果の確率と逆行列Ｓとを計算することは容易な作業である。
【０１０６】
２つのシンプルなケースが考えられる。４つのサンプルのノイズが無相関である場合、共分散行列とそれの逆行列との双方は対角となり、所望される尤度はガウス関数の積となる。他の興味深いケースは、
【０１０７】
【数４８】
2016506267.tif 000122
の形式の共分散行列についてである。幾何学的には、この共分散行列は、２つの画像位置における画像共分散が同じであり（すなわち、ｃ１１＝ｃ２２，ｃ３３＝ｃ４４）、画像間の共分散が２つの画像位置（ｃ１３＝ｃ２４，ｃ１４＝ｃ２３）に関して対称的であるケースをモデル化する。典型的な状況では、これらの関係は、再構成されたＣＴ画像における分散は空間においてゆっくりとしか変化しないため、ほぼ充足される。
【０１０８】
この増加した対称性はまた、逆行列Ｓにおいて反映され、それは、
【０１０９】
【数４９】
2016506267.tif 000123
を読む。
【０１１０】
明らかに、共分散行列の対称性がまたそれの逆行列に存在し、従って、
【０１１１】
【数５０】
2016506267.tif 000124
として記述できる。この仮定の下、ベクトル（ａ１，ａ２，ｂ１，ｂ２）を測定するための所望の尤度は、
【０１１２】
【数５１】
2016506267.tif 000125
である。
【０１１３】
フィルタの式（１６）に挿入されるような式（１４）に与えられる結合尤度は、従来技術において行われるような２つの画像に対する導出された独立した類似度の積の形態を有さないが、追加的な項、すなわち、差分の積ΔａΔｂがあることに留意されたい。
【０１１４】
結論として、２つの画像ａ及びｂに対する提案されたフィルタは、ボクセルに固有の正規化ファクタである
【０１１５】
【数５２】
2016506267.tif 000126
によって
【０１１６】
【数５３】
2016506267.tif 000127
として記述でき、ｐは式（１４）の尤度であり、ｗｄ（ｉ，ｊ）はバイラテラルフィルタにおける通常の距離重み付け（式（５）及び（６）において
【０１１７】
【数５４】
2016506267.tif 000128
として参照される）である。２つの和は、ボクセルｉの近傍における全てのボクセルに対して実行される。
【０１１８】
要約すると、画像の既知のバイラテラルフィルタリングは、ローカルな平均化のためのボクセルに固有のウェイトを利用する。
【０１１９】
【数５５】
2016506267.tif 000129
ここで、ｗｄは距離に依存する重み付け（ガウスなど）であり、σは画像ａにおける（ローカルな）ノイズレベルである。
【０１２０】
画像ａの既知のジョイントバイラテラルフィルタリングは、ローカルな平均化のためのボクセルに固有のウェイトを利用する。
【０１２１】
【数５６】
2016506267.tif 000130
ここで、本発明との相違を明確に示すため、式（４）においてすでに説明された類似度の指標Ｌ１及びＬ２としてノイズレベルσａ及びσｂによるシンプルなガウス類似指標が慎重に利用される。
【０１２２】
他方、画像ａの提案されたバイラテラルフィルタリングは、全ての画像からのエッジに関する情報とフル２次ノイズ情報とを利用して、
【０１２３】
【数５７】
2016506267.tif 000131
を取得する。ここで、Δａ及びΔｂはそれぞれａｉ−ａｊ及びｂｉ−ｂｊの短縮形であり、Ｓｉｊは画像サンプル（ａｉ，ａｊ，ｂｉ，ｂｊ）のノイズ共分散行列の逆行列である。
【０１２４】
本発明が図面及び上記説明において詳細に図示及び説明されたが、このような図示及び説明は、例示的又は一例であり、限定的でないとみなされるべきであり、本発明は開示された実施例に限定されるものでない。開示された実施例に対する他の変形は、図面、開示及び添付した請求項の考察から請求された発明を当業者が実施することによって理解及び実現可能である。
【０１２５】
請求項において、“有する”という単語は他の要素又はステップを排除せず、“ある”という不定冠詞は複数を排除するものでない。単一の要素又は他のユニットは、請求項に記載された複数のアイテムの機能を実現してもよい。特定の手段が互いに異なる従属形式の請求項に記載されているという事実は、これらの手段の組み合わせが効果的には利用可能でないことを示すものでない。
【０１２６】
コンピュータプログラムは、他のハードウェアと一緒に又は一部として供給される光記憶媒体又はソリッドステート媒体などの適切な非一時的な媒体に格納／配布されてもよいが、また、インターネットや他の有線若しくは無線通信システムなどを介し他の形態により配布されてもよい。
【０１２７】
請求項における参照符号は、範囲を限定するものとして解釈されるものでない。
【国際調査報告】
2016506267.tif 000132
2016506267.tif 000133