JP2007523402A - 内部個別要素を用いるメッシュモデル - Google Patents

内部個別要素を用いるメッシュモデル Download PDF

Info

Publication number
JP2007523402A
JP2007523402A JP2006548466A JP2006548466A JP2007523402A JP 2007523402 A JP2007523402 A JP 2007523402A JP 2006548466 A JP2006548466 A JP 2006548466A JP 2006548466 A JP2006548466 A JP 2006548466A JP 2007523402 A JP2007523402 A JP 2007523402A
Authority
JP
Japan
Prior art keywords
internal
mesh
processing system
image processing
individual elements
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.)
Pending
Application number
JP2006548466A
Other languages
English (en)
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 JP2007523402A publication Critical patent/JP2007523402A/ja
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation
    • G06T17/205Re-meshing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20Finite element generation, e.g. wire-frame surface description, tesselation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Analysis (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Image Processing (AREA)
  • Image Generation (AREA)

Abstract

画像処理システムは、表面個別の要素(T)及び内部個別の要素(TH)を有する非構造化変形可能メッシュモデルを用いて、注目すべき対象物のセグメント化の画像データ処理手段を有し、表面個別の要素のサイズの局所的変化に従って、内部個別の要素のサイズを自動的且つ動的に適合させることにより非構造化変形可能メッシュモデルを改善する手段を更に有する。このシステムは、内部個別の要素に割り当てられる最適なサイズを評価するように表面個別の要素に関するサイズ情報(L)を取得するための、及び、改善処理中に新たな内部個別の要素が生成される間に、表面個別の要素から内部個別の要素にこのサイズ情報を伝えるための画像データ処理手段を更に有する一方、新たな内部個別の要素(TH)は、内部個別の要素の内部の新しい頂点の挿入による改善処理中に与えられる。

Description

本発明は、注目すべき対象物に可変メッシュモデルをマッピングする操作を有する二次元画像又は三次元画像において、注目すべき対象物のセグメント化のために画像データ処理手段を有する画像処理システムに関する。本発明は、異常又は病状を研究又は検出するために、体の組織又は体の流体流のような対象物のセグメント化のために、処理システムにより処理される医療用二次元又は三次元画像を生成するための医療用検査装置に更に関する。本発明は、医療用撮像方法、そのプログラムプロダクト及び装置又はシステムの分野における特定のアプリケーションに適する。
三次元において、四面体メッシュ、即ち、四面体を有するボリュームメッシュは、血管系中の血液流のような三次元対象物における物理量をモデル化するために主に用いられる。メッシュ要素の形状の適合度は、それが計算の正確度及び安定度に非常に影響するため、重要である。理想的な要素の形状は、等辺の面及び同じ辺の長さを有する規則的な四面体である。
四面体メッシュは三角形を構成する表面メッシュにより形成される。三角形メッシュは3Dの対象物の表面を表す一方、その四面体メッシュは同じ3Dの対象物内のボリュームを表す。それら両方のタイプのメッシュは同じ表面三角測量を共有する。
四面体メッシュの生成は、主に、所謂、ドロネー四面体測量方法に基づいている。そのドロネー方法は、例えば、“Reasonably efficient Delaunay based mesh generator in three dimensins”と題され、H.Borouchaki,F.Hecht,E.Salteloyobi P.L.George著の1995年8月23日に出版された文献(INRIA,Domaine de Rocquencourt,BP 105,78153 Le Chesnay Cedex FRANCE,EUROPE)において開示されている。
この方法に従って、四面体要素が、改善される必要がある四面体内図において、ドロネー基準による新たな頂点を追加的に挿入することにより与えられる。その方法は三角形から成るメッシュ表面から始まり、表面メッシュの頂点を有する一般的な頂点を有する四面体を用いて、粗いボリュームメッシュを更に生成する。次いで、このボリュームメッシュは、最適な要素サイズが得られるまで、ドロネー方法を用いて追加的に改善される。
その問題点は、最適な形状及びサイズを実際に示す要素の規定の仕方である。高速で且つ単純過ぎる方法は、ボリュームメッシュの各々の四面体要素に同じサイズを与えるようになる。しかしながら、この方法は、表面の三角形メッシュの局所的なサイズ変化を考慮していないため、非常に限定的であり、そのことは形の悪いボリューム要素をもたらす。
H.Borouchaki,F.Hecht,E.Salteloyobi P.L.George,"Reasonably efficient Delaunay based mesh generator in three dimensins",August 23,1995,INRIA,Domaine de Rocquencourt,BP 105,78153 Le Chesnay Cedex FRANCE,EUROPE
本発明の目的は、完全自動方法を実行するために画像データ処理手段を有する画像処理システムであって、3D画像におけるボリュームメッシュモデルか又は2D画像における内部メッシュモデルかどちらかを生成することができる、画像処理システムを提供することである。このボリュームメッシュモデルは、三角形から成る表面メッシュから生成される四面体要素を有し、表面三角形の局所的なサイズの変化に従った四面体要素サイズに自動的且つ動的に適合する。内部メッシュモデルは、線分から成る輪郭メッシュから与えられる三角形要素を有し、輪郭線分のサイズの局所的変化に三角形要素を適合させる。ボリューム四面体要素及び内部三角形要素は個別内部要素と更に呼ばれる一方、表面三角形要素及び輪郭線分要素は個別表面要素と呼ばれる。
本発明の目的は、個別内部要素のメッシュ品質を評価するように画像データ処理手段を有する画像処理システムを提供することである。本発明に従って、メッシュモデルは、前記個別内部要素の内側の新たな頂点の挿入により改善される。本発明のシステムは、その処理を改善するための処理手段であって、
個別内部要素に関連付けられるように最適サイズを評価するように、個別表面要素により規定されるサイズ情報を取得するための手段、及び
改善処理中に新たな個別内部要素が与えられる間に、個別表面要素から個別内部要素にこのサイズ情報を伝えるための手段、
を有する、処理手段を有する。
本発明の目的はまた、このシステムを操作するための段階を有する画像処理方法を提供することを目的とする。本発明はまた、3D画像処理のためにこのシステムに結合された医療用診断撮像装置に関する。医療用撮像装置は、MRI医療用検査装置、X線医療用検査装置又は何れの他の3D医療用撮像装置であることが可能である。本発明は、画像処理方法を実行するためのプログラムパッケージ又はプログラムプロダクトに更に関する。
本発明については、以下、図を参照して詳述する。
本発明は、研究される注目すべき対象物を表す医療用画像の改善に関する。注目すべき対象物は、例えば、二次元又は三次元医療用画像で表された腹部大動脈瘤(AAA)を研究するための腹部大動脈のような血管であることが可能である。
それらの画像は、特定の患者用の血流の演算流体動的(CFD)シミュレーション及び血流に対する血管系の短期間及び長期間反応による心臓血管疾患の研究及び検出のために用いられることが可能である。これに関連して、CFDシミュレーションは、有限要素法(FEM)により、血管成分についての幾何学的及び機械的情報をモデル化することにある。その幾何学的情報は、三次元表面メッシュ(ボクセル分類)の形で医療用画像のセグメント化からもたらされる。FEMについて、不可欠の段階は、有限ボリューム要素を構成するボリュームメッシュにおける表面メッシュのモザイク化である。この操作はボリュームメッシュ生成と呼ばれる。
三次元においては、有限ボリューム要素は、通常、四面体タイプ及び六面体タイプと呼ばれる2つの可能なタイプを有し、それらの各々は点及びそれらの点の間の接続の集合として表される。
有限ボリューム要素が六面体タイプである場合、構造化メッシュと呼ばれるボリュームメッシュモデルのタイプは要素タイプに関連付けられる。構造化メッシュは、各々の点で点と規則的な接続(即ち、一定隣接数、例えば、常に3つの隣接要素であって、それ以上でもそれ以下でもない)の集合を有する。
本発明は、六面体形状として知られている可能な形状には関連しない。その代わり、本発明に従って、有限ボリューム要素は六面体タイプを有する。六面体タイプの場合、非構造化メッシュと呼ばれるボリュームメッシュモデルは要素タイプに関連する。各々の点の接続は規則的でない(例えば、その数を変わることが可能であり、3つ、4つ、5つ又はそれ以上の隣接要素を求めることが可能である)。
非構造化メッシュの有利点は、四面体要素が適切な正確度を有する不規則な境界に適合するようにするそれらのフレキシビリティである。非構造化メッシュの他の有利点は、メッシュ適合要求を満足するそれらの能力である。実際には、正確度と計算時間との間のトレードオフを可能にするように、メッシュが制御可能であることがしばしば要求される。この場合、要素密度は、局所的な正確度の要求に応じて変化する必要があり、この変化はスムーズである必要がある。このことはメッシュ適合化と呼ばれる。非構造化要素を用いる場合、要素サイズ及び密度の変化は、接続性は制限されないために、制御されることが可能である。四面体要素については、演算の最高の正確度が規則的な四面体により得られる。十分な正確度を保証するように、そのメッシュは、要素の幾何学的形状の品質を測定する品質基準の最適条件、例えば、最小条件を満足する必要がある。
本発明は、三次元の個別の変形可能なボリュームメッシュモデルを用いて、三次元画像に表される注目すべき対象物を自動的にセグメント化するための画像処理システムの第1実施形態に関する。セグメント化のボリュームメッシュモデルの表面Sは前記三次元の対象物の表面に適合され、そのモデルのボリュームメッシュVは表面Sのメッシュに適合される。本発明に従って、四面体メッシュ、即ち、四面体を有するボリュームメッシュは、複数の三角形を有する表面メッシュから生成される。三角形メッシュは、注目すべき3Dの対象物の表面の表現である一方、四面体メッシュは、同じ3Dの対象物内のボリュームの表現である。両方のタイプのメッシュにおいては、同じ表面三角測量が共通である。理想的な要素形状は等辺の面及び同じ辺の長さを有する規則的な四面体である。
本発明は、二次元の個別の変形可能メッシュモデルを用いて、二次元画像で表される注目すべき対象物をセグメント化するための画像処理システムの他の実施形態に更に関する。このシステムは手段を有し、それにより、セグメント化の変形可能な2Dのメッシュモデルの輪郭Sのセグメントは、2Dの画像における前記対象物の境界において適合され、その輪郭の内部の三角形メッシュVはその輪郭のセグメントのサイズに適合される。注目すべき対象物は二次元の医療用画像で表される器官の断面である。
本発明に従って、三角形メッシュ、即ち、輪郭Sの内部メッシュはセグメントから成る輪郭から生成される。セグメント化された輪郭メッシュは、2D画像で表される注目すべき2Dの対象物の表面の表現である一方、三角形メッシュから成る2D領域は同じ2Dの対象物の輪郭内の領域の表現である。理想的な内部の要素の形状は正三角形である。
実際には、本発明は、三次元画像又は二次元画像において同じ問題点を解決するための手段を有する。本発明は、対象物の表面のセグメント化された輪郭に対して個別の内部要素のサイズの自動的最適化のための画像データセグメントか手段を有する画像処理システムを提供する。それらの個別内部要素は、三角形で構成される3Dセグメント化表面に対する3D四面体か又は、セグメントで構成される2Dセグメントの輪郭に対する2Dの正三角形のどちらかである。
第1実施形態は、ボリュームメッシュモデルを用いて3D対象物をモデル化することについて表現される。図1Aは、この第1実施形態に関する本発明のシステムの手段を表す図である。本発明に従って処理される画像は、図5A乃至図5C、図6A及び図6Bに示されている。図5Aは三角形から成る球の表面メッシュSを表し、図5Bは同じ球の内部のボリュームメッシュVを表し、それら両方の画像は、その球の内側を見るように面Pに対して切り取られている。複数の三角形を構成するセグメント化表面Sが先ず、利用可能である。3D対象物の前記表面メッシュSから、同じ対象物の最初のボリュームメッシュVが生成される。最初のボリュームメッシュVの四面体の頂点全ては表面Sの三角形の頂点である。それ故、四面体要素は対象物の表面に全て接続される。図5Bに示すように、このようなボリュームメッシュは非常に平坦な四面体要素であり、それらの四面体要素は四面体の低い形状品質をもたらし、四面体のサイズ及びボリュームにおいて大きい変化をもたらす。
図1Aを参照するに、本発明の自動システムは、先ず、非構造化ボリュームメッシュモデルを自動的且つ動的に構成するためのデータ処理手段を有するシステムであって、
1)3Dセグメント化表面Sの個別の表面要素を生成するための演算手段1Aであって、図5Aに示されるように、それらの表面要素は、辺で隣接するSにおける頂点で規定される三角形Tから構成される、演算手段1A;及び
2)最初のボリューム要素Vを生成するための演算要素2Aであって、それらのボリューム要素はTHで表される四面体であり、その四面体の4つの頂点はSの頂点であり、このことは、図5Bに示すように、そのような四面体は非常に平坦であることが可能であることを確実にし、そのようなことが粗いと言われる理由である、演算要素2A;
を有する。
自動システムはまた、最初のボリューム要素を改善するための演算手段を有し、その自動システムは:
表面要素のサイズ情報を取得するための評価手段3A;
Sの表面要素Tに関連するサイズ情報を用いて、Vのボリューム要素THに関連付けられる最適サイズを評価するための評価手段4A、5A;及び
改善処理中に新しい要素THが生成される間に、表面SからボリュームVにこのサイズ情報を伝えるための改善手段6A;
を有する。
本発明に従って、ボリューム要素は、表面要素のサイズ情報を考慮して、最初のボリューム要素の内側の新しい頂点の挿入により改善される。図1Aを参照するに、それらの改善処理手段は、好適には、次のように詳細に説明する構成要素を有する。
3)個別の要素の各々の頂点に関連する重みパラメータLを規定するための処理手段。例えば、頂点Bの集合の各々の頂点について、個別の要素の集合を表している図3を参照するに、隣の頂点に前記頂点を接続する辺の異なる長さ、例えば、セグメントの長さBA、BK、BG、BF、BE、BC、BDが演算され、距離Lで表される。それ故、最適距離Lと呼ばれる重みパラメータはBについて、及び、個別の要素の集合の他の頂点Jについて更に、演算される。重みパラメータは、好適には、Bに関連する、それ故、前記他の頂点Jに関連する異なる距離Lの平均値である。頂点Jに関連付けられる重みパラメータLを演算するこの操作は、Vを構成する図5Bの粗い四面体THで表される最初のボリューム要素の頂点に対して3Dの状態で適用され、前記頂点全ては三角形Tを構成するセグメントの表面S上にある。
4)各々の四面体要素THに関連する最適なボリュームVを演算するための処理手段4A。3Dの最初の四面体要素はそれぞれの3D表面メッシュの頂点に基づいている。3Dにおいて、四面体要素THは3D表面Sの4つの頂点に基づいていて、各々の頂点は、予め演算された最適な距離Lにより構成されるそれぞれの重みパラメータにより成るそれぞれの重みパラメータに関連している。前記距離LをメッシュVの頂点に割り当てることは、表面メッシュSおよびボリュームメッシュVは同じ表面三角測量を共有する。規則的な四面体である最適な要素形状、即ち、最適なボリュームVは、要素を有する頂点の4つの最適距離LJの平均に等しい辺の長さを有する規則的な四面体のボリュームである。ボリュームVは、次式(1a)で与えられる。
Figure 2007523402
5)各々の最初の四面体要素の真のボリュームvRJを演算するための処理手段5A。
6)実際のボリュームvRJと最適なボリュームvを比較するため、従って、検討中の四面体THの改善を開始するための処理手段6Aであって、a)実際のボリュームvRJが最適なボリュームvより大きい場合、本発明に従って、即ち、更なる処理手段7Aが用いて、検討中の四面体要素の改善を操作し、そうでない場合、b)ボリュームVの他の四面体にスキップされ、そして、c)それ以上、改善されるべき四面体がない場合、改善することを停止する。
7)実際のボリュームvRJが最適なボリュームvより大きい四面体において幾つかの頂点挿入の位置を選択するための処理手段7A。新しい頂点を挿入するための、幾つかの可能な位置は、図2Aで示されているような辺の一の中央、図2Bに示されているような面の中心、図2Cで示されている四面体の中心、及び図2Dで示されているような外接球の中心である。
8)選択された位置に依存して、新たに挿入された頂点に関連付けられる最適距離Lと呼ばれるパラメータを演算するための処理手段8A。その選択位置が、その四面体の辺の一の中央(図2A)にある場合、新たに挿入されるように割り当てられる最適距離は、予め演算され且つその辺の先端における2つの頂点に割り当てられた2つの最適な距離の平均であり、その四面体の面の一の中心(図2B)にある場合、新たに挿入されるように割り当てられる最適距離は、予め演算され且つその面の3つの先端に割り当てられた2つの最適な距離の平均であり、その四面体の中心(図2C)にある場合、新たに挿入されるように割り当てられる最適距離は、予め演算され且つその四面体の4つの頂点に割り当てられた4つの最適な距離の平均であり、その外接円の中心(図2D)にある場合、最適距離は四面体の4つの頂点の4つの最適距離の平均である。
9)挿入される頂点の各々の選択肢により与えられる形状を評価し、比較するように、各々の四面体の形状に基づいて、四面体の形状品質の指標qを演算するための指標手段9A。上記の四面体において頂点を挿入する4つの可能性の間の差異を評価するように、第1の有効な基準が次式により与えられ、
Figure 2007523402
ここで、ρは四面体要素の内接球の直径であり、hは四面体要素の最大の辺の長さである。
他の形状品質の簡単な基準qは次式のようであり、
Figure 2007523402
ここで、ρは四面体に内接された球の直径であり、dは外接球の直径である。生成される最も適切でない要素の最も適切な品質を与える挿入の位置は維持される。
10)選択された位置における新たな頂点の挿入によりメッシュを改善するための処理手段。四面体を改善することはそのボリュームの内側の最適なサイズの情報を伝えることを可能にし、最初の四面体を置き換えるための幾つかの小さい四面体を生成する。四面体を生成している間にまた、ドロネー妥当性基準が適用される。
ドロネー妥当性基準については、下で説明する。外接球、即ち、四面体の4つの点により規定される球が他のメッシュの点を有しない場合及びその場合に限り、四面体は、所謂“妥当なドロネー”と呼ばれる。拡大解釈すれば、全てのメッシュが妥当なドロネーである場合及びその場合にのみ、メッシュは妥当なドロネーである。この基準については、2D画像に関連して、図3B及び3Cにより示されている。新しい頂点Oが三角形ACDに挿入されている。しかしながら、この新しい頂点は、三角形ABCの外接円Φ1及び三角形CDEの外接円Φ3の内側にある。それ故、線分AC及びCDは削除され、新たな線分OA、OB、OC、OD、OEが形成される。これは、新たな三角形AOB、BOC、COE、DOE、AODの形成を可能にする。
図3B及び図3Cは、2D三角形メッシュにおいて新たな点がどのように挿入されるかを示し、その処理は、同じ方法で、3Dに拡張される。各々の新たに挿入された頂点は、メッシュに接続されるために必要である。その点に接続するために、外接球がその点で重なり合う、即ち、もはや妥当なドロネーではなく、メッシュからそれらの妥当なドロネーを取り除くように、全ての四面体を位置付ける。これは、図3Cにおいて2Dで示されているように、周囲を囲われた凸状キャビティを規定する。それ故、新たな四面体の構築が、キャビティ面に対して点を接続することにより実施される。
それ故、本発明の処理手段を用いて、完全自動方式が適用され、その完全自動方式は、表面の三角形のサイズの局所的変化に従って、四面体要素のサイズを動的に適合させる。
本発明の手段は2D画像に適用するために十分に適切である。第2実施形態については、2Dの変形可能メッシュモデルを用いて2Dの対象物をセグメント化するために説明する。図1Bは、この第2実施形態に関する本発明のシステムの手段を表す図である。2Dセグメント化については図4A乃至図4Cに示している。図4Aに示すように、注目すべき対象物は、線分ESを構成する輪郭メッシュSに従ってセグメント化される。2Dの対象物の前記輪郭メッシュSから、同じ対象物の最初の内部の2DのメッシュVはITにより表される三角形が与えられる。最初の内部の2DメッシュVの三角形の頂点の全ては輪郭Sの線分の頂点である。更に、このような最初の内部の三角形Vの集合は、輪郭メッシュSの頂点以外の頂点は含まない。それ故、三角形ITは全て、対象物の輪郭に接続される。図4Bに示すように、その結果は非常に粗いものである。Vの全ての頂点はまた、Sの頂点であるため、この内部メッシュは形状が決して等辺形状でない三角形要素を有し、それにより、メッシュの低い形状品質が得られる。
図1Bを参照するに、本発明の自動的及び動的システムは、最初に、2D輪郭メッシュS及び内部の個別の要素Vを自動的に構成するためのデータ処理手段を有する。それらの手段は図1Aの手段に相当し、
1)2Dの個別の輪郭要素Sを与えるための計算手段1Bであって、それらの輪郭要素は、図4Aに示すように、それらの辺により隣接する、Sにおける頂点A′、B′、C′、...、K′で規定される線分ESJを有する、計算手段1Bと、
2)最初の内部の個別の要素Vを与えるための計算手段2Bであって、複数の三角形を有し、図4Bに示すように、それら三角形の3つの頂点は、例えばA′B′D′のようなSの頂点である、計算手段2Bと、
を有する。
最初の内部の個別の要素Vを改善するための手段3B乃至11Bは、
輪郭要素のサイズ情報を取得するための評価手段3Bと、
Vの個別の内部要素に割り当てられるように最適なサイズを評価するためにSの輪郭要素により規定されるからサイズ情報を用いる計算手段4B、5Bと、
その改善処理中に新たな内部要素が与えられる間に、輪郭Sから内部領域Vにサイズ情報を伝えるための改善手段6B乃至11Bと
を有する。
本発明に従って、内部要素は、輪郭要素のサイズ情報を考慮して、最初の三角形要素の内側における新しい頂点の挿入により改善される。図1Bに示すように、それらの改善処理手段は、好適には且つ詳細には、次の構成要素を有する。
3)セグメント化された輪郭Sの各々の頂点に割り当てられた重みパラメータLを規定するための処理手段3B。
4)各々の三角形要素ITに割り当てられた最適な割り当て表面sを計算するための処理手段4B。2Dにおいて、ITで表される三角形要素は2D輪郭Sの3つの頂点に基づいていて、各々の頂点は、予め計算された最適な距離Lで構成されるそれぞれの重みパラメータを割り当てられる。正三角形、即ち、最適な領域sである最適な要素形状は、その要素を成す頂点の3つの最適な距離Lの平均に等しい辺の長さを有する正三角形の表面である。
5)各々の最初の三角形要素の実際の面積を計算するための処理手段5B。
6)実際の領域sRJと最適な領域sを比較するため、従って、検討中の三角形ITの改善を開始するための処理手段6Bであって、a)三角形ITの実際の領域sRJが最適な表面sより大きい場合、本発明に従って、即ち、更なる処理手段7Bを用いて、検討中の三角形要素の改善を操作し、そうでない場合、b)内部領域Vの他の三角形をスキップし、そしてc)改善されるべき三角形がもはやない場合、停止する。
7)実際の領域sRJが最適な領域sより大きい三角形への頂点の挿入の幾つかの位置を選択するための処理手段7B。新しい頂点を挿入するために、幾つかの有効な位置は、複数の辺の位置の中央、三角形の中心又は外接円の中心にある。
8)選択された位置に応じて、新たに挿入された頂点に割り当てられるように最適距離Lと呼ばれるパラメータを演算するための処理手段8B。選択された位置は、辺の位置の中央にある場合、新たに挿入される頂点に割り当てられる最適距離はその辺の先端の2つの頂点に対して割り当てられ且つ予め計算された2つの最適な距離の平均であり、三角形の中心である場合、挿入される頂点に割り当てられる最適距離はその三角形の3つの頂点に割り当てられ且つ予め計算された3つの最適距離の平均であり、外接円の中心にある場合。
9)頂点を挿入する各々のオプションにより与えられる形状を評価し且つ計算するように、三角形の形状品質の指標qを計算するための指標手段9B。
10)選択された位置における新しい頂点の挿入によりメッシュを改善するための処理手段10B。
それ故、本発明の処理手段を用いて、完全自動方法を適用し、その完全自動方法は、輪郭の辺のサイズの局所的変化に従った三角形の要素のサイズに動的に適合する。
医療用検査装置及び観察システム
上記の手段は、本発明の観察システム内に含まれるかその観察システムに結合されている。図7は、医療用検査装置内に組み込まれている本発明に従った画像観察システムの実施形態の基本構成要素を示している。医療用撮像装置100は、CTスキャナ又はX線装置又は超音波装置のような他の医療用撮像装置であることが可能である。その装置100により得られる画像データは汎用コンピュータのようなデータ処理手段70に供給され、そのデータ処理手段は本発明のインタラクティブな適合手段を構成するユーザ制御手段及び計算手段を有する。データ処理手段70は、典型的には、モニタ60のような視覚化装置と、ユーザがそのシステムと対話することができるようにユーザが操作可能であるキーボード又はマウス71、ポインティング装置70等のような入力装置72とに結合されている。データ処理装置70は、本発明に従って医療用画像データを処理するための処理手段を実施するようにプログラムされている。特に、データ処理装置70は、図1及び図4に関連して述べた操作を実行するために必要な計算手段及びメモリ手段を有する。それらの操作を実行するために予めプログラムされた命令を有するコンピュータプログラムプロダクトがまた、実行されることができる。
図面及びそれら図面に関連する上記の説明は本発明を限定するのではなく、例示である。特許請求の範囲における範囲に網羅される多くの変形が存在することが理解されるであろう。更に、本発明は、ディスプレイのための画像データを生成することに関して述べたが、本発明は、表示装置及びプリンティングを含む画像データの視覚化の何れの形式を実質的にカバーすることが意図されているが、それらに限定されるものではない。
ボリュームメッシュモデル及び表面メッシュモデルを用いてそれぞれ、対象物をモデル化するための本発明のシステムの手段を表している図である。 ボリュームメッシュモデル及び表面メッシュモデルを用いてそれぞれ、対象物をモデル化するための本発明のシステムの手段を表している図である。 変形可能ボリュームメッシュモデルにおいて四面体を改善するために頂点挿入について示す図であって、辺の中央に頂点を挿入するオプションを示す図である。 変形可能ボリュームメッシュモデルにおいて四面体を改善するために頂点挿入について示す図であって、四面体の三角形の面の中心点に頂点を挿入するオプションを示す図である。 変形可能ボリュームメッシュモデルにおいて四面体を改善するために頂点挿入について示す図であって、四面体の中心点に頂点を挿入するオプションを示す図である。 変形可能ボリュームメッシュモデルにおいて四面体を改善するために頂点 挿入について示す図であって、四面体の外接球の中心点に頂点を挿入するオプションを示す図である。 メッシュを改善するためのドロネー基準の適用を示す図であって、頂点の重みの決定を示す図である。 メッシュを改善するためのドロネー基準の適用を示す図であって、三角形内への頂点の挿入を示す図である。 メッシュを改善するためのドロネー基準の適用を示す図であって、三角形内への頂点の挿入を示す図である。 2D対象物のセグメント化された輪郭を示す図である。 輪郭から構成される2Dの個別の内部要素を示す図である。 改善された2Dの個別の内部要素を示す図である。 表面メッシュが三角形の集合から成る場合、3Dの対象物のセグメント化された表面を示す図である。 粗い四面体から構成されるボリュームメッシュモデルを示す図であって、そのボリュームメッシュモデルの頂点はこの表面メッシュモデルの頂点である、図である。 このボリュームメッシュモデルを示す図であって、そのボリュームメッシュモデルの四面体は本発明に従って改善され、表面の三角形に適合された小さい四面体を有する、図である。 表面メッシュが他の三角形の集合を有するセグメント化された表面を示す図である。 本発明に従ってボリュームメッシュが改善されるセグメント化ボリュームであって、それらの他の表面の三角形に適合された小さい四面体を有する、セグメント化ボリュームを示す図である。 医療用検査装置に結合された医療用観察システムを示す図である。

