JP5231210B2 - 多次元データセットにおいて表面をセグメント化する方法、システム及びコンピュータプログラム - Google Patents

多次元データセットにおいて表面をセグメント化する方法、システム及びコンピュータプログラム Download PDF

Info

Publication number
JP5231210B2
JP5231210B2 JP2008506040A JP2008506040A JP5231210B2 JP 5231210 B2 JP5231210 B2 JP 5231210B2 JP 2008506040 A JP2008506040 A JP 2008506040A JP 2008506040 A JP2008506040 A JP 2008506040A JP 5231210 B2 JP5231210 B2 JP 5231210B2
Authority
JP
Japan
Prior art keywords
features
selection factor
segmentation
data set
matching
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.)
Expired - Fee Related
Application number
JP2008506040A
Other languages
English (en)
Other versions
JP2008535614A5 (ja
JP2008535614A (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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
Koninklijke Philips Electronics NV
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 Koninklijke Philips NV, Koninklijke Philips Electronics NV filed Critical Koninklijke Philips NV
Publication of JP2008535614A publication Critical patent/JP2008535614A/ja
Publication of JP2008535614A5 publication Critical patent/JP2008535614A5/ja
Application granted granted Critical
Publication of JP5231210B2 publication Critical patent/JP5231210B2/ja
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/149Segmentation; Edge detection involving deformable models, e.g. active contour models
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • 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/10072Tomographic images
    • 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/30048Heart; Cardiac

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Analysis (AREA)
  • Image Processing (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Description

本発明は、複数の画像を含む多次元データセットにおいて表面をセグメント化する方法に関する。
本発明は、更に、複数の画像を含む多次元データセットにおいて表面をセグメント化するシステムに関する。
本発明は、更に、ワークステーションに関する。
本発明は、更に、ビューイングステーションに関する。
本発明は、更に、複数の画像を含む多次元データセットにおいて表面をセグメント化するためのコンピュータプログラムに関する。
本発明は、更に、複数の画像を含む多次元データセットにおいて表面をセグメント化するためのユーザインタフェースに関する。
冒頭に記載した方法の実施例は、J.Weese他による刊行物「Shape Constrained Deformable Models for 3D Medical Image Segmentation」(Proc. IPMI 380-387, 2001)から知られている。既知の方法は、特に、セグメント化されることが企図されるオブジェクトの以前の形状情報を使用して、画像セグメンテーション方法のロバストネスを改善するように構成されている。既知の方法において、形状情報は、弾力的に変形可能な表面形状モデルに組み込まれており、画像への適応は、局所的な表面検出から得られる外部エネルギー及び内部エネルギーによって支配され、これは、変形可能な表面が、形状モデルによって規定されるサブスペースの近傍にとどまるように制約する。
既知の方法の不利益は、その信頼性が、動きうる身体の動きを記述するように構成される弾力性のある制約の精度に大きく依存することである。実際に、イメージングアーチファクト又は不完全な特徴抽出のため、画像から抽出される特徴のうちのいくつかは、事実、オブジェクトに属さず、又はオブジェクトの形状が、形状モデルによって与えられる変形レンジから局所的な逸脱を有することがある。従って、このような「アウトライア」の抽出される特徴が存在する場合、セグメンテーションの結果は不正確になることがある。
本発明の目的は、画像セグメンテーション方法の精度を改善することである。
このために、本発明による方法は、
−前記画像内の複数の特徴を入手するステップと、
−アプリオリに決定された選択ファクタを入手するステップと、
−多次元データセット内の前記複数の特徴をマッチングすることによって、表面を与えるように複数の画像をセグメント化するステップであって、各特徴について、個々のマッチングエラーが計算され、特徴のうち、マッチングエラーの大きさが上位である可変の割合の特徴が、選択ファクタによって与えられる限界まで無視される、ステップと、
を含む。
本発明の技術的な方策は、例えばマッチングエラーの大きさが上位である特徴の所定割合のような、選択ファクタを規定することによって、及び適切なセグメンテーションアルゴリズムにそれを組み込むことによって、セグメント化された表面の精度が改善されるという洞察に基づく。好適には、選択されることが企図される特徴は、セグメント化される表面の位置についての候補ロケーションに似ており、言い換えると、表面は、それが選択ファクタに従って除外された特徴を避け、すべての選択された特徴を通るように変形されるべきである。
予め設定された割合のアウトライアを無視することに基づく画像の位置合わせの概念は、それ自体、D.M. Mount他による「Efficient algorithms for robust feature matching」(Pattern Recogn. 32, p.17, 1999)から知られていることに留意されたい。知られている方法において、2つの時間的に不連続の画像が、計算技法を使用してマッチングされる。予め決められたマッチング精度が、画像内のいくつかの特徴を無視することによって保たれ、これは、画像内の他の特徴に関して顕著なマッチング精度を示す。
本発明による方法において、最適な割合のアウトライア特徴を識別する準備ステップにおいて、セグメント化された表面の可変割合の潜在的なロケーションを無視し、それによって正確なセグメンテーションの結果を与えることから利益が得られることを認識されたい。無視されるべき可変の割合の特徴に対応するあらゆる画像セグメンテーションは、好適にはスコアを割り当てられる。セグメンテーションは、個別のスコアに関して比較され、最善のものが選択される。選択ファクタの予め決められた値は、経験的に決定されることができる。
特に心臓のような動きうる身体の画像について使用されることが意図される適用例の場合、選択ファクタを数パーセントに設定することが有用であり、例えば10−15%の設定値が、信頼できる結果を提供するために示されている。選択ファクタの値があまりに低い場合、いくつかのアウトライア特徴がモデルマッチングの際に考慮され、結果として得られるセグメント化された表面の局所解剖学的な精度を低減させることに留意されたい。選択ファクタの値があまりに高い場合、いくつかの正しい特徴が、誤って除外されることがあり、これは、実際のオブジェクトの一部が正しくセグメント化されないという結果をもたらす。
本発明による方法の実施例において、方法は更に、
−多次元データセットにおいて表面を再構成するステップと、
−表面を表示手段に表示するステップと、
を含み、又は、
−多次元データセットにおいて表面を再構成するステップと、
−表面を表示手段に表示するステップと、
−他の選択ファクタを与えるように、選択ファクタを変更するステップと、
−複数の画像をセグメント化しなおすステップであって、画像特徴のうち、個々のマッチングエラーの大きさが上位である或る割合の画像特徴が、他の選択ファクタに従って無視される、ステップと、
−多次元データセットにおいて他の表面を再構成するステップと、
−他の表面を表示手段に表示するステップと、
を含む。
各々の画像セグメンテーションステップについて、選択ファクタの適応的な最適化を可能にすることが特に有利であることが分かった。方法は、半自動的なアプローチを使用することができ、オペレータは、セグメント化された表面と画像内のオブジェクトとの間の局所解剖学的な適合の程度を評価することが可能である。特に選択ファクタの初期値が深遠な経験的知識に基づく最適な状況において、更なる繰り返しは必要でない。このように、オペレータは、好適には元のデータに重ねられる表示される表面を調べ、セグメンテーションの結果を受け入れる。セグメンテーションの結果は、例えばアーカイビングの目的で保存され、及び/又はリモート解析のためにエクスポートされることもできる。オペレータが、より良好な局所解剖学的な適合が便利であると思う場合、オペレータは、選択ファクタを適切に変更することができ(増加させ又は減少させる)、本発明による方法は、他の選択ファクタに基づいて他の表面をセグメント化しなおし、好適には元のデータに再構成された表面を重ねることによって再構成された他の表面を適切な表示手段上に視覚化することによって、進む。オペレータが満足する場合、オペレータは、セグメンテーションルーチンを出る。そうでない場合、オペレータは、選択ファクタを変更し、プロセスは、オペレータが満足するまで繰り返す。我々の経験から、セグメンテーションは、ほんの数回の繰り返しのうちに収束することが分かっている。本発明による方法の本実施例によれば、高速且つロバストな自動化されたセグメンテーションが、最小限のオペレータ対話によって可能にされ、それによって、画像セグメンテーション方法は、未知の割合の画像アーチファクト及びモデルからのオブジェクトの形状ずれにも適応できる。
本発明による方法の他の実施例において、選択ファクタを変更するステップは、予め決められたアルゴリズムに従って自動的に実施される。
選択ファクタの反復的な適応化のプロセスを自動化することが好ましいことが分かる。好適には、適切なルックアップテーブルが、選択ファクタを変更するために前もって準備される。代替例として、予め選択された関数に従って選択ファクタを変更することも可能である。更なる代替例として、表面ベースの方策の代わりに又はそれに加えて、セグメンテーション品質の領域ベースの方策を使用することも可能である。セグメンテーション品質の適切な領域ベースの方策の例は、例えばセグメンテーションを表す表面によって描かれる領域(例えば、表面の内側及び外側)のような、予め規定された領域内の画像のテクスチャの均一性の計算である。更なる代替例として、ユーザは、ユーザの満足の程度を示すことによって、変更プロセスを速めるためのオプションを与えられることができる。満足の程度が低い場合、アルゴリズムは、ルックアップテーブルの他のエントリに又は関数従属性の他のファクタにジャンプすることができる。変更アルゴリズムのさまざまな変形が可能である。例えばモデルベースの特に医用画像セグメンテーションの適用例の場合、モデルは、一般に、例えば専門家によって実行される注意深い手動のトレーシングプロシージャを使用して得られる基準セグメンテーションの組から構築される。初期の選択ファクタP_start及び選択ファクタを変更するためのインクリメントP_stepの2つのパラメータが、モデル構築プロセスの副産物として評価されることができる。基準セグメンテーションが利用できるので、自動化された最適化プロシージャが、各々の画像データセットについてアウトライアの最適な割合を与える最適な最終の選択ファクタP=P_finalを見つけるために使用されることができる。セグメンテーションプロセスの最適化のために、適切な複数の候補セグメンテーションが、さまざまな選択ファクタに対応して予め計算される。ユーザは、対応するセグメント化された表面間をスクロールして、最適の表面を選択することができる。更に、ユーザは、すでに計算されたセグメンテーションから好適なセグメンテーションを選択することができるとともに、選択ファクタを適切に対話的に更に調整するための始点として、それを使用することができる。ある経験的実践に基づいて、P_finalについてのヒストグラムが、選択ファクタP_finalの或る値ごとに最終選択の数を示して構築されることができる。このヒストグラムに基づいて、選択ファクタの予め記憶された値であるP_startが、新しい検査されていないデータセットについて予め設定される。インクリメントP_stepは、或る画像カテゴリに関するヒストグラムの幅に依存して選択されることができる。このように、ロバストな自動化された画像セグメンテーション方法が構築され、その結果、画像セグメンテーションステップの局所解剖学的な精度が改善される。好適には、適切な画像セグメンテーションルーチンのために、変形可能なモデルに基づく画像セグメンテーションアルゴリズムが選択される。
本発明によるシステムは、
−前記画像内の複数の特徴を入手し、アプリオリに決定された選択ファクタを入手する入力部と、
−多次元データセット内の前記複数の特徴をマッチングすることによって、表面を与えるように前記複数の画像をセグメント化する計算手段であって、各特徴について、個々のマッチングエラーが計算され、特徴のうち、個々のマッチングエラーの大きさが上位である可変の割合の特徴が、選択ファクタによって与えられる限界まで無視される、計算手段と、
を有する。
本発明によるシステムは、オフラインのデータ解析を可能にするワークステーションとして実現されることが有利である。好適には、本発明によるシステムは、多次元データセットにおいて表面を再構成する再構成ユニットと、表面を表示する表示手段と、を更に有する。更に好適には、本発明によるシステムが、オフラインの対話的なデータ解析を可能にするビューイングステーションとして実現されることが有利である。更に好適には、本発明によるシステムは、多次元データセットを取得するデータ取得ユニットを有する。
本発明によるコンピュータプログラムは、
−前記画像内の複数の特徴を入手するステップと、
−アプリオリに決定された選択ファクタを入手するステップと、
−多次元データセット内の前記複数の特徴をマッチングすることによって、表面を与えるように前記複数の画像をセグメント化するステップであって、各々の特徴について、個々のマッチングエラーが計算され、特徴のうち、個々のマッチングエラーの大きさが上位である可変の割合の特徴が、選択ファクタによって与えられる限界まで無視される、ステップと、
をプロセッサに実行させるための命令を含む。
本発明によるユーザインタフェースは、
−再構成された表面を表示手段に表示する処理であって、前記再構成された表面が、複数の特徴を含む多次元データセット内の個々の画像部分のセグメンテーションに基づいて計算され、特徴のうち、マッチングエラーの大きさが上位である可変の割合の特徴が、予め決められた選択ファクタによって設定される限界まで無視される、処理と、
−他の選択ファクタを与えるように、選択ファクタを変更する処理と、
−表面を与える複数の画像をセグメント化しなおす処理であって、特徴のうち、個々のマッチングエラーの大きさが上位である或る割合の特徴が、他の選択ファクタに従って無視される、処理と、
−多次元データセットにおいて他の表面を再構成する処理と、
−他の表面を表示手段に表示する処理と、
を行うように構成される。
好適には、ユーザインタフェースは、更に、各々の繰り返しごとに各々の個々の他の選択ファクタを記憶するように構成される。この場合、選択ファクタがその最適値を越えて変更され、それゆえ局所解剖学的な適合が、以前の繰り返しよりも低下する状況において、適切なバックトラッキングが可能である。この場合、ユーザは、選択ファクタの以前の値へ戻ることが可能にされ、その値から始まって、ユーザは、選択ファクタのインクリメントP−stepの値を変更することができる。このようにして、選択ファクタの微調整された最適化が可能にされる。好適には、ユーザインタフェースは、更に、選択ファクタのさまざまな異なる値についてのセグメンテーションの事前の、特にバックグランド又はオフラインの複数の計算に応じて、複数の表面セグメンテーションを表示するように構成される。こうして、ユーザは、最も良好なものを選択するためにセグメント化しなおすことを待つ必要がない。この技術的な方策は、ワークフローを改善する。
本発明のこれら及び他の見地は、図面を参照してより詳細に議論される。
図1は、本発明による方法の実施例の概略図を示している。本発明による方法1は、複数の画像を含む多次元データセットにおいて表面をセグメント化するように構成されている。好適には、画像は、準備ステップ2において、適切なデータ取得ユニットを使用して取得される。データ処理ステップ及びデータ取得ステップは、時間的に又は地理的に離されることが可能であり、それによりステップ4において、適切なデータセグメンテーションステップ6の結果が入手されることが可能であり、前記結果は、本発明の方法を使用して表面をセグメント化するためにのちに使用される画像の部分を含む。画像セグメンテーションのステップは、リアルタイムに実施されることができ、又は前もって実行されてもよいことに留意されたい。後者の場合、初期データセグメンテーション6の結果は、適切なコンピュータ手段によって入手される。次に、ステップ8において、セグメント化されようとする表面の予想される空間的形状と類似する適切な複数の画像特徴が選択され、入手される。このような選択は、事前の知識に基づいて、又は識別可能な特徴の全体量のパーセンテージのような或る基準に従って、実施されることができる。前記パーセンテージが100%にセットされてもよいことに留意されたい。適切な量の特徴を選択する際、それらの特徴は、ステップ10において、すべての画像部分についてマッチングされる。特徴マッチングのプロシージャは、それ自体、データ処理の分野において知られており、ここでは詳しく説明されない。しかしながら、ステップ10のマッチングプロシージャの間、前記複数の特徴の各特徴について、マッチングエラーが割り当てられることが極めて重要である。ステップ14において、予め規定された選択ファクタが入手される。複数の選択ファクタ14aが予め記憶されることが可能であり、それによって、実際の選択ファクタが、前記予め記憶された複数の選択ファクタから選択されることが可能である。ステップ14において入手される選択ファクタに従って、特徴のうち、マッチングエラーの大きさが上位である或る割合の特徴が除外され、求められる表面のセグメンテーションが、ステップ16において実施され、それによって、除外された特徴は、表面の候補変形の適合の品質を評価するための考慮に入れられない。結果として得られる表面は、ステップ18において適切な表示手段に表示され、好適には、ユーザによる解析のために、表面は、元の、特に診断データに重ね合わせられる。ユーザがセグメンテーションの結果に満足する場合、セグメント化された表面は、ステップ19において記憶される。ユーザがセグメンテーションの結果に満足しない場合、プロセスは、ステップ14に戻され、他の選択ファクタが入手される。他の選択ファクタは、ユーザによる入力によってステップ11において変更されることができ、又は代替として、或る変更アルゴリズム13によって供給されることもできる。簡単な変更アルゴリズムの例は、ルックアップテーブルの複数の予め記憶された値の間のループである。代替例として、関数従属性が使用されることができ、又は画像均一性ファクタに基づく洗練されたアルゴリズムが使用されることができる。前記画像均一性ファクタは、予め決められた領域内の画像のテクスチャの均一性に基づく。他の選択ファクタは、リアルタイムに計算され、画像均一性ファクタに基づいて(自動的に)最適化される。言い換えると、最も正確なセグメンテーションのための最適な選択ファクタが、ユーザによって指示される代わりに自動的に計算される。
図2は、本発明によるシステムの実施例の概略図を示している。本発明によるシステムは、複数の画像を含む多次元データセットにおいて表面をセグメント化するように構成される。システム30は、前記画像内の複数の特徴を入手するための入力部32を有する。特徴は、適切な準備の画像セグメンテーションステップの結果であってもよいことに留意されたい。特徴は、局所的に又は遠隔的に位置しうる適切な記憶ユニット(図示せず)から入手されてもよい。代替として及び/又は付加的に、入力部32は、適切な内部データセグメンテーションアルゴリズム35dから前記特徴を受け取るように構成されることができる。この場合、多次元データセットが、入力部32によって入手され、セグメンテーションアルゴリズム35dを使用して準備の画像セグメンテーションを実施するために、システム30の計算ユニット35に対して利用可能にされることに留意されたい。
システム30のコア部は、入力部32、計算ユニット35、作業メモリ36及びバックグラウンド記憶ユニット38であるシステム30の構成要素を動作させるように構成されるプロセッサ34によって形成される。適切なプロセッサ34の例は、通常のマイクロプロセッサ又は信号プロセッサ、バックグラウンド記憶装置38(一般にはハードディスクに基づく)及び作業メモリ36(一般にはRAMに基づく)である。バックグラウンド記憶装置38は、処理されていないときには適切なデータセット(又はその一部)を記憶するために使用されることができ、画像セグメンテーションステップ、適切な特徴及び選択ファクタを選択するステップ及び任意の他の適切な中間又は最終の計算ステップ、の結果を記憶するために使用されることができる。作業メモリ36は、一般に、処理されているデータセット及び表面の一部のセグメンテーションの結果を保有する。計算ユニット35は、好適には、適切な数の実行可能なサブルーチン35a、35b、35c及び35dを含む。サブルーチン35aは、画像の一部分の中の適切な複数の特徴の選択を実施するように構成される。サブルーチン35bは、多次元データセット内の選択された複数の特徴をマッチングするように構成され、それによって、各々の特徴について、マッチングエラーが割り当てられる。サブルーチン35cは、選択ファクタPの実際の値を入手するように及び/又はそれを計算するように構成される。サブルーチン35dは、特徴のうち、マッチングエラーの大きさが上位である可変の割合の特徴を除外することによって表面をセグメント化するように構成され、この場合、このような特徴の最大量は、実際の選択ファクタPによって支配される。結果として得られるサブセグメンテーションは、スコアを割り当てられ、それによってサブセグメンテーションは、個別のスコアに関して相互に比較される。最適なサブセグメンテーションが、求められる表面を示す結果として得られたセグメンテーションであるように選択される。
本発明によるシステム30は、セグメント化された表面と元のデータ、特に診断画像との適切な重ね合わせのレンダリングを生成するように構成される重ね合わせコーダ37を更に有する。好適には、計算された重ね合わせは、ファイル37aに記憶される。好適には、重ね合わせコーダ37、計算ユニット35及びプロセッサ34は、好適にはメモリ38に記憶されるコンピュータプログラム33によって動作可能である。出力部39は、セグメント化された表面の適切なレンダリングを重ね合わせた心臓の解剖学的構造を表す重ね合わせられた画像データのような、処理の結果を出力するように使用される。
図3は、本発明による装置の他の実施例の概略図を示している。システム40は、例えば時間的に間隔のあいた複数の心臓画像を含む多次元データセットにおいて表面をセグメント化するように構成される。好適には、システム40は、特に多次元データセットを取得するための磁気共鳴撮影装置、トモグラフィユニット、超音波装置又はX線ユニットのようなデータ取得装置41を有する。通常、データは、適切にコード化された信号Sによって、データ取得装置41からプロセッサ42に転送されることが企図される。図2を参照して説明したように、プロセッサは、適切なデータセグメンテーションを実施し、その出力部において、さまざまな可能なデータが生成されうる。例えば、データ42aは、第1の選択ファクタに基づく第1の潜在的な表面のセグメンテーションを含み、データ42bは、第2の選択ファクタに基づく第2の潜在的な表面のセグメンテーションを提供する、等が可能である。好適には、プロセッサ42は、ワークステーション44aに組み込まれる。
データ42a、42bのいずれか又はそれらの適切な組み合わせが、適切なビューワ43の他の入力部45に利用可能にされる。好適には、他の入力部45は、ユーザインタフェース48を制御するように適応されるプログラム46を使用して、適切なインタフェースを動作させるように構成される適切な他のプロセッサを有し、それにより解剖学的データの画像が、セグメンテーションステップの結果、特にデータ42a、42bと適切に重ね合わせられ、結果的に画像部分48a、48bを与える。好適には、ユーザの便宜のために、ビューワ43は、高解像度表示手段47を備えており、ユーザインタフェースは、例えばマウス、キーボード又は他の任意の適切なユーザの入力装置のような、適切なインタラクティブ手段49によって動作可能である。好適には、ユーザインタフェースは、実際の選択ファクタPを変更するために、ユーザが画像と対話することを可能にし、実際の選択ファクタPは、表面の他のセグメンテーションを計算するためにシステム40によって使用される。適切なグラフィックユーザ入力が、コンピュータプログラム46によって変数に変えられる。このオプションは、予め記憶された又は計算された選択ファクタのレンジが満足のいくセグメンテーション結果を供給しないとき、表面の正確なセグメンテーションを可能にする。好適には、プロセッサ42及びビューワ43は、ビューイングステーション44bを形成するように構成される。
図4は、本発明によるコンピュータプログラムのフローチャートの実施例を図式的に示している。本発明によるコンピュータプログラム50は、複数の画像を含む多次元データセットにおいて表面をセグメント化するように構成される。好適には、画像は、準備ステップ52において、適切なデータ取得ユニットを使用して取得される。データ処理ステップ及びデータ取得ステップは、時間的に又は地理的に離されており、ステップ54において、適切なデータセグメンテーションステップ56の結果が、適切なデータアクセスサブルーチンを使用して入手され、それによって個々の画像特徴を与えることが可能である。これらの特徴は、本発明の方法を使用して、表面をセグメント化するためにのちに使用される画像の部分を含むことが企図される。画像セグメンテーション56のステップは、リアルタイムに実施されることができ、又は前もって実行されることができることに留意されたい。前者の場合、適切なそれ自体知られているデータセグメンテーションアルゴリズムが使用されうる。後者の場合、初期のデータセグメンテーション56の結果が、適切なコンピュータ入力/出力データアクセスプロトコルによって入手される。次に、ステップ58において、セグメント化されようとする表面の予想される空間的形状と類似する適切な複数の画像特徴が、選択される。このステップは、或る特徴を画像から抽出するように構成される適切な画像解析アルゴリズムによって実現されることができる。特徴抽出の方法は、それ自体、当分野において知られており、ここでは説明されない。更に、特徴選択は、事前の知識に基づいて又は識別可能な特徴の全体量のパーセンテージのような或る基準に従って、実施されることもできる。前記パーセンテージは、100%に設定されてもよいことに留意されたい。
適切な量の特徴を選択する際、それらの特徴は、適切なデータマッチングサブルーチンを使用して、ステップ60において、すべての画像部分についてマッチングされる。特徴マッチングのプロシージャは、それ自体、データ処理の分野において知られており、ここで詳しく説明されない。しかしながら、ステップ60のマッチングプロシージャの間に、前記複数の特徴の各特徴について、マッチングエラーが割り当てられることが極めて重要である。ステップ64において、予め規定された選択ファクタは、適切なデータアクセスサブルーチンによって入手される。複数の選択ファクタ64aが、予め記憶され、それによって、実際の選択ファクタが、前記予め記憶された複数の選択ファクタから選択されることが可能である。ステップ64において入手される選択ファクタに従って、特徴のうち、マッチングエラーの大きさが上位である可変の割合の特徴が、選択ファクタによって与えられる最大限まで除外され、それによって複数の個々のサブセグメンテーションを与える。各々のサブセグメンテーションは、スコアを割り当てられ、前記スコアに関して最善のものが選択され、ステップ66において求められる表面が与えられる。除外された特徴は、スコアを計算するための考慮に入れられない。結果として得られた表面は、ステップ68において適切な表示手段に表示される。好適には、前記表面は、ユーザによる解析のために、元の特に診断データに重ね合わせられる。ステップ68の実現のために、本発明によるコンピュータプログラムは、適切なグラフィックユーザインタフェースを制御するように構成されることが好ましい。ユーザがセグメンテーションの結果に満足する場合、ユーザは、適当なコンピュータコマンドによってこのイベントを知らせ、そののち、セグメント化された表面がステップ69において記憶される。ユーザがセグメンテーションの結果に満足しない場合、プロセスは、ステップ64に戻され、他の選択ファクタが入手される。他の選択ファクタは、ステップ61においてユーザ入力によって変更されることが可能であり、又は代替として或る変更アルゴリズム63によって提供されることができる。簡単な変更アルゴリズムの例は、ルックアップテーブルの複数の予め記憶された値の間のループである。代替例として、関数従属性が使用されることができ、又は画像均一性ファクタに基づく洗練されたアルゴリズムが使用されてもよい。前記画像均一性ファクタは、予め決められた領域内の画像のテクスチャの均一性に基づく。他の選択ファクタは、リアルタイムに計算され、画像均一性ファクタに基づいて最適化される。
図5は、本発明によるユーザインタフェースの実施例を図式的に示している。本発明によるユーザインタフェース70の実施例は、英数字情報及び画像情報が投影されることが可能なグラフィックインタフェース74を含む。このために、グラフィックインタフェース74は、78a、78bのような作動可能なボタンを有するインタラクティブウィンドウ78を含む。オペレータが、「表示」ボタン78aを始動すると、適切な複数の画像74a、74bが、グラフィックスクリーン74の画像スクリーン内に投影されることができる。この例において、2つの画像が表示されており、これらの画像は、セグメント化された表面77b、79bが重ねられた診断情報77a、79aを含んでいる。本実施例において、セグメント化された表面は、異なる選択ファクタに関して除外された異なる割合の特徴に対応する。代替例として、各々の画像74a又は74bは、予め計算されたセグメンテーションを通じたムービールーピングを含むことができ、それによって、ユーザは、予め計算されたセグメンテーションのうちの1つを承認することができ、又は手動で選択ファクタを変更することができる。後者の場合、新しいセグメンテーションの結果が、それに応じてスクリーン72に投影される。
本発明による方法の実施例の概略図。 本発明によるシステムの実施例を図式的に示す図。 本発明によるシステムの他の実施例を図式的に示す図。 本発明によるコンピュータプログラムのフローチャートの実施例を図式的に示す図。 本発明によるユーザインタフェースの実施例を図式的に示す図。

Claims (13)

  1. 複数の画像を含む多次元データセットにおいて表面をセグメント化する方法であって、
    前記画像内の、セグメント化しようとする表面の予想される空間的形状と類似する複数の特徴を入手するステップと、
    アプリオリに決定された選択ファクタを入手するステップであって、前記選択ファクタは、前記複数の特徴のうち、前記セグメント化の実施において無視されるべき特徴の割合を規定するものである、ステップと、
    前記多次元データセット内の前記複数の特徴を形状モデルとマッチングして、各特徴について個々のマッチングエラーを計算し、前記マッチングエラーを考慮して当該形状モデルを変形させることにより、前記表面をセグメント化するステップであって、前記セグメント化は、前記複数の特徴のうち、マッチングエラーの大きい上位のものを前記選択ファクタによって規定される割合無視して当該形状モデルを変形させることにより実施され、前記選択ファクタの値は、前記セグメント化の異なる結果を与えるように変更可能である、ステップと、
    を含む方法。
  2. 前記多次元データセットにおいて前記表面を再構成するステップと、
    前記表面を表示手段に表示するステップと、
    を更に含み、又は
    前記多次元データセットにおいて前記表面を再構成するステップと、
    前記表面を表示手段に表示するステップと、
    他の選択ファクタを与えるように、前記選択ファクタを変更するステップと、
    前記表面をセグメント化しなおすステップであって、前記セグメント化が、前記複数の特徴のうち、マッチングエラーの大きい上位のものを前記他の選択ファクタによって規定される割合無視して当該形状モデルを変形させることにより実施される、ステップと、
    前記多次元データセットにおいて他の表面を再構成するステップと、
    前記他の表面を前記表示手段に表示するステップと、
    を更に含む、請求項1に記載の方法。
  3. 前記選択ファクタを変更する前記ステップが、予め決められたアルゴリズムに従って自動的に実施される、請求項2に記載の方法。
  4. 前記アルゴリズムが、予め決められた領域内の前記画像のテクスチャの均一性に基づいて前記選択ファクタを最適化する、請求項3に記載の方法。
  5. 前記アルゴリズムは、前記他の選択ファクタの最適化を可能にするように構成される、請求項3又は請求項4に記載の方法。
  6. 複数の画像を含む多次元データセットにおいて表面をセグメント化するシステムであって、
    前記画像内の、セグメント化しようとする表面の予想される空間的形状と類似する複数の特徴を入手するとともに、アプリオリに決定された選択ファクタを入手する入力部であって、前記選択ファクタは、前記複数の特徴のうち、前記セグメント化の実施において無視されるべき特徴の割合を規定するものである、入力部と、
    前記多次元データセット内の前記複数の特徴を形状モデルとマッチングして、各特徴について個々のマッチングエラーを計算し、前記マッチングエラーを考慮して当該形状モデルを変形させることにより、前記表面をセグメント化する計算手段であって、前記セグメント化は、前記複数の特徴のうち、マッチングエラーの大きい上位のものを前記選択ファクタによって規定される割合無視して当該形状モデルを変形させることにより実施され、前記選択ファクタの値は、前記セグメント化の異なる結果を与えるように変更可能である、計算手段と、
    を有するシステム。
  7. 前記多次元データセットにおいて前記表面を再構成する再構成ユニットと、
    前記表面を表示する表示手段と、
    を更に有する、請求項6に記載のシステム。
  8. 前記多次元データセットを取得するデータ取得ユニットを更に有する、請求項6又は請求項7に記載のシステム。
  9. 請求項7に記載のシステムを有するワークステーション。
  10. 請求項8に記載のシステムを有するビューイングステーション。
  11. 複数の画像を含む多次元データセットにおいて表面をセグメント化するためのコンピュータプログラムであって、
    前記画像内の、セグメント化しようとする表面の予想される空間的形状と類似する複数の特徴を入手するステップと、
    アプリオリに決定された選択ファクタを入手するステップであって、前記選択ファクタは、前記複数の特徴のうち、前記セグメント化の実施において無視されるべき前記特徴の割合を規定するものである、ステップと、
    前記多次元データセット内の前記複数の特徴を形状モデルとマッチングして、各特徴について個々のマッチングエラーを計算し、前記マッチングエラーを考慮して当該形状モデルを変形させることにより、前記表面をセグメント化するステップであって、前記セグメント化は、前記複数の特徴のうち、マッチングエラーの大きい上位のものを前記選択ファクタによって規定される割合無視して当該形状モデルを変形させることにより実施され、前記選択ファクタの値は、前記セグメント化の異なる結果を与えるように変更可能である、ステップと、
    をプロセッサに実行させるための命令を含む、コンピュータプログラム。
  12. 複数の画像を含む多次元データセットにおいてセグメント化された表面を表示手段に表示するためのユーザインタフェースを提供するコンピュータプログラムであって、
    再構成された表面を前記表示手段に表示する処理であって、前記画像内の、セグメント化しようとする表面の予想される空間的形状と類似する複数の特徴を入手し、アプリオリに決定された選択ファクタを入手し、前記複数の特徴を形状モデルとマッチングして各特徴について個々のマッチングエラーを計算し、前記マッチングエラーを考慮して当該形状モデルを変形させることにより、前記表面をセグメント化することに基づいて、前記再構成された表面を計算することを含み、前記選択ファクタは、前記複数の特徴のうち、前記セグメント化の実施において無視されるべき特徴の割合を規定し、前記セグメント化は、前記複数の特徴のうち、マッチングエラーの大きい上位のものを前記選択ファクタによって規定される割合無視して当該形状モデルを変形させることにより実施される、処理と、
    他の選択ファクタを与えるように、前記選択ファクタを変更する処理と、
    前記表面をセグメント化しなおす処理であって、前記複数の特徴のうち、マッチングエラーの大きい上位のものを前記他の選択ファクタによって規定される割合無視して当該形状モデルを変形させることにより実施される、処理と、
    前記多次元データセットにおいて前記他の表面を再構成する処理と、
    前記他の表面を前記表示手段に表示する処理と、
    をプロセッサに実施させるためのコンピュータプログラム。
  13. 前記コンピュータプログラムは更に、各々の繰り返しについて各々の個々の他の選択ファクタを記憶するように構成される、請求項12に記載のコンピュータプログラム。
JP2008506040A 2005-04-13 2006-04-13 多次元データセットにおいて表面をセグメント化する方法、システム及びコンピュータプログラム Expired - Fee Related JP5231210B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP05102894 2005-04-13
EP05102894.2 2005-04-13
PCT/IB2006/051145 WO2006109269A1 (en) 2005-04-13 2006-04-13 A method, a system and a computer program for segmenting a surface in a multi-dimensional dataset

Publications (3)

Publication Number Publication Date
JP2008535614A JP2008535614A (ja) 2008-09-04
JP2008535614A5 JP2008535614A5 (ja) 2013-03-07
JP5231210B2 true JP5231210B2 (ja) 2013-07-10

Family

ID=36822416

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008506040A Expired - Fee Related JP5231210B2 (ja) 2005-04-13 2006-04-13 多次元データセットにおいて表面をセグメント化する方法、システム及びコンピュータプログラム

Country Status (5)

Country Link
US (1) US8144987B2 (ja)
EP (1) EP1875435B1 (ja)
JP (1) JP5231210B2 (ja)
CN (1) CN101189639B (ja)
WO (1) WO2006109269A1 (ja)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8503761B2 (en) * 2009-11-12 2013-08-06 Harris Corporation Geospatial modeling system for classifying building and vegetation in a DSM and related methods
CN103648385B (zh) * 2011-04-28 2016-03-16 皇家飞利浦有限公司 具有用于选择候选分割图像的独立按钮的医学成像设备
WO2012176100A1 (en) * 2011-06-22 2012-12-27 Koninklijke Philips Electronics N.V. System and method for processing a medical image
US9629598B2 (en) * 2014-02-27 2017-04-25 Impac Medical Systems, Inc. System and method for auto-contouring in adaptive radiotherapy
US10055671B2 (en) * 2014-06-26 2018-08-21 Siemens Aktiengesellschaft Automatic assessment of perceptual visual quality of different image sets

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5999653A (en) * 1996-01-19 1999-12-07 Xerox Corporation Fast techniques for searching images using the Hausdorff distance
SE510664C2 (sv) * 1996-10-29 1999-06-14 Ericsson Telefon Ab L M Metoder och anordning för meddelandehantering i ett kommunikationssystem
DE19705017A1 (de) * 1997-02-10 1998-08-13 Delphi Systemsimulation Gmbh Verfahren zum Segmentieren eines digitalen Bildes
US6591004B1 (en) 1998-09-21 2003-07-08 Washington University Sure-fit: an automated method for modeling the shape of cerebral cortex and other complex structures using customized filters and transformations
US6975755B1 (en) * 1999-11-25 2005-12-13 Canon Kabushiki Kaisha Image processing method and apparatus
US6728424B1 (en) * 2000-09-15 2004-04-27 Koninklijke Philips Electronics, N.V. Imaging registration system and method using likelihood maximization
US6816607B2 (en) * 2001-05-16 2004-11-09 Siemens Corporate Research, Inc. System for modeling static and dynamic three dimensional anatomical structures by 3-D models
DE10144004A1 (de) 2001-09-07 2003-03-27 Philips Corp Intellectual Pty Verfahren zur Messung geometrischer Größen einer in einem Bild enthaltenen Struktur
US20030132936A1 (en) * 2001-11-21 2003-07-17 Kevin Kreeger Display of two-dimensional and three-dimensional views during virtual examination
CN100541541C (zh) * 2001-12-29 2009-09-16 田捷 基于单层表面跟踪的超大规模医学图像表面重建方法
WO2004036500A2 (en) * 2002-10-16 2004-04-29 Koninklijke Philips Electronics N.V. Hierarchical image segmentation
EP1437685A3 (de) * 2003-01-13 2004-07-21 Philips Intellectual Property & Standards GmbH Verfahren zum Segmentieren einer dreidimensionalen Struktur
WO2004068401A2 (en) 2003-01-30 2004-08-12 Philips Intellectual Property & Standards Gmbh Method for the interactive segmentation of an object in a three-dimensional data record
DE602004028471D1 (de) * 2003-06-13 2010-09-16 Philips Intellectual Property Dreidimensionale bildsegmentierung
US7103399B2 (en) * 2003-09-08 2006-09-05 Vanderbilt University Apparatus and methods of cortical surface registration and deformation tracking for patient-to-image alignment in relation to image-guided surgery
CN1265316C (zh) * 2004-06-11 2006-07-19 浙江大学 球形生物物料图像的灰度变换方法
US8295541B2 (en) * 2004-06-30 2012-10-23 Vision Fire & Security Pty Ltd System and method for detecting a change in an object scene

Also Published As

Publication number Publication date
CN101189639A (zh) 2008-05-28
CN101189639B (zh) 2010-09-29
US8144987B2 (en) 2012-03-27
EP1875435B1 (en) 2013-10-02
JP2008535614A (ja) 2008-09-04
WO2006109269A1 (en) 2006-10-19
US20090060319A1 (en) 2009-03-05
EP1875435A1 (en) 2008-01-09

Similar Documents

Publication Publication Date Title
Zhao et al. An overview of interactive medical image segmentation
JP5660432B2 (ja) 領域データ編集装置、領域データ編集方法、プログラム、及び、記録媒体
Hamarneh et al. 3D live-wire-based semi-automatic segmentation of medical images
JP6877534B2 (ja) 三次元医用データの対話型輪郭描出のための画像処理システム及び方法
EP2916738A1 (en) Lung, lobe, and fissure imaging systems and methods
WO2012017375A2 (en) In-plane and interactive surface mesh adaptation
US10275114B2 (en) Medical image display system and method for providing user interface enabling three-dimensional mesh to be edited on three-dimensional volume
JP2006506164A (ja) 3次元メッシュモデルをオブジェクトの3次元表面に自動的に適応させる画像処理システム
JP5231210B2 (ja) 多次元データセットにおいて表面をセグメント化する方法、システム及びコンピュータプログラム
US20090010505A1 (en) Method, a System and a Computer Program for Segmenting a Structure in a Dataset
US7668349B2 (en) Three-dimensional segmentation using deformable surfaces
US20060210158A1 (en) Object-specific segmentation
WO2005078666A1 (en) Real time user interaction with deformable surfaces of segmentation
JP2008535613A (ja) 多次元データセットにおける生体構造をセグメント化する方法、装置及びコンピュータプログラム
JP2008535614A5 (ja)
EP1851721B1 (en) A method, a system and a computer program for segmenting a surface in a multidimensional dataset
US8805073B2 (en) Image segmentation by assigning classes to adaptive mesh primitives
US20220375099A1 (en) Segmentating a medical image
Vanacken et al. Force feedback to assist active contour modelling for tracheal stenosis segmentation

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090410

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110517

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20110812

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20110819

A524 Written submission of copy of amendment under section 19 (pct)

Free format text: JAPANESE INTERMEDIATE CODE: A524

Effective date: 20110829

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20111020

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20120120

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20120127

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120419

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120719

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20121018

A602 Written permission of extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A602

Effective date: 20121025

A524 Written submission of copy of amendment under section 19 (pct)

Free format text: JAPANESE INTERMEDIATE CODE: A524

Effective date: 20130121

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: 20130219

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20130321

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20160329

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20160329

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees