JP2009515585A - ボリュメトリック位置合わせのための方法、システムおよびコンピュータプログラム - Google Patents
ボリュメトリック位置合わせのための方法、システムおよびコンピュータプログラム Download PDFInfo
- Publication number
- JP2009515585A JP2009515585A JP2008539557A JP2008539557A JP2009515585A JP 2009515585 A JP2009515585 A JP 2009515585A JP 2008539557 A JP2008539557 A JP 2008539557A JP 2008539557 A JP2008539557 A JP 2008539557A JP 2009515585 A JP2009515585 A JP 2009515585A
- Authority
- JP
- Japan
- Prior art keywords
- image
- parameter set
- floating image
- limited
- floating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000004590 computer program Methods 0.000 title claims abstract description 20
- 238000012549 training Methods 0.000 claims abstract description 45
- 230000009466 transformation Effects 0.000 claims abstract description 37
- 238000005457 optimization Methods 0.000 claims description 39
- 238000009826 distribution Methods 0.000 claims description 21
- 238000005070 sampling Methods 0.000 claims description 10
- 238000012545 processing Methods 0.000 claims description 6
- 238000011156 evaluation Methods 0.000 claims description 3
- 201000010099 disease Diseases 0.000 abstract description 5
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 abstract description 5
- 210000004556 brain Anatomy 0.000 description 6
- 230000007170 pathology Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 2
- 238000013138 pruning Methods 0.000 description 2
- 244000141353 Prunus domestica Species 0.000 description 1
- 238000000137 annealing Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/14—Transformations for image registration, e.g. adjusting or mapping for alignment of images
- G06T3/153—Transformations for image registration, e.g. adjusting or mapping for alignment of images using elastic snapping
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/35—Determination of transform parameters for the alignment of images, i.e. image registration using statistical methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/38—Registration of image sequences
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Image Analysis (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本発明は、浮動画像の基準画像とのボリュメトリック位置合わせのための方法に関する。ステップ2'において、浮動画像及び基準画像がアクセスされる。ステップ4及びステップ6において、変換関数T及び類似性関数Sがアクセスされる。本発明による方法は、先験的知識、特にステップ3においてアクセスされる制限されたパラメータ組を使用する。好適には、制限されたパラメータ組は、訓練画像の組の適切なボリュメトリック位置合わせを実施することによって得られる。訓練組は、好適には、各々の臨床アプリケーションについて、浮動画像及び基準画像のシーケンスを含む。同様に、訓練組は、疾患、年齢、性別、人種等の特定のグループを表す患者グループの画像からなりうる。本発明は更に、ボリュメトリック位置合わせを可能にするシステム及びコンピュータプログラムに関する。
Description
本発明は、浮動画像の基準画像とのボリュメトリックなレジストレーション(以下、位置合わせとする)のための方法であって、
−浮動画像及び基準画像にアクセスするステップと、
−浮動画像を空間的にワープする(ゆがめる)ためのパラメータ化された幾何学的変換関数を選択するステップと、
−ワープされた浮動画像及び基準画像の間の類似性を量的に評価するための類似性関数を選択するステップと、
を含む方法に関する。
−浮動画像及び基準画像にアクセスするステップと、
−浮動画像を空間的にワープする(ゆがめる)ためのパラメータ化された幾何学的変換関数を選択するステップと、
−ワープされた浮動画像及び基準画像の間の類似性を量的に評価するための類似性関数を選択するステップと、
を含む方法に関する。
本発明は更に、浮動画像の基準画像とのボリュメトリックな位置合わせのためのシステムに関する。
本発明は更に、浮動画像の基準画像とのボリュメトリック画像位置合わせのためのコンピュータプログラムであって、
−浮動画像及び基準画像にアクセスするステップと、
−浮動画像を空間的にワープするためのパラメータ化された幾何学的変換関数を選択するステップと、
−ワープされた浮動画像及び基準画像の間の類似性を量的に評価するための類似性関数を選択するステップと、
をプロセッサに実行させる命令を含む、コンピュータプログラムに関する。
−浮動画像及び基準画像にアクセスするステップと、
−浮動画像を空間的にワープするためのパラメータ化された幾何学的変換関数を選択するステップと、
−ワープされた浮動画像及び基準画像の間の類似性を量的に評価するための類似性関数を選択するステップと、
をプロセッサに実行させる命令を含む、コンピュータプログラムに関する。
位置合わせ方法はそれ自体、米国特許第5,633,951号明細書から、表面に基づく方法の形で知られている。冒頭の段落に示されている方法の一実施例は、文献「Adaptive Search Space Scaling in Digital Image Registration」(Venkat R. Mandava et al, IEEE Transactions on Medical Imaging, Vol.8, No. 3, September 1989)から知られている。知られている方法は、浮動画像の基準画像との適当な幾何学的な整合を見つける前に、浮動画像をワープすることによって画像のボリュメトリック位置合わせを実行するように構成される。探索空間の剪定を実現可能なものにするための知られている方法において、探索空間は、N次元のハイパーキューブに離散化される。遺伝的アルゴリズム又は擬似焼きなまし法のような最適化方法は、位置合わせの間、ワープされた浮動画像及び基準画像の間の幾何学的整合の品質を表す類似性関数の最適解、特に大域的最大値、を探索するために使用される。知られている方法において、これらの最適化アルゴリズムは、低減された探索空間が、動きアーチファクトの複雑さの或る先験的情報に基づいて人間オペレータによって選択される場合、ボリュメトリック位置合わせのために使用される。代替例として、知られている方法は、適応的な探索空間スケーリングを可能にするように構成され、その場合、人間オペレータが、最初に関心領域を規定する必要がある。
多数のパラメータが調整されなければならず、最適化アルゴリズムの多次元の探索空間をもたらすことは、知られている方法の不利益である。この多次元性は、多くの局所的最適解と相関し、このことは、最適化アルゴリズムをトラップし、それが大域的最適解を求めること妨げることがあり、それによってボリュメトリック位置合わせの品質を低減させる。
本発明の目的は、信頼性がありロバストであるボリュメトリック位置合わせのための方法を提供することである。
このために、本発明による方法は、
−先験的知識に基づいて、パラメータ化された幾何学的変換関数のための制限されたパラメータ組を選択するステップと、
−制限されたパラメータ組を有するパラメータ化された幾何学的変換関数を使用して浮動画像を空間的にワープし、ワープされた浮動画像を与えるステップと、
−類似性関数の最適解を見つけるために、制限されたパラメータ組を最適化するステップと、
を含む。
−先験的知識に基づいて、パラメータ化された幾何学的変換関数のための制限されたパラメータ組を選択するステップと、
−制限されたパラメータ組を有するパラメータ化された幾何学的変換関数を使用して浮動画像を空間的にワープし、ワープされた浮動画像を与えるステップと、
−類似性関数の最適解を見つけるために、制限されたパラメータ組を最適化するステップと、
を含む。
本発明による方法は、最適化方法が、パラメータ空間にわたる類似性関数の大域的最適解、特に最大値、を見つけるための効率的なツールを提示し、スピード及びロバストネスに関するそれらの性能が、探索空間及び収縮戦略の規定に依存するという洞察に基づく。増加する数のパラメータは、ワープ変換の柔軟性を増やすだけでなく、探索されなければならない探索空間をも増やし、最適化アルゴリズムが大域的最適解に達することを妨げうる増加する数の局所的最適解をまねくことが理解される。制限されたパラメータ組は、大域的最適解を見つけるために最も重要であるパラメータに、パラメータ組を剪定する。この剪定は、浮動画像及び基準画像を表す訓練画像の組から導き出される先験的知識に基づく。
本発明による方法の一実施例において、制限されたパラメータ組は、浮動画像及び基準画像を表す訓練画像のボリュメトリック位置合わせの結果を解析することによって得られ、前記ボリュメトリック位置合わせは、拡大されたパラメータ組を有するパラメータ化された幾何学的変換関数を使用して実施される。
訓練組は、好適には、各々の臨床アプリケーションについて、浮動画像及び基準画像のシーケンスを含む。例えば、特定の解剖学的部位について又は特定のタイプの病理について、ボリュメトリック位置合わせが首尾よく実施されることができる。同様に、訓練組は、疾患、年齢、性別、人種等の特定のグループを表す患者グループの画像からなりうる。上述のカテゴリのすべての可能な組み合わせも同様に企図される。本発明の方法によれば、予備知識は、先験的に準備された訓練組から導き出される。
この技術的な方策のため、以下の有利な効果が達成される。すなわち、最適化アルゴリズム、特に確率論的な探索戦略又はランダムにリスタートされる決定論的最適化方法が、考察中の画像のタイプに対応する類似性関数の大域的最適解を見つけるために自動的に同調される。例えば特定の解剖学的部位の医用画像のような、対象アプリケーションを表すk個の訓練画像の組は、拡大されたパラメータ組を有するパラメータ化された幾何学的変換関数を使用する計算集約的な位置合わせアルゴリズムを使用して、事前に位置合わせされる。幾何学的変換関数は、例えばアプリケーション領域に特徴的なパラメータ化された動きモデルに基づいて、構成されることができる。この処理は自動的に実施されることが好ましい。結果として得られる拡大されたパラメータ組(p1...pn)は、例えば探索空間の新しいパラメタリゼーションを規定するパラメータの制限された組Q=(q1...qm)をもたらす主成分分析(PCA)を受ける。ここで、好適にはm<nである。このようにして、新しい制限されたパラメータ組Qは、十分な精度をもって初期の訓練組の位置合わせのための最も関連する変形パターンを近似するのに十分である。組Qは、こうして、最適化アルゴリズムを低次元のサブ空間を通じて大域的最適解に向かうように導く主探索方向を規定する。こうして、訓練画像の拡大されたパラメータ組から最も重要なパラメータが得られると、これらのパラメータは、訓練組と同じ画像クラスからの臨床的な位置合わせの幾何学的変換関数のために使用される。例えば、アプリケーション「脳」は、訓練組「脳」からの制限されたパラメータ組を使用する等である。最適化関数は、制限されたパラメータ組を最適化し、ゆえに最適化された制限されたパラメータ組はより効果的であり且つ正確である。最適化アルゴリズムのために、任意の局所的な又は大域的な収束最適化方法が適切である。当分野においてそれ自体知られている確率論的最適化方法が選択される場合、付加の利点が達成される。すなわち、制限されたパラメータ組は、大域的最適解の方へ探索を向かわせる戦略を提供する。付加の利点として、探索空間にわたる変形パターンの確率分布は、訓練組から評価されることができる。Qによってスパンされるサブ空間へのこの分布の投影は、Qによって与えられるサブ空間探索方向に沿って最適サンプリング戦略を決定するために使用されることができる最も関連する変形パターンの周辺分布を提供する。従って、パラメタリゼーションQは更に、探索空間の範囲をよりよく定めるとともに、確率論的最適化アルゴリズムとともに適用される場合に確率論的なサンプリング戦略のための最適密度分布を改善することによって、低減された空間次元において機能する本発明による方法の性能を改善する。
本発明による方法の更なる実施例において、制限されたパラメータ組は、拡大されたパラメータ組の解析から導き出されるフィーチャに基づいて得られる。
フィーチャの適切な例は、低減された数のパラメータ、低減された数の座標軸、個々のパラメータの許容レンジ、パラメータの密度分布を含む。
本発明による方法の更なる実施例において、方法は更に、
−フィーチャの実質的なずれを検出するステップと、
−予備知識を更新するステップと、
−更新された予備知識から、更新された制限されたパラメータ組を導き出すステップと、を含む。
−フィーチャの実質的なずれを検出するステップと、
−予備知識を更新するステップと、
−更新された予備知識から、更新された制限されたパラメータ組を導き出すステップと、を含む。
本発明による方法は、有利には、自己適応的である。この技術的方策によれば、ボリュメトリック位置合わせ方法は、浮動画像の変形パターンの十分な変化を自動的に追跡し、シグナリングすることができる。選択された訓練組が、もはや臨床的なケースを表さない場合、実質的なずれの信号が、新しい最新の訓練組を編集し、それに適切な洗練されたボリュメトリック位置合わせ技法を受けさせ、臨床的なケースをより正確に位置合わせするために使用されることができる更新された制限されたパラメータ組を導き出すことによって、先験的知識を更新するために使用される。
本発明による方法の更なる実施例において、方法は、
−サブ空間にわたる浮動画像の変形パターンの期待される確率分布を評価するために予備知識を使用するステップと、
−前記評価からサブ空間内の最適化関数のためのサンプリング戦略を決定するステップと、を含む。
−サブ空間にわたる浮動画像の変形パターンの期待される確率分布を評価するために予備知識を使用するステップと、
−前記評価からサブ空間内の最適化関数のためのサンプリング戦略を決定するステップと、を含む。
すでに説明されたように、制限されたパラメータ組(q1...qm)は、最適化アルゴリズム、特に確率論的最適化を、低次元のサブ空間を通じて大域的最適解のほうへ導く主探索方向を規定する。更に、探索空間にわたる変形パターンの確率分布は、訓練組から評価されることができ、付加の予備知識を与える。この分布のサブ空間(q1...qm)への投影は、(q1...qm)によって与えられるサブ空間探索方向に沿って最適サンプリング戦略を決定するために使用されることができる最も関連する変形パターンの周辺分布を提供する。
本発明によるシステムは、
−浮動画像及び基準画像にアクセスし、
−浮動画像を空間的にワープするためのパラメータ化された幾何学的変換関数を選択し、
−ワープされた浮動画像及び基準画像の間の類似性を量的に評価するための類似性関数を選択し、
−先験的知識に基づいてパラメータ化された幾何学的変換関数のための制限されたパラメータ組を選択する、
ための入力部と、
−制限されたパラメータ組を有するパラメータ化された幾何学的変換関数を使用して浮動画像を空間的にワープし、ワープされた浮動画像を与え、
−類似性関数の最適解を見つけるために制限されたパラメータ組を最適化する、
ための処理手段と、を有する。
−浮動画像及び基準画像にアクセスし、
−浮動画像を空間的にワープするためのパラメータ化された幾何学的変換関数を選択し、
−ワープされた浮動画像及び基準画像の間の類似性を量的に評価するための類似性関数を選択し、
−先験的知識に基づいてパラメータ化された幾何学的変換関数のための制限されたパラメータ組を選択する、
ための入力部と、
−制限されたパラメータ組を有するパラメータ化された幾何学的変換関数を使用して浮動画像を空間的にワープし、ワープされた浮動画像を与え、
−類似性関数の最適解を見つけるために制限されたパラメータ組を最適化する、
ための処理手段と、を有する。
本発明によるシステムは、有利には、先験的知識が使用される画像のボリュメトリック位置合わせを可能にする。好適には、制限されたパラメータ組は、浮動画像及び基準画像を表す訓練画像のボリュメトリック位置合わせの結果を解析することによって、処理手段によって得られる。前記ボリュメトリック位置合わせは、拡大されたパラメータ組を有するパラメータ化された幾何学的変換関数を使用して実施される。訓練組は、各々の臨床アプリケーションについて、浮動画像及び基準画像のシーケンスを含むことが好ましい。例えば、特定の解剖学的部位について又は特定のタイプの病理について、ボリュメトリック位置合わせが首尾よく実施されることができる。同様に、訓練組は、疾患、年齢、性別、人種等の特定のグループを表す患者グループの画像からなりうる。上述のカテゴリのすべての可能な組み合わせが同様に企図される。更に好適には、本発明によるシステムの処理手段は、拡大されたパラメータ組の解析から導き出されるフィーチャに基づいて、制限されたパラメータ組を得るようにも構成される。適切なフィーチャの例は、低減された数のパラメータ、個々のパラメータの許容レンジ、パラメータの密度分布を含む。更に好適には、本発明によるシステムの処理手段は、フィーチャの実質的なずれを検出し、先験的知識を更新し、更新された先験的知識から更新された制限されたパラメータ組を導き出すようにも構成される。このようにして、訓練組がその表現の実際のデータとの整合性に基づいて更新される自己学習システムが、提供される。本発明によるシステムの更なる利点は図2に記述される。
本発明によるコンピュータプログラムは、
−先験的知識に基づいて、パラメータ化された幾何学的変換関数のための制限されたパラメータ組を選択するステップと、
−制限されたパラメータ組を有するパラメータ化された幾何学的変換関数を使用して浮動画像を空間的にワープし、ワープされた浮動画像を与えるステップと、
−類似性関数の最適解を見つけるための制限されたパラメータ組を最適化するステップと、
をプロセッサに実行させるための命令を含む。
−先験的知識に基づいて、パラメータ化された幾何学的変換関数のための制限されたパラメータ組を選択するステップと、
−制限されたパラメータ組を有するパラメータ化された幾何学的変換関数を使用して浮動画像を空間的にワープし、ワープされた浮動画像を与えるステップと、
−類似性関数の最適解を見つけるための制限されたパラメータ組を最適化するステップと、
をプロセッサに実行させるための命令を含む。
本発明のコンピュータプログラムによれば、先験的知識が、好適には、浮動画像及び基準画像を表す訓練画像のボリュメトリック位置合わせの結果の解析を含み、前記ボリュメトリック位置合わせが、拡大されたパラメータ組を有するパラメータ化された幾何学的変換関数を使用して実施される、改善されたボリュメトリック位置合わせが可能にされる。本発明によるシステムの更なる有利な実施形態は、請求項13−15に与えられる。本発明によるコンピュータプログラムの処理は、図3を参照してより詳しく述べられる。
本発明のこれら及び他の詳細は、図面を参照して更に詳しく記述される。
図1は、本発明による方法のブロックスキームの一実施例を概略的に示している。本発明による方法1のステップ2'において、入力処理が実行される。例えば、ステップ2において、ボリュメトリックに位置合わせされることが企図される2つの浮動画像及び基準画像が、ロードされる。ステップ4において、浮動画像(F)を空間的にワープするための変換関数Tが、アクセスされる。ステップ6において、ワープされた浮動画像(F')及び基準画像(R)の間の類似性を量的に評価するための類似性関数(S)が、アクセスされる。本発明による方法は、先験的知識、特に制限されたパラメータ組、を使用する。制限されたパラメータ組は、ステップ3においてアクセスされる。好適には、制限されたパラメータ組は、訓練画像の組の適切なボリュメトリック位置合わせを実施することによって得られる。訓練組は、好適には、各々の臨床アプリケーションについて、浮動画像及び基準画像のシーケンスを含む。例えば、特定の解剖学的部位について又は特定のタイプの病理について、ボリュメトリック位置合わせは首尾よく実施されることができる。同様に、訓練組は、疾患、年齢、性別、人種等の特定のグループを表す患者グループの画像から構成されうる。上述のカテゴリのすべての可能な組み合わせが同様に企図される。本発明の方法によれば、予備知識は、先験的に準備された訓練組から導き出される。
この技術的な方策のため、以下の有利な効果が達成される。すなわち、最適化アルゴリズム、特に確率論的探索戦略又はランダムにリスタートされる決定論的最適化方法が、考察中の画像のタイプに対応する類似性関数の大域的最適解を見つけるために自動的に調整される。例えば特定の解剖学的部位の医用画像のような、対象アプリケーションを表すk個の訓練画像の組は、拡大されたパラメータ組を有するパラメータ化された幾何学的変換関数を使用する計算集約的な位置合わせアルゴリズムを使用して、事前に位置合わせされる。幾何学的変換関数は、例えば、アプリケーション領域に特徴的なパラメータ化された動きモデルに基づいて、構成されることができる。好適には、この処理は自動的に実施される。結果として得られる拡大されたパラメータ組(p1...pn)は、探索空間の新しいパラメタリゼーションを規定するパラメータQ=(q1...qm)の制限された組をもたらす主成分分析(PCA)を受けさせられる。ここで、好適にはm<nである。こうして、新しい制限されたパラメータ組Qは、十分な精度をもって初期の訓練組の位置合わせのための最も関連する変形パターンを近似するに十分である。こうして、組Qは、最適化アルゴリズムを低次元のサブ空間を通じて大域的最適解のほうへ導く主探索方向を規定する。こうして、訓練画像の拡大されたパラメータ組から最も重要なパラメータが得られると、これらのパラメータは、訓練組と同じ画像のクラスからの臨床的な位置合わせの幾何学的変換関数のために使用される。例えば、アプリケーション「脳」は、訓練組「脳」の制限されたパラメータ組を使用する等である。最適化関数は、制限されたパラメータ組を最適化し、ゆえに最適化された制限されたパラメータ組は、より効果的であり且つ正確である。好適には、訓練画像は、浮動及び基準画像に関して同じデータ取得ユニットを使用してステップ7において取得される。
本発明による方法1は、次のように動作する。ステップ8において、浮動画像(F)は、ステップ3において得られた制限されたパラメータ組を有する変換関数を使用してワープされる。ステップ10において、制限されたパラメータ組は、類似性関数12の大域的最適解を見つける目的で最適化される。類似性関数12は、ワープされた浮動画像14及び基準画像16の間の類似性の程度を評価する。類似性のこの程度が品質基準を満たす場合、最適化は、ステップ17において止められ、方法は、次のステップ19へ進む。他の場合、制限されたパラメータ組は、ステップ10において更に最適化され、ループ10−12−17が続く。最適化アルゴリズムのために、いかなる局所的に又は大域的に収束する最適化方法も適切である。当分野においてそれ自体知られている確率論的最適化方法が選択される場合、付加の利点が達成される。すなわち、制限されたパラメータ組は、大域的最適解の方へ探索を導く戦略を提供する。付加な利点として、探索空間にわたる変形パターンの確率分布が、訓練組から評価されることができる。Qによってスパンされるサブ空間へのこの分布の投影は、Qによって提供されるサブ空間探索方向に沿って最適サンプリング戦略を決定するために使用されることができる最も関連する変形パターンの周辺分布を提供する。従って、パラメタリゼーションQは更に、探索空間の範囲をよりよく定め、確率論的最適化アルゴリズムとともに適用される場合に確率論的サンプリング戦略のための最適密度分布を改善することによって、低減された空間次元で動作する本発明による方法の性能を改善する。
図2は、本発明によるシステムの一実施例を示している。本発明によるシステム20は、浮動画像、基準画像、及び図1を参照して議論された方法に従って得られる制限されたパラメータ組(図示せず)にアクセスするように構成される入力部22を具えるコンピュータ23を有する。入力部22は更に、浮動画像を空間的にワープするためのパラメータ化された幾何学的変換関数にアクセスするように構成される。入力部22は更に、ワープされた入力画像及び基準画像の間の類似性を量的に評価するための類似性関数にアクセスするように構成される。入力部22は更に、制限されたパラメータ組にわたる類似性関数の大域的最適解を見つけるための適切な最適化アルゴリズムにアクセスするように構成される。本発明によるシステム20のコンピュータ23は更に、類似性関数の大域的最適解を見つけるために、制限されたパラメータ組を最適化するように構成されるコンピュータ手段24を有する。本発明によるシステムの動作は、好適には、図1を参照して説明されたような方法1のステップをプロセッサ26に実行させる命令を含むコンピュータプログラムによって制御される。本発明によるシステム20は更に、好適には、少なくとも基準画像を取得するように構成されるデータ取得ユニット21を有する。データ取得ユニットの適切な例は、MRユニット、CTユニット又はX線ユニット、超音波装置等を含む。システム20は更に、好適には、適切なディスプレイ28に投影されるボリュメトリック位置合わせ28の結果を表示するように構成されるビューワ25を有する。データ取得ユニット21、コンピュータ23及び表示ユニット15は、互いにかなり離れて位置してもよく、この場合、好適にはインターネットのようなデータ伝送リンクによって又は適切なワイヤレスデータ通信によって接続される。
図3は、本発明によるコンピュータプログラムのフローチャートの一実施例を概略的に示している。本発明によるコンピュータプログラム31の命令32'に従って、入力処理が実行される。例えば、命令32は、プロセッサ(図示せず)に、ボリュメトリックに位置合わせされることが企図される浮動画像及び基準画像をロードさせる。命令34に従って、浮動画像(F)を空間的にワープするための変換関数Tが、アクセスされ、命令36に従って、ワープされた浮動画像(F')及び基準画像(R)の間の類似性を量的に評価するための類似性関数(S)がアクセスされる。本発明によるコンピュータプログラムは、先験的知識、特に制限されたパラメータ組を使用する。制限されたパラメータ組は、命令33に従ってアクセスされる。好適には、制限されたパラメータ組は、一組の訓練画像の適切なボリュメトリック位置合わせを実施することによって得られる。訓練組は、好適には、各々の臨床アプリケーションについて、浮動画像及び基準画像のシーケンスを含む。例えば、特定の解剖学的部位について又は特定のタイプの病理について、ボリュメトリック位置合わせが首尾よく実施されることができる。同様に、訓練組は、疾患、年齢、性別、人種等の特定のグループを表す患者グループの画像からなりうる。上述のカテゴリのすべての可能な組み合わせが同様に企図される。本発明のコンピュータプログラムに従って、先験的知識は、先験的に準備される訓練組から導き出される。
この技術的な方策のため、以下の有利な効果が達成される。すなわち、最適化アルゴリズム、特に確率論的探索戦略又はランダムにリスタートされる決定論的最適化方法が、考察中の画像のタイプに対応する類似性関数の大域的最適解を見つけるために自動的に調整される。例えば特定の解剖学的部位の医用画像である、対象アプリケーションを表すk個の訓練画像の組は、拡大されたパラメータ組を有するパラメータ化された幾何学的変換関数を使用する計算集約的な位置合わせアルゴリズムを使用して、事前に位置合わせされる。幾何学的変換関数は、例えばアプリケーション領域に特徴的なパラメータ化された動きモデルに基づいて構成されることができる。好適には、この処理は自動的に実施される。結果として得られる拡大されたパラメータ組(p1...pn)は、例えば、探索空間の新しいパラメタリゼーションを規定するパラメータの制限された組Q=(q1...qm)をもたらす主成分分析(PCA)を受けさせられる。ここで、好適にはm<nである。このようにして、新しい制限されたパラメータ組Qは、十分な精度をもって初期の訓練組の位置合わせのための最も関連する変形パターンを近似するに十分である。組Qは、こうして、低次元のサブ空間を通じて大域的最適解のほうへ最適化アルゴリズムを導く主探索方向を規定する。こうして、訓練画像の拡大されたパラメータ組から最も重要なパラメータが得られると、これらのパラメータは、訓練組と同じ画像のクラスからの臨床的な位置合わせの幾何学的変換関数のために使用される。例えば、アプリケーション「脳」は、訓練組「脳」からの制限されたパラメータ組を使用する等である。最適化関数は、制限されたパラメータ組を最適化し、ゆえに最適化された制限されたパラメータ組は、より効果的であり且つ正確である。好適には、訓練画像は、浮動画像及び基準画像に関して同じデータ取得ユニットを使用して適切な命令37によって取得される。
本発明によるコンピュータプログラム31は、次のように動作する。命令38は、命令33により得られた制限されたパラメータ組を有する変換関数を使用して浮動画像(F)をワープする。命令40は、ワープされた浮動画像44及び基準画像46の間の類似性の程度を評価する類似性関数42の大域的最適解を見つける目的で、プロセッサ(図示せず)に制限されたパラメータ組を最適化させる。類似性のこの程度が、品質基準をみたす場合、最適化は、命令47に従って止められ、コンピュータプログラムは、次の命令49へ進む。他の場合、制限されたパラメータ組は、命令40に戻ることによって更に最適化され、ループ40−42−47が続く。最適化アルゴリズムのために、いかなる局所的に又は大域的に収束する最適化方法も適切である。当分野においてそれ自体知られている確率論的最適化方法が選択される場合、付加の利点が達成される。すなわち、制限されたパラメータ組は、大域的最適解の方へ探索を導く戦略を提供する。付加の利点として、探索空間にわたる変形パターンの確率分布が、訓練組から評価されることができる。Qによってスパンされるサブ空間へのこの分布の投影は、Qによって提供されるサブ空間探索方向に沿って最適サンプリング戦略を決定するために使用されることができる最も関連する変形パターンの周辺分布を提供する。従って、パラメタリゼーションQは更に、低減された空間次元において動作し、従って、探索空間の範囲をよりよく定めるとともに、確率論的最適化アルゴリズムとともに適用される場合に確率論的サンプリング戦略のための最適密度分布を改善することによって、ボリュメトリック位置合わせ方法の性能を改善する。
Claims (15)
- 浮動画像の基準画像とのボリュメトリックな位置合わせのための方法であって、
前記浮動画像及び前記基準画像にアクセスするステップと、
前記浮動画像を空間的にワープするためのパラメータ化された幾何学的変換関数を選択するステップと、
ワープされた浮動画像及び前記基準画像の間の類似性を量的に評価するための類似性関数を選択するステップと、
先験的知識に基づいて、前記パラメータ化された幾何学的変換関数のための制限されたパラメータ組を選択するステップと、
前記制限されたパラメータ組を有する前記パラメータ化された幾何学的変換関数を使用して前記浮動画像を空間的にワープし、ワープされた浮動画像を与えるステップと、
前記類似性関数の最適解を見つけるために前記制限されたパラメータ組を最適化するステップと、
を含む方法。 - 前記制限されたパラメータ組は、前記浮動画像及び前記基準画像を表す訓練画像のボリュメトリックな位置合わせの結果を解析することによって得られ、前記ボリュメトリックな位置合わせは、拡大されたパラメータ組を有する前記パラメータ化された幾何学的変換関数を使用して実施される、請求項1に記載の方法。
- 前記制限されたパラメータ組は、前記拡大されたパラメータ組の解析から導き出されるフィーチャに基づいて得られる、請求項2に記載の方法。
- 前記フィーチャが、低減された数のパラメータを含む、請求項3に記載の方法。
- 前記フィーチャが、低減された数の座標軸を含む、請求項3に記載の方法。
- 前記フィーチャが、個々のパラメータの許容レンジを含む、請求項3に記載の方法。
- 前記フィーチャが、前記パラメータの密度分布を含む、請求項3に記載の方法。
- 前記フィーチャの実質的なずれを検出するステップと、
前記先験的知識を更新するステップと、
前記更新された先験的知識から、更新された制限されたパラメータ組を導き出すステップと、
を含む、請求項4乃至請求項7のいずれか1項に記載の方法。 - サブ空間にわたる前記浮動画像の変形パターンの期待される確率分布を評価するために前記先験的知識を使用するステップと、
前記評価から前記サブ空間内の前記最適化関数のためのサンプリング戦略を決定するステップと、
を更に含む、請求項1に記載の方法。 - 浮動画像の基準画像とのボリュメトリックな位置合わせのためのシステムであって、
前記浮動画像及び前記基準画像にアクセスし、
前記浮動画像を空間的にワープするためのパラメータ化された幾何学的変換関数を選択し、
ワープされた浮動画像及び前記基準画像の間の類似性を量的に評価するための類似性関数を選択し、
先験的知識に基づいて、前記パラメータ化された幾何学的変換関数のための制限されたパラメータ組を選択する、
入力部と、
前記制限されたパラメータ組を有する前記パラメータ化された幾何学的変換関数を使用して前記浮動画像を空間的にワープし、ワープされた浮動画像を与え、
前記類似性関数の最適解を見つけるために前記制限されたパラメータ組を最適化する、
処理手段と、
を有するシステム。 - 少なくとも前記基準画像を取得するように構成されるデータ取得ユニットを更に有する、請求項10に記載のシステム。
- 浮動画像の基準画像とのボリュメトリックな位置合わせのためのコンピュータプログラムであって、
前記浮動画像及び前記基準画像にアクセスするステップと、
前記浮動画像を空間的にワープするためのパラメータ化された幾何学的変換関数を選択するステップと、
ワープされた浮動画像及び前記基準画像の間の類似性を量的に評価するための類似性関数を選択するステップと、
先験的知識に基づいて、前記パラメータ化された幾何学的変換関数のための制限されたパラメータ組を選択するステップと、
前記制限されたパラメータ組を有する前記パラメータ化された幾何学的変換関数を使用して前記浮動画像を空間的にワープし、ワープされた浮動画像を与えるステップと、
前記類似性関数の最適解を見つけるために前記制限されたパラメータ組を最適化するステップと、
をプロセッサに実行させる命令を含むコンピュータプログラム。 - 拡大されたパラメータ組の解析から導き出されるフィーチャに基づいて、前記制限されたパラメータ組を得るステップを、前記プロセッサに実行させる命令を更に含む、請求項12に記載のコンピュータプログラム。
- 前記フィーチャの実質的なずれを検出するステップと、
前記先験的知識を更新するステップと、
前記更新された先験的知識から、更新された制限されたパラメータ組を導き出すステップと、
を前記プロセッサに実行させる命令を更に含む、請求項13に記載のコンピュータプログラム。 - サブ空間にわたる前記浮動画像の変形パターンの期待される確率分布を評価するために前記先験的知識を使用するステップと、
前記評価から前記サブ空間内の前記最適化関数のためのサンプリング戦略を決定するステップと、
を前記プロセッサに実行させる命令を更に含む、請求項12乃至請求項14のいずれか1項に記載のコンピュータプログラム。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP05110705 | 2005-11-14 | ||
PCT/IB2006/054078 WO2007054864A1 (en) | 2005-11-14 | 2006-11-03 | A method, a system and a computer program for volumetric registration |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2009515585A true JP2009515585A (ja) | 2009-04-16 |
Family
ID=37772573
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008539557A Withdrawn JP2009515585A (ja) | 2005-11-14 | 2006-11-03 | ボリュメトリック位置合わせのための方法、システムおよびコンピュータプログラム |
Country Status (5)
Country | Link |
---|---|
US (1) | US20080310760A1 (ja) |
EP (1) | EP1952341A1 (ja) |
JP (1) | JP2009515585A (ja) |
CN (1) | CN101310302A (ja) |
WO (1) | WO2007054864A1 (ja) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016002214A (ja) * | 2014-06-16 | 2016-01-12 | キヤノン株式会社 | 画像処理装置、画像処理方法およびプログラム |
JP2019005566A (ja) * | 2017-06-21 | 2019-01-17 | インターナショナル・ビジネス・マシーンズ・コーポレーションInternational Business Machines Corporation | 自動画像レジストレーション品質評価 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8385687B1 (en) * | 2006-12-06 | 2013-02-26 | Matrox Electronic Systems, Ltd. | Methods for determining a transformation between images |
US8588488B2 (en) * | 2009-02-11 | 2013-11-19 | Koninklijke Philips N.V. | Group-wise image registration based on motion model |
WO2010118629A1 (en) * | 2009-04-17 | 2010-10-21 | The Hong Kong University Of Science And Technology | Method, device and system for facilitating motion estimation and compensation of feature-motion decorrelation |
US20130072797A1 (en) * | 2010-05-31 | 2013-03-21 | Samsung Medison Co., Ltd. | 3d ultrasound apparatus and method for operating the same |
BR112013012981A2 (pt) * | 2010-11-26 | 2016-12-06 | Koninkl Philips Electronics Nv | aparelho de processamento de imagem, estação de trabalho, aparelhos de imagem, método de criação e produto de programa de computador |
US20130163854A1 (en) * | 2011-12-23 | 2013-06-27 | Chia-Ming Cheng | Image processing method and associated apparatus |
KR20140105103A (ko) * | 2013-02-21 | 2014-09-01 | 삼성전자주식회사 | 장기의 움직임을 추적하는 방법, 장치 및 의료 영상 시스템 |
CN109785225B (zh) * | 2017-11-13 | 2023-06-16 | 虹软科技股份有限公司 | 一种用于图像矫正的方法和装置 |
CN109785390B (zh) * | 2017-11-13 | 2022-04-01 | 虹软科技股份有限公司 | 一种用于图像矫正的方法和装置 |
CN111260700B (zh) * | 2020-01-09 | 2023-05-30 | 复旦大学 | 一种多参数磁共振图像全自动配准与分割方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0602730B1 (en) * | 1992-12-18 | 2002-06-19 | Koninklijke Philips Electronics N.V. | Registration of Volumetric images which are relatively elastically deformed by matching surfaces |
-
2006
- 2006-11-03 JP JP2008539557A patent/JP2009515585A/ja not_active Withdrawn
- 2006-11-03 WO PCT/IB2006/054078 patent/WO2007054864A1/en active Application Filing
- 2006-11-03 EP EP06821305A patent/EP1952341A1/en not_active Withdrawn
- 2006-11-03 CN CNA2006800423349A patent/CN101310302A/zh active Pending
- 2006-11-03 US US12/092,991 patent/US20080310760A1/en not_active Abandoned
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016002214A (ja) * | 2014-06-16 | 2016-01-12 | キヤノン株式会社 | 画像処理装置、画像処理方法およびプログラム |
JP2019005566A (ja) * | 2017-06-21 | 2019-01-17 | インターナショナル・ビジネス・マシーンズ・コーポレーションInternational Business Machines Corporation | 自動画像レジストレーション品質評価 |
JP7034013B2 (ja) | 2017-06-21 | 2022-03-11 | インターナショナル・ビジネス・マシーンズ・コーポレーション | 自動画像レジストレーション品質評価 |
Also Published As
Publication number | Publication date |
---|---|
CN101310302A (zh) | 2008-11-19 |
US20080310760A1 (en) | 2008-12-18 |
EP1952341A1 (en) | 2008-08-06 |
WO2007054864A1 (en) | 2007-05-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2009515585A (ja) | ボリュメトリック位置合わせのための方法、システムおよびコンピュータプログラム | |
US8437521B2 (en) | Systems and methods for automatic vertebra edge detection, segmentation and identification in 3D imaging | |
US20080317382A1 (en) | Adaptive Point-Based Elastic Image Registration | |
US7450780B2 (en) | Similarity measures | |
US20200334825A1 (en) | Method and apparatus for medical imaging | |
CN113298830B (zh) | 一种基于自监督的急性颅内ich区域图像分割方法 | |
CN110288615B (zh) | 一种基于深度学习的倾斜定位框定位方法 | |
CN115345938B (zh) | 基于全局到局部的头影标记点定位方法、设备及介质 | |
US20060242218A1 (en) | Prior-constrained mean shift analysis | |
CN109949288A (zh) | 肿瘤类型确定系统、方法及存储介质 | |
CN114332132A (zh) | 图像分割方法、装置和计算机设备 | |
US7983464B2 (en) | System and method for corpus callosum segmentation in magnetic resonance images | |
CN113066111B (zh) | 一种基于ct影像的心脏二尖瓣顶点的自动定位方法 | |
CN116113986A (zh) | 用于病理器官分割的用户交互快速标注的用户引导域自适应 | |
Pei et al. | Learning spatiotemporal probabilistic atlas of fetal brains with anatomically constrained registration network | |
JP2022509126A (ja) | 対話的な反復的画像注釈付け | |
US20240055104A1 (en) | Method for analyzing output of neural network, and system therefor | |
US20220261985A1 (en) | System for determining the presence of features in a dataset | |
CN114581340A (zh) | 一种图像校正方法及设备 | |
Hoogendoorn et al. | A groupwise mutual information metric for cost efficient selection of a suitable reference in cardiac computational atlas construction | |
US9135401B2 (en) | Information-theoretic view of the scheduling problem in whole-body computer aided detection/diagnosis (CAD) | |
US11721022B2 (en) | Apparatus and method for automated analyses of ultrasound images | |
CN110738231A (zh) | 改进S-DNet神经网络模型的乳腺X射线图像的分类方法 | |
CN113223104B (zh) | 一种基于因果关系的心脏mr图像插补方法及系统 | |
Aranguren et al. | Segmentation of magnetic resonance brain images through the self-adaptive differential evolution algorithm and the minimum cross-entropy criterion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20100105 |