JP4170096B2 - 対象の3次元表面上にマップされた3次元メッシュモデルの適合性評価のための画像処理装置 - Google Patents
対象の3次元表面上にマップされた3次元メッシュモデルの適合性評価のための画像処理装置 Download PDFInfo
- Publication number
- JP4170096B2 JP4170096B2 JP2002578256A JP2002578256A JP4170096B2 JP 4170096 B2 JP4170096 B2 JP 4170096B2 JP 2002578256 A JP2002578256 A JP 2002578256A JP 2002578256 A JP2002578256 A JP 2002578256A JP 4170096 B2 JP4170096 B2 JP 4170096B2
- Authority
- JP
- Japan
- Prior art keywords
- gradient
- mesh model
- image processing
- processing apparatus
- color
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite 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)
- Processing Or Creating Images (AREA)
- Image Processing (AREA)
Description
本発明は、3次元物体上に3次元メッシュモデルをマッピングする処理を含む、3次元画像における3次元物体分割用画像処理装置に関する。本発明は、また、3次元メッシュモデルによりマップされた分割3次元物体を、当該物体に対する当該メッシュモデルの適合性の可視的な指示と共に、表示するための画像処理装置に関する。本発明は、更に、器官の病状を検出若しくは研究すべく体の器官である物体の分割のための医療用イメージング装置若しくはシステム、及び、これらの装置若しくはシステムにより生成される医療用3次元画像を処理するプログラムに関する。
3次元物体のモデル化技術は、既に、“processing of the International Conference on Computer Vision and Pattern Recognition (CVPR’94), 20−24 June 1994, Seattle, USA”において“Simplex Meshes: a General Representation for 3D shape Reconstruction”という刊行物にH.DELINGETTEにより開示されている。この文献では、3次元物体を回復するための物理的ベースのアプローチが示されている。このアプローチは、“(シンプレックス・メッシュ)Simplex Meshes”の形状に基づく。メッシュの弾性挙動は、各頂点(メッシュの節点(ノード))で抽出されるシンプレックス角度により平均曲率を調整する局所安定化関数によって、モデル化される。この関数は、視点が一定で、固有なもので、且つ、スケールに影響されやすいものである。規則的なグリッドに定義される変形可能な表面と異なり、シンプレックス・メッシュは、非常に適応性のある構造である。大きく湾曲した部分や不精密な部分でメッシュの分解能を増加させるための細分化(改善)処理も開示されている。複雑なモデルを回復するためにシンプレックス・メッシュを接続する処理は、より単純な形状を有する部分を使用して実行されてよい。
本発明は、例えばグレーレベルで表現された3次元デジタル画像に適用される画像処理方法に関する。画像は、対象物体と称する器官の荒い3次元表面を表現してよい。例えば背景に対して、対象物体のより良好なビュー(表示)をユーザに提供すべく、この対象は分割される。分割画像により、ユーザが器官の異常や病気をより良好に研究し若しくは検出することができるようになる。本画像処理方法は、以下の幾つかのステップを含む。
3次元画像を取得する方法は、本発明の部分でない。分割方法は、当業者に知られる超音波システムやX線装置や他のシステムによって取得できる器官の3次元デジタル画像に適用できうる。3次元物体は、図1Aに示される。本例では、この対象は、立方体Cである。この立方体表面は、本方法が非常に多数の異なる且つ複雑な表面に適用されてよいことを実証すべく、選ばれている。対象の3次元物体を表現する3次元画像の取得後、当該画像が分割される。先行技術として上で引用した刊行物に関連して説明された分割技術は、ロバストで良好な結果を与えるゆえに、採用される。それは、シンプレックス・メッシュモデルと称される離散的モデルを用いて対象物を表現できる反復法則である。
3次元デジタルシンプレックス・メッシュモデルは、図1Bに示されるように生成される。本例では、それは、セルと称される小さな3次元離散曲面F0,F1,F2,....のセットから形成される単純な球M0であり、セルは、メッシュモデルのエッジと称される境界によりリンクしており、メッシュモデルの頂点と称される共通のノードを有する。分割処理は、本例では立方体である図1Aの3次元対象物体C上に、図1Bの3次元シンプレックス・メッシュモデルM0をマッピングする処理を含む。立方体表面の分割は、立方体表面が鋭い角を有する故に、3次元シンプレックス・メッシュモデルに対しては困難な図形であることを意味する。
この分割処理は、シンプレックス・メッシュモデルの元の球形状M0を対象物体C上にマップするために、即ち、その表面を可能な限り対象物体Cの表面に近づけるために、その球形状M0を変形させることを含む。この処理は、引用刊行物によって教示されるような反復法則に従った反復ステップによって実行される。この法則は、対象Cの表面に向けた、モデルのセルF0,F1,F2,....の引っ張りの第1の力である外力と、即ち、対象表面にセル表面を近づけさせる外力と、メッシュモデルの一般面を円滑化するための調整力である内力との間の釣り合いを確立する。
図2A及び図3Aは、上で引用した反復法則に従って実行されたある数の反復ステップ後の変形したメッシュモデルM0を示す。この段階でのメッシュモデルの新たな形状は、M1により指示される。初期のメッシュモデルM0のセルの表面は、外力の作用により対象Cの表面に引き付けられる一方、内力が、メッシュモデル表面を円滑化して、メッシュモデルM1の形状を、対象の形状により近づけるようにする。ユーザは、対象Cの形状に対するメッシュモデルM1の新たな形状の適合性を、当該対象Cと当該メッシュモデルM1の重畳画像を表示することによってのみ、評価することができる。更に、Cが高濃度の画像である際、可視的な評価が、一連の二次元スライスから通常的に実行できるだけである。従って、定量化された評価が、当該適合性をより良く且つより迅速に認識すべく必要となる。本方法は、3Dでの当該適合性の可視的な定量化を提供するための自動化技術を提案する。
セルの位置での勾配ベクトルの分布に基づく統計結果、
勾配ベクトルの長さでなく向き、若しくは、
勾配ベクトルのべき関数(パワー関数)に基づく。
b)平行四辺形分解を用いて三角形に沿って積分する。サンプリングステップは、既知のボクセルサイズよりも大きな副セルを有しないように選択される。分解のための技術は、次のステップを含む。
b.1)線分T1T2を要素ベクトル
b.2)aとbの間の所与のduに対して、適合する平行四辺形の最大数が求められ、即ち、bb’間隔が、duがT1T2から求められたのと同様の態様で、要素dvを決定すべく分割される。
b.3)各平行四辺形に対して、
と共に可変の“流れ(Flow)”を増加することによって維持される。所与の点に対する勾配値は補間されない。代わりに、最も近い隣接値が取られる。したがって、セルの領域(Area)により正規化された算出流れは、次の式で与えられる。
ステップ4.5)、4.6)及び4.7)を実行することによる上述の如く適合性の最初の推定後、ユーザは、この適合性をより良好にすべく反復ステップを続けることを決めても良い。
Claims (15)
- 3次元画像中の3次元物体の分割のための画像処理装置であって、該3次元物体上に3次元メッシュモデルをマッピングする処理を含む画像処理装置において、
分割されるべき対象物体の3次元画像を取得する手段と、
3角形セルに分解できるセルからなる、メッシュモデルを生成する手段と、
前記対象物体上に前記メッシュモデルをマップするためにメッシュモデルを変形させる手段と、
前記メッシュモデルの所定数のセルのセル表面領域を通る勾配ベクトル場の勾配流値若しくは勾配導出測度レベルを推定する推定手段と、
前記勾配流値若しくは勾配導出測度レベルが適合性閾値と称する少なくとも所定レベルに到達するセルの割合に基づいて前記メッシュモデルの適合性の良好性を評価する評価手段とを含む、画像処理装置。 - 所定の色が所与の勾配流値若しくは勾配導出測度レベルに対応する色分けテーブルを構築する手段と、
前記メッシュモデルの所与のセルの勾配流値若しくは勾配導出測度レベルを、該勾配流値若しくは勾配導出測度レベルに対応する前記色分けテーブルにより与えられる色に関連付ける手段とを更に含む、請求項1記載の画像処理装置。 - 前記所与のセルに、該セルの勾配流値若しくは勾配導出測度レベルに対応する前記色分けテーブルから決定される色を割り当てることによって、色分け処理を実行する手段と、
前記色分け処理により色付けられたセルを有する前記メッシュモデルの画像を表示する手段とを更に含む、請求項2記載の画像処理装置。 - 前記色分け処理は、全てのセル若しくは所定数のセルに対して実行される、請求項3記載の画像処理装置。
- 前記色分けテーブルにおいて、所与の色は、勾配流値若しくは勾配導出測度レベルの範囲に対応して作成されるか、若しくは、色又は色相の所定のスケール内にある色を持つセルの割合に対応して作成される、請求項3記載の画像処理装置。
- 前記色分けテーブルにおいて、所与の色相は、勾配流値若しくは勾配導出測度レベルの範囲の副分割に対応して作成される、請求項2乃至4の何れか記載の画像処理装置。
- 前記評価手段の結果に基づいて対象物体上へのメッシュモデルのマッピング処理を終了することを決定する手段を更に含む、請求項1乃至6の何れか記載の画像処理装置。
- メッシュモデルの所定数のセル若しくは所定割合のセルが、閾値と称される勾配流値若しくは勾配導出測度レベルの所与のレベルに達していない間、対象物体上へのメッシュモデルのマッピング処理を改善することを決定し、該閾値に達したときに処理を終了することを決定する手段を更に含む、請求項1乃至6の何れか記載の画像処理装置。
- メッシュモデルの所定数のセル若しくは所定割合のセルが、閾値と称される勾配流値若しくは勾配導出測度レベルの所定のレベルに対応する所与の色範囲内で表示されていない間、対象物体上へのメッシュモデルのマッピング処理を改善することを決定し、該所与の色範囲内に達したときに処理を終了することを決定する手段を更に含む、請求項3乃至6の何れか記載の画像処理装置。
- 前記改善する処理は、セルを2つに分割する処理を含む、請求項8又は9記載の画像処理装置。
- 勾配導出測度は、前記勾配ベクトル場の分散に基づく統計結果、若しくは、勾配ベクトルの長さでなく向き、又は、勾配のべき関数に基づく、請求項1乃至10の何れか記載の画像処理装置。
- 前記推定手段は、該セルを三角形に副分割し、平行四辺形分解を用いて前記三角形に沿った積分を実行し、前記セル領域に比例する前記勾配流値若しくは勾配導出測度レベルを供給する手段を含む、請求項1乃至11の何れか記載の画像処理装置。
- 請求項1乃至12の何れか記載の画像処理装置の各手段を実現するように画像データを処理するように構成された、適切にプログラムされたコンピューター、若しくは、回路手段を有する特殊用途用のプロセッサを含む、システム。
- 体の器官の3次元画像を取得する手段と、
請求項1乃至12の何れか記載の画像処理装置によりシンプレックス・セルから形成されるシンプレックス・メッシュモデルを生成し、対象物体上に前記メッシュモデルをマップするためにメッシュモデルを変形させ、前記メッシュモデルの所定数のセルのセル表面領域を通る勾配ベクトル場の勾配流値若しくは勾配導出測度レベルを推定し、前記勾配流値若しくは勾配導出測度レベルが適合性閾値と称する少なくとも所定レベルに到達するセルの割合に基づいて前記メッシュモデルの適合性の良好性を評価することを含む、前記3次元の器官画像上に3次元メッシュモデルをマッピングする処理を含む3次元画像内の器官の分割方法を実行する処理手段と、
医療用デジタル画像及び処理されたデジタル画像を表示する表示手段とを有する、医療診断用イメージング装置。 - コンピューターをして請求項1乃至12の何れか記載の装置を実現させるための指令のセットを含むコンピュータープログラム。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| EP01400817 | 2001-03-29 | ||
| PCT/IB2002/000894 WO2002080110A1 (en) | 2001-03-29 | 2002-03-19 | Image processing method for fitness estimation of a 3d mesh model mapped onto a 3d surface of an object |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2004521425A JP2004521425A (ja) | 2004-07-15 |
| JP4170096B2 true JP4170096B2 (ja) | 2008-10-22 |
Family
ID=8182666
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002578256A Expired - Fee Related JP4170096B2 (ja) | 2001-03-29 | 2002-03-19 | 対象の3次元表面上にマップされた3次元メッシュモデルの適合性評価のための画像処理装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US20020172406A1 (ja) |
| EP (1) | EP1374181A1 (ja) |
| JP (1) | JP4170096B2 (ja) |
| WO (1) | WO2002080110A1 (ja) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1588325B1 (en) * | 2003-01-15 | 2006-10-11 | Koninklijke Philips Electronics N.V. | Image processing method for automatic adaptation of 3-d deformable model onto a subtantially tubular surface of a 3-d object |
| US8276091B2 (en) * | 2003-09-16 | 2012-09-25 | Ram Consulting | Haptic response system and method of use |
| JP2007512042A (ja) * | 2003-11-17 | 2007-05-17 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 診断画像形成における特徴検出のためにユーザ入力を利用する方法 |
| WO2005059831A1 (en) * | 2003-12-11 | 2005-06-30 | Philips Intellectual Property & Standards Gmbh | Elastic image registration |
| EP1828991A1 (en) * | 2004-12-17 | 2007-09-05 | Koninklijke Philips Electronics N.V. | High quality accurate surface triangulation from a simplex mesh |
| ATE520101T1 (de) * | 2005-02-11 | 2011-08-15 | Koninkl Philips Electronics Nv | Bildverarbeitungsvorrichtung und -verfahren |
| US9129387B2 (en) * | 2005-06-21 | 2015-09-08 | Koninklijke Philips N.V. | Progressive model-based adaptation |
| US7817841B2 (en) * | 2005-11-12 | 2010-10-19 | General Electric Company | Time-lapse cell cycle analysis of unstained nuclei |
| JP5586953B2 (ja) * | 2006-09-29 | 2014-09-10 | コーニンクレッカ フィリップス エヌ ヴェ | 解剖学的形状情報を用いた医学画像データベースへのアクセス |
| CN101903911A (zh) * | 2007-12-20 | 2010-12-01 | 皇家飞利浦电子股份有限公司 | 身体和身体轮廓的3d重建 |
| US20110227923A1 (en) * | 2008-04-14 | 2011-09-22 | Xid Technologies Pte Ltd | Image synthesis method |
| JP5451135B2 (ja) * | 2009-03-26 | 2014-03-26 | キヤノン株式会社 | 画像処理装置、画像処理方法 |
| US9508162B2 (en) | 2009-08-01 | 2016-11-29 | Tsinghua University | Method and system for rapidly vectorizing image by gradient meshes based on parameterization |
| JP6232054B2 (ja) * | 2012-05-31 | 2017-11-15 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | 品質メトリックを決定するための方法、システム、コンピュータプログラム及び記憶媒体 |
| RU2015125701A (ru) * | 2012-11-30 | 2017-01-10 | Конинклейке Филипс Н.В. | Количественное определение шероховатости поверхности ткани на основе данных изображения и определение присутствия заболевания на их основе |
| CN107452047B (zh) * | 2017-08-09 | 2018-10-12 | 安徽中科超安科技有限公司 | 一种自适用分辨率的蒙特卡罗几何截面可视化方法 |
| CN119043434A (zh) * | 2024-07-23 | 2024-11-29 | 武汉大学 | 基于Argo浮标与海底地形数据的洋流流量方法及装置 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5189626A (en) * | 1991-03-27 | 1993-02-23 | Caterpillar Inc. | Automatic generation of a set of contiguous surface patches on a computer modeled solid |
| US6295464B1 (en) * | 1995-06-16 | 2001-09-25 | Dimitri Metaxas | Apparatus and method for dynamic modeling of an object |
| US6106466A (en) * | 1997-04-24 | 2000-08-22 | University Of Washington | Automated delineation of heart contours from images using reconstruction-based modeling |
| US6271861B1 (en) * | 1998-04-07 | 2001-08-07 | Adobe Systems Incorporated | Smooth shading of an object |
-
2002
- 2002-03-19 EP EP02713109A patent/EP1374181A1/en not_active Withdrawn
- 2002-03-19 JP JP2002578256A patent/JP4170096B2/ja not_active Expired - Fee Related
- 2002-03-19 WO PCT/IB2002/000894 patent/WO2002080110A1/en not_active Ceased
- 2002-03-28 US US10/109,371 patent/US20020172406A1/en not_active Abandoned
Also Published As
| Publication number | Publication date |
|---|---|
| EP1374181A1 (en) | 2004-01-02 |
| JP2004521425A (ja) | 2004-07-15 |
| US20020172406A1 (en) | 2002-11-21 |
| WO2002080110A1 (en) | 2002-10-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4170096B2 (ja) | 対象の3次元表面上にマップされた3次元メッシュモデルの適合性評価のための画像処理装置 | |
| US7015907B2 (en) | Segmentation of 3D medical structures using robust ray propagation | |
| EP1565880B1 (en) | Image processing system for automatic adaptation of a 3-d mesh model onto a 3-d surface of an object | |
| Tang et al. | Curvature-augmented tensor voting for shape inference from noisy 3d data | |
| EP1514229B1 (en) | Statistical model | |
| US8031919B2 (en) | Computerised cortex boundary extraction from MR images | |
| US7990379B2 (en) | System and method for coronary segmentation and visualization | |
| EP0526881B1 (en) | Three-dimensional model processing method, and apparatus therefor | |
| US8384716B2 (en) | Image processing method | |
| US7043062B2 (en) | Image processing method for displaying an image sequence of a deformable 3-D object with indications of the object wall motion | |
| JP2008511366A (ja) | 距離座標を使った特徴重み付け医療オブジェクト輪郭検出 | |
| CN107708550A (zh) | 用于检测和测量解剖异常的分段回声结构的表面建模 | |
| US20140161352A1 (en) | Iterative method for determining a two-dimensional or three-dimensional image on the basis of signals arising from x-ray tomography | |
| Ballester et al. | Segmentation and measurement of brain structures in MRI including confidence bounds | |
| Duan et al. | A subdivision-based deformable model for surface reconstruction of unknown topology | |
| CN111369662B (zh) | Ct图像中血管的三维模型重建方法及系统 | |
| Brett et al. | A Method of 3D Surface Correspondence for Automated Landmark Generation. | |
| Djuricic et al. | 3d central line extraction of fossil oyster shells | |
| JP2003123057A (ja) | 3次元形状モデルの生成方法および装置 | |
| He et al. | Brain image analysis using spherical splines | |
| EP4521355A1 (en) | Correcting topological defects on a surface mesh representing an organ | |
| Karbowski | Edge detection in reverse engineering system | |
| Fripp et al. | A subdivision-based parametric deformable model for surface extraction and statistical shape modeling of the knee cartilages | |
| Ličev et al. | Search and Implementation of Optimization Algorithms in Analysis of Ultrasonic Pictures in Neurology | |
| Patten | Three-dimensional modeling of point clouds |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050316 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20071211 |
|
| A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20080310 |
|
| A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20080317 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080611 |
|
| 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: 20080708 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080806 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110815 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |