Approximate Implicitization Using Monoid Curves and Surfaces
Approximate Implicitization Using Monoid Curves and Surfaces
Approximate Implicitization Using Monoid Curves and Surfaces
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
2 <strong>Approximate</strong> implicitization method<br />
Dokken developed a general method to approximately implicitize g-manifolds represented by parametric equations<br />
[11]. Here, we outline his approach.<br />
Let p(s) be a parametric || say, Bezier || curve or surface. We would like to nd an approximate<br />
low-degree implicit equation q(x) =0forp(s). In general, the distance between the parametric curve or surface<br />
<strong>and</strong> its implicit approximation depends not only on the value of q p(s), but also on the gradient of q(x).<br />
However, in practice, the value of jq p(s)j alone can serve as a measure of the goodness of an approximation,<br />
provided the coecients of q are normalized. This reduces the implicitization problem to nding a degree m<br />
algebraic hypersurface q(x) =0with normalized coecients such that jq p(s)j, where is a prescribed<br />
tolerance.<br />
Let S be an arbitrary simplex containing p(s), <strong>and</strong> q(x) = P I b I m!<br />
I! xI be represented in barycentric coordinates<br />
x over S. Then the composition of q(x) = 0 <strong>and</strong> p(s) can be expressed as<br />
q p(s) =(Db) T (s); (1)<br />
where D is a constant matrix, called composite matrix, b is a vector consisting of the coecients of q, <strong>and</strong> (s)<br />
is a vector whose components are basis functions of the composition.<br />
If the components of (s) form a unit partition for s 2 , the following properties hold:<br />
Property 1. max<br />
s2<br />
jq p(s)jkDbk 2<br />
Property 2.<br />
min<br />
kbk 2=1<br />
kDbk 2 = 1 , where 1 0 is the smallest singular value of D.<br />
Thus we can obtain q(x) by setting b to the eigenvector of D D that corresponds to the eigenvalue 2 1 . This<br />
will ensure q(x) satises jq p(s)j when 1