JP2004133550A - 非多様体の陰関数表現方法と陰関数曲面の直接描画方法及びそれらのプログラム - Google Patents
非多様体の陰関数表現方法と陰関数曲面の直接描画方法及びそれらのプログラム Download PDFInfo
- Publication number
- JP2004133550A JP2004133550A JP2002295332A JP2002295332A JP2004133550A JP 2004133550 A JP2004133550 A JP 2004133550A JP 2002295332 A JP2002295332 A JP 2002295332A JP 2002295332 A JP2002295332 A JP 2002295332A JP 2004133550 A JP2004133550 A JP 2004133550A
- Authority
- JP
- Japan
- Prior art keywords
- implicit function
- manifold
- volume
- implicit
- value
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/30—Polynomial surface description
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Computer Graphics (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Image Generation (AREA)
Abstract
【解決手段】非多様体の陰関数場をボリュームデータの形で保持し、格子点の間の点における陰関数の値を補間により決定し、かつ補間の対象となる2つの互いに隣接するボクセルでの符号付距離の差が一定幅より大きい場合には、このボクセル間には表面を生成しない。また入力曲面を表裏の判別が出来る曲面パッチに分解し、その表裏にそれぞれ番号を振ってこれらを区別し、最近傍点の面の番号を用いて空間を複数領域に分類する。
【選択図】 図3
Description
【発明の属する技術分野】
本発明は、非多様体の陰関数表現方法と陰関数曲面の直接描画方法及びそれらのプログラムに関する。
【0002】
【従来の技術】
計算機で3次元物体の表面形状を表現する際には、物体の表面の座標値で形状を定義する、パラメトリック表現が広く用いられている。この例として、CAD(Computer−Aided Design)で利用される自由曲面表現、仮想現実感モデルや娯楽の分野で利用されるポリゴン表現などが挙げられる。パラメトリック表現はデータ構造が直感的で分かりやすく、比較的少ないデータ量で形状を表現でき、また高速な描画ハードウェアの利用が容易であるという特長がある。
【0003】
一方、表面形状の変形や集合演算といった処理を行う場合には、陰形式の関数を用いて間接的に曲面を表現する、陰関数表現を用いるほうが便利な場合が多く、モデリングや形状処理、物理シミュレーション等の分野で使われている。この表現形式は、データ構造が単純なため複雑な処理を容易に行うことができ、処理の階層化や並列化に向いているという特長がある。曲面のパラメトリック表現と陰関数表現にはそれぞれ利点と欠点があり、処理に応じて形式を選択できることが望ましく、その際には両者の形式を変換する必要が生じる。
【0004】
なお、本発明に関連する先行技術文献として、以下のものが挙げられる。
【0005】
【非特許文献1】
J. Bloomenthal, Polygonization of implicit surfaces, Computer Aided Geometric Design, 5:341−355,1988.
【非特許文献2】
J. Bloomenthal, Introduction to Implicit Surface, Morgan Kaufmann Publishers, Inc., 1997.
【非特許文献3】
J. Bloomenthal and K. Ferguson, Polygonization of non−manifold implicitsurfaces, In SIGGRAPH ’95, pages 309−316, 1995.
【非特許文献4】
M. Brady, K. Jung, H. T. Nguyen and T Nguyen, Two−phase perspective raycasting for interactive volume navigation, In Visualization 97, pages 243−56,1997.
【非特許文献5】
K. Engel, M. Kraus and T. Ertl, High−quality pre−integrated volume rendering using hardware−accelerated pixel shading, In Eurographics/SIGGRAPH
Workshop on Graphics Hardware ’01,
pages 9−16,2001.
【非特許文献6】
H. C. Hege, M. Seebas, D. Stalling and M. Zockler, A generalized marching cubes algorithm based on non−binary classifications, Technical report, Konrad−Zuse−Zentrum fur Information stechnik(ZIB),1997.
【非特許文献7】
Philippe Lacroute and Marc Levoy, Fast volume rendering using a shear−warp factorization of the viewing transformation, In SIGGRAPH ’94, pages
451−458, 1994.
【非特許文献8】
W. E. Lorensen and H. E. Cline, Marching cubes: a high resolution 3d surface reconstruction algorithm, In SIGGRAPH ’87, pages 163−169,1987.
【非特許文献9】
J. Rossignac and M. O’Connor, SGC:A dimension independent model for pointsets with internal structures and incomplete boundaries, In GeometricModeling for Product Engineering, pages 145−180,1990.
【非特許文献10】
H. Tuy and L. Tuy, Direct 2d display of 3d objects, IEEE Mag, ComputerGraphics and Applications,4(10):29−33,1984.
【非特許文献11】
A. P. Witkin and P. S. Heckbert, Using particles to sample and control implicit surfaces, In SIGGRAPH ’94, pages 269−278,1994.
【非特許文献12】
Shuntaro Yamazaki, Kiwamu Kase and Katsushi Ikeuchi, Nonmanifold implicit surfaces based on discontinuous implicitization and polygonization, In Geometric Modeling and Processing, pages 138−146. IEEE, July 2002.
【0006】
【発明が解決しようとする課題】
陰関数表現をパラメトリック表現に変換する方法に関しては多くの方法が提案されており、特に三角形メッシュへの変換は安定にかつ高速に行うことが可能である。ところが、パラメトリック表現を陰関数表現に変換する方法に関しては必ずしも正しい結果が得られるとは限らない。これは従来の陰関数表現では、関数場の定義に連続な実数値関数を用いているため扱える表面の種類が2−多様体に限られており、入力に非多様体的な特徴がある場合、例えば、曲面に境界がある場合と表面に分岐がある場合には、それらが失われて形状が大きく変わってしまう問題点があった。
【0007】
一方、陰関数表現はそれ自体が曲面を表すのに十分な情報を持っているが、データとして保持する値が曲面の幾何形状と直接結びついていないため、実際の形状を確認するためには表面を可視化する必要がある。しかし、陰関数曲面を用いて形状設計を行う際には、表面の混合や変形などの処理の過程で逐次更新される変形などの処理によって関数場が変更されるに再メッシュ化を行う必要があり、処理に時間がかかる。また、表面形状が複雑になるほど生成されるポリゴン数が増え、描画速度が遅くなるという問題があった。
【0008】
本発明は、上述した問題点を解決し要望を満たすために創案されたものである。すなわち本発明の第1の目的は、入力に非多様体的な特徴がある場合、例えば、曲面に境界や分岐がある場合でも、非多様体曲面をパラメトリック表現から陰関数表現へ形状を正しく保持しながら変換する非多様体の陰関数表現方法とそのプログラムを提案することにある。
また本発明の第2の目的は、変形などによって陰関数曲面が更新された場合にも高速に処理でき、形状が複雑な場合でも高速描画ができ、これにより、陰関数場を対話的な速度で表示できる陰関数曲面の直接描画方法とそのプログラムを提案することにある。
【0009】
【課題を解決するための手段】
本発明によれば、非多様体の陰関数場をボリュームデータの形で保持し、格子点の間の点における陰関数の値を補間により決定し、かつ補間の対象となる2つの互いに隣接するボクセルでの符号付距離の差が一定幅より大きい場合には、このボクセル間には表面を生成しない、ことを特徴とする非多様体の陰関数表現方法が提案される。
この方法により、曲面に境界がある場合に、補間によって不要な面が生成されるのを防止することができる。
【0010】
本発明の好ましい実施形態によれば、任意の2つのサンプル点の間隔をw、値をそれぞれu、v(u≦v)とするとき、
u∈(−∞,t)・・・(2)
v∈[t,∞)・・・(3)
0<((−u)−t)+(v−t)<αw・・・(4)
ただしα(≧1)、
がすべて成立する場合にのみ、これら2つの点の間に表面が存在するものとする。
この方法により、補間の対象となる2つの互いに隣接するボクセルでの符号付距離の差が一定幅より大きい場合には、このボクセル間に表面を生成しないようにすることができる。
【0011】
また、q(∈[0,1])は表面の位置が、q=0の時に値がuの格子点上に、q=1の時に値がvの格子点上になるように正規化し、表面が存在する位置qを、
q=(t−u)/(v−u)・・・(5)
で求める。
この方法により、曲面に境界がある場合でも、表面の位置を誤りなく求めることができる。
【0012】
また本発明によれば、入力曲面を表裏の判別が出来る曲面パッチに分解し、その表裏にそれぞれ番号を振ってこれらを区別し、最近傍点の面の番号を用いて空間を複数領域に分類する、ことを特徴とする非多様体の陰関数表現方法が提供される。
この方法により、表面に分岐がある場合でも、分岐部分において接続している表面の一部が除去されてしまうのを防ぐことができる。
【0013】
本発明の好ましい実施形態によれば、1.入力の非多様体曲面を分岐線上に沿って分割し、それぞれ分岐を持たない曲面パッチに分解し、
2.得られたパッチに順に番号iを振り、パッチの表裏を区別し、表面の番号をi+、裏面の番号をi−とし、
3.空間を格子点pでサンプリングし、格子点に曲面へのEuclid距離dE(p)と最近傍点の面の番号i(p)を割り当て、
4.各格子点pに対して、6隣接点pnでi(pn)を調べ、i(p)≠i(pn) であるような(i(p),i(pn))の組を列挙し、
5.上で作られた番号の組を別の新たな番号で置換し、ただし置換の結果、最初にi+とi−であった番号が同じ番号になる場合にはその組み合わせに対する置換を行わず、これにより、最終的に番号が0から順に並ぶようにし、
6.置換表に従い各格子点pでの領域番号i(p)を書き換え、
得られたボリュームの領域番号i(p)と、各ボクセルにおける表面へのEuclid距離dE(p)から、実数値の陰関数ボリュームを構成する。
【0014】
この方法により、最近傍点の面の番号を用いて空間を複数領域に分類し、分岐のある曲面を陰関数表現することが可能になる。
【0015】
また、距離を表す実数空間を領域分割し、分割された領域の幅をDsとする時、領域iに含まれる距離di sは
di s∈[Dsi,Ds(i+1))・・・(6)
であり、各ボクセルの位置pで、dE(p)とi(p)から領域付距離fs(p)を次の式で計算し、
fs(p)=min(dE、2B−ε)+2Bi(p)・・・(7)
ε(>0)はfs(p)が式(6)の半開区間に含まれるようにdE(p)を切り捨てるための微小な正の実数に設定する。
この方法により、曲面に分岐がある場合でも、表面の位置を誤りなく求めることができる。
【0016】
また、任意の2つのサンプル点の間隔をw、値をそれぞれu、v(u≦v)とするとき、
u∈(2Bi,2B(i+1))・・・(8)
v∈[2Bj,2B(j+1))・・・(9)
0<(u−2Bi)+(v−2Bj)<αw・・・(10)
ただしi、j(0≦i≦j≦n−1)、α(≧1)、
がすべて成立する場合にのみ、これら2つの点の間に表面が存在するものとする。
この方法により、補間の対象となる2つの互いに隣接するボクセルでの領域付距離関数の差が一定幅より大きい場合には、このボクセル間に表面を生成しないようにすることができる。
【0017】
また、q(∈[0,1])は表面の位置が、q=0の時に値がuの格子点上に、q=1の時に値がvの格子点上になるように正規化し、表面が存在する位置qを、
q=(u−2Bi)/((u−2Bi)+(v−2Bj))・・・(11)
で求める。
この方法により、曲面に境界がある場合でも、表面の位置を誤りなく求めることができる。
【0018】
さらに本発明によれば、スライス前面のボリュームの値を表すテクスチャTfrontとスライス後面のボリュームの値を表すテクスチャTbackを用いてスライス前面とスライス後面で囲まれる領域のボリュームの値を補間して表示する、ことを特徴とする陰関数曲面の直接描画方法が提供される。
この方法により、スライスが存在しない領域も含めてボリューム全体の影響を描画結果に反映させることができ、モアレを除去することができる。
【0019】
また、本発明の好ましい実施形態によれば、視線とスライス前面及びスライス後面の交点を計算して、スライス前面のテクセル値tfrontとスライス後面のテクセル値tbackとから、スライス間で視線上にあるボリュームがこの位置で観測される色と透明度に与える影響を計算し、ポリゴン上に表示する。
この方法により、少ない枚数のスライスを用いた描画でも正しい結果が得られる。
【0020】
また、テクセル値tfront、tbackの組から、観測される色と透明度を計算する処理を前もって行っておき、結果を早見表として2次元テクスチャとしてグラフィクスカード内に保存し、描画時にテクスチャ合成機能を利用して参照する。
この方法により、計算量を減らし、高速に描画処理できる。
【0021】
また、領域付距離場ボリュームで表された陰関数曲面を3Dテクスチャとして使える形に変換し、テクセル値tfront、tbackからなる任意の領域付距離の組に対して、それらの間の観測される色と透明度を計算する処理を前もって行って早見表を作成し、これにより描画色を決定する。
この方法により、正しい表示結果を、高速に描画処理できる。
【0022】
更に本発明によれば、計算機に前記請求項1乃至13の方法を実行させる、ことを特徴とするコンピュータプログラムが提供される。
このコンピュータプログラムを用いることにより、変形などによって陰関数曲面が更新された場合にも高速に処理でき、形状が複雑な場合でも高速描画ができ、これにより、陰関数場を対話的な速度で表示できる。
【0023】
【発明の実施の形態】
以下、本発明の好ましい実施形態を図面を参照して説明する。
1. 本発明は、領域付距離関数を用いて陰関数ボリュームを生成し、この値を不連続性関数を用いて補間するように拡張することで、非多様体のパラメトリック曲面を陰関数表現する方法を提案する。
また、本発明の方法で定義された陰関数場に対してボリュームレンダリングの方法を応用し、格子点サンプリングされた陰関数場から曲面を直接描画する方法を提案する。この方法を用いると、多様体または非多様体の陰関数曲面を、表面形状の複雑さに依存することなく高速に描画できる。
また具体的なアルゴリズムとして、三角形メッシュを入力とし、陰関数場を格子点サンプリングして得られるボリュームデータを出力する方法について説明する。
この方法はグラフィクスカードのプログラム可能なテクスチャ合成機能を利用して効率的に実装することができ、 対話的な速度で高画質の曲面を表示できる。またその速度は表面の形状によらず一定であり、多様体、及び非多様体曲面を扱うことが出来る。
【0024】
まとめると、本発明の主な特徴は次の2点である。
(1)パラメトリック曲面を陰関数表現に変換する方法を提案する。特に、任意の三角形メッシュからボリューム化された陰関数場を生成し、そこから陰関数曲面を抽出するアルゴリズムを提案する。
(2)多様体、非多様体の陰関数曲面を、ボリュームレンダリングの方法を応用して描画する方法を提案する。
【0025】
2.関連研究
2.1 非多様体曲面の陰関数表現
陰関数表現はもともと多様体曲面のみを扱うことを主眼としているが、Rossignacらはベクトル値関数を用いることで非多様体曲面を陰関数表現することが可能であることを指摘した[非特許文献9]。これ基づきBloomenthalらは曲面を陰関数曲面として表現する方法を提案した[非特許文献3]。
しかし、一般にパラメトリック表現で与えられる曲面からこの方法が扱える形式の陰関数場を構成する方法は明らかではない。本発明では非多様体曲面の定義を明確にし、この定義に含まれる曲面が十分に実用的であることを述べた後、この定義に含まれる任意のパラメトリック曲面を陰関数表現に変換する方法を提案する。
【0026】
2.2 陰関数曲面の可視化
陰関数曲面は光線追跡法(ray tracing)による可視化と相性が良く、この方法を用いると高い画質で表面を可視化することが可能である[非特許文献2]。しかしこの方法は、処理に時間がかかるため対話的な速度で描画することは難しい。
表示を簡略化することにより描画速度を大幅に向上する方法として、Witkinらはparticle systemを用いた表示方法を提案している[非特許文献11]。しかし、この方法は陰関数が曲面上で微分可能であること仮定しているため、非多様体へ拡張された陰関数曲面に対しては適用できない。
陰関数曲面を三角形メッシュに変換することにより、描画処理を高速化することが可能である。Bloomenthalらは曲面を表す陰関数場を四面体分割し、各四面体ごとに三角形メッシュ化する方法を提案しており[非特許文献1]、この方法は一部の非多様体形状に対して拡張されている[非特許文献3]。
一方Lorensenらは陰関数場の六面体分割を用いて陰関数場を三角形メッシュ化する方法を提案しており[非特許文献8]、後に非多様体形状に拡張されている[非特許文献6]、[非特許文献12]。
これらの方法では、描画処理の前に関数場からメッシュを生成する必要があるため、変形などの処理によって関数場が変更されるに再メッシュ化を行う必要があり、処理に時間がかかる。また、表面形状が複雑になるほど生成されるポリゴン数が増え、描画速度が遅くなるという問題がある。
【0027】
本発明の方法では陰関数場から直接、表面を可視化するため、変形などによって陰関数場が更新された場合にも高速に処理できる。また形状の複雑度に無関係に描画することが出来るため、描画速度は一定である。
【0028】
3.非多様体曲面の陰関数表現
3.1 陰関数曲面
3次元空間内の曲面S⊂R3に対して、実数tと実数値関数f(p)の組で次を満たすものを考える。
p∈S ⇔ f(p)=t・・・(1)
ただしp∈R3は3次元空間内の点である。fはSへの距離に応じた値を返す関数と考えることができ、Sはfの値がtの等値面である。曲面Sの式(1)右辺の形式による定義を「陰関数表現」と呼び、このときSを「陰関数曲面」と呼ぶ。
【0029】
実数値関数fで定義される関数場を計算機で処理する時には、fをサンプリングする必要がある。そこで、あらかじめf格子点でサンプリングしてボリュームデータの形で保持することにより、陰関数場に対する処理を効率的に行うことが出来る。本発明ではこのボリュームデータを「陰関数ボリューム」と呼び、陰関数場はボリュームで定義されるものとする。陰関数場がボリュームで与えられる場合、格子点の間の点における陰関数の値はサンプル点での値を補間して決定する。即ち、表面での関数値tと陰関数ボリュームで陰関数曲面を定義するためには、陰関数場を定義する距離関数と、サンプル点から実際に表面の位置を決定する補間関数の両者が必要である。
【0030】
ある曲面の陰関数表現の定義は一つではないが、サンプリング幅に対して曲面の形状が複雑な場合にも安定して陰関数場を構成するために、符号付距離関数を用いて陰関数場を生成することが多い。この関数は、曲面へのEuclid距離に、最近傍点が面の表側である場合には正の符号を、裏面にある場合には負の符号を付けた値を返す。この関数を用いて生成された陰関数ボリュームでは、格子点間の値は周辺の格子点での値から線形補間して得ることができる。この様子を2次元で模式化した例を図1及び図9に示す。図中では格子の辺上でのみ関数値を補間して表面の位置を決定し、それらを繋ぐことにより表面を再構成している。
【0031】
3.2 非多様体曲面
3次元空間で、「集合Sが2−多様体である」とは、S上の任意の点における無限に小さい近傍が2次元の円盤と位相同相であることをいう。以降では、「多様体」という言葉を用いた場合、2−多様体を指すものとする。陰関数場が、サンプリングされた値の線形補間で与えられる場合には、空間内の任意の点で関数値が微分可能である。従って等値面上の任意の点pの近傍は法線が∇fの微小平面近似することができ、このとき陰関数曲面は多様体である。
一方、「Sが非多様体である」とは、表面上のある点の近傍領域において、座標軸が1つ以下しかとれないか、もしくは2次元局所座標系のとり方が3通り以上あることをいう。
【0032】
本発明で扱う形状はパラメトリック曲面、特に三角形メッシュで表現できる曲面であり、個々の三角形は自己交差を考えない限り、境界以外は多様体である。従って曲面上で多様体の性質を満たさなくなる点が存在するのは「三角形の境界のみ」であり、これは次の2つの場合に分類できる。
【0033】
(1)表面に境界がある場合:ある三角形の境界上の点が、別のどの三角形の境界上にもない場合。2次元で模式化された例を図2Aに示す。
(2)表面に分岐がある場合:ある三角形の境界上の点が、別の2つ以上の異なる三角形の境界上にある場合。2次元で模式化された例を図2Bに示す。
ただし入力の三角形の集合内の異なる2つの三角形が、その境界以外の点を共有する場合には、共有している点を含む直線に沿ってこれらの三角形を分割することにより、三角形同士の交差は生じないものとする。
図2にある通り、従来の方法を用いてこれらの非多様体形状を陰関数曲面に変換すると形状が大きく変わってしまう。そこで、非多様体曲面を正しく陰関数表現するためには、これらの特徴を陰関数場の中で保持し、また表面として再構成するアルゴリズムを構築する必要がある。
【0034】
3.3 距離関数と補間関数の設計
3.1節で述べたように、陰関数ボリュームを用いて曲面を陰関数表現するためには、陰関数場を構成するための距離関数と、表面の位置を決定するための補間関数を定義する必要がある。そこでこれらをうまく定めることにより、3.2節で挙げた2つの非多様体の特長を陰関数曲面の枠組みの中で表現する方法を考える。
【0035】
まず曲面に境界がある場合に関して、従来の陰関数表現では境界部分が延長されて、もともと存在しなかった表面が過剰に生成されてしまう(図2A)。これは、従来の陰関数を用いて定義された陰関数ボリューム中にも正しい表面の位置が定義されているが、補間によって不要な面が生成されていることを意味する。そこで元の曲面が存在しない位置には表面を生成しないように補間関数を拡張することで、表面の境界を表現することが可能である(図3A、図10A、図11)。
【0036】
距離関数として符号付距離を用いた場合には、補間関数として線形補間を用いて表面の位置を決定していた。そこでこれを拡張し、補間の対象となる2つの互いに隣接するボクセルでの符号付距離の差が一定幅より大きい場合には、このボクセル間には表面を生成しないようにする。
【0037】
即ち任意の2つのサンプル点の間隔をw、値をそれぞれu、v(u≦v)とするとき、これら2つの点の間に表面が存在するのは
u∈(−∞,t)・・・(2)
v∈[t,∞)・・・(3)
0<((−u)−t)+(v−t)<αw・・・(4)
が成立する場合のみとする。ただしα(≧1)は曲面の陰関数場への変換が誤差を含む場合にも表面を正しく発生させるために導入する係数である。
【0038】
陰関数場が格子点サンプリングされた値を線形補間した値で与えられる場合、空間内の点における関数値は最大格子点間隔の半分程度の誤差を含むと考えられるため、実験ではα=1.5を使用した。
【0039】
上の3つの式を満たすとき、表面が存在する位置は
q=(t−u)/(v−u)・・・(5)
で求められる。ただしq(∈[0,1])は表面の位置が、q=0の時に値がuの格子点上に、q=1の時に値がvの格子点上になるように正規化されている。
【0040】
次に、表面に分岐がある場合に関して、従来の陰関数表現では分岐部分において接続している表面の一部が除去されてしまう(図2B)。これは、連続な実数値関数を用いて定義された陰関数ボリューム中では分岐を持つ表面の位置を正しく表現できないことが原因である。
【0041】
曲面が多様体の場合、曲面上の任意の点の近傍において面の表裏を定義することが出来る。したがって表面上の任意の点は、最近傍点がその表側にある領域と裏側にある領域の2つの領域に挟まれる部分に存在し、このような領域分割を符号付距離関数を用いて表現することが出来る。ところが表面に分岐があり非多様体になっている場合には面の裏と表を定義することは出来ないため、特に分岐線上の点の周辺領域を最近傍点の面の方向から2種類に分類することはできない。
【0042】
入力の曲面を表裏の判別が出来る曲面パッチに分解しその表裏にそれぞれ番号を振ってこれらを区別することにより、最近傍点の面の番号を用いて空間を複数領域に分類すると、分岐のある曲面を陰関数表現することが可能になる(図3B、図10B、図12)。以下にその領域分割のアルゴリズムを示す。
【0043】
1.入力の非多様体曲面を分岐線上に沿って分割し、それぞれ分岐を持たない曲面パッチに分解する。
2.得られたパッチに順に番号iを振る。パッチの表裏を区別し、表面の番号をi+、裏面の番号をi−とする。
3.空間を格子点pでサンプリングし、格子点に曲面へのEuclid距離dE(p)と最近傍点の面の番号i(p)を割り当てる。
4.各格子点pに対して、6隣接点pnでi(pn)を調べ、i(p)≠i(pn) であるような(i(p),i(pn))の組を列挙する。
5.上で作られた番号の組を別の新たな番号で置換する。ただし置換の結果、最初にi+とi−であった番号が同じ番号になる場合にはその組み合わせに対する置換を行わない。最終的に番号が0から順に並ぶようにする。
6.置換表に従い各格子点pでの領域番号i(p)を書き換える。
【0044】
このようにして得られたボリュームの領域番号i(p)と、各ボクセルにおける表面へのEuclid距離dE(p)から、実数値の陰関数ボリュームを構成する。そのために、符号付距離を拡張した、「領域付距離」を定義する。領域付距離における「領域」とは符号付距離における正と負の2種類の「符号」を任意の数に拡張したものであり、距離を表す実数空間を領域分割することにより実現する。
【0045】
分割された領域の幅をDsとする時、領域iに含まれる距離di sは
di s∈[Dsi,Ds(i+1))・・・(6)
である。
以上のような領域付距離は以下の手順で計算できる。
まずDs=2Bと置き、B=Bmax−[log2n]から領域の大きさを決定する。ただしnは領域の数である。Bmaxは十分に大きい数であり、計算機上で表現する場合には領域付距離を表す変数のビット数として与えることができる。
【0046】
このとき各ボクセルの位置pで、dE(p)とi(p)から領域付距離fs(p)を次の式で計算する。
fs(p)=min(dE,2B−ε)+2Bi(p)・・・(7)
ただしε(>0)はfs(p)が式(6)の半開区間に含まれるようにdE(p) を切り捨てるための微小な正の実数である。実験ではεを格子点間隔、Bmax=8とした。
【0047】
領域付距離関数は連続関数ではないため、その値を線形に補間することは出来ない(図13参照)。そこで、与えられた領域付距離場中の領域数をn、距離場の2つのサンプル点の間隔をw、値をそれぞれu、v(u≦v)とするとき、これら2つの点の間に表面が存在するのは
u∈(2Bi,2B(i+1))・・・(8)
v∈[2Bj,2B(j+1))・・・(9)
0<(u−2Bi)+(v−2Bj)<αw・・・(10)
が全て成立するi、j(0≦i≦j≦n−1)が存在する場合のみとする。α≧1 は式(5)で用いた係数である。このとき表面が存在する位置は
q=(u−2Bi)/((u−2Bi)+(v−2Bj))・・・(11)
で求められる。qの定義は式(5)と同様である。図14は上述した領域付距離関数とその条件を示す模式図である。
【0048】
4.陰関数曲面の可視化
4.1 スライスを用いたボリュームレンダリング法
陰関数ボリュームに対してray casting法[非特許文献10]によるボリュームレンダリングの方法を適用することで陰関数曲面を直接、描画することが可能である。ray casting法では高画質な描画が可能である。
【0049】
ray castingは任意の視線上での関数値のサンプリングが必要であるため比較的コストの高い計算であるが、Lacrouteらはボリュームをx、 y、z軸と垂直な3組のスライスの集合で表現して(図4A)、ボリュームのサンプリングをスライス平面上のみで行うことで計算効率を上げて描画を高速化する方法を提案した[非特許文献7]。
【0050】
この方法はポリゴンとテクスチャマップを用いて効率的に実装することが可能である[非特許文献4]。この方法では、ボリュームのスライスを透明度付き2Dテクスチャとしてポリゴンにマップし、得られた半透明ポリゴンの集合を視線に対して奥の方から順にαブレンドしながら表示する。この表示処理ではグラフィクスカードの描画加速を利用できるため、高速に描画できる。またグラフィクスカードを使う場合には、3Dテクスチャ機能を利用することで、スライスの方向を固定せずに常に視線方向と垂直なポリゴンを発生させて描画することで(図4B)、多少の速度の減少と引き換えにより高画質な描画を行うことも可能である。
【0051】
4.2 Pre−integration法[非特許文献5]
スライスを使ったボリュームレンダリングの大きな問題点は、ボリュームの値の補間が、スライス上では任意の解像度で行われるのに対し、スライスと垂直な方向では固定のスライス間隔でしか行われないため、補間精度の不一致によりモアレと呼ばれる縞状の不具合や亀裂が観測される点である。この問題は、発生させるスライスの枚数を、必要な補間精度に応じて変更することにより回避可能であるが、描画速度とスライス枚数は反比例の関係にあり、多くの場合モアレが気にならない程度にスライスを発生させることは現実的ではない。
【0052】
この問題に対してEngelらは、各スライスに、そのスライス上でのボリュームの値だけでなく、隣接するスライス上でのボリュームの値も同時にマップし、スライス間に存在するボリュームの影響も考慮することでモアレを除去する、pre−integration法を提案した。
【0053】
図5Aのようにボリュームをスライス集合で表現し、frontで示されるスライスSfrontをポリゴンで描画する状況を考える。従来のスライス法ではこのポリゴンに対し、スライス上でのボリュームの値を表すテクスチャTfrontをマップして表示を行うため、スライス間のボリュームの値は描画結果に影響しない。
【0054】
一方、pre−integration法ではTfrontと同時にbackで示される隣接スライスSback上のテクスチャTbackもマップし、2つのテクスチャを用いてSfrontとSbackで囲まれる領域のボリュームの値を補間して表示する。その結果、スライスが存在しない領域も含めてボリューム全体の影響を描画結果に反映させることが出来る。実際にある視線での見え方を描画する時には、視線とSfront、Sbackの交点を計算して、front側のテクセル値tfrontとback側のテクセル値tbackから、スライス間で視線上にあるボリュームがこの位置で観測される色と透明度に与える影響を計算して、ポリゴン上に表示する。
【0055】
例えば図5Bにあるように、スライスを用いたボリュームレンダリングの方法で、表面を表示する状況を考える。図の中央の視線での見え方を考えるとき、従来のスライスを用いた方法ではスライスが表面上にはないため表面が描画されない。一方、pre−integration法ではスライス間のボクセル値の変化を考慮することが出来るため、少ない枚数のスライスを用いた描画でも正しい結果が得られる。
【0056】
tfront、tbackの組から、観測される色と透明度を計算する処理を前もって行っておき、結果をpre−integration tableと呼ばれる早見表として保存しておくことで計算量を減らすことが可能である。またその表は2次元テクスチャとしてグラフィクスカード内に保存しておき、描画時にテクスチャ合成機能を利用して参照することにより高速に描画処理できる。
【0057】
これを実現するためのグラフィクスカードの設定では、使用するハードウェアはnVidia社のGeForce3 GPUで、このグラフィクスカードのテクスチャ合成機能であるTexture hader機能とRegister Combiners機能を、NVParseを使ってプログラムした。
【0058】
この設定では図4Bにあるスクリーンに垂直なスライスに3Dテクスチャマップを行って描画する。また、ボリュームはフォン・シェーディングによって照光処理される。そのため、ボリュームとして与えられるスカラー場の勾配をあらかじめ計算しておく必要がある。各ボクセルにおけるスカラー値sと曲面の勾配(ux、uy、uz)が与えられるとき、ここからピクセル値(s、uy、uz、ux)を持つ3Dテクスチャを生成し、これを読み込む。
【0059】
描画時にはマルチテクスチャ機能を利用し、3Dテクスチャをテクスチャユニット0と1に読み込み、ユニット0からtbackが、ユニット1からtfrontが取り出せるようにテクスチャ座標を与える。pre−integrationtableは2Dテクスチャで表現してtfront、tbackで与えられる任意のスカラー値の組に対する色をRGBチャネルに、透明度をAチャネルに割り当てる。
【0060】
ただし照光計算の際にスライス間で勾配を補間するために、Bチャネルに勾配を補間する係数を割り当て、色は実際にはRGの2チャネルのみで表現する。このテクスチャをユニット3に読み込むと、このユニットに対するテクスチャ座標がグラフィクスハードウェアによって(tfront、tback)に設定され、pre−integration tableが参照される。
【0061】
4.3 陰関数曲面の可視化方法
領域付距離場ボリュームで表された陰関数曲面をpre−integration法で描画するためには、陰関数ボリュームを3Dテクスチャとして使える形に変換し、tfront、tbackからなる任意の領域付距離の組に対して、それらの間の描画色を決定するためのpre−integration tableを作ればよい。
【0062】
まず陰関数ボリュームに関して、領域付距離は一般に実数で与えられるが、現行のグラフィクスカードでテクスチャとして使うために距離を8bit整数(0−255)で量子化する。またシェーディングに必要な勾配は、陰関数場の生成の際に元の曲面の勾配から計算し、同様に量子化しておく。
【0063】
pre−integration tableは256×256のサイズの2次元のRGBAテクスチャとして作成する。RGの2つのチャネルに割り当てる色は、領域付距離の距離の組ごとに任意の色を設定でき、表面を均一に描画する時には単一色を与える。Aチャネルには、スライス間に面が存在しない場合にはA=0(透明)、存在する場合にはA=1(不透明)とする。最後にBチャネルに割り当てる補間係数として式(11)から計算される、q∈[0,1]を8bit整数で量子化した値を用いる。
【0064】
例として図6Aに示す非多様体曲面を陰関数表現して可視化することを考える。このモデルは球面が平面にはめ込まれた形の非多様体形状であり、空間は3つの領域に分類できる。これに対応するpre−integrationtableは図6Bの通りである。表の横軸がtfront∈[0,255]、縦軸がtback∈[0,255]に対応し、表の各点の色が描画される表面の色を、色がない部分は表面が存在せずに透明色で描画される組み合わせに対応する。色の濃さは勾配の補間のための係数に対応している。なお、図15、16は別のpre−integrationtableを例示している。
【0065】
非多様体の陰関数曲面をpre−integration法で可視化する際には、テクスチャとして与えられるボクセル値の領域付距離が線形補間されないようにする必要がある。例えば図6の非多様体曲面を陰関数ボリュームに変換し、ボクセル値の線形補間を行って可視化した結果が図7Aである。線形補間によって、領域境界付近に本来存在しない領域に属する距離が生成され、本来存在しない表面が現れている。この不具合は、ボクセル値の線形補間を行わないことで解決することができ、結果として図7Bのように正しい表示結果を得ることが出来る。
【0066】
5.実験結果
上述した本発明の方法に基づき、三角形メッシュを陰関数ボリューム化するソフトウェア、及びそのボリュームを陰関数曲面表示するソフトウェアを実装し、複数のモデルに対して実験を行った。可視化の実験に利用した計算機はCPU Pentium4 1.7GHz、主記憶1.0GB、VPU GeForce4 Ti4600、VRAM 128MBのPCである。可視化のソフトウェアとして描画面に水平なスライスを用いてpre−integrationを行うボリュームレンダラを作成した。
【0067】
図8は非多様体表面モデルを陰関数表現に変換し、可視化した結果である。上段が入力として与えたパラメトリック表現された表面モデルであり、三角形メッシュを使用した。これらのモデルを領域付距離場を用いて2563サイズの陰関数ボリュームに変換し、可視化した結果をそれぞれ下段に示す。
【0068】
図8A,Bに示すモデルには部分的に穴が開いており、表面に境界が存在する非多様体である。従来の陰関数表現では扱うことの出来なかった境界のある曲面を、本発明の方法では表現できている。ただし陰関数場をボリュームでサンプリングしているため、ボクセル幅以下の形状の特徴を正確に表現することは出来ず、特に境界付近の形状は多少変化している。この問題はボリュームのサンプリング密度を高くすることで解消できる。
【0069】
図8Cは、平面と半球面からなる表面形状に対する実験結果である。この表面は境界線と表面の分岐をもつ非多様体であるが、本発明の方法では境界や分岐が正しく表現できることがわかる。距離関数の値と勾配を8ビットで量子化して最近傍点補間しているが、描画色は2つのテクセル値の組み合わせから決定するため、描画面が極端に大きくない限りは結果画像中の量子化ノイズはほとんど目立たない。描画時の色は、表面が存在する領域の対ごとに、表面と裏面を区別して設定することが可能である。この例では全ての組み合わせに異なる色を設定している。
図17は、本発明による可視化した図を示す図8と同様のディスプレイ上の画像である。
【0070】
また1283と2563サイズのボリュームデータに対し可視化処理を行った際の描画速度を表1に示す。発生させるスライス枚数を変更することによって速度は変化するが、多くの実験結果で対話的な速度に必要とされる30fpsを達成できた。
【0071】
【表1】
【0072】
上述したように、本発明では、陰関数曲面を定義する距離関数と補間関数として不連続関数を用いることにより、陰関数表現で非多様体曲面を扱えるように拡張した。またこのように定義された多様体・非多様体陰関数曲面を高速に描画する方法を提案した。
【0073】
実験で用いた入力のパラメトリック曲面は三角形メッシュのみであるが、本発明の方法は自由曲面にも適用可能である。ただし自由曲面を多様体パッチに分割するアルゴリズムなどの細かい変更は必要であり今後、実験する必要がある。また形状の変形や集合演算を行う際には曲面の陰関数表現が有効であるが、本発明の方法で定義された非多様体陰関数曲面表現に対してもこれらの処理が適用可能である。
【0074】
従来の陰関数曲面では取り扱うことの出来る形状に制限が多く、表面の境界と表面の分岐という大きな問題に関して、本発明の方法で解決することが可能である。
【0075】
また、上述した実施例に示すように、本発明の方法に基づき、三角形メッシュを陰関数ボリューム化するソフトウェア、及びそのボリュームを陰関数曲面表示するソフトウェアをコンピュータプログラムとして実装することにより、変形などによって陰関数曲面が更新された場合にも高速に処理でき、形状が複雑な場合でも高速描画ができ、これにより、陰関数場を対話的な速度で表示できる。
【0076】
【発明の効果】
上述したように、本発明の非多様体の陰関数表現方法と陰関数曲面の直接描画方法及びそれらのプログラムは、入力に非多様体的な特徴がある場合、例えば、曲面に境界や分岐がある場合でも、非多様体曲面をパラメトリック表現から陰関数表現へ形状を正しく保持しながら変換することができ、かつ、変形などによって陰関数曲面が更新された場合にも高速に処理でき、形状が複雑な場合でも高速描画ができ、これにより、陰関数場を対話的な速度で表示できる、等の優れた効果を有する。
【図面の簡単な説明】
【図1】符号付距離場と線形補間で再構築された表面を示す2次元模式図である。
【図2】従来法で陰関数表現された非多様体形状を示す2次元模式図である。
【図3】本発明の方法で陰関数表現された非多様体形状を示す2次元模式図である。
【図4】ボリュームの軸に沿ったスライス生成(A)と、描画面と垂直なスライス生成(B)の模式図である。
【図5】隣接する2つのスライス上のテクスチャを用いた描画の図である。
【図6】非多様体モデルに対する領域分け(A)と、対応するpre−integration table(B)のディスプレイ上の画像である。
【図7】領域付距離場の値を線形補間した場合の不具合(A)と、最近傍点をとることでその不具合を防いだ場合(B)を示すディスプレイ上の画像である。
【図8】入力の三角形メッシュモデル(上段)と、本発明による可視化した図(下段)を示すディスプレイ上の画像である。
【図9】符号付距離場と線形補間で再構築された表面を示す図1と同様の2次元模式図である。
【図10】本発明の方法で陰関数表現された非多様体形状を示す図3と同様の2次元模式図である。
【図11】表面に境界がある場合の距離関数とその条件を示す説明図である。
【図12】表面に分岐がある場合の条件を示す説明図である。
【図13】本発明による領域付距離関数を示す模式図である。
【図14】本発明による領域付距離関数とその条件を示す模式図である。
【図15】線形補間(A)と境界のある表面に対応するpre−integrationtable(B)のディスプレイ上の画像である。
【図16】線形補間(A)と分岐・境界のある表面に対応するpre−integration table(B)のディスプレイ上の画像である。
【図17】本発明による可視化した図を示す図8と同様のディスプレイ上の画像である。
Claims (13)
- 非多様体の陰関数場をボリュームデータの形で保持し、格子点の間の点における陰関数の値を補間により決定し、かつ補間の対象となる2つの互いに隣接するボクセルでの符号付距離の差が一定幅より大きい場合には、このボクセル間には表面を生成しない、ことを特徴とする非多様体の陰関数表現方法。
- 任意の2つのサンプル点の間隔をw、値をそれぞれu、v(u≦v)とするとき、
u∈(−∞,t)・・・(2)
v∈[t,∞)・・・(3)
0<((−u)−t)+(v−t)<αw・・・(4)
ただしα(≧1)、
がすべて成立する場合にのみ、これら2つの点の間に表面が存在するものとする、ことを特徴とする請求項1に記載の非多様体の陰関数表現方法。 - q(∈[0,1])は表面の位置が、q=0の時に値がuの格子点上に、q=1の時に値がvの格子点上になるように正規化し、表面が存在する位置qを、
q=(t−u)/(v−u)・・・(5)
で求める、ことを特徴とする請求項2に記載の非多様体の陰関数表現方法。 - 入力曲面を表裏の判別が出来る曲面パッチに分解し、その表裏にそれぞれ番号を振ってこれらを区別し、最近傍点の面の番号を用いて空間を複数領域に分類する、ことを特徴とする非多様体の陰関数表現方法。
- 1.入力の非多様体曲面を分岐線上に沿って分割し、それぞれ分岐を持たない曲面パッチに分解し、
2.得られたパッチに順に番号iを振り、パッチの表裏を区別し、表面の番号をi+、裏面の番号をi−とし、
3.空間を格子点pでサンプリングし、格子点に曲面へのEuclid距離dE(p)と最近傍点の面の番号i(p)を割り当て、
4.各格子点pに対して、6隣接点pnでi(pn)を調べ、i(p)≠i(pn) であるような(i(p),i(pn))の組を列挙し、
5.上で作られた番号の組を別の新たな番号で置換し、ただし置換の結果、最初にi+とi−であった番号が同じ番号になる場合にはその組み合わせに対する置換を行わず、これにより、最終的に番号が0から順に並ぶようにし、
6.置換表に従い各格子点pでの領域番号i(p)を書き換え、
得られたボリュームの領域番号i(p)と、各ボクセルにおける表面へのEuclid距離dE(p)から、実数値の陰関数ボリュームを構成する、ことを特徴とする請求項4に記載の非多様体の陰関数表現方法。 - 距離を表す実数空間を領域分割し、分割された領域の幅をDsとする時、領域iに含まれる距離di sは
di s∈[Dsi,Ds(i+1))・・・(6)
であり、各ボクセルの位置pで、dE(p)とi(p)から領域付距離fs(p)を次の式で計算し、
fs(p)=min(dE,2B−ε)+2Bi(p)・・・(7)
ε(>0)はfs(p)が式(6)の半開区間に含まれるようにdE(p)を切り捨てるための微小な正の実数に設定する、ことを特徴とする請求項4に記載の非多様体の陰関数表現方法。 - 任意の2つのサンプル点の間隔をw、値をそれぞれu、v(u≦v)とするとき、
u∈(2Bi,2B(i+1))・・・(8)
v∈[2Bj,2B(j+1))・・・(9)
0<(u−2Bi)+(v−2Bj)<αw・・・(10)
ただしi、j(0≦i≦j≦n−1)、α(≧1)、
がすべて成立する場合にのみ、これら2つの点の間に表面が存在するものとする、ことを特徴とする請求項4に記載の非多様体の陰関数表現方法。 - q(∈[0,1])は表面の位置が、q=0の時に値がuの格子点上に、q=1の時に値がvの格子点上になるように正規化し、表面が存在する位置qを、
q=(u−2Bi)/((u−2Bi)+(v−2Bj))・・・(11)
で求める、ことを特徴とする請求項7に記載の非多様体の陰関数表現方法。 - スライス前面のボリュームの値を表すテクスチャTfrontとスライス後面のボリュームの値を表すテクスチャTbackを用いてスライス前面とスライス後面で囲まれる領域のボリュームの値を補間して表示する、ことを特徴とする陰関数曲面の直接描画方法。
- 視線とスライス前面及びスライス後面の交点を計算して、スライス前面のテクセル値tfrontとスライス後面のテクセル値tbackとから、スライス間で視線上にあるボリュームがこの位置で観測される色と透明度に与える影響を計算し、ポリゴン上に表示する、ことを特徴とする請求項9に記載ぼ陰関数曲面の直接描画方法。
- テクセル値tfront、tbackの組から、観測される色と透明度を計算する処理を前もって行っておき、結果を早見表として2次元テクスチャとしてグラフィクスカード内に保存し、描画時にテクスチャ合成機能を利用して参照する、ことを特徴とする請求項9又は10に記載ぼ陰関数曲面の直接描画方法。
- 領域付距離場ボリュームで表された陰関数曲面を3Dテクスチャとして使える形に変換し、テクセル値tfront、tbackからなる任意の領域付距離の組に対して、それらの間の観測される色と透明度を計算する処理を前もって行って早見表を作成し、これにより描画色を決定する、ことを特徴とする請求項9又は10に記載ぼ陰関数曲面の直接描画方法。
- 計算機に前記請求項1乃至13の方法を実行させる、ことを特徴とするコンピュータプログラム。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002295332A JP4165695B2 (ja) | 2002-10-08 | 2002-10-08 | 非多様体の陰関数表現方法のコンピュータプログラム |
US10/615,782 US7158131B2 (en) | 2002-10-08 | 2003-07-10 | Implicit function rendering method of nonmanifold, direct drawing method of implicit function curved surface and programs thereof |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002295332A JP4165695B2 (ja) | 2002-10-08 | 2002-10-08 | 非多様体の陰関数表現方法のコンピュータプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004133550A true JP2004133550A (ja) | 2004-04-30 |
JP4165695B2 JP4165695B2 (ja) | 2008-10-15 |
Family
ID=32040761
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002295332A Expired - Fee Related JP4165695B2 (ja) | 2002-10-08 | 2002-10-08 | 非多様体の陰関数表現方法のコンピュータプログラム |
Country Status (2)
Country | Link |
---|---|
US (1) | US7158131B2 (ja) |
JP (1) | JP4165695B2 (ja) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006013813A1 (ja) * | 2004-08-02 | 2006-02-09 | Kyoto University | 情報処理装置およびプログラム |
JP2006048518A (ja) * | 2004-08-06 | 2006-02-16 | Toshiba Corp | 画像処理方法、画像処理プログラム、及び画像処理コンピュータ |
WO2007015365A1 (ja) * | 2005-08-01 | 2007-02-08 | National University Corporation NARA Institute of Science and Technology | 情報処理装置およびプログラム |
JP2007094628A (ja) * | 2005-09-28 | 2007-04-12 | Japan Science & Technology Agency | 粒子系を用いた陰関数曲面のポリゴン化方法 |
JP2014021659A (ja) * | 2012-07-17 | 2014-02-03 | Fujitsu Ltd | 描画処理プログラム、描画処理方法及び描画処理装置 |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ATE358358T1 (de) * | 2003-08-25 | 2007-04-15 | Sony Ericsson Mobile Comm Ab | Antennenschalterstruktur für ein mobiles endgerät in einem drahtlosen kommunikationssystem |
US7443393B2 (en) * | 2006-01-19 | 2008-10-28 | International Business Machines Corporation | Method, system, and program product for re-meshing of a three-dimensional input model using progressive implicit approximating levels |
GB0818279D0 (en) * | 2008-10-06 | 2008-11-12 | Advanced Risc Mach Ltd | Graphics processing systems |
US9123161B2 (en) | 2010-08-04 | 2015-09-01 | Exxonmobil Upstream Research Company | System and method for summarizing data on an unstructured grid |
US20140267266A1 (en) * | 2013-03-14 | 2014-09-18 | Nvidia Corporation | Generating anti-aliased voxel data |
US11024004B2 (en) * | 2018-08-31 | 2021-06-01 | International Business Machines Corporation | Reconstructing missing slices in slice-based scans of 3D objects |
CN110033519B (zh) * | 2019-04-23 | 2022-06-28 | 中南大学 | 基于隐式函数的三维建模方法、装置、系统及存储介质 |
CN116630568B (zh) * | 2023-07-24 | 2023-11-03 | 航天宏图信息技术股份有限公司 | 基于点、线、面结构的地质三维模型生成方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0681275B2 (ja) * | 1985-04-03 | 1994-10-12 | ソニー株式会社 | 画像変換装置 |
-
2002
- 2002-10-08 JP JP2002295332A patent/JP4165695B2/ja not_active Expired - Fee Related
-
2003
- 2003-07-10 US US10/615,782 patent/US7158131B2/en not_active Expired - Fee Related
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006013813A1 (ja) * | 2004-08-02 | 2006-02-09 | Kyoto University | 情報処理装置およびプログラム |
US8149237B2 (en) | 2004-08-02 | 2012-04-03 | National University Corporation NARA Institute of Science and Technology | Information processing apparatus and program |
JP2006048518A (ja) * | 2004-08-06 | 2006-02-16 | Toshiba Corp | 画像処理方法、画像処理プログラム、及び画像処理コンピュータ |
JP4557629B2 (ja) * | 2004-08-06 | 2010-10-06 | 株式会社東芝 | 画像処理方法、画像処理プログラム、及び画像処理コンピュータ |
WO2007015365A1 (ja) * | 2005-08-01 | 2007-02-08 | National University Corporation NARA Institute of Science and Technology | 情報処理装置およびプログラム |
US8149236B2 (en) | 2005-08-01 | 2012-04-03 | National University Corporation NARA Institute of Science and Technology | Information processing apparatus and program |
JP2007094628A (ja) * | 2005-09-28 | 2007-04-12 | Japan Science & Technology Agency | 粒子系を用いた陰関数曲面のポリゴン化方法 |
JP4699152B2 (ja) * | 2005-09-28 | 2011-06-08 | 独立行政法人科学技術振興機構 | 粒子系を用いた陰関数曲面のポリゴン化方法 |
JP2014021659A (ja) * | 2012-07-17 | 2014-02-03 | Fujitsu Ltd | 描画処理プログラム、描画処理方法及び描画処理装置 |
Also Published As
Publication number | Publication date |
---|---|
JP4165695B2 (ja) | 2008-10-15 |
US20040068530A1 (en) | 2004-04-08 |
US7158131B2 (en) | 2007-01-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9367943B2 (en) | Seamless fracture in a production pipeline | |
KR20080022551A (ko) | 프로시저 기하 개체의 삼각화 | |
JPH0950537A (ja) | ボリューム・レンダリング装置及び方法 | |
JP4165695B2 (ja) | 非多様体の陰関数表現方法のコンピュータプログラム | |
CN1655191A (zh) | 基于可编程图形硬件的多边形网格模型的快速体素化方法 | |
US20120038640A1 (en) | Spatial decomposition methods using bit manipulation | |
US7777740B2 (en) | Spatial decomposition methods using bit manipulation | |
US11955246B2 (en) | Unified anisotropic volume and surface mesh storage | |
Scholz et al. | Real‐time isosurface extraction with view‐dependent level of detail and applications | |
Lee et al. | Adaptive synthesis of distance fields | |
Zhou et al. | Pixel-exact rendering of spacetime finite element solutions | |
Priyadarshi et al. | Finding mold-piece regions using computer graphics hardware | |
Luo et al. | Gpu-based multi-slice per pass algorithm in interactive volume illumination rendering | |
Brasher et al. | Rendering planar cuts through quadratic and cubic finite elements | |
Shakaev et al. | View-Dependent Level of Detail for Real-Time Rendering of Large Isosurfaces | |
CN117152334B (zh) | 一种基于电波、气象云图大数据实现的三维仿真模拟方法 | |
Bolla | High quality rendering of large point-based surfaces | |
US11869123B2 (en) | Anti-aliasing two-dimensional vector graphics using a compressed vertex buffer | |
Sud | Efficient computation of discrete voronoi diagram and homotopy-preserving simplified medial axis of a three-dimensional polyhedron | |
Harbinson et al. | Hybrid polygon-point rendering of singular and non-manifold implicit surfaces | |
Zhou et al. | A Novel Pointcloud-based Isosurface Extraction | |
Heinonen | Simplification of 3D Computer-aided Design Models to Improve Rendering Performance | |
He et al. | A Rendering Method for Visualization of Medical Data | |
Zhang et al. | Graphics hardware based volumetric medical dataset visualization and classification | |
Johansson | Accelerating isosurface extraction by caching cell topology with graphics hardware |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
RD03 | Notification of appointment of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7423 Effective date: 20040322 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050927 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20080512 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080516 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080701 |
|
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: 20080723 |
|
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: 20080723 |
|
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: 20110808 Year of fee payment: 3 |
|
LAPS | Cancellation because of no payment of annual fees |