10.07.2015 Views

To the Graduate Council: I am submitting herewith a thesis written by ...

To the Graduate Council: I am submitting herewith a thesis written by ...

To the Graduate Council: I am submitting herewith a thesis written by ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Chapter 4: Algorithm Overview 44where N is <strong>the</strong> normal vector field. The real valued function κ where κ(s) = || α '(s) ||,s <strong>the</strong> curvature function of α and completely describes <strong>the</strong> shape of α in R 2 , up to atranslation and rotation. This curvature function is what we would like to exploit todefine shape information of a curve. We would like to formulate <strong>the</strong> task of curvatureestimation on discrete s<strong>am</strong>ples of such curves. For a planar curve α, we have s<strong>am</strong>plesα j =α(s j ). We assume that uniform s<strong>am</strong>pling across <strong>the</strong> arc length of <strong>the</strong> curve suchthat ∆s = s j - s j-1 is a constant. This approach leads to N s<strong>am</strong>ples over <strong>the</strong> curve α.Since we have uniform s<strong>am</strong>pling along <strong>the</strong> curve κ j is directly proportional to <strong>the</strong>turning angle θ j formed <strong>by</strong> <strong>the</strong> line segments from end point α j-1 to α j and from α j toα j+1 .With 2D curves <strong>the</strong> definition and hence <strong>the</strong> computation of curvature isstraightforward while its extension to 3D surfaces require some concepts indifferential geometry.On a smooth surface S, we can define normal curvature as a starting point. ConsiderFigure 4.4, <strong>the</strong> point p lies on a smooth surface S, and we specify <strong>the</strong> orientation of Sat p with <strong>the</strong> unit-length normal N. We define S as a manifold embedded in R 3 . We canconstruct a plane Π p that contains p and N such that <strong>the</strong> intersection of Π p with Sforms a contour α. As before, we can arbitrarily par<strong>am</strong>eterize α(s) <strong>by</strong> arc length swhere α(0) = p and α’(0)= T. The normal curvature κ p (T) in <strong>the</strong> direction of T is thusα’’ (0) = κ p (T)N .This single κ p (T) does not specify <strong>the</strong> surface curvature of S at psince Π p is not a unique plane. If we rotate Π p around N, we form a new contour on Swith its own normal curvature. We can see that we actually have an infinite set of<strong>the</strong>se normal curvatures around p in every direction. Fortunately, herein enters <strong>the</strong>elegance of surface curvature. For this infinite set, we can construct an orthonormalbasis {T 1 , T 2 } that completely describes <strong>the</strong> set. The natural choice for this basis is <strong>the</strong>tangent vectors associated with <strong>the</strong> maximum and minimum normal curvatures at psince <strong>the</strong> directions of <strong>the</strong>se curvatures are always orthogonal.(a)(b)Figure 4.4: Illustration to understand curvature of a surface.

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

Saved successfully!

Ooh no, something went wrong!