Claims (14)

  1. 表面個別要素及び内部個別要素を有する非構造化変形可能メッシュモデルを用いて、注目すべき対象物のセグメント化の画像データ処理手段を有し、前記表面個別要素のサイズの局所的変化に従って、前記内部個別要素のサイズを自動的且つ動的に適合させることにより非構造化変形可能メッシュモデルを改善する手段を更に有することを特徴とする画像処理システム。
  2. 請求項1に記載の画像処理システムであって、前記内部個別要素に割り当てられる最適なサイズを評価するように前記表面個別要素に関するサイズ情報を取得するための、及び、改善処理中に新たな内部個別要素が生成される間に、前記表面個別要素から前記内部個別要素にこのサイズ情報を伝えるための画像データ処理手段を更に有する、ことを特徴とする画像処理システム。
  3. 請求項2に記載の画像処理システムであって、新たな内部個別要素は、前記内部個別要素の内部の新しい頂点の挿入による改善処理中に与えられる、ことを特徴とする画像処理システム。
  4. 請求項3に記載の画像処理システムであって、前記内部個別要素のメッシュ品質を評価するように及び前記評価されたメッシュ品質に基づいて非構造化メッシュモデルを改善するように画像データ処理手段を有する、ことを特徴とする画像処理システム。
  5. 請求項4に記載の画像処理システムであって、前記非構造化メッシュモデルは、三角形を有する表面個別要素及び四面体を有する内部個別要素を有する三次元メッシュモデルであり、又は、非構造化メッシュモデルは、輪郭線分を有する表面個別要素及び三角形を有する内部個別要素を有する表面個別要素を有する二次元メッシュモデルである、ことを特徴とする画像処理システム。
  6. 請求項5に記載の画像処理システムであって、三次元においては、内部四面体は、はじめに表面の三角形の頂点に基づいて構成され、次いで、四面体の辺の中央、四面体の面の中心、四面体の中心又は四面体の外接円の中心の何れに頂点を挿入することにより改善され、又は、二次元においては、内部三角形は、はじめに輪郭要素の頂点に基づいて構成され、次いで、三角形の辺の中央、三角形の面の中心又は三角形の中心の何れに頂点を挿入することにより改善される、ことを特徴とする画像処理システム。
  7. 請求項5又は6に記載の画像処理システムであって、隣の頂点に前記頂点を結合する辺の長さの平均に基づいて、前記個別要素の各々の頂点に関連する重みパラメータを評価するための画像データ処理手段と、各々の内部個別要素に関連する最適なボリューム又は表面であって、最適な内部個別要素の形状は規則的な四面体又は三角形及び各々のはじめの内部個別要素の実際のボリューム又は表面である、最適なボリューム又は表面と、前記の最適なボリューム又は表面と前記実際のボリューム又は表面をそれぞれ比較するため、従って、前記内部個別要素の前記実際のボリューム又は表面が最適なボリューム又は表面より大きい場合、検討中の内部個別要素の改善を開始するための画像データ処理手段と、を有する、ことを特徴とする画像処理システム。
  8. 請求項7に記載の画像処理システムであって、有効性基準を評価するための画像データ処理手段であって、その有効性基準に従って、外接球又は外接円がメッシュの他の頂点を囲んでいない場合、新しい内部要素が有効である、画像データ処理手段を有する、ことを特徴とする画像処理システム。
  9. 請求項7に記載の画像処理システムであって、前記内部個別要素のメッシュ品質を評価するための前記画像データ処理手段は、外接球又は外接円の直径及び前記内部個別要素の辺の長さに基づく基準と、前記内部個別要素のボリューム又は表面に基づく基準とを有する、ことを特徴とする画像処理システム。
  10. 請求項1乃至9に記載の画像処理システムであって、処理された画像を表示するための視覚化手段を更に有する、ことを特徴とする画像処理システム。
  11. 請求項1乃至10に記載の画像処理システムであって、メッシュ品質の所定の閾値が適合されるとき、内部個別要素の改善を停止させるための手段を更に有する、ことを特徴とする画像処理システム。
  12. 回路手段を有する特殊目的のプロセッサ又は適切にプログラムされたコンピュータを有する医療用撮像システムであって、医療用画像データを処理するように、請求項1乃至11の何れ一項に記載の画像処理システムを構成するように備えられている、ことを特徴とする医療用撮像システム。
  13. 医療用検査撮像装置であって:
    体の器官の三次元画像を取得するための手段;及び
    請求項1乃至12の何れ一項に記載のシステム;
    を有することを特徴とする医療用検査撮像装置
  14. 請求項1乃至12の何れ一項に記載のシステムで用いられる指令の集合を有することを特徴とするコンピュータプログラム。
JP2006548466A 2004-01-13 2005-01-05 内部個別要素を用いるメッシュモデル Pending JP2007523402A (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP04300019 2004-01-13
PCT/IB2005/000051 WO2005069227A1 (en) 2004-01-13 2005-01-05 Mesh models with internal discrete elements

Publications (1)

Publication Number Publication Date
JP2007523402A true JP2007523402A (ja) 2007-08-16

Family

ID=34778232

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2006548466A Pending JP2007523402A (ja) 2004-01-13 2005-01-05 内部個別要素を用いるメッシュモデル

Country Status (5)

Country Link
US (1) US20070165948A1 (ja)
EP (1) EP1706848A1 (ja)
JP (1) JP2007523402A (ja)
CN (1) CN1910625A (ja)
WO (1) WO2005069227A1 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008134224A (ja) * 2006-08-31 2008-06-12 Sierra Nevada Corp 三次元レーダー像のレンダリングのための装置及び方法
JP2011500187A (ja) * 2007-10-19 2011-01-06 ヴァスコップス 管状構造用の自動的な幾何学的及び機械的解析方法及びシステム
JP2011138350A (ja) * 2009-12-28 2011-07-14 National Institute Of Information & Communication Technology 対象体の解剖構造解析方法及び対象体解剖構造の表示方法並びに対象体解剖構造表示装置
JP5022508B1 (ja) * 2011-06-02 2012-09-12 三菱重工業株式会社 形状計測方法
JP2013254445A (ja) * 2012-06-08 2013-12-19 Fujitsu Ltd 描画プログラム、描画装置、および描画方法

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3840173B2 (ja) * 2002-11-15 2006-11-01 キヤノン株式会社 三次元解析用メッシュ生成方法、三次元解析用メッシュ生成装置、プログラムおよび記憶媒体
GB0305315D0 (en) * 2003-03-07 2003-04-09 Weber Martin Image processing system
US7702490B1 (en) * 2004-11-16 2010-04-20 Livermore Software Technology Corporation Method and system for adaptive mesh-free shell structures
TWI275971B (en) * 2005-10-27 2007-03-11 Coretech Sys Co Ltd Automated meshes creation method
JP4999522B2 (ja) * 2007-04-06 2012-08-15 株式会社日立製作所 解析メッシュ生成装置
US8209157B2 (en) * 2009-05-06 2012-06-26 Livermore Software Technology Corporation Methods and systems for simulating beam-to-surface contacts in finite element analysis
WO2011127940A1 (en) 2010-04-13 2011-10-20 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Method for segmentation of objects from a three dimensional image data set with a deformable model
CN102034272B (zh) * 2010-09-29 2012-07-18 浙江大学 个性化颌面部软组织六面体网格的生成方法
JP6049272B2 (ja) * 2012-02-24 2016-12-21 キヤノン株式会社 メッシュ生成装置、方法およびプログラム
FR2999325A1 (fr) * 2012-12-07 2014-06-13 Commissariat Energie Atomique Procede iteratif de determination d'une image en deux dimensions ou trois dimensions a partir de signaux issus de tomographie par rayons x
US9633459B2 (en) * 2013-08-23 2017-04-25 Disney Enterprises, Inc. Methods and systems for creating a hull that may have concavities
US9773331B2 (en) 2013-08-23 2017-09-26 Disney Enterprises, Inc. Methods and systems for efficient graphics rendering
CN105014971B (zh) * 2015-07-31 2017-02-01 创生医疗器械(中国)有限公司 一种3d打印方法
EP3353748B1 (en) * 2015-09-23 2022-08-03 Koninklijke Philips N.V. Generation of triangle mesh for a three dimensional image
US20170116779A1 (en) * 2015-10-26 2017-04-27 Microsoft Technology Licensing, Llc Volumetric representation of objects
CN105427381B (zh) * 2015-11-16 2018-06-19 中国石油大学(华东) 空间球外联的岩心骨架提取方法
CN107423498B (zh) * 2017-07-13 2020-03-10 山东大学 一种高致密度离散颗粒多相体系的建模方法
KR101991666B1 (ko) * 2017-12-06 2019-10-01 주식회사 카이 구형 영상을 생성하는 방법, 구형 영상을 재생하는 방법 및 그 장치들
CN108520138B (zh) * 2018-04-02 2022-03-11 重庆大学 基于渐进有限元网格和分割法的块体离散元模型生成方法
ES2913434T3 (es) * 2018-12-03 2022-06-02 Siemens Ag Reconociendo desviaciones en el comportamiento de la seguridad de unidades automatizadas
US11995854B2 (en) * 2018-12-19 2024-05-28 Nvidia Corporation Mesh reconstruction using data-driven priors
CN110115840B (zh) * 2019-05-16 2021-10-29 腾讯科技(深圳)有限公司 虚拟场景中的物体移动控制方法、装置、终端及存储介质
CN113405966B (zh) * 2021-06-08 2022-08-23 浙江广天构件集团股份有限公司 一种水泥基材料颗粒堆积体系孔径分布计算方法
CN114036773B (zh) * 2021-11-25 2023-06-30 河海大学 一种块体离散元数值模型外轮廓几何快速检索方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6879324B1 (en) * 1998-07-14 2005-04-12 Microsoft Corporation Regional progressive meshes
US6816820B1 (en) * 1999-09-24 2004-11-09 Moldflow Ireland, Ltd. Method and apparatus for modeling injection of a fluid in a mold cavity
US6560570B1 (en) * 1999-11-22 2003-05-06 Sandia Corporation Method and apparatus for connecting finite element meshes and performing simulations therewith
WO2001048699A1 (en) * 1999-12-27 2001-07-05 Alcoa Nederland B.V. Mesh generator for and method of generating meshes in an extrusion process
US20020177985A1 (en) * 2001-04-23 2002-11-28 Kraft Joseph Anthony Computer system and method for radial cooled bucket optimization
WO2002097735A1 (en) * 2001-05-31 2002-12-05 Kent Ridge Digital Labs System and method of anatomical modeling
US6987511B2 (en) * 2002-10-17 2006-01-17 International Business Machines Corporation Linear anisotrophic mesh filtering
CN1764917A (zh) * 2003-02-05 2006-04-26 莫尔德弗洛爱尔兰有限公司 采用混合模型进行过程模拟的设备和方法
US20050143965A1 (en) * 2003-03-14 2005-06-30 Failla Gregory A. Deterministic computation of radiation doses delivered to tissues and organs of a living organism
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
US20050096515A1 (en) * 2003-10-23 2005-05-05 Geng Z. J. Three-dimensional surface image guided adaptive therapy system
US7454319B2 (en) * 2003-11-19 2008-11-18 Jun Wan System, method, and computer program product for determining wall thickness in graphic model

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008134224A (ja) * 2006-08-31 2008-06-12 Sierra Nevada Corp 三次元レーダー像のレンダリングのための装置及び方法
JP2011500187A (ja) * 2007-10-19 2011-01-06 ヴァスコップス 管状構造用の自動的な幾何学的及び機械的解析方法及びシステム
JP2011138350A (ja) * 2009-12-28 2011-07-14 National Institute Of Information & Communication Technology 対象体の解剖構造解析方法及び対象体解剖構造の表示方法並びに対象体解剖構造表示装置
JP5022508B1 (ja) * 2011-06-02 2012-09-12 三菱重工業株式会社 形状計測方法
WO2012164987A1 (ja) * 2011-06-02 2012-12-06 三菱重工業株式会社 形状計測方法
JP2013254445A (ja) * 2012-06-08 2013-12-19 Fujitsu Ltd 描画プログラム、描画装置、および描画方法

Also Published As

Publication number Publication date
US20070165948A1 (en) 2007-07-19
CN1910625A (zh) 2007-02-07
EP1706848A1 (en) 2006-10-04
WO2005069227A1 (en) 2005-07-28

Similar Documents

Publication Publication Date Title
JP2007523402A (ja) 内部個別要素を用いるメッシュモデル
US8384716B2 (en) Image processing method
US7215810B2 (en) Method for creating single 3D surface model from a point cloud
US7773786B2 (en) Method and apparatus for three-dimensional interactive tools for semi-automatic segmentation and editing of image objects
US20110262015A1 (en) Image processing apparatus, image processing method, and storage medium
US20050018885A1 (en) System and method of anatomical modeling
JP6385318B2 (ja) 3d医用画像中の対象物をセグメンテーションするための3d対象物の変換
JP2007528529A (ja) 3dデータ集合のサーフェスを識別するための方法及びシステム(「ボクセル分割」)
JP2004518473A (ja) 変形する3d物体の画像シーケンスを該物体の壁の動きの指示を伴って表示する画像処理方法
JP2007537811A (ja) 結節の検出方法
JP2015531946A (ja) 物体の3dモデリング
JP2005198708A (ja) 血管狭窄率解析装置及び血管狭窄率解析方法
JP2007518484A (ja) 変形可能な表面のセグメント化のリアルタイムなユーザ対話処理
US9836891B2 (en) Shape data generation method and apparatus
EP1685534B1 (en) Three-dimensional segmentation using deformable surfaces
JP6643821B2 (ja) 画像処理装置
JP4170096B2 (ja) 対象の3次元表面上にマップされた3次元メッシュモデルの適合性評価のための画像処理装置
JP4910478B2 (ja) モデリング装置、モデリング方法およびプログラム
JP2007267979A (ja) 生物の臓器形態解析方法と生物の臓器形態解析システム
EP2876607B1 (en) Shape data generation method and device
Caliskan et al. Three-dimensional modeling in medical image processing by using fractal geometry
CN111833448A (zh) 源自局部体积的腔室重建
EP4238534A1 (en) Method for determining object area from three-dimensional model, and three-dimensional model processing device
JP7382059B2 (ja) 画像処理装置、画像処理方法、およびプログラム
EP2991033A1 (en) System and method for three-dimensional shape generation from closed curves