Yamyam: 1版をインポートしました

2018-06-28T17:45:29Z

1版をインポートしました

2010年4月3日 (土) 12:34にwiki>Terrywallworkによる

2010-04-03T12:34:51Z

新規ページ

=MetaMesh=

MetaMesh project is part of my thesis. Purpose of this project is to create simple and easy way of [[BlenderDev/MetaBall|MetaBall]] editting. Non-manifold Mesh will be directing structure for MetaBalls.

[[Image:meta-mesh.jpg]]

== Data structure==

*visualisation of [[BlenderDev/MetaVert|MetaVert]] data structure: [[Image:MetaVert.png|visualisation of MetaVert data structure]]

*visualisation of [[BlenderDev/MetaEdge|MetaEdge]] data structure: [[Image:MetaEdge.png|visualisation of MetaEdge data structure]]

*visualisation of [[BlenderDev/MetaFace|MetaFace]] data strucure: [[Image:MetaFace.jpg|visualisation_of_MetaFace_data_structure]]

== Theory of MetaBalls==

=== Implicit Function===

Polygonised surface of [[BlenderDev/MetaBalls|MetaBalls]] is just set of points in space with the same value of implicit function:

[[Image:implicit_function.png|implicit function]]

where F(P) is value of implicit funciton. Points at polygonised surface has zero value. c parameter is "stiffness" of each MetaElem. F blending function and r is distance from each MetaElem.

=== Blending Function===

We use Blending function with this form:

{|
|valign=top|[[image:Dev-blending_function.png|blending_function]]
|}

r parameter is distance from [[BlenderDev/MetaElem|MetaElem]] and R is radius of MetaElem. [[BlenderDev/MetaElem|MetaElem]] behind this radius has no influence at polygonised surface.

=== Computation of distance from directing structures===

==== Point (Vertex)====

Distance from point is very easy to computate:

{|
|valign=top|[[image:Dev-point_length.png|distance_from_point]]
|}

==== Segment (Edge)====

Computation of distance from segment isn't so easy as computation distance from a point. I will show you how hard it is if we computate it 3D space. I don't do it so stupid, of course ;-)

Segment is part of the line. Line is identified by two points A[a0, a1, a2], B[b0, b1, b2] in space. Parametric simultaneous sequations of line has these form:

{|
|valign=top|[[image:Dev-line_equations.png|line_equations]]
|}

At the first time we have to computate parameter tI of rectangular projection of point P[p0, p1, p2] to the line.

{|
|valign=top|[[image:Dev-rectangular_projection.png|rectangular_projection]]
|}

If parameter tI is inbetween <0,1> then we can computate position C[c0, c1, c2] of this projection using previous simultaneous equations. If parameter tI is bigger then 1 or smaller then 0 value, then it means, then rectangular projection doesn't intersects segment.

If rectangular projection intersects segment then distance from this segment is:

{|
|valign=top|[[image:Dev-two_points_length.png|two_points_length]]
|}

It rectangular projection doesn't intersect segment, then we have to computate distance from point P[p0, p1, p2] to the both end points of segment. Minimal distance is distance from segment.

Computation distance from segment using these math formulas isn't too much smart. We use different way. We pre-computate inverse transformation matrix for each segment(edge), which "transfer" an segment to the origin of coordinates. Computation of distance from such segment is much more simple.

Ilustration of transformation in 2D:

{|
|valign=top|[[image:Dev-metaball-segment-distance.jpg|metaball-segment-distance]]
|}

We don't have transform segment with inverse trnasform matrix, because we know, that transformed segment will lay in the origin of coordinates. We have to transform only point P=X[x,y,z,w] w (weight of point is 1):

{|
|valign=top|[[image:Dev-transform_matrix.jpg|transform_matrix]]
|}

This transformation could be computate using modern instruction set (MMX/SSE) of our processors.

Computation of distance transformed point from segment laing in th origin of coordinates is very easy then:

<source lang="cpp">
/* simplificated densfunc() in mball.c
* vec ... coordinates of point X
* ball ... segment primitive */
float densfunc(MetaElem *ball, float *vec)
{
float dist2;

Mat4MulVecfl(ball->imat, vec);
if(ball->type==MB_EDGE) {
if( vec[0] > ml->expx) vec[0]-= ml->expx;
else if(vec[0]< -ml->expx) vec[0]+= mll->expx;
else vec[0]= 0.0;
}
/* only powers of distance squared are used in blending function */
dist2= (vec[0]*vec[0] + vec[1]*vec[1] + vec[2]*vec[2]);
/* calculation of blending function */
dist2= 1.0f-(dist2/ball->rad2);
if(dist2 < 0.0) return 0.5f;
return 0.5f-ball->s*dist2*dist2*dist2;
}</source>

==== Triangle (Face)====

I captured several demonstraiting videos of my devel version:

*http://e-learning.vslib.cz/~hnidek/captured-videos/meta-mesh.avi
*http://e-learning.vslib.cz/~hnidek/captured-videos/point-dist.avi
*http://e-learning.vslib.cz/~hnidek/captured-videos/meta-edge-ellipsoid.avi

[[User:JiriHnidek|Jiří Hnídek]] 01 Nov 2004

[[Category:Script]]

← 古い版	2018年6月28日 (木) 17:45時点における版
(相違点なし)

Dev:Source/Modeling/MetaMesh - 版の履歴

Yamyam: 1版 をインポートしました

2010年4月3日 (土) 12:34にwiki>Terrywallworkによる

Yamyam: 1版をインポートしました