21.01.2015 Views

Approximate Implicitization Using Monoid Curves and Surfaces

Approximate Implicitization Using Monoid Curves and Surfaces

Approximate Implicitization Using Monoid Curves and Surfaces

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!