(57)【要約】画像処理システムＩＰＳは、イメージング機器ＸＡによって異なる方向に沿って取得された可動物体ＰＡＴの複数の投影画像からの投影画像を受け取る入力インターフェースＩＮを含む。複数の投影画像は、検出器Ｄの放射線感受性面が広がる投影領域において規定される。システムは、取得中の物体の動きによって投影画像が破壊されているか否かを判断するために、投影領域内で作動する動きチェッカＭＣを含む。
2022502108.tif 000002
【特許請求の範囲】
【請求項１】
イメージング機器によって異なる方向に沿って取得された可動物体の複数の投影画像であって、検出器の放射線感受性面に広がる投影領域内に規定された前記複数の投影画像からの投影画像を受け取るための入力インターフェースと、
前記投影画像が、取得中の前記可動物体の動きによって破損されているかどうかを判断するために前記投影領域内で動作する動きチェッカと、
前記動きチェッカによる判断に基づいて、前記可動物体の断面情報を含む画像を再構成する再構成器と、
を含む、画像処理システム。
【請求項２】
前記再構成器は、前記複数の投影画像から前記可動物体の前記断面情報を含む前記画像を再構成し、前記投影画像は、前記動きチェッカが動き破損があると判断した場合かつその場合にのみ、前記再構成器によって無視される、請求項１に記載の画像処理システム。
【請求項３】
前記動きチェッカが動き破損があると判断した場合かつその場合にのみ前記動きの推定値を計算する動き推定器を含む、請求項１に記載の画像処理システム。
【請求項４】
前記動き推定値に基づいて前記動きを補償して、動き補償された投影画像を得る動き補償器を含む、請求項３に記載の画像処理システム。
【請求項５】
前記再構成器は、少なくとも前記動き補償された投影画像から前記画像を再構成する、請求項４に記載の画像処理システム。
【請求項６】
前記動きチェッカは、少なくとも２つの投影画像に基づいて、前記イメージング機器のイメージングジオメトリに関連する整合性スコアを計算し、前記動きチェッカによる前記判断は、前記整合性スコアに基づいている、請求項１から５のいずれか一項に記載の画像処理システム。
【請求項７】
前記整合性スコア値を表示デバイス上で視覚化するビジュアライザを含む、請求項１から６のいずれか一項に記載の画像処理システム。
【請求項８】
前記可動物体は、異なる物理的状態となることができ、前記物理的状態間で移行することができる、請求項１から７のいずれか一項に記載の画像処理システム。
【請求項９】
前記動きチェッカは、前記投影画像が前記移行中に取得された場合に、動き破損があると判断する、請求項８に記載の画像処理システム。
【請求項１０】
前記動き補償器は、前記補償を行うために補正動作を適用し、前記動き補償器は、前記可動物体が前記物理的状態のいずれか１つにある間に取得された前記投影画像又は他の１つ以上の投影画像に同じ補正動作を適用する、請求項８に記載の画像処理システム。
【請求項１１】
前記イメージング機器は、コーンビームコンピュータ断層撮影方式機器である、請求項１から１０のいずれか一項に記載の画像処理システム。
【請求項１２】
イメージング機器によって異なる方向に沿って取得された可動物体の複数の投影画像であって、検出器の放射線感受性面に広がる投影領域内に規定された前記複数の投影画像からの投影画像を受け取るステップと、
前記投影画像が、取得中の前記可動物体の動きによって破損されているかどうかを、前記投影領域内の動作に基づいて、判断するステップと、
前記判断するステップによる判断に基づいて、前記可動物体の断面情報を含む画像を再構成するステップと、
を含む、画像処理方法。
【請求項１３】
イメージング機器と、請求項１から１１のいずれか一項に記載の画像処理システムと、を含む、画像処理装置。
【請求項１４】
少なくとも１つの処理ユニットによって実行されると、前記少なくとも１つの処理ユニットに請求項１２に記載の方法を行わせる、コンピュータプログラム。
【請求項１５】
請求項１４に記載のコンピュータプログラムを記憶した、コンピュータ可読媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理システム、画像処理方法、イメージング装置、コンピュータプログラム要素、及びコンピュータ可読媒体に関する。
【背景技術】
【０００２】
今日では、コーンビームＣＴ（コンピュータ断層撮影）スキャン中に患者の動きを検出する方法は、主に外部センサに基づいている。光学カメラベースのシステムで使用されている技法には、患者の皮膚に取り付けられる受動マーカ及び能動マーカが含まれる（ただし、これらに限定されない）。
【０００３】
患者の動きは、特にコーンビームＣＴのイメージングを妨害する。これは、スキャン速度がファンビームＣＴシステムと比較して依然として遅いことによる。測定データに対する最初の再構成の剛体３Ｄ／２Ｄ位置合わせは、すべての投影画像に対して剛体運動を推定するためのデータ駆動型方法である。このデータ駆動型アプローチは、３Ｄ順投影を必要とし、この動作は計算コストが高いことから、かなりの時間を要する。例えば、Ａ．Ｓｉｓｎｉｅｇａ他の「Ｍｏｔｉｏｎ  ｃｏｍｐｅｎｓａｔｉｏｎ  ｉｎ  ｅｘｔｒｅｍｉｔｙ  ｃｏｎｅ−ｂｅａｍ  ＣＴ  ｕｓｉｎｇ  ａ  ｐｅｎａｌｉｚｅｄ  ｉｍａｇｅ  ｓｈａｒｐｎｅｓｓ  ｃｒｉｔｅｒｉｏｎ」（Ｐｈｙｓ．Ｍｅｄ．Ｂｉｏｌ．、第６２巻、２０１７年、３７１２頁）を参照されたい。
【発明の概要】
【発明が解決しようとする課題】
【０００４】
したがって、より反応が早いＣＴイメージングが求められている。
【課題を解決するための手段】
【０００５】
本発明の目的は、独立請求項の主題によって解決され、さらなる実施形態は従属請求項に組み込まれる。なお、以下に説明する本発明の態様は、画像処理方法、イメージング装置、コンピュータプログラム要素、及びコンピュータ可読媒体に等しく適用されることに留意されたい。
【０００６】
本発明の第１の態様によれば、画像処理システムが提供される。画像処理システムは、
イメージング機器によって異なる方向に沿って取得された可動物体の複数の投影画像からの投影画像を受け取るための入力インターフェースと、
取得中の物体の動きによって投影画像が破損されているかどうかを判断する動きチェッカと、
動きチェッカによる判断に基づいて、物体の断面情報を含む画像を再構成する再構成器とを含む。
【０００７】
実施形態では、動きチェッカは、投影画像が規定される投影領域において、好ましくは排他的に動作する。投影領域には、検出器の放射線感受性面が広がる。投影領域は好ましくは２Ｄである。これにより、３Ｄでの複雑な計算が不要であるため、判断を迅速に行える。動きチェッカの動作は、尺度又はスコア、すなわち、数を計算することを含み、この尺度は、決断のために１つ以上の閾値と比較され得る。システムが再構成器を含むことが好ましいが、これは必須ではなく、また、教示目的のため、又は患者の動き挙動を分析するときなど、必ずしも再構成がなくても、動きチェッカの出力が考慮される代替案が想定される。
【０００８】
実施形態では、尺度は、データ整合性を評価する。実施形態では、基礎となるイメージングジオメトリの仮定に由来するデータ冗長性が利用される。データ冗長性に対する妨害に敏感な数値尺度が構築され、そのような妨害は患者の動きによって引き起こされたものとみなされる。
【０００９】
一実施形態では、再構成器は、複数の投影画像から物体の断面情報を含む画像（３Ｄボリュームなど）を再構成するが、動きチェッカが動き破損があると判断した場合にのみ、投影画像は再構成器によって無視される。
【００１０】
一実施形態では、システムは、動きチェッカが動き破損があると判断した場合にのみ、動きの推定値を計算する動き推定器を含む。
【００１１】
一実施形態では、システムは、動き推定値に基づいて動きを補償して、動き補償された投影画像を得る動き補償器を含む。動き補償された画像を形成することは、（ピクセルで符号化された）画像情報を変更することを含み得るが、代わりに、投影画像に割り当てられたイメージングジオメトリのみを変更することを含んでもよい。他の実施形態では、イメージングジオメトリ及びピクセルの両方が変更されてもよい。
【００１２】
好ましい実施形態では、動きチェッカの計算量は、動き推定器及び動き補償器の計算量よりも少ない。これにより、動き推定及び補償は動き破損された投影画像に対してのみ行われるので、計算時間をさらに短縮できる。好ましくは、システムはモジュラー式であり、動きチェッカは、動き推定器及び／又は動き補償器とは異なる処理段として配置される。好ましくは、動き推定及び動き補償も別個の段として配置され、これらを使用できる。
【００１３】
これらの実施形態では、再構成器は、少なくとも動き補償された投影画像から画像を再構成する。画像は、３Ｄボリュームであってよい。
【００１４】
一実施形態では、動きチェッカは、少なくとも２つの投影画像に基づいて、イメージング機器のイメージングジオメトリに関連する整合性スコアを計算し、動きチェッカによる判断は、上記スコアに基づいている。
【００１５】
一実施形態では、システムは、スコア値を表示デバイス上で視覚化するビジュアライザを含む。
【００１６】
一実施形態では、物体は、異なる物理的状態となることができ、また、物理的状態間で移行することができる、
【００１７】
一実施形態では、動きチェッカは、投影画像が移行中に取得された場合に、動き破損があると判断する。
【００１８】
一実施形態では、動き補償器は、補償を行うために補正動作を適用し、動き補償器は、物体が物理的状態のいずれか１つにある間に取得された投影画像又は他の１つ以上の投影画像に同じ補正動作を適用する。これにより、動き推定を投影画像に対して個別に新たに行う必要がなく、再利用できるため、さらに一層計算時間を削減できる。
【００１９】
さらなる態様によれば、画像処理方法が提供される。方法は、
イメージング機器によって異なる方向に沿って取得された可動物体の複数の投影画像からの投影画像を受け取るステップと、
投影画像が、取得中の物体の動きによって破損されているかどうかを判断するステップと、
上記判断ステップによる判断に基づいて、物体の断面情報を含む画像を再構成するステップとを含む。
【００２０】
投影画像は、検出器の放射線感受性面が広がる投影領域内で規定され、投影画像が動き破損されているかどうかを判断するステップは、実施形態では、投影領域内で行われる１つ以上の動作に基づいている。ここでも、システムについて上述したように、実施形態では、再構成ステップを省略してもよい。
【００２１】
さらなる態様によれば、イメージング機器と、上記の実施形態のいずれか１つによるシステムとを含む画像処理装置が提供される。
【００２２】
さらなる態様によれば、少なくとも１つの処理ユニットによって実行されると、処理ユニットに方法を行わせるコンピュータプログラム要素が提供される。
【００２３】
さらなる態様によれば、プログラム要素を格納したコンピュータ可読媒体が提供される。
【００２４】
動き補正を動き破損されている投影に限定したことにより、動き補正の全体的なロバスト性が向上される。さらに、動き補正によるアイソセンタ又は向きの変化などの影響が回避される。
【００２５】
実施形態では、動きによって実際に影響を受けているビューに対してのみ、順投影及び３Ｄ／２Ｄ位置合わせを行うことによって、第２のパスの動き補償された再構成を行う効率的なやり方が提案される。
【図面の簡単な説明】
【００２６】
本発明の例示的な実施形態について、次の図面を参照しながら説明する。
【００２７】
【図１】図１は、イメージング装置の概略ブロック図を示す。
【図２】図２は、コーンビーム用の例示的なイメージングジオメトリを示す。
【図３】図３は、画像処理システムの概略ブロック図を示す。
【図４】図４は、動き破損を確定するためのスコア値の曲線を示す。
【図５】図５は、画像処理方法のフロー図を示す。
【発明を実施するための形態】
【００２８】
図１を参照すると、イメージング装置ＩＡの概略ブロック図が示されている。
【００２９】
イメージング装置１Ａは、特に、Ｘ線イメージング機器ＸＡ（「イメージャ」）、好ましくはコーンビームタイプのＣＴ(コンピュータ断層撮影）スキャナと、画像処理システムＩＰＳとを含む。イメージング機器は、Ｃアーム又はＵアームスキャナとして可動又は固定ユニットとして配置されるか、又はドーナツ型ガントリ内に配置され得る。
【００３０】
画像処理システムＩＰＳは、画像取得中にＸ線イメージャＸＡによって取得された投影画像πを処理する。画像処理システムＩＰＳは、汎用コンピュータなどの計算ユニットＰＵ上のハードウェア回路及び／又はソフトウェア構成要素として構成され得る。画像処理システムは、通常、１つ以上のＸ線イメージャＸＡに関連付けられているワークステーション上に実装され得る。実施形態では、画像処理システムＩＰＳが１つ以上のイメージャＸＡのために機能するクラウドベースのアーキテクチャも想定される。画像プロセッサは、集積回路、ＦＰＧＡ、ＡＳＩＣなどの専用ハードウェア構成要素として少なくとも部分的に実装され得る。ＩＰＳは、Ｘ線イメージャＸＡに（完全に）統合されていてもよい。Ｘ線イメージャＸＡから画像プロセッサＩＰＳへの投影画像πの通信は、有線又は無線通信ネットワークを介して行われ得る。
【００３１】
イメージャＸＡを使用するとき、患者ＰＡＴの関心領域（ＲＯＩ）は、イメージャＸＡのＸ線源ＸＳとＸ線感受性検出器Ｄとの間の検査領域に配置される。画像取得の間、Ｘ線源ＸＳは、関心領域の周りの軌道又は軌跡Ｔを辿るように移動する。実施形態では、検出器Ｄも、Ｘ線源ＸＳと協働して上記軌跡を辿り、関心領域の周りのＸ線源の移動中、Ｘ線源ＸＳの反対側にあり続ける。好ましくは、必須ではないが、軌跡は、平面円又は楕円の少なくとも一部である。軌跡は、必ずしも閉じた形である必要はない。
【００３２】
あるいは、検出器Ｄが、患者ＰＡＴがその中に位置する管又はリングとして本質的に配置され、Ｘ線源だけが上記管又はリング内の軌跡Ｔを追跡する、より高い世代のイメージャＸＡもまた、本明細書で想定される。軌跡Ｔに沿ったＸ線源ＸＳの移動中、少なくとも関心領域を照射するようにＸ線ビームがＸ線源ＸＳによって生成される。関心領域は、頭部、心臓、胸部、脚部、膝などの患者ＰＡＴの任意の解剖学的部分であってよい。
【００３３】
本明細書で好ましく想定されるコーンビーム装置のため、ビームＸＢの空間形状は、図１の破線の楕円によって示されるような空間コーンの空間形状である。コーンは、円形又は楕円形の断面を有する。コーンビームタイプは、ビームが発散平面扇形の形を有するか、又は平行ビームジオメトリが使用される他のタイプのＣＴイメージング機器とは対照的である。これらの他のタイプのいずれにおいても、ビームは実質的に２Ｄ構造であるのに対して、コーンビームは３Ｄ構造である。
【００３４】
Ｘ線ビームＸＢは、患者ＰＡＴの組織と相互作用する。その結果、ビームＸＢは吸収、位相シフト／偏差、及び種々の散乱のために、特に減衰を経験する。この相互作用のために、Ｘ線ビームＸＢは変更され、（Ｘ線源の位置に対して）遠い側で患者から出て、検出器Ｄの検出器感知表面に衝突する。検出器表面は２Ｄ構造である。
【００３５】
行及び列に配列された検出器ピクセルは、検出器表面をレイアウトするか又は他の方法で形成し、変更された放射線を検出する。検出器ピクセルによる検出動作は、変更された放射線を電気信号に変換することを含む。電気信号は、データ取得ユニットＤＡＱによって、デジタル投影画像のインスタンスを形成する一組の数（「投影画像値」又は「ｐ値」）に変換される。線源ＸＳが、患者ＰＡＴに対して異なる空間位置をとるように軌跡Ｔに沿って適切な増分で移動するにつれて、そのような位置ごとに、異なる投影画像（又はフレーム）が検出器によって検出される。このようにして、（投影）画像取得の終了時には、投影画像の集合（軌跡Ｔ上の各位置に対して１つ以上）が得られる。
【００３６】
本明細書で主に（ただし、排他的にではなく）想定されるコーンビームのジオメトリのために、各投影画像は、好ましくは行列（ｕ，ｖ）に配列され格納されるｐ値の２次元アレイである。
【００３７】
Ｘ線源ＸＳが患者の周りで完全な曲線（例えば、円）を辿ることが可能であるが、必ずしも必要ではない。通常、１８０°以下の円弧で十分である。
【００３８】
投影画像の集合は、画像処理システムＩＰＳに転送され得る。投影画像は、画像取得中に、それらが取得されるときにそれぞれ転送されても、一括して、一度に異なるセットで転送されてもよい。あるいは、すべての投影画像を、画像取得の終了時にまとめて転送してもよい。
【００３９】
画像処理システムは特に、再構成器ＲＥＣＯＮを含む。再構成器ＲＥＣＯＮは、コーンビームタイプのイメージングに適切に調整された再構成アルゴリズムを実行するソフトウェア構成要素又はハードウェア構成要素である。大まかに言えば、再構成器は、投影画像を処理して、患者ＰＡＴの体内の内部構造（器官、組織界面、病変など）の点を特に表す、ボクセルＶｘと呼ばれる画像値から構成される３Ｄ画像ボリュームを生成する。
【００４０】
次いで、このボリュームＶを通る特定のビュー、レンダリング、又は断面が、適切なリフォーマッタ又はレンダラによって生成されて、ビジュアライザＶＩＺに渡される。ビジュアライザＶＩＺは、画像値をグレー値又はカラー値にマッピングする。ビジュアライザＶＩＺは、グラフィックス回路と対話して、表示デバイスＤＤを制御する。次いで、３Ｄボリュームの要求された断面画像又はレンダリングなどの視覚的表現が表示デバイスＤＤ上に表示され、診断、ナビゲーション、介入計画又は他のタスクにおいてユーザ（例えば、放射線科医）を支援することができる。
【００４１】
ここで、再構成器ＲＥＣＯＮをより詳細に参照すると、この構成要素は、投影画像から３Ｄ画像ボリュームへのマッピングを実施する。上から分かるように、投影画像は、検出器Ｄの放射線感受性面を形成するピクセルによって規定される。この面は通常平面であるが、曲面検出器も本明細書において想定されていることから、すべての実施形態において必ずしもそうとは限らない。投影画像が規定されるこの検出器面を、投影領域と呼ぶ。これと区別するために、また、これに加えて、Ｘ線源及び好ましくは検出器が軌跡を追跡するときにその周りを回転する検査領域の３Ｄ空間内のグリッド点によって形成される第２の領域もある。再構成画像、すなわち、ＲＯＩの断面画像を含む３Ｄボリュームは、再構成の結果として再構成器ＲＥＣＯＮによってグリッド点に割り当てられるボクセル値によって規定される。したがって、再構成演算は、投影画像πを入力とし、画像値を計算し、３Ｄ内のボクセル位置に割り当てる。したがって、ボクセル位置における画像値は合わせて計算されて、再構成３ＤボリュームＶが形成される。コーンビームは、ファンビーム又は平行ジオメトリイメージングとは異なる。従来のビームジオメトリでは、Ｘ線源は、複数の断面画像を取得するために患者の長手方向軸に沿って前進する必要がある。これらの断面画像は画像ボリュームとなるようにスタックされる。コーンビームＣＴでは、イメージングが真に又は本質的に３Ｄであり、個々の断面画像を再構成し、これらをスタックする必要がない。コーンビームＣＴでは、完全な３Ｄボリュームが再構成されると、リフォーマットなどの比較的簡単な演算によって断面画像を計算することができる。
【００４２】
ファンビーム又は平行ジオメトリスキャナで取得した投影画像では、従来のフィルタ付き逆投影（ＦＢＰ）アルゴリズムを使用して再構成を行うことができる。コーンビームジオメトリでは、発散ビームジオメトリのために、これは通常当てはまらない。Ｌ．Ａ．Ｆｅｌｄｋａｍｐ他の「Ｐｒａｃｔｉｃａｌ  ｃｏｎｅ−ｂｅａｍ  ａｌｇｏｒｉｔｈｍ」（Ｊ．Ｏｐｔ．Ｓｏｃ．Ａｍ．Ａ／第１巻、第６号、１９８４年６月）の専用コーンビーム再構成アルゴリズムなどの専用のコーンビーム再構成アルゴリズムが過去に開発されている。ファンビームのジオメトリも発散型であるが、この発散は平面に限定され、一方で、複数の発散平面が存在するコーンビームには当てはまらないことが理解されるであろう。
【００４３】
ファンビームや平行ジオメトリビームのより単純なビーム構造は、それぞれの再構成アルゴリズムを導出する際に、特定の単純化仮定を行うことを可能にする。ＦＢＰなどのこれらのアルゴリズムは、フーリエ変換及びラドン変換の理論的操作を通して到達される。これらの単純化仮定は、コーンビームの場合でもはや行うことができない。しかしながら、コーンビーム状況をファンビーム又は平行ジオメトリの場合のそれに少なくとも近似的に変換するように、特定の変換Ｇ（以下により詳細に説明する）を投影データに適用できることが観察され、報告されている。Ｂｒｕｃｅ  Ｓｍｉｔｈがこのことについて書いている。例えば「Ｃｏｎｅ−ｂｅａｍ  ｔｏｍｏｇｒａｐｈｙ：  ｒｅｃｅｎｔ  ａｄｖａｎｃｅｓ  ａｎｄ  ａ  ｔｕｔｏｒｉａｌ  ｒｅｖｉｅｗ」（ＯＰＴＩＣＡＬ  ＥＮＧＩＮＥＥＲＩＮＧ、第２９巻（５）、５２４〜５３４頁）を参照されたい。
【００４４】
再構成画像Ｖ（π）の精度の大部分は、基礎となるイメージングジオメトリに関連する特定のモデル仮定に基づいている。図２に、これらの仮定をよりうまく説明するために、１つの基礎となるイメージングジオメトリを概略的に示す。図２は、具体的にはＸ線源の平面軌跡Ｔが位置する平面Π（Ｔ）を示す。実施形態では、この軌跡Ｔは、円、楕円又は他の閉曲線の一部である。非平面のＴもまた、実施形態において想定され得る。Ｔの曲線自体は閉じている場合もあるが、軌跡Ｔ自体は通常そうではなく、代わりに円弧などの曲線の一部のみを形成する。
【００４５】
軌跡の平面Π（Ｔ）に対して、検出器Ｄの平面Π（Ｄ）が垂直である。２つの平面の交点は、３Ｄ座標系（Ｕ，Ｖ，Ｗ）の１次元を規定する。別の次元Ｗは、平面Π（Ｄ）に対して垂直であり得る。したがって、検出器平面内の点は、局所座標（ｕ，ｖ）によって記述され得る。座標系の原点「０」は、Ｕ及びＶと、軌跡の平面Π（Ｔ）に垂直でかつその内に延在し、また、Ｘ線源ＸＳの位置に向かって延在する第３の次元ｗとの交点で規定され得る。
【００４６】
座標系全体は、周囲の破線の円上の矢印によって示されるように、Ｘ線源ＸＳと共に移動（例えば、回転）される。具体的には、図２は、軌跡Ｔ上のＸ線源ＸＳの所与の位置についての構成を示す。平面Π（Ｄ）は関心領域を通過し、それと交差する中心平面を表す。ビームのコーンビーム構造により、複数のそのような平面が存在することが理解されるであろう。これらの平面は、中心平面Π（Ｄ）に平行であり、中心平面に対して−ｖ及び＋ｖ方向にずれている。
【００４７】
例示的なボクセル点Ｖｘが、３Ｄ空間内のグリッド点上に位置するものとして示されている。より良く例示するために、この例示的な場合では、ボクセルＶｘは平面Π（Ｔ）外にある。幾何学的光線が、光線に沿ったすべてのボクセルＶｘ（１つのみを示す）の像点Ｐ（Ｖｘ）を規定し、Ｘ線源ＸＳの特定の位置αλについて検出器Ｄによって記録される。線源ＸＳの位置は、位置パラメータλによってパラメータ化されたベクトルαによってパラメータ化され得る。パラメータは角度であってもよい。検出器表面Π（Ｄ）上の像点Ｐ（Ｖｘ）の位置（ｕ，ｖ）は、上記幾何学的光線と検出器表面Π（Ｄ）との交点である。像点Ｐ（Ｖｘ）の画像値は、上記幾何学的光線上にあるすべてのボクセル（Ｖｘを含む）の線積分である。幾何学的光線は、Ｘ線ビームの一部が進む経路をマークする。すべてのボクセルＶｘに対するすべての幾何学的光線（１つのみを示す）は合わされて空間コーンを形成する。図２の構成は、イメージングジオメトリを規定する。図２に示すイメージングジオメトリは、検出器とＸ線源との間の空間的関係、及びイメージング物体（この例では患者の解剖学的構造又は関心領域）に対するそれらの空間的関係を規定する。
【００４８】
基礎となるジオメトリに基づくこのモデル化の仮定は、例えば、各投影画像πが軌跡Ｔ上の事前に規定された線源ＸＳ位置から取得されることを必要とすることがある。したがって、投影画像πと線源ＸＳの位置パラメータλとの間に関連付けが存在する。さらに、軌跡Ｔは、画像取得全体を通して不変と仮定される事前に規定された形状を有すると仮定される。図２の例示的な実施形態の場合、Ｔの形状は、円形であるか、少なくとも円弧であると仮定される。事前に規定された軌跡を仮定することは、線源と患者ＰＡＴとの間、又は少なくとも線源ＸＳと関心領域ＲＯＩとの間の不変の空間的関係を暗示する。
【００４９】
残念なことに、これらのモデル仮定は、患者の動きを含む外部因子により妨害されることがある。例えば、脳スキャンの間に、患者は、咳払いしなければならないことがある。これは、Ｘ線源ＸＳが事前に規定された軌跡Ｔ上でその経路を辿る間に、画像取得中に頭部の動きを引き起こす。動きが生じたとき、基礎となるジオメトリの仮定はもはや適用されないか、又は取得全体を通して有効でなくなり、考慮されなければ、イメージングアーチファクトがもたらされる。例えば、関心領域の動きがある場合、相対軌跡Ｔはもはや円形ではなく、より複雑な仮想形状を有するボクセル位置がある。したがって、基礎となるモデル化の仮定は壊れ、再構成アルゴリズムは、少なくとも局所的に、画質（「ＩＱ」）を損なうアーチファクトをもたらす可能性がある。
【００５０】
この状況を軽減するために、提案される画像処理システムＩＰＳは、画像取得中に生じる可能性がある動き破損を補償又は相殺することを可能にする動き合成機能を含む。
【００５１】
しかしながら、現存の動き補償方式とは異なり、提案される画像処理システムは選択的である。つまり、出願人によって発見されたように、画像取得中に、すべての投影画像が動きによって破損されるわけではない。したがって、すべての投影画像が動き補償されなければならないわけではない。例えば、上記の脳イメージングの例では、外乱は、画像取得の初めから終わりまでショートバーストで発生する。例えば、咳発作によって引き起こされる頭部の動きは、時間的に分離した「特異点」であり、非常に短い期間に急激な動きのピークを引き起こす。したがって、本出願人は、動き破損によって実際に影響を受けた投影画像に対してのみ動き補償を行う、的を絞った又は選択的な動き補償を提案する。したがって、動きによって影響を受けていない残りの投影画像は保持されるか、又は、ある場合に、少なくとも、動き補償が、すべての投影画像の集合のうちのいくつかのフレームについて単純化されてもよい。この選択的アプローチにより、計算時間が削減され、再構成の反応をより早めることができる。
【００５２】
画像処理システムＩＰＳの動作をより詳細に説明する前に、以下の観察が役立つであろう。図２のイメージングジオメトリの上記の説明から、投影画像は、以下の情報成分を含むことが明らかになるであろう。すなわち、投影画像は、検出器ピクセルによって記録された行及び列（ｕ，ｖ）にされた画像値Ｐ（Ｖｘ）を含むが、基礎となるイメージングジオメトリにおける各像点Ｐ（Ｖｘ）の予想される空間位置λ、特に、線源ＸＳに対するその関係、及び／又は線源からＲＯＩを通って像点Ｐ（Ｖｘ）までの光線に沿った関連するボクセル位置Ｖｘも含む。さらなるジオメトリ情報は、検出器の位置、検出器の向き、患者テーブルの位置、又はアイソセンタ位置のうちの任意の１つを含み得るが、これらに限定されない。したがって、投影画像は、π＝（ｕ，ｖ，λ）と書くことができ、（ｕ，ｖ）は画像ピクセルであり、λはイメージングジオメトリを記述する３Ｄ空間座標である。イメージングジオメトリは、特に、取得時の３Ｄ空間における線源ＸＳの位置について記述する。したがって、ある線源位置αλで取得された所与の瞬間における任意の１つの投影画像は、画像情報（ｕ，ｖ）と、本明細書でλと呼ぶそのイメージングジオメトリの空間座標とを含む多成分構造である。「λ」は、一般的な意味で使用され、必ずしも図２の例示的な円形ジオメトリに関係するとは限らない。空間情報、特に空間座標λは、例えば、対応する投影画像などのヘッダファイル内に符号化されることによって、各投影画像に適切に関連付けられる。あるいは、空間情報は、較正データファイルなどのファイルとしてメモリに別個に格納されてもよい。
【００５３】
座標λは、空間における対応する投影画像とＲＯＩとの予想される空間構成を指定する。特に、座標は、ＲＯＩとＸ線源ＸＳ及び／又は検出器Ｄとの予想される空間コンステレーション（配置）を指定する。
【００５４】
ここで図３を参照すると、本明細書で想定される動き補償機能を有する画像処理システムＩＰＳの概略ブロック図が示されている。
【００５５】
画像処理システムは、入力ポートＩＮを含み、ここで、好ましくはＡＤ変換後に、動き破損された可能性のある投影画像πがイメージャＸＡから受信される。πはまた、ＰＡＣＳ又は他のメモリなどの記憶装置からイメージング後に取得されてもよい。
【００５６】
各投影画像は、数の行列として表現可能である。各投影画像は、上述したように、基礎となるイメージングジオメトリに対応する特定の空間座標に関連付けられている。
【００５７】
投影画像のいくつか、又は好ましくはすべてが、動きチェッカＭＣＨによって処理される。動きチェッカＭＣＨは、投影画像の所与のフレームが動きによって破損されているかどうかを確定するためにスコアを計算することを可能にする動き検出尺度ｍｅａｓ（π）を使用する。一実施形態では、コーンビームデータ整合性尺度ｍｅａｓ（・）＝ＤＣＣ（・）が、例えば「４ｔｈ  Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｃｏｎｆｅｒｅｎｃｅ  ｏｎ  Ｉｍａｇｅ  Ｆｏｒｍａｔｉｏｎ  ｉｎ  Ｘ−Ｒａｙ  Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ」（２０１６年７月１８〜２２日、ドイツ、バンベルグ）で発表されたＪｅｒｏｍｅ  Ｌｅｓａｉｎｔ他による「Ｔｗｏ  ｃｏｎｅ−ｂｅａｍ  ｃｏｎｓｉｓｔｅｎｃｙ  ｃｏｎｄｉｔｉｏｎｓ  ｆｏｒ  ａ  ｃｉｒｃｕｌａｒ  ｔｒａｊｅｃｔｏｒｙ」において報告されているように、かつ、以下で図５でより詳細に説明されるように、使用される。実施形態では、ＤＤＣは、投影画像を構成するいくつかの又はすべての他のフレームとのフレームの整合性を評価する。実施形態では、基礎となるイメージングジオメトリ及び／又は静止している被験者の仮定との画像πの整合性が定量化される。
【００５８】
所与のπについての動き破損尺度ｍｅａｓ（π）は、所与の投影画像πの前及び／又は後に取得された追加の投影画像の関数であってよい。好ましくは、少なくとも２つの投影画像が処理されて、所与の投影画像πに対するスコアｍｅａｓ（π）が生成される。例えば、一実施形態では、所与の投影画像πλ１の場合、所与の画像πλ１のスコアを計算するために、後の時点λ１＜λ２で取得された後続の画像πλ２も使用される。３つ以上の後続の投影画像及び／又は１つ以上の前のフレームも、所与の時点πλ１についての計算に含められ得る。本明細書で使用される場合、「前の」又は「後の」フレームは、投影画像が取得されたときにＸ線源ＸＳがある位置を規定する位置パラメータλに対応し得る。
【００５９】
投影画像又はそのセットは、物体が動いている間に取得された場合、動き破損されているとみなされる。さらに、少なくとも２つの投影画像からなるセットは、それらが物体が２つの異なる物理的状態にあるときに取得された場合、動き破損されているか、又は不整合であるとみなされ得る。したがって、１つの状態で取得された１つのフレームは、物体の異なる状態で取得されたフレームと比較して、動き破損されてる可能性がある。本明細書で使用される「状態」の概念は、後に図４Ｂにおいて詳述される。
【００６０】
まず、画像処理システムＩＰＳをより詳細に参照すると、選択的な動き補償の様々な実施形態が想定される。
【００６１】
これらの実施形態のうちの１つでは、画像処理システムＩＰＳは、動きチェッカＭＣＨをフィルタとして機能させることによって実装される。例えば、一実施形態では、動き破損されていない投影画像のみが、画像処理システムの次の段に転送される。患者の動きに起因する破損情報を含むものとして動きチェッカによって判断された投影画像は破棄され、ＩＰＳによってそれ以上処理されない。具体的には、ボリュームＶを再構成する際に再構成器によって使用されない。この実施形態では、動きチェッカＭＣＨは、投影画像のセット全体から動き破損されていない画像のみを実質的にフィルタリングし、これらのみが、３ＤボリュームＶに再構成されるように、再構成器ＲＥＣＯＮ（図３の点線）に渡される。この実施形態では、再構成器ＲＥＣＯＮは、全投影画像の元の集合のうちの真のサブセットからのみ再構成する。次いで、再構成されたボリュームは、出力ポートＯＵＴで出力された後、記憶されても、表示デバイス上での視覚化のために処理されても、他の方法で処理されてもよい。
【００６２】
しかしながら、好ましい実施形態では、動き破損された投影画像は破棄されずに、動き推定段ＭＥに渡されることによってさらに解析される。動き破損があるかどうかの指示を提供する動きチェッカＭＣＨとは対照的に、動き推定器は、考慮される投影画像に関連する経験された動きの推定値
【数１】
2022502108.tif 000003
を実際に計算する。
【００６３】
次いで、この動き推定値は、動き補償器ＭＣに転送され、そこで、動き破損された投影画像の補正動作が行われる。
【００６４】
これらの実施形態のいずれにおいても、動きチェッカの動作は、閾値に基づく。つまり、すべての動き破損が等しく重大であるとみなされるわけではない。閾値に基づいて、物体の十分に大きい動きの間に取得された投影画像のみが、動き破損されているとフラグが立てられる。
【００６５】
閾値化は、動き補償器によって使用され得る。動きチェッカによって計算された動き尺度は、閾値と比較される。例えば（ただし、これは慣例ではあるが）、特定の閾値を超えた場合にのみ、イメージング物体が十分に動いているときに起きたとみなされる。
【００６６】
動き補償器によって行われる補正動作は、動き破損されていることがわかったフレーム内の画像情報を変更することを含む。他の好ましい実施形態では、フレーム内の画像情報は変更されないが、推定された動き
【数２】
2022502108.tif 000004
に基づいて、そのフレームの関連する座標に対して座標変換
【数３】
2022502108.tif 000005
が行われる。このようにして、考慮される投影画像には、イメージングされた物体及び軌跡Ｔに関して空間内の新しい位置が割り当てられる。あるいは、
【数４】
2022502108.tif 000006
であり、したがって、ピクセル画像情報（ｕ，ｖ）は維持され、空間イメージングジオメトリ座標λのみがλ’に変更される。いくつかの実施形態では、変換
【数５】
2022502108.tif 000007
は、投影画像のうちの１つ以上の投影画像のジオメトリを再割り当てするように、画像を再分類する。頭部イメージングなどの一実施形態では、２つ以上の投影画像が
【数６】
2022502108.tif 000008
によって同じジオメトリに再割り当てされる。これは、例えば脳イメージングにおいて、例えば画像取得中に、患者ＰＡＴが頭部を左から右に向けたときに起こる可能性がある。したがって、２つの投影画像をλにより異なる方向から取得することができるが、これらは同じ情報、すなわちＲＯＩのビューを表す。
【００６７】
次いで、保持されたフレーム（もしあれば）と共に動き補正された投影画像の集合が、再構成器ＲＥＣＯＮに転送された後、出力ポートＯＵＴで出力される動き補償された３Ｄボリュームに再構成される。
【００６８】
いくつかの実施形態では、動き推定器ＭＥ、動き補償器ＭＣ及び再構成器ＲＥＣＯＮ間にフィードバックループがあってもよい。例えば、一実施形態では、推定器ＭＥによる動き推定は、３Ｄ／２Ｄ位置合わせを含む。この実施形態では、第１のパスで画像ボリュームＶ’を得るために、結果が動き破損され得ることを知りつつ、第１の再構成が行われる。次いで、中間ボリュームＶ’を順投影し、投影領域における投影画像と比較する。偏差がある場合、投影画像を再整列させるために位置合わせが行われ、そのように再整列された投影画像は次に、第２のパスで再構成器に再び転送される。次に、第２のパスにおいて、再構成器ＲＥＣＯＮは、補正された、すなわち再整列された投影画像から、残りの破損されていない投影画像（ある場合に）と共に、動き補償された最終３ＤボリュームＶを再構成する。
【００６９】
ここで、提案される画像処理システムＩＰＳの動作の実施形態を、図４Ａ、図４Ｂの曲線をさらに参照して説明する。図４Ａ、図４Ｂのグラフにおいて、それぞれのスコアＤＣＣは、例示的な番号付け１〜６００で示されるそれぞれの投影画像フレームにわたって捕捉される。投影画像の番号付けは、縦軸ＤＣＣに沿ったスコア値に対する横軸に沿って示される。円錐ビームＣＴ頭部スキャンのための中心平面のデータ整合性尺度が取られる。
【００７０】
図４Ａは、一実施形態による動きチェッカと動き推定器との相互作用を示す。動きチェッカは、いくつかの又は各フレームπについて、尺度ＤＣＣ（π）を計算する。尺度は、動き破損に敏感である。好ましくは、尺度ＤＣＣは、投影領域全体で計算される。すなわち、投影領域が埋め込まれている周囲の３Ｄ空間の「外へ出る（stepping outside）」ことは不要である。特に、スコアを計算するための３Ｄ−２Ｄ位置合わせは不要である。より一般的には、実施形態では、スコアは、ジオメトリ情報λ以外の、周囲の３Ｄ空間からの情報を使用せずに計算できる。
【００７１】
各投影画像フレームπについて計算されたスコアは、図４Ａの右斜線の幅広間隔のハッチングで設定されたエリアによって示されるように、閾値に対して処理され得る。スコアが閾値を超えるフレームπは次に、動き破損されているとしてフラグが立てられ、次に、これらのみが動き推定される。当然ながら、他の実施形態では、使用される特定の尺度に応じて、閾値を下回る低下が動き破損を示す逆の状況があってもよい。
【００７２】
動き破損されていることがわかっており、個々の動き推定が行われたフレームが、図４Ａの下部に示されている。このようなフレームは、左斜線幅広ハッチングによって示されている。これらのフレームについてのみ、それぞれの動き推定値に基づいて動き補償が行われる。残りの投影画像（右斜線幅細ハッチングで示される）は、変更されないままである。
【００７３】
図４Ａによるアプローチは有用であり、本明細書において、実施形態において想定されているが、図４Ａの実施形態の改良版と考えられる図４Ｂに示される第２の実施形態も想定される。この第２の実施形態では、動き推定は、動きによって破損されていることが分かっている特定のフレームに限定される。動きは、所与のフレームについて推定されると、他のフレームについて再利用され得る。具体的には、この実施形態は、イメージングされる物体が、左、真っ直ぐ、右などの人間の頭部の異なる姿勢又は向きなどの異なる物理的状態を取ることができる場合に特に有用である。物体は、ある状態になると、その状態に特定の期間の間留まる。例えば、頭部は、左又は右の特定の位置に静止する。ある状態になると、その状態に留める期間には動きはない。しかしながら、１つの状態から別の状態に移ると、例えば、頭部が左からまっすぐ又は右へと向けられると、動きが生じる。それぞれの状態は、動き破損スコアＤＣＣ曲線におけるプラトー領域ＰＬ１及びＰＬ２によって示される。プラトー領域は、異なるレベルを有していても有していなくてもよい。
【００７４】
図４Ｂは、物体がなり得る２つの異なる物理状態又はフェーズＰＨ１及びＰＨ２を示す例示的な状況を示す。あるフェーズになると、１つ、２つ、又はそれ以上の移行領域によって示されるように、別のフェーズに戻る、又は別のフェーズに移るまで、特定の期間の間、そのあるフェーズにおける動きはない。図４Ｂの例では、３つの移行領域ＴＲ１、ＴＲ２、ＴＲ３が示されている。
【００７５】
ここでは、このような移行期間ＴＲ１、ＴＲ又はＴＲ３の間に取得された投影画像に対してのみ、個々の動き推定を適用することが提案される。物体が（それぞれのプラトーＰＬ１及びＰＬ２によって示される）所与の同じ状態に留まっている間に取得された投影画像については、個々の動き推定は不要である。そのスコアがプラトーＰＬ１、ＰＬ２を形成する投影画像は、真っ直ぐのハッチングのセクションによって下部に示されている。
【００７６】
例えば、先行する移行フェーズＴＲ又はＴＲ２からの投影画像のうちの１つについて取得された動き推定値を、各後続のプラトー領域内のすべての投影画像に使用してもよい。より具体的には、また、プラトー領域ＰＬ２を参照すると、フレーム毎の個々の動き推定は、先行する移行フェーズＴＲ２における投影画像に対してのみ行われ得る。動き補償は移行領域について行われるが、動き推定は、後続プラトー領域ＰＬ２内のすべての投影画像に対してそのまま再利用され得る。これは、プラトー領域において個々の推定を行う必要がないため、計算時間を節約できる。ＲＯＩが所与の物理状態ＰＨ１、ＰＨ２にある間に取得された投影画像は、先行する移行領域で計算された１つの再利用動き推定に基づく同じ補正動作を使用して動き補償される。同様の考察は、ＰＬ１などの他のプラトー領域にも当てはまる。これらにおいても、先行する移行領域ＴＲ２において取得された投影画像に基づいて計算されたものと同じ動き推定値が使用される。フェーズのうちの１つ（ＰＨ１）は、動き推定を行う必要がない基準又は基底フェーズであることが理解されよう。基底フェーズＰＨ１で取得された投影画像は、図４Ｂの前述のように変更されないままにされる。移行領域からの動き推定値を再利用する場合、移行フェーズからの最新又は最後のそのような投影画像に対する推定値を使用することが好ましい。上記は、３つ以上のフェーズ及び任意の２つのこのようなフェーズの間の移行に容易に一般化され得ることが理解される。
【００７７】
実施形態では、一連の投影画像にわたって計算されたスコア値自体を表示デバイスＤＤ上に視覚化し得る。実施形態では、視覚化は、図４Ａ、図４Ｂに示されるように構成され得る。任意選択で、関与する閾値についての視覚的指示もまた、図４Ａ、図４Ｂに示されるように示される。スコア値の視覚化は、取得中に発生した動きの量に関するフィードバックをユーザに提供する。本明細書で提案されるシステムは、ユーザ対話を必要とせずに完全自動であるが、ユーザは依然として、適切なユーザインターフェースによって介入することができる。例えば、視覚化されたスコア値曲線に基づいて、ユーザは、場合によってはシステムＩＰＳを無効にして、特定のフレームに個別の動き補正を適用することを決定することができる。
【００７８】
ここで、画像処理システムの動作の基礎となる画像処理方法のフロー図を示す図５を参照する。しかしながら、以下のステップは、それら自体で教示として理解され、図２の画像処理システムＩＰＳのアーキテクチャに必ずしも結びついていないことが理解されるであろう。
【００７９】
ステップＳ５１０において、Ｘ線イメージング機器から投影画像が受信される。これらは、関心領域に対して異なる方向から取得される。あるいは、投影画像は、取得後にメモリに記憶され、後にそこから取り出されてもよい。しかしながら、好ましくは、投影画像は、イメージャＸＡから適切なインタフェースＩＮを介してリアルタイムで受信される。投影画像は、ピクセル（ｕ，ｖ）での画像情報を含み、また、取得に使用されたイメージングジオメトリに関する空間情報λを含む。所与の投影画像の空間座標は、領域に対するＸ線源の空間における位置を表す。これは、投影画像が取得されたときにＸ線源がとっていた位置である。空間情報は、各投影画像に関連する空間座標を含み得る。空間座標は、取得時間における各投影画像の位置とグリッド点（ボクセル位置）との間の空間関係を記述し得る。コーンビームＣＴのイメージングでは、図２に示すジオメトリが使用されてもよい。患者がグリッド点に対して静止したままであるという仮定もあり得る。少なくとも、例えば、患者を患者テーブルと共にＸ線システムに対して移動させるなど、起こるあらゆる動きは、事前に規定されているか又は知られているという仮定があり得る。これらの仮定のいずれかが満たされない場合、いくつかの又はすべての投影画像は、患者の動きに起因して動き破損されている可能性がある。
【００８０】
ステップＳ５２０において、投影画像が処理され、いくつかの又は各投影画像について尺度又はスコアが計算される。スコアは、所与の投影画像の動き破損があったか否か、すなわち、各投影画像が取得されたときに患者の動きがあったかどうかを判断するために使用される。実施形態では、動き破損があったか否かの判断は、スコアを１つ以上の閾値と比較することによる閾値化に基づいている。尺度を計算するステップＳ５２０は、好ましくは、ステップＳ５１０において受信された全投影画像のサブセットに対して行われる。
【００８１】
実施形態では、１つ以上の投影画像を使用して、所与の１つの投影画像の動き破損スコアを計算し得る。１つ以上の投影画像は、所与の投影画像が参照画像と整合するかどうか、したがって、所与の画像が動き破損されているか否かを判断するために、参照画像を提供できる。実施形態では、上述のＬｅｓａｉｎｔの論文など、他の場所で報告されているように、スコアデータ整合性条件ＤＣＣを使用することが想定される。ＤＣＣは、Ｌｅｓａｉｎｔなどによって示されるように、再構築（reformulation）を使用することによって効率的に評価できる。
【００８２】
ステップＳ５２０において、動きに敏感である限り、任意の適切な尺度を使用してスコアを計算することができる。そのような適切な尺度の１つのクラスは、上述のＬｅｓａｉｎｔの参照文献にあるように、過去に提案されたＤＣＣであることがわかっている。しかしながら、例示的である一実施形態では、以下のように動機付けされ得る以下の尺度Ｇを使用してもよい：ｇは、減衰物体ｆの線源軌道λに沿って測定されたＸ線コーンビーム変換を示すとする。関数ｆは、３Ｄスカラー場であり、材料密度、散乱、回折（位相コントラスト）などの、再構成が望まれる（一般に）空間依存品質を記述する。線源位置λに対して測定されるコーンビームデータπ＝ｇは、取得の終わりにおいて、次式で与えられる：
【数７】
2022502108.tif 000009
（１）において、
【数８】
2022502108.tif 000010
は線源ＸＳの位置であり、
【数９】
2022502108.tif 000011
は線源ＸＳの位置から検出器平面Π（Ｄ）までの線であり、それに沿って物体ｆ（ＲＯＩを表す）の線積分が取られる。
【００８３】
実施形態では、冗長情報に基づき、冗長性が維持されているかどうかを確立することを可能にする尺度が使用される。より詳細には、冗長情報は、より複雑な形式で、例えば、測定されたコーンビーム投影ｇのフィルタリングされた２Ｄラドン変換として表すことができる。本明細書では、実施形態において、Ｇと呼ばれる関数表現を用いて、イメージングジオメトリの仮定に起因する情報の冗長性が確実に維持されているかどうかを測定することを提案する。実施形態では、同一の情報を含むと予想されるラドン平面内の情報が測定される。情報が同一ではなく、逸脱している場合、これはデータ（不）整合性の尺度とみなされる。ラドン平面は、ラドン変換が画像領域からの点をマッピングする平面である。言い換えれば、ラドン平面は、ラドン変換の関数画像である。
【００８４】
より詳細には、実施形態において、データ整合性は、関数Ｇ（Ｌｅｓａｉｎｔでは同様の関数が中間関数と呼ばれている）を使用して測定され得る：
【数１０】
2022502108.tif 000012
ここで、
【数１１】
2022502108.tif 000013
は、中心平面などの平面を決定する３Ｄ空間における単位球面Ｓ２上のベクトルである。上述のＬｅｓａｉｎｔの論文によれば、
【数１２】
2022502108.tif 000014
と３Ｄラドン変換
【数１３】
2022502108.tif 000015
との間の関係は、次式：
【数１４】
2022502108.tif 000016
によって与えられる。
【００８５】
次いで、畳み込みカーネルｈは、例えば、ランプフィルタ又は微分フィルタとして選択され得る：
【数１５】
2022502108.tif 000017
【００８６】
少なくとも２つの線源位置
【数１６】
2022502108.tif 000018
及び
【数１７】
2022502108.tif 000019
を含む、原点に対するその法線
【数１８】
2022502108.tif 000020
及びその距離ｓによって定義される平面が、２つの投影
【数１９】
2022502108.tif 000021
及び
【数２０】
2022502108.tif 000022
の整合性尺度を提供する。これは、
【数２１】
2022502108.tif 000023
及び式（２．１）から、以下の恒等式：
【数２２】
2022502108.tif 000024
を導出することができるからである。
【００８７】
したがって、恒等式（４）を使用して、
【数２３】
2022502108.tif 000025
のいずれか１つなど、必要に応じて、乗法的又は加法的な形成で、動きに敏感な尺度又はスコアを構成することができる。
【００８８】
（５）又は（６）の組み合わせ、又は（４）から導出可能な他の関数式も想定される。（６）では、適切なノルム||.||を使用することができる。定数係数だけ（５）、（６）と異なる任意の尺度は、本明細書では、代替実施形態において想定されていると理解される。関係（４）、したがって（５）、（６）は、適切な検出器座標変換を使用することによって数値的に効率的に評価され得る。特に、Ｇは、上述の参照文献の４３２頁のセクションＡ）においてＬｅｓａｉｎｔによって報告されているように、投影データの事前重み付け、２Ｄラドン変換及び畳み込みを含む３ステップ処理として、再構築され得る。（４）、（５）、（６）のとおり、計算は、λ１における所与の投影画像のスコアを計算するために、１つ以上の投影画像（λ２）において必要となり得る。しかしながら、所与の投影画像、例えば第１の投影画像が参照画像として選択され、式（４）〜（６）は、１変数依存性に縮小することができる。参照投影画像は、図４Ｂのプラトー領域を表す静止フェーズ／状態の間に取得されたものであってもよい。
【００８９】
式（４）〜（６）によるＤＣＣは、単に１つの例示的な実施形態であり、他の式において及び／又は円形以外のジオメトリについての他のＤＣＣも、本明細書では想定されることが理解されるであろう。特に、他の実施形態では、非ラドン平面ベースのデータ整合性尺度を代わりに使用してもよい。また、本明細書で想定されるすべての尺度が冗長情報に基づいているわけではない。
【００９０】
一般に、ＤＣＣは、ＤＣＣが正常値又は値の範囲から逸脱する場合、動き破損を検出するために有用であり得る。上記正常値又は範囲に基づいて適切に構成された閾値を使用して、投影画像のセットについて計算されたスコアに基づいてそのような偏差を検出することができる。ＤＣＣ又は他の尺度の正常値は、即座には知られていないが、１つ以上の較正スキャンで得ることができる。次いで、ＤＣＣ関数によるスコア値を一連の投影画像にわたって取得して、ＤＣＣ値がどこに落ち着くかの見当を得ることができ、これは、動きがないことを示す正常値を表すものとみなされる。実際に動きがないことを確かめるために、較正においてファントムを使用することができる。より具体的には、Ｓｈｅｐｐ−Ｌｏｇａｎファントムを使用することができる。
【００９１】
ステップＳ５２０におけるスコアの計算は、投影領域における情報に基づいて排他的に行われ得る。このステップでは、周囲の３Ｄ空間における演算は不要である。より具体的には、このステップでは、３Ｄ−２Ｄ位置合わせは不要である。
【００９２】
実施形態では、ステップＳ５２０で計算された動き破損の量を表すスコア値は、例えば、図４Ａ、図４Ｂに示されるように、適切な表示デバイスにおいてユーザにグラフ表現として表示されるように使用され得る。
【００９３】
ステップＳ５２０における動き破損があるかどうかの判断は、後続の処理フローを決定する。一実施形態では、（閾値と比較された）スコアが、動き破損があったことを示す場合、対応する投影画像は破棄され、ステップＳ５５０における後続の再構築において使用されなくなる。この実施形態では、ＲＯＩの３Ｄボリュームを再構成するために、再構成ステップＳ５５０において、利用可能な全投影画像のうちのサブセットのみが使用される。
【００９４】
代替実施形態では、投影画像は、動き破損されていることがわかっていても依然として使用される。これらの実施形態では、投影画像は、動きの推定が行われるステップＳ５３０に渡される。動き破損されていないことが分かっている投影画像は、そのまま保持され、ステップＳ５５０の再構成までリープフロッグ（飛び越し）される。
【００９５】
ステップＳ５３０において、破損を引き起こした動きの推定値が計算される。このステップには、動きの大きさ及び／又は方向を計算することが含まれ得る。動きは、剛体変換、特に３Ｄ剛体変換を計算することによって推定され得る。
【００９６】
ステップＳ５３０における動き推定は、ステップ５２０における判断が動き破損があるとの場合にのみ、所与の投影画像に適用される。
【００９７】
このようにして得られた動き推定値は、後続のステップＳ５４０において使用されて、実際の動き補償が実行される。これについては、多くの異なる実施形態が想定され、以下でより詳細に説明する。
【００９８】
ステップＳ５４０は動き補償された投影画像をもたらし、この投影画像は、再構築のためにステップＳ５５０に転送される。当初に受信された投影画像の全セットは、ステップＳ５２０において十分な動き破損が見つからなかったためにそのままにされた投影画像と、ステップＳ５４０において動き補償ステップが実行された投影画像との２つのグループに属する。次に、２つのグループの投影画像、すなわち、元の投影画像と動き補償された投影画像は共に、再構成ステップＳ５５０において使用されて３ＤボリュームＶが生成され、これは、次いで、記憶、表示、又はさらなる処理のために、ステップＳ５６０で出力され得る。
【００９９】
上述のステップＳ５２０〜Ｓ５４０は、ステップＳ５１０で受信された異なる方向からの投影画像の一部又は全部について繰り返されると理解される。
【０１００】
ここで、動き補償ステップＳ５４０をより詳細に再び参照すると、上述のいくつかの実施形態についてより詳細に説明する。実施形態では、動き補償における関連する補正動作は、投影画像におけるピクセル画像情報（ｕ，ｖ）を変更することを含み得る。具体的には、画像情報は、例えば、オプティカルフロー法などによって、補償を行うように変更される。しかしながら、他の実施形態では、各投影画像に関連付けられた空間イメージングジオメトリ座標λのみが変更される。この変更は、上述したように、患者の動きを考慮するために、座標変換
【数２４】
2022502108.tif 000026
で行われてよい。特に、１つ以上の投影画像の変換された空間イメージングジオメトリ座標に基づく一実施形態では、仮想軌跡Ｔ’が計算され、再構成はその新しい仮想軌跡に基づいて行われる。実際の線源軌跡ＴがＲＯＩに対して円弧だったとしても、動きがあることにより、軌跡はもはや円形を記述せず、より複雑な形状を有し、この複雑な形状が、ステップＳ５５０において、より正確な再構成をもたらすために、円形の代わりに使用され得る。他の実施形態では又はそれに加えて、２つのパスを有する３Ｄ−２Ｄ位置合わせが、動き破損されている投影画像に対して行われる。
【０１０１】
好ましくは、尺度ｍｅａｓ（．）に基づいて、動き破損があったかどうかを判断するためのステップＳ５２０は、好ましくは、後続の動き推定ステップＳ５３０及び／又は動き補償ステップＳ５４０ののいずれよりも計算コストが安価である。特に、動きチェッカステップＳ５２０は主に、３Ｄ計算を含まないものとして、特に、３Ｄ順投影が不要であるものとして想定されている。このようにして、計算コストの高い動き推定ステップＳ５３０及び／又は動き補償ステップＳ５４０が、動き破損によって実際に十分に深刻な影響を受けている投影画像にのみ行われるので、計算時間が大幅に削減される。さらに、ＤＣＣ条件式（４）〜（６）などの尺度ｍｅａｓ（）に関する計算の一部は、予想イメージングジオメトリのみに依存し、事前に分かっているので、事前に計算できる。この場合、投影データにｍｅａｓ（．）を適用するステップＳ５２０を迅速に行うことができ、その結果、判断ステップＳ５２０でｍｅａｓ（．）を計算するときの応答が早くなり、ＣＰＵ時間が短縮される。
【０１０２】
実施形態では、また、図４Ａで説明したように、ステップＳ５２０において閾値化が行われて、どの投影画像に対して動き推定Ｓ５３０及び動き補償Ｓ５４０を行うべきかが判断される。しかしながら、この実施形態は、患者又はＲＯＩの様々な状態を識別することを可能にする、２つ以上の複数の閾値を使用することによって、さらに改良できる。具体的には、関心の状態が仮定されると、患者は、比較的穏やかなままであり、これらの状態の間の移行時にのみ、動き推定、したがって動き補償が実行される。この実施形態では、取得された投影画像の時系列のスコア値を分析して、静止状態を表すプラトー領域を特定する必要がある。静止状態は、実質的にゼロ勾配のセクションによって示され得る。異なるプラトーレベルに対応する２つ以上の異なる状態があってよい。プラトーレベルは、プラトーを見つけるために時系列のスコア値に対して自動識別器を実行することによって見つけることができる。
【０１０３】
次に、ステップＳ５３０における動き推定が、移行を表す非消失勾配があるところに行われる。図４Ｂによる本実施形態のさらなる改良において、プラトー領域に関連する投影画像に適用される実際の動き補償Ｓ５４０は、個別に処理される必要がない。言い換えると、同じ動き補償Ｓ５４０を、プラトー領域に関連するすべての投影画像に適用できる。所与のプラトー領域に関連する投影画像は、それら自体間で整合性がある。ただし、異なるプラトー領域について取られた投影画像についてはそのような整合性はない場合がある。したがって、１つのプラトーＰＬ１についての投影画像は、別のプラトー領域ＰＬ２についての投影画像に対して、動き破損されていると判断され得る。この場合、１つのプラトー／状態についての投影画像に動き補償を適用することにより、別の（参照）プラトー又は状態からの投影画像との整合性を確立することが可能になる。
【０１０４】
移行中に取得された投影画像について計算された動き推定値のいずれも、その移行の直後に続くプラトー領域に使用してもよい。より具体的には、移行領域において計算された最新の動き推定値を実質的にすべてのプラトー投影画像に使用してもよい。この場合、動き補償も、同じ動き推定値に基づいて、プラトー投影画像について同じであり得る。この実施形態は、図４Ａのように、動き破損の影響を受けていることがわかった各投影画像に対して、個別の動き推定／補償が行われる実施形態と比較して、さらに一層計算時間が削減できる。図４Ｂの実施形態では、（移行領域内の投影画像について推定される）同じ動き推定が、後続のプラトー領域のすべての投影画像に使用され得る。個別の推定を行う必要がなく、前の推定が再利用される。プラトー領域の投影画像には、動き補償だけを適用する必要があるが、これは同じ動き推定に基づいているので、同じ動き補償である。
【０１０５】
実施形態では、提案される方法及びシステムは、神経イメージングにおける心臓バリスティック効果（Ｃａｒｄｉｏ−ｂａｌｌｉｓｔｉｃ  ｅｆｆｅｃｔ：ＣＢＥ）を検出及び補正するために使用され得る。ＣＢＥは、心拍活動及び脈拍動によって引き起こされる動きに関連する。上述の原理をＣＢＥタイプの動きに適用できる。動き推定及び補償は、非剛体変換を含んでもよい。
【０１０６】
投影画像πがさらに、切り詰められた投影も含む場合、ステップＳ５２０における整合性計算は、動きによって引き起こされる外乱ではなく、切り詰めから生じる追加の不整合性効果を回避するために、切り詰められていない平面にのみ適用され得る。良好な実行のために、本明細書では、これらの２つのタイプの不整合性の混合を回避することが提案される。あるいは、切り詰めによって引き起こされる不整合性も考慮するために、整合性尺度が増強される。
【０１０７】
提案される方法はまた、例えば（比較的少ない）投影からの３Ｄサーフェイスモデリングのために、非断層撮影再構成において整合性のある投影を特定するために使用されてもよい。
【０１０８】
ＣＢＣＴについて特に参照したが、提案される選択的動き補正のコーンビーム以外の他のＸ線イメージングモダリティへの適用も想定され、本明細書に説明される原理はそれへの直接的な適用を許容する。
【０１０９】
本明細書で開示される１つ以上の特徴は、コンピュータ可読媒体内にコード化された回路として／それを用いて、及び／又はそれらの組合せとして構成又は実装され得る。回路は、ディスクリート及び／又は集積回路、特定用途向け集積回路(ＡＳＩＣ)、システムオンチップ(ＳＯＣ)、及びそれらの組合せ、機械、コンピュータシステム、プロセッサ及びメモリ、コンピュータプログラムを含み得る。
【０１１０】
本発明の別の例示的な実施形態では、適切なシステム上で、前述の実施形態のうちの１つによる方法の方法ステップを実行するように構成されることを特徴とするコンピュータプログラム又はコンピュータプログラム要素が提供される。
【０１１１】
したがって、上記コンピュータプログラム要素は、本発明の実施形態の一部であり得るコンピュータユニットに格納されてもよい。この計算ユニットは、上述の方法のステップを実行し、又は該ステップの実行を誘導するように構成され得る。また、上述の装置の各構成要素を動作させるようにしてもよい。計算ユニットは、自動的に動作するように、及び／又はユーザの指令を実行するように構成されてもよい。コンピュータプログラムは、データプロセッサの作業メモリにロードされてもよい。したがって、データプロセッサは、本発明の方法を実行するように装備され得る。
【０１１２】
本発明のこの例示的な実施形態は、最初から本発明を使用するコンピュータプログラムと、アップデートによって既存のプログラムを本発明を使用するプログラムに変えるコンピュータプログラムとの両方を対象とする。
【０１１３】
さらに、コンピュータプログラム要素は、上述の方法の例示的な実施形態の手順を満たすために必要なすべてのステップを提供することができる。
【０１１４】
本発明のさらなる例示的な実施形態によれば、ＣＤ‐ＲＯＭなどのコンピュータ可読媒体が提示される。コンピュータ可読媒体は、その上に格納されたコンピュータプログラム要素を有し、コンピュータプログラム要素は、前述の段落によって説明されている。
【０１１５】
コンピュータプログラムは、他のハードウェアと共に又はその一部として供給される光記憶媒体又はソリッドステート媒体などの適切な媒体(特に、必ずしもそうである必要はないが、非一時的媒体)上に記憶及び/又は配布することができるが、インターネット又は他の有線もしくは無線電気通信システムなどを介して他の形態で配布することもできる。
【０１１６】
しかしながら、コンピュータプログラムはまた、ワールドワイドウェブのようなネットワークを介して提示されてもよく、このようなネットワークからデータプロセッサの作業メモリにダウンロードすることができる。本発明のさらなる例示的な実施形態によれば、ダウンロードのためにコンピュータプログラム要素を利用可能にするための媒体が提供され、このコンピュータプログラム要素は、本発明の前述の実施形態のうちの１つによる方法を行うように構成される。
【０１１７】
本発明の実施形態は、異なる主題を参照して説明されていることに留意されたい。特に、方法タイプの請求項を参照して説明される実施形態もあれば、デバイスタイプの請求項を参照して説明される実施形態もある。しかしながら、当業者は、上記及び以下の説明から、別段の通知がない限り、１つのタイプの主題に属する特徴の任意の組合せに加えて、異なる主題に関する特徴間の任意の組合せも、本出願で開示されると考えられることを理解するであろう。しかしながら、すべての特徴を組み合わせて、特徴の単なる寄せ集め以上の相乗効果を提供できる。
【０１１８】
本発明を、図面及び前述の説明において詳細に図示及び説明したが、このような図示及び説明は、例示的であり、限定的ではないと考えられるべきである。本発明は、開示された実施形態に限定されない。開示された実施形態に対する他の変形は、図面、開示及び従属請求項の精査から、特許請求されている発明を実施する際に当業者によって理解され、達成され得る。
【０１１９】
特許請求の範囲において、「含む」という用語は、他の要素又はステップを排除するものではなく、単数形は、複数を排除するものではない。単一のプロセッサ又は他のユニットが、特許請求の範囲に記載されるいくつかのアイテムの機能を果たすことができる。特定の手段が相互に異なる従属請求項に記載されているという単なる事実は、これらの手段の組合せを有利に使用することができないことを示すものではない。特許請求の範囲におけるいかなる参照符号も、範囲を限定するものとして解釈されるべきではない。
【図１】
2022502108.tif 000027
【図２】
2022502108.tif 000028
【図３】
2022502108.tif 000029
【図４Ａ）】
2022502108.tif 000030
【図４Ｂ）】
2022502108.tif 000031
【図５】
2022502108.tif 000032
【国際調査報告】
2022502108.tif 000033
2022502108.tif 000034
2022502108.tif 000035