画像登録方法及び前記画像登録方法を実施するためのアルゴリズム及び前記アルゴリズムを用いて画像を登録するためのプログラム及び対象の移動による画像アーチファクトを低減するための生体医学画像の取り扱い方法
少なくとも第一と第二デジタル画像又は同じ対象の一連の断面画像のような生物医学画像を登録する方法において、第一画像又は画像セットの中で、一定数のランドマーク、いわゆる特徴が、一定数のピクセル又はボクセルを選択することで個別化される。特徴として選択された各ピクセル又はボクセルの位置は、オプティカルフローベクトルを決定することによって第一から第二画像又は画像セット内で追跡される。第一と第二画像又は画像セットは、逆オプティカルフローベクトルを第二画像又は画像セットのピクセル又はボクセルに適用することで登録される。本発明は、第一画像又は第一の一連の断面画像の各ピクセル又はボクセルの周囲近傍のピクセル又はボクセルを定義することに存する自動追跡可能なランドマーク選択ステップを提供する。各目標ピクセル又はボクセルは、この目標ピクセル又はボクセルと、ウインドウの各ピクセル又はボクセル又はその一部の外観を記述するパラメータの関数として、そして数字で表したマトリックス又は前記ウインドウのピクセル又はボクセルを表現する前記数字で表したマトリックスの変形の一つ以上の特徴的なパラメータの関数として計算される一つ以上の特徴的なパラメータを決定する。有効に追跡可能なランドマークと同時に起こるピクセル又はボクセルが、前記目標ピクセル又はボクセルの前記特徴的なパラメータの関数として決定される。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、
a)ピクセル又はボクセルからなる2次元又は3次元アレイで構成された、少なくとも第一と第二デジタル画像又は同じ対象の一連の断面画像を提供するステップと、
b)第一画像又は画像セットの中に、ランドマークとして設定された一定数のピクセル又はボクセルを選択し、追跡すべき前記特徴のリストを生成することで、一定数のいわゆる特徴であるランドマークを定義するステップと、
c)第一画像から第二画像又は、第一から第二の画像位置間のオプティカルフローベクトルによって決定された後の時点で得られた画像又は画像セット又は、特徴として選択された各ピクセル又はボクセルの後の時点で得られた画像又は画像セットから、特徴として選択された各ピクセル又はボクセルの位置を追跡するステップと、
d)第一と第二画像又は、逆オプティカルフローベクトルを第二画像又は画像セットのピクセル又はボクセルの位置に適用することによって後の時点で得られた画像又は画像セットを登録するステップと、
からなる対象の移動によって生じる画像アーチファクトを低減する生物医学画像の登録方法に関する。
【背景技術】
【0002】
この種類の時間系列の一部としての画像登録方法は、例えば米国特許6、553、152によって知られている。
【0003】
米国特許6、553、152は、視覚的ランドマーク識別と、人間オペレータによる二つの画像の対応する画素のマーキングとによって、二つの画像に対する画像登録方法を提供している。
【0004】
この方法は、ランドマークの選択を完全にオペレーターの技能にゆだねる。そして、顕著な又は一意の認識可能な構造が画像内に存在していなければ、実際に実行することは非常に困難である。例えば磁気共鳴映像(MRI)や超音波やレントゲン写真画像といった診断画像を考えると、画像の構造的特徴によってランドマークの視覚的識別をすることは、極めて困難であり、大きなエラーが生ずる可能性がある。
【0005】
種々の異なるアルゴリズムが輪郭を識別するために試みられており、例えばFischer, B.;Modersitzki, J. Curvature Based Registration with Applications to MRMammography; Springer-Verlag Berlin Heidelberg: ICCS 2002, LNCS 2331, 2002; pp 202- 206に開示されている。また、2回以上測定されたか又は異なるモダリティで調査された組織/領域の放射線画像内のランドマークの識別に関しては、例えばShi J, Tomasi C. Good Features to Track. 1994 IEEE Conference on Computer Vision and Pattern Recognition (CVPR'94) 1994; 593- 600や、Tommasini T, Fusiello A, Trucco E, Roberto V. Making Good Features Track Better Proc. IEEE Int. Conf. on Computer Vision and Pattern Recognition (CVPR '98) 1998; 178- 183や、Ruiter, N. V.; Muller, T. 0.; Stotzka, R.; Kaiser, W. A. Elastic registration of x ray mammograms and three-dimensional magnetic resonance imaging data. J Digit Imaging 2001, 14, 52-55 が開示されている。
【0006】
異なる時間に撮影された同一対象の二枚の平面画像を登録して、時間に対する対称のう動きを計算することは、二次元登録と呼ばれる。同様のことを画像の三次元セット、例えばMRI又はコンピュータ断層撮影(CT)スキャナの断面画像セットに対して行うアルゴリズムは、三次元と称される。ランドマークや対象の移動は、これら画像セットのあらゆる方向に生じる場合がある。
【0007】
比較される画像の異なるモダリティの数に依存して、アルゴリズムは単一様式の、そして複数様式の登録アルゴリズムに分割される。
【0008】
二枚の連続したMRI画像を比較することは、単一様式の登録である。
【0009】
通常、単一様式登録アルゴリズムは、両画像内の明確なランドマークや複雑な形状を識別しようとし、ほとんどの場合、二次元上で動作する。二番目のステップにおいて、両画像の対応するランドマークの位置が比較され、それらの運動ベクトルが計算される。この情報は、動きアーチファクトを除去するために、二番目の画像の全ての画素を新しい位置にシフトバックするために用いられる。
【0010】
「剛体レジストレーション」は、単一ユニットとしての二次元又は三次元画像/ボリュームを一定の方向に移動させる。「非剛体レジストレーション」は、一定のエリア/ボリュームの単一又は複数の画素を異なる方向へ移動させる。
【0011】
例えば非常に弾力性のある乳房組織に関し、動きアーチファクトを除去するために非剛体レジストレーションが必要である。
【0012】
骨のような一定の構造をとらない特別な柔らかさゆえに、乳房に直接隣接した部分は、異なる方向に移動可能である。従って、どの登録アルゴリズムも、他の既存のアプローチも、乳房組織全体に広がるランドマークの識別が必要である。
【0013】
アルゴリズムの一つのグループは、限定された数の発見すべきランドマークを定義する。これらは、その数価によって分類される。乳房のMRIイメージに関し、MRIにおいて脂肪と空気の間に高いコントラストがあるため、これらのアプローチはしばしば乳房の外側境界や乳房構造において多数のランドマークを見つける(Huwer, S.; Rahmel, J.; Wangenheim, A. v. Data-Driven Registration for Lacal Deformations. Pattern Recognition Letters 1996, 17, 951-957)。
【0014】
圧縮パネルを用いるか、少なくとも両乳房を腹臥位に横たえて固定することで、通常圧縮装置に直接接する皮膚のすべての部分は良く固定される(図1)。対照的に、心臓の鼓動や呼吸によって、両乳房の内部領域は僅かに動き続ける。外側の皮膚と空気の境界と比較して、脂肪と柔組織間のコントラストは非常に低いので、乳房の中心にある有効なランドマークの数は、これらのアルゴリズムにおいて低いままである。
【0015】
ランドマーク選択における一般的な問題は、選択したランドマークが非常に少ないことが、不十分な又は不正確な登録につながるという事実に存する。しかし、人間の観察者もプログラムも、例えば乳房中央の異なる組織タイプのより低いコントラストによって、より再割付けが不確実であるランドマークを含めざるを得ないので、追加のランドマークを選択することは、必ずしも正確な登録を保証しない。ランドマークの数を増やすことは、常に登録の計算の複雑さを著しく増加させる。
【0016】
いわゆる非剛体レジストレーションに特に適した特徴選択と特徴追跡アルゴリズムは、いわゆるLucas & Kanadeアルゴリズムと、Lucas BD, Kanade T. An Iterative Stereo Vision, IJCAI81; 674-679 1981 及び付録1として下に同封したBouguet J-Y. Pyramidal Implementation of the Lucas Kanade Feature Tracker, TechReport as Part of Open Computer Vision Library Documentation, Intel Corporationに詳細が開示された、特にピラミッド形の実施である。Lucas & Kanade技術とそのピラミッド形の特徴追跡実施は、画像内の信頼できるランドマークの自動認識や、異なる時間に撮影された二枚の画像間の追跡に、特に適している。選択された特徴の置換又は移動は、シフトとして又はいわゆるオプティカルフローベクトルとして求められる。この技術は、例えばロボット工学におけるコンピュータ・ビジョンに広く適用され、デジタル画像処理といくつかの基礎コースで教えられるコンピュータビジョンにおける当業者の一般的な知識の一部である。
【発明の開示】
【発明が解決しようとする課題】
【0017】
しかしLucas & Kanade技術を実用的に特に乳房のデジタルの又はデジタル処理された診断画像に適用することを考慮すると、前記技術は興味の範囲外も含めた非常に多くの特徴を識別し、例えば乳房まわりの空気のノイズの影響を受ける。
【0018】
以下の例は、上記問題をより詳細に明らかにする。
【0019】
コントラストを強調した(CE)磁気共鳴乳房X線撮影(MRM)は、浸潤性腫瘍に対する鋭敏な感度によって有用で重要な検査である。それは乳房柔組織の密度と独立な高い陰性適中率を示す。最近のマルチセンターな研究は、画像解析技術に従って、感度と特異性が異なることを示す。1.5Tスキャナを用いて、96、93、又は86%の感度と30、50、又は70%の特異性がそれぞれ計算された。従って、MRMの主な欠点は、偽の確定診断が高い割合で残存していることである。
【0020】
通常MRイメージングは、患者が腹臥位の状態で、製造業者の二重乳房コイル(米国の場合、一重乳房コイル装置も多い)を用いて1.0又は1.5テスラの撮像装置で実施される。このプロトコルはあらゆる方向において30から180秒の間の収集時間で、少なくともダイナミックコントラスト強調されたT1加重したシーケンスを含む。gadopentetate dimeglumine又は他の常磁性MR造影剤のボーラス注入の前に、少なくとも一回測定が実施され、その後に数回の測定が実施される。
【0021】
乳房組織の領域への実際の造影剤の摂取を計算するために、造影剤の使用前の各ボクセルの信号強度は、造影剤の使用後の信号強度から差し引かれなければならない。MR画像のスライス厚さによる部分容積効果と同様に、呼吸と心臓の鼓動によって生じる両乳房の最低限の動きによって、異なる時間に撮影された同じ画像位置のボクセルは、乳房組織の正確に同じボリュームを示さない。従って(腫瘍の場合を除いて)引き算した画像は完全に黒くはならない。小さな移動の効果は、乳房の外側境界で最も良く現れる。脂肪組織の信号強度は大きく、周囲の空気の信号強度はほとんどゼロであることから、非常に小さな移動が脂肪を先の画像の空気の位置に置換する。脂肪の大きな信号強度から空気のゼロの信号強度を引き算することで、CM用途における大きな信号上昇をシミュレートする。結果として、引き算画像が、画像化された乳房の少なくとも部分まわりの強い造影剤の摂取を表す太く白い線を示す(図2−乳房の白い外側境界が表現されている)
【0022】
デカルト座標系によって定義された三次元又は容積測定用の画像収集を考慮すると、xとy方向のあらゆる移動に加えて、(胸に向かっての)z方向の幾らかの移動が常に見られる。これは、一連の検査中の患者の弛緩によって生じる。
【0023】
アーチファクトを低減するために、乳房は通常何らかの圧縮方法(製造業者に依存する)によって乳房コイルの内側に固定される。この事実にもかかわらず、最小限の移動は常に発生する。従って、乳房を固定しているにもかかわらず、画像には常にx、y、x方向に僅かな動きアーチファクトが見られる。この動きアーチファクトは、明るい画素として引き算画像の中で視覚化される。動きアーチファクトがなければ、腫瘍性組織が占める領域を除き、引き算画像は絶対的に黒くなる。
【0024】
周囲の空気のノイズ入りの信号内で、一定数のランドマークが乳房の外側で検出される。各特徴とみなされるものは追跡しなければならず、コンピュータに不必要な負荷を発生させるという事実の他にも、ノイズは変量効果であり、ノイズの中で見つけられたすべての特徴又はランドマークは、二番目の又はそれに続く画像において対応するランドマークを有さない。しかし実は、ノイズの中から発見されたランドマークに対応する二番目の、又はそれに続く画像又は画像セットにおけるランドマークを調査するアルゴリズムは、一定量のランドマークを第二画像又は画像セット内に再割付けする。これが、画像の不正確な登録につながり、結果に影響を与える。
【0025】
現在知られている方法によれば、「特徴選択」は以下のように実行できる。
B1)第一画像又は第一画像セットの各ピクセル又はボクセル周りに隣接するピクセル又はボクセルを定義する。ここで、前記隣接するピクセル又はボクセルは、限られた数のピクセル又はボクセルからなり、単一画像の場合には二次元隣接部又はボクセルが選択され、断面画像セットの場合には、三次元隣接部又はボクセルが選択される。
B2)各ピクセル又はボクセルに対して、前記隣接部又はボクセルのすべてのピクセル又はボクセルの平均信号強度を定義する。
B3)平均信号強度閾値を定義する。
B4)ステップB2で各隣接部又はボクセルに対して定義した平均信号強度を、前記平均信号強度閾値と比較する。
B5)前記隣接部又はボクセルの平均信号強度がステップB4の閾値より大きい場合に、そのピクセル又はボクセルは追跡すべき特徴として定義され、追跡すべき特徴のリストに追加される。
【0026】
平均信号強度閾値は、調整可能であり、変更可能である。経験的又は実験的基準は、例えばテストセットやサンプル画像に対しこの方法を実行し、異なる閾値に対して得られた結果を比較することによって、このような閾値を求めるために適用可能である。
【0027】
自動特徴選択が、例えば二次元画像又は二次元の場合の一般化としての三次元画像に対していわゆるLukas&Kanadeアルゴリズムを応用することによって、現在用いられている。良く知られたLukas&Kanadeアルゴリズムによる自動特徴選択法は、Lucas BD, Kanade T. An Iterative Image Registration Technique with an Application to Stereo Vision, IJCAI81; 674-679 に非常に詳細に記載されている。
【0028】
三次元画像や三次元データボリュームを表す断面画像セットに用いられる自動特徴選択は、Lukas&Kanadeの二次元特徴検出アルゴリズムの実質的な拡張であって、以下により詳細に説明する。
【0029】
本明細書と請求項において、「特徴」という用語は「有効に追跡可能なランドマーク」という用語と等価である。
【0030】
追跡すべきランドマーク又は特徴を自動的に選択するための上記方法を用いて、対象の画像の登録が実行可能である。登録のおかげで、動きの差異、又は動きアーチファクトを補償することができる。しかし、登録された画像シーケンスは、十分な外観を有さない。上記方法で画像を処理した後、画像はシャープさが減ったように見え、又は他の効果が見える。
【0031】
さらに、有効なランドマーク又は特徴の自動選択のための方法は、手作業で、即ち人間オペレータの技能に基づいて実行されるのを避けて、客観的な基準に従って、画像内で追跡されるべき有効なランドマークを選択することを可能にする有用なツールではある。しかし、有効なランドマークの識別は、画像の各ピクセルを目標ピクセルとして考慮し、目標ピクセルの周囲の一定数のピクセルを備えた各目標ピクセルまわりの窓を定義することによって求められた勾配行列の固有値の関数であるところの特定の基準にのみ委ねられている。
【0032】
この方法ステップは、一種の境界検出であって、境界の特徴を有するとして識別された目標ピクセルの中で、幾つかのみが登録プロセスにおいて有効に追跡可能なランドマークとして選択される。
【0033】
本発明の目的は、特に二次元又は三次元の生物医学的診断画像に対して画像を登録する方法を提供することに存する。この方法は、二次元及び三次元で動きを追跡可能にすることによって、また特にデジタル画像又はデジタル処理された画像内の対象領域の外側のノイズによるアーチファクトを低減することによって、また隣接する組織タイプの信号強度における高低から独立な、乳房全体にわたりランドマークを広げることによって、既知の登録方法の欠点を効率的に克服している。
【0034】
本発明は、特に画像内の追跡すべき有効なランドマークを選択するステップを向上させることを目的としている。
【課題を解決するための手段】
【0035】
本発明による画像を登録するための方法は、
a)組織又は組織領域又は解剖した領域のMRI,超音波又はレントゲン写真の走査により得られた第一のデジタル又はデジタル処理された画像又は断面画像セットを提供し、
ここで前記画像はピクセル又はボクセルの二次元又は三次元アレイによって形成され、
同じ組織又は組織領域又は同じ解剖した領域のMRI,超音波又はレントゲン写真の走査により得られた少なくとも第二デジタル又はデジタル処理された画像又は同じ解剖した領域の断面画像セットを提供するステップと、
b)ランドマーク又は特徴として設定された一定数のピクセル又はボクセルを選択し、追跡すべき前記特徴のリストを生成することによって、第一画像又は画像セットの中に、一定数のランドマーク、いわゆる特徴を定義するステップと、
c)特徴として選択された各ピクセル又はボクセルに対する第一から第二画像又は画像セットのオプティカルフローベクトルを求めることによって、第一から第二画像又は画像セット間で、特徴としての各ピクセル又はボクセルの位置を追跡するステップと、
d)逆オプティカルフローベクトルを第二画像又は画像セットのピクセル又はボクセルに適用することによって、第一と第二画像又は画像セットを登録するステップと、からなる。
【0036】
ここにおいて、
B1)第一画像又は第一の一連の断面画像の各ピクセル又はボクセル周りに隣接した、限られた数で構成されるピクセル又はボクセルを定義するステップと、
B2)各目標ピクセル又はボクセルに対し、
外観を記載する数値パラメータ、すなわち前記ピクセル又はボクセル窓の前記目標ピクセル又はボクセルとピクセル又はボクセルの部分とからなるいわゆる数値外観パラメータの関数として計算され、又は、
前記窓のピクセル又はボクセルを表現する数値パラメータのマトリックス又は数値パラメータの前記マトリックスの変形で構成される一つ以上の特徴パラメータの関数として計算される、一つ以上の特徴のあるパラメータを定義するステップと、
B3)前記目標ピクセル又はボクセルの特徴的なパラメータの関数としての有効に追跡可能なランドマーク又は特徴に存するピクセル又はボクセルを定義するステップと、
からなる自動追跡可能なランドマーク選択ステップが実行されることを特徴とする。
【0037】
本発明の第一実施例によれば、目標ピクセル又はボクセルが追跡すべき有効なランドマーク又は特徴であるかを求めるための関数は、それぞれの値又は目標ピクセル又はボクセルの特徴パラメータと閾値との比較に存する。
【0038】
閾値は、実験によって求められ、テストセット又はサンプル画像に対して方法を実行し、異なる閾値に対して得られた結果を比較して求められる。
【0039】
目標ピクセル又はボクセルが追跡すべき有効なランドマーク又は特徴であるかを求めるための代替的な方法は、分類アルゴリズムの入力ベクトルであるところのベクトル形式の目標ピクセル又はボクセルの特徴パラメータを組織化することに存する。
【0040】
コンポーネントとして目標ピクセル又はボクセルの一定数の特徴パラメータを備えたベクトルを処理するために、前記分類アルゴリズムは、前記タスクを実行するために具体的に訓練される必要がある。
【0041】
従って、上記代替方法による本発明の方法は、
一定数の有効なランドマークが、有効に追跡可能なランドマーク又は特徴として識別された既知の場合の一定数の画像を提供するステップと、
上記ステップB1)とB2)に存する自動追跡可能ランドマーク選択ステップを適用することによって、一定数の画像の中で有効に追跡可能であると識別された前記ランドマークに対応して前記ピクセル又はボクセルに対する特徴パラメータセットを求めるステップと、
有効に追跡可能であると識別された各ランドマークに一意的に関連付けられ、成分として前記有効に追跡可能なランドマークと一致する前記ピクセル又はボクセルの特徴パラメータを備えたベクトルを生成するステップと、
変数の所定の数値によって有効に追跡可能な質を表現し、
前記数値を、有効に追跡可能なランドマークと一致する各ピクセル又はボクセルをコード化するベクトルに関連付けるステップと、
有効に追跡可能なランドマークと一致する各ピクセル又はボクセルをコード化する各ベクトルが、分類アルゴリズムのためのトレーニングデータベースの登録を形成するステップと、
前記データベースによって分類アルゴリズムをトレーニングするステップと、
鍛えられた分類アルゴリズムの入力に前記目標ピクセル又はボクセルの前記特徴パラメータを備えたベクトルを提供することによって、目標ピクセル又はボクセルの有効に追跡可能なランドマークの質を求めるステップと、を備える。
【0042】
改善した本発明によれば、前記トレーニングデータベースは、有効に追跡可能ではないと知られたランドマークと一致する画像のピクセル又はボクセルのベクトル形式に組織化され、有効に追跡可能でないランドマークの質を、追跡可能なランドマークの質を表現する前記変数の値とは異なる変数の所定の数値によって表現する特徴パラメータも備え、一方で有効に追跡可能でないランドマークの質を表現する前記数値を、既知の有効に追跡可能でないランドマークと一致した各ピクセル又はボクセルをコード化するベクトルに加える。
【0043】
分類アルゴリズムを用いた有効に追跡可能なランドマークを求めるためのこの種の自動的方法は、例えば適応度や他の同様な統計指標のような分類の信頼性を示す統計的尺度をも提供するという更なる長所を有する。
【0044】
分類アルゴリズムとして、あらゆる種類のこれらアルゴリズムを用いることができる。分類アルゴリズムの特別な場合としては、いわゆるクラスタ化アルゴリズム又は人工の神経ネットワークがある。
【0045】
目標ターゲットピクセル又はボクセルや、前記ピクセル又はボクセル窓のピクセル又はボクセルのそれぞれ、又は一部の外観を表現する数値外観パラメータの関数として、又は前記窓のピクセル又はボクセルを表現する数値行列又は前記数値行列の転換の一つ以上の特徴パラメータの関数として計算される特徴数値パラメータに関し、多くの異なる関数が多くの異なる組み合わせの中で使用可能である。
【0046】
特徴又は有効に追跡可能なランドマークに対応する目標ピクセル又はボクセルの特徴パラメータを求めるための第一の関数は、前記隣接部又はボクセルのすべてのピクセル又はボクセルの平均強度値に存する。
【0047】
目標ピクセル又はボクセルのもう一つの特徴パラメータは、ステップB1)で定義された前記窓のピクセル又はボクセル、即ち目標ピクセル又はボクセルに新設するピクセル又はボクセルの一定数を表現する数値行列の特徴パラメータに存する。前記特徴パラメータは、前記窓のピクセル又はボクセルの画像データを備えた数値行列の特異値である。
【0048】
代替的に、又は上述の特徴パラメータとの組合わせにおいて、目標ピクセル又はボクセルのもう一つの特徴パラメータが、前記窓のピクセル又はボクセルを表現する数値行列の勾配行列の固有値に存する。
【0049】
目標ピクセル又はボクセルのもう一つの特徴パラメータが、前記窓のピクセル又はボクセルを表現する数値行列のヘッシアン行列の固有値に存する。
【0050】
代替的に又は上記開示されたものとの組合わせで提供された目標ピクセル又はボクセルの特徴パラメータのためのもう一つの選択が、前記窓のピクセル又はボクセルに対応するデータの数値行列のウェーブレット変換係数の一つ以上に存する。
【0051】
この場合、いくつかのウェーブレットを基礎にした関数が、代替的に又は組み合わせとして用いるために選択可能である。
【0052】
ウェーブレット変換のより詳細な記載は、付録1で与えられている。この付録は、インターネット上で利用可能なように公開されており、"Wavelet for Kids, A tutorial introduction"というタイトルで、Duke大学のBrandi Vidakovic とPeter Muellerが著している。この文献では、ウェーブレット理論が要約され、議論され、画像処理への幾つかの応用が開示されている。画像処理でのウェーブレットの実行を開示した章から明らかに、ウェーブレット分解はパラメータを得ることを可能にする。分解の各レベルに対して、ウェーブレット変換は、平均を表す一つの行列と、いわゆる詳細を表す3つの行列とを生成する。上記行列の一つ以上から、例えば行列の要素の平均を取る以外にも、行列の特異値を取ることなどのよって、幾つかのパラメータを抽出することが可能である。これらのパラメータのすべて又は幾つかは、目標ピクセル又はボクセルの特徴パラメータとして用いることができ、また、選択された窓に含まれる周辺ピクセル又はボクセルと関係する各目標ピクセル又はボクセルを表すベクトルの成分を形成するために用いることができる。
【0053】
代替的に又は上記開示されたものとの組合わせで提供された目標ピクセル又はボクセルの特徴パラメータのためのもう一つの選択が、前記窓のピクセル又はボクセルを表す数値行列の自己相関変換係数の一つ以上に存する。
【0054】
画像処理における自己相関は、一般的に画像テクスチャを特徴づけるための道具として用いられ、二つの画像の数学的評価から成る。この二つの画像は、異なる時刻に撮影されるか、又は第一画像から場所を移動した結果としての第二画像として生成される。自己相関は、これら二つの画像間の関係を測定する。この数学的評価は、目標ピクセル又はボクセルのコード化において、考慮されるパラメータの数を減らす可能性を提供する。
【0055】
目標ピクセル又はボクセルのもう一つの特徴パラメータは、前記窓のピクセル又はボクセルを表す前記数値行列の同時生起行列(又はその特異値)に存する。
【0056】
前記同時生起行列は、固定された空間的関係によって分離された一対のピクセルに対する中間調の二次元ヒストグラムである。この行列は、一対のピクセルの共有の確率分布を近似する。一般的にコントラストやエントロピーのような、テクスチャ測定を計算するために用いられる。
【0057】
目標ピクセル又はボクセルの特徴パラメータの上記例によれば、各目標ピクセル又はボクセルは、前記窓のピクセル又はボクセルを表現する数値行列の固有値又は特異値及び/又は、前記窓のピクセル又はボクセルを表現する数値行列の勾配行列又はヘッシアン行列の固有値及び/又は、ウェーブレット変換係数の一つ以上の係数及び/又は、自己相関変換の一つ以上の係数及び/又は、同時生起行列の一つ以上の入力又は特異値及び/又は前記窓のピクセル又はボクセルと前記目標ピクセル又はボクセルの平均強度値、の二つ以上の組み合わせ又はサブコンビネーションによって記載される。
【0058】
既に開示したように、前記特徴パラメータは、前記特徴パラメータのそれぞれに対し、又はその組合わせに対し閾値を定義することによって、画像における目標ピクセル又はボクセルが、有効に追跡可能なランドマークとして考慮可能かどうかを判断するために用いることができる。
【0059】
もう一つの基準は、前記特徴パラメータをベクトル形式に組織化することに存し、そこで各パラメータは前記ベクトルの成分であり、前記ベクトルの係数を求め、この値を前記係数に対する閾値と比較することで、組織化する。
【0060】
代替的に、前記ベクトルは既に上で開示した分類アルゴリズムによって処理することもできる。
【0061】
各ピクセル又はボクセルと、選択されたピクセル又はボクセル窓によって定義された選択された隣接部又はボクセルによって記載された数値間の関係は、単に各ピクセル又はボクセルの外観を表す値に存するオリジナルの数値行列の前記異なる変換の前記特異値及び/又は固有値によって、及び/又は、ウェーブレット変換の一つ以上の係数及び/又は自己相関変換の一つ以上の係数及び/又は同時生起行列の一つ以上の入力又は特異値
によって要約される。そして、各ベクトルにおいて、この関係は強調表示に特に適した、又は選択された窓内の画像のピクセル又はボクセル間の関係の特定の種類に感受性のある数値によって定義される。
【0062】
従って、ピクセル又はボクセルに対する有効に追跡可能なランドマークの質を求めるためのこれら特徴パラメータの一部又は全てを考慮することは、このピクセル又はボクセルと残りの画像との関係に関しより良くより完全に考慮する助けになる。
【0063】
すでに述べたように、上記方法は磁気共鳴画像(MRI)又はX線CTにおける同一の又は異なる種類のスキャナで撮影された二枚の診断画像又は3次元画像ボリュームを登録するのに用いることができる。
【0064】
上述の登録方法は、例えばMRIや超音波画像や特にコントラストを強調した磁気共鳴マンモグラフィー(MRM)といった造影剤強調診断画像との組み合わせによっても提供される。この方法において、乳房組織領域の実際の造影剤摂取を計算するために、造影剤使用前の各ボクセルの信号強度が、造影剤使用後の信号強度から差し引かれる必要がある。このように、画像化された組織に造影剤の存在しない時に撮影された第一画像と、画像化された組織に造影剤が存在しかん流した第二の時間又は後の時間とを登録することは、登録が実行されなかった場合に患者の動きにより引き算画像に現れるアーチファクトを非常に成功裏に抑制するのに有用である。前記二枚の画像の引き算は、造影剤が溜まった組織領域に対応するピクセル又はボクセルの信号強度を高めることに繋がる。造影剤の溜まっていない他の領域は、低い信号強度のピクセル又はボクセルであって、特に暗い灰色か黒に見える。
【発明を実施するための最良の形態】
【0065】
画像の登録は、特に医療診断の分野において、幾つかの種類の調査の実行を可能にする重要なプロセスである。
【0066】
いわゆる画像化された領域のかん流挙動の測定において、登録が重要であることが以下の例により明らかになる。この場合、異なる時刻に撮影された身体上の領域の一連の画像が得られる。最初の画像は、造影剤が身体上の領域にない状態で撮影される。続く画像は、造影剤が身体上の領域に達した状態で撮影される。造影剤は、身体上の領域の周辺組織に対して血管やリンパ管の流れの応答を強調し、この流れがより良く見えるようにする。腫瘍性組織や他の病巣のような病変によって影響された組織は、強化された血管新生によって特徴付けられるので、画像化された組織のかん流挙動の測定は、この病理組織を個性化し、画像においてより良く見えるようにレンダリングするための手段である。画像化された身体上の領域の幾つかの領域における強化された血管新生は、この領域において造影剤が他の領域より集中するという結果をもたらす。また、造影剤は投与した時から血管やリンパ管のシステムによって組織に拡散されるまでに時間を要するので、強化された血管新生は、この強化された血管新生が存在する領域において、造影剤をより迅速に供給し、同時に、画像化された身体上の領域の通常に血管新生された領域に比較して、造影剤の最大ピークに達した後に、造影剤の集中のより迅速な減少を提供する。
【0067】
比較的長い時間間隔に広がった画像の時間系列が得られるので、画像化された身体上の領域の動きは、一枚の画像とそれに続く画像間で生じることができる。この動きは、患者の動きによって、また例えば心臓の鼓動や呼吸といった生理上の理由によって決定される。
【0068】
かん流挙動は、造影剤の画像化された身体上の領域への投与後に得られた時間系列の各画像から前記画像化された身体上の領域に造影剤の存在しない時点での画像を引き算することによって求められる。図2に示すように、造影剤が存在する領域が明るいままであるのに対して、この引き算が実質的に造影剤のかん流の生じない領域のピクセルを黒かほとんど黒色にするという事実のほかに、(上で定義した)各引き算画像に対して平均強度を求め、グラフに一連の画像の各引き算画像の平均強度をプロットすることで、造影剤が検出された領域のかん流挙動の尺度である、いわゆるかん流曲線を引くことができる。
【0069】
診断に関するフォローアップ又は治療が施された際の画像化された領域の病状の進展の
フォローアップを実行するために、画像化が異なる時間間隔で実行された場合に、登録も重要である。ここで、画像を撮影するための装置に対して患者はいつでも位置決めしなおす必要があるので、先のかん流調査の場合よりも、一つの画像とそれに続く画像間の時間間隔はより長くなり、一つの画像とその他の画像間の動きは、よりドラマチックになる。
【0070】
他の場合において、異なる画像化技術で撮影された身体上の領域の画像は登録される必要がある。相関画像又は組み合わせ画像がこの対象に関するより完全で詳細な情報を、又はより良いより明確な組合わせ画像を提供できるように、異なる画像技術はしばしば身体上の領域の画像化された対象の異なる様相や質を強調する又は暴露するので、これは重要である。
【0071】
強固な構造が画像化される場合には、登録は一連の画像の一枚の他に対する単純な平行移動又は回転平行移動とみなすことができる。画像化された対象が柔組織からなるか、柔組織に埋め込まれている場合、画像化された対象や画像化された領域に歪みが生じるので、登録にはこの効果も考慮することが必要である。
【0072】
診断画像の時間系列を登録する一般的なケースは、かん流挙動測定法を適用することで実行される乳房中央を調査する場合である。
【0073】
図1に患者の乳房を固定するための一般的な装置を示す。通常、MR造影は1.0又は1.5テスラの撮像装置に製造業者の二重乳房コイル(米国の場合、単一乳房コイルである場合も多い)を用いて行われ、その間患者はうつ向けになっている。このプロトコルはあらゆる方向において30から180秒の間の収集時間で、少なくともダイナミックコントラスト強調されたT1加重したシーケンスを含む(個々のプロトコルは異なっていてもよいし、追加的に脂肪抑圧を用いることもできる)。0.1 mmol/kgから0.2 mmol/kgの gadopentetate dimeglumine又は他の常磁性MR造影剤のボーラス注入の前に、少なくとも一回測定が実施され、その後に数回の測定が実施される。
【0074】
少なくともMRI画像の最初のボリュームが造影剤が肘正中静脈に注射されて乳房組織に到達する前の時間t0に撮影される。少なくとも画像のもう一つのボリュームが造影剤注射後の時間t1に直接撮影される。通常、時間遅れのある連続したボリュームが撮影される(例えば、造影剤投与後、6回以下の同一ボリュームが連続する)。
【0075】
図2に、断面画像の3次元ボリュームから取り出された2枚の画像を示す。一枚は造影剤投与前の時刻t0に撮影されたものであり、もう一枚は造影剤投与後の時刻t1に撮影されたものである。引き算画像、即ち第二画像から第一画像をピクセル毎に引き算した結果の像も、同様に示す。ここで登録は行われておらず、例えば矢印で示された領域のように、造影剤の摂取されていない領域のピクセルは、輝いたままであることは明らかである。
【0076】
図3に、従来の登録方法とも共通の、登録プロセスの主なステップを示す。
【0077】
立方体で表現された画像ボリュームが異なる時刻t0とt1に画像化される。ここでは三角と円で表現された二つの対象が、ランドマークとして選択される。動きのせいで、画像ボリュームt0での二つの対象の位置は、t1での画像の同じ対象の位置とは一致しない。
【0078】
従って、第一ステップは、画像の中で追跡すべきランドマークを個別化することである。三角と円で表現された追跡すべきランドマークがひとたび個別化されると、いわゆる動きベクトルの測定を可能にする追跡ステップが実行される。
【0079】
これらベクトルによって、画像化された身体構造上の領域の動きを大域的運動として、また画像化された身体構造上の領域の一定の領域の局部的運動として表現する、画像中に拡がるベクトル場を生成することが可能である。逆動きベクトル場を画像に適用すると、第一画像の時刻t0と第二画像の時刻t1との間で生じる動きの効果は、補償される。乳房のかん流試験の例におけるこの登録の効果は、図19で見ることができる。ここで、左側の連続画像は、時刻t0とt1で撮影された乳房画像が登録されていない場合に関し、下側の画像は、図2の一つと対応する引き算画像である。右側の連続画像は、時刻t0とt1に撮影された2枚の画像を示すが、時刻t1に撮影された画像は、登録方法で処理済みのものである。引き算画像に対応する最後の画像において、左側の引き算画像ではピクセルが輝いていた多くの領域が消失し、又は対応するピクセルが完全に黒か灰色の様相を呈していることが明らかに見ることができる。
【0080】
本発明の目的は、以下のステップでもクリティカルな、画像の登録プロセスの第一ステップを改良することである。画像内の信頼できるランドマークを定義することは、より良い登録結果を得るためにクリティカルである。したがって、有効に追跡可能なランドマークとなる可能性の非常に高い画像のピクセル又はボクセルの特徴的性質を求めるためのより精密で完全な方法が本発明の目的である。
【0081】
図4を参照して、単純化のために二次元の場合について、デジタル又はデジタル処理された画像が、n, m = 1, ... , NであるピクセルアレイP(n, m)によって形成される。各ピクセルは、一定の外観を有するモニター又は画面又は印刷された紙上の小さなドットである。ドットの外観は、数値データによって表現可能である。いわゆるグレースケール画像において、各ドットは前記グレースケールの一定のグレーレベルに対応した一定の強度を有する。カラー画像の場合には、ピクセルの外観を数値データで完全に表現するために、より多くのパラメータが通常必要である。カラー画像におけるピクセルの外観を一義的に表現するために、幾つかのシステムが定義されている。そのうち一つのシステムは、いわゆる既知のHSV(色調、彩度、明暗度)又はいわゆるRGB(赤、緑、青)システムである。もしカラー画像を処理した場合、これらの値が対応する数値データに置き換える必要があることが当業者には自明であるため、本単純化した例においては強度I(n, m)のみが示される。
【0082】
従って、図5に示されたn, m = 1, ... , 10の数値データのアレイI(n, m)は、図1のピクセルアレイの概要図としてのn, m = 1, ... , 10のピクセルP(n, m)で構成される画像に対応し、そして、数値データのアレイは行列である。
【0083】
各ピクセルP(n, m)は、例えばグレースケールデジタル画像のピクセルのグレーレベルのように、数値的にその外観を表現した数値データI(n, m)に一義的に関連する。
【0084】
隣接する断面(二次元)画像の一定数によって形成されるボリューム画像として、三次元画像を定義することができる。この場合、画像はいわゆるボクセルで形成され、そしてアレイは画像データの数値行列も三次元であるように、三次元である。実際には、断面画像セットは定義された厚みを有し、従ってこれら画像は、純粋な二次元画像よりもよりスライスされている。この観点から、常にボクセルと称することがより正確である。それにもかかわらず、モニター画面上の断面画像の視覚的表現は、純粋に二次元画像であって、スライスの厚みは考慮されず、この場合、ピクセルが画像ドットの適切な定義である。
【0085】
図5、6、7に示すように、画像の目標ピクセルは、その成分がその目標ピクセル周りのピクセルに関する情報を含んだベクトルによってコード化可能である。通常、前記周辺領域は目標ピクセルを中心とした(2? + 1)( (2? + 1)ピクセルの次元の窓から成る。ここで、?は(? = 1, 2, ..., N)でユーザーによって任意に定義される整数であり、ここでは「ステップ」と称す。
【0086】
周囲のピクセルで形成された窓の上記定義は、勾配xの周辺ピクセルの定義と等価であり、ここでxは整数であって、この整数はステップにおける目標ピクセルから隣接部までの距離を示す。目標ピクセルとして前記窓の中心にあるピクセルを考慮すれば、勾配1の周辺ピクセルを備えた前記窓は、目標ピクセルに直接隣接したピクセルのシェルであり、勾配2の周辺ピクセルは、目標ピクセルの周囲の二番目に近いピクセルのシェルを備え、画像を形成するピクセルアレイの各方向において、勾配1は目標ピクセルから1距離ステップに対応し、勾配2は、目標ピクセルから2距離ステップに対応する。この定義は、ボクセルアレイで形成された3次元画像に対して準用する。
【0087】
前記窓の最小寸法は、3×3寸法のピクセルアレイに存し、その中心ピクセルが目標ピクセルであり、上で勾配として定義されたステップ?は、1に等しい。例えば5×5又は7×7ピクセル窓といった、より大きな窓を選択することもできて、それぞれのステップ又は勾配は、?=2と?=3である。簡単にするために、本例においては、目標ピクセルが中心にある3×3ピクセルアレイに対応した窓を選択する。
【0088】
画像のピクセルを表現する数値データ行列を参照して、図5にこの窓を示す。窓Wは九つのピクセルを備え、そのうち一つが目標ピクセルである。図5に示された窓は、目標ピクセルP(2, 2)を中心とし、ピクセルP(1, 1), P(1, 2), P(1, 3), P(2, 1), P(2, 2), P(2, 3), P(3, 1), P(3, 2), P(3, 3)を備えており、各ピクセルは対応する数値、即ち前記ピクセルの信号強度に関連する数値パラメータI(1, 1), I(1, 2), I(1, 3), I(2, 1), I(2, 2), I(2, 3), I(3, 1), I(3, 2), I(3, 3)で表現される。
【0089】
前記目標ピクセルと画像内の特定の周辺領域との関係が考慮されるように、目標ピクセルP(2, 2)は、窓Wの隣接画素に関する情報も用いてコード化される。前記目標ピクセルP(2, 2)と上で定義した周辺ピクセルとの関係を表すベクトルの成分として、前記ピクセルの強度は、目標ピクセルの強度とともに測定される。
【0090】
前記ベクトルは、図6に図式的に示される。明確に表すために、目標ピクセルを表現する各ベクトルは、九つの成分を有する。
【0091】
図7に、目標ピクセルをコード化するための前記ベクトルの一般的表現を示す。この場合、目標ピクセルはその強度によって表現され、強度I(n, m)を有するピクセルP(n, m)として定義される。
【0092】
図8に図7のベクトルを参照する窓Wと、いわゆる勾配行列におけるその変換を示す。
【0093】
当業者には自明のことであろうが、二方向が定義可能で、対象の窓における各ピクセルに対し、この二方向において勾配が評価可能であるように、ピクセルを表現する数値アレイと、この場合窓W内のピクセルの強度データアレイI(m-1, n-1), I(m-1, n), I(m-1, n+1), I(m, n-1), I(m, n), I(m, n+1), I(m+1, n-1), I(m, n), I(m+1, n+1) とは、二次元オブジェクトである。
【0094】
勾配行列は、固有値?p(p = 1, 2)で、この場合?1? ?2で表現可能なように、対角化可能である。
【0095】
選択された窓Wのピクセルの強度Iのオリジナル行列は、さらに加工可能であり、前記オリジナル行列に対し、いわゆるヘッシアン行列が計算可能である。また、この場合ヘッシアン行列は、その固有値?pによって表現可能である。
【0096】
選択された窓の強度値Iの3×3行列を考慮した場合、通常、前記行列は対角化可能でなく、固有値は上述のように有意義なものではない。これを考慮して、より一般的な場合において、図8に示すように選択された窓のピクセルの強度Iのオリジナル行列は、いわゆる特異値siで表現可能である。
【0097】
選択された窓に対応する強度行列の特異値又は代替的に、勾配行列又はヘッシアン行列の固有値と組合わせて用いることで、対応した目標ピクセルを一義的にコード化するためのベクトルを生成することが可能である。前記ベクトルの成分は、選択された窓のピクセルと、選択された窓のピクセルに対応する強度値の前記行列を処理することで得られる勾配行列又はヘッシアン行列の固有値とに対応した強度値の行列の前記特異値に存する。
【0098】
上で開示したように、図示された例は、選択された窓のピクセルに対応する強度値のオリジナル行列の変換の特別な場合に限定されるが、その他の変換も代替的に又は組み合わせとして適用可能である。例えば、ウェーブレット分解が実行可能であり、目標ピクセルのコード化ベクトルの成分として、それらのすべて又は少なくとも一部に対し、平均値と詳細値が利用可能である。
【0099】
代替的に、又は組合わせにおいて、選択された窓のピクセルに対応する強度値のオリジナル行列の自己相関変換が利用可能であり、目標ピクセルのコード化ベクトルの成分として、それらのすべて又は少なくとも一部に対し、得られたパラメータが利用可能である。
【0100】
ウェーブレット変換も前処理ステップとして利用可能である。得られた数値データの行列は、対応する勾配行列及び/又は対応するヘッシアン行列を用いて処理可能であり、その特異値又は固有値は、代替的に又は組合わせにおいて求めることができる。
【0101】
前記窓のピクセル又はボクセルを表現する行列のいわゆる同時生起変換に存する、ピクセルの窓を表現するオリジナル行列の別な変換が利用可能である。この場合、この同時生起行列の特異値又は固有値が、目標ピクセルを表すパラメータとして利用可能である。
【0102】
一定数の隣接部からなる選択された窓における目標ピクセルを特徴付けるための幾つかのパラメータの上記記載は、単に単純化のために二次元の場合に限定されているが、上述の選択は、三次元の場合にも適用可能であることを再び強調しておく必要がある。
【0103】
上記結果として、画像の各目標ピクセル又はボクセルは、上記パラメータの一つ以上によって表現可能である。これらのパラメータは、暗示的に目標ピクセル又はボクセルと選択された窓のピクセル又はボクセルとの関係を示し、それぞれ選択された窓に対応する画像の一定の特徴を現すことができる。従って、これらパラメータの二つ以上、又はすべてを、ピクセル又はボクセルが画像内の有効に追跡可能なランドマークと一致するかどうかを求めるための数値としてみなすことで、非常に信頼性のあるランドマークの識別を提供することになる。
【0104】
目標ピクセル又はボクセルと関連付けられた少なくとも上述のパラメータの一部を利用して有効に追跡可能なランドマークを選択するための第一の方法によれば、例えば閾値の数値の定義を実験データに基づかせることで、閾値を求めることができる。
【0105】
いったん閾値が上述の各目標ピクセル又はボクセルに対するベクトルの成分である各パラメータに対して定義されると、前記パラメータのそれぞれ又はパラメータのグループの組み合わせと、対応する閾値とを比較することによって、目標ピクセル又はボクセルが有効に追跡可能なランドマークと対応するかどうかを求めるための基準を設定する関数を定義することができる。
【0106】
非常に多数の関数が異なる基準に対して定義可能であり、異なる可能性が考慮される。
【0107】
最も単純だが非常に時間を要する解は、ピクセル又はボクセルをコード化するベクトルの成分である各パラメータの値に対して閾値を定義することであり、各パラメータを関連する閾値と比較して、パラメータが閾値より高いか低いかを調べることである。各パラメータに対する比較結果は、パラメータが関連する閾値より低いか高いか、又は比較結果が有効に追跡可能なランドマークかそうでないかを、例えば0か1かの数値によって表現可能である。有効に追跡可能なランドマークに対して、特徴パラメータの幾つかは負の比較結果、即ち、前記特徴パラメータがランドマークが有効に追跡可能でないことを示す比較結果を有することもあるので、特徴パラメータと関連する閾値との各個別の比較結果の数値からなる成分を有するベクトルは、各特徴パラメータとその閾値との比較結果を要約する数値からなる包括的評価パラメータを得るために、さらに処理可能である。
【0108】
このことは、特徴パラメータの一つ又はサブセットに対する比較結果がランドマークが有効に追跡可能ではないと示しているにもかかわらず、有効に追跡可能なランドマークを識別することを可能にする。
【0109】
上記特徴パラメータと関連する閾値との比較結果の包括的評価は、各比較結果の関数として、又は特徴パラメータの関数として包括的評価パラメータを算出することで求めることができる。
【0110】
例えば包括的評価パラメータは、特徴パラメータのセットの平均値として算出することができる。この包括的評価パラメータは、包括的閾値と比較することができる。この場合、包括的評価パラメータに対する包括的閾値は、各特徴パラメータと関連する各単一閾値の関数として、例えば閾値の平均値として、求めることができる。
【0111】
包括的評価パラメータを求める他の方法は、特徴パラメータにより形成された成分からなるベクトルを定義し、次に前記ベクトルの係数、即ちその長さを計算することである。包括的評価パラメータである前記ベクトルの係数と比較される包括的閾値は、各単一特徴パラメータと関連する閾値によって形成された成分からなる閾値ベクトルの係数として算出することができる。
【0112】
また別な方法は、対応する特徴パラメータの対応する閾値との比較結果が有効に追跡可能なランドマークを示すか否かに依存する、例えば0又は1といった2値を前提とする、比較評価変数によって、各特徴パラメータの対応する閾値との比較結果をパラメータ化することである。
【0113】
このように、特徴パラメータと対応する閾値との比較結果が数値的に表現可能であり、そして、包括的結果は、その成分が比較評価パラメータ値であるベクトル形式で表現可能である。包括評価パラメータが二つの離散値0と1のうち一つを有することで定義される上記例に関し、比較評価ベクトルの成分は、0と1の数値のセットから成る。
【0114】
上記オプションと同様に、包括的評価パラメータは、この比較評価ベクトルの係数として、又はその成分の平均値として計算可能である。
【0115】
上記方法の別の特徴として、特徴パラメータの一つ以上のサブセットが定義可能であり、それは前記サブセットの特徴パラメータの関数としての二次特徴パラメータを計算するために用いられる。
【0116】
あらゆる種類の関数が、二次特徴パラメータを求め、又は算出するために利用可能である。
【0117】
二次特徴パラメータを求めるための特別な方法によれば、特徴パラメータの各サブセットは、目標ピクセル又はボクセル又は前記ピクセル又はボクセル窓のピクセル又はボクセルのそれぞれ又は一部を表す数値的外観の同じ種類の関数に関する特徴パラメータと、
窓のピクセル又はボクセルを表現する数値行列又は同一の数値行列の同一の変換についての一つ以上の特徴パラメータの同じ種類の関数に関する特徴パラメータと、を備える。
【0118】
これによって、特徴パラメータのセット全体が、各サブセットの特徴パラメータの関数として求められた二次特徴パラメータによって要約される特徴パラメータのサブセットに分割される。
【0119】
実際の例によれば、目標ピクセル又はボクセルを備えた一定の窓内のピクセル又はボクセルの信号強度と、前記目標ピクセル又はボクセルの信号強度の値は、特徴パラメータのサブセットを形成することとして定義される。これらの値を要約する二次特徴パラメータは、目標ピクセル又はボクセル又は窓に含まれるピクセル又はボクセルの信号強度の平均信号強度として計算可能である。
【0120】
同一の平均値関数が、窓のピクセル又はボクセルと目標ピクセル又はボクセルを表す数値画像データの行列の特異値から、及び/又は勾配行列の固有値から、又はヘッシアン行列の固有値から、又は自己相関行列の特異値から、二次特徴パラメータを計算するために利用可能である。
【0121】
上記ステップは、一定数のパラメータを有する一次特徴パラメータセットを、より少ない数のパラメータを有する二次特徴パラメータセットに減じることに存する。
【0122】
上述のように一次特徴パラメータとして定義可能な特徴パラメータのセットの評価と同様に、二次パラメータは、ピクセル又はボクセルが有効に追跡可能なランドマークに対応するか、有効に追跡可能でないランドマークに対応するかを決定するために用いることができる。前記二次パラメータは、一次特徴パラメータに適用したものと同じステップによって処理可能である。従って、二次パラメータをその二次閾値と比較するステップが実行可能であるのと同様に、二次閾値は定義可能である。またこの場合に、包括的閾値と比較される包括的比較評価パラメータも求めることができる。
【0123】
また上述の基準に従って特徴パラメータをグルーピングすることによって得られた前記二次特徴パラメータの場合において、単一二次特徴パラメータの組み合わせと対応する閾値の組み合わせは、例えば加重和又はその他の加重した組み合わせによって作成され、そして、そのように定義された包括的二次特徴パラメータと包括的二次特徴閾値のレベルとで比較が実行可能である。
【0124】
原理的に、一次特徴パラメータとその対応する閾値との比較を実行する上述のそれぞれの方法は、二次特徴パラメータのセットに対しても適用できる。特に、例えば二次特徴パラメータが閾値より大きいか小さいかによって0と1を取るような、各特徴パラメータとその関連する閾値との間の比較に関する異なる離散値を含む比較結果ベクトルに対して提供することにメリットがある。
【0125】
次に、包括的結果を和として、特に比較結果ベクトルの成分の加重和として求めることができる。
【0126】
パラメータ数が減ることによって、計算と評価がより簡単になることが長所である。さらに、二次特徴パラメータを求めるために選択された関数の種類が、正しくない又は統計的許容誤差の外にある一つ以上の前記パラメータの値の範囲を制限することを可能にする。
【0127】
比較のために包括的評価パラメータを求める場合において、ピクセル又はボクセルが有効に追跡可能なランドマークかどうかを求める場合の妥当性に相応じて、各一次又は二次特徴パラメータを加重する加重値を定義することが可能である。図20は、上述のごとく定義された前記ピクセル又はボクセルの特徴を考慮することによって、ピクセル又はボクセルが有効に追跡可能なランドマークと対応するか否かを評価するための上述のステップを表現した概要図である。
【0128】
図20で示されたケースは、一次特徴としてここで定義された目標ピクセル又はボクセルの特徴パラメータの選択されたサブセットの関数としての、各目標ピクセル又はボクセルに対して二次特徴パラメータを求めた場合を指す。
【0129】
以下で説明することではあるが、前記処理ステップのメカニズムも、第一特徴パラメータに適用することができる。
【0130】
図4ないし8に関連して、先に説明したように、第一特徴パラメータのセットが求められる。目標ピクセルは、I(m, n)で表される。隣接するピクセルは、3×3ピクセル寸法を有する窓に含まれるものとして選択される。一次特性は、選択された窓内の隣接部の強度Iと、目標ピクセルの強度と、目標ピクセルの一つを含む選択された窓に含まれるピクセルの信号強度の行列の特異値siと、目標ピクセルの一つを含む選択された窓に含まれるピクセルの信号強度の行列の勾配行列の固有値と、目標ピクセルの一つを含む選択された窓に含まれるピクセルの信号強度の行列のヘッシアン行列の固有値と、で構成される。
【0131】
上述のように求められたあらゆる別の特徴パラメータは、図示された特徴パラメータに付加したり、置き換えたりすることができるので、図20の例で選択された特徴パラメータの種類に、本発明は限定されるものではないことを強調する必要がある。
【0132】
一次特徴パラメータの前記種類のそれぞれは、特徴パラメータのサブセットとして定義される。これは、中括弧によって図20に示される。二次特徴パラメータは、関数1として示された関数によって、一次特徴パラメータの各サブセットから求められる。対応する一次特徴パラメータのサブセットから二次特徴パラメータを求めるための関数は、各サブセット間で同一又は異なるものであってよい。
【0133】
関数の例は、二乗関数又は係数関数等の和の平均値関数又は加重関数に存することができる。
【0134】
他の種類の関数は、サブセットの一次特徴パラメータの少なくとも幾つか、又は一つだけ除いた全部を除去するために、閾値との比較を提供する。各サブセットに対し、二次特徴パラメータはそのように求めることができ、図20において、この二次特徴パラメータは、P1,P2,P3,P4で示される。ここで、目標ピクセルが有効に追跡可能なランドマークと認められるか否かを求めるために、上述のすべてのステップが実行可能である。各二次特徴パラメータに対して個々の閾値を求めることができ、そして各二次特徴パラメータと対応する閾値との個々の比較が実行可能である。ピクセルが有効に追跡可能なランドマークか否かの決定は、上述の包括的比較評価プロセスによるこれらの比較結果に基づいて行われる。
【0135】
本実施例において、包括的二次特徴パラメータが図20の関数2によって示される二次特徴パラメータの関数として求められる。関数2は、関数1又は関数1に対して示された例のうち一つと等しくすることができる。
【0136】
上述のすでに開示した方法のうち一つに従って求められた包括的閾値は、例えば各二次特徴パラメータに対する個々の閾値を求め、次に、関数2と同じ又は異なる関数によって、包括的二次特徴パラメータをもとめるために、この個々の閾値のセットから包括的閾値を求めることによって求められる。
【0137】
包括的二次特徴パラメータと包括的閾値との比較は、比較関数によって実行され、そして、比較結果はピクセルI(m, n)が有効に追跡可能なランドマークと一致するか否かの評価を与える。
【0138】
代替案として、全ての一次特徴パラメータから包括的評価パラメータを求めるために関数を選ぶことができ、そして、閾値から求めるための別の関数又は同じ関数が、各一次特徴パラメータに対して包括的閾値を定義することを認めなければならない。比較関数は、次に、この一次特徴パラメータに対する包括的閾値に適用され、また包括的評価特徴パラメータに適用される。
【0139】
好ましい実施例において、関数1と関数2として示された異なる関数が、線形結合であって、好ましくは前記一次又は二次特徴パラメータの加重線形結合であって、例えば加重和である。同様に、第二特徴パラメータに対する閾値と包括的閾値とは、各一次パラメータに対する閾値と、各二次パラメータに対する閾値それぞれの加重和として計算可能である。
【0140】
特別に長所を有する実施例は、以下のステップからなる。各一次特徴パラメータとそれらの閾値とを比較するステップと、各成分が各一次パラメータと対応する閾値との比較結果からなる比較結果ベクトルを定義するステップ。もし一次パラメータが閾値より小さいとき、値0とし、もし一次パラメータが対応する閾値より大きいとき、値を1とするように比較結果を定義するステップ。比較結果ベクトルを値0と1の連続で形成するステップ。比較結果ベクトルの成分を組合わせるステップ。特に比較結果ベクトルの成分の加重和を求めるステップと、そして、この結果を包括的閾値と比較するステップ。包括的閾値は、実験又は各一次特徴パラメータに対する閾値の正規化された加重和によって求めることができる。比較結果ベクトルの成分の加重和と包括的閾値との比較結果に関連して、ピクセル又はボクセル又はピクセル又はボクセルのグループが有効に追跡可能なランドマークであるかどうかを決定するステップ。
【0141】
代替的に、上述のように二次パラメータを形成するために一次パラメータがグループ化された場合に、同じステップが二次特徴パラメータに対して適用可能である。この場合、比較結果ベクトルは、その成分が各二次特徴パラメータと対応する閾値との比較結果として定義される。また、ここで数値、特に0と1が、二次パラメータが対応する閾値より高いか低いかという条件にそれぞれ対応する比較結果を表現するために選択可能である。次に、ピクセル又はボクセル又はピクセル又はボクセルのグループが追跡すべき良いランドマークであるか否かということに関連するピクセル又はボクセル又はピクセル又はボクセルのグループの有効性が、比較結果ベクトルの成分の加重和と包括的閾値との比較によって求められる。
【0142】
特徴パラメータは、前記特徴パラメータの種類又はその数学的意味に依存した異なる方法によっても評価可能である。
【0143】
従って、例えば一次特徴パラメータのリストに含まれるヘッシアン行列の固有値を考慮すると、前記ヘッシアン行列の固有値の値から、ピクセル又はボクセルの重要性を評価することができる。
【0144】
ウェーブレット変換係数を考慮することによっても、画像のピクセル又はボクセルの重要性を評価することができる。この技術は、文献「Wavelet-based Salient Points: Applications to Using Color and Texture Features, E. Loupias, N. Sebe, Visual' 00, pp. 223-232, Lyon, France」と「Image Retrieval using Wavelet-based Salient Points, Q. Tian, N. Sebe, M.S. Lew, E. Loupias, T.S. Huang, Journal of Electronic Imaging, Vol. 10, No.4, pp. 835-849, October, 2001」とにより詳細に開示されている。
【0145】
従って、上述の開示内容に加えて、特徴がピクセル又はボクセルのものかを評価するための閾値が、追跡すべき特徴又はランドマーク即ち画像のピクセル又はボクセルの一つに対応し、幾つかの特徴パラメータは、ヘッシアン行列の固有値とウェーブレット変換係数に関する代替的な方法によって評価することができる。
【0146】
図21に、有効に追跡可能なランドマーク又は有効に追跡可能なランドマークでないものの質を評価する別の実施例を示す。ここで、いわゆる分類アルゴリズム、例えばクラスタ化アルゴリズム又は予測アルゴリズム、例えば人工の神経回路網又は統計アルゴリズムが用いられる。
【0147】
この実施例を実行可能にするために、第一に、有効に追跡可能なランドマークか否かが知られたランドマークに対応したピクセルのデータベースを構築する必要がある。ここで、i=1, ..., n の各ランドマークLiに対し、成分が前記ランドマークLiが特定された画像のピクセル又はボクセルの前記方法ステップによって求められた一次又は二次特徴パラメータであるベクトルが構築される。
【0148】
ベクトルの成分は、更にそれぞれ有効に追跡可能なランドマークか否かの質を示す二つの異なる値を推測可能な有効な変数を備える。
【0149】
前記データベースは、図21に概略的に示される。ここで、Liはランドマーク、Viは対応する特性のベクトルであり、「有効性」は、ランドマークLi又はそれに一致するピクセルが有効に追跡可能か否かを表現するパラメータである。データベースのレコードは、通常の手段で分類アルゴリズムを鍛えるために用いられる。ひとたびアルゴリズムが鍛えられると、ピクセルI(m, n)の特徴パラメータは、それによって処理され、そして分類アルゴリズムは出力として有効な変数の値と、ピクセルI(m, n)が有効に追跡可能なランドマーク又は有効でない追跡可能なランドマークと一致するか否かの予測又は決定を与える。
【0150】
本実施例及び図20の一つ前の実施例に従う評価プロセスの組み合わせも適用可能であることを強調しておくことは重要である。この場合、一次特徴パラメータを用いる代わりに、目標ピクセル又はボクセルと、鍛えられたデータベースのランドマークのための入力データベクトルを構築するために二次特徴パラメータのみが用いられる。
【0151】
上述の異なる評価ステップに対して、特に有効に追跡可能であると識別された二つのランドマークが互いに一定の距離にある場合に、別なランドマークを確認する又は廃棄するステップが実行可能であることも考慮されるべきである。これは、前記二つのランドマークの一次特徴パラメータの少なくとも一つ、又は一部又は全部、又は二次特徴パラメータの少なくとも一つ、又は一部又は全部を比較することによって、そして、選択関数を二つのランドマーク間の選択基準を表す比較結果に適用することによって、実行可能である。
【0152】
以下において、自動ランドマーク識別の例を図11ないし14を用いて示す。これらの図式的な例は、理解を容易にするために、二次元画像を用いている。三次元画像への方法の一般化は、ここでは簡潔に記述された数学的形式主義によって、概要が与えられる。
【0153】
三次元と、三次元ボリュームデータを表現する断面画像セットにおける自動特徴選択作業は、Lucas & Kanadeの二次元特徴検出アルゴリズムの実質的拡張である。
【0154】
時間の僅かに異なる間に記録された同一患者の二枚の画像ボリュームを、IとJと定義する。
【0155】
第一の画像ボリュームIにおいて、二つのボリューム間で追跡すべき特徴とみなされた一定数の選択されたボクセルが識別される必要がある。
【0156】
Iをオリジナルの第一画像ボリュームとし、Jを次のボリュームとする。ここで、対応する特徴は見出されていなければならない。
【0157】
第一ステップは、第一画像ボリュームにおける最小固有値の計算で構成される。これは、以下のように実行される。
【0158】
第一画像ボリュームは、三つの勾配ボリュームIx, Iy, Iz(画像ボリュームを表現するデカルト座標系のx-, y-, z方向)に変換され、x、y、zは各シングルボクセルの座標である。
【0159】
勾配ボリュームは、一般的にx方向に(2?x + 1)、y方向に(2?y + 1)、z方向に(2?z + 1)(ここで?x, ?y, ?z = 1, 2, 3, 4, ..., nボクセル)のサイズのいわゆる窓に含まれる隣接するボクセルに対する各ボクセルの強度勾配から算出される。
【0160】
目標ボクセルを中心とする隣接ボクセルの3×3×3アレイのサイズを考えると、3つの勾配ボリュームは、以下のように定義される。
a)x方向の勾配ボリューム
【数1】
b)y方向の勾配ボリューム
【数2】
c)z方向の勾配ボリューム:
【数3】
次のステップとして、各ボクセルに対し、三つすべての先に算出した勾配ボリュームを用いて、勾配行列Gが設定される。勾配行列は、以下のように生成される。
【数4】
各ボクセルの各勾配行列に対して、以下のように勾配行列から最小固有値?minが算出される。
定義:Graphics Gems IV, Paul S. Heckbert, 1994, Academic Press, S 193-98
【数5】
ここで
【数6】
なぜなら、Gは中心モーメントの行列だからである。
勾配行列Gの固有値は、次式で算出される。
【数7】
次に、固有値?mの最小値に対する閾値が定義される。
【0161】
次に、ボクセルが追跡可能な特徴を表現しているか否かを検討するために、以下の基準が適用される。
1.もし?mが閾値より大きい場合、ボクセルの実際の位置が、追跡すべき特徴のリストに 保存される。
2.もし?mが固有値の全ての最小値?mの最大値のパーセンテージより小さい場合、追跡 すべき特徴のリストから落とされる。
3.もし実際の特徴周りの定義された領域(調整可能、例えば3×3×3)において、他 の特徴が存在する場合、固有値のより小さな最小値?mを有する特徴が、追跡すべき特 徴のリストから落とされる。
4.もし特徴位置周りの3次元隣接ボクセルの平均信号強度(例えば? = 1)が、ゼロに近 い非常に低い値(例えば10)より小さいか等しい場合、その特徴は、周囲の空気の 中で乳房の外側にあるとみなされ、追跡すべき特徴のリストから廃棄される。
5.最後に、特徴の定義可能な最大数のみがキープされる。もしそれ以上の特徴が追跡す べき特徴のリストに存在する場合、より小さな最小値?mを有する特徴が、そのリスト から落とされる。
【0162】
これらのステップは、二次元の場合について図11ないし14で視覚的に説明されている。この場合、ピクセルのアレイが小さな四角のアレイとして示されている。第一の基準を満足する固有値の最小値を有する選択された特徴の位置に対応した異なるピクセルが示されている。これらの特徴と位置は、図3において対応するピクセルに異なる外観を与えることによって強調表示されている。図4において、5×5隣接部が、その5×5隣接部の内側を刻む円によって強調表示されている。
【0163】
図11ないし18の例において、グリッドはボクセルアレイを表現している。グリッドの各要素は、画像を形成するボクセルである。
【0164】
図11の例において、4種類のボクセルがランドマークとして区別され、対応する勾配行列の最小固有値に従って分類されている。102で示されたボクセルは、最大の固有値を有し、202で示されたボクセルは、大きな固有値を有し、302で示されたボクセルは、小さな固有値を有し、402で示されたボクセルは、最小の固有値を有する。
【0165】
図12において、各ボクセルに対して、直径5ピクセル/ボクセルのユークリッド距離内の隣接部が特徴として選択された中央ピクセルの周りに定義される。前記隣接部選択は、図4に円で示されている。
【0166】
図12において、このイベントの第一のケースが、502’と記されたより太い円で強調表示されている。ここで、特徴として選択された二つのピクセルが、ピクセルのユークリッド距離として上で定義された同一の隣接部内に存在する。これらの二つのピクセルは202’と302’で示され、上記定義に従って、特徴202’として選択されたピクセルは、第二のピクセル302’より大きな固有値を備えた勾配行列を有する。よって、この第二のピクセル302’は廃棄される。
【0167】
図13には次のステップを示し、図11で特徴として選択されたピクセル302’はもはや存在しない。図13において、隣接部を示す円を備えたピクセル202’’が設けられている。図12と同様の状況がここでも存在し、ピクセル202’’の隣接部に、ピクセル202’’より小さな固有値を有する勾配行列を持ったピクセル302’’が存在する。
【0168】
従って、この場合においても、特徴として選択されたピクセル302’’は追跡すべき特徴のリストから廃棄される。
【0169】
最後に、図14のピクセルアレイにおいて、元の特徴として選択されたピクセルの幾つかのみが存在している。ここで、単一の隣接部内において、一つの特徴以上は見られなくなっている。従って、追跡可能な特徴は、画像/ボリューム中に散開している。
【0170】
有効に追跡可能なピクセルを定義するための数学的形式化と基準と同様に、図11ないし14に示された本例は、目標ピクセルと選択された窓のピクセルとの強度信号に存し、また前記窓の画像データの勾配行列の固有値に存する特徴パラメータのセットに限られる。
【0171】
本発明の次のステップによれば、ピクセルの強度はさらに調査され、そして、所定値より低い強度を有する特徴として選択されたピクセルも、追跡すべき特徴に対応するピクセルのリストから廃棄される。ほとんどゼロに近い非常に低い閾値が、周囲空気内で見つけられた特性に用いられる。このステップは、図示しないが、当業者には自明である。
【0172】
ひとたび2次元又は3次元画像の特徴、即ち有効に追跡可能なランドマークとみなされたピクセル又はボクセルが求められると、前記特徴の追跡が実行される。本発明の好ましい実施例によれば、これはいわゆる「Lucas Kanade特徴追跡のピラミッド形の実施」を用いて実行される。詳細な記載は、Jean Yves Bouguetが記した記事「Pyramidal Implementation of the Lucas Kanade Feature Tracker, Description of the Algorithm」Intel Corporation Microprocessor Research Labs.発行、により与えられる。
【0173】
上述の出版物によって二次元の例について開示された理論を、以下で3次元の場合に更に適用する。上記方法とアルゴリズムとによって個別化された追跡すべき特徴を表わす各ボクセルに対し、以下のステップを実行する必要がある。
【0174】
ボリュームIにおけるポイントとしてuを、ボリュームJにおける対応するポイントとしてvを定義する。第一ステップにおいて、レベルL 0, ..., mにおけるボリュームを表現するILとJLとを備え、元のボリュームを表わす基底部(L = 0)を備えたボリュームIとJに対応する二つのピラミッドが設けられる必要がある。各以下のフロアのボリュームは、一つの平均値に対する直接の隣接部における一定量のピクセルを組合わせることによって寸法が削減される。追加のフロアの数は、ボリュームのサイズに依存する。
【0175】
4.2.以下のステップにおいて、いわゆる包括的ピラミッド形動きベクトルgが最上階で初期化される。
【数8】
全てのレベルLに対して、結果として生じる動きベクトルdLが次に計算される。
これは、以下のステップに従って実行される。
a)各ボリュームILに対して、ポイントuLの位置が次式によって計算される。
【数9】
ここでpは実際のボリューム位置である。
b)ボリュームのピラミッド形の表現の各レベルLにおいて、デカルト座標x、y、zの各方向において勾配が計算される。デカルト座標xに対し、ボリューム勾配が以下の式によって計算される。この等式は、他の二つのデカルト座標y、zに対しても適用される。
【数10】
c)次に、ボリューム勾配が以下の式による勾配行列を計算するために用いられる。
【数11】
ここで値?は、エリアサイズ又はボクセルに影響を与える隣接部を定義する。行列の要素は、上述の定義によるものと同様にして得られる。即ち、
【数12】
d)各レベルLに対して、反復ベクトル?が初期化される。
【数13】
e)次に、k=1から最大カウントKまで、又は最小置換
【数14】
まで、以下のボリューム差を計算する。
【数15】
以下の等式に従って、不整合なベクトルを計算する。
【数16】
次に、次式によってオプティカルフローベクトルが定義される。
【数17】
上で定義したオプティカルフローベクトルを用いて、反復動きベクトルが計算可能である。
【数18】
そしてこれをレベルLに対する最後のオプティカルフローと等しく設定する。
【数19】
上述のステップを各レベルに対し、最後のレベルL=0になるまで繰り返す。それによって、次のレベルL−1に対して、包括的オプティカルフローが上記計算値から算出される。
【数20】
f)最終的なオプティカルフローベクトルは、次式となる。
【数21】
g)次に、ボリュームIのポイントuに対応するボリュームJのポイントvの座標は、以下のように計算できる。
v = u + d
【0176】
図15ないし18は、明らかに実際のものの簡単表現である上述の追跡方法の効果の図による概略的な発想を与えようとするものである。
【0177】
図15ないし18において、異なる時間T0とT1に撮影された二枚の画像IとJが、図11ないし14の例のように、ピクセルの二次元アレイとして表現されている。
【0178】
追跡すべき特徴に対応し、上述の方法によって選択されたピクセル20が、画像Iの二次元ピクセルアレイにおいて示されている。画像Jにおいて同一の位置のピクセル20’’が、画像化された対象この場合乳房の動きが無い場合におけるピクセル20に対応している。ピクセル20’’の一定の隣接ピクセルが、ピクセル20’’を中心とした7×7ピクセルアレイに対応した灰色の影つきの四角で示されており、隣接部は番号10が付されている。
【0179】
画像Iの特徴20の新しい位置が20’で示されている。これは、ピクセル20の位置における乳房組織が、ピクセル20’の位置に時間とともに移動したことを意味している。明らかに、前記ピクセル20’は、追跡が実行される方法を説明するためのみに示されている。実際の場合において、このピクセルは未知であり、さもなければ追跡は不必要である。
【0180】
既に上で述べた数学的形式による選択された特徴に対応するピクセルを追跡する一般的な方法について、以下に図で与えれ、簡単のために二次元の場合に限定された実施例によって記載する。しかし、追跡の一般的数学的記述を踏まえた当業者には、それでも3次元の場合について追跡する方法を理解し、評価することができる。
【0181】
追跡の初めに、画像I内のピクセル20が、同じ位置、即ち画像Jのピクセル20’’にあるものと仮定する。このピクセル20’’を、以下では「追跡されたピクセル」と呼び、一方で特徴として選択された画像Iのピクセル20を「特徴ピクセル」と呼ぶこととする。まず追跡されたピクセル20’’の隣接部に対応した一定の領域の反転勾配行列にアルゴリズムを適用し、本例においては画像Jのピクセル20’’を中心とした前記7×7ピクセルアレイ(この値は調整可能)が計算される。この隣接部が、いわゆる追跡領域に相当する。Lucas & Kanadeの特徴追跡アルゴリズムのピラミッド形の実施を適用することで、これは3次元画像ボリュームに対しても適応されて、詳細は後述するが、画像J上の追跡された特徴の新しい位置を見つけることができる。
【0182】
画像J内の追跡された特徴20’の新しい位置は、オプティカルフローベクトルを決定する。これは、第一の反復ステップにおいて計算されて定義される。図15と16の比較から明らかなように、追跡領域は、追跡領域の中心が同じ組織領域20’の新しい位置により近くなるようにシフトされる。
【0183】
追跡は、画像I内の特徴を表わすピクセル20と画像J内の追跡領域の中心にある新しい識別された位置との間の不整合ベクトルの計算を繰り返すことによって、さらに反復的して適用される。これは、新しいオプティカルフローベクトルを定義することになり、図17に示すように、ピクセル隣接部又は追跡領域の新たなシフトを生じさせる。
【0184】
画像I内の特徴ピクセル20と追跡領域の中心間の不整合ベクトルを求めるステップは、不整合ベクトルが調整可能な最小値に達するまで続けられる。
【0185】
この状況を図18に示す。図15ないし17の図的表現に関連して、ピクセル20’は追跡領域がシフトするという事実を強調するために位置を変えておらず、追跡領域は最終的に画像Iのピクセル20としての乳房の同じ位置を表現する正しい新しい位置20’を中心とする位置に動いている。
【0186】
上記例は、ピクセルによって表現されたただ一つの特徴に対して実行されている。上記ステップは、画像内で選択された各特徴に対して実行する必要がある。
【0187】
すべての特徴がその新しい位置にまで追跡された後に、画像/ボリューム全体にわたって、ただし乳房組織を表現する領域内でのみ、異なる位置にあるオプティカルフローベクトルの一定量が存在する。他の登録アルゴリズムと比較した最大の長所は、画像/ボリューム全体にわたってオプティカルフローベクトルが散会することである。次に、逆オプティカルフローベクトルを画像/ボリュームJ内の各特徴と、その周囲の隣接部とに適用することによって、モーフィングが実行される。これは、新しい仮想画像/ボリュームJ'をもたらす。
【0188】
図19に、実際の診断画像に対する前記方法の効果を示す。左側の第一の画像セットは、造影剤投与無しでの時刻t0における乳房検査からの代表的なスライスを示す。左側の第二画像は、造影剤投与後の時刻t1における同じスライス位置を撮影したものである。左側の第三画像は、時刻t1のピクセルの信号強度から時刻t0のものをピクセル毎に引き算した結果を視覚化したものであり、いわゆる引き算画像である。時刻t1とt0間で乳房の相対的移動が無い場合、造影剤強調が生じない中で、結果の画像では黒色を示すゼロ信号強度のピクセルが多数を占める。時刻t1とt0間での移動は、画像/ボリュームIとJの同じ位置のピクセル単位でアーチファクトが発生するという効果を有し、これは組織の同じ領域には対応していない。本例においては、引き算画像は小さな矢印で示されたコントラスト強調された腫瘍を含む乳房を表わしている。乳房のほとんどの種類の悪性および良性腫瘍は、増加した血管新生によって生じる追加的な造影剤摂取を示す。両乳房の中央の白く残った構造は、移動と組み合わされた正常な乳房柔組織の弱い造影剤摂取の結果である。
【0189】
図19の右側は、前記登録アルゴリズムが適用された対応する画像を示す。
【0190】
時刻t0における画像はアルゴリズムによって修正されておらず、従って同じものを留めている。しかし、三次元特徴追跡後は、新たな仮想画像ボリュームが算出される。時刻t1’での右側の第二画像は、登録適用後の同じスライス位置の対応する仮想画像を示す。右側の最後の画像は、時刻t1’から時刻t0でのピクセルの信号強度を引く、ピクセル毎の引き算結果を視覚化したものであり、新たな(仮想)引き算画像である。一つの乳房は、小さな矢印で示されたコントラスト強調された腫瘍を表わしている。正常な乳房柔組織を表わす乳房内の領域は、幾分抑制された造影剤摂取を示し、動きアーチファクトの減少により、僅かに少ない輝きを見せている。両乳房上部の白い構造の大部分は、動きアーチファクトの補償によって消失している。
【0191】
一枚は造影剤投与前、一枚は投与後の二枚の連続した画像のみが検討されているので、図19の例は単純化されている。
【0192】
実際には、一連の断面画像(例えば50枚以上)が、造影剤の投与前後の一定時間毎に撮影され、引き算画像が連続した画像の各二枚間で算出される。
【0193】
1組の中に一枚以上の断面画像を確保するのに加えて、2組以上(例えば5組以上)の画像を次々に確保することもできる。これは、複数の引き算画像/ボリュームをもたらす。もちろん、開示したアルゴリズムは各組のすべての組み合わせに適用できるので、造影剤投与の前か後かに撮影されたにかかわらず、その間に生じた動きを等しく削除することができる。
【0194】
本発明に係る上述の方法は、コンピュータで実行されるプログラムとしてコード化されたアルゴリズムによって実行可能である。
【0195】
このプログラムは、一枚以上のフレキシブルディスクや、一枚以上の光磁気ディスク、一枚以上のCD−ROM又はDVD−ROMその他携帯用記憶装置や、安定したハードディスク等に保存できる。代替的に、ペンドライブ等の固体電子装置にも保存できる。
【図面の簡単な説明】
【0196】
【図1】胸部MRIで使用される乳房固定装置の概要図である。
【図2】引き算プロセスと登録ステップがない場合に生じる動きアーチファクトを示す。上部左右には、同じMRスキャナ位置の乳房の二枚の断面画像を示す。第二画像は第一画像の約2.5分後の造影剤投与後に測定した図である。両画像は、断面画像の全体的な3次元セットの代表的な断面図である。
【図3】本発明に係る登録ステップを示す概要図である。
【図4】10×10画素アレイで形成されたデジタルの又はデジタル処理された画像の概要図である。
【図5】図4の画像に対応する画像データアレイの概要図であり、そこにおいて各画素の視覚的外観、例えばグレースケールデジタル画像の画素の強度は数値で表現される。
【図6】ピクセルの外観、特に選択された窓のピクセルの強度(図5に示す)を表す数値によって目標ピクセルをコード化するために用いられた最先端技術に従うベクトルの第一実施例を示す概要図である。
【図7】図6に従うピクセルコード化方法の概要図を示し、そこにおいてベクトルの成分はより一般化された方法で表現されている。
【図8】座標m、nの目標ピクセルを中心とする一般的な3×3窓と、一般的な勾配行列と、数値的な偏導関数の計算の一つの方法を示す図であり、例えばSobel演算子やFrei-Chen演算子(Digital Image Processing, W. K. Pratt, John Wiley & Sons, Inc. New York, NY, USA, 1991)のような、より複雑な演算子を用いても導関数を計算することができる。
【図9】一般的なヘッシアン行列と数値的二次偏導関数の一つの計算方法を示す。
【図10】本発明に係る目標ピクセルを表現するベクトルの概要図である。
【図11】有効に追跡可能なランドマークの自動測定法の一例の概要図である。
【図12】有効に追跡可能なランドマークの自動測定法の一例の概要図である。
【図13】有効に追跡可能なランドマークの自動測定法の一例の概要図である。
【図14】有効に追跡可能なランドマークの自動測定法の一例の概要図である。
【図15】追跡がどのように実行されるかを示す概要図である。
【図16】追跡がどのように実行されるかを示す概要図である。
【図17】追跡がどのように実行されるかを示す概要図である。
【図18】追跡がどのように実行されるかを示す概要図である。
【図19】造影剤投与前後の異なる時間に撮影された二枚の画像と、特徴追跡と画像登録を実行した後の、両方の画像の引き算で得られた画像とを、図の左側に示した図2の例と比較して示すことによる本発明の結果を表した図である。図2から図19までのすべての画像は、登録以外の他の処理は一切しないで同一のMR試験で撮影されたものである。
【図20】対応する特徴パラメータによって、ピクセルの有効に又は非有効に追跡可能な特徴の質を測定するためのステップからなる本発明にかかる方法の第一実施例を示す図である。
【図21】対応する特徴パラメータによって、ピクセルの有効に又は非有効に追跡可能な特徴の質を測定するためのステップからなる本発明にかかる方法の第二の実施例を示す図である。
【技術分野】
【0001】
本発明は、
a)ピクセル又はボクセルからなる2次元又は3次元アレイで構成された、少なくとも第一と第二デジタル画像又は同じ対象の一連の断面画像を提供するステップと、
b)第一画像又は画像セットの中に、ランドマークとして設定された一定数のピクセル又はボクセルを選択し、追跡すべき前記特徴のリストを生成することで、一定数のいわゆる特徴であるランドマークを定義するステップと、
c)第一画像から第二画像又は、第一から第二の画像位置間のオプティカルフローベクトルによって決定された後の時点で得られた画像又は画像セット又は、特徴として選択された各ピクセル又はボクセルの後の時点で得られた画像又は画像セットから、特徴として選択された各ピクセル又はボクセルの位置を追跡するステップと、
d)第一と第二画像又は、逆オプティカルフローベクトルを第二画像又は画像セットのピクセル又はボクセルの位置に適用することによって後の時点で得られた画像又は画像セットを登録するステップと、
からなる対象の移動によって生じる画像アーチファクトを低減する生物医学画像の登録方法に関する。
【背景技術】
【0002】
この種類の時間系列の一部としての画像登録方法は、例えば米国特許6、553、152によって知られている。
【0003】
米国特許6、553、152は、視覚的ランドマーク識別と、人間オペレータによる二つの画像の対応する画素のマーキングとによって、二つの画像に対する画像登録方法を提供している。
【0004】
この方法は、ランドマークの選択を完全にオペレーターの技能にゆだねる。そして、顕著な又は一意の認識可能な構造が画像内に存在していなければ、実際に実行することは非常に困難である。例えば磁気共鳴映像(MRI)や超音波やレントゲン写真画像といった診断画像を考えると、画像の構造的特徴によってランドマークの視覚的識別をすることは、極めて困難であり、大きなエラーが生ずる可能性がある。
【0005】
種々の異なるアルゴリズムが輪郭を識別するために試みられており、例えばFischer, B.;Modersitzki, J. Curvature Based Registration with Applications to MRMammography; Springer-Verlag Berlin Heidelberg: ICCS 2002, LNCS 2331, 2002; pp 202- 206に開示されている。また、2回以上測定されたか又は異なるモダリティで調査された組織/領域の放射線画像内のランドマークの識別に関しては、例えばShi J, Tomasi C. Good Features to Track. 1994 IEEE Conference on Computer Vision and Pattern Recognition (CVPR'94) 1994; 593- 600や、Tommasini T, Fusiello A, Trucco E, Roberto V. Making Good Features Track Better Proc. IEEE Int. Conf. on Computer Vision and Pattern Recognition (CVPR '98) 1998; 178- 183や、Ruiter, N. V.; Muller, T. 0.; Stotzka, R.; Kaiser, W. A. Elastic registration of x ray mammograms and three-dimensional magnetic resonance imaging data. J Digit Imaging 2001, 14, 52-55 が開示されている。
【0006】
異なる時間に撮影された同一対象の二枚の平面画像を登録して、時間に対する対称のう動きを計算することは、二次元登録と呼ばれる。同様のことを画像の三次元セット、例えばMRI又はコンピュータ断層撮影(CT)スキャナの断面画像セットに対して行うアルゴリズムは、三次元と称される。ランドマークや対象の移動は、これら画像セットのあらゆる方向に生じる場合がある。
【0007】
比較される画像の異なるモダリティの数に依存して、アルゴリズムは単一様式の、そして複数様式の登録アルゴリズムに分割される。
【0008】
二枚の連続したMRI画像を比較することは、単一様式の登録である。
【0009】
通常、単一様式登録アルゴリズムは、両画像内の明確なランドマークや複雑な形状を識別しようとし、ほとんどの場合、二次元上で動作する。二番目のステップにおいて、両画像の対応するランドマークの位置が比較され、それらの運動ベクトルが計算される。この情報は、動きアーチファクトを除去するために、二番目の画像の全ての画素を新しい位置にシフトバックするために用いられる。
【0010】
「剛体レジストレーション」は、単一ユニットとしての二次元又は三次元画像/ボリュームを一定の方向に移動させる。「非剛体レジストレーション」は、一定のエリア/ボリュームの単一又は複数の画素を異なる方向へ移動させる。
【0011】
例えば非常に弾力性のある乳房組織に関し、動きアーチファクトを除去するために非剛体レジストレーションが必要である。
【0012】
骨のような一定の構造をとらない特別な柔らかさゆえに、乳房に直接隣接した部分は、異なる方向に移動可能である。従って、どの登録アルゴリズムも、他の既存のアプローチも、乳房組織全体に広がるランドマークの識別が必要である。
【0013】
アルゴリズムの一つのグループは、限定された数の発見すべきランドマークを定義する。これらは、その数価によって分類される。乳房のMRIイメージに関し、MRIにおいて脂肪と空気の間に高いコントラストがあるため、これらのアプローチはしばしば乳房の外側境界や乳房構造において多数のランドマークを見つける(Huwer, S.; Rahmel, J.; Wangenheim, A. v. Data-Driven Registration for Lacal Deformations. Pattern Recognition Letters 1996, 17, 951-957)。
【0014】
圧縮パネルを用いるか、少なくとも両乳房を腹臥位に横たえて固定することで、通常圧縮装置に直接接する皮膚のすべての部分は良く固定される(図1)。対照的に、心臓の鼓動や呼吸によって、両乳房の内部領域は僅かに動き続ける。外側の皮膚と空気の境界と比較して、脂肪と柔組織間のコントラストは非常に低いので、乳房の中心にある有効なランドマークの数は、これらのアルゴリズムにおいて低いままである。
【0015】
ランドマーク選択における一般的な問題は、選択したランドマークが非常に少ないことが、不十分な又は不正確な登録につながるという事実に存する。しかし、人間の観察者もプログラムも、例えば乳房中央の異なる組織タイプのより低いコントラストによって、より再割付けが不確実であるランドマークを含めざるを得ないので、追加のランドマークを選択することは、必ずしも正確な登録を保証しない。ランドマークの数を増やすことは、常に登録の計算の複雑さを著しく増加させる。
【0016】
いわゆる非剛体レジストレーションに特に適した特徴選択と特徴追跡アルゴリズムは、いわゆるLucas & Kanadeアルゴリズムと、Lucas BD, Kanade T. An Iterative Stereo Vision, IJCAI81; 674-679 1981 及び付録1として下に同封したBouguet J-Y. Pyramidal Implementation of the Lucas Kanade Feature Tracker, TechReport as Part of Open Computer Vision Library Documentation, Intel Corporationに詳細が開示された、特にピラミッド形の実施である。Lucas & Kanade技術とそのピラミッド形の特徴追跡実施は、画像内の信頼できるランドマークの自動認識や、異なる時間に撮影された二枚の画像間の追跡に、特に適している。選択された特徴の置換又は移動は、シフトとして又はいわゆるオプティカルフローベクトルとして求められる。この技術は、例えばロボット工学におけるコンピュータ・ビジョンに広く適用され、デジタル画像処理といくつかの基礎コースで教えられるコンピュータビジョンにおける当業者の一般的な知識の一部である。
【発明の開示】
【発明が解決しようとする課題】
【0017】
しかしLucas & Kanade技術を実用的に特に乳房のデジタルの又はデジタル処理された診断画像に適用することを考慮すると、前記技術は興味の範囲外も含めた非常に多くの特徴を識別し、例えば乳房まわりの空気のノイズの影響を受ける。
【0018】
以下の例は、上記問題をより詳細に明らかにする。
【0019】
コントラストを強調した(CE)磁気共鳴乳房X線撮影(MRM)は、浸潤性腫瘍に対する鋭敏な感度によって有用で重要な検査である。それは乳房柔組織の密度と独立な高い陰性適中率を示す。最近のマルチセンターな研究は、画像解析技術に従って、感度と特異性が異なることを示す。1.5Tスキャナを用いて、96、93、又は86%の感度と30、50、又は70%の特異性がそれぞれ計算された。従って、MRMの主な欠点は、偽の確定診断が高い割合で残存していることである。
【0020】
通常MRイメージングは、患者が腹臥位の状態で、製造業者の二重乳房コイル(米国の場合、一重乳房コイル装置も多い)を用いて1.0又は1.5テスラの撮像装置で実施される。このプロトコルはあらゆる方向において30から180秒の間の収集時間で、少なくともダイナミックコントラスト強調されたT1加重したシーケンスを含む。gadopentetate dimeglumine又は他の常磁性MR造影剤のボーラス注入の前に、少なくとも一回測定が実施され、その後に数回の測定が実施される。
【0021】
乳房組織の領域への実際の造影剤の摂取を計算するために、造影剤の使用前の各ボクセルの信号強度は、造影剤の使用後の信号強度から差し引かれなければならない。MR画像のスライス厚さによる部分容積効果と同様に、呼吸と心臓の鼓動によって生じる両乳房の最低限の動きによって、異なる時間に撮影された同じ画像位置のボクセルは、乳房組織の正確に同じボリュームを示さない。従って(腫瘍の場合を除いて)引き算した画像は完全に黒くはならない。小さな移動の効果は、乳房の外側境界で最も良く現れる。脂肪組織の信号強度は大きく、周囲の空気の信号強度はほとんどゼロであることから、非常に小さな移動が脂肪を先の画像の空気の位置に置換する。脂肪の大きな信号強度から空気のゼロの信号強度を引き算することで、CM用途における大きな信号上昇をシミュレートする。結果として、引き算画像が、画像化された乳房の少なくとも部分まわりの強い造影剤の摂取を表す太く白い線を示す(図2−乳房の白い外側境界が表現されている)
【0022】
デカルト座標系によって定義された三次元又は容積測定用の画像収集を考慮すると、xとy方向のあらゆる移動に加えて、(胸に向かっての)z方向の幾らかの移動が常に見られる。これは、一連の検査中の患者の弛緩によって生じる。
【0023】
アーチファクトを低減するために、乳房は通常何らかの圧縮方法(製造業者に依存する)によって乳房コイルの内側に固定される。この事実にもかかわらず、最小限の移動は常に発生する。従って、乳房を固定しているにもかかわらず、画像には常にx、y、x方向に僅かな動きアーチファクトが見られる。この動きアーチファクトは、明るい画素として引き算画像の中で視覚化される。動きアーチファクトがなければ、腫瘍性組織が占める領域を除き、引き算画像は絶対的に黒くなる。
【0024】
周囲の空気のノイズ入りの信号内で、一定数のランドマークが乳房の外側で検出される。各特徴とみなされるものは追跡しなければならず、コンピュータに不必要な負荷を発生させるという事実の他にも、ノイズは変量効果であり、ノイズの中で見つけられたすべての特徴又はランドマークは、二番目の又はそれに続く画像において対応するランドマークを有さない。しかし実は、ノイズの中から発見されたランドマークに対応する二番目の、又はそれに続く画像又は画像セットにおけるランドマークを調査するアルゴリズムは、一定量のランドマークを第二画像又は画像セット内に再割付けする。これが、画像の不正確な登録につながり、結果に影響を与える。
【0025】
現在知られている方法によれば、「特徴選択」は以下のように実行できる。
B1)第一画像又は第一画像セットの各ピクセル又はボクセル周りに隣接するピクセル又はボクセルを定義する。ここで、前記隣接するピクセル又はボクセルは、限られた数のピクセル又はボクセルからなり、単一画像の場合には二次元隣接部又はボクセルが選択され、断面画像セットの場合には、三次元隣接部又はボクセルが選択される。
B2)各ピクセル又はボクセルに対して、前記隣接部又はボクセルのすべてのピクセル又はボクセルの平均信号強度を定義する。
B3)平均信号強度閾値を定義する。
B4)ステップB2で各隣接部又はボクセルに対して定義した平均信号強度を、前記平均信号強度閾値と比較する。
B5)前記隣接部又はボクセルの平均信号強度がステップB4の閾値より大きい場合に、そのピクセル又はボクセルは追跡すべき特徴として定義され、追跡すべき特徴のリストに追加される。
【0026】
平均信号強度閾値は、調整可能であり、変更可能である。経験的又は実験的基準は、例えばテストセットやサンプル画像に対しこの方法を実行し、異なる閾値に対して得られた結果を比較することによって、このような閾値を求めるために適用可能である。
【0027】
自動特徴選択が、例えば二次元画像又は二次元の場合の一般化としての三次元画像に対していわゆるLukas&Kanadeアルゴリズムを応用することによって、現在用いられている。良く知られたLukas&Kanadeアルゴリズムによる自動特徴選択法は、Lucas BD, Kanade T. An Iterative Image Registration Technique with an Application to Stereo Vision, IJCAI81; 674-679 に非常に詳細に記載されている。
【0028】
三次元画像や三次元データボリュームを表す断面画像セットに用いられる自動特徴選択は、Lukas&Kanadeの二次元特徴検出アルゴリズムの実質的な拡張であって、以下により詳細に説明する。
【0029】
本明細書と請求項において、「特徴」という用語は「有効に追跡可能なランドマーク」という用語と等価である。
【0030】
追跡すべきランドマーク又は特徴を自動的に選択するための上記方法を用いて、対象の画像の登録が実行可能である。登録のおかげで、動きの差異、又は動きアーチファクトを補償することができる。しかし、登録された画像シーケンスは、十分な外観を有さない。上記方法で画像を処理した後、画像はシャープさが減ったように見え、又は他の効果が見える。
【0031】
さらに、有効なランドマーク又は特徴の自動選択のための方法は、手作業で、即ち人間オペレータの技能に基づいて実行されるのを避けて、客観的な基準に従って、画像内で追跡されるべき有効なランドマークを選択することを可能にする有用なツールではある。しかし、有効なランドマークの識別は、画像の各ピクセルを目標ピクセルとして考慮し、目標ピクセルの周囲の一定数のピクセルを備えた各目標ピクセルまわりの窓を定義することによって求められた勾配行列の固有値の関数であるところの特定の基準にのみ委ねられている。
【0032】
この方法ステップは、一種の境界検出であって、境界の特徴を有するとして識別された目標ピクセルの中で、幾つかのみが登録プロセスにおいて有効に追跡可能なランドマークとして選択される。
【0033】
本発明の目的は、特に二次元又は三次元の生物医学的診断画像に対して画像を登録する方法を提供することに存する。この方法は、二次元及び三次元で動きを追跡可能にすることによって、また特にデジタル画像又はデジタル処理された画像内の対象領域の外側のノイズによるアーチファクトを低減することによって、また隣接する組織タイプの信号強度における高低から独立な、乳房全体にわたりランドマークを広げることによって、既知の登録方法の欠点を効率的に克服している。
【0034】
本発明は、特に画像内の追跡すべき有効なランドマークを選択するステップを向上させることを目的としている。
【課題を解決するための手段】
【0035】
本発明による画像を登録するための方法は、
a)組織又は組織領域又は解剖した領域のMRI,超音波又はレントゲン写真の走査により得られた第一のデジタル又はデジタル処理された画像又は断面画像セットを提供し、
ここで前記画像はピクセル又はボクセルの二次元又は三次元アレイによって形成され、
同じ組織又は組織領域又は同じ解剖した領域のMRI,超音波又はレントゲン写真の走査により得られた少なくとも第二デジタル又はデジタル処理された画像又は同じ解剖した領域の断面画像セットを提供するステップと、
b)ランドマーク又は特徴として設定された一定数のピクセル又はボクセルを選択し、追跡すべき前記特徴のリストを生成することによって、第一画像又は画像セットの中に、一定数のランドマーク、いわゆる特徴を定義するステップと、
c)特徴として選択された各ピクセル又はボクセルに対する第一から第二画像又は画像セットのオプティカルフローベクトルを求めることによって、第一から第二画像又は画像セット間で、特徴としての各ピクセル又はボクセルの位置を追跡するステップと、
d)逆オプティカルフローベクトルを第二画像又は画像セットのピクセル又はボクセルに適用することによって、第一と第二画像又は画像セットを登録するステップと、からなる。
【0036】
ここにおいて、
B1)第一画像又は第一の一連の断面画像の各ピクセル又はボクセル周りに隣接した、限られた数で構成されるピクセル又はボクセルを定義するステップと、
B2)各目標ピクセル又はボクセルに対し、
外観を記載する数値パラメータ、すなわち前記ピクセル又はボクセル窓の前記目標ピクセル又はボクセルとピクセル又はボクセルの部分とからなるいわゆる数値外観パラメータの関数として計算され、又は、
前記窓のピクセル又はボクセルを表現する数値パラメータのマトリックス又は数値パラメータの前記マトリックスの変形で構成される一つ以上の特徴パラメータの関数として計算される、一つ以上の特徴のあるパラメータを定義するステップと、
B3)前記目標ピクセル又はボクセルの特徴的なパラメータの関数としての有効に追跡可能なランドマーク又は特徴に存するピクセル又はボクセルを定義するステップと、
からなる自動追跡可能なランドマーク選択ステップが実行されることを特徴とする。
【0037】
本発明の第一実施例によれば、目標ピクセル又はボクセルが追跡すべき有効なランドマーク又は特徴であるかを求めるための関数は、それぞれの値又は目標ピクセル又はボクセルの特徴パラメータと閾値との比較に存する。
【0038】
閾値は、実験によって求められ、テストセット又はサンプル画像に対して方法を実行し、異なる閾値に対して得られた結果を比較して求められる。
【0039】
目標ピクセル又はボクセルが追跡すべき有効なランドマーク又は特徴であるかを求めるための代替的な方法は、分類アルゴリズムの入力ベクトルであるところのベクトル形式の目標ピクセル又はボクセルの特徴パラメータを組織化することに存する。
【0040】
コンポーネントとして目標ピクセル又はボクセルの一定数の特徴パラメータを備えたベクトルを処理するために、前記分類アルゴリズムは、前記タスクを実行するために具体的に訓練される必要がある。
【0041】
従って、上記代替方法による本発明の方法は、
一定数の有効なランドマークが、有効に追跡可能なランドマーク又は特徴として識別された既知の場合の一定数の画像を提供するステップと、
上記ステップB1)とB2)に存する自動追跡可能ランドマーク選択ステップを適用することによって、一定数の画像の中で有効に追跡可能であると識別された前記ランドマークに対応して前記ピクセル又はボクセルに対する特徴パラメータセットを求めるステップと、
有効に追跡可能であると識別された各ランドマークに一意的に関連付けられ、成分として前記有効に追跡可能なランドマークと一致する前記ピクセル又はボクセルの特徴パラメータを備えたベクトルを生成するステップと、
変数の所定の数値によって有効に追跡可能な質を表現し、
前記数値を、有効に追跡可能なランドマークと一致する各ピクセル又はボクセルをコード化するベクトルに関連付けるステップと、
有効に追跡可能なランドマークと一致する各ピクセル又はボクセルをコード化する各ベクトルが、分類アルゴリズムのためのトレーニングデータベースの登録を形成するステップと、
前記データベースによって分類アルゴリズムをトレーニングするステップと、
鍛えられた分類アルゴリズムの入力に前記目標ピクセル又はボクセルの前記特徴パラメータを備えたベクトルを提供することによって、目標ピクセル又はボクセルの有効に追跡可能なランドマークの質を求めるステップと、を備える。
【0042】
改善した本発明によれば、前記トレーニングデータベースは、有効に追跡可能ではないと知られたランドマークと一致する画像のピクセル又はボクセルのベクトル形式に組織化され、有効に追跡可能でないランドマークの質を、追跡可能なランドマークの質を表現する前記変数の値とは異なる変数の所定の数値によって表現する特徴パラメータも備え、一方で有効に追跡可能でないランドマークの質を表現する前記数値を、既知の有効に追跡可能でないランドマークと一致した各ピクセル又はボクセルをコード化するベクトルに加える。
【0043】
分類アルゴリズムを用いた有効に追跡可能なランドマークを求めるためのこの種の自動的方法は、例えば適応度や他の同様な統計指標のような分類の信頼性を示す統計的尺度をも提供するという更なる長所を有する。
【0044】
分類アルゴリズムとして、あらゆる種類のこれらアルゴリズムを用いることができる。分類アルゴリズムの特別な場合としては、いわゆるクラスタ化アルゴリズム又は人工の神経ネットワークがある。
【0045】
目標ターゲットピクセル又はボクセルや、前記ピクセル又はボクセル窓のピクセル又はボクセルのそれぞれ、又は一部の外観を表現する数値外観パラメータの関数として、又は前記窓のピクセル又はボクセルを表現する数値行列又は前記数値行列の転換の一つ以上の特徴パラメータの関数として計算される特徴数値パラメータに関し、多くの異なる関数が多くの異なる組み合わせの中で使用可能である。
【0046】
特徴又は有効に追跡可能なランドマークに対応する目標ピクセル又はボクセルの特徴パラメータを求めるための第一の関数は、前記隣接部又はボクセルのすべてのピクセル又はボクセルの平均強度値に存する。
【0047】
目標ピクセル又はボクセルのもう一つの特徴パラメータは、ステップB1)で定義された前記窓のピクセル又はボクセル、即ち目標ピクセル又はボクセルに新設するピクセル又はボクセルの一定数を表現する数値行列の特徴パラメータに存する。前記特徴パラメータは、前記窓のピクセル又はボクセルの画像データを備えた数値行列の特異値である。
【0048】
代替的に、又は上述の特徴パラメータとの組合わせにおいて、目標ピクセル又はボクセルのもう一つの特徴パラメータが、前記窓のピクセル又はボクセルを表現する数値行列の勾配行列の固有値に存する。
【0049】
目標ピクセル又はボクセルのもう一つの特徴パラメータが、前記窓のピクセル又はボクセルを表現する数値行列のヘッシアン行列の固有値に存する。
【0050】
代替的に又は上記開示されたものとの組合わせで提供された目標ピクセル又はボクセルの特徴パラメータのためのもう一つの選択が、前記窓のピクセル又はボクセルに対応するデータの数値行列のウェーブレット変換係数の一つ以上に存する。
【0051】
この場合、いくつかのウェーブレットを基礎にした関数が、代替的に又は組み合わせとして用いるために選択可能である。
【0052】
ウェーブレット変換のより詳細な記載は、付録1で与えられている。この付録は、インターネット上で利用可能なように公開されており、"Wavelet for Kids, A tutorial introduction"というタイトルで、Duke大学のBrandi Vidakovic とPeter Muellerが著している。この文献では、ウェーブレット理論が要約され、議論され、画像処理への幾つかの応用が開示されている。画像処理でのウェーブレットの実行を開示した章から明らかに、ウェーブレット分解はパラメータを得ることを可能にする。分解の各レベルに対して、ウェーブレット変換は、平均を表す一つの行列と、いわゆる詳細を表す3つの行列とを生成する。上記行列の一つ以上から、例えば行列の要素の平均を取る以外にも、行列の特異値を取ることなどのよって、幾つかのパラメータを抽出することが可能である。これらのパラメータのすべて又は幾つかは、目標ピクセル又はボクセルの特徴パラメータとして用いることができ、また、選択された窓に含まれる周辺ピクセル又はボクセルと関係する各目標ピクセル又はボクセルを表すベクトルの成分を形成するために用いることができる。
【0053】
代替的に又は上記開示されたものとの組合わせで提供された目標ピクセル又はボクセルの特徴パラメータのためのもう一つの選択が、前記窓のピクセル又はボクセルを表す数値行列の自己相関変換係数の一つ以上に存する。
【0054】
画像処理における自己相関は、一般的に画像テクスチャを特徴づけるための道具として用いられ、二つの画像の数学的評価から成る。この二つの画像は、異なる時刻に撮影されるか、又は第一画像から場所を移動した結果としての第二画像として生成される。自己相関は、これら二つの画像間の関係を測定する。この数学的評価は、目標ピクセル又はボクセルのコード化において、考慮されるパラメータの数を減らす可能性を提供する。
【0055】
目標ピクセル又はボクセルのもう一つの特徴パラメータは、前記窓のピクセル又はボクセルを表す前記数値行列の同時生起行列(又はその特異値)に存する。
【0056】
前記同時生起行列は、固定された空間的関係によって分離された一対のピクセルに対する中間調の二次元ヒストグラムである。この行列は、一対のピクセルの共有の確率分布を近似する。一般的にコントラストやエントロピーのような、テクスチャ測定を計算するために用いられる。
【0057】
目標ピクセル又はボクセルの特徴パラメータの上記例によれば、各目標ピクセル又はボクセルは、前記窓のピクセル又はボクセルを表現する数値行列の固有値又は特異値及び/又は、前記窓のピクセル又はボクセルを表現する数値行列の勾配行列又はヘッシアン行列の固有値及び/又は、ウェーブレット変換係数の一つ以上の係数及び/又は、自己相関変換の一つ以上の係数及び/又は、同時生起行列の一つ以上の入力又は特異値及び/又は前記窓のピクセル又はボクセルと前記目標ピクセル又はボクセルの平均強度値、の二つ以上の組み合わせ又はサブコンビネーションによって記載される。
【0058】
既に開示したように、前記特徴パラメータは、前記特徴パラメータのそれぞれに対し、又はその組合わせに対し閾値を定義することによって、画像における目標ピクセル又はボクセルが、有効に追跡可能なランドマークとして考慮可能かどうかを判断するために用いることができる。
【0059】
もう一つの基準は、前記特徴パラメータをベクトル形式に組織化することに存し、そこで各パラメータは前記ベクトルの成分であり、前記ベクトルの係数を求め、この値を前記係数に対する閾値と比較することで、組織化する。
【0060】
代替的に、前記ベクトルは既に上で開示した分類アルゴリズムによって処理することもできる。
【0061】
各ピクセル又はボクセルと、選択されたピクセル又はボクセル窓によって定義された選択された隣接部又はボクセルによって記載された数値間の関係は、単に各ピクセル又はボクセルの外観を表す値に存するオリジナルの数値行列の前記異なる変換の前記特異値及び/又は固有値によって、及び/又は、ウェーブレット変換の一つ以上の係数及び/又は自己相関変換の一つ以上の係数及び/又は同時生起行列の一つ以上の入力又は特異値
によって要約される。そして、各ベクトルにおいて、この関係は強調表示に特に適した、又は選択された窓内の画像のピクセル又はボクセル間の関係の特定の種類に感受性のある数値によって定義される。
【0062】
従って、ピクセル又はボクセルに対する有効に追跡可能なランドマークの質を求めるためのこれら特徴パラメータの一部又は全てを考慮することは、このピクセル又はボクセルと残りの画像との関係に関しより良くより完全に考慮する助けになる。
【0063】
すでに述べたように、上記方法は磁気共鳴画像(MRI)又はX線CTにおける同一の又は異なる種類のスキャナで撮影された二枚の診断画像又は3次元画像ボリュームを登録するのに用いることができる。
【0064】
上述の登録方法は、例えばMRIや超音波画像や特にコントラストを強調した磁気共鳴マンモグラフィー(MRM)といった造影剤強調診断画像との組み合わせによっても提供される。この方法において、乳房組織領域の実際の造影剤摂取を計算するために、造影剤使用前の各ボクセルの信号強度が、造影剤使用後の信号強度から差し引かれる必要がある。このように、画像化された組織に造影剤の存在しない時に撮影された第一画像と、画像化された組織に造影剤が存在しかん流した第二の時間又は後の時間とを登録することは、登録が実行されなかった場合に患者の動きにより引き算画像に現れるアーチファクトを非常に成功裏に抑制するのに有用である。前記二枚の画像の引き算は、造影剤が溜まった組織領域に対応するピクセル又はボクセルの信号強度を高めることに繋がる。造影剤の溜まっていない他の領域は、低い信号強度のピクセル又はボクセルであって、特に暗い灰色か黒に見える。
【発明を実施するための最良の形態】
【0065】
画像の登録は、特に医療診断の分野において、幾つかの種類の調査の実行を可能にする重要なプロセスである。
【0066】
いわゆる画像化された領域のかん流挙動の測定において、登録が重要であることが以下の例により明らかになる。この場合、異なる時刻に撮影された身体上の領域の一連の画像が得られる。最初の画像は、造影剤が身体上の領域にない状態で撮影される。続く画像は、造影剤が身体上の領域に達した状態で撮影される。造影剤は、身体上の領域の周辺組織に対して血管やリンパ管の流れの応答を強調し、この流れがより良く見えるようにする。腫瘍性組織や他の病巣のような病変によって影響された組織は、強化された血管新生によって特徴付けられるので、画像化された組織のかん流挙動の測定は、この病理組織を個性化し、画像においてより良く見えるようにレンダリングするための手段である。画像化された身体上の領域の幾つかの領域における強化された血管新生は、この領域において造影剤が他の領域より集中するという結果をもたらす。また、造影剤は投与した時から血管やリンパ管のシステムによって組織に拡散されるまでに時間を要するので、強化された血管新生は、この強化された血管新生が存在する領域において、造影剤をより迅速に供給し、同時に、画像化された身体上の領域の通常に血管新生された領域に比較して、造影剤の最大ピークに達した後に、造影剤の集中のより迅速な減少を提供する。
【0067】
比較的長い時間間隔に広がった画像の時間系列が得られるので、画像化された身体上の領域の動きは、一枚の画像とそれに続く画像間で生じることができる。この動きは、患者の動きによって、また例えば心臓の鼓動や呼吸といった生理上の理由によって決定される。
【0068】
かん流挙動は、造影剤の画像化された身体上の領域への投与後に得られた時間系列の各画像から前記画像化された身体上の領域に造影剤の存在しない時点での画像を引き算することによって求められる。図2に示すように、造影剤が存在する領域が明るいままであるのに対して、この引き算が実質的に造影剤のかん流の生じない領域のピクセルを黒かほとんど黒色にするという事実のほかに、(上で定義した)各引き算画像に対して平均強度を求め、グラフに一連の画像の各引き算画像の平均強度をプロットすることで、造影剤が検出された領域のかん流挙動の尺度である、いわゆるかん流曲線を引くことができる。
【0069】
診断に関するフォローアップ又は治療が施された際の画像化された領域の病状の進展の
フォローアップを実行するために、画像化が異なる時間間隔で実行された場合に、登録も重要である。ここで、画像を撮影するための装置に対して患者はいつでも位置決めしなおす必要があるので、先のかん流調査の場合よりも、一つの画像とそれに続く画像間の時間間隔はより長くなり、一つの画像とその他の画像間の動きは、よりドラマチックになる。
【0070】
他の場合において、異なる画像化技術で撮影された身体上の領域の画像は登録される必要がある。相関画像又は組み合わせ画像がこの対象に関するより完全で詳細な情報を、又はより良いより明確な組合わせ画像を提供できるように、異なる画像技術はしばしば身体上の領域の画像化された対象の異なる様相や質を強調する又は暴露するので、これは重要である。
【0071】
強固な構造が画像化される場合には、登録は一連の画像の一枚の他に対する単純な平行移動又は回転平行移動とみなすことができる。画像化された対象が柔組織からなるか、柔組織に埋め込まれている場合、画像化された対象や画像化された領域に歪みが生じるので、登録にはこの効果も考慮することが必要である。
【0072】
診断画像の時間系列を登録する一般的なケースは、かん流挙動測定法を適用することで実行される乳房中央を調査する場合である。
【0073】
図1に患者の乳房を固定するための一般的な装置を示す。通常、MR造影は1.0又は1.5テスラの撮像装置に製造業者の二重乳房コイル(米国の場合、単一乳房コイルである場合も多い)を用いて行われ、その間患者はうつ向けになっている。このプロトコルはあらゆる方向において30から180秒の間の収集時間で、少なくともダイナミックコントラスト強調されたT1加重したシーケンスを含む(個々のプロトコルは異なっていてもよいし、追加的に脂肪抑圧を用いることもできる)。0.1 mmol/kgから0.2 mmol/kgの gadopentetate dimeglumine又は他の常磁性MR造影剤のボーラス注入の前に、少なくとも一回測定が実施され、その後に数回の測定が実施される。
【0074】
少なくともMRI画像の最初のボリュームが造影剤が肘正中静脈に注射されて乳房組織に到達する前の時間t0に撮影される。少なくとも画像のもう一つのボリュームが造影剤注射後の時間t1に直接撮影される。通常、時間遅れのある連続したボリュームが撮影される(例えば、造影剤投与後、6回以下の同一ボリュームが連続する)。
【0075】
図2に、断面画像の3次元ボリュームから取り出された2枚の画像を示す。一枚は造影剤投与前の時刻t0に撮影されたものであり、もう一枚は造影剤投与後の時刻t1に撮影されたものである。引き算画像、即ち第二画像から第一画像をピクセル毎に引き算した結果の像も、同様に示す。ここで登録は行われておらず、例えば矢印で示された領域のように、造影剤の摂取されていない領域のピクセルは、輝いたままであることは明らかである。
【0076】
図3に、従来の登録方法とも共通の、登録プロセスの主なステップを示す。
【0077】
立方体で表現された画像ボリュームが異なる時刻t0とt1に画像化される。ここでは三角と円で表現された二つの対象が、ランドマークとして選択される。動きのせいで、画像ボリュームt0での二つの対象の位置は、t1での画像の同じ対象の位置とは一致しない。
【0078】
従って、第一ステップは、画像の中で追跡すべきランドマークを個別化することである。三角と円で表現された追跡すべきランドマークがひとたび個別化されると、いわゆる動きベクトルの測定を可能にする追跡ステップが実行される。
【0079】
これらベクトルによって、画像化された身体構造上の領域の動きを大域的運動として、また画像化された身体構造上の領域の一定の領域の局部的運動として表現する、画像中に拡がるベクトル場を生成することが可能である。逆動きベクトル場を画像に適用すると、第一画像の時刻t0と第二画像の時刻t1との間で生じる動きの効果は、補償される。乳房のかん流試験の例におけるこの登録の効果は、図19で見ることができる。ここで、左側の連続画像は、時刻t0とt1で撮影された乳房画像が登録されていない場合に関し、下側の画像は、図2の一つと対応する引き算画像である。右側の連続画像は、時刻t0とt1に撮影された2枚の画像を示すが、時刻t1に撮影された画像は、登録方法で処理済みのものである。引き算画像に対応する最後の画像において、左側の引き算画像ではピクセルが輝いていた多くの領域が消失し、又は対応するピクセルが完全に黒か灰色の様相を呈していることが明らかに見ることができる。
【0080】
本発明の目的は、以下のステップでもクリティカルな、画像の登録プロセスの第一ステップを改良することである。画像内の信頼できるランドマークを定義することは、より良い登録結果を得るためにクリティカルである。したがって、有効に追跡可能なランドマークとなる可能性の非常に高い画像のピクセル又はボクセルの特徴的性質を求めるためのより精密で完全な方法が本発明の目的である。
【0081】
図4を参照して、単純化のために二次元の場合について、デジタル又はデジタル処理された画像が、n, m = 1, ... , NであるピクセルアレイP(n, m)によって形成される。各ピクセルは、一定の外観を有するモニター又は画面又は印刷された紙上の小さなドットである。ドットの外観は、数値データによって表現可能である。いわゆるグレースケール画像において、各ドットは前記グレースケールの一定のグレーレベルに対応した一定の強度を有する。カラー画像の場合には、ピクセルの外観を数値データで完全に表現するために、より多くのパラメータが通常必要である。カラー画像におけるピクセルの外観を一義的に表現するために、幾つかのシステムが定義されている。そのうち一つのシステムは、いわゆる既知のHSV(色調、彩度、明暗度)又はいわゆるRGB(赤、緑、青)システムである。もしカラー画像を処理した場合、これらの値が対応する数値データに置き換える必要があることが当業者には自明であるため、本単純化した例においては強度I(n, m)のみが示される。
【0082】
従って、図5に示されたn, m = 1, ... , 10の数値データのアレイI(n, m)は、図1のピクセルアレイの概要図としてのn, m = 1, ... , 10のピクセルP(n, m)で構成される画像に対応し、そして、数値データのアレイは行列である。
【0083】
各ピクセルP(n, m)は、例えばグレースケールデジタル画像のピクセルのグレーレベルのように、数値的にその外観を表現した数値データI(n, m)に一義的に関連する。
【0084】
隣接する断面(二次元)画像の一定数によって形成されるボリューム画像として、三次元画像を定義することができる。この場合、画像はいわゆるボクセルで形成され、そしてアレイは画像データの数値行列も三次元であるように、三次元である。実際には、断面画像セットは定義された厚みを有し、従ってこれら画像は、純粋な二次元画像よりもよりスライスされている。この観点から、常にボクセルと称することがより正確である。それにもかかわらず、モニター画面上の断面画像の視覚的表現は、純粋に二次元画像であって、スライスの厚みは考慮されず、この場合、ピクセルが画像ドットの適切な定義である。
【0085】
図5、6、7に示すように、画像の目標ピクセルは、その成分がその目標ピクセル周りのピクセルに関する情報を含んだベクトルによってコード化可能である。通常、前記周辺領域は目標ピクセルを中心とした(2? + 1)( (2? + 1)ピクセルの次元の窓から成る。ここで、?は(? = 1, 2, ..., N)でユーザーによって任意に定義される整数であり、ここでは「ステップ」と称す。
【0086】
周囲のピクセルで形成された窓の上記定義は、勾配xの周辺ピクセルの定義と等価であり、ここでxは整数であって、この整数はステップにおける目標ピクセルから隣接部までの距離を示す。目標ピクセルとして前記窓の中心にあるピクセルを考慮すれば、勾配1の周辺ピクセルを備えた前記窓は、目標ピクセルに直接隣接したピクセルのシェルであり、勾配2の周辺ピクセルは、目標ピクセルの周囲の二番目に近いピクセルのシェルを備え、画像を形成するピクセルアレイの各方向において、勾配1は目標ピクセルから1距離ステップに対応し、勾配2は、目標ピクセルから2距離ステップに対応する。この定義は、ボクセルアレイで形成された3次元画像に対して準用する。
【0087】
前記窓の最小寸法は、3×3寸法のピクセルアレイに存し、その中心ピクセルが目標ピクセルであり、上で勾配として定義されたステップ?は、1に等しい。例えば5×5又は7×7ピクセル窓といった、より大きな窓を選択することもできて、それぞれのステップ又は勾配は、?=2と?=3である。簡単にするために、本例においては、目標ピクセルが中心にある3×3ピクセルアレイに対応した窓を選択する。
【0088】
画像のピクセルを表現する数値データ行列を参照して、図5にこの窓を示す。窓Wは九つのピクセルを備え、そのうち一つが目標ピクセルである。図5に示された窓は、目標ピクセルP(2, 2)を中心とし、ピクセルP(1, 1), P(1, 2), P(1, 3), P(2, 1), P(2, 2), P(2, 3), P(3, 1), P(3, 2), P(3, 3)を備えており、各ピクセルは対応する数値、即ち前記ピクセルの信号強度に関連する数値パラメータI(1, 1), I(1, 2), I(1, 3), I(2, 1), I(2, 2), I(2, 3), I(3, 1), I(3, 2), I(3, 3)で表現される。
【0089】
前記目標ピクセルと画像内の特定の周辺領域との関係が考慮されるように、目標ピクセルP(2, 2)は、窓Wの隣接画素に関する情報も用いてコード化される。前記目標ピクセルP(2, 2)と上で定義した周辺ピクセルとの関係を表すベクトルの成分として、前記ピクセルの強度は、目標ピクセルの強度とともに測定される。
【0090】
前記ベクトルは、図6に図式的に示される。明確に表すために、目標ピクセルを表現する各ベクトルは、九つの成分を有する。
【0091】
図7に、目標ピクセルをコード化するための前記ベクトルの一般的表現を示す。この場合、目標ピクセルはその強度によって表現され、強度I(n, m)を有するピクセルP(n, m)として定義される。
【0092】
図8に図7のベクトルを参照する窓Wと、いわゆる勾配行列におけるその変換を示す。
【0093】
当業者には自明のことであろうが、二方向が定義可能で、対象の窓における各ピクセルに対し、この二方向において勾配が評価可能であるように、ピクセルを表現する数値アレイと、この場合窓W内のピクセルの強度データアレイI(m-1, n-1), I(m-1, n), I(m-1, n+1), I(m, n-1), I(m, n), I(m, n+1), I(m+1, n-1), I(m, n), I(m+1, n+1) とは、二次元オブジェクトである。
【0094】
勾配行列は、固有値?p(p = 1, 2)で、この場合?1? ?2で表現可能なように、対角化可能である。
【0095】
選択された窓Wのピクセルの強度Iのオリジナル行列は、さらに加工可能であり、前記オリジナル行列に対し、いわゆるヘッシアン行列が計算可能である。また、この場合ヘッシアン行列は、その固有値?pによって表現可能である。
【0096】
選択された窓の強度値Iの3×3行列を考慮した場合、通常、前記行列は対角化可能でなく、固有値は上述のように有意義なものではない。これを考慮して、より一般的な場合において、図8に示すように選択された窓のピクセルの強度Iのオリジナル行列は、いわゆる特異値siで表現可能である。
【0097】
選択された窓に対応する強度行列の特異値又は代替的に、勾配行列又はヘッシアン行列の固有値と組合わせて用いることで、対応した目標ピクセルを一義的にコード化するためのベクトルを生成することが可能である。前記ベクトルの成分は、選択された窓のピクセルと、選択された窓のピクセルに対応する強度値の前記行列を処理することで得られる勾配行列又はヘッシアン行列の固有値とに対応した強度値の行列の前記特異値に存する。
【0098】
上で開示したように、図示された例は、選択された窓のピクセルに対応する強度値のオリジナル行列の変換の特別な場合に限定されるが、その他の変換も代替的に又は組み合わせとして適用可能である。例えば、ウェーブレット分解が実行可能であり、目標ピクセルのコード化ベクトルの成分として、それらのすべて又は少なくとも一部に対し、平均値と詳細値が利用可能である。
【0099】
代替的に、又は組合わせにおいて、選択された窓のピクセルに対応する強度値のオリジナル行列の自己相関変換が利用可能であり、目標ピクセルのコード化ベクトルの成分として、それらのすべて又は少なくとも一部に対し、得られたパラメータが利用可能である。
【0100】
ウェーブレット変換も前処理ステップとして利用可能である。得られた数値データの行列は、対応する勾配行列及び/又は対応するヘッシアン行列を用いて処理可能であり、その特異値又は固有値は、代替的に又は組合わせにおいて求めることができる。
【0101】
前記窓のピクセル又はボクセルを表現する行列のいわゆる同時生起変換に存する、ピクセルの窓を表現するオリジナル行列の別な変換が利用可能である。この場合、この同時生起行列の特異値又は固有値が、目標ピクセルを表すパラメータとして利用可能である。
【0102】
一定数の隣接部からなる選択された窓における目標ピクセルを特徴付けるための幾つかのパラメータの上記記載は、単に単純化のために二次元の場合に限定されているが、上述の選択は、三次元の場合にも適用可能であることを再び強調しておく必要がある。
【0103】
上記結果として、画像の各目標ピクセル又はボクセルは、上記パラメータの一つ以上によって表現可能である。これらのパラメータは、暗示的に目標ピクセル又はボクセルと選択された窓のピクセル又はボクセルとの関係を示し、それぞれ選択された窓に対応する画像の一定の特徴を現すことができる。従って、これらパラメータの二つ以上、又はすべてを、ピクセル又はボクセルが画像内の有効に追跡可能なランドマークと一致するかどうかを求めるための数値としてみなすことで、非常に信頼性のあるランドマークの識別を提供することになる。
【0104】
目標ピクセル又はボクセルと関連付けられた少なくとも上述のパラメータの一部を利用して有効に追跡可能なランドマークを選択するための第一の方法によれば、例えば閾値の数値の定義を実験データに基づかせることで、閾値を求めることができる。
【0105】
いったん閾値が上述の各目標ピクセル又はボクセルに対するベクトルの成分である各パラメータに対して定義されると、前記パラメータのそれぞれ又はパラメータのグループの組み合わせと、対応する閾値とを比較することによって、目標ピクセル又はボクセルが有効に追跡可能なランドマークと対応するかどうかを求めるための基準を設定する関数を定義することができる。
【0106】
非常に多数の関数が異なる基準に対して定義可能であり、異なる可能性が考慮される。
【0107】
最も単純だが非常に時間を要する解は、ピクセル又はボクセルをコード化するベクトルの成分である各パラメータの値に対して閾値を定義することであり、各パラメータを関連する閾値と比較して、パラメータが閾値より高いか低いかを調べることである。各パラメータに対する比較結果は、パラメータが関連する閾値より低いか高いか、又は比較結果が有効に追跡可能なランドマークかそうでないかを、例えば0か1かの数値によって表現可能である。有効に追跡可能なランドマークに対して、特徴パラメータの幾つかは負の比較結果、即ち、前記特徴パラメータがランドマークが有効に追跡可能でないことを示す比較結果を有することもあるので、特徴パラメータと関連する閾値との各個別の比較結果の数値からなる成分を有するベクトルは、各特徴パラメータとその閾値との比較結果を要約する数値からなる包括的評価パラメータを得るために、さらに処理可能である。
【0108】
このことは、特徴パラメータの一つ又はサブセットに対する比較結果がランドマークが有効に追跡可能ではないと示しているにもかかわらず、有効に追跡可能なランドマークを識別することを可能にする。
【0109】
上記特徴パラメータと関連する閾値との比較結果の包括的評価は、各比較結果の関数として、又は特徴パラメータの関数として包括的評価パラメータを算出することで求めることができる。
【0110】
例えば包括的評価パラメータは、特徴パラメータのセットの平均値として算出することができる。この包括的評価パラメータは、包括的閾値と比較することができる。この場合、包括的評価パラメータに対する包括的閾値は、各特徴パラメータと関連する各単一閾値の関数として、例えば閾値の平均値として、求めることができる。
【0111】
包括的評価パラメータを求める他の方法は、特徴パラメータにより形成された成分からなるベクトルを定義し、次に前記ベクトルの係数、即ちその長さを計算することである。包括的評価パラメータである前記ベクトルの係数と比較される包括的閾値は、各単一特徴パラメータと関連する閾値によって形成された成分からなる閾値ベクトルの係数として算出することができる。
【0112】
また別な方法は、対応する特徴パラメータの対応する閾値との比較結果が有効に追跡可能なランドマークを示すか否かに依存する、例えば0又は1といった2値を前提とする、比較評価変数によって、各特徴パラメータの対応する閾値との比較結果をパラメータ化することである。
【0113】
このように、特徴パラメータと対応する閾値との比較結果が数値的に表現可能であり、そして、包括的結果は、その成分が比較評価パラメータ値であるベクトル形式で表現可能である。包括評価パラメータが二つの離散値0と1のうち一つを有することで定義される上記例に関し、比較評価ベクトルの成分は、0と1の数値のセットから成る。
【0114】
上記オプションと同様に、包括的評価パラメータは、この比較評価ベクトルの係数として、又はその成分の平均値として計算可能である。
【0115】
上記方法の別の特徴として、特徴パラメータの一つ以上のサブセットが定義可能であり、それは前記サブセットの特徴パラメータの関数としての二次特徴パラメータを計算するために用いられる。
【0116】
あらゆる種類の関数が、二次特徴パラメータを求め、又は算出するために利用可能である。
【0117】
二次特徴パラメータを求めるための特別な方法によれば、特徴パラメータの各サブセットは、目標ピクセル又はボクセル又は前記ピクセル又はボクセル窓のピクセル又はボクセルのそれぞれ又は一部を表す数値的外観の同じ種類の関数に関する特徴パラメータと、
窓のピクセル又はボクセルを表現する数値行列又は同一の数値行列の同一の変換についての一つ以上の特徴パラメータの同じ種類の関数に関する特徴パラメータと、を備える。
【0118】
これによって、特徴パラメータのセット全体が、各サブセットの特徴パラメータの関数として求められた二次特徴パラメータによって要約される特徴パラメータのサブセットに分割される。
【0119】
実際の例によれば、目標ピクセル又はボクセルを備えた一定の窓内のピクセル又はボクセルの信号強度と、前記目標ピクセル又はボクセルの信号強度の値は、特徴パラメータのサブセットを形成することとして定義される。これらの値を要約する二次特徴パラメータは、目標ピクセル又はボクセル又は窓に含まれるピクセル又はボクセルの信号強度の平均信号強度として計算可能である。
【0120】
同一の平均値関数が、窓のピクセル又はボクセルと目標ピクセル又はボクセルを表す数値画像データの行列の特異値から、及び/又は勾配行列の固有値から、又はヘッシアン行列の固有値から、又は自己相関行列の特異値から、二次特徴パラメータを計算するために利用可能である。
【0121】
上記ステップは、一定数のパラメータを有する一次特徴パラメータセットを、より少ない数のパラメータを有する二次特徴パラメータセットに減じることに存する。
【0122】
上述のように一次特徴パラメータとして定義可能な特徴パラメータのセットの評価と同様に、二次パラメータは、ピクセル又はボクセルが有効に追跡可能なランドマークに対応するか、有効に追跡可能でないランドマークに対応するかを決定するために用いることができる。前記二次パラメータは、一次特徴パラメータに適用したものと同じステップによって処理可能である。従って、二次パラメータをその二次閾値と比較するステップが実行可能であるのと同様に、二次閾値は定義可能である。またこの場合に、包括的閾値と比較される包括的比較評価パラメータも求めることができる。
【0123】
また上述の基準に従って特徴パラメータをグルーピングすることによって得られた前記二次特徴パラメータの場合において、単一二次特徴パラメータの組み合わせと対応する閾値の組み合わせは、例えば加重和又はその他の加重した組み合わせによって作成され、そして、そのように定義された包括的二次特徴パラメータと包括的二次特徴閾値のレベルとで比較が実行可能である。
【0124】
原理的に、一次特徴パラメータとその対応する閾値との比較を実行する上述のそれぞれの方法は、二次特徴パラメータのセットに対しても適用できる。特に、例えば二次特徴パラメータが閾値より大きいか小さいかによって0と1を取るような、各特徴パラメータとその関連する閾値との間の比較に関する異なる離散値を含む比較結果ベクトルに対して提供することにメリットがある。
【0125】
次に、包括的結果を和として、特に比較結果ベクトルの成分の加重和として求めることができる。
【0126】
パラメータ数が減ることによって、計算と評価がより簡単になることが長所である。さらに、二次特徴パラメータを求めるために選択された関数の種類が、正しくない又は統計的許容誤差の外にある一つ以上の前記パラメータの値の範囲を制限することを可能にする。
【0127】
比較のために包括的評価パラメータを求める場合において、ピクセル又はボクセルが有効に追跡可能なランドマークかどうかを求める場合の妥当性に相応じて、各一次又は二次特徴パラメータを加重する加重値を定義することが可能である。図20は、上述のごとく定義された前記ピクセル又はボクセルの特徴を考慮することによって、ピクセル又はボクセルが有効に追跡可能なランドマークと対応するか否かを評価するための上述のステップを表現した概要図である。
【0128】
図20で示されたケースは、一次特徴としてここで定義された目標ピクセル又はボクセルの特徴パラメータの選択されたサブセットの関数としての、各目標ピクセル又はボクセルに対して二次特徴パラメータを求めた場合を指す。
【0129】
以下で説明することではあるが、前記処理ステップのメカニズムも、第一特徴パラメータに適用することができる。
【0130】
図4ないし8に関連して、先に説明したように、第一特徴パラメータのセットが求められる。目標ピクセルは、I(m, n)で表される。隣接するピクセルは、3×3ピクセル寸法を有する窓に含まれるものとして選択される。一次特性は、選択された窓内の隣接部の強度Iと、目標ピクセルの強度と、目標ピクセルの一つを含む選択された窓に含まれるピクセルの信号強度の行列の特異値siと、目標ピクセルの一つを含む選択された窓に含まれるピクセルの信号強度の行列の勾配行列の固有値と、目標ピクセルの一つを含む選択された窓に含まれるピクセルの信号強度の行列のヘッシアン行列の固有値と、で構成される。
【0131】
上述のように求められたあらゆる別の特徴パラメータは、図示された特徴パラメータに付加したり、置き換えたりすることができるので、図20の例で選択された特徴パラメータの種類に、本発明は限定されるものではないことを強調する必要がある。
【0132】
一次特徴パラメータの前記種類のそれぞれは、特徴パラメータのサブセットとして定義される。これは、中括弧によって図20に示される。二次特徴パラメータは、関数1として示された関数によって、一次特徴パラメータの各サブセットから求められる。対応する一次特徴パラメータのサブセットから二次特徴パラメータを求めるための関数は、各サブセット間で同一又は異なるものであってよい。
【0133】
関数の例は、二乗関数又は係数関数等の和の平均値関数又は加重関数に存することができる。
【0134】
他の種類の関数は、サブセットの一次特徴パラメータの少なくとも幾つか、又は一つだけ除いた全部を除去するために、閾値との比較を提供する。各サブセットに対し、二次特徴パラメータはそのように求めることができ、図20において、この二次特徴パラメータは、P1,P2,P3,P4で示される。ここで、目標ピクセルが有効に追跡可能なランドマークと認められるか否かを求めるために、上述のすべてのステップが実行可能である。各二次特徴パラメータに対して個々の閾値を求めることができ、そして各二次特徴パラメータと対応する閾値との個々の比較が実行可能である。ピクセルが有効に追跡可能なランドマークか否かの決定は、上述の包括的比較評価プロセスによるこれらの比較結果に基づいて行われる。
【0135】
本実施例において、包括的二次特徴パラメータが図20の関数2によって示される二次特徴パラメータの関数として求められる。関数2は、関数1又は関数1に対して示された例のうち一つと等しくすることができる。
【0136】
上述のすでに開示した方法のうち一つに従って求められた包括的閾値は、例えば各二次特徴パラメータに対する個々の閾値を求め、次に、関数2と同じ又は異なる関数によって、包括的二次特徴パラメータをもとめるために、この個々の閾値のセットから包括的閾値を求めることによって求められる。
【0137】
包括的二次特徴パラメータと包括的閾値との比較は、比較関数によって実行され、そして、比較結果はピクセルI(m, n)が有効に追跡可能なランドマークと一致するか否かの評価を与える。
【0138】
代替案として、全ての一次特徴パラメータから包括的評価パラメータを求めるために関数を選ぶことができ、そして、閾値から求めるための別の関数又は同じ関数が、各一次特徴パラメータに対して包括的閾値を定義することを認めなければならない。比較関数は、次に、この一次特徴パラメータに対する包括的閾値に適用され、また包括的評価特徴パラメータに適用される。
【0139】
好ましい実施例において、関数1と関数2として示された異なる関数が、線形結合であって、好ましくは前記一次又は二次特徴パラメータの加重線形結合であって、例えば加重和である。同様に、第二特徴パラメータに対する閾値と包括的閾値とは、各一次パラメータに対する閾値と、各二次パラメータに対する閾値それぞれの加重和として計算可能である。
【0140】
特別に長所を有する実施例は、以下のステップからなる。各一次特徴パラメータとそれらの閾値とを比較するステップと、各成分が各一次パラメータと対応する閾値との比較結果からなる比較結果ベクトルを定義するステップ。もし一次パラメータが閾値より小さいとき、値0とし、もし一次パラメータが対応する閾値より大きいとき、値を1とするように比較結果を定義するステップ。比較結果ベクトルを値0と1の連続で形成するステップ。比較結果ベクトルの成分を組合わせるステップ。特に比較結果ベクトルの成分の加重和を求めるステップと、そして、この結果を包括的閾値と比較するステップ。包括的閾値は、実験又は各一次特徴パラメータに対する閾値の正規化された加重和によって求めることができる。比較結果ベクトルの成分の加重和と包括的閾値との比較結果に関連して、ピクセル又はボクセル又はピクセル又はボクセルのグループが有効に追跡可能なランドマークであるかどうかを決定するステップ。
【0141】
代替的に、上述のように二次パラメータを形成するために一次パラメータがグループ化された場合に、同じステップが二次特徴パラメータに対して適用可能である。この場合、比較結果ベクトルは、その成分が各二次特徴パラメータと対応する閾値との比較結果として定義される。また、ここで数値、特に0と1が、二次パラメータが対応する閾値より高いか低いかという条件にそれぞれ対応する比較結果を表現するために選択可能である。次に、ピクセル又はボクセル又はピクセル又はボクセルのグループが追跡すべき良いランドマークであるか否かということに関連するピクセル又はボクセル又はピクセル又はボクセルのグループの有効性が、比較結果ベクトルの成分の加重和と包括的閾値との比較によって求められる。
【0142】
特徴パラメータは、前記特徴パラメータの種類又はその数学的意味に依存した異なる方法によっても評価可能である。
【0143】
従って、例えば一次特徴パラメータのリストに含まれるヘッシアン行列の固有値を考慮すると、前記ヘッシアン行列の固有値の値から、ピクセル又はボクセルの重要性を評価することができる。
【0144】
ウェーブレット変換係数を考慮することによっても、画像のピクセル又はボクセルの重要性を評価することができる。この技術は、文献「Wavelet-based Salient Points: Applications to Using Color and Texture Features, E. Loupias, N. Sebe, Visual' 00, pp. 223-232, Lyon, France」と「Image Retrieval using Wavelet-based Salient Points, Q. Tian, N. Sebe, M.S. Lew, E. Loupias, T.S. Huang, Journal of Electronic Imaging, Vol. 10, No.4, pp. 835-849, October, 2001」とにより詳細に開示されている。
【0145】
従って、上述の開示内容に加えて、特徴がピクセル又はボクセルのものかを評価するための閾値が、追跡すべき特徴又はランドマーク即ち画像のピクセル又はボクセルの一つに対応し、幾つかの特徴パラメータは、ヘッシアン行列の固有値とウェーブレット変換係数に関する代替的な方法によって評価することができる。
【0146】
図21に、有効に追跡可能なランドマーク又は有効に追跡可能なランドマークでないものの質を評価する別の実施例を示す。ここで、いわゆる分類アルゴリズム、例えばクラスタ化アルゴリズム又は予測アルゴリズム、例えば人工の神経回路網又は統計アルゴリズムが用いられる。
【0147】
この実施例を実行可能にするために、第一に、有効に追跡可能なランドマークか否かが知られたランドマークに対応したピクセルのデータベースを構築する必要がある。ここで、i=1, ..., n の各ランドマークLiに対し、成分が前記ランドマークLiが特定された画像のピクセル又はボクセルの前記方法ステップによって求められた一次又は二次特徴パラメータであるベクトルが構築される。
【0148】
ベクトルの成分は、更にそれぞれ有効に追跡可能なランドマークか否かの質を示す二つの異なる値を推測可能な有効な変数を備える。
【0149】
前記データベースは、図21に概略的に示される。ここで、Liはランドマーク、Viは対応する特性のベクトルであり、「有効性」は、ランドマークLi又はそれに一致するピクセルが有効に追跡可能か否かを表現するパラメータである。データベースのレコードは、通常の手段で分類アルゴリズムを鍛えるために用いられる。ひとたびアルゴリズムが鍛えられると、ピクセルI(m, n)の特徴パラメータは、それによって処理され、そして分類アルゴリズムは出力として有効な変数の値と、ピクセルI(m, n)が有効に追跡可能なランドマーク又は有効でない追跡可能なランドマークと一致するか否かの予測又は決定を与える。
【0150】
本実施例及び図20の一つ前の実施例に従う評価プロセスの組み合わせも適用可能であることを強調しておくことは重要である。この場合、一次特徴パラメータを用いる代わりに、目標ピクセル又はボクセルと、鍛えられたデータベースのランドマークのための入力データベクトルを構築するために二次特徴パラメータのみが用いられる。
【0151】
上述の異なる評価ステップに対して、特に有効に追跡可能であると識別された二つのランドマークが互いに一定の距離にある場合に、別なランドマークを確認する又は廃棄するステップが実行可能であることも考慮されるべきである。これは、前記二つのランドマークの一次特徴パラメータの少なくとも一つ、又は一部又は全部、又は二次特徴パラメータの少なくとも一つ、又は一部又は全部を比較することによって、そして、選択関数を二つのランドマーク間の選択基準を表す比較結果に適用することによって、実行可能である。
【0152】
以下において、自動ランドマーク識別の例を図11ないし14を用いて示す。これらの図式的な例は、理解を容易にするために、二次元画像を用いている。三次元画像への方法の一般化は、ここでは簡潔に記述された数学的形式主義によって、概要が与えられる。
【0153】
三次元と、三次元ボリュームデータを表現する断面画像セットにおける自動特徴選択作業は、Lucas & Kanadeの二次元特徴検出アルゴリズムの実質的拡張である。
【0154】
時間の僅かに異なる間に記録された同一患者の二枚の画像ボリュームを、IとJと定義する。
【0155】
第一の画像ボリュームIにおいて、二つのボリューム間で追跡すべき特徴とみなされた一定数の選択されたボクセルが識別される必要がある。
【0156】
Iをオリジナルの第一画像ボリュームとし、Jを次のボリュームとする。ここで、対応する特徴は見出されていなければならない。
【0157】
第一ステップは、第一画像ボリュームにおける最小固有値の計算で構成される。これは、以下のように実行される。
【0158】
第一画像ボリュームは、三つの勾配ボリュームIx, Iy, Iz(画像ボリュームを表現するデカルト座標系のx-, y-, z方向)に変換され、x、y、zは各シングルボクセルの座標である。
【0159】
勾配ボリュームは、一般的にx方向に(2?x + 1)、y方向に(2?y + 1)、z方向に(2?z + 1)(ここで?x, ?y, ?z = 1, 2, 3, 4, ..., nボクセル)のサイズのいわゆる窓に含まれる隣接するボクセルに対する各ボクセルの強度勾配から算出される。
【0160】
目標ボクセルを中心とする隣接ボクセルの3×3×3アレイのサイズを考えると、3つの勾配ボリュームは、以下のように定義される。
a)x方向の勾配ボリューム
【数1】
b)y方向の勾配ボリューム
【数2】
c)z方向の勾配ボリューム:
【数3】
次のステップとして、各ボクセルに対し、三つすべての先に算出した勾配ボリュームを用いて、勾配行列Gが設定される。勾配行列は、以下のように生成される。
【数4】
各ボクセルの各勾配行列に対して、以下のように勾配行列から最小固有値?minが算出される。
定義:Graphics Gems IV, Paul S. Heckbert, 1994, Academic Press, S 193-98
【数5】
ここで
【数6】
なぜなら、Gは中心モーメントの行列だからである。
勾配行列Gの固有値は、次式で算出される。
【数7】
次に、固有値?mの最小値に対する閾値が定義される。
【0161】
次に、ボクセルが追跡可能な特徴を表現しているか否かを検討するために、以下の基準が適用される。
1.もし?mが閾値より大きい場合、ボクセルの実際の位置が、追跡すべき特徴のリストに 保存される。
2.もし?mが固有値の全ての最小値?mの最大値のパーセンテージより小さい場合、追跡 すべき特徴のリストから落とされる。
3.もし実際の特徴周りの定義された領域(調整可能、例えば3×3×3)において、他 の特徴が存在する場合、固有値のより小さな最小値?mを有する特徴が、追跡すべき特 徴のリストから落とされる。
4.もし特徴位置周りの3次元隣接ボクセルの平均信号強度(例えば? = 1)が、ゼロに近 い非常に低い値(例えば10)より小さいか等しい場合、その特徴は、周囲の空気の 中で乳房の外側にあるとみなされ、追跡すべき特徴のリストから廃棄される。
5.最後に、特徴の定義可能な最大数のみがキープされる。もしそれ以上の特徴が追跡す べき特徴のリストに存在する場合、より小さな最小値?mを有する特徴が、そのリスト から落とされる。
【0162】
これらのステップは、二次元の場合について図11ないし14で視覚的に説明されている。この場合、ピクセルのアレイが小さな四角のアレイとして示されている。第一の基準を満足する固有値の最小値を有する選択された特徴の位置に対応した異なるピクセルが示されている。これらの特徴と位置は、図3において対応するピクセルに異なる外観を与えることによって強調表示されている。図4において、5×5隣接部が、その5×5隣接部の内側を刻む円によって強調表示されている。
【0163】
図11ないし18の例において、グリッドはボクセルアレイを表現している。グリッドの各要素は、画像を形成するボクセルである。
【0164】
図11の例において、4種類のボクセルがランドマークとして区別され、対応する勾配行列の最小固有値に従って分類されている。102で示されたボクセルは、最大の固有値を有し、202で示されたボクセルは、大きな固有値を有し、302で示されたボクセルは、小さな固有値を有し、402で示されたボクセルは、最小の固有値を有する。
【0165】
図12において、各ボクセルに対して、直径5ピクセル/ボクセルのユークリッド距離内の隣接部が特徴として選択された中央ピクセルの周りに定義される。前記隣接部選択は、図4に円で示されている。
【0166】
図12において、このイベントの第一のケースが、502’と記されたより太い円で強調表示されている。ここで、特徴として選択された二つのピクセルが、ピクセルのユークリッド距離として上で定義された同一の隣接部内に存在する。これらの二つのピクセルは202’と302’で示され、上記定義に従って、特徴202’として選択されたピクセルは、第二のピクセル302’より大きな固有値を備えた勾配行列を有する。よって、この第二のピクセル302’は廃棄される。
【0167】
図13には次のステップを示し、図11で特徴として選択されたピクセル302’はもはや存在しない。図13において、隣接部を示す円を備えたピクセル202’’が設けられている。図12と同様の状況がここでも存在し、ピクセル202’’の隣接部に、ピクセル202’’より小さな固有値を有する勾配行列を持ったピクセル302’’が存在する。
【0168】
従って、この場合においても、特徴として選択されたピクセル302’’は追跡すべき特徴のリストから廃棄される。
【0169】
最後に、図14のピクセルアレイにおいて、元の特徴として選択されたピクセルの幾つかのみが存在している。ここで、単一の隣接部内において、一つの特徴以上は見られなくなっている。従って、追跡可能な特徴は、画像/ボリューム中に散開している。
【0170】
有効に追跡可能なピクセルを定義するための数学的形式化と基準と同様に、図11ないし14に示された本例は、目標ピクセルと選択された窓のピクセルとの強度信号に存し、また前記窓の画像データの勾配行列の固有値に存する特徴パラメータのセットに限られる。
【0171】
本発明の次のステップによれば、ピクセルの強度はさらに調査され、そして、所定値より低い強度を有する特徴として選択されたピクセルも、追跡すべき特徴に対応するピクセルのリストから廃棄される。ほとんどゼロに近い非常に低い閾値が、周囲空気内で見つけられた特性に用いられる。このステップは、図示しないが、当業者には自明である。
【0172】
ひとたび2次元又は3次元画像の特徴、即ち有効に追跡可能なランドマークとみなされたピクセル又はボクセルが求められると、前記特徴の追跡が実行される。本発明の好ましい実施例によれば、これはいわゆる「Lucas Kanade特徴追跡のピラミッド形の実施」を用いて実行される。詳細な記載は、Jean Yves Bouguetが記した記事「Pyramidal Implementation of the Lucas Kanade Feature Tracker, Description of the Algorithm」Intel Corporation Microprocessor Research Labs.発行、により与えられる。
【0173】
上述の出版物によって二次元の例について開示された理論を、以下で3次元の場合に更に適用する。上記方法とアルゴリズムとによって個別化された追跡すべき特徴を表わす各ボクセルに対し、以下のステップを実行する必要がある。
【0174】
ボリュームIにおけるポイントとしてuを、ボリュームJにおける対応するポイントとしてvを定義する。第一ステップにおいて、レベルL 0, ..., mにおけるボリュームを表現するILとJLとを備え、元のボリュームを表わす基底部(L = 0)を備えたボリュームIとJに対応する二つのピラミッドが設けられる必要がある。各以下のフロアのボリュームは、一つの平均値に対する直接の隣接部における一定量のピクセルを組合わせることによって寸法が削減される。追加のフロアの数は、ボリュームのサイズに依存する。
【0175】
4.2.以下のステップにおいて、いわゆる包括的ピラミッド形動きベクトルgが最上階で初期化される。
【数8】
全てのレベルLに対して、結果として生じる動きベクトルdLが次に計算される。
これは、以下のステップに従って実行される。
a)各ボリュームILに対して、ポイントuLの位置が次式によって計算される。
【数9】
ここでpは実際のボリューム位置である。
b)ボリュームのピラミッド形の表現の各レベルLにおいて、デカルト座標x、y、zの各方向において勾配が計算される。デカルト座標xに対し、ボリューム勾配が以下の式によって計算される。この等式は、他の二つのデカルト座標y、zに対しても適用される。
【数10】
c)次に、ボリューム勾配が以下の式による勾配行列を計算するために用いられる。
【数11】
ここで値?は、エリアサイズ又はボクセルに影響を与える隣接部を定義する。行列の要素は、上述の定義によるものと同様にして得られる。即ち、
【数12】
d)各レベルLに対して、反復ベクトル?が初期化される。
【数13】
e)次に、k=1から最大カウントKまで、又は最小置換
【数14】
まで、以下のボリューム差を計算する。
【数15】
以下の等式に従って、不整合なベクトルを計算する。
【数16】
次に、次式によってオプティカルフローベクトルが定義される。
【数17】
上で定義したオプティカルフローベクトルを用いて、反復動きベクトルが計算可能である。
【数18】
そしてこれをレベルLに対する最後のオプティカルフローと等しく設定する。
【数19】
上述のステップを各レベルに対し、最後のレベルL=0になるまで繰り返す。それによって、次のレベルL−1に対して、包括的オプティカルフローが上記計算値から算出される。
【数20】
f)最終的なオプティカルフローベクトルは、次式となる。
【数21】
g)次に、ボリュームIのポイントuに対応するボリュームJのポイントvの座標は、以下のように計算できる。
v = u + d
【0176】
図15ないし18は、明らかに実際のものの簡単表現である上述の追跡方法の効果の図による概略的な発想を与えようとするものである。
【0177】
図15ないし18において、異なる時間T0とT1に撮影された二枚の画像IとJが、図11ないし14の例のように、ピクセルの二次元アレイとして表現されている。
【0178】
追跡すべき特徴に対応し、上述の方法によって選択されたピクセル20が、画像Iの二次元ピクセルアレイにおいて示されている。画像Jにおいて同一の位置のピクセル20’’が、画像化された対象この場合乳房の動きが無い場合におけるピクセル20に対応している。ピクセル20’’の一定の隣接ピクセルが、ピクセル20’’を中心とした7×7ピクセルアレイに対応した灰色の影つきの四角で示されており、隣接部は番号10が付されている。
【0179】
画像Iの特徴20の新しい位置が20’で示されている。これは、ピクセル20の位置における乳房組織が、ピクセル20’の位置に時間とともに移動したことを意味している。明らかに、前記ピクセル20’は、追跡が実行される方法を説明するためのみに示されている。実際の場合において、このピクセルは未知であり、さもなければ追跡は不必要である。
【0180】
既に上で述べた数学的形式による選択された特徴に対応するピクセルを追跡する一般的な方法について、以下に図で与えれ、簡単のために二次元の場合に限定された実施例によって記載する。しかし、追跡の一般的数学的記述を踏まえた当業者には、それでも3次元の場合について追跡する方法を理解し、評価することができる。
【0181】
追跡の初めに、画像I内のピクセル20が、同じ位置、即ち画像Jのピクセル20’’にあるものと仮定する。このピクセル20’’を、以下では「追跡されたピクセル」と呼び、一方で特徴として選択された画像Iのピクセル20を「特徴ピクセル」と呼ぶこととする。まず追跡されたピクセル20’’の隣接部に対応した一定の領域の反転勾配行列にアルゴリズムを適用し、本例においては画像Jのピクセル20’’を中心とした前記7×7ピクセルアレイ(この値は調整可能)が計算される。この隣接部が、いわゆる追跡領域に相当する。Lucas & Kanadeの特徴追跡アルゴリズムのピラミッド形の実施を適用することで、これは3次元画像ボリュームに対しても適応されて、詳細は後述するが、画像J上の追跡された特徴の新しい位置を見つけることができる。
【0182】
画像J内の追跡された特徴20’の新しい位置は、オプティカルフローベクトルを決定する。これは、第一の反復ステップにおいて計算されて定義される。図15と16の比較から明らかなように、追跡領域は、追跡領域の中心が同じ組織領域20’の新しい位置により近くなるようにシフトされる。
【0183】
追跡は、画像I内の特徴を表わすピクセル20と画像J内の追跡領域の中心にある新しい識別された位置との間の不整合ベクトルの計算を繰り返すことによって、さらに反復的して適用される。これは、新しいオプティカルフローベクトルを定義することになり、図17に示すように、ピクセル隣接部又は追跡領域の新たなシフトを生じさせる。
【0184】
画像I内の特徴ピクセル20と追跡領域の中心間の不整合ベクトルを求めるステップは、不整合ベクトルが調整可能な最小値に達するまで続けられる。
【0185】
この状況を図18に示す。図15ないし17の図的表現に関連して、ピクセル20’は追跡領域がシフトするという事実を強調するために位置を変えておらず、追跡領域は最終的に画像Iのピクセル20としての乳房の同じ位置を表現する正しい新しい位置20’を中心とする位置に動いている。
【0186】
上記例は、ピクセルによって表現されたただ一つの特徴に対して実行されている。上記ステップは、画像内で選択された各特徴に対して実行する必要がある。
【0187】
すべての特徴がその新しい位置にまで追跡された後に、画像/ボリューム全体にわたって、ただし乳房組織を表現する領域内でのみ、異なる位置にあるオプティカルフローベクトルの一定量が存在する。他の登録アルゴリズムと比較した最大の長所は、画像/ボリューム全体にわたってオプティカルフローベクトルが散会することである。次に、逆オプティカルフローベクトルを画像/ボリュームJ内の各特徴と、その周囲の隣接部とに適用することによって、モーフィングが実行される。これは、新しい仮想画像/ボリュームJ'をもたらす。
【0188】
図19に、実際の診断画像に対する前記方法の効果を示す。左側の第一の画像セットは、造影剤投与無しでの時刻t0における乳房検査からの代表的なスライスを示す。左側の第二画像は、造影剤投与後の時刻t1における同じスライス位置を撮影したものである。左側の第三画像は、時刻t1のピクセルの信号強度から時刻t0のものをピクセル毎に引き算した結果を視覚化したものであり、いわゆる引き算画像である。時刻t1とt0間で乳房の相対的移動が無い場合、造影剤強調が生じない中で、結果の画像では黒色を示すゼロ信号強度のピクセルが多数を占める。時刻t1とt0間での移動は、画像/ボリュームIとJの同じ位置のピクセル単位でアーチファクトが発生するという効果を有し、これは組織の同じ領域には対応していない。本例においては、引き算画像は小さな矢印で示されたコントラスト強調された腫瘍を含む乳房を表わしている。乳房のほとんどの種類の悪性および良性腫瘍は、増加した血管新生によって生じる追加的な造影剤摂取を示す。両乳房の中央の白く残った構造は、移動と組み合わされた正常な乳房柔組織の弱い造影剤摂取の結果である。
【0189】
図19の右側は、前記登録アルゴリズムが適用された対応する画像を示す。
【0190】
時刻t0における画像はアルゴリズムによって修正されておらず、従って同じものを留めている。しかし、三次元特徴追跡後は、新たな仮想画像ボリュームが算出される。時刻t1’での右側の第二画像は、登録適用後の同じスライス位置の対応する仮想画像を示す。右側の最後の画像は、時刻t1’から時刻t0でのピクセルの信号強度を引く、ピクセル毎の引き算結果を視覚化したものであり、新たな(仮想)引き算画像である。一つの乳房は、小さな矢印で示されたコントラスト強調された腫瘍を表わしている。正常な乳房柔組織を表わす乳房内の領域は、幾分抑制された造影剤摂取を示し、動きアーチファクトの減少により、僅かに少ない輝きを見せている。両乳房上部の白い構造の大部分は、動きアーチファクトの補償によって消失している。
【0191】
一枚は造影剤投与前、一枚は投与後の二枚の連続した画像のみが検討されているので、図19の例は単純化されている。
【0192】
実際には、一連の断面画像(例えば50枚以上)が、造影剤の投与前後の一定時間毎に撮影され、引き算画像が連続した画像の各二枚間で算出される。
【0193】
1組の中に一枚以上の断面画像を確保するのに加えて、2組以上(例えば5組以上)の画像を次々に確保することもできる。これは、複数の引き算画像/ボリュームをもたらす。もちろん、開示したアルゴリズムは各組のすべての組み合わせに適用できるので、造影剤投与の前か後かに撮影されたにかかわらず、その間に生じた動きを等しく削除することができる。
【0194】
本発明に係る上述の方法は、コンピュータで実行されるプログラムとしてコード化されたアルゴリズムによって実行可能である。
【0195】
このプログラムは、一枚以上のフレキシブルディスクや、一枚以上の光磁気ディスク、一枚以上のCD−ROM又はDVD−ROMその他携帯用記憶装置や、安定したハードディスク等に保存できる。代替的に、ペンドライブ等の固体電子装置にも保存できる。
【図面の簡単な説明】
【0196】
【図1】胸部MRIで使用される乳房固定装置の概要図である。
【図2】引き算プロセスと登録ステップがない場合に生じる動きアーチファクトを示す。上部左右には、同じMRスキャナ位置の乳房の二枚の断面画像を示す。第二画像は第一画像の約2.5分後の造影剤投与後に測定した図である。両画像は、断面画像の全体的な3次元セットの代表的な断面図である。
【図3】本発明に係る登録ステップを示す概要図である。
【図4】10×10画素アレイで形成されたデジタルの又はデジタル処理された画像の概要図である。
【図5】図4の画像に対応する画像データアレイの概要図であり、そこにおいて各画素の視覚的外観、例えばグレースケールデジタル画像の画素の強度は数値で表現される。
【図6】ピクセルの外観、特に選択された窓のピクセルの強度(図5に示す)を表す数値によって目標ピクセルをコード化するために用いられた最先端技術に従うベクトルの第一実施例を示す概要図である。
【図7】図6に従うピクセルコード化方法の概要図を示し、そこにおいてベクトルの成分はより一般化された方法で表現されている。
【図8】座標m、nの目標ピクセルを中心とする一般的な3×3窓と、一般的な勾配行列と、数値的な偏導関数の計算の一つの方法を示す図であり、例えばSobel演算子やFrei-Chen演算子(Digital Image Processing, W. K. Pratt, John Wiley & Sons, Inc. New York, NY, USA, 1991)のような、より複雑な演算子を用いても導関数を計算することができる。
【図9】一般的なヘッシアン行列と数値的二次偏導関数の一つの計算方法を示す。
【図10】本発明に係る目標ピクセルを表現するベクトルの概要図である。
【図11】有効に追跡可能なランドマークの自動測定法の一例の概要図である。
【図12】有効に追跡可能なランドマークの自動測定法の一例の概要図である。
【図13】有効に追跡可能なランドマークの自動測定法の一例の概要図である。
【図14】有効に追跡可能なランドマークの自動測定法の一例の概要図である。
【図15】追跡がどのように実行されるかを示す概要図である。
【図16】追跡がどのように実行されるかを示す概要図である。
【図17】追跡がどのように実行されるかを示す概要図である。
【図18】追跡がどのように実行されるかを示す概要図である。
【図19】造影剤投与前後の異なる時間に撮影された二枚の画像と、特徴追跡と画像登録を実行した後の、両方の画像の引き算で得られた画像とを、図の左側に示した図2の例と比較して示すことによる本発明の結果を表した図である。図2から図19までのすべての画像は、登録以外の他の処理は一切しないで同一のMR試験で撮影されたものである。
【図20】対応する特徴パラメータによって、ピクセルの有効に又は非有効に追跡可能な特徴の質を測定するためのステップからなる本発明にかかる方法の第一実施例を示す図である。
【図21】対応する特徴パラメータによって、ピクセルの有効に又は非有効に追跡可能な特徴の質を測定するためのステップからなる本発明にかかる方法の第二の実施例を示す図である。
【特許請求の範囲】
【請求項1】
a)ピクセル又はボクセルからなる2次元又は3次元アレイで構成された、少なくとも第一と第二デジタル画像又は同じ対象の一連の断面画像を提供するステップと、
b)第一画像又は画像セットの中に、ランドマークとして設定された一定数のピクセル又はボクセルを選択し、追跡すべき前記特徴のリストを生成することで、一定数のいわゆる特徴であるランドマークを定義するステップと、
c)第一画像から第二画像又は、第一から第二の画像位置間のオプティカルフローベクトルによって決定された後の時点で得られた画像又は画像セット又は、特徴として選択された各ピクセル又はボクセルの後の時点で得られた画像又は画像セットから、特徴として選択された各ピクセル又はボクセルの位置を追跡するステップと、
d)第一画像と第二画像或いは、逆オプティカルフローベクトルを前記第二画像又は前記画像セットのピクセル又はボクセルの位置に適用することによって後の時点で得られた画像又は画像セット、を登録するステップと、
からなる、対象の移動によって生じる画像アーチファクトを低減する、生物医学画像の登録方法であって、
B1)第一画像又は第一の一連の断面画像の各ピクセル又はボクセル周りに隣接した、限られた数で構成されるピクセル又はボクセルを定義するステップと、
B2)各目標ピクセル又はボクセルに対して、外観を記載する数値パラメータ、すなわち前記ピクセル又はボクセル窓の前記目標ピクセル又はボクセルとピクセル又はボクセルの部分とからなるいわゆる数値外観パラメータの関数として計算され、或いは、前記窓のピクセル又はボクセルを表現する数値パラメータのマトリックス又は数値パラメータの前記マトリックスの変形で構成される一つ以上の特徴パラメータの関数として計算される、一つ以上の特徴のあるパラメータを定義するステップと、
B3)前記目標ピクセル又はボクセルの特徴的なパラメータの関数としての有効に追跡可能なランドマークに存するピクセル又はボクセルを定義するステップと、
からなる自動追跡可能なランドマーク選択ステップが実行されることを特徴とする登録方法。
【請求項2】
もし前記目標ピクセル又はボクセルが追跡すべき有効なランドマークであった場合に、定義するための関数が、前記目標ピクセル又はボクセルの特徴パラメータの各一つ又は組み合わせの値と閾値とを比較することであることを特徴とする請求項1記載の方法。
【請求項3】
前記各特徴パラメータに対する閾値を求め、前記各特徴パラメータを対応する閾値と比較して、単一の比較結果の関数としての全体的評価パラメータから、有効に追跡可能なランドマークと有効でない追跡可能なランドマークの質を求めることを特徴とする請求項2に記載の方法。
【請求項4】
前記各特徴パラメータに対して閾値の関数としての全体的閾値が求められ、前記全体的評価パラメータと前記全体的閾値とを比較することによって、有効に追跡可能なランドマークと有効でない追跡可能なランドマークの質を求めることを特徴とする請求項3に記載の方法。
【請求項5】
有効に追跡可能なランドマークか否かに関して比較した各結果を数値で表現するために可変有効性が定義され、
部分又は各特性と対応する閾値との比較の前記可変有効性の値の関数として、全体的評価パラメータが求められ、そして
全体的な評価パラメータの関数として、有効に追跡すべきランドマークと有効でない追跡すべきランドマークの質が求められる
ことを特徴とする請求項3又は4に記載の方法。
【請求項6】
包括的評価パラメータから有効に追跡可能なランドマークと有効でない追跡可能なランドマークを求めるための関数が、前記包括的評価パラメータと包括的閾値との比較であることを特徴とする請求項5に記載の方法。
【請求項7】
特徴パラメータセットが一定数のサブセットに細分化されて、各サブセットに対し二次特徴パラメータがそのサブセットの特徴パラメータの関数として求められ、そして、
各二次特徴パラメータに対し閾値が定義されて、各二次特徴パラメータに対し有効に追跡可能なランドマークと有効でない追跡可能なランドマークの質が、前記二次特徴パラメータと対応する閾値との関数によって求められる
ことを特徴とする請求項1ないし6のいずれか1項に記載の方法。
【請求項8】
各二次特徴パラメータと対応する閾値が、特徴パラメータの対応するサブセットの閾値の関数として求められることを特徴とする請求項7に記載の方法。
【請求項9】
前記サブセットの特徴パラメータが、少なくとも共通の特徴を有し、
特に前記目標ピクセル又はボクセル又は前記ピクセル又はボクセル窓のピクセル又はボクセルのそれぞれ又は一部の外観を記載する数値的外観パラメータからなり、そして
前記窓のピクセル又はボクセルを表現する数値行列の一つ以上の特徴パラメータの同じ関数によって、又は前記数値行列の同じ変換によって求められる
ことを特徴とする請求項7又は8に記載の方法。
【請求項10】
前記包括的評価パラメータが二次特徴パラメータの関数として定義され、そして、
有効に追跡可能なランドマークと有効でない追跡可能なランドマークの質が、前記包括的評価パラメータの関数と、前記二次特徴パラメータに対する閾値の関数として存する包括的閾値と、によって求められることを特徴とする請求項7ないし10のいずれか1項に記載の方法。
【請求項11】
有効に追跡可能なランドマークと有効でない追跡可能なランドマークの質を求めるために、各特徴パラメータ及び/又は各二次特徴パラメータ及び/又は各包括的閾値が加重され、そして、その加重が特徴パラメータ又は二次特徴パラメータの関連性に従って求められることを特徴とする請求項2ないし10のいずれか1項に記載の方法。
【請求項12】
前記包括的評価パラメータ及び/又は前記二次特徴パラメータ及び/又は包括的閾値を求めるための関数が、一つ以上の特徴パラメータの係数の線形結合か、二乗値の合計か、前記特徴パラメータ又は前記第二パラメータ又は前記閾値が処理される非線形関数に存することを特徴とする請求項1ないし11のいずれか1項に記載の方法。
【請求項13】
前記包括的評価パラメータ及び/又は前記二次特徴パラメータ及び/又は前記包括的閾値を求めるための関数が、各特徴パラメータ又は第二パラメータ又は閾値を加重するためにも用いられることを特徴とする請求項12に記載の方法。
【請求項14】
前記目標ピクセル又はボクセルの有効に追跡可能なランドマークと有効でない追跡可能なランドマークの質が、分類アルゴリズムを用いて画像の各ピクセル又はボクセルの特徴パラメータを処理することによって求められることを特徴とする請求項1に記載の方法。
【請求項15】
前記各目標ピクセルのIDが、ベクトル処理されるためにコード化され、前記ベクトルの成分が、前記ピクセル又はボクセルの一つ以上の特徴パラメータに存することを特徴とする請求項14に記載の方法。
【請求項16】
特徴パラメータセットが一定数のサブセットに細分化され、
前記特徴パラメータの代わりに、二次特徴パラメータが前記分類アルゴリズムによって処理され、そして、前記二次特徴パラメータが、特徴パラメータのサブセットの関数として求められることを特徴とする請求項14又は15に記載の方法。
【請求項17】
各ピクセル又はボクセルに対して前記分類アルゴリズムによって処理された入力ベクトル成分が、各特徴パラメータと対応する閾値との比較結果、又は各二次特徴パラメータと対応する閾値との比較結果を数値的に表現する可変有効性に存することを特徴とする請求項14又は15に記載の方法。
【請求項18】
前記分類アルゴリズムが、クラスタ化アルゴリズム又は予測アルゴリズムであることを特徴とする請求項14ないし17のいずれか1項に記載の方法。
【請求項19】
前記分類アルゴリズムが、人工の神経ネットワークであることを特徴とする請求項19記載の方法。
【請求項20】
一定数の有効なランドマークが、有効に追跡可能なランドマークとして識別された既知の場合の一定数の画像を提供するステップと、
上記ステップB1)とB2)に存する自動追跡可能ランドマーク選択ステップを適用することによって、一定数の画像の中で有効に追跡可能であると識別された前記ランドマークに対応して前記ピクセル又はボクセルに対する特徴パラメータセットを求めるステップと、
有効に追跡可能であると識別された各ランドマークに一意的に関連付けられ、成分として前記有効に追跡可能なランドマークと一致する前記ピクセル又はボクセルの特徴パラメータを備えたベクトルを生成するステップと、
変数の所定の数値によって有効に追跡可能な質を表現して、前記数値を、有効に追跡可能なランドマークと一致する各ピクセル又はボクセルをコード化するベクトルに関連付けるステップと、
有効に追跡可能なランドマークと一致する各ピクセル又はボクセルをコード化する各ベクトルが、分類アルゴリズムのためのトレーニングデータベースの登録を形成するステップと、
前記データベースによって分類アルゴリズムをトレーニングするステップと、
鍛えられた分類アルゴリズムの入力に前記目標ピクセル又はボクセルの前記特徴パラメータを備えたベクトルを提供することによって、目標ピクセル又はボクセルの有効に追跡可能なランドマークの質を求めるステップと、
からなることを特徴とする請求項14ないし19のいずれか1項に記載の方法。
【請求項21】
前記トレーニングデータベースが、有効に追跡可能ではないと知られたランドマークと一致する画像のピクセル又はボクセルのベクトル形式に組織化され、
有効に追跡可能でないランドマークの質を、追跡可能なランドマークの質を表現する前記変数の値とは異なる変数の所定の数値によって表現する特徴パラメータが備えられ、
一方で有効に追跡可能でないランドマークの質を表現する前記数値が、既知の有効に追跡可能でないランドマークと一致した各ピクセル又はボクセルをコード化するベクトルに加えられたことを特徴とする請求項20に記載の方法。
【請求項22】
画像のピクセル又はボクセルが、いわゆるLucas & Kanade自動特徴追跡法又はアルゴリズムによって、事前に又は平行して処理されることを特徴とする請求項1ないし22のいずれか1項に記載の方法。
【請求項23】
前記Lucas & Kanadeアルゴリズムが、
a)ボクセルの3次元アレイによって形成された第一3次元画像ボリュームIを提供するステップと、
b)目標ボクセルを取り囲み、例えば前記目標ボクセルを中心とした3×3×3アレイのような調整可能なサイズを有する隣接ボクセルを各目標ボクセルに対して定義するステップと、
c)デカルト座標系を定義し、前記デカルト座標系の各軸に対し、隣接ボクセルに対して第一画像ボリュームの各ボクセルの強度勾配によって形成されるいわゆる勾配ボリュームを計算するステップと、
ここで傾斜ボリュームは次式で定義され、
ア)x方向の勾配ボリューム:
【数1】
イ)y方向の勾配ボリューム:
【数2】
ウ)z方向の勾配ボリューム:
【数3】
d)各目標ボクセルに対し、次式で定義されるいわゆる勾配行列を計算するステップと、
【数4】
e)前記3次元画像ボリュームの各目標ボクセルの各勾配行列に対して、以下の等式を適用することで最小固有値?mを計算するステップと、
【数5】
ここで
【数6】
に対して、前記勾配行列の最小固有値を計算する
【数7】
f)前記最小固有値?mの閾値を定義するステップと、
g)以下の基準を満足する、対応する勾配行列が最小固有値?mを有する画像Iの追跡すべき特徴としての各ボクセルを選択するステップと、
i)?mが前記閾値より大きく、
ii)?mが固有値の全ての最小値?mの最大値の割合より小さくなく、
iii)もし特徴として選択された別のボクセルが、選択された特徴として維持されている目標ボクセル回りの定義された隣接ボクセルに存在する場合、勾配行列がより大きな固有値の最小値?mを有するボクセルの一つのみが特徴として選択され、他の一つは追跡可能な特徴のリストから廃棄され、
iv)特徴として選択されたボクセルの回りの3次元隣接ボクセルの平均信号強度が、調整
可能な平均信号強度閾値より大きい、
以上の各ステップを備えることを特徴とする請求項22に記載の方法。
【請求項24】
特徴として選択されたボクセル回りの隣接ボクセルの平均信号強度に対する閾値が約10に設定されることを特徴とする請求項23に記載の方法。
【請求項25】
前記特徴パラメータが、選択窓のピクセル又はボクセルの強度値に存し、また目標ピクセル又はボクセルの強度値に存することを特徴とする請求項1ないし24のいずれか1項に記載の方法。
【請求項26】
前記特徴パラメータが、選択された窓のピクセル又はボクセルの画像データを備えた数値行列の特異値に代替的に又は付加的に存することを特徴とする請求項1ないし25のいずれか1項に記載の方法。
【請求項27】
前記特徴パラメータが、選択された窓のピクセル又はボクセルを表現する数値行列の勾配行列の固有値に代替的に又は付加的に存することを特徴とする請求項1ないし26のいずれか1項に記載の方法。
【請求項28】
前記特徴パラメータが、選択された窓のピクセル又はボクセルを表現する数値行列のヘッシアン行列の固有値に代替的に又は付加的に存することを特徴とする請求項1ないし27のいずれか1項に記載の方法。
【請求項29】
前記特徴パラメータが、一つ以上のウエーブレットに基づく関数を用いて前記数値行列を処理することで得られた前記窓のピクセル又はボクセルを表現する数値行列のウエーブレット変換係数の一つ以上又はその組み合わせに代替的に又は付加的に存することを特徴とする請求項1ないし28のいずれか1項に記載の方法。
【請求項30】
前記特徴パラメータが、いわゆる行列の同時生起変換に代替的に又は付加的に存することを特徴とする請求項1ないし29のいずれか1項に記載の方法。
【請求項31】
前記特徴パラメータが、前記窓のピクセル又はボクセルを表現する数値行列の自己相関変換の一つ以上の係数に代替的に又は付加的に存することを特徴とする請求項1ないし30のいずれか1項に記載の方法。
【請求項32】
前記特徴パラメータが、
前記窓のピクセル又はボクセルを表現する数値行列の固有値又は特異値及び/又は、
前記窓のピクセル又はボクセルを表現する数値行列の勾配行列又はヘッシアン行列の固有値及び/又は、
ウェーブレット変換係数の一つ以上の係数及び/又は、
自己相関変換の一つ以上の係数及び/又は、
前記窓のピクセル又はボクセルを表現する数値行列の同時生起行列の特異値
の組み合わせに存することを特徴とする請求項1ないし31のいずれか1項に記載の方法。
【請求項33】
同一対象の少なくとも二つ以上の画像が得られるか、又は異なる時間に得られた少なくとも二つ以上の画像が提供され、
第一画像ボリュームIにおいて有効に追跡可能なランドマークと一致して区別された各ボクセルに対して、前記特徴が、第二の遅い時間に得られた同一対象の追加の画像Jのボクセルアレイにおけるその位置に対して追跡され、
前記追跡が、いわゆる「Lucas&Kanade特徴追跡」のピラミッド形の実施によって行われることを特徴とする請求項1ないし32のいずれか1項に記載の方法。
【請求項34】
前記いわゆる「Lucas&Kanade特徴追跡」のピラミッド形の実施が、
画像Iにおける特徴として選択されたボクセルの座標をそれぞれ有する、画像Iのボクセルの3次元アレイに対応する画像ボリュームIにおけるポイントとしてのポイントuと、画像Jのボクセルの3次元アレイに対応する画像ボリュームJにおける対応するポイントvとを定義するステップと、
元のボリュームを表現する下部構造(L=0)と共に、レベルL 0,...,m におけるボリュームを表現するILとJLとを用いて画像ボリュームIとJの二つの理想的なピラミッドを構築するステップと、
最上階Lmにおいて初期化されるいわゆる包括的ピラミッド形の運動ベクトルgを定義するステップと、
【数8】
i)uLとしてポイントuの位置を定義し、以下の式で前記位置を算出し、ここでLは
ピラミッドの実際の階又はレベルの値を有し、
【数9】
ここでpは実際のボリューム位置であり、
ii)ボリューム勾配に対する以下の式に従って、デカルト座標x、y、zの各方向にお
けるボリューム勾配を計算し、
【数10】
iii)以下の等式に従って、勾配行列を計算するために勾配ボリュームを使用し、
【数11】
ここで
【数12】
ここで数値?は、エリアサイズ又は追跡される特徴を表現する目標ボクセルに影響を与える隣接ボクセルを定義し、
iv)各レベルLに対し、次式で定義される反復ベクトルを初期化し、
【数13】
v)k=1から最大カウントKまで、又は最小置換
【数14】
まで、以下のボリューム差を計算し、
【数15】
vi)以下の等式に従って、不整合なベクトルを計算し、
【数16】
vii)以下の式に従って、オプティカルフローベクトルを求め、
【数17】
ここで、G-1はステップiii)で求めた行列Gの逆勾配行列であり、
viii)上で定義されたオプティカルフローベクトルを用いて以下のように反復運動ベクト
ルを計算し、
【数18】
そしてこれをレベルLに対する最後のオプチカルフローと等しくし、
【数19】
ix)最後のレベルL=0が次式で定義される最終オプチカルフローベクトルに達するま
でステップi)からviii)を各レベルに対し繰り返し、
【数20】
x)以下の式でボリュームにおける追跡すべき特徴を表現するポイントuに対応するボ
リュームJにおけるポイントvの座標を求め、
v=u+d
xi)画像Iにおいて選択された特徴に対応する各ボクセルに対して上記方法を繰り返
す、
上記各ステップを備えることを特徴とする請求項33記載の方法。
【請求項35】
画像Iにおける選択された特徴に対応するボクセルのポイントUに対応する特徴を表現するボクセルに対応すると認定された画像Jにおける各ポイントvに反転オプティカルフローベクトルを適用するか、又は
画像Iにおいて選択された特徴であると認定されたボクセルにそれぞれ対応するポイントuにオプティカルフローベクトルを適用することによって、
二つの画像IとJを登録することを特徴とする請求項34に記載の方法。
【請求項36】
a)MRI又は超音波で得られた同一対象の少なくとも第一と第二ディジタル画像又はディジタル処理された画像又は断面画像セットを提供するステップと、
ここで前記画像はピクセル又はボクセルの二次元又は三次元アレイによって形成され、
ここで同一組織又は組織領域又は同一の身体構造上の領域の操作が、前記同一の組織又は組織領域又は前記身体構造上の領域に造影剤が存在する状態で、二度目又はそれ以降の時間で実施され、
b)ランドマークとして設定されたピクセル又はボクセルの一定数を選択し、追跡すべき前記特徴のリストを生成することによって、第一画像又は画像セット内に一定数のランドマーク、いわゆる特徴を定義するステップと、
c)特徴として選択された各ピクセル又はボクセルに対する第一から第二画像又は画像のセットからオプティカルフローベクトルを求めることによって、第一から第二画像又は画像のセットから特徴として選択された各ピクセル又はボクセルの位置を追跡するステップと、
d)第二画像又は画像セットのピクセル又はボクセルに反転オプティカルフローを適用することによって、第一と第二画像又は画像セットを登録するステップと、
からなる造影剤で強調した画像診断法、特に造影剤で強調したMRI又は超音波画像診断法との組合わせで提供されることを特徴とする請求項33ないし35のいずれか1項に記載の方法。
【請求項37】
前記特徴選択ステップの後で、前記特徴追跡ステップを実行する前に、
B1)第一画像又は第一の断面画像セットの各ピクセル又はボクセルの回りに隣接するピクセル又はボクセルを定義するステップと、
ここで前記隣接するピクセル又はボクセルは、制限された数のピクセル又はボクセルを備え、単一画像の場合に二次元隣接部又はボクセルが選択され、断面画像セットの場合に三次元隣接部又はボクセルが選択され、
B2)各ピクセル又はボクセルに対し、前記隣接部又はボクセルの全ピクセル又はボクセルの平均信号強度を求めるステップと、
B3)平均信号強度閾値を定義するステップと、
B4)ステップB2で各隣接部又はボクセルに対して求められた平均信号強度値と、前記平均信号強度値とを、前記平均信号強度閾値に対して比較するステップと、
B5)前記隣接部又はボクセルの前記平均信号強度値が前記ステップB4での閾値より大きい場合、前記ピクセル又はボクセルを追跡すべき特徴として定義することにより、追跡すべき特徴のリストに加えるステップと、
からなる追加の自動特徴選択ステップを実行することを特徴とする請求項33ないし36のいずれか1項に記載の方法。
【請求項1】
a)ピクセル又はボクセルからなる2次元又は3次元アレイで構成された、少なくとも第一と第二デジタル画像又は同じ対象の一連の断面画像を提供するステップと、
b)第一画像又は画像セットの中に、ランドマークとして設定された一定数のピクセル又はボクセルを選択し、追跡すべき前記特徴のリストを生成することで、一定数のいわゆる特徴であるランドマークを定義するステップと、
c)第一画像から第二画像又は、第一から第二の画像位置間のオプティカルフローベクトルによって決定された後の時点で得られた画像又は画像セット又は、特徴として選択された各ピクセル又はボクセルの後の時点で得られた画像又は画像セットから、特徴として選択された各ピクセル又はボクセルの位置を追跡するステップと、
d)第一画像と第二画像或いは、逆オプティカルフローベクトルを前記第二画像又は前記画像セットのピクセル又はボクセルの位置に適用することによって後の時点で得られた画像又は画像セット、を登録するステップと、
からなる、対象の移動によって生じる画像アーチファクトを低減する、生物医学画像の登録方法であって、
B1)第一画像又は第一の一連の断面画像の各ピクセル又はボクセル周りに隣接した、限られた数で構成されるピクセル又はボクセルを定義するステップと、
B2)各目標ピクセル又はボクセルに対して、外観を記載する数値パラメータ、すなわち前記ピクセル又はボクセル窓の前記目標ピクセル又はボクセルとピクセル又はボクセルの部分とからなるいわゆる数値外観パラメータの関数として計算され、或いは、前記窓のピクセル又はボクセルを表現する数値パラメータのマトリックス又は数値パラメータの前記マトリックスの変形で構成される一つ以上の特徴パラメータの関数として計算される、一つ以上の特徴のあるパラメータを定義するステップと、
B3)前記目標ピクセル又はボクセルの特徴的なパラメータの関数としての有効に追跡可能なランドマークに存するピクセル又はボクセルを定義するステップと、
からなる自動追跡可能なランドマーク選択ステップが実行されることを特徴とする登録方法。
【請求項2】
もし前記目標ピクセル又はボクセルが追跡すべき有効なランドマークであった場合に、定義するための関数が、前記目標ピクセル又はボクセルの特徴パラメータの各一つ又は組み合わせの値と閾値とを比較することであることを特徴とする請求項1記載の方法。
【請求項3】
前記各特徴パラメータに対する閾値を求め、前記各特徴パラメータを対応する閾値と比較して、単一の比較結果の関数としての全体的評価パラメータから、有効に追跡可能なランドマークと有効でない追跡可能なランドマークの質を求めることを特徴とする請求項2に記載の方法。
【請求項4】
前記各特徴パラメータに対して閾値の関数としての全体的閾値が求められ、前記全体的評価パラメータと前記全体的閾値とを比較することによって、有効に追跡可能なランドマークと有効でない追跡可能なランドマークの質を求めることを特徴とする請求項3に記載の方法。
【請求項5】
有効に追跡可能なランドマークか否かに関して比較した各結果を数値で表現するために可変有効性が定義され、
部分又は各特性と対応する閾値との比較の前記可変有効性の値の関数として、全体的評価パラメータが求められ、そして
全体的な評価パラメータの関数として、有効に追跡すべきランドマークと有効でない追跡すべきランドマークの質が求められる
ことを特徴とする請求項3又は4に記載の方法。
【請求項6】
包括的評価パラメータから有効に追跡可能なランドマークと有効でない追跡可能なランドマークを求めるための関数が、前記包括的評価パラメータと包括的閾値との比較であることを特徴とする請求項5に記載の方法。
【請求項7】
特徴パラメータセットが一定数のサブセットに細分化されて、各サブセットに対し二次特徴パラメータがそのサブセットの特徴パラメータの関数として求められ、そして、
各二次特徴パラメータに対し閾値が定義されて、各二次特徴パラメータに対し有効に追跡可能なランドマークと有効でない追跡可能なランドマークの質が、前記二次特徴パラメータと対応する閾値との関数によって求められる
ことを特徴とする請求項1ないし6のいずれか1項に記載の方法。
【請求項8】
各二次特徴パラメータと対応する閾値が、特徴パラメータの対応するサブセットの閾値の関数として求められることを特徴とする請求項7に記載の方法。
【請求項9】
前記サブセットの特徴パラメータが、少なくとも共通の特徴を有し、
特に前記目標ピクセル又はボクセル又は前記ピクセル又はボクセル窓のピクセル又はボクセルのそれぞれ又は一部の外観を記載する数値的外観パラメータからなり、そして
前記窓のピクセル又はボクセルを表現する数値行列の一つ以上の特徴パラメータの同じ関数によって、又は前記数値行列の同じ変換によって求められる
ことを特徴とする請求項7又は8に記載の方法。
【請求項10】
前記包括的評価パラメータが二次特徴パラメータの関数として定義され、そして、
有効に追跡可能なランドマークと有効でない追跡可能なランドマークの質が、前記包括的評価パラメータの関数と、前記二次特徴パラメータに対する閾値の関数として存する包括的閾値と、によって求められることを特徴とする請求項7ないし10のいずれか1項に記載の方法。
【請求項11】
有効に追跡可能なランドマークと有効でない追跡可能なランドマークの質を求めるために、各特徴パラメータ及び/又は各二次特徴パラメータ及び/又は各包括的閾値が加重され、そして、その加重が特徴パラメータ又は二次特徴パラメータの関連性に従って求められることを特徴とする請求項2ないし10のいずれか1項に記載の方法。
【請求項12】
前記包括的評価パラメータ及び/又は前記二次特徴パラメータ及び/又は包括的閾値を求めるための関数が、一つ以上の特徴パラメータの係数の線形結合か、二乗値の合計か、前記特徴パラメータ又は前記第二パラメータ又は前記閾値が処理される非線形関数に存することを特徴とする請求項1ないし11のいずれか1項に記載の方法。
【請求項13】
前記包括的評価パラメータ及び/又は前記二次特徴パラメータ及び/又は前記包括的閾値を求めるための関数が、各特徴パラメータ又は第二パラメータ又は閾値を加重するためにも用いられることを特徴とする請求項12に記載の方法。
【請求項14】
前記目標ピクセル又はボクセルの有効に追跡可能なランドマークと有効でない追跡可能なランドマークの質が、分類アルゴリズムを用いて画像の各ピクセル又はボクセルの特徴パラメータを処理することによって求められることを特徴とする請求項1に記載の方法。
【請求項15】
前記各目標ピクセルのIDが、ベクトル処理されるためにコード化され、前記ベクトルの成分が、前記ピクセル又はボクセルの一つ以上の特徴パラメータに存することを特徴とする請求項14に記載の方法。
【請求項16】
特徴パラメータセットが一定数のサブセットに細分化され、
前記特徴パラメータの代わりに、二次特徴パラメータが前記分類アルゴリズムによって処理され、そして、前記二次特徴パラメータが、特徴パラメータのサブセットの関数として求められることを特徴とする請求項14又は15に記載の方法。
【請求項17】
各ピクセル又はボクセルに対して前記分類アルゴリズムによって処理された入力ベクトル成分が、各特徴パラメータと対応する閾値との比較結果、又は各二次特徴パラメータと対応する閾値との比較結果を数値的に表現する可変有効性に存することを特徴とする請求項14又は15に記載の方法。
【請求項18】
前記分類アルゴリズムが、クラスタ化アルゴリズム又は予測アルゴリズムであることを特徴とする請求項14ないし17のいずれか1項に記載の方法。
【請求項19】
前記分類アルゴリズムが、人工の神経ネットワークであることを特徴とする請求項19記載の方法。
【請求項20】
一定数の有効なランドマークが、有効に追跡可能なランドマークとして識別された既知の場合の一定数の画像を提供するステップと、
上記ステップB1)とB2)に存する自動追跡可能ランドマーク選択ステップを適用することによって、一定数の画像の中で有効に追跡可能であると識別された前記ランドマークに対応して前記ピクセル又はボクセルに対する特徴パラメータセットを求めるステップと、
有効に追跡可能であると識別された各ランドマークに一意的に関連付けられ、成分として前記有効に追跡可能なランドマークと一致する前記ピクセル又はボクセルの特徴パラメータを備えたベクトルを生成するステップと、
変数の所定の数値によって有効に追跡可能な質を表現して、前記数値を、有効に追跡可能なランドマークと一致する各ピクセル又はボクセルをコード化するベクトルに関連付けるステップと、
有効に追跡可能なランドマークと一致する各ピクセル又はボクセルをコード化する各ベクトルが、分類アルゴリズムのためのトレーニングデータベースの登録を形成するステップと、
前記データベースによって分類アルゴリズムをトレーニングするステップと、
鍛えられた分類アルゴリズムの入力に前記目標ピクセル又はボクセルの前記特徴パラメータを備えたベクトルを提供することによって、目標ピクセル又はボクセルの有効に追跡可能なランドマークの質を求めるステップと、
からなることを特徴とする請求項14ないし19のいずれか1項に記載の方法。
【請求項21】
前記トレーニングデータベースが、有効に追跡可能ではないと知られたランドマークと一致する画像のピクセル又はボクセルのベクトル形式に組織化され、
有効に追跡可能でないランドマークの質を、追跡可能なランドマークの質を表現する前記変数の値とは異なる変数の所定の数値によって表現する特徴パラメータが備えられ、
一方で有効に追跡可能でないランドマークの質を表現する前記数値が、既知の有効に追跡可能でないランドマークと一致した各ピクセル又はボクセルをコード化するベクトルに加えられたことを特徴とする請求項20に記載の方法。
【請求項22】
画像のピクセル又はボクセルが、いわゆるLucas & Kanade自動特徴追跡法又はアルゴリズムによって、事前に又は平行して処理されることを特徴とする請求項1ないし22のいずれか1項に記載の方法。
【請求項23】
前記Lucas & Kanadeアルゴリズムが、
a)ボクセルの3次元アレイによって形成された第一3次元画像ボリュームIを提供するステップと、
b)目標ボクセルを取り囲み、例えば前記目標ボクセルを中心とした3×3×3アレイのような調整可能なサイズを有する隣接ボクセルを各目標ボクセルに対して定義するステップと、
c)デカルト座標系を定義し、前記デカルト座標系の各軸に対し、隣接ボクセルに対して第一画像ボリュームの各ボクセルの強度勾配によって形成されるいわゆる勾配ボリュームを計算するステップと、
ここで傾斜ボリュームは次式で定義され、
ア)x方向の勾配ボリューム:
【数1】
イ)y方向の勾配ボリューム:
【数2】
ウ)z方向の勾配ボリューム:
【数3】
d)各目標ボクセルに対し、次式で定義されるいわゆる勾配行列を計算するステップと、
【数4】
e)前記3次元画像ボリュームの各目標ボクセルの各勾配行列に対して、以下の等式を適用することで最小固有値?mを計算するステップと、
【数5】
ここで
【数6】
に対して、前記勾配行列の最小固有値を計算する
【数7】
f)前記最小固有値?mの閾値を定義するステップと、
g)以下の基準を満足する、対応する勾配行列が最小固有値?mを有する画像Iの追跡すべき特徴としての各ボクセルを選択するステップと、
i)?mが前記閾値より大きく、
ii)?mが固有値の全ての最小値?mの最大値の割合より小さくなく、
iii)もし特徴として選択された別のボクセルが、選択された特徴として維持されている目標ボクセル回りの定義された隣接ボクセルに存在する場合、勾配行列がより大きな固有値の最小値?mを有するボクセルの一つのみが特徴として選択され、他の一つは追跡可能な特徴のリストから廃棄され、
iv)特徴として選択されたボクセルの回りの3次元隣接ボクセルの平均信号強度が、調整
可能な平均信号強度閾値より大きい、
以上の各ステップを備えることを特徴とする請求項22に記載の方法。
【請求項24】
特徴として選択されたボクセル回りの隣接ボクセルの平均信号強度に対する閾値が約10に設定されることを特徴とする請求項23に記載の方法。
【請求項25】
前記特徴パラメータが、選択窓のピクセル又はボクセルの強度値に存し、また目標ピクセル又はボクセルの強度値に存することを特徴とする請求項1ないし24のいずれか1項に記載の方法。
【請求項26】
前記特徴パラメータが、選択された窓のピクセル又はボクセルの画像データを備えた数値行列の特異値に代替的に又は付加的に存することを特徴とする請求項1ないし25のいずれか1項に記載の方法。
【請求項27】
前記特徴パラメータが、選択された窓のピクセル又はボクセルを表現する数値行列の勾配行列の固有値に代替的に又は付加的に存することを特徴とする請求項1ないし26のいずれか1項に記載の方法。
【請求項28】
前記特徴パラメータが、選択された窓のピクセル又はボクセルを表現する数値行列のヘッシアン行列の固有値に代替的に又は付加的に存することを特徴とする請求項1ないし27のいずれか1項に記載の方法。
【請求項29】
前記特徴パラメータが、一つ以上のウエーブレットに基づく関数を用いて前記数値行列を処理することで得られた前記窓のピクセル又はボクセルを表現する数値行列のウエーブレット変換係数の一つ以上又はその組み合わせに代替的に又は付加的に存することを特徴とする請求項1ないし28のいずれか1項に記載の方法。
【請求項30】
前記特徴パラメータが、いわゆる行列の同時生起変換に代替的に又は付加的に存することを特徴とする請求項1ないし29のいずれか1項に記載の方法。
【請求項31】
前記特徴パラメータが、前記窓のピクセル又はボクセルを表現する数値行列の自己相関変換の一つ以上の係数に代替的に又は付加的に存することを特徴とする請求項1ないし30のいずれか1項に記載の方法。
【請求項32】
前記特徴パラメータが、
前記窓のピクセル又はボクセルを表現する数値行列の固有値又は特異値及び/又は、
前記窓のピクセル又はボクセルを表現する数値行列の勾配行列又はヘッシアン行列の固有値及び/又は、
ウェーブレット変換係数の一つ以上の係数及び/又は、
自己相関変換の一つ以上の係数及び/又は、
前記窓のピクセル又はボクセルを表現する数値行列の同時生起行列の特異値
の組み合わせに存することを特徴とする請求項1ないし31のいずれか1項に記載の方法。
【請求項33】
同一対象の少なくとも二つ以上の画像が得られるか、又は異なる時間に得られた少なくとも二つ以上の画像が提供され、
第一画像ボリュームIにおいて有効に追跡可能なランドマークと一致して区別された各ボクセルに対して、前記特徴が、第二の遅い時間に得られた同一対象の追加の画像Jのボクセルアレイにおけるその位置に対して追跡され、
前記追跡が、いわゆる「Lucas&Kanade特徴追跡」のピラミッド形の実施によって行われることを特徴とする請求項1ないし32のいずれか1項に記載の方法。
【請求項34】
前記いわゆる「Lucas&Kanade特徴追跡」のピラミッド形の実施が、
画像Iにおける特徴として選択されたボクセルの座標をそれぞれ有する、画像Iのボクセルの3次元アレイに対応する画像ボリュームIにおけるポイントとしてのポイントuと、画像Jのボクセルの3次元アレイに対応する画像ボリュームJにおける対応するポイントvとを定義するステップと、
元のボリュームを表現する下部構造(L=0)と共に、レベルL 0,...,m におけるボリュームを表現するILとJLとを用いて画像ボリュームIとJの二つの理想的なピラミッドを構築するステップと、
最上階Lmにおいて初期化されるいわゆる包括的ピラミッド形の運動ベクトルgを定義するステップと、
【数8】
i)uLとしてポイントuの位置を定義し、以下の式で前記位置を算出し、ここでLは
ピラミッドの実際の階又はレベルの値を有し、
【数9】
ここでpは実際のボリューム位置であり、
ii)ボリューム勾配に対する以下の式に従って、デカルト座標x、y、zの各方向にお
けるボリューム勾配を計算し、
【数10】
iii)以下の等式に従って、勾配行列を計算するために勾配ボリュームを使用し、
【数11】
ここで
【数12】
ここで数値?は、エリアサイズ又は追跡される特徴を表現する目標ボクセルに影響を与える隣接ボクセルを定義し、
iv)各レベルLに対し、次式で定義される反復ベクトルを初期化し、
【数13】
v)k=1から最大カウントKまで、又は最小置換
【数14】
まで、以下のボリューム差を計算し、
【数15】
vi)以下の等式に従って、不整合なベクトルを計算し、
【数16】
vii)以下の式に従って、オプティカルフローベクトルを求め、
【数17】
ここで、G-1はステップiii)で求めた行列Gの逆勾配行列であり、
viii)上で定義されたオプティカルフローベクトルを用いて以下のように反復運動ベクト
ルを計算し、
【数18】
そしてこれをレベルLに対する最後のオプチカルフローと等しくし、
【数19】
ix)最後のレベルL=0が次式で定義される最終オプチカルフローベクトルに達するま
でステップi)からviii)を各レベルに対し繰り返し、
【数20】
x)以下の式でボリュームにおける追跡すべき特徴を表現するポイントuに対応するボ
リュームJにおけるポイントvの座標を求め、
v=u+d
xi)画像Iにおいて選択された特徴に対応する各ボクセルに対して上記方法を繰り返
す、
上記各ステップを備えることを特徴とする請求項33記載の方法。
【請求項35】
画像Iにおける選択された特徴に対応するボクセルのポイントUに対応する特徴を表現するボクセルに対応すると認定された画像Jにおける各ポイントvに反転オプティカルフローベクトルを適用するか、又は
画像Iにおいて選択された特徴であると認定されたボクセルにそれぞれ対応するポイントuにオプティカルフローベクトルを適用することによって、
二つの画像IとJを登録することを特徴とする請求項34に記載の方法。
【請求項36】
a)MRI又は超音波で得られた同一対象の少なくとも第一と第二ディジタル画像又はディジタル処理された画像又は断面画像セットを提供するステップと、
ここで前記画像はピクセル又はボクセルの二次元又は三次元アレイによって形成され、
ここで同一組織又は組織領域又は同一の身体構造上の領域の操作が、前記同一の組織又は組織領域又は前記身体構造上の領域に造影剤が存在する状態で、二度目又はそれ以降の時間で実施され、
b)ランドマークとして設定されたピクセル又はボクセルの一定数を選択し、追跡すべき前記特徴のリストを生成することによって、第一画像又は画像セット内に一定数のランドマーク、いわゆる特徴を定義するステップと、
c)特徴として選択された各ピクセル又はボクセルに対する第一から第二画像又は画像のセットからオプティカルフローベクトルを求めることによって、第一から第二画像又は画像のセットから特徴として選択された各ピクセル又はボクセルの位置を追跡するステップと、
d)第二画像又は画像セットのピクセル又はボクセルに反転オプティカルフローを適用することによって、第一と第二画像又は画像セットを登録するステップと、
からなる造影剤で強調した画像診断法、特に造影剤で強調したMRI又は超音波画像診断法との組合わせで提供されることを特徴とする請求項33ないし35のいずれか1項に記載の方法。
【請求項37】
前記特徴選択ステップの後で、前記特徴追跡ステップを実行する前に、
B1)第一画像又は第一の断面画像セットの各ピクセル又はボクセルの回りに隣接するピクセル又はボクセルを定義するステップと、
ここで前記隣接するピクセル又はボクセルは、制限された数のピクセル又はボクセルを備え、単一画像の場合に二次元隣接部又はボクセルが選択され、断面画像セットの場合に三次元隣接部又はボクセルが選択され、
B2)各ピクセル又はボクセルに対し、前記隣接部又はボクセルの全ピクセル又はボクセルの平均信号強度を求めるステップと、
B3)平均信号強度閾値を定義するステップと、
B4)ステップB2で各隣接部又はボクセルに対して求められた平均信号強度値と、前記平均信号強度値とを、前記平均信号強度閾値に対して比較するステップと、
B5)前記隣接部又はボクセルの前記平均信号強度値が前記ステップB4での閾値より大きい場合、前記ピクセル又はボクセルを追跡すべき特徴として定義することにより、追跡すべき特徴のリストに加えるステップと、
からなる追加の自動特徴選択ステップを実行することを特徴とする請求項33ないし36のいずれか1項に記載の方法。
【図1】
【図3】
【図4】
【図5−7】
【図8】
【図9−10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図20】
【図21】
【図2】
【図19】
【図3】
【図4】
【図5−7】
【図8】
【図9−10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図20】
【図21】
【図2】
【図19】
【公表番号】特表2009−512527(P2009−512527A)
【公表日】平成21年3月26日(2009.3.26)
【国際特許分類】
【出願番号】特願2008−537088(P2008−537088)
【出願日】平成18年10月24日(2006.10.24)
【国際出願番号】PCT/EP2006/067722
【国際公開番号】WO2007/048794
【国際公開日】平成19年5月3日(2007.5.3)
【出願人】(503427359)ブラッコ イメージング ソチエタ ペル アチオニ (19)
【氏名又は名称原語表記】BRACCO IMAGING S.P.A.
【Fターム(参考)】
【公表日】平成21年3月26日(2009.3.26)
【国際特許分類】
【出願日】平成18年10月24日(2006.10.24)
【国際出願番号】PCT/EP2006/067722
【国際公開番号】WO2007/048794
【国際公開日】平成19年5月3日(2007.5.3)
【出願人】(503427359)ブラッコ イメージング ソチエタ ペル アチオニ (19)
【氏名又は名称原語表記】BRACCO IMAGING S.P.A.
【Fターム(参考)】
[ Back to top ]