(57)【要約】時系列における第１および第２の画像間の動き補正のための方法は、それぞれ低減された解像度の第１および第２の画像を与えるために縮小処理において第１および第２の画像を処理するステップと；それぞれ低減された解像度の第１および第２の画像における画像データのそれぞれのラプラシアンを誘導するステップを含むそれぞれ第１および第２の特徴マップを、それぞれ低減された解像度の第１および第２の画像から誘導するステップと；変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって第１および第２の特徴マップを処理することによって変位フィールドを誘導するステップと；上記変位フィールドによって第２の画像をワープするステップと、を含む。
2008511395.tif 000002
【特許請求の範囲】
【請求項１】
時系列における第１および第２の画像間の動き補正のための方法であって、前記方法は、
前記第１および第２の画像からそれぞれ第１および第２の特徴マップを誘導するステップと、
変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって、前記第１および第２の特徴マップを処理することによって変位フィールドを誘導するステップと、
前記変位フィールドによって前記第２の画像をワープするステップと、を備える動き補正のための方法。
【請求項２】
変位フィールドを誘導する前記ステップは、前記変位フィールドを誘導するために前に誘導された変位フィールドを利用する、請求項１に記載の動き補正のための方法。
【請求項３】
変位フィールドを誘導する前記ステップは、前記変位フィールドの初期誘導のために所与の初期変位フィールドを利用する、請求項２に記載の動き補正のための方法。
【請求項４】
前記所与の初期変位フィールドのためのデフォルト条件は、空集合である、請求項３に記載の動き補正のための方法。
【請求項５】
前記所与の初期変位フィールドは、患者の動きについての事前の知識を考慮している、請求項３に記載の動き補正のための方法。
【請求項６】
変位フィールドを誘導するステップを繰り返すステップを備え、各繰り返しは前の繰り返しのときより高い解像度を有する前記第１および第２の画像に対応する第１および第２の特徴マップに対して実行される、請求項２に記載の動き補正のための方法。
【請求項７】
各繰り返しは、直前のステップからの前に誘導された変位フィールドを利用する、請求項２に記載の動き補正のための方法。
【請求項８】
変位フィールドを誘導する前記ステップは、前記変位フィールドの初期誘導のために所与の変位フィールドを利用する、請求項７に記載の動き補正のための方法。
【請求項９】
変位フィールドを誘導する前記ステップは、前記前に誘導された変位フィールドを前記増加する解像度バージョンの解像度レベルに拡張するステップを備える、請求項７に記載の動き補正のための方法。
【請求項１０】
動き修正済みの画像を生成するように前記第２の画像をワープするために存在する最高の解像度で誘導された変位フィールドを利用するステップを含む、請求項７に記載の動き補正のための方法。
【請求項１１】
それぞれ第１および第２の特徴マップを誘導する前記ステップは、前記それぞれ低減された解像度の第１および第２の画像における画像データのそれぞれのラプラシアンを誘導するステップを備える、請求項７に記載の動き補正のための方法。
【請求項１２】
時系列における第１および第２の画像間の動き補正のための方法であって、
それぞれ低減された解像度の第１および第２の画像を与えるために縮小処理において前記第１および第２の画像を処理するステップと、
前記それぞれ低減された解像度の第１および第２の画像における画像データのそれぞれのラプラシアンを誘導するステップを含むそれぞれ第１および第２の特徴マップを前記それぞれ低減された解像度の第１および第２の画像から誘導するステップと、
変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって、前記第１および第２の特徴マップを処理することによって変位フィールドを誘導するステップと、
前記変位フィールドによって前記第２の画像をワープするステップと、を備える動き補正のための方法。
【請求項１３】
時系列における第１および第２の画像間の動き補正のための方法であって、前記方法は、
それぞれ低減された解像度の第１および第２の画像を与えるために縮小処理において前記第１および第２の画像を処理するステップと、
前記それぞれ低減された解像度の第１および第２の画像からそれぞれ第１および第２の特徴マップを誘導するステップと、
変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって、前記第１および第２の特徴マップを処理することによって変位フィールドを誘導するステップと、
前記変位フィールドによって前記第２の画像をワープするステップと、を備える動き補正のための方法。
【請求項１４】
それぞれ第１および第２の特徴マップを誘導する前記ステップは、前記それぞれ低減された解像度の第１および第２の画像における画像データのそれぞれのラプラシアンを誘導するステップを備える、請求項１３に記載の動き補正のための方法。
【請求項１５】
変位フィールドを誘導する前記ステップは、初期変位フィールドを利用する、請求項１３に記載の動き補正のための方法。
【請求項１６】
前記初期変位フィールドのためのデフォルト条件は、空集合である、請求項１５に記載の動き補正のための方法。
【請求項１７】
前記初期変位フィールドは、患者の動きについての事前の知識を考慮している、請求項１５に記載の動き補正のための方法。
【請求項１８】
見当合わせアルゴリズムにしたがって前記第１および第２の特徴マップを処理する前記ステップは、前記第２の画像の変位フィールドに整合する解像度レベルに前記変位フィールドを拡張するステップを備える、請求項１３に記載の動き補正のための方法。
【請求項１９】
時系列における基準画像と浮動画像との間の動き補正のための方法であって、前記方法は、
前記基準画像から漸進的に低減する解像度の第１のセットの修正された画像を誘導するステップと、
前記浮動画像から漸進的に低減する解像度の第２のセットの修正された画像を誘導するステップと、
前記第１のセットの修正された画像から第１のセットの特徴マップを誘導するステップと、
前記第２のセットの修正された画像から第２のセットの特徴マップを誘導するステップと、
変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって、それぞれ前記第１および第２のセットの特徴マップの各々の最低解像度メンバーと所与の初期変位フィールドから第１の変位フィールドを誘導するステップと、
前記見当合わせアルゴリズムにしたがって、それぞれ前記第１および第２のセットの特徴マップの各々の最低の次の解像度メンバーと前のステップで得られた前記第１の変位フィールドから第２の変位フィールドを誘導するステップと、
それぞれ前記第１および第２のセットの特徴マップの各々の、何か存在すれば、順次に高くなる解像度メンバーのために前のステップを繰り返し、各場合に前記浮動画像の解像度が到達されるまで現在のステップの前のステップで得られた変位フィールドを使用し、またより高い解像度のメンバーが存在しなければ直ぐ次のステップに進むステップと、
最後に得られた変位フィールドによって前記浮動画像をワープするステップと、を備える動き補正のための方法。
【請求項２０】
特徴マップを誘導する前記ステップは、それぞれの画像における画像データのラプラシアンを誘導するステップを備える、請求項１９に記載の動き補正のための方法。
【請求項２１】
変位フィールドを誘導する前記ステップは、拡張が必要とされれば、次に高い解像度レベルの変位フィールドに整合する解像度に前記変位フィールドを拡張するステップを備える、請求項１９に記載の動き補正のための方法。
【請求項２２】
前記所与の初期変位フィールドのためのデフォルト条件は、空集合である、請求項１９に記載の動き補正のための方法。
【請求項２３】
前記所与の初期変位フィールドは、患者の動きについての事前の知識を考慮している、請求項１９に記載の動き補正のための方法。
【請求項２４】
時系列における基準画像と浮動画像との間の動き補正のための方法であって、前記方法は、
（ａ）本明細書では、最低の解像度レベルに関してレベル０と呼ばれ、次に高い解像度レベルに関してレベル１、第２の次に高い解像度レベルに関してレベル２、そして任意の既存のより高い解像度レベルに関して同様に呼ばれるそれぞれ複数の順次に低減された解像度の基準画像と浮動画像とを与えるためにそれぞれ第１および第２の複数の縦続する解像度低減処理において前記基準画像と浮動画像を処理するステップと、
（ｂ）解像度レベルＬ０、Ｌ１、Ｌ２などの前記それぞれ複数の順次に低減された解像度の基準画像と浮動画像からそれぞれ複数の基準および浮動特徴マップを誘導するステップと、
（ｃ）第１の変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して所与の初期変位フィールドに基づいて、局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがってレベルＬ０に対応する基準特徴マップと浮動特徴マップとを処理することによって第１の変位フィールドを誘導するステップと、
（ｄ）拡張された第１の変位フィールドを与えるために解像度Ｌ１の変位フィールドに整合する解像度レベルに前記第１の変位フィールドを拡張するステップと、
（ｅ）第２の変位フィールドを形成するそれぞれのベクトルを誘導するように前記拡張された第１の変位フィールドに基づいて、前記見当合わせアルゴリズムにしたがって解像度レベルＬ１に対応する基準特徴マップと浮動特徴マップとを処理することによって第２の変位フィールドを誘導するステップと、
（ｆ）拡張された第２の変位フィールドを与えるために解像度Ｌ１の変位フィールドに整合する解像度レベルに前記第２の変位フィールドを拡張するステップと、
（ｇ）前記第２の変位フィールドを形成するそれぞれのベクトルを誘導するように前記拡張された第２の変位フィールドに基づいて、前記見当合わせアルゴリズムにしたがって解像度レベルＬ２に対応する基準特徴マップと浮動特徴マップとを処理することによって第３の変位フィールドを誘導するステップと、
（ｈ）Ｌ２が前記基準および浮動画像の解像度レベルであれば、前記第３の変位フィールドを利用することによって前記浮動画像をワープして、終了するステップと、を備える、またそうでなければ、
（ｉ）拡張された第３の変位フィールドを与えるために次に高い解像度レベルに整合する解像度レベルに前記第３の変位フィールドを拡張するステップと、
（ｊ）前記基準および浮動画像の解像度レベルが到達されるまで、最後の、前に得られて拡張された変位フィールドに基づいて段々高くなるレベルに対応する変位フィールドを取得するために解像度レベルの適当な修正をしながらステップ（ｇ）から始まる一連のステップを繰り返し、前記基準および浮動画像の解像度レベルが到達されたとき、得られた最後の変位フィールドによって前記浮動画像をワープして、終了するステップと、を備える動き補正のための方法。
【請求項２５】
それぞれ複数の順次に低減された解像度の基準および浮動画像を与える前記ステップは、前記それぞれ複数の順次に低減された解像度の基準および浮動画像における画像データのそれぞれのラプラシアンを誘導するステップを備える、請求項２４に記載の動き補正のための方法。
【請求項２６】
第１の変位フィールドを誘導する前記ステップは、所与の初期変位フィールドを入力するステップを備える、請求項２４に記載の動き補正のための方法。
【請求項２７】
時系列における第１および第２の画像間の動き補正のための方法であって、前記方法は、
それぞれ低減された解像度の第１および第２の画像からそれぞれ第１および第２の特徴マップを誘導するステップと、
変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって、前記第１および第２の特徴マップを処理することによって第１の変位フィールドを誘導するステップと、
前記第１および第２の画像の解像度に対応する前記第１の変位フィールドを拡張するステップと、
前記第１の変位フィールドに基づいて、前記見当合わせアルゴリズムにしたがって前記第１および第２の画像を処理することによって第２の変位フィールドを誘導するステップと、
動き修正済み画像を得るように前記第２の変位フィールドによって前記第２の画像をワープするステップと、を備える動き補正のための方法。
【請求項２８】
第１の変位フィールドを誘導する前記ステップは、所与の初期変位フィールドを入力するステップを備える、請求項２７に記載の動き補正のための方法。
【請求項２９】
前記所与の初期変位フィールドのためのデフォルト条件は、空集合である、請求項２８に記載の動き補正のための方法。
【請求項３０】
前記所与の初期変位フィールドは、患者の動きについての事前の知識を考慮している、請求項２８に記載の動き補正のための方法。
【請求項３１】
それぞれ第１および第２の特徴マップを誘導する前記ステップは、前記それぞれ低減された解像度の第１および第２の画像における画像データのそれぞれのラプラシアンを誘導するステップを備える、請求項２７に記載の動き補正のための方法。
【請求項３２】
時系列における第１および第２の画像間の動き補正のための方法であって、前記方法は、
漸進的に低減された解像度を有する第１の系列の画像を前記第１の画像から誘導するステップと、
漸進的に低減された解像度を有する第２の系列の画像を前記第２の画像から誘導するステップと、
前記第１の系列の画像から第１の系列の特徴マップを誘導するステップと、
前記第２の系列の画像から第２の系列の特徴マップを誘導するステップと、
変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える、それぞれの変位フィールドを与えるための見当合わせアルゴリズムにしたがって、同じ解像度に関する、前記第２の系列の特徴マップの対応する特徴マップによって前記第１の系列の特徴マップの特徴マップを処理することによって１系列の変位フィールドを誘導し、最低の解像度レベルにおける第１の変位フィールドを誘導するために使用される所与の変位フィールドから始め、その後に前記系列における次に高い解像度に対応する次の変位フィールドを誘導するために直前の変位フィールドを利用するステップと、
前記第１および第２の画像の解像度レベルにおいて誘導された最後の変位フィールドを利用して前記第２の画像をワープするステップと、を備える動き補正のための方法。
【請求項３３】
前記所与のフィールドと次に高い解像度に対応する前記最後のフィールドとを除いて前記各解像度のフィールドを拡張するステップを備える、請求項３２に記載の動き補正のための方法。
【請求項３４】
時系列における基準画像と浮動画像との間の動き補正のための方法であって、前記方法は、
（ａ）本明細書では、最低の解像度レベルに関してレベル０と呼ばれ、次に高い解像度レベルに関してレベル１、第２の次に高い解像度レベルに関してレベル２、そして任意の既存のより高い解像度レベルに関して同様に呼ばれるそれぞれ複数の順次に低減された解像度の基準画像と浮動画像を与えるためにそれぞれ第１および第２の複数の縦続する解像度低減処理において前記基準画像と浮動画像とを処理するステップと、
（ｂ）前記基準および浮動画像の解像度レベルであるＬ２によって増加する解像度の順に、解像度レベルＬ０、Ｌ１、Ｌ２における前記それぞれ複数の順次に低減された解像度の基準画像と浮動画像からそれぞれ複数の基準および浮動特徴マップを誘導するステップと、
（ｃ）第１の変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して所与の初期変位フィールドに基づいて、局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがってレベルＬ０に対応する基準特徴マップと浮動特徴マップとを処理することによって第１の変位フィールドを誘導するステップと、
（ｄ）拡張された第１の変位フィールドを与えるために解像度Ｌ１の変位フィールドに整合する解像度レベルに前記第１の変位フィールドを拡張するステップと、
（ｅ）第２の変位フィールドを形成するそれぞれのベクトルを誘導するように前記拡張された第１の変位フィールドに基づいて、前記見当合わせアルゴリズムにしたがって解像度レベルＬ１に対応する基準特徴マップと浮動特徴マップとを処理することによって第２の変位フィールドを誘導するステップと、
（ｆ）拡張された第２の変位フィールドを与えるために解像度Ｌ２の変位フィールドに整合する解像度レベルに前記第２の変位フィールドを拡張するステップと、
（ｇ）最後の変位フィールドを形成するそれぞれのベクトルを誘導するように前記拡張された第２の変位フィールドに基づいて、前記見当合わせアルゴリズムにしたがって解像度レベルＬ２に対応する基準特徴マップと浮動特徴マップとを処理することによって最後の変位フィールドを誘導するステップと、
（ｈ）前記最後の変位フィールドを利用することによって前記浮動画像をワープするステップと、を備える動き補正のための方法。
【請求項３５】
前記所与の初期変位フィールドのためのデフォルト条件は、空集合である、請求項３４に記載の動き補正のための方法。
【請求項３６】
前記所与の初期変位フィールドは、患者の動きについての事前の知識を考慮している、請求項３４に記載の動き補正のための方法。
【請求項３７】
特徴マップを誘導する前記ステップは、前記それぞれ低減された解像度の第１および第２の画像における画像データのそれぞれのラプラシアンを誘導するステップを備える、請求項３４に記載の動き補正のための方法。
【請求項３８】
第１の変位フィールドを誘導するステップ（ｃ）は、
その座標がベクトルｘj∈Ｒ3によって与えられるインデックスｊによって各ボクセルにラベル付けするステップと、
前記アルゴリズムのｋ番目の繰り返しの後に点ｘjにおいて回収される変位ベクトルとしてｐjk∈Ｒ3を定義するステップと、
１セットの変位ベクトルｐjkが与えられたとして、Ｊ2,k：
【数１】
2008511395.tif 000023
によって繰り返しｋにおける前記浮動特徴マップの修正済みバージョンを示すステップと、
各ボクセルｊにおいてｓjkによる更新ルールｐjk+1＝ｐjk＋ｓjkを非線形最小二乗平均問題の解：
【数２】
2008511395.tif 000024
（ｚik（ｓ）＝Ｊ1（ｘi）−Ｊ2,k（ｘi＋ｓ）であり、Ｇσ（ｘ）は標準偏差σの等方性３次元ガウシアン・カーネル：
【数３】
2008511395.tif 000025
である）として定義するステップと、
下記の基準：
【数４】
2008511395.tif 000026
（
【数５】
2008511395.tif 000027
は
【数６】
2008511395.tif 000028
として定義される）を最小化することによって前記非線形最小二乗
平均問題の線形化バージョンを解くステップと、
最適性の必要条件：
【数７】
2008511395.tif 000029
が閉じた形式の解ｓjk：
【数８】
2008511395.tif 000030
を与えるようにＥ（ｓjk）の第１の変動を計算するステップと、を含む、請求項３４に記載の動き補正のための方法。
【請求項３９】
プログラムと他のデータとを記憶するためのメモリ装置と、
前記メモリ装置と通信し、また時系列における第１および第２の画像間の動き補正のための方法を実行するために前記プログラムによって動作するプロセッサと、を備えた、画像の動き補正を実行するためのシステムであって、前記方法は、
それぞれ低減された解像度の第１および第２の画像を与えるために縮小処理において前記第１および第２の画像を処理するステップと、
前記それぞれ低減された解像度の第１および第２の画像からそれぞれ第１および第２の特徴マップを誘導するステップと、
変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって、前記第１および第２の特徴マップを処理することによって変位フィールドを誘導するステップと、
前記変位フィールドによって前記第２の画像をワープするステップと、を備える、画像の動き補正を実行するためのシステム。
【請求項４０】
時系列における第１および第２の画像間の動き補正のための方法によって画像の動き補正を実行するためのプログラムコードのためにコンピュータプログラム論理を媒体上に記録したコンピュータ使用可能媒体を備えるコンピュータプログラム製品であって、前記方法は、
それぞれ低減された解像度の第１および第２の画像を与えるために縮小処理において前記第１および第２の画像を処理するステップと、
前記それぞれ低減された解像度の第１および第２の画像からそれぞれ第１および第２の特徴マップを誘導するステップと、
変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって、前記第１および第２の特徴マップを処理することによって変位フィールドを誘導するステップと、
前記変位フィールドによって前記第２の画像をワープするステップと、を備える、コンピュータプログラム製品。
【請求項４１】
時系列における第１および第２の画像間の動き補正のための方法であって、前記方法は、
変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムを含む処理にしたがって前記第１および第２の画像を処理することによって変位フィールドを誘導するステップと、
前記変位フィールドによって前記第２の画像をワープするステップと、を備える動き補正のための方法。
【請求項４２】
処理にしたがって前記第１および第２の画像を処理する前記ステップは、前記第１および第２の画像からそれぞれ第１および第２の特徴マップを誘導して、前記特徴マップにそれぞれ前記見当合わせアルゴリズムを適用するステップを含む、請求項４１に記載の動き補正のための方法。
【請求項４３】
変位フィールドを誘導する前記ステップは、前記変位フィールドを誘導するために前に誘導された変位フィールドを利用する、請求項４１に記載の動き補正のための方法。
【請求項４４】
変位フィールドを誘導する前記ステップは、前記変位フィールドの初期誘導のために所与の初期変位フィールドを利用する、請求項４１に記載の動き補正のための方法。
【発明の詳細な説明】
【技術分野】
【０００１】
関連出願の相互参照
優先権
ここにおいて、本出願の発明者であるＣＨＲＩＳＴＯＰＨＥ  ＣＨＥＦＤ’ＨＯＴＥＬの名前で２００４年８月３１に出願された、「Ａ  ＧＡＵＳＳＩＡＮ  ＷＥＩＧＨＴＥＤ  ＬＥＡＳＴ  ＭＥＡＮ  ＳＱＵＡＲＥ  ＲＥＧＩＳＴＲＡＴＩＯＮ  ＡＬＧＯＲＩＴＨＭ  ＦＯＲ  ＢＲＥＡＳＴ  ＭＲ  ＭＯＴＩＯＮ  ＣＯＲＲＥＣＴＩＯＮ」と題する米国仮特許出願第６０／６０５，７５９号、代理人整理番号第２００４Ｐ１４８２２ＵＳ号に対して特定の参照が行われ、その優先権の恩典が請求されており、その開示内容は参照により本明細書に組み込まれる。
【０００２】
本出願は、一般的には画像形成における動き修正に関し、また特定的には例えば医療画像形成プロセスによって取得され得るような時系列的な一連の画像において発生し得るような被験者の動きの修正に関する動き補正に関する。
【背景技術】
【０００３】
医療画像形成技法は、例えば患者の癌または前癌状態の検出時を含む多くの医療処置において使用される。重要な用途は、乳癌における腫瘍または潜在腫瘍の検出にある。潜在腫瘍は、検出が困難である。潜在的に有用な情報を与える利用可能な技法の間で、例えばこのような腫瘍関連組織が典型的に、造影剤のより迅速な取り入れ（ウォッシュイン）ならびに隣接する非腫瘍組織よりも迅速なウォッシュアウト（洗い流し）を示し得ることが知られている。これらおよびその他のような特性は、造影剤のウォッシュインおよび／またはウォッシュアウトのような処置の前後で行われる患者の画像の比較による懸念組織の検出と組織特性の識別を含むある幾つかの診断に有用であり得る。磁気共鳴映像法（ＭＲＩ）のような画像形成技法によって作られたこのような時系列的画像を使用して、取得されたＭＲボリュームの異なる部位によって示される挙動の差異を検出するために、画像間の比較が行われ得る。
【０００４】
有利にはこの検出を実行するための技法は、このようなボリュームの時系列における単一ボクセルの強度を追跡するための技法を必要とする。しかしながら困難は、患者が典型的には連続的な画像取得の間で動き、それによって、取得された画像間での動き関連差異を導入し、それによって空間内の単一の点は、動き修正が実施されなければ、もはや追跡され得ないという点において発生する。本明細書で使用されるように空間内の点は、古典的な幾何学的に定義された寸法を持たない点を意味することを意図しておらず、むしろディジタル化された画像を構成するようになる要素の小さな寸法を有するディジタル化手順の結果から得られる点を意味する。
【０００５】
過去におけるこの問題を解決するための従来技術のアプローチは、任意の一方の画像がその系列の画像の間で基準として選択される二つの画像間における光学的流れを計算している。例えばこれら二つの画像は、ラプラシアン・ピラミッドを計算することによって取得された画像から取得できる。光学的流れは、例えばこれら二つのラプラシアン画像の間のポイントツーポイント差の最小化問題を解くことによって計算され得る。
【発明の開示】
【０００６】
例えば造影剤の迅速な取り入れ（ウォッシュイン）ならびに迅速なウォッシュアウト（洗い流し）を有する組織として検出される潜在的腫瘍の胸部ＭＲ検出において有利な仕方で動き修正または補正問題を解決することは本発明の目的である。
【０００７】
本発明の一態様によれば患者の第１および第２の画像は、ＭＲＩ装置のような画像形成プロセスから時系列的に取得され、この場合第２の遅い方の画像は第１の早い方の画像が取られた後の例えば患者の動きに起因する第１の早い方の画像からの差異を含む可能性がある。これらの画像は、縮小処理において処理され、その結果、元の画像と比較して低い解像度のそれぞれ第１および第２の画像を生じる。低い解像度の第１および第２の画像における画像データのそれぞれのラプラシアンを誘導することによって、低い解像度の第１および第２の画像からそれぞれ第１および第２の特徴マップが誘導される。所与の初期変位フィールドに基づいてこれらの特徴マップは、見当合わせアルゴリズムにしたがって見当合わせ処理において処理される。このアルゴリズムは、歪みをモデル化するために高濃度の変位フィールドを形成するそれぞれのベクトルを誘導するように各ボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える。この変位フィールドは、動き修正された第２の画像を取得するように第２の画像をワープするために利用される。
【０００８】
上記の所与の初期変位フィールドのためのデフォルト条件は、例えば空集合またはゼロ変位であるか、あるいは上記の所与の初期変位フィールドは、患者の動きについての事前の知識を考慮している可能性がある。
【０００９】
本発明のもう一つの態様によれば、時系列における第１および第２の画像間の動き補正のための方法は、第１および第２の画像からそれぞれ第１および第２の特徴マップを誘導するステップと；変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって第１および第２の特徴マップを処理することによって変位フィールドを誘導するステップと；この変位フィールドによって第２の画像をワープするステップと、を備える。
【００１０】
本発明のもう一つの態様によれば、変位フィールドを誘導するステップは、変位フィールドを誘導するために前に誘導された変位フィールドを利用し、また変位フィールドを誘導するステップは、変位フィールドの初期誘導のために所与の初期変位フィールドを利用する。
【００１１】
所与の初期変位フィールドのためのデフォルト条件は空集合であるか、あるいは例えば所与の初期変位フィールドは、患者の動きについての事前の知識を考慮している可能性がある。
【００１２】
本発明のもう一つの態様によれば、変位フィールドを誘導するステップは繰り返えされ、各繰り返しは前の繰り返しのときより高い解像度を有する第１および第２の画像に対応する第１および第２の特徴マップに対して実行される。
【００１３】
本発明のもう一つの態様によれば各繰り返しは、直前のステップからの前に誘導された変位フィールドを利用する。
【００１４】
本発明のもう一つの態様によれば変位フィールドを誘導するステップは、変位フィールドの初期誘導のために所与の変位フィールドを利用する。
【００１５】
本発明のもう一つの態様によれば変位フィールドを誘導するステップは、前に誘導された変位フィールドを増加する解像度バージョンの解像度レベルに拡張するステップを備える。
【００１６】
本発明のもう一つの態様によれば動き補正のための方法は、動き修正済みの画像を生成するように第２の画像をワープするために存在する最高の解像度で誘導された変位フィールドを利用するステップを含む。
【００１７】
本発明のもう一つの態様によればそれぞれ第１および第２の特徴マップを誘導するステップは、それぞれ低減された解像度の第１および第２の画像における画像データのそれぞれのラプラシアンを誘導するステップを備える。
【００１８】
本発明のもう一つの態様によれば、時系列における第１および第２の画像間の動き補正のための方法は、それぞれ低減された解像度の第１および第２の画像を与えるために縮小処理において第１および第２の画像を処理するステップと；それぞれ低減された解像度の第１および第２の画像から、それぞれ低減された解像度の第１および第２の画像における画像データのそれぞれのラプラシアンを誘導するステップを備えたそれぞれ第１および第２の特徴マップを誘導するステップと；変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって第１および第２の特徴マップを処理することによって変位フィールドを誘導するステップと；この変位フィールドによって第２の画像をワープするステップと、を備える。
【００１９】
本発明のもう一つの態様によれば、時系列における第１および第２の画像間の動き補正のための方法は、それぞれ低減された解像度の第１および第２の画像を与えるために縮小処理において第１および第２の画像を処理するステップと；それぞれ低減された解像度の第１および第２の画像からそれぞれ第１および第２の特徴マップを誘導するステップと；変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって第１および第２の特徴マップを処理することによって変位フィールドを誘導するステップと；この変位フィールドによって第２の画像をワープするステップと、を備える。
【００２０】
本発明のもう一つの態様によれば、時系列における基準画像と浮動画像との間の動き補正のための方法は、基準画像から漸進的に低減する解像度の第１のセットの修正された画像を誘導するステップと；浮動画像から漸進的に低減する解像度の第２のセットの修正された画像を誘導するステップと；第１のセットの修正された画像から第１のセットの特徴マップを誘導するステップと；第２のセットの修正された画像から第２のセットの特徴マップを誘導するステップと；変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって、それぞれ第１および第２のセットの特徴マップの各々の最低の解像度のメンバーと所与の初期変位フィールドから第１の変位フィールドを誘導するステップと；上記の見当合わせアルゴリズムにしたがって、それぞれ第１および第２のセットの特徴マップの各々の最低の次の解像度のメンバーと前のステップで得られた第１の変位フィールドから第２の変位フィールドを誘導するステップと；それぞれ第１および第２のセットの特徴マップの各々の、何か存在すれば、順次に高い解像度のメンバーのために前のステップを繰り返し、各場合に上記浮動画像の解像度が到達されるまで現在のステップの前のステップで得られた変位フィールドを使用し、またより高い解像度のメンバーが存在しなければ直ぐ次のステップに進むステップと；最後に得られた変位フィールドによって上記浮動画像をワープするステップと、を備える。
【００２１】
本発明のもう一つの態様によれば、時系列における基準画像と浮動画像との間の動き補正のための方法は、（ａ）本明細書では、最低の解像度レベルに関してレベル０と呼ばれ、次に高い解像度レベルに関してレベル１、第２の次に高い解像度レベルに関してレベル２、そして任意の既存のより高い解像度レベルに関して同様に呼ばれるそれぞれ複数の順次に低減された解像度の基準画像と浮動画像を与えるためにそれぞれ第１および第２の複数の縦続する解像度低減処理において基準画像と浮動画像を処理するステップと；（ｂ）解像度レベルＬ０、Ｌ１、Ｌ２などのそれぞれ複数の順次に低減された解像度の基準画像と浮動画像からそれぞれ複数の基準および浮動特徴マップを誘導するステップと；（ｃ）第１の変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して所与の初期変位フィールドに基づいて、局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがってレベルＬ０に対応する基準特徴マップと浮動特徴マップとを処理することによって第１の変位フィールドを誘導するステップと；（ｄ）拡張された第１の変位フィールドを与えるために解像度Ｌ１の変位フィールドに整合する解像度レベルに上記第１の変位フィールドを拡張するステップと；（ｅ）第２の変位フィールドを形成するそれぞれのベクトルを誘導するように上記拡張された第１の変位フィールドに基づいて上記見当合わせアルゴリズムにしたがって解像度レベルＬ１に対応する基準特徴マップと浮動特徴マップとを処理することによって第２の変位フィールドを誘導するステップと；（ｆ）拡張された第２の変位フィールドを与えるために解像度Ｌ１の変位フィールドに整合する解像度レベルに第２の変位フィールドを拡張するステップと；（ｇ）第２の変位フィールドを形成するそれぞれのベクトルを誘導するように上記拡張された第２の変位フィールドに基づいて上記見当合わせアルゴリズムにしたがって解像度レベルＬ２に対応する基準特徴マップと浮動特徴マップとを処理することによって第３の変位フィールドを誘導するステップと；（ｈ）Ｌ２が上記基準および浮動画像の解像度レベルであれば、第３の変位フィールドを利用することによって上記浮動画像をワープして、終了するステップと、を備える、またそうでなければ、（ｉ）拡張された第３の変位フィールドを与えるために次に高い解像度レベルに整合する解像度レベルに第３の変位フィールドを拡張するステップと、（ｊ）上記基準および浮動画像の解像度レベルが到達されるまで、最後の、前に得られて拡張された変位フィールドに基づいて段々に高くなる解像度レベルに対応する変位フィールドを取得するために解像度レベルの適当な修正をしながらステップ（ｇ）から始まる一連のステップを繰り返し、上記基準および浮動画像の解像度レベルが到達されたとき、得られた最後の変位フィールドによって上記浮動画像をワープして、終了するステップと、を備える。
【００２２】
本発明のもう一つの態様によれば、それぞれ複数の順次に低減された解像度の基準および浮動画像を与えるステップは、それぞれ複数の順次に低減された解像度の基準および浮動画像における画像データのそれぞれのラプラシアンを誘導するステップを備える。
【００２３】
本発明のもう一つの態様によれば、第１の変位フィールドを誘導するステップは、所与の初期変位フィールドを入力するステップを備える。
【００２４】
本発明のもう一つの態様によれば、時系列における第１および第２の画像間の動き補正のための方法は、それぞれ低減された解像度の第１および第２の画像からそれぞれ第１および第２の特徴マップを誘導するステップと；変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって、第１および第２の特徴マップを処理することによって第１の変位フィールドを誘導するステップと；第１および第２の画像の解像度に対応するように第１の変位フィールドを拡張するステップと；第１の変位フィールドに基づいて上記見当合わせアルゴリズムにしたがって、第１および第２の画像を処理することによって第２の変位フィールドを誘導するステップと；動き修正済みの画像を得るように第２の変位フィールドによって第２の画像をワープするステップと、を備える。
【００２５】
本発明のもう一つの態様によれば、時系列における第１および第２の画像間の動き補正のための方法は、漸進的に低減された解像度を有する第１の系列の画像を第１の画像から誘導するステップと；漸進的に低減された解像度を有する第２の系列の画像を第２の画像から誘導するステップと；第１の系列の画像から第１の系列の特徴マップを誘導するステップと；第２の系列の画像から第２の系列の特徴マップを誘導するステップと；変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える、それぞれの変位フィールドを与えるための見当合わせアルゴリズムにしたがって、同じ解像度に関する、第２の系列の特徴マップの対応する特徴マップによって第１の系列の特徴マップの特徴マップを処理することによって１系列の変位フィールドを誘導し、最低の解像度レベルにおける第１の変位フィールドを誘導するために使用される所与の変位フィールドから始め、その後に上記系列における次に高い解像度に対応する次の変位フィールドを誘導するために直前の変位フィールドを利用するステップと；上記第１および第２の画像の解像度レベルにおいて最後の変位フィールドを利用して第２の画像をワープするステップと、を備える。
【００２６】
本発明のもう一つの態様によれば、時系列における基準画像と浮動画像との間の動き補正のための方法は、（ａ）本明細書では、最低の解像度レベルに関してレベル０と呼ばれ、次に高い解像度レベルに関してレベル１、第２の次に高いレベルに関してレベル２、そして任意の既存のより高い解像度レベルに関して同様に呼ばれるそれぞれ複数の順次に低減された解像度の基準画像と浮動画像とを与えるためにそれぞれ第１および第２の複数の縦続する解像度低減処理において基準画像と浮動画像とを処理するステップと；（ｂ）上記基準および浮動画像の解像度レベルであるＬ２によって解像度の増加する順に、解像度レベルＬ０、Ｌ１、Ｌ２における上記それぞれ複数の順次に低減された解像度の基準画像と浮動画像からそれぞれ複数の基準および浮動特徴マップを誘導するステップと；（ｃ）第１の変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して所与の初期変位フィールドに基づいて、局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがってレベルＬ０に対応する基準特徴マップと浮動特徴マップとを処理することによって第１の変位フィールドを誘導するステップと；（ｄ）拡張された第１の変位フィールドを与えるために解像度Ｌ１の変位フィールドに整合する解像度レベルに第１の変位フィールドを拡張するステップと；（ｅ）第２の変位フィールドを形成するそれぞれのベクトルを誘導するように上記拡張された第１の変位フィールドに基づいて上記見当合わせアルゴリズムにしたがって解像度レベルＬ１に対応する基準特徴マップと浮動特徴マップとを処理することによって第２の変位フィールドを誘導するステップと；（ｆ）拡張された第２の変位フィールドを与えるために解像度Ｌ２の変位フィールドに整合する解像度レベルに第２の変位フィールドを拡張するステップと；（ｇ）最後の変位フィールドを形成するそれぞれのベクトルを誘導するように上記拡張された第２の変位フィールドに基づいて上記見当合わせアルゴリズムにしたがって解像度レベルＬ２に対応する基準特徴マップと浮動特徴マップとを処理することによって最後の変位フィールドを誘導するステップと；（ｈ）上記最後の変位フィールドを利用することによって浮動画像をワープするステップと、を備える。
【００２７】
本発明のもう一つの態様によれば、画像の動き補正を実行するためのシステムは、プログラムと他のデータとを記憶するためのメモリ装置と；メモリ装置と通信し、また時系列における第１および第２の画像間の動き補正のための方法を実行するために上記プログラムによって動作するプロセッサと、を備え、上記方法は、それぞれ低減された解像度の第１および第２の画像を与えるために縮小処理において第１および第２の画像を処理するステップと；それぞれ低減された解像度の第１および第２の画像からそれぞれ第１および第２の特徴マップを誘導するステップと；変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均
問題を解くステップを備える見当合わせアルゴリズムにしたがって第１および第２の特徴マップを処理することによって変位フィールドを誘導するステップと；上記変位フィールドによって第２の画像をワープするステップと、を備える。
【００２８】
本発明のもう一つの態様によれば、コンピュータプログラム製品は、時系列における第１および第２の画像間の動き補正のための方法であって、それぞれ低減された解像度の第１および第２の画像を与えるために縮小処理において第１および第２の画像を処理するステップと；それぞれ低減された解像度の第１および第２の画像からそれぞれ第１および第２の特徴マップを誘導するステップと；変位フィールドを形成するそれぞれのベクトルを誘導するように各画素またはボクセルに関して局所的ガウス重み付け最小二乗平均問題を解くステップを備える見当合わせアルゴリズムにしたがって第１および第２の特徴マップを処理することによって変位フィールドを誘導するステップと；上記変位フィールドによって第２の画像をワープするステップと、を備える方法によって画像の動き補正を実行するためのプログラムコードのためにコンピュータプログラム論理を媒体上に記録したコンピュータ使用可能媒体を備える。
【００２９】
本発明は、図面と関連付けて上記の例示的実施形態を含む詳細な説明から更に十分に理解される。
【発明を実施するための最良の形態】
【００３０】
本発明の一実施形態によれば、患者の基準画像は、ＭＲＩ装置を利用するような医療画像形成手順によって取得され、また浮動画像は典型的には、この基準画像より後の時刻に取得される。浮動画像は、例えば第１の画像が取られた後の患者の動きによる基準画像との差を含み得る。これらの画像は、元のそれぞれの画像と比較して低い解像度のそれぞれの基準画像と浮動画像という結果をもたらす縮小処理において処理される。それぞれの基準画像および浮動画像の特徴マップは、より低い解像度の基準および浮動画像における画像データのそれぞれのラプラシアンを誘導することによって低い解像度の基準および浮動画像から誘導される。所与の初期変位フィールドに基づいて特徴マップは、見当合わせ（レジストレーション）アルゴリズムにしたがって見当合わせ手順において処理される。このアルゴリズムは、歪みをモデル化するための高濃度の変位フィールドを形成するそれぞれのベクトルを誘導するように局所的ガウス重み付け最小二乗平均問題を各ボクセルに関して解くステップを備える。この変位フィールドは、動き修正された第２の画像を得るように浮動画像をワープするために利用される。上記のように所与の初期変位フィールドに関する可能なデフォルト条件は、空集合であるか、あるいは例えば所与の初期変位フィールドは、患者の動きについての事前知識を考慮している可能性がある。
【００３１】
胸部ＭＲ動き修正のためのガウス重み付け最小二乗平均見当合わせアルゴリズムは、腫瘍関連組織の検出時に動き修正を与えるために適用され、また胸部ＭＲ動き修正のための柔軟な見当合わせアルゴリズムは、腫瘍関連組織の検出時に動き修正を与えるために適用される。歪みは、今後説明されるように一連の局所的ガウス重み付け最小二乗平均問題を解くことによって得られるベクトルで形成される高濃度変位フィールドとしてモデル化される。
【００３２】
この見当合わせアルゴリズムは、体積測定的画像間で高濃度変位フィールドを推定するために利用される。目的は、磁気共鳴（ＭＲ）画像の時系列における動きの人為構成物を補正することである。典型的には動きの人為構成物は主として、呼吸、心臓拍動および患者の動きに起因する。本発明に関する主要な応用分野は、胸部ＭＲ研究における強度変化の研究である。胸部組織における造影剤取り入れ量の時間的変化は、病変検出時に有用な情報を与えるために使用され得る。このような変化を計算することは画像取得間における胸部組織の正確な空間的整合を必要とする。
【００３３】
本明細書では、高速フィルタリング技法を使用して実現されるガウス重み付けの導入が、少なくとも二つの利点を提供すること、すなわち：
変位フィールドの平滑さを保証し、それによって特異性を避けて無条件的歪みモデルの規則性を与えるための簡単で計算的に効率的な方法を提供することと、
最小二乗平均問題が適当に公式化されることを保証することと、を提供することが認められる。
【００３４】
更に本発明の原理による画像類似性の尺度は、画像強度を使用することよりむしろ画像ラプラシアンの計算に依存している。本発明は、このアプローチが造影剤取り入れに起因する強度変化に対処する際に有効であることを見出している。このアプローチの捕捉範囲と速度と強さとを改善するために、多解像度戦略が使用される。
【００３５】
本発明は、例示的な例によって更に詳細に説明される。本発明の本例示的例の目的のために見当合わせされる１対の画像の強度は二つの実数値関数Ｉ1：
【数９】
2008511395.tif 000003
およびＩ2：
【数１０】
2008511395.tif 000004
、Ω⊂Ｒ3の規則的な格子上でサンプリングされたバージョン
であることが仮定されており、ここでこれらの記号は通常の意味を有し、したがって本例ではＩ2は入力Ωと出力Ｒのマップレットを備える関数であり、ΩはＲ3の適当なサブセットであり、以下同様である。本発明者らはここで体積測定的画像の場合を扱っているが、引き続いて生じる展開は任意の次元に直ちに拡張することに留意すること。
【００３６】
下記において本発明者らは、基準画像としてＩ1を指し、浮動画像としてＩ2を指す。本発明者らは、Ｉ1と浮動画像Ｉ2ο（ｉｄ＋ｐ）（ここでｉｄはアイデンティティ・マップであり、οはコンポジション演算子である）のワープバージョンとの間の類似性を最大にする変位フィールドｐ：
【数１１】
2008511395.tif 000005
（格子点において推定される）を見出すステップとして見当合わせ問題を定義する。
【００３７】
画像強度は、必ずしも見当合わせの仕事を実行するための最善の特徴ではない。この見通しから、他の画像の特徴（例えば画像ラプラシアンまたは画像階調度のノルム）を表す、Ｉ1およびＩ2から計算される補助関数Ｊ1：
【数１２】
2008511395.tif 000006
およびＪ2：
【数１３】
2008511395.tif 000007
を研究対象にすることが時には有用である。特に胸部ＭＲへの応用においては
、画像ラプラシアンを計算することは、実際上極めて有効となっている。
【００３８】
このアルゴリズムは反復性である。これは、「真の」動きフィールドに向かって収束する一連の変位フィールドｐ0、ｐ1、・・・、ｐkを構築するステップを備える。最後の「動き修正済み」画像は、Ｉ2ο（ｉｄ＋ｐk）によって与えられる。この見当合わせ技法の核心を形成する変位フィールドの反復構造は、下記に説明される。
【００３９】
反復手順は、粗大から微細への多解像度戦略と組み合わせて適用される。このアプローチの利点は、
局所的最小値にトラップされる危険を減らすことによって捕捉範囲が改善されることと、
より低い解像度でより少ないデータを処理することによって計算コストが減らされることである。
【００４０】
この多解像度戦略は、下記のような例によって説明できる。
平均化とサブサンプリングとによって基準画像と浮動画像とに関する二つの多解像度ピラミッドを構築する。このタイプのピラミッド構造は、画像処理の標準的技法であって、詳細については例えば２．Ｒ．ＧｏｎｚａｌｅｓとＲ．Ｗｏｏｄｓによるテキストブック、Ｄｉｇｉｔａｌ  Ｉｍａｇｅ  Ｐｒｏｃｅｓｓｉｎｇ、２ｎｄ  Ｅｄｉｔｉｏｎ、Ｐｒｅｎｔｉｃｅ−Ｈａｌｌ，２００２が参照される。
Ｉ1lおよびＩ2lを結果として得られる解像度ｌの画像とする。
Ｊ1lおよびＪ2lとラベル付けされた対応する特徴画像を計算する。
ｐk,lを解像度ｌでｋ番目の反復の後に得られる変位フィールドとする。ｐ0,0は最低解像度における初期変位フィールド（一般に空フィールド）を示す。以下更に説明されるように、解像度ｌで一定回数の反復のために見当合わせアルゴリズムを実行する。トライリニアー（ｔｒｉ−ｌｉｎｅａｒ）内挿法と拡大縮小とによって次の解像度に変位フィールドを拡張する。この結果を初期変位フィールドｐ0,l+1として使用する。
元の解像度が到達されるまで上記の演算を実行する。
【００４１】
これらのステップは、図１に概略図的形式で例示されている。図１を参照すると解像度は、「ｌ」として示されている。この図に描かれている３レベルは、例示的実施形態によるもので、例示目的のためであり、本発明の範囲内でより多い、またはより少ない拡張ステップが考えられる。図１において、解像度レベルｌ＝ｌmax＝２においてそれぞれｌ1およびｌ2で示される基準画像と浮動画像は、第１の縮小処理を受け、その結果として解像度ｌ＝１のそれぞれの画像Ｉ11およびＩ21を与え、続いて更なる縮小処理を受け、その結果得られる画像は、それぞれ解像度レベルｌ＝０のＩ10、Ｉ20となる。特徴マップＪ10およびＪ20は、それぞれ画像Ｉ10およびＩ20から誘導される。特徴マップＪ10およびＪ20は、所与の初期変位フィールドｐ0,0を使用する、本明細書に説明されたような見当合わせ処理を受ける。この初期変位フィールドは、任意のものであってよく、ゼロであってもよい。この見当合わせ処理は、初期変位フィールドと比較して改善された変位フィールドと考えられ得る変位フィールドｐk,0という結果をもたらす。それから変位フィールドｐk,0は、画像解像度ｌ＝１に対応する拡張あるいは膨張を受け、その結果、拡張された解像度フィールドｐ0,1となる。
【００４２】
特徴マップＪ11およびＪ21は、それぞれ画像Ｉ11およびＩ21から誘導される。特徴マップＪ11およびＪ21は、ｐ0,2に拡張される変位フィールドｐk,1という結果をもたらすように前に誘導された拡張済み変位フィールドｐ0,1を使用する、本明細書で説明されるように見当合わせ処理を受ける。本明細書で説明されたようなＩ2ο（ｉｄ＋ｐ）によって表される動き修正済み画像を生成するために、画像Ｉ2の最後のワープを実行するために使用される最後の歪みフィールドという結果をもたらす画像Ｉ12およびＩ22に対して前と同様なステップが実行される。
【００４３】
このアルゴリズムの固有パラメータを下記のもの：
ピラミッドレベルのレベル数と、
各レベルにおける反復の回数と、
変位フィールドの平滑さを制御する正規化パラメータσ（このパラメータは以下に論じられる）と、に設定することだけが必要であることは留意される。
【００４４】
本発明によるガウス重み付け最小二乗平均アルゴリズムの更に詳細な説明が次に述べられ、続いて本発明の原理による特徴マップ計算の説明が行われる。
【００４５】
最初に、特に見当合わせ技法の核心をなす考えが説明される。下記の反復手順は、所与の解像度で実行される。下記の記号において解像度レベルｌは省略されている。
【００４６】
Ｇσ（ｘ）を標準偏差σの等方性３次元ガウスカーネル：
【数１４】
2008511395.tif 000008
であるとする。
各ボクセルは、インデックスｊによってラベル付けされている。その座標は、ベクトルｘj∈Ｒ3によって与えられる。本発明者らは、ｐjk∈Ｒ3をアルゴリズムのｋ番目の反復の後に点ｘjで回収された変位ベクトルと呼ぶ。１セットの変位ベクトルｐjkが与えられれば、反復ｋにおける浮動特徴マップの修正済みバージョンはＪ2,kと表される。
【数１５】
2008511395.tif 000009
各ボクセルｊにおいて、非線形最小二乗平均問題：
【数１６】
2008511395.tif 000010
のｓjk解によって更新ルールｐjk+1＝ｐjk＋ｓjkを定義することが提案されている。ガウス重み付けは、変位を推定するために使用される近隣サイズの制御を（暗に）可能にする。副産物としてこれは、歪みフィールドの全体的平滑さも保証する。
【００４７】
実際に前述の問題の線形化バージョンは、下記の基準：
【数１７】
2008511395.tif 000011
（
【数１８】
2008511395.tif 000012
は
【数１９】
2008511395.tif 000013
として定義される）を最小化することによって解決される。
【００４８】
Ｅ（ｓjk）の第１の変動を計算することによって最適性の必要条件：
【数２０】
2008511395.tif 000014
は、ｓjkに関して閉じた形式の解：
【数２１】
2008511395.tif 000015
を与える。
【００４９】
マトリックス
【数２２】
2008511395.tif 000016
は対称であって、６個の係数によって特徴付けられ得ることに
留意すること。ベクトル
【数２３】
2008511395.tif 000017
は、３個の係数によって定義される。
【００５０】
一旦これらの係数が計算されると、変位更新を推定することは、ガウスカーネルによってマトリックスとベクトル係数との畳み込みを実行することと、各格子点における３×３対称線形問題を解くことと、になる。
【００５１】
畳み込みは、刊行物、Ｒ．Ｄｅｒｉｃｈｅの「Ｆａｓｔ  ａｌｇｏｒｉｔｈｍｓ  ｆｏｒ  ｌｏｗ−ｌｅｖｅｌ  ｖｉｓｉｏｎ」、ＩＥＥＥ  Ｔｒａｎｓａｃｔｉｏｎｓ  ｏｎ  Ｐａｔｔｅｒｎ  Ａｎａｌｙｓｉｓ  ａｎｄ  Ｍａｃｈｉｎｅ  Ｉｎｔｅｌｌｉｇｅｎｃｅ，１（１２），１９９０，７８−８８に提案されているような第２次再帰フィルタを使用して近似される。これらのフィルタは、ガウス平滑化ならびにその誘導体を近似するために必要とされる計算的努力を軽減する。このアプローチによって演算は、フィルタのサイズとは独立に１出力点ごとに一定回数の乗算と加算とによって実行される。
【００５２】
３×３線形問題のために閉じた形式の解が使用される。
【数２４】
2008511395.tif 000018
が対称マトリックス
：
【数２５】
2008511395.tif 000019
によって与えられれば、その逆数は
【数２６】
2008511395.tif 000020
となる。ここで
【数２７】
2008511395.tif 000021
【００５３】
前述の全ての計算において、格子点以外で値が必要とされるときにはトライリニアー内挿方式が使用されている。
【００５４】
次に特徴マップ計算の説明が続けられる。胸部ＭＲにおいて造影剤取り入れは、被観測組織の強度値を大きく、また局所的に変化させ得る。見当合わせのために強度に依存するよりむしろ、元のデータのラプラシアンとして関数Ｊ1およびＪ2の値が取られる。この考えは、造影剤取り入れ量によって影響されることの少ない画像のエッジ情報（ラプラシアン特徴付けエッジのゼロ交差点）に集中することである。このアプローチは、実際に極めて有効であることが分かった。ここでデカルト座標系を使用してボクセルを識別するために整数のトリプレット（ｉ1，ｉ2，ｉ3）が使用され、ラプラシアンのために下記の有限差分方式が適用される。
【数２８】
2008511395.tif 000022
【００５５】
ラプラシアン演算子を推定するために使用される中心差分方式がボクセルの異方性を考慮していることに留意すること。係数λは、第３の軸に沿ったボクセル間隔（スライス間の距離）と第１の２本の軸に沿ったボクセル間隔（本明細書では同じであると仮定されている）との間の比率に対応する。
【００５６】
Ｉ1からＪ1の計算に同じ方程式が適用される。
【００５７】
明らかになるように本発明は、プログラムされたディジタルコンピュータの使用と応用によって実現されるように意図されている。図２は、双方向データ通信のために入力装置と、出力装置と、プログラムおよびその他のデータを記憶するためのメモリ装置と、に接続されたディジタルプロセッサを基本的概略形式で示している。入力装置は、本発明による処理のために適当な一つまたは複数の画像を与えるための装置のような幅広い言葉で表される。例えば入力は、ＣＡＴＳＣＡＮ、Ｘ線装置、ＭＲＩまたは他の装置に組み込まれた装置のような画像形成装置から、または記憶された画像からのものである可能性があり、あるいは直接接続による他のコンピュータまたは装置との通信によるもの、または例えばワールドワイドウェブまたはインターネットまたは他の任意の適当なこのようなデータ源を経由したような、変調された赤外線ビーム、無線、地上回線、ファクシミリまたは衛星との通信によるものである可能性がある。出力装置は、陰極線受像管、プラズマディスプレイ、液晶ディスプレイなどのような任意の適当な装置を使用するコンピュータタイプの表示装置を含むことができ、あるいは画像を表現するための装置を含むことも含まないこともあり、また好都合であり得るように更なる処理、観察または評価のために画像を記憶するためのメモリ装置を、または図２のメモリ装置の一部を含む可能性もあり、あるいはこれは、入力装置に関連して前に述べられたようなものを含む接続または連結を利用することもできる。プロセッサは、本発明のステップを実現するために本発明にしたがってセットアップされたプログラムによって動作する。このようなプログラムされたコンピュータは、画像データの取得と伝送のために通信を経由して直ちに地上回線、無線、インターネットなどにインタフェースすることができる。
【００５８】
本発明は、ソフトウエアメモリ装置に少なくとも部分的に直ちに実現可能であって、ソフトウエア製品としてこの形式にパッケージされ得る。これは、本発明の方法を利用して画像の動き補正を実行するためのプログラムコードのために媒体上にコンピュータプログラム論理を記録したコンピュータ利用可能媒体を備えるコンピュータプログラム製品の形式であり得る。
【００５９】
本発明は、人体胸部の潜在的腫瘍のＭＲＩ検出時の時系列的画像における動き補正に関する例証的例示的実施形態を使用して例によって説明されてきたが、本発明はまた一般に、ＰＥＴ−ＣＴ見当合わせの例のような、しかしこれに限定されない他の分野において空間的整列を必要とする問題の解にも適用可能である。
【００６０】
例示的実施形態による説明が限定することを意図していないことと、本明細書で明示的に説明されていない種々の変更と代用が特許請求の範囲にその範囲が記載された本発明の精神から逸脱せずに実施可能であることとは、理解される。
【図面の簡単な説明】
【００６１】
【図１】本発明の原理による１時系列の画像からの動き修正を流れ図形式で示す図である。
【図２】本発明の実現のためのプログラム可能ディジタルコンピュータの応用を概略形式で示す図である。
【図１】
2008511395.tif 000031
【図２】
2008511395.tif 000032
【国際調査報告】
2008511395.tif 000033
2008511395.tif 000034
2008511395.tif 000035