JP6781417B2 - 患者間脳レジストレーション - Google Patents

患者間脳レジストレーション Download PDF

Info

Publication number
JP6781417B2
JP6781417B2 JP2019547070A JP2019547070A JP6781417B2 JP 6781417 B2 JP6781417 B2 JP 6781417B2 JP 2019547070 A JP2019547070 A JP 2019547070A JP 2019547070 A JP2019547070 A JP 2019547070A JP 6781417 B2 JP6781417 B2 JP 6781417B2
Authority
JP
Japan
Prior art keywords
registration
image
diffeomorphic
demon
rigid registration
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.)
Active
Application number
JP2019547070A
Other languages
English (en)
Other versions
JP2019535472A (ja
Inventor
チェン−ルイ チョウ
チェン−ルイ チョウ
ミン−チャン リウ
ミン−チャン リウ
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sony Corp
Original Assignee
Sony Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sony Corp filed Critical Sony Corp
Publication of JP2019535472A publication Critical patent/JP2019535472A/ja
Application granted granted Critical
Publication of JP6781417B2 publication Critical patent/JP6781417B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/32Determination of transform parameters for the alignment of images, i.e. image registration using correlation-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/02Affine transformations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10004Still image; Photographic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30016Brain
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30096Tumor; Lesion

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Medical Treatment And Welfare Office Work (AREA)

Description

本発明は、画像分析に関する。特に、本発明は、腫瘍画像分析に関する。
並列計算技術により、コンピュータは、人口分析によって、脳の異常を発見することができる。患者画像を同一の座標空間内に配置するデータ正規化は、人口分析の成功にとって重要である。しかしながら、現在利用可能な脳人口分析のための正規化技術は、まだ、非常に基本的なものである。脳正規化は、現在、軟部組織の変形や腫瘍の存在の不一致にかかわらず、線形変換を用いて、被変形画像を目標画像にマッピングするアフィンレジストレーションに基づいている。
データ正規化のための患者間脳レジストレーション方法は、腫瘍が存在する場合でも、異なる患者から取得される2つの脳画像を変形可能に位置合わせする。
一態様では、装置の非一時的メモリにプログラムされる方法は、アフィンレジストレーションを行って、被変形画像及び目標画像を、おおよそ同一の座標空間に位置決めするステップであって、前記被変形画像及び前記目標画像は、各々、脳構造を含む、ステップと、腫瘍適応指標を用いた微分同相デモンズ非剛体レジストレーションを実行して、前記脳構造の境界を精細化するステップと、を含む。前記アフィンレジストレーションは、点、直線及び平面を保存することによって、画像を位置合わせすることを含む。前記微分同相デモンズ非剛体レジストレーションは、ベクトル場及びベクトルの大きさを利用する。前記微分同相デモンズ非剛体レジストレーションは、オプティカルフローから導出された変形力及び制御可能な正則化を利用する。前記微分同相デモンズ非剛体レジストレーションは、変形場において、画像強度のマッチング誤差を最小化することを含む。前記腫瘍適応指標は、腫瘍領域の誤差の寄与を除外する。前記微分同相デモンズ非剛体レジストレーションは、レジストレーション時間中、滑らかな変形軌跡を用いることを含む。
別の態様では、装置は、アプリケーションを記憶するための非一時的メモリであって、前記アプリケーションは、アフィンレジストレーションを行って、被変形画像及び目標画像を、おおよそ同一の座標空間に位置決めするステップであって、前記被変形画像及び前記目標画像は、各々、脳構造を含む、ステップと、腫瘍適応指標を用いた微分同相デモンズ非剛体レジストレーションを実行して、前記脳構造の境界を精細化するステップと、を行うためのものである、非一時的メモリと、前記メモリに結合され、前記アプリケーションを処理するように構成される処理構成要素と、を備える。前記アフィンレジストレーションは、点、直線及び平面を保存することによって、画像を位置合わせすることを含む。前記微分同相デモンズ非剛体レジストレーションは、ベクトル場及びベクトルの大きさを利用する。前記微分同相デモンズ非剛体レジストレーションは、オプティカルフローから導出された変形力及び制御可能な正則化を利用する。前記微分同相デモンズ非剛体レジストレーションは、変形場において、画像強度のマッチング誤差を最小化することを含む。前記腫瘍適応指標は、腫瘍領域の誤差の寄与を除外する。前記微分同相デモンズ非剛体レジストレーションは、レジストレーション時間中、滑らかな変形軌跡を用いることを含む。
別の態様では、システムは、磁気共鳴撮像装置と、計算装置であって、前記計算装置は、アフィンレジストレーションを行って、被変形画像及び目標画像を、おおよそ同一の座標空間に位置決めするステップであって、前記被変形画像及び前記目標画像は、各々、脳構造を含む、ステップと、腫瘍適応指標を用いた微分同相デモンズ非剛体レジストレーションを実行して、前記脳構造の境界を精細化するステップと、を行うように構成される、計算装置と、前記メモリに結合され、前記アプリケーションを処理するように構成される処理構成要素と、を備える。前記アフィンレジストレーションは、点、直線及び平面を保存することによって、画像を位置合わせすることを含む。前記微分同相デモンズ非剛体レジストレーションは、ベクトル場及びベクトルの大きさを利用する。前記微分同相デモンズ非剛体レジストレーションは、オプティカルフローから導出された変形力及び制御可能な正則化を利用する。前記微分同相デモンズ非剛体レジストレーションは、変形場において、画像強度のマッチング誤差を最小化することを含む。前記腫瘍適応指標は、腫瘍領域の誤差の寄与を除外する。前記微分同相デモンズ非剛体レジストレーションは、レジストレーション時間中、滑らかな変形軌跡を用いることを含む。
いくつかの実施形態による患者間脳レジストレーション方法を示す図である。 いくつかの実施形態によるデモンズ非剛体レジストレーションを示す図である。 いくつかの実施形態による第1の画像の第2の画像への変形を示す図である。 いくつかの実施形態によるデモンズ非剛体レジストレーションのための腫瘍適応指標を示す図である。 いくつかの実施形態による正規SSD指標と適応SSD指標との間の差を示す図である。 いくつかの実施形態による微分同相デモンズレジストレーションを示す図である。 非剛体レジストレーション後の改善された脳表面及び脳構造の位置合わせの例を示す図である。 非剛体レジストレーション後の改善された脳表面及び脳構造の位置合わせの例を示す図である。 非剛体レジストレーション後の改善された脳表面及び脳構造の位置合わせの例を示す図である。 いくつかの実施形態による患者間脳レジストレーション方法を実行するように構成される例示的な計算装置のブロック図である。 いくつかの実施形態による磁気共鳴撮像(MRI)システムを示す図である。
患者間脳レジストレーション方法は、2つの画像をおおよそ同一の座標空間に配置するアフィンレジストレーションから開始する。次に、この方法は、脳構造の境界を更に精細化する新規の腫瘍適応指標を用いた微分同相デモンズ非剛体レジストレーションを用いる。この腫瘍適応指標は、「自動3D脳腫瘍セグメント化及び分類(AUTOMATIC 3D BRAIN TUMOR SEGMENTATION AND CLASSIFICATION)」という名称の米国特許代理人整理番号Sony−68500に記載の方法を用いてセグメント化される腫瘍位置で実施される。尚、この米国特許代理人整理番号Sony−68500は、全ての目的に対してその全体が引用によりここに組み込まれる。
図1は、いくつかの実施形態による患者間脳レジストレーション方法の図を示す。ステップ100では、アフィンレジストレーションが実行される。アフィンレジストレーションは、点、直線及び平面を保存することによって、画像を位置合わせしようとする。例えば、目標画像及び被変形画像(移動画像又は移動被変形画像と呼ばれる場合もある)が、アフィンレジストレーションを用いて、位置合わせされる。ステップ102では、非剛体レジストレーションが実行される。非剛体レジストレーションは、レジストレーション処理において、ベクトル場及びベクトルの大きさを利用する。非剛体レジストレーションは、また、デモンズレジストレーション、腫瘍適応指標及び微分同相性を用いて、画像の適合を改善する。非剛体レジストレーションは、構造境界の変形精細化を可能にする。例えば、目標画像及び被変形画像が、アフィンレジストレーションを用いて、位置合わせされた後、非剛体レジストレーションは、被変形画像を目標画像に又はその逆に変形することなどによって、位置合わせに改良を加えることができる。いくつかの実施形態では、より少ない又は追加のステップが実行される。例えば、MRI情報を取得するステップ及び/又は照合された画像を表示するステップを含むことができる。いくつかの実施形態では、ステップの順序は変更される。
図2は、いくつかの実施形態によるデモンズ非剛体レジストレーションの図を示す。セット200内の画像は、アフィンレジストレーション後であり、セット202内の画像は、非剛体レジストレーション後である。セット202内の画像は、非剛体レジストレーションを用いて、より良好に位置合わせされる。
図3は、いくつかの実施形態による第1の画像の第2の画像への変形の図を示す。Imoving画像300(被変形画像とも呼ばれる)は、ここで説明されるように、Itarget画像302に変形される。
デモンズ非剛体レジストレーションは、オプティカルフローから導出された変形力及び制御可能な正則化を利用する。オプティカルフローから導出された変形力は、以下のように計算される。
速度は、空間勾配に沿って、不一致の強度に比例する。
制御可能な正則化は、以下の通りである。
ここで、αは、正則化の優先度、すなわち、滑らかさ又はデータの一致を制御する。
は、変形(弾性)のガウシアン平滑化である。
この変形を行うため、マッチング誤差(画像強度の差の二乗和(SSD))が、変形場において、下式のように、最小化される。
SSD = Γx(Itarget(x) - Imoving(x))2
図4は、いくつかの実施形態によるデモンズ非剛体レジストレーションのための腫瘍適応指標の図を示す。腫瘍の存在は、レジストレーションの最適化の誤差に寄与するが、腫瘍領域の誤差の寄与を除外する腫瘍適応誤差指標は、この問題に対処する。その結果は、以下の通りである。
SSDadaptive = Γx0tumor(Itarget(x) - Imoving(x))2
すなわち、最小化は、腫瘍であると判定された領域以外の画像に対して行われる。腫瘍は、画像比較(例えば、腫瘍の無い脳画像を目標画像又は被変形画像と比較し、何らかの有意差がある場合、その領域を腫瘍として表示する)、局所極小検出、及び/又は画像において腫瘍を判定するための他の任意の画像処理技術などのいずれかの方法で判定することができる。Imoving画像400は、ここで説明されるように、Itarget画像402に変形される。
図5は、いくつかの実施形態による正規SSD指標と適応SSD指標との間の差の図を示す。正規SSD指標500は、制限された変形を行う。適応SSD指標502は、腫瘍領域の誤差の寄与を除外することによって、より大きな変形を可能にする。
図6は、いくつかの実施形態による微分同相デモンズレジストレーションの図を示す。Imoving画像600は、ここで説明されるように、レジストレーション時間中、滑らかな変形軌跡を用いて、Itarget画像602に変形される。患者間の解剖学的構造の大きな変化は、大変形を可能にするレジストレーションを用いるので、最終的な変形場(弾性)の代わりに、経時変化する速度場(流体状)を平滑化することによって、「微分同相」(逆一貫性のある滑らかな)変形を計算することによって、変形軌跡を滑らかにすることが重要である。
図7〜図9は、非剛体レジストレーション後の改善された脳表面及び脳構造の位置合わせの例を示す。
図10は、いくつかの実施形態による患者間脳レジストレーション方法を実行するように構成される例示的な計算装置のブロック図を示す。計算装置1000を用いて、画像及びビデオなどの情報を、取得、記憶、計算、処理、通信、及び/又は表示することができる。一般に、計算装置1000を実現するのに適したハードウェア構造は、ネットワークインタフェース1002、メモリ1004、プロセッサ1006、入出力装置1008、バス1010、及び記憶装置1012を含む。プロセッサの選択は、十分な速度を有する好適なプロセッサが選ばれる限り重要ではない。メモリ1004は、当該技術分野で公知の任意の従来のコンピュータメモリとすることができる。記憶装置1012は、ハードドライブ、CDROM、CDRW、DVD、DVDRW、高精細ディスク/ドライブ、超高精細ドライブ、フラッシュメモリカード、又は他の任意の記憶装置を含むことができる。計算装置1000は、1つ又はそれ以上のネットワークインタフェース1002を含むことができる。ネットワークインタフェースの一例は、イーサネット又は他の種類のLANに接続されるネットワークカードを含む。入出力装置1008は、以下のもの、すなわち、キーボード、マウス、モニタ、スクリーン、プリンタ、モデム、タッチスクリーン、ボタンインタフェース、及び他のデバイスのうちの1つ又はそれ以上を含むことができる。患者間脳レジストレーション方法を実行するのに用いられる患者間脳レジストレーションアプリケーション1030は、記憶装置1012及びメモリ1004に記憶され、アプリケーションが通常処理されるように処理される可能性が高い。計算装置1000は、図10に示すより多い又は少ない構成要素を含むことができる。いくつかの実施形態では、患者間脳レジストレーション方法のハードウェア1020が含まれる。図10の計算装置1000は、患者間脳レジストレーション方法のためのアプリケーション1030及びハードウェア1020を含むが、患者間脳レジストレーション方法は、ハードウェア、ファームウェア、ソフトウェア、又はそれらの任意の組み合わせとして、計算装置上で実現することができる。例えば、いくつかの実施形態では、患者間脳レジストレーション方法のアプリケーション1030は、メモリ内にプログラムされ、プロセッサを用いて実行される。別の実施例では、いくつかの実施形態において、患者間脳レジストレーション方法のハードウェア1020は、患者間脳レジストレーション方法を実行するように専用に設計されるゲートを含む、プログラムされたハードウェアロジックである。
いくつかの実施形態では、患者間脳レジストレーション方法のアプリケーション1030は、いくつかのアプリケーション及び/又はモジュールを含む。いくつかの実施形態では、モジュールは、1つ又はそれ以上のサブモジュールも含む。いくつかの実施形態では、より少ない又は追加のモジュールを含むことができる。
好適な計算装置の実施例は、パーソナルコンピュータ、ラップトップコンピュータ、コンピュータワークステーション、サーバ、メインフレームコンピュータ、ハンドヘルドコンピュータ、携帯情報端末、セルラ電話/携帯電話、スマート家電、ゲーム機、デジタルカメラ、デジタルカムコーダ、カメラ付き携帯電話、スマートフォン、携帯音楽プレーヤー、タブレットコンピュータ、移動体デバイス、ビデオプレーヤー、ビデオディスクライター/プレーヤー(例えば、DVDライター/プレーヤー、高精細ディスクライター/プレーヤー、超高精細ディスクライター/プレーヤー)、テレビジョン、家庭用娯楽システム、拡張現実デバイス、仮想現実デバイス、スマートジュエリー(例えば、スマートウォッチ)、又は他の任意の好適な計算装置を含む。
図11は、いくつかの実施形態による磁気共鳴撮像(MRI)システムの図を示す。MRI装置1100を用いて、脳スキャンなどのMRI画像を取得する。計算装置1000は、MRI装置1100から(例えば、クラウドでダウンロード又はアクセスすることによって)MRI情報を受信し、計算装置1000は、患者間脳レジストレーション方法を実行する。
本明細書で説明する患者間脳レジストレーション方法を利用するため、この患者間脳レジストレーション方法を用いて、MRI情報を分析する。この分析に基づいて、脳の異常を識別することができる。識別された異常に基づいて、外科医は、適切な処置をとることができる。
動作時、患者間脳レジストレーション方法は、アフィンレジストレーションと腫瘍適応指標を用いた非剛体レジストレーションとに基づいた脳正規化を用いて、人口分析によって、脳の異常を発見することができる。この改良の結果として、単にアフィンレジストレーションを用いるよりも良好な脳構造の位置合わせが得られる。
本明細書では、脳腫瘍に関連して、患者間脳レジストレーション方法を説明したが、患者間脳レジストレーション方法は、他の任意の腫瘍又は疾患に適合することができる。
患者間脳レジストレーションのいくつかの実施形態
1.装置の非一時的メモリにプログラムされる方法であって、前記方法は、
アフィンレジストレーションを行って、被変形画像及び目標画像を、おおよそ同一の座標空間に位置決めするステップであって、前記被変形画像及び前記目標画像は、各々、脳構造を含む、ステップと、
腫瘍適応指標を用いた微分同相デモンズ非剛体レジストレーションを実行して、前記脳構造の境界を精細化するステップと、
を含むことを特徴とする方法。
2.前記アフィンレジストレーションは、点、直線及び平面を保存することによって、画像を位置合わせすることを含むことを特徴とする、第1項に記載の方法。
3.前記微分同相デモンズ非剛体レジストレーションは、ベクトル場及びベクトルの大きさを利用することを特徴とする、第1項に記載の方法。
4.前記微分同相デモンズ非剛体レジストレーションは、オプティカルフローから導出された変形力及び制御可能な正則化を利用することを特徴とする、第1項に記載の方法。
5.前記微分同相デモンズ非剛体レジストレーションは、変形場において、画像強度のマッチング誤差を最小化することを含むことを特徴とする、第1項に記載の方法。
6.前記腫瘍適応指標は、腫瘍領域の誤差の寄与を除外することを特徴とする、第1項に記載の方法。
7.前記微分同相デモンズ非剛体レジストレーションは、レジストレーション時間中、滑らかな変形軌跡を用いることを含むことを特徴とする、第1項に記載の方法。
8.アプリケーションを記憶するための非一時的メモリであって、前記アプリケーションは、
アフィンレジストレーションを行って、被変形画像及び目標画像を、おおよそ同一の座標空間に位置決めするステップであって、前記被変形画像及び前記目標画像は、各々、脳構造を含む、ステップと、
腫瘍適応指標を用いた微分同相デモンズ非剛体レジストレーションを実行して、前記脳構造の境界を精細化するステップと、を行うためのものである、
非一時的メモリと、
前記メモリに結合され、前記アプリケーションを処理するように構成される処理構成要素と、
を備えることを特徴とする装置。
9.前記アフィンレジストレーションは、点、直線及び平面を保存することによって、画像を位置合わせすることを含むことを特徴とする、第8項に記載の装置。
10.前記微分同相デモンズ非剛体レジストレーションは、ベクトル場及びベクトルの大きさを利用することを特徴とする、第8項に記載の装置。
11.前記微分同相デモンズ非剛体レジストレーションは、オプティカルフローから導出された変形力及び制御可能な正則化を利用することを特徴とする、第8項に記載の装置。
12.前記微分同相デモンズ非剛体レジストレーションは、変形場において、画像強度のマッチング誤差を最小化することを含むことを特徴とする、第8項に記載の装置。
13.前記腫瘍適応指標は、腫瘍領域の誤差の寄与を除外することを特徴とする、第8項に記載の装置。
14.前記微分同相デモンズ非剛体レジストレーションは、レジストレーション時間中、滑らかな変形軌跡を用いることを含むことを特徴とする、第8項に記載の装置。
15.磁気共鳴撮像装置と、
計算装置であって、前記計算装置は、
アフィンレジストレーションを行って、被変形画像及び目標画像を、おおよそ同一の座標空間に位置決めするステップであって、前記被変形画像及び前記目標画像は、各々、脳構造を含む、ステップと、
腫瘍適応指標を用いた微分同相デモンズ非剛体レジストレーションを実行して、前記脳構造の境界を精細化するステップと、を行うように構成される、
計算装置と、
前記メモリに結合され、前記アプリケーションを処理するように構成される処理構成要素と、
を備えることを特徴とするシステム。
16.前記アフィンレジストレーションは、点、直線及び平面を保存することによって、画像を位置合わせすることを含むことを特徴とする、第15項に記載のシステム。
17.前記微分同相デモンズ非剛体レジストレーションは、ベクトル場及びベクトルの大きさを利用することを特徴とする、第15項に記載のシステム。
18.前記微分同相デモンズ非剛体レジストレーションは、オプティカルフローから導出された変形力及び制御可能な正則化を利用することを特徴とする、第15項に記載のシステム。
19.前記微分同相デモンズ非剛体レジストレーションは、変形場において、画像強度のマッチング誤差を最小化することを含むことを特徴とする、第15項に記載のシステム。
20.前記腫瘍適応指標は、腫瘍領域の誤差の寄与を除外することを特徴とする、第15項に記載のシステム。
21.前記微分同相デモンズ非剛体レジストレーションは、レジストレーション時間中、滑らかな変形軌跡を用いることを含むことを特徴とする、第15項に記載のシステム。
本発明の構成及び動作の原理の理解を容易にするために、詳細内容を組み込んだ特定の実施形態に関して本発明を説明してきた。このような本明細書における特定の実施形態及びその詳細内容への言及は、本明細書に添付される特許請求の範囲を限定することを意図するものではない。特許請求の範囲によって規定される本発明の精神及び範囲から逸脱することなく、例示のために選択された実施形態に、他の様々な修正を行うことができることは、当業者に容易に理解されるであろう。
100 ステップ
102 ステップ
200 セット
202 セット
300 Imoving画像
302 Itarget画像
400 Imoving画像
402 Itarget画像
500 正規SSD指標
502 適応SSD指標
600 Imoving画像
602 Itarget画像
1000 計算装置
1002 ネットワークインタフェース
1004 メモリ
1006 プロセッサ
1008 入出力装置
1010 バス
1012 記憶装置
1020 患者間脳レジストレーション方法のハードウェア
1030 患者間脳レジストレーションアプリケーション
1100 MRI装置

Claims (18)

  1. 装置の非一時的メモリにプログラムされる方法であって、
    アフィンレジストレーションを行って、被変形画像及び目標画像を、おおよそ同一の座標空間に位置決めするステップであって、前記被変形画像及び前記目標画像は、各々、脳構造を含む、ステップと、
    腫瘍適応指標を用いた微分同相デモンズ非剛体レジストレーションを実行して、前記脳構造の境界を精細化するステップと、
    を含み、
    前記微分同相デモンズ非剛体レジストレーションは、オプティカルフローから導出された変形力及び制御可能な正則化を利用するものであり、前記制御可能な正則化は、滑らかさ及びデータの一致のうちの一方を選択することによって前記滑らかさ又はデータの一致に基づいて優先順位付けされ、前記オプティカルフローから導出された変形力は、空間勾配に沿って不一致の強度に比例する速度を用いて計算される
    ことを特徴とする方法。
  2. 前記アフィンレジストレーションは、点、直線及び平面を保存することによって、画像を位置合わせすることを含むことを特徴とする、請求項1に記載の方法。
  3. 前記微分同相デモンズ非剛体レジストレーションは、ベクトル場及びベクトルの大きさを利用することを特徴とする、請求項1に記載の方法。
  4. 前記微分同相デモンズ非剛体レジストレーションは、変形場において、画像強度のマッチング誤差を最小化することを含むことを特徴とする、請求項1に記載の方法。
  5. 前記腫瘍適応指標は、腫瘍領域の誤差の寄与を除外することを特徴とする、請求項1に記載の方法。
  6. 前記微分同相デモンズ非剛体レジストレーションは、レジストレーション時間中、滑らかな変形軌跡を用いることを含むことを特徴とする、請求項1に記載の方法。
  7. アプリケーションを記憶するための非一時的メモリであって、
    アフィンレジストレーションを行って、被変形画像及び目標画像を、おおよそ同一の座標空間に位置決めするステップであって、前記被変形画像及び前記目標画像は、各々、脳構造を含む、ステップと、
    腫瘍適応指標を用いた微分同相デモンズ非剛体レジストレーションを実行して、前記脳構造の境界を精細化するステップと、を行うためのものである、
    非一時的メモリと、
    前記メモリに結合され、前記アプリケーションを処理するように構成される処理構成要素と、
    を備え、
    前記微分同相デモンズ非剛体レジストレーションは、オプティカルフローから導出された変形力及び制御可能な正則化を利用するものであり、前記制御可能な正則化は、滑らかさ及びデータの一致のうちの一方を選択することによって前記滑らかさ又はデータの一致に基づいて優先順位付けされ、前記オプティカルフローから導出された変形力は、空間勾配に沿って不一致の強度に比例する速度を用いて計算される
    ことを特徴とする装置。
  8. 前記アフィンレジストレーションは、点、直線及び平面を保存することによって、画像を位置合わせすることを含むことを特徴とする、請求項7に記載の装置。
  9. 前記微分同相デモンズ非剛体レジストレーションは、ベクトル場及びベクトルの大きさを利用することを特徴とする、請求項7に記載の装置。
  10. 前記微分同相デモンズ非剛体レジストレーションは、変形場において、画像強度のマッチング誤差を最小化することを含むことを特徴とする、請求項7に記載の装置。
  11. 前記腫瘍適応指標は、腫瘍領域の誤差の寄与を除外することを特徴とする、請求項7に記載の装置。
  12. 前記微分同相デモンズ非剛体レジストレーションは、レジストレーション時間中、滑らかな変形軌跡を用いることを含むことを特徴とする、請求項7に記載の装置。
  13. 磁気共鳴撮像装置と、
    計算装置であって、
    アフィンレジストレーションを行って、被変形画像及び目標画像を、おおよそ同一の座標空間に位置決めするステップであって、前記被変形画像及び前記目標画像は、各々、脳構造を含む、ステップと、
    腫瘍適応指標を用いた微分同相デモンズ非剛体レジストレーションを実行して、前記脳構造の境界を精細化するステップと、を行うように構成される、
    計算装置と、
    前記メモリに結合され、前記アプリケーションを処理するように構成される処理構成要素と、
    を備え、
    前記微分同相デモンズ非剛体レジストレーションは、オプティカルフローから導出された変形力及び制御可能な正則化を利用するものであり、前記制御可能な正則化は、滑らかさ及びデータの一致のうちの一方を選択することによって前記滑らかさ又はデータの一致に基づいて優先順位付けされ、前記オプティカルフローから導出された変形力は、空間勾配に沿って不一致の強度に比例する速度を用いて計算される
    ことを特徴とするシステム。
  14. 前記アフィンレジストレーションは、点、直線及び平面を保存することによって、画像を位置合わせすることを含むことを特徴とする、請求項13に記載のシステム。
  15. 前記微分同相デモンズ非剛体レジストレーションは、ベクトル場及びベクトルの大きさを利用することを特徴とする、請求項13に記載のシステム。
  16. 前記微分同相デモンズ非剛体レジストレーションは、変形場において、画像強度のマッチング誤差を最小化することを含むことを特徴とする、請求項13に記載のシステム。
  17. 前記腫瘍適応指標は、腫瘍領域の誤差の寄与を除外することを特徴とする、請求項13に記載のシステム。
  18. 前記微分同相デモンズ非剛体レジストレーションは、レジストレーション時間中、滑らかな変形軌跡を用いることを含むことを特徴とする、請求項13に記載のシステム。
JP2019547070A 2016-11-16 2017-08-15 患者間脳レジストレーション Active JP6781417B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US15/353,679 2016-11-16
US15/353,679 US10176570B2 (en) 2016-11-16 2016-11-16 Inter-patient brain registration
PCT/US2017/046974 WO2018093430A1 (en) 2016-11-16 2017-08-15 Inter-patient brain registration

Publications (2)

Publication Number Publication Date
JP2019535472A JP2019535472A (ja) 2019-12-12
JP6781417B2 true JP6781417B2 (ja) 2020-11-04

Family

ID=59791128

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019547070A Active JP6781417B2 (ja) 2016-11-16 2017-08-15 患者間脳レジストレーション

Country Status (6)

Country Link
US (1) US10176570B2 (ja)
EP (1) EP3526769A1 (ja)
JP (1) JP6781417B2 (ja)
KR (1) KR102227713B1 (ja)
CN (1) CN109952597B (ja)
WO (1) WO2018093430A1 (ja)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102019107952B4 (de) * 2019-03-27 2023-08-10 Volume Graphics Gmbh Computer-implementiertes Verfahren zur Analyse von Messdaten eines Objekts
CN112419377B (zh) * 2020-11-20 2021-10-08 推想医疗科技股份有限公司 配准图像确定方法及配准图像确定装置
CN115005802B (zh) * 2022-07-21 2023-04-28 首都医科大学宣武医院 脑网络疾病发作起始部位定位方法、系统及装置

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001022352A1 (en) * 1999-09-22 2001-03-29 Miller Michael I Method and apparatus for image registration using large deformation diffeomorphisms on a sphere
US8571278B2 (en) 2005-06-24 2013-10-29 The University Of Iowa Research Foundation System and methods for multi-object multi-surface segmentation
US7876938B2 (en) 2005-10-06 2011-01-25 Siemens Medical Solutions Usa, Inc. System and method for whole body landmark detection, segmentation and change quantification in digital images
US20070116381A1 (en) * 2005-10-19 2007-05-24 Ali Khamene Method for deformable registration of images
EP2081494B1 (en) 2006-11-16 2018-07-11 Vanderbilt University System and method of compensating for organ deformation
EP2369551B1 (en) 2010-03-25 2019-10-30 Emory University Imaging system and method
US9129382B2 (en) 2010-06-25 2015-09-08 Institute Of Automation, Chinese Academy Of Sciences Method and system for brain tumor segmentation in multi-parameter 3D MR images via robust statistic information propagation
CN103098090B (zh) 2011-12-21 2015-01-21 中国科学院自动化研究所 多参数三维磁共振图像脑肿瘤分割方法
CN102722890B (zh) * 2012-06-07 2014-09-10 内蒙古科技大学 基于光流场模型的非刚性心脏图像分级配准方法
CN103854276B (zh) * 2012-12-04 2018-02-09 东芝医疗系统株式会社 图像配准及分割装置和方法,以及医学图像设备
US9747421B2 (en) 2014-02-10 2017-08-29 Picofemto LLC Multi-factor brain analysis via medical imaging decision support systems and methods
CN104091337B (zh) * 2014-07-11 2017-07-14 北京工业大学 一种基于PCA及微分同胚Demons的变形医学图像配准方法

Also Published As

Publication number Publication date
KR102227713B1 (ko) 2021-03-15
WO2018093430A1 (en) 2018-05-24
CN109952597B (zh) 2023-03-31
EP3526769A1 (en) 2019-08-21
JP2019535472A (ja) 2019-12-12
CN109952597A (zh) 2019-06-28
US10176570B2 (en) 2019-01-08
US20180137617A1 (en) 2018-05-17
KR20190057396A (ko) 2019-05-28

Similar Documents

Publication Publication Date Title
US10846853B2 (en) Medical image processing apparatus, medical image processing method, and medical image processing program
US10614599B2 (en) Anatomy segmentation through low-resolution multi-atlas label fusion and corrective learning
JP7350878B2 (ja) 画像解析方法、装置、プログラム
JP6819898B2 (ja) 自動3d脳腫瘍セグメンテーション及び分類
US10424067B2 (en) Image processing apparatus, image processing method and storage medium
JP6781417B2 (ja) 患者間脳レジストレーション
JP2009520558A (ja) ポイント・ベースの適応的弾性画像登録
JP2002269559A (ja) 画像のテンプレートマッチング方法及び画像処理装置
US20180253861A1 (en) Information processing apparatus, method and non-transitory computer-readable storage medium
JP2005321902A (ja) 画像処理装置および方法、記録媒体、並びにプログラム
CN104508703B (zh) 用于配准图像序列的系统和方法
EP2624009B1 (en) Method and apparatus for obtaining symmetry information of objects
US20160275357A1 (en) Method and system for tracking a region in a video image
JP6824845B2 (ja) 画像処理システム、装置、方法およびプログラム
JP2014174174A (ja) 医用画像処理装置、医用画像診断装置及び画像処理方法
CN107274450A (zh) 信息处理装置及其控制方法
JP2013037539A (ja) 画像特徴量抽出装置およびそのプログラム
US11847786B2 (en) Motion learning without labels
CN110197100B (zh) 物体定位方法和装置
US9792533B2 (en) Processing apparatus, processing method, and program for determining an image taking position
US20210248769A1 (en) Array-based depth estimation
Yanli et al. A robust and accurate non-rigid medical image registration algorithm based on multi-level deformable model
Chen et al. Automatic fiducial localization in brain images
Gadermayr et al. Feature extraction with intrinsic distortion correction in celiac disease imagery: no need for rasterization
KR101627974B1 (ko) 영상 블러에 강인한 영상 특징 기술자 생성 방법 및 장치

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20190515

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200415

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200612

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20200916

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200929

R151 Written notification of patent or utility model registration

Ref document number: 6781417

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151