(57)【要約】  非常に多くの場合、画像レジストレーションは、手動で実行されなければならない退屈な作業であった。本発明の１つの実施例によれば、画像時系列のレジストレーションは、薬物動態学的モデルに基づいて関心領域の代替変換列が互いに比較される、薬物動態学的モデルに基づいて実行され、最高の変換ベクトル列は、画像レジストレーションのために用いられる。有利なことに、これは、たとえ解剖学的コントラストが全く又はほとんどない場合でも、臓器運動の効果的な補償を可能にすることができる。
2007536054.tif 000002
【特許請求の範囲】
【請求項１】
少なくとも第１の画像及び第２の画像を有する画像の時系列をレジスタする方法において、
前記第１の画像中で第１の関心領域を選択するステップと；
前記第２の画像中で第２の関心領域を選択するステップと；
薬物動力学モデルに基づいて前記第１の関心領域から前記第２の関心領域への第１の変換を決定するステップと；
前記第１の変換に基づいて前記第１の画像及び前記第２の画像をレジスタするステップと；
を有し、前記第１の関心領域は、前記第２の関心領域に対応する、方法。
【請求項２】
請求項１に記載の方法において、前記第１の変換の前記決定ステップは、
前記第１の関心領域から前記第２の画像中で選択される第３の関心領域への第２の変換に基づいて、且つ、前記薬物動力学モデルに基づいて、第１の薬物動力学パラメータを決定するステップと；
前記第１の関心領域から前記第２の画像中で選択される第４の関心領域への第３の変換に基づいて、且つ、前記薬物動力学モデルに基づいて、第２の薬物動力学パラメータを決定するステップであって、前記第３の変換は、前記第２の変換のバリエーションであり、前記第３の及び前記第４の関心領域は、前記第１の関心領域に対応する、ステップと；
第１のコンパートメントパラメータ及び第２のコンパートメントパラメータの品質決定を実行し、第１の品質値及び第２の品質値を結果として生じるステップと；
前記第１の及び第２の品質値のうちどちらがより良い品質値であるかを決定するステップと；
前記第２の変換及び前記第３の変換から変換を選択するステップであって、当該選択される変換は、前記より良い品質値に対応する変換である、ステップと；
を有し、前記選択される変換は前記第１の変換である、方法。
【請求項３】
請求項２に記載の方法において、前記品質決定は、
前記第１の薬物動力学パラメータの推定及び前記第２の薬物動力学パラメータの推定のうちの少なくとも１つに基づく統計パラメータの決定と、
薬物動力学パラメータのライブラリと、
前記第１の薬物動力学パラメータ推定及び前記第２の薬物動力学パラメータ推定の整合性と、
のうちの少なくとも１つに基づいて実行される、方法。
【請求項４】
請求項１に記載の方法において、前記第１の及び第２の関心領域のうちの少なくとも１つ並びに前記薬物動力学モデルは、選択肢の既定の組からインタラクティブに選択される、方法。
【請求項５】
請求項２に記載の方法において、前記第１の品質値及び前記第２の品質値のうちの少なくとも１つが予め設定された閾値を上回るまで、反復的に繰り返される方法。
【請求項６】
請求項１に記載の方法において、医療イメージングにおいてＣＴデータセット、ＭＲＩデータセット、ＰＥＴデータセット、ＳＰＥＣＴデータセット及び超音波イメージングデータセットのうちの１つに適用される、方法。
【請求項７】
第１の画像及び第２の画像をレジスタするための画像処理装置において、
前記第１の画像及び前記第２の画像を有する多次元データセットを記憶するためのメモリと；
前記多次元データセットをロードする動作と、
前記第１の画像中で第１の関心領域を選択する動作と、
前記第２の画像中で第２の関心領域を選択する動作と、
薬物動力学モデルに基づいて前記第１の関心領域から前記第２の関心領域への第１の変換を決定する動作と、
前記第１の変換に基づいて前記第１の画像及び前記第２の画像をレジスタする動作と、
を実行するように適応され、前記第１の関心領域は、前記第２の関心領域に対応する、画像プロセッサと；
を有する装置。
【請求項８】
第１の画像及び第２の画像をレジスタするためのスキャナシステムにおいて、
前記第１の画像及び前記第２の画像を有する多次元データセットを記憶するためのメモリと；
前記多次元データセットをロードする動作と、
前記第１の画像中で第１の関心領域を選択する動作と、
前記第２の画像中で第２の関心領域を選択する動作と、
薬物動力学モデルに基づいて前記第１の関心領域から前記第２の関心領域への第１の変換を決定する動作と、
前記第１の変換に基づいて前記第１の画像及び前記第２の画像をレジスタする動作と、
を実行するように適応され、前記第１の関心領域は、前記第２の関心領域に対応する、前記第１の画像及び前記第２の画像のレジストレーションを実行するように適応された画像プロセッサと；
を有するシステム。
【請求項９】
ＣＴスキャナシステム、ＭＲＩスキャナシステム、ＰＥＴスキャナシステム、ＳＰＥＣＴスキャナシステム及び超音波イメージングシステムのうちの１つである請求項８に記載のスキャナシステム。
【請求項１０】
第１の画像及び第２の画像をレジスタするためのコンピュータプログラムであって、
当該コンピュータプログラムが画像プロセッサにおいて実行されると、当該コンピュータプログラムは当該画像プロセッサに、
前記第１の画像及び前記第２の画像を有する多次元データセットをロードする動作と、
前記第１の画像中の第１の関心領域を選択する動作と、
前記第２の画像中の第２の関心領域を選択する動作と、
薬物動力学モデルに基づいて前記第１の関心領域から前記第２の関心領域への第１の変換を決定する動作と、
前記第１の変換に基づいて前記第１の画像及び前記第２の画像をレジスタする動作と、
を実行させ、
前記第１の関心領域は、前記第２の関心領域に対応する、コンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、デジタルイメージングの分野に関し、例えばメディカルイメージングの分野に関する。特に、本発明は、少なくとも第１の画像及び第２の画像を有する画像の時系列（time series）をレジスタ（register）する方法、画像処理装置、スキャナシステム、並びに、第１の画像及び第２の画像をレジスタするためのコンピュータプログラムに関する。
【背景技術】
【０００２】
同じ物体の２つの画像が、異なった投影（projections）、異なった時点、又は、関心物体の異なった移動段階の最中に、撮られなければならないとき、画像をレジスタすることは、非常に望ましい可能性がある。
【０００３】
画像をレジスタすることは、当該画像の幾何学的特性の共通基準系への統合を意味する。通常、これは、対応する解剖学的グレー値構造が何らかの類似性測度（相互相関、相互情報量等）とマッチするまで、画像を歪める（warp）ことによってなされる。核医療画像のレジストレーションは、トレーサの特異性が増加するにつれて、より困難になる。トレーサが特異的になるほど、その組織内取込み（tissue up-take）は普遍的でなくなり、従って、該トレーサが生じる解剖学的背景コントラストはより低くなり（解剖学的コントラストが失われるという問題）、従来型のグレー値に基づくレジストレーションを信頼できないものにする。
【発明の開示】
【発明が解決しようとする課題】
【０００４】
本発明の目的は、改良された画像レジストレーションを提供することである。
【課題を解決するための手段】
【０００５】
請求項１に記載の本発明の１つの実施例によれば、上記目的は、第１の画像及び第２の画像の画像時系列をレジスタする方法であって、前記第１の画像中で第１の関心領域（region of interest）が選択され、前記第２の画像中で第２の関心領域が選択される方法によって解決されることができる。更に、第１の関心領域の第２の関心領域への第１の変換（translation）は、薬物動力学モデル（pharmacokinetic model）に基づいて決定される。第１の画像及び第２の画像は第１の変換に基づいてレジスタされ、ここで、第１の関心領域は第２の関心領域に対応する。
【０００６】
例えば、本発明のこの実施例によれば、少なくとも２つの画像を有する画像時系列の局所的レジストレーションは、特定の期間に亘る特定の関心領域（例えばガン組織）の変換を決定することによって実行される。有利なことに、本発明のこの実施例によれば、変換は薬物動力学モデルに基づいて決定され、従って、解剖学的コントラストがほとんど又は全くない画像（高度に特異的なトレーサが核医学／メディカルイメージングにおいて用いられるときに頻繁に起こる）についての関心領域のトラッキングを可能にする。
【０００７】
請求項２に記載の本発明の他の実施例によれば、第１の変換は、第１の関心領域から第２の画像中で選択される第３の関心領域への第２の変換に基づいて、且つ、薬物動力学モデルに基づいて、第１の薬物動力学パラメータを決定することによって識別される。更に、第２の薬物動力学パラメータは、第１の関心領域から第２の画像中で選択される第４の関心領域への第３の変換に基づいて、且つ、薬物動力学モデルに基づいて、決定される。第３の変換は、第２の変換のバリエーションであり、第３の及び第４の関心領域は、第１の関心領域に対応する。更に、第１のコンパートメントパラメータ及び第２のコンパートメントパラメータの品質決定が実行され、第１の品質値及び第２の品質値を結果として生じる。次に、本発明のこの実施例によれば、これら第１の及び第２の品質値は比較され、第１の及び第２の品質値のうちどちらがより良い品質値であるかが決定される。変換は、第２の変換及び第３の変換から選択される。選択される変換はより良い品質値に対応し、ここで、選択される変換は第１の変換である。
【０００８】
有利なことに、本発明のこの実施例によれば、関心物体の異なった（代替の）変換が、対応する薬物動力学パラメータを、異なった変換に基づいて薬物動力学モデルを利用することによって推定することによって比較される。次に薬物動力学パラメータは品質化され、結果として対応する品質値になる。このとき、「より良い」変換（「より良い」品質値に対応する変換）が、画像表示のために用いられる。
【０００９】
有利なことに、これは改良された画像レジストレーションを可能にする可能性がある。
【００１０】
本発明の他の実施例は、請求項３に記載されており、ここで、品質決定は、第１の薬物動力学パラメータ推定及び第２の薬物動力学パラメータ推定のうちの少なくとも１つに基づく統計的品質測度の決定と、薬物動力学パラメータのライブラリと、第１の薬物動力学パラメータ推定及び第２の薬物動力学パラメータ推定の整合性とのうちの少なくとも１つに基づいて実行される。
【００１１】
有利なことに、これは、速いか、効率的であるか、又は自動的でありさえする品質決定を許可しうる。
【００１２】
請求項４に記載された本発明の他の実施例によれば、第１の及び第２の関心領域のうちの少なくとも１つ並びに薬物動力学モデルは、選択肢の既定の組からインタラクティブに選択される。従って、本発明のこの実施例によれば、ユーザは、画像取得の少し後に、又は画像取得の最中に、病巣候補（candidate lesion）及び適用されるべき薬物動力学モデルを選択することができる。
【００１３】
有利なことに、これは速くてユーザーフレンドリーなインタラクティブ画像レジストレーションを可能にすることができる。
【００１４】
請求項５に記載された本発明の他の実施例によれば、画像レジストレーションの方法は、第１の品質値及び第２の品質値のうちの少なくとも１つが予め設定された閾値を上回るまで、反復的に繰り返される。
【００１５】
請求項６に記載された本発明の他の実施例によれば、画像レジストレーションの方法は、医療イメージングにおいてＣＴデータセット、ＭＲＩデータセット、ＰＥＴデータセット、ＳＰＥＣＴデータセット及び超音波イメージングデータセットデータセットのうちの１つに適用される。
【００１６】
請求項７に記載された本発明の他の実施例によれば、第１の画像及び第２の画像をレジスタするための画像処理装置が提供され、該装置は、第１の画像及び第２の画像を有する多次元データセットを記憶するためのメモリと；前記多次元データセットをロードする動作と、前記第１の画像中で第１の関心領域を選択する動作と、前記第２の画像中で第２の関心領域を選択する動作と、薬物動力学モデルに基づいて前記第１の関心領域から前記第２の関心領域への第１の変換を決定する動作と、前記第１の変換に基づいて前記第１の画像及び前記第２の画像をレジスタする動作とを実行するように適応された画像処理プロセッサとを有する。第１の関心領域は、第２の関心領域に対応する。
【００１７】
有利なことに、本発明のこの実施例による画像処理装置は、改善された画像レジストレーション速度及び高いレジストレーション精度を可能にすることができる。
【００１８】
本発明は、第１の画像及び第２の画像を有する多次元データセットを記憶するためのメモリと、第１の画像及び第２の画像のレジストレーションを実行するように適応された画像処理プロセッサとを有するスキャナシステムにも関する。本発明の一側面によれば、スキャナシステムは、ＣＴスキャナシステム、ＭＲＩスキャナシステム、ＰＥＴスキャナシステム、ＳＰＥＣＴスキャナシステム及び超音波イメージングシステムのうちの１つである。本発明によるスキャナシステムは、請求項８及び９に記載される。
【００１９】
有利なことに、これは、本発明によるスキャナシステムによって得られる画像の時系列の改良された画像レジストレーションを可能にすることができる。
【００２０】
本発明は、例えば、画像処理プロセッサ等のプロセッサにおいて実行されることができるコンピュータプログラムにも関する。この種のコンピュータプログラムは、例えば、ＣＴスキャナシステム、ＭＲＩスキャナシステム、ＰＥＴスキャナシステム、ＳＰＥＣＴスキャナシステム又は超音波システムの一部であってもよい。本発明の１つの実施例によるコンピュータプログラムは、請求項１０に記載される。好適には、これらのコンピュータプログラムは、画像処理プロセッサの作業メモリにロードされることができる。これにより、画像処理プロセッサには、本発明の実施例を実行するための能力が与えられる。コンピュータプログラムは、ＣＤ−ＲＯＭ等のコンピュータ可読媒体に保存されることができる。コンピュータプログラムは、また、ＷｏｒｌｄＷｉｄｅＷｅｂ等のネットワークを通じて与えられることができ、この種のネットワークから画像処理プロセッサの作業メモリにダウンロードされることができる。本発明のこの実施例によるコンピュータプログラムは、いかなる適切なプログラム言語によっても（例えばＣ＋＋）書き込まれることができる。
【００２１】
以下は、本発明の１つの実施例の要旨とみなされることができる。即ち、画像時系列のレジストレーションは、コンパートメントモデル等の薬物動力学モデルに基づいて実行され、当該モデルにおいて、画像取得の間、例えば患者移動に起因して移動する関心領域の代替の変換シーケンスは、薬物動力学モデルに基づいて互いに比較され、最高の変換ベクトルが画像レジストレーションのために用いられる。有利なことに、これはたとえ解剖学的コントラストが全く又はほとんどない場合でも、臓器運動の有効な補償を可能にすることができる。有利なことに、本発明の１つの実施例によれば、閾値が導入されることができ、異なった可能な変換ベクトルを比較する工程が、それぞれの変換ベクトルの品質に対応する品質値が閾値を上回るまで、反復的に繰り返されることができる。
【００２２】
本発明のこれらの及び他の側面は、以下で説明される実施例を参照して説明され明らかにされる。
【発明を実施するための最良の形態】
【００２３】
本発明の実施例が以下で、図面を参照して説明される。
【００２４】
図１は、本発明の方法の１つの実施例を実行するための、本発明による画像処理装置の１つの実施例を示す。図１に記載される画像処理装置は、中央処理装置（ＣＰＵ）、画像処理プロセッサ１５１を有し、該プロセッサ１５１は、関心領域（例えばガン組織）を有する内臓等の関心物体を描く画像を有する多次元データセットを記憶するためのメモリ１５２に接続される。画像処理プロセッサ１５１は、複数の入出力ネットワークデバイス（ＭＲ装置又はＣＴ装置等）に接続されていてよい。画像処理プロセッサは更に、画像プロセッサ１５１において計算又は適応される情報又は画像を表示するためのコンピュータのための例えば表示装置１５４に接続されている。オペレータは、キーボード１５５及び／又は他の出力又は入力デバイス（例えば図１には記載されないコンピュータマウス）を介して画像処理プロセッサ１５１と相互作用することができる。更に、また、バスシステム１５３を介して、画像処理・制御プロセッサ１５１を、例えば、関心物体の移動をモニタするモーションモニタに接続することが可能である。例えば、患者の肺がイメージングされる場合には、移動センサは呼気センサであってもよい。心臓がイメージングされる場合には、移動センサは心電図であってもよい。
【００２５】
図２は、本発明の１つの実施例による画像レジストレーションの方法の１つの実施例のフローチャートを示す。本方法はステップＳ０で開始し、その後、多次元データセットの取得は、ステップＳ１で実行され、これは例えば多色ビームを生成する電磁気放射線の多色源によって、そして、多色ビームを検出する放射線検出器（例えばＣＴイメージングの場合）によって、実行される。
【００２６】
取得される多次元データセットは、二次元データセット又は三次元データセットの時系列であってもよいことに注意されたい。更に、データセットは、追加情報を有することができ、例えば画像取得の間に得られる心電図データによる周期的な動きに関する情報を有することができる。
【００２７】
心電図データは、（例えば患者の心臓がイメージングされる場合）心拍数に基づいて、動き前補償（pre-motion compensation）を（画像レジストレーションが実行される前に）実行するために用いられることができる。更に、例えば、患者の肺がイメージングされる場合、多次元データセットは呼気センサで測定されるデータを有することができる。
【００２８】
図２に表される方法は、いわゆる薬物動力学モデル（コンパートメントモデルである）に基づく核医学画像の時系列の局所的レジストレーションのための、モデルに基づく機構の１つの実施例である。
【００２９】
薬物動力学モデルは、管理から最終的な排出までの人間又は動物の体の中での具体的な薬物の流れを説明する。コンパートメントモデルは、この種の薬物動力学モデルの特別な数学的表現である。これは、体を通じる薬物の流れを、コンパートメントごとの物質の正味流入及び正味流出を説明する多くの交換レートによって結合される幾つかのコンパートメント又は貯蔵部間のこの物質の交換にまで、帰納的レベルで簡略化する。コンパートメント間で流速のバランスをとることは、交換レート（ｋパラメータ）の観点で各コンパートメント中の物質の時間依存含量を記述することを許可する。この種の薬物動力学モデル又はコンパートメントモデルは、よく知られた技術であるので、これ以上詳細には説明されない。
【００３０】
病巣の悪性及び治療に対するこれら病巣の応答を決定するための重要な診断ツールであるトレーサ摂取量の時間依存性を検出するために、核医学画像の時系列が取得され、病巣候補による特定の摂取値の時間変動が測定される。薬物動力学モデルは、いわゆるコンパートメント間のトレーサフローをモデル化する特性パラメータｋ１、…、ｋｎの観点で、特定の摂取値の時間依存性を説明する（コンパートメントモデル）。
【００３１】
正確な特定の摂取値を得て、特性ｋパラメータの推定値を導くためには、時間経過全体に亘って特定の摂取値が同じ解剖学的領域から得られることを保証するために、病巣候補は患者及び臓器の移動を補償するために時系列の全体を通じてトラッキングされなければならない。
【００３２】
ステップＳ２において、多次元データセットからの画像スライス又は体積データ（三次元画像）の時系列の選択が実行される。その後、ステップＳ３で、第１の関心領域（例えば病巣候補）が、データセットの第１の画像中で選択される。次に、ステップＳ４で、データセットの第２の画像中で第２の関心領域が選択される。第２の画像は、時間に関して第１の画像に続く（第１の画像が得られた後に得られることを意味する）画像であってもよい。第２の関心領域は、以下のように第１の関心領域に対応する。即ち、第２の関心領域は、同じ病巣候補であるが異なった時点における病巣候補であり、例えば臓器運動のため、空間中の異なった位置にある、病巣候補に対応する。
【００３３】
その後、ステップＳ５で、第１の変換が決定され、これは、第１の画像中の第１の関心領域の第２の画像中の第２の関心領域への変換を説明する。更に、ステップＳ６で、第１の画像中の第１の関心領域から第２の画像中の第３の関心領域への変換を説明する代替の第２の変換が決定され、第３の関心領域は、第２の関心領域に対して僅かにシフトされているものである。ステップＳ７において、（後の時点の）他の画像が、多次元データセットから選択されることができ（例えば第３の画像及び第４の画像）、病巣候補の変換は、画像ごとにトラッキングされる。これは、病巣候補の第１の画像からの、第２の画像、第３の画像及び第４の画像への移動を記述する変換ベクトルの第１の列を結果として生じ、また、第１の画像からの、第２の画像、第３の画像及び第４の画像への関心領域（病巣候補）の代替トラックを記述する、変換ベクトルの第２の列（第１の変換ベクトル列の僅かに変化したものである）を結果として生じる（ステップＳ８）。従って、４つの画像の時系列の場合には、第１の及び第２の変換ベクトル列は、それぞれ、３つの変換ベクトルを有し、各変換ベクトルは、入力画像の次元に応じて２又は３の次元を有する。一般に、ｎ個の画像の時系列は、ｎ−１個の２又は３次元変換ベクトルの変換ベクトル列を結果として生じる。
【００３４】
その後、ステップＳ９で、関心コンパートメントモデルの全てのｋパラメータを有する第１のコンパートメントベクトルＫ１（ベクトル（ｋ１，１、ｋ１，２、ｋ１，３、ｋ１，４）である）は、第１の画像の第１の関心領域から、例えば、（第２の画像の第２の関心領域及び第３の画像の第３の関心領域を介した）第４の画像の第４の関心領域への病巣の変換を記述する第１の変換ベクトル列に基づいて決定される。更に、このコンパートメントベクトルＫ１は、薬物動力学モデル（コンパートメントモデルである）に基づいて決定される。
【００３５】
更に、第２のコンパートメントベクトルＫ２＝（ｋ２，１、ｋ２，２、ｋ２，３、ｋ２，４）は、第１の画像中の第１の関心領域の第４の画像中の第５の関心領域への変換を記述する第２の変換ベクトル列に基づいて得られ、第４の画像の第５の関心領域は、第４の画像中の第４の関心領域に対して僅かにシフトされたものである（第１の変換ベクトルによって記述される）。従って、第２の変換ベクトル列は、第１の変換ベクトル列のバリエーションであり、画像の時系列を通じる病巣の第２のトラックを記述する。
【００３６】
コンパートメントモデルに基づいてコンパートメントパラメータを決定した後に、第１のコンパートメントベクトル及び第２のコンパートメントベクトルの品質決定が実行され、これは、対応する第１の品質値及び対応する第２の品質値を結果として生じる（ステップＳ１０）。
【００３７】
第１の及び第２の品質値を結果として生じる品質決定は、ｋ推定の統計品質（例えばそれらの統計分散）に基づいて、又は、イメージングされる病巣及び解剖学的組織についてのｋ値の可能な組の所与のライブラリに基づいて、実行されることができ、これは、特徴空間中におけるマッチするｋ値（ｋベクトル）の最も近い組の距離になる。特徴空間におけるこの距離は、品質値として用いられることができる。更に、本発明の１つの実施例によれば、品質決定は、画像系列を時間的に先方に（Ｉ１→Ｉ２→Ｉ３→Ｉ４）及び後方に（Ｉ４→Ｉ３→Ｉ２→Ｉ１）レジスタすることから得られるｋ推定の整合性に基づいて実行されることができる。
【００３８】
しかし、品質決定工程の他の選択肢が可能である点に留意する必要がある。
【００３９】
その後、ステップＳ１１で、第１の及び第２の品質値は互いに比較され、これらの値のうちどちらが「より良い」品質値であるかが決定される。この「より良い」品質値は、更に処理される。ステップＳ１２において、「より良い」品質値が予め設定された閾値を上回るかどうかが決定される。換言すれば、「より良い」品質値が特定の閾値基準を満たすどうかが決定される。ステップＳ１２で、「より良い」品質値が、予め設定された閾値基準（これはユーザによって手動で又は自動的にソフトウェア側からセットされることができるものである）を満たさないと決定される場合、本方法は、ステップＳ６にジャンプバックし、ここで、他の代替変換ベクトル列が決定される。この代替の他の変換ベクトル列もまた、２つの変換ベクトル列のいずれが「より良い」品質値を生じたかに応じて、第１の変換ベクトル列の他のバリエーションであるか、又は、第２の変換ベクトル列のバリエーションであってよい。
【００４０】
ステップＳ１２で、閾値基準が満たされたと決定される場合、本方法は、ステップＳ１４に進み、４つの画像のレジストレーションが実行される。レジストレーションの後、本方法はステップＳ１５で終了する。
【００４１】
本発明の１つの実施例によれば、関心領域が選択肢の既定の組から選択されることができる点に留意する必要がある。更に、コンパートメントモデルは、コンパートメントモデルの既定の組から選択されることができる。有利なことに、選択はインタラクティブに実行されることができ、データ取得の間又はこれより少し後のユーザ入力を許可する。
【００４２】
図３は、本発明の１つの実施例による第１の画像及び第２の画像並びに画像をレジスタするために実行される動作の概略図を示す。第１の画像３０１は、多次元データセットからの画像スライスであり、例えば、ＭＲＩスキャナシステム又は超音波イメージングシステムによって得られるものである。画像スライス３０１は、例えば病巣を規定する関心領域３０３を視覚化する。特定の時間後、第２の画像スライス３０２が得られ、病巣を規定する関心領域が識別される。臓器運動のため、画像スライス３０２の病巣の位置は、画像スライス３０１の病巣の位置３０３に対してシフトされうる。解剖学的コントラストがほとんど又は全くないため、関心物体（病巣）の識別は、可能ではない。従って、本発明の１つの実施例によれば、画像スライス３０２の関心物体の可能な位置の全部の組が識別され、この組は領域３０４、領域３０５及び領域３０６を有する。更に、第１の関心領域３０３から関心領域３０４への変換が識別される（画像スライス３０１を参照）。この変換は、変換３０７である。関心領域３０３の関心領域３０５への変換を指す他の変換３０８が識別され、更に、関心領域３０３の関心領域３０６への変換を指す第３の変換３０９が識別される。３つの関心領域３０４、３０５、３０６が互いに対して僅かにシフトされているので、変換３０７、３０８、３０９も同様に僅かに変化する。このとき、本発明の１つの実施例によれば、これらの変換は、３つの変換の各々についての特定の品質値を導くことによって、互いに比較される。その後、３つの変換のうち「最高」のものが、２つの画像スライス３０１、３０２のレジストレーションのために選択される。
【００４３】
本発明は、関心領域が例えばガン組織に関するものであるメディカルイメージングを参照して説明されるが、本発明はまた、例えば、実際の製品の動く画像がレジスタされるべき、材料試験又は品質管理といった非医療用アプリケーションに利用されることができることに留意されたい。
【００４４】
図４は、本発明の１つの実施例による画像の時系列をレジスタするために実行される動作の概略図を示す。
【００４５】
画像４１０は、関心領域３０３の可能な移動を表す。図４に記載される本発明の実施例によれば、関心領域運動３０３の移動についての変換ベクトルの３つの異なった軌跡は、選択されて、互いに比較される。第１の変換ベクトルは、関心領域３０３の位置３０４への変換を変換３０７によって、位置４０２への変換を変換４０７によって、そして、位置４０５への変換を変換４１２によって、記述する。第２の変換ベクトルは、関心領域３０３の、位置３０５への変換、位置４０３への変換、そして位置４０６への変換を、対応する変換３０８、４０８及び４１３によって記述する。第３の変換ベクトルは、関心領域３０３の、位置３０６への変換、位置４０１への変換、そして位置４０４への変換を、対応する変換３０９、４０９及び４１１によって記述する。位置３０４、３０５及び３０６は、２度目に取得される第２の画像スライスに関連するものであり、これは、位置３０３を有する第１の画像スライスが得られる第１回目のときよりも後であることに注意されたい。更に、位置４０１、４０２及び４０３は、第２の時間より後である第３の時間に得られる第３の画像スライスに関するものであり、位置４０４、４０５及び４０６は、第４の（最後の）時間に得られる第４の画像スライスに関するものである。
【００４６】
本発明の一側面によれば、領域３０３を領域４０５、４０４又は４０６のうちの１つにリンクする３つの変換ベクトル列の決定の後、それぞれがコンパートメントパラメータを有する対応するコンパートメントベクトルが、３つの変換ベクトル列、対応する具体的な摂取値（各画像において、各関心領域について測定される）、及び、ユーザによってインタラクティブに選択されることができる対応する薬物動力学モデル又はコンパートメントモデルに基づいて決定される。コンパートメントベクトルに基づいて、品質決定が実行され、３つの変換ベクトルのうちの１つが選択され、これは、品質決定から導かれる「最高の」品質値に関連する。本発明の１つの実施例によって、それぞれの「最高の」品質値が、予め設定された閾値を上回れば、対応する変換ベクトルは、４つの画像のレジストレーションのために用いられる。
【００４７】
換言すれば、Ｒが、疑わしい病巣を囲む関心領域を示し、Ｉ１、…、Ｉｍが、時間ｔ１、…、ｔｍに取得された核医療画像を示し、Ｔｋが、時系列中で画像Ｉｋから次の画像Ｉｋ＋１に進むときの関心領域Ｒの変換を示すこととする。変換Ｔｉ＝（Ｔｉ，１、…、Ｔｉ，ｍ−１）の全ての列について、コンパートメントパラメータｋｉ，１、…、ｋｉ，ｍの推定が得られ、これらの推定の質は、何らかの品質基準Ｑに関して決定される（ここでｉ＝１、…、ｒであり、ｒは、異なった変換の列又は変換ベクトルの数である）。時系列は、コンパートメントパラメータｋｊ，１、…、ｋｊ，ｍの推定の最高の品質を生じる変換Ｔｊのその特定の列についてのＲに関して局所的にレジスタされる。
【００４８】
品質基準Ｑの可能な選択肢は、
−ｋ推定の統計品質（例えば統計分散）
−イメージングされる病巣及び解剖学的組織についてのｋ値の可能な組のライブラリを与えられた場合には、ｋ値の最も近いマッチングセットからの特徴空間の距離、又は、
−画像系列を時間的に先方及び後方にレジスタすることにより得られるｋ推定の整合性
である。
【００４９】
本方法は、薬物動力学モデルが、用いられるトレーサ及びイメージングされる解剖学的領域で利用可能であるという条件の下で、２次元又は３次元核医療画像データセットのいかなる時系列に対しても適用されることができる。トレーサ摂取量についての正確で再現性の良い尺度を得ることは、疑わしい病巣の悪性、治療の応答性及びガンの再発の早期発見を決定するために決定的に重要である。この技術のアプリケーション領域は、分子のイメージングの分野における急速な進歩及び特異的トレーサの激増によって、急速に成長することになる。
【００５０】
有利なことに、本発明は、関心のある病巣の形状及びサイズを保存しながらも、核画像の時系列を通じた関心領域のトラッキングを可能にする。更に、トレーサ摂取量を特徴付けているパラメータについて可能な限り最高の推定を得るために、患者又は臓器移動の補償を可能にする可能性がある。従って、ほとんど又は全く解剖学的コントラストのない画像（しばしば、非常に特異的なトレーサについて起こる）が、首尾よくレジスタされることができる。提案された方法は、単に局所的に規定された関心領域を分析するだけで、増加したレジストレーション速度を可能にすることができる。
【００５１】
従って、本発明は、時系列中のコンパートメントパラメータの推定の大域的な整合を可能にし、従って、改善された画像レジストレーションを可能にする。なぜなら、トレーサ摂取の時間依存性を説明するソリッドモデルと画像の全部の組（時系列）（例えば、２つの画像及び対応するグレー値構造だけではない）が、レジストレーションのために用いられるからである。
【図面の簡単な説明】
【００５２】
【図１】本発明の方法の１つの実施例を実行するための、本発明による画像処理装置の１つの実施例を示す。
【図２】本発明による画像をレジスタする方法の１つの実施例のフローチャートを示す。
【図３】本発明の１つの実施例による画像をレジスタするために実行される画像及び動作の概略図を示す。
【図４】本発明の他の実施例による画像をレジスタするために実行される動作の概略図を示す。
【図１】
2007536054.tif 000003
【図２】
2007536054.tif 000004
【図３】
2007536054.tif 000005
【図４】
2007536054.tif 000006
【国際調査報告】
2007536054.tif 000007
2007536054.tif 000008
2007536054.tif 000009