930 likes | 1.11k Views
Degree reduction of Bézier curves. Lizheng Lu lulz_zju@yahoo.com.cn Mar. 8, 2006. Outline. Overview Recent developments Our work. Problem formulation. Problem I: Given a curve of degree n in , to find a curve of degree m , such that,. An example.
E N D
Degree reduction of Bézier curves Lizheng Lu lulz_zju@yahoo.com.cn Mar. 8, 2006
Outline • Overview • Recent developments • Our work
Problem formulation Problem I: Given a curve of degree n in , to find a curve of degree m, such that,
An example Degree from 7 to 4
Two kinds of methods • Component-wise • Vector decomposition • Degree reduction at each decomposition • Combining all the components • Euclidean [Brunnett et al., 1996] • Consider all the components together
Constrained degree reduction Problem II: Given a curve of degree n in , to find a curve of degree m, such that, I) II)
Metric choice • -norms on C[0,1] • Weighted -norms • Others • Control points perturbing
Approximation Theory • L1-norm • Chebyshev polynomials of second kind • L2-norm • Legendre polynomials • L∞-norm • Chebyshev polynomials of first kind
Lp-norms • L1-norm • Kim and Moon, 1997 • L2-norm • Ahn et al., 2004; Chen and Wang, 2002; Eck, 1995; Zheng and Wang, 2003; Zhang and Wang, 2005; • L∞-norm • Eck, 1993; Ahn, 2003
Present status • Unconstrained • Solved and very mature • Constrained (Optimal approximation) • Solved for L2-norm • Unsolved for L1-norm and L∞-norm • Some methods have been proposed, but not optimal
Outline • Overview • Recent developments • Our work
Multiple degree reduction and elevation of Bézier curves using Jacobi-Bernstein basis transformations Rababah, A., Lee, B.G., Yoo, J. Submitted to Applied Numerical Mathematics
Main contribution • Unified matrix representations • An unconstrained component-wise method • Explicit • Simple and efficient • Explicit approximating error • Include three previous methods • Lp-norm, p =1, 2, ∞
Jacobi polynomials [Szegö, 1975] Orthonormal on [0,1] with
Special kinds of orthonormal polynomials • α=β=-1/2 • Chebyshev polynomials of second kind • α=β=0 • Legendre polynomials • α=β=1/2 • Chebyshev polynomials of first kind
Jacobi-Bernstein basis transformation [Rababah, 2004] Jacobi Bernstein
Degree elevation (1) (2) (3) (4)
Special cases forJacobi-weighted L2-norm • α=β=-1/2 • L∞ -approximation • α=β=0 • L2 -approximation • α=β=1/2 • L1 -approximation
Degree reduction by Jacobi polynomials ? (1) (2)
Degree reduction by Jacobi polynomials ? (1) (2) (3)
Degree reduction by Jacobi polynomials (1) (2) (3) (4)
Summary:Advantages • Unified matrix representations • Explicit approximating error • Optimal approximation under unconstrained degree reduction • Include three previous methods • Lp-norm, p=1, 2, ∞
Summary:Disadvantages • Constrained degree reduction • Unsolved • Challenge to unified representation ?
Outline • Overview • Recent developments • Our work
Optimal multi-degree reduction of Bézier curves with G1-continuity Lizheng Lu and Guozhao Wang To be published in JZUS
Motivation:parametric representations Optimal approximation not unique
Motivation:Geometric Hermite Interpolation • Theorem. [Boor et al., 1987] If the curvature at one endpoint is not vanished, a planar curve can be interpolated by cubic spline with G 2-continuity and that the approximation order is 6. BHS method. More methods about GHI. [Degen, 2005]
Main contributions • Multi-degree reduction • G1: position and tangent direction • Minimize Euclidean distance between control points • Optimal approximation
An example Ahn et al., 2004 Ours
Problem Problem III: Given a curve of degree n in , to find a curve of degree m, such that, I) Gk-continuous: II)
Special cases • G0-continuity • Endpoint interpolation • G1-continuity • Position and tangent direction • G2-continuity • G1 + curvature
Main challenges • Consider all the components together • Control points free moving • How to be optimal? • Error estimating • Numerical problems • Convergence, uniqueness, stability, etc.
Algorithm overview • G1 condition • Discrete coefficient norm • Through degree elevation • Solution and improvement • Numerical methods
G1-continuous Goal Given
G1-continuous Goal Given G1 continuous at both endpoints