(57)【要約】  画像訂正に使用する少なくとも１つのアーチファクトのテンプレートを生成するシステム及び方法が開示される。そのアーチファクトを含む画像は、少なくとも２つの一様な照射を用いて生成される。各照射は、異なるディテクタ動作温度で生成される。その画像における各ピクセル位置でのグレー値の局所分散が計算される。すると、画像における各ピクセルが分類される。その分類に基づき二値画像が生成される。その後、その二値画像とそのアーチファクトを含む画像データとの両方に基づき、テンプレートが形成される。
2009531109.tif 000002
【特許請求の範囲】
【請求項１】
画像訂正に使用する少なくとも１つのアーチファクトのテンプレートを生成する方法において、
少なくとも２つのゲインマップからピクセラライゼーション画像を生成するステップであって、各ゲインマップが、異なるディテクタ温度で生成される、ステップと、
前記ピクセラライゼーション画像での各ピクセル位置でグレー値の局所分散を計算するステップと、
前記ピクセラライゼーション画像における、ピクセル状の、ピクセラライゼーション領域及び非ピクセラライゼーション領域を決定するステップと、
決定されたピクセラライゼーション領域と非ピクセラライゼーション領域とに基づき、ピクセラライゼーション画像の二値画像を生成するステップと、
前記ピクセラライゼーション画像と前記二値画像とに基づきテンプレートを生成するステップとを有する、方法。
【請求項２】
前記局所分散を計算する前に、欠陥マップに基づき、前記ピクセラライゼーション画像における欠陥ピクセルの排除を行うステップを含む、請求項１に記載の方法。
【請求項３】
前記ピクセラライゼーション画像における欠陥訂正を実行するステップを含む、請求項１に記載の方法。
【請求項４】
前記テンプレートを生成するステップが、欠陥であると決定されるか、又は非ピクセラライゼーション領域として指定された各ピクセルに対して、前記ピクセラライゼーション画像におけるピクセル値をnullにセットするステップを含む、請求項１に記載の方法。
【請求項５】
前記形成されたテンプレートを前記画像訂正に使用する分離したピクセラライゼーションテンプレートに更に分割するステップを含む、請求項１に記載の方法。
【請求項６】
画像訂正の方法において、
訂正を要する画像における少なくとも１つのアーチファクトのテンプレートを用いるステップと、
前記テンプレートと前記画像とに基づきスカラー積を生成するステップと、
前記スカラー積に基づき前記テンプレートと前記画像との間の適合度を取得するステップと、
前記取得された適合度に基づき前記画像を訂正するステップとを有する、方法。
【請求項７】
前記テンプレートを用いるステップが、前記テンプレートを生成するステップ、又はデータベースから前記テンプレートにアクセスするステップを有する、請求項６に記載の方法。
【請求項８】
前記スカラー積を生成するステップの前に、前記テンプレート又は前記画像における変調効果を除去するため、前記テンプレート及び前記画像の少なくとも１つを正規化するステップを有する、請求項６に記載の方法。
【請求項９】
前記テンプレートと前記画像との間の前記適合度を決定した後、前記変調効果を回復するステップを含む、請求項８に記載の方法。
【請求項１０】
イメージング手順で用いるコンピュータ可読媒体であって、
少なくとも２つのゲインマップからピクセラライゼーション画像を生成するよう構成されるコードと、
前記ピクセラライゼーション画像における各ピクセル位置でグレー値の局所分散を計算するよう構成されるコードと、
前記ピクセラライゼーション画像における、ピクセル状のピクセラライゼーション領域及び非ピクセラライゼーション領域を決定するよう構成されるコードと、
前記決定されるピクセラライゼーション領域と非ピクセラライゼーション領域とに基づき、前記ピクセラライゼーション画像の二値画像を生成するよう構成されるコードと、
前記ピクセラライゼーション画像と前記二値画像とに基づきテンプレートを形成するよう構成されるコードとを有する、コンピュータ可読媒体。
【請求項１１】
イメージング手順で用いるコンピュータ可読媒体であって、
訂正を要する画像における少なくとも１つのアーチファクトのテンプレートを用いるよう構成されるコードと、
前記画像と前記テンプレートとに基づきスカラー積を生成するよう構成されるコードと、
前記ベクトル積に基づき前記テンプレートと前記画像との間の適合度を決定するよう構成されるコードと、
前記適合度に基づき前記画像を訂正するよう構成されるコードとを有する、コンピュータ可読媒体。
【請求項１２】
画像訂正に使用するテンプレートを生成するシステムであって、
少なくとも２つのゲインマップに基づきピクセラライゼーション画像を生成する手段と、
前記ピクセラライゼーション画像での各ピクセル位置でグレー値の局所分散を計算する手段と、
各ピクセル位置での成分をピクセラライゼーション領域又は非ピクセラライゼーション領域へと分類する手段と、
前記分類されたピクセラライゼーション領域と非ピクセラライゼーション領域とに基づき、二値画像を生成する手段と、
前記二値画像と前記ピクセラライゼーション画像とに基づきテンプレートを生成する手段とを有する、システム。
【請求項１３】
画像を訂正するシステムであって、
訂正を要する画像における少なくとも１つのアーチファクトのテンプレートを用いる手段と、
前記テンプレートと前記画像とに基づきスカラー積を生成する手段と、
前記スカラー積に基づき前記テンプレートと前記画像との間の適合度を取得する手段と、
前記適合度に基づき前記画像を訂正する手段とを有する、システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、X線イメージングの分野に関し、より詳細には、フラットパネルディテクタを用いるデジタルイメージングの分野に関する。
【背景技術】
【０００２】
イメージング技術は従来、所望の注目領域から画像データをキャプチャするため、写真用フィルムを画像データレセプタとして用いてきた。近年、画像レセプタに関連する技術における重要なシフトが行われており、それは、アナログ技術の使用からデジタル技術の使用へと進化させるものである。
【０００３】
フラットパネルディテクタ(FPD)技術が、イメージング分野において広く実現されており、アナログイメージングからデジタルイメージングへの遷移に影響を与えている。イメージングデータを素早く取得できることとは別に、FPDは、イメージインテンシファイアといった従来利用可能な技術と比べて、そのコンパクトなサイズとより長寿命である点とが広く知られている。人間医学における利用増大が見られるだけでなく、歯科、非破壊検査、及び獣医学を含む他の分野においても、FPDにより提供される利点を理解し始めている。例えば、時間浪費する処理を省略できる点や、フィルムを保管するコストを減らせる点である。
【０００４】
FPDは、１つ又は複数のフォトディテクタと共にシンチレータを用いる。放射線が、通常アクティブ物質としてヨウ化セシウム(CsI)を有するシンチレータにあたるとき、シンチレータは、入射放射線を光の光子に変換する。これらの光子がフォトディテクタにあたるとき、光子の数に比例して多数の電子が生成される。こうして、画像情報が、更なる処理のための電気信号へと変換される。
【発明の開示】
【発明が解決しようとする課題】
【０００５】
FPDは、感度変動が原因で貧弱な画像品質を生み出す。これは、様々な理由によるものである。例えば、FPDとシンチレータとの間の気泡の存在、温度変動、FPDの製造及び組み立て時における欠陥等の理由がある。FPDがアクティブクーリングを持つ場合、即ち、温度が一定に保たれるとき、局所ピクセル強度における変動は、普通に知られるゲイン訂正法を用いて訂正されることができる。しかしながら、FPDがアクティブクーリング機構を持たない場合、FPDのフォトダイオードとシンチレータとの間の光学接触が弱い及び／又は不安定であることによりもたらされる画像アーチファクトが、ディテクタの温度と共にドリフトする可能性がある。
【０００６】
ディテクタの温度がゲイン較正の間に使用される温度と異なるとき、この従来使用されるゲイン訂正法は、画像アートファクトを訂正する際実効的ではなくなる。ドリフト感度は、結果として、貧弱な画像品質を生じさせる。従って、温度依存効果を訂正するイメージング手順を持つことが有利であろう。
【課題を解決するための手段】
【０００７】
従って、画像訂正に使用する少なくとも１つのアーチファクトのテンプレートを生成するシステム及び方法が本書で開示される。そのアーチファクトを含む画像は、少なくとも２つの一様な照射(exposure)を用いて生成される。各照射は、異なるディテクタ動作温度で生成される。その画像における各ピクセル位置でのグレー値の局所分散が計算される。その後、画像における各ピクセルが分類される。その分類に基づき二値画像が生成される。その後、その二値画像とそのアーチファクトを含む画像データとの両方に基づき、テンプレートが形成される。
【０００８】
更に、画像訂正のシステム及び方法が本書に開示される。訂正が必要な画像におけるアーチファクトを示す連続ピクセルの少なくとも1つの領域のテンプレートが使用される。そのテンプレートと画像とに基づかれるスカラー積が生成される。その後、そのスカラー積に基づいて決定される適合度に基づき、画像が訂正される。
【０００９】
更に、上述した方法を実現するコードを含む、コンピュータ可読媒体が開示される。
【発明を実施するための最良の形態】
【００１０】
図面を参照して以下の詳細な説明を読めば、様々な特徴、側面及び利点が明らかとなるであろう。
【００１１】
本書における説明において、ピクセラライゼーション(pixellisation)は、FPDのフォトダイオードとシンチレータとの間の弱くかつ不安定な光学接触によりもたらされる画像アーチファクトに対して与えられる用語である。光学接触は、温度の関数として変化する。ピクセラライゼーションは、検出器感度を低下させる。アーチファクト性質上、感度の空間変動は高いものとなる。温度もピクセラライゼーション効果を持つ。なぜなら、感度の空間変動は温度の関数としてドリフトするからである。ピクセラライズドされた領域、即ちピクセラライゼーションを示す領域における密度分布は、ランダムなノイズ信号の密度分布と同様に振舞う。ピクセラライゼーションの別の特徴は、温度の関数として、ピクセラライズド領域における平均強度の変化が、非ピクセラライズド領域、即ちピクセラライゼーションを示さない領域の平均強度変化と類似し、それに匹敵することである。これは、ピクセラライゼーション領域の内側及び外側における感度の全体の変化が似ていることを意味する。しかしながら、感度ドリフトの空間分散は、ピクセラライゼーションにおいてより高い。感度ドリフトは、特定の温度での感度と、その特定の温度とは異なる温度での感度との比として規定されることができる。ピクセラライゼーションは、空間周波数スペクトルの高端部における画像を劣化させる。
【００１２】
高い空間分散、即ち、ピクセラライズド領域において感度分布が確率的な(ランダムな)様相を呈していることは、ピクセラライゼーションが、臨床画像における臨床情報から容易に区別されることを可能にする。感度分布のランダムな様相は、訂正法が、画像内容にほとんど影響を与えないであろうことを意味する。
【００１３】
ある実施形態において、画像訂正がサンプル画像に実行されることができる前に、ピクセラライゼーションテンプレートが生成される必要がある。このテンプレートは、サンプル画像におけるピクセラライゼーションを訂正する指紋として機能する。従って、テンプレートを生成するステップは、画像訂正に関する較正段階と考えられることができる。図１は、テンプレートを生成する方法を説明するフローチャートを示す。以下本書のセクションにおいて、FPDは、ディテクタと呼ばれることになる。
【００１４】
テンプレートのサイズは、そのテンプレートが、ディテクタにおける温度非一様性に対して不変であるよう選択される。これは、テンプレートのサイズが適切に選択される場合、ディテクタ方向の関数としてテンプレートにおける温度変動が低いであろうことによる。
【００１５】
第１のステップ１００は、ピクセラライゼーション画像を生成することである。このステップは、ピクセラライゼーション画像としても知られる、浮動小数点画像を得るためピクセル状に(pixelwise)分割されることができる少なくとも２つのゲインマップを用いる。X線量子ノイズの残りとは別に、ピクセラライゼーション画像は理想的にはフラット画像である。即ち、異なるピクセルにわたりグレーレベル値に変動がないことを意味する。
【００１６】
ピクセラライゼーション画像を生成するために、少なくとも２つの一様な照射又は画像取得が、少なくとも２つのディテクタ動作温度でそのディテクタに関して実行される。一様な照射は、X線量子ノイズを減らすため、系列状の一様な照射の時間平均である。これが、「ゲインマップ」と呼ばれる画像を各ゲイン較正に対して生成する。ある実施形態において、その較正は、最低及び最高のディテクタ動作温度で実行されることができる。これにより、特徴的なゲインドリフトパターンが、最大振幅を持つことになる。しかしながら、許容可能な最低及び最高ディテクタ動作温度間の任意の適切な温度範囲が、適切に選択されることができる。ゲイン較正は、ピクセラライゼーションパターンにおける平面X線量子ノイズの量を減らすため、高線量で複数の画像の平均に基づき実行されることができる。オプションで、ディテクタ感度変動の規則的な訂正を実行するための個別のゲイン較正が実行されることができる。
【００１７】
オプションで、ディテクタの欠陥マップが利用可能であるとき、ディテクタにある欠陥ピクセルにより画像統計が乱されることを避けるため、ピクセラライゼーション画像に関して欠陥排除ステップ１５０が実行されることができる。例えば、20%を超える感度変動が除去されることができる。別の実施形態では、欠陥ピクセルの検出を改善するため、追加的な事前フィルタリングが実行されることができる。更に、ピクセルにおける任意の欠陥を低減又は除去するため、ディテクタに関して欠陥訂正ステップ１６０が実行されることもできる。所望の効果を実現するため、従来より知られるディテクタ欠陥訂正法が採用されることもできる。
【００１８】
次のステップ２００は、ピクセラライゼーション画像における各ピクセル位置でグレー値の局所分散を決定することである。(実行されるとき)欠陥検出から廃棄されるピクセルは、このステップから排除されるべきである。ある実施形態では、ピクセラライゼーション画像の5x5のサブウィンドウにおける各ピクセル位置で局所分散が計算される。合理的な量の空間分解能を持つ一方、ピクセラライゼーションの信頼性の高い検出を可能にするよう、サブウィンドウサイズの選択がなされることができる。
【００１９】
次のステップ３００は、ピクセラライゼーション画像からピクセラライズド領域及び非ピクセラライゼーション領域を決定することである。このため、適切な閾値レベルが決定されなければならない。ある実施形態では、ランクオーダフィルタリングが使用されることができる。例えば、1400 x 1440ピクセルを持つディテクタのエッジから40ピクセルを超えないピクセラライゼーションを想定して考える。この場合、ピクセラライゼーションピクセルの最悪な量は、想定されるディテクタサイズに対して12%より少ないものとして計算されることができる。ピクセラライゼーションが、ディテクタのエッジの周囲にのみ存在すると想定される一方、ピクセラライゼーションはディテクタにわたりいずれの場所又はすべてに存在すると想定することも可能である。ピクセラライゼーションピクセルの量を表す割合は、ディテクタのサイズと存在するピクセル総数とに関して変化することになる。以下の議論において、明確さのため、ピクセラライゼーションは、ディテクタのエッジから40ピクセルを超えることはないものと想定して考える。
【００２０】
ピクセルの80%が低い分散レベルを持つ分散レベルを見つけることにより、ピクセラライゼーションのないピクセルを見つけることができる。その後これらのピクセルは、非ピクセラライズドピクセルの統計量を決定するのに使用されることができる。ある実施形態では、分散画像の位置座標(i,j)でのピクセルが、値Pvar(i,j)を持つ場合、標準偏差の空間平均値は、以下の式から決定されることができる。即ち、標準偏差の空間平均値(Mσ)は、Pvar(i,j) < Pvar(80%)であり、かつ
2009531109.tif 000003
を満たすすべてのピクセルに対して、
2009531109.tif 000004
である。ここで、nは、Pvar(80%)より少ない分散を持つ非欠陥ピクセルの数である。
【００２１】
局所標準偏差値の標準偏差は、同じ条件が適用されるすべてのピクセルに対して、
2009531109.tif 000005
として計算される。ここで、Mσは、標準偏差の空間平均値である。
【００２２】
ピクセラライズド領域と非ピクセラライズド領域との間を区別する適切なピクセラライゼーション閾値は、
2009531109.tif 000006
を用いて見つけられることができる。ここで、Fは、標準偏差値に対する乗数である。正規分布に対して、乗数Fの適切な値は、3から4の範囲にある。他の実施形態では、閾値は、当業者に知られる技術を用いて適切に発見されることができる。
【００２３】
一旦閾値が見つかると、二値画像が生成される(ステップ４００)。これは、各ピクセルがNULL値又は非NULL値のいずれかを持つことを意味する。ある実施形態において、非NULLピクセルは、ピクセラライズドピクセルとして表される。他の実施形態では、NULLピクセルが、ピクセラライズドピクセルとして表されることもできる。しかしながら、その式に対する適切な修正がなされる必要がある。二値画像において、すべてのピクセルは閾値に基づきマークされることに留意されたい。ピクセラライゼーションを示すピクセルは、ピクセラライゼーションを示さないピクセルとは異なってマークされる。ピクセラライゼーションの有無に基づくピクセルのマーキングの態様は、いずれかの特定の方法に限定されるものではない。任意の適切な識別マークが使用されることができる。
【００２４】
オプションで、最大分散閾値法４５０も実行されることができる。ディテクタの高品質欠陥マップが利用可能であるとき、このステップは必要とされることはない。このステップの１つの利点は、非現実的に高いピクセラライゼーション分散を持つピクセルが除外されることができる点にある。これは、追加の欠陥検出マップと考えられることができる。
【００２５】
次のステップ５００において、接続されたピクセラライゼーション領域の数を減らすために得られた二値画像が処理される。ある実施形態では、その処理は、狭い水平ピクセラライゼーション領域(y方向におけるオープニング処理の間)と狭い垂直ピクセラライゼーション領域(x方向におけるオープニング処理の間)とを除去するための形態学的処理を含むことができる。更に、形態学的処理は、x-y方向に沿ったオープニング処理(opening operation)も含むことができる。他のタイプの形態学的処理が実行されることもできる。例えば、小さなテンプレート穴を除去するため、及び／又はテンプレート等高線付近の１つ若しくは複数のピクセルを追加するため、拡張処理(dilation operation)が実行されることができる。当業者には理解されるように、二値領域を滑らかにする様々な方法が存在する。形態学的処理は、斯かる方法の１つを表すにすぎない。しかしながら、任意の他の適切な方法が採用されることができる。
【００２６】
ある実施形態では、形態学的処理が完了した後、欠陥マップがピクセラライゼーション画像と融合されることができる。すべての欠陥ピクセル及びピクセラライゼーション検出ステップから除外されるすべてのピクセル(前述されたもの)は、NULLにセットされる。これは、ピクセラライゼーション訂正処理が、斯かるピクセルを訂正処理から除外することを可能にする。斯かるピクセルを除外する1つの利点は、斯かるピクセルが訂正処理に影響を与える(dominate)ことを除外する点にある。その後、ピクセラライゼーション画像は、二値画像と組み合わされることができる。
【００２７】
別の実施形態では、すべての形態学的処理が完了したとき、二値画像がピクセラライゼーション画像と組み合わされることができる。２つの画像を結合する1つの方法は、二値画像における対応するピクセルがNULLであるような箇所で、ピクセラライゼーション画像におけるすべてのピクセルをNULLにセットすることを含む。これにより、結果として生じるピクセラライゼーション画像におけるすべての非NULLピクセルが、ピクセラライゼーションテンプレートを形成するのに必要なデータを持つ。
【００２８】
その後、結果として生じるピクセラライゼーション画像は、ピクセラライゼーションテンプレートを形成するよう更に分割される。より一般的に本書では以下テンプレートと呼ぶ。ある実施形態では、テンプレートを生成する方法の結果が、単一のテンプレートを生じさせることができる。他の実施形態では、複数のテンプレートが生成されることができる。更に別の実施形態では、生成されるテンプレートのセットが、単一のテンプレートを形成するよう結合されることもできる。
【００２９】
図２は、１つ又は複数の画像アーチファクトを持つ、例えば臨床画像といった画像の画像訂正の方法のフローチャートを示す。以下のセクションにおいて、画像訂正の方法が、臨床画像に対して説明される。しかしながら、当業者であれば、必要な任意の画像からのアーチファクトを訂正又は除去するのに、本方法を適用することができる。本方法は、臨床画像において発見される１つ又は複数の画像アーチファクトに対して１つ又は複数のテンプレートを用いる。ある実施形態では、本書で上述した方法により、１つ又は複数のテンプレートが生成されることができる。他の実施形態では、そのテンプレートは既に存在しており、その方法は、単にこれらのテンプレートにアクセスし、画像訂正のためそのテンプレートを用いることができる。
【００３０】
ここでの議論では、画像訂正の方法が、テンプレートを用いるステップ６００を説明することになる。「用いる」という用語は、画像訂正の方法の処理の間そのテンプレートを生成すること、又は以前に生成され、データベースに格納されたテンプレートにアクセスすることのいずれかを意味することができると解釈されなければならない。
【００３１】
一旦テンプレートがアクセスされると、そのテンプレートは、大きな平滑化カーネル(smoothing kernel)を備える畳み込み処理６５０の支配下に置かれることができる。それは、局所平均化処理と考えられることができる。畳み込みステップは、低域通過フィルタリングステップとして機能し、テンプレートデータに存在する場合がある任意の故意でないオフセット又は弱いグラジエントを削除するのに使用される。オフセット及びグラジエント情報のみを備える低域通過フィルタリングの結果は、画像訂正に使用されることになる所望のテンプレート(TAC)を与えるためテンプレートデータから減算される。しかしながら、そのテンプレートは、大きなカーネルを備える任意の畳み込みを行うことなく、画像訂正に対する所望のテンプレート(TAC)として直接使用されることができる。
【００３２】
一旦所望のテンプレート(TAC)が得られると、臨床画像だけでなくテンプレートに関して高域通過フィルタリングが行われる(ステップ７００)。前述したように、ピクセラライゼーションは、周波数スペクトルの高端部における主要な貢献を持つ。高域通過フィルタリングは、画像訂正法における画像成分の影響を減らす。ある実施形態では、高域通過フィルタリングは、臨床画像に関してのみ実行されることができる。他の実施形態では、高域通過フィルタリングは、臨床画像とテンプレートとの両方に行われることができる。このステップの結果は、ピクセラライゼーションの信号遷移が、臨床画像だけでなくテンプレートの処理においても同一になるであろうことを示す。
【００３３】
画像訂正法は、更に、正規化ステップ８００も含む。このステップ８００の利点は、テンプレートの位置で画像コントラストが高いとき、訂正の有効性が増加されることができる点にある。感度ドリフトは、倍数的な現象であることが知られている。例えば、画像の明るい領域において、ピクセラライゼーションの振幅は高くなるであろう。その振幅は、その明るい領域における局所画像強度で線形にスケール化する。そのテンプレートは通常何の変調も含まない。なぜなら、それが一様に照射された画像から生成されているからである。従って、臨床画像を一様なレベルのテンプレートにするために、臨床画像におけるピクセラライゼーションの変調が除去される必要がある。局所画像強度に比例する信号により臨床画像を割ることで、これが実現される。臨床画像データの低域通過フィルタリングは、信号(PLP)を提供する。テンプレート生成の間、斯かるタイプの低域通過フィルタリングが、テンプレートにおける任意の非一様性を訂正するためそのテンプレートに関して実行されることもできる。
【００３４】
ここではテンプレート及び臨床画像は、同じサイズの２つの１次元ベクトルと考えられる。テンプレートベクトルのベクトル長でそのテンプレートのすべてのピクセルを割ることが有利であろう。言い換えると、テンプレートベクトルは、単位長へ正規化される。テンプレートベクトルと臨床画像ベクトルとを用いてスカラー積が計算されると、画像データの長さとテンプレートベクトルの長さとを適合させるため、スカラー積計算の結果とテンプレートベクトルとが掛けられる。
【００３５】
テンプレートベクトル長は、
2009531109.tif 000007
であるすべてのテンプレートピクセルに対して、式
2009531109.tif 000008
により計算されることができる。スカラー積は、テンプレートベクトルを単位ベクトルに正規化する計算を含み、
2009531109.tif 000009
であるすべてのテンプレートピクセルに対して、式
2009531109.tif 000010
により計算されることができる。
【００３６】
ここで、PHP(i,j)は、高域通過フィルタリングされた臨床画像における位置座標(i,j)での対応するピクセルである。テンプレート及び画像データが無相関であるとき、スカラー積の結果はほとんどゼロであろう。なぜなら、両方のベクトルは、平均値ゼロの符号付き数字からなるからである。スカラー積計算は、テンプレートと画像との間の相関を決定する１つの可能な方法である。
【００３７】
ここで、正規化されたテンプレートベクトルの各ピクセルは、テンプレートのピクセラライゼーションレベルを臨床画像のピクセラライゼーションレベルに適合させる係数αによりスケール化されることができる。スカラー積計算ステップ９００の間、画像ベクトル長による除算は実行されない。従って、このベクトル長との積算の必要はない。係数αは、正規化されたベクトルTHPNをベクトルPHPに適合させるのにフィットされる係数である。ベクトルTHPに適用されるようフィットされる係数Fは、αを追加項|THP|で割ることにより得られることができる。この決定の利点は、係数Fが、
2009531109.tif 000011
であるすべてのテンプレートピクセルに対して、平方根を取ることなく、式
2009531109.tif 000012
により計算されることができる点にある。
【００３８】
元のテンプレートデータTから高域通過フィルタリングされたデータTHPにするのに、及び同様に元の画像データPから高域通過フィルタリングされた画像データPHPにするのに、線形フィルタリング処理だけが実行されたので、同じ係数Fが、臨床画像Pのピクセラライゼーションレベルにもっていくために、テンプレートデータTからのピクセラライゼーションデータに適用可能でなければならない(ステップ１０００)。オフセット及びグラジエントエラーから解放されている、テンプレートTACを得るためフィルタリングに使用された比較的大きなカーネルのため、係数Fも、フィルタリングされたテンプレートデータTACに安全に適用されることができる。従って、Fは、臨床画像に存在する変調を含むテンプレートに基づき、テンプレートピクセラライゼーション画像を提供する訂正を実行するために使用されることになる係数である。
【００３９】
テンプレートピクセラライゼーション画像は、それに関連付けられる変調を持たない。前述されたように、臨床画像におけるピクセラライゼーションの実際のレベルは、ピクセラライゼーションの局所強度により変調される。正規化処理の間、まず、この変調が、低域通過フィルタリング段の結果を用いてピクセル状分割により除去される。テンプレートピクセラライゼーション画像に関する変調を回復するため、テンプレートピクセラライゼーション画像のピクセル状の乗法が実行される。積算されるピクセル状係数は、以前に得られる低域通過情報PLPである。このステップ１１００は、変調を回復するより前に実行される正規化処理とは反対のものとして考えられることができる。
【００４０】
次のステップ１２００は、臨床画像の画像訂正である。このステップは、臨床画像からピクセラライゼーションを除去する。なぜなら、臨床画像は、ピクセラライゼーションを伴う画像成分を含む一方、テンプレートピクセラライゼーション画像は、ピクセラライゼーションのみを含み、画像成分を含まないからである。ある実施形態では、画像訂正が、ピクセル状減算によりなされることができる。テンプレートデータと画像データとにおけるピクセラライゼーション間の相関が高いとき、ピクセラライゼーションレベルは、概略的な図の最終画像におけるものとして可視閾値未満のレベルにまで減少することになるであろう。そこでは、訂正が実際に行われる。
【００４１】
ある実施形態では、前述の方法は、X線イメージングシステムといったシステムで実現されることができる。そのシステムは、本方法の各機能を実現する分離したモジュールを持つことにより、その機能を実現する手段を含むことができる。他の実施形態では、様々な機能が、１つ又は少数のモジュールで実現されることができる。そのシステムは、オペレータが、システムにコマンド又は命令を提供することを可能にするオペレータワークステーションも含むことができる。そのコマンド又は命令は、臨床画像の所望の訂正が実現されるとき、訂正処理を開始及び終了するためのものである。そのシステムは、マイクロプロセッサ及びディスプレイデバイスも含むことができる。別の実施形態では、前述の方法が、スタンドアロンシステムで実現されることができる。斯かるスタンドアロンシステムは、イメージングシステム又は取得画像を含むデータベースに接続されることができる。
【００４２】
別の実施形態では、画像訂正を実行するシステムが、上述の画像訂正を実行する手段を含む。ある実施形態では、そのシステムはデータベースから１つ又は複数のテンプレートにアクセスすることができる。他の実施形態では、そのシステムは、画像訂正に使用する１つ又は複数のテンプレートを生成する手段を含むことができる。別の可能な実施形態では、その１つ又は複数のテンプレートを生成するシステムが、画像訂正のためのシステムの一部を形成することができる。その逆も成り立つ。
【００４３】
当業者であれば理解されるであろうが、画像訂正に使用する１つ又は複数のテンプレートを生成に関し実現される方法及び画像訂正に関し実現される方法は、例えばコンピュータコードの形式のプログラム命令を用いて実現されることができる。斯かるコードは、コンピュータ可読媒体に含まれることができる。１つの可能な実施形態では、そのコードは上述された方法を実現するシステムに直接格納されることができる。別の実施形態では、そのコードが具体的な媒体に含まれ、システムに供給されることができる。その媒体は、そのコードが適切に格納されることができる光学又は磁気媒体を含むことができる。斯かるコンピュータ媒体の例は、CDROM、DVD、フラッシュメモリカード、コンピュータハードディスク、フロッピー（登録商標）ディスク等を含む。
【図面の簡単な説明】
【００４４】
【図１】画像訂正において使用する少なくとも１つのアーチファクトのテンプレートを生成する方法の実現を表すフローチャートである。
【図２】画像訂正の方法の実現を表すフローチャートである。
【図１】
2009531109.tif 000013
【図２】
2009531109.tif 000014
【国際調査報告】
2009531109.tif 000015
2009531109.tif 000016