Next Article in Journal
A Network-Based Analysis for Evaluating Conditional Covariance Estimates
Next Article in Special Issue
Polynomial-Based Non-Uniform Ternary Interpolation Surface Subdivision on Quadrilateral Mesh
Previous Article in Journal
Some Characterizations of Complete Hausdorff KM-Fuzzy Quasi-Metric Spaces
Previous Article in Special Issue
Generalized de Boor–Cox Formulas and Pyramids for Multi-Degree Spline Basis Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Control Polygons of Planar Sextic Pythagorean Hodograph Curves

1
Information Institute, Zhejiang University of Finance and Economics Dongfang College, Jiaxing 314408, China
2
School of Information Management and Artificial Intelligence, Zhejiang University of Finance and Economics, Hangzhou 310018, China
3
School of Mathematical Sciences, Zhejiang University, Hangzhou 310058, China
4
School of Mathematical Sciences, Xiamen University, Xiamen 361005, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(2), 383; https://doi.org/10.3390/math11020383
Submission received: 11 December 2022 / Revised: 3 January 2023 / Accepted: 4 January 2023 / Published: 11 January 2023
(This article belongs to the Special Issue Computer-Aided Geometric Design)

Abstract

:
In this paper, we analyze planar parametric sextic curves to determine conditions for Pythagorean hodograph (PH) curves. By expressing the curves to be analyzed in the complex form, the analysis is conducted in algebraic form. Since sextic PH curves can be classified into two classes according to the degrees of their derivatives’ factors, we introduce auxiliary control points to reconstruct the internal algebraic structure for both classes. We prove that a sextic curve is completely characterized by the lengths of legs and angles formed by the legs of their Bézier control polygons. As such conditions are invariant under rotations and translations, we call them the geometric characteristics of sextic PH curves. We demonstrate that the geometric characteristics form the basis for an easy and intuitive method for identifying sextic PH curves. Benefiting from our results, the computations of the parameters of cusps and/or inflection points can also be simplified.

1. Introduction

This paper describes an intuitive geometric method for determining whether a planar parametric sextic curve is a Pythagorean hodograph (PH) curve [1]. We call this determination the geometric characteristics of sextic PH curves. A curve is a PH curve if the length of its hodograph coincides with the square of a real polynomial. PH curves have wide-ranging applications due to the special algebraic structures of hodographs. For instance, in numerical control (NC) machining operation, many algorithms rely on the computation of arc lengths or offsets of curves. The arc length of a PH curve can be computed without numerical integration, which speeds up the algorithms for NC machining [2]. The offset of a PH curve can be represented exactly instead of being approximated [3].
PH curves in both polynomial and rational forms have been extensively studied in 2- or 3-dimensional Euclidean spaces [4,5,6,7,8]. Among different representations of PH curves, planar polynomial PH curves are the most widely used due to their simplicity. Hence, the identification of planar polynomial PH curves becomes an important problem in most related applications, e.g., reverse engineering [9].
A fair amount of literature has been focused on developing methods for identifying planar polynomial PH curves of different degrees. They can be classified into algebraic methods and geometric methods. In these methods, PH curves are usually expressed in the Bézier form due to its popularity in CAD systems. Algebraic methods provide conditions for a curve to be a PH curve according to the algebraic properties of PH curves. For example, Farouki et al. determined whether or not a polynomial curve is a PH curve based on the satisfaction of a system of algebraic constraints on its Bézier control polygon legs [9]. Although algebraic methods are straightforward, they are less intuitive in identifying PH curves in practical applications. Geometric methods provide a more flexible and efficient alternative construction or identification of PH curves according to the geometric properties of PH curves. For example, an intuitive geometric condition on the control polygon for a cubic Bézier curve to be a PH curve was proposed in [1]. To obtain more geometric freedoms for practical design problems, the study of geometric constraints on control polygons has proceeded to PH curves of higher degrees [10,11,12,13]. The conditions for control polygons of cubics and quintics in Euclidean and Minkowski spaces have been well studied [14], and the relationship between these curves in Euclidean and Minkowski spaces have been further investigated [15]. Šír [16] studied the classification of planar PH curves with respect to similarities and linear reparameterizations, and he gave general expressions for planar PH curves of low degrees. The work was then employed to discuss the classification of polynomial minimal surfaces [17].
The most-relevant previous work is [13], in which necessary and sufficient conditions for sextic PH curves were provided. Based on these conditions, three methods for the identification of planar sextic PH curves are proposed, where the first method relies on solving a system of degree eleven, the second method can only deal with convex control polygons, and the third method is limited to some special cases of planar sextic PH curves. In this paper, we provide a more succinct and complete description of the geometric characteristics of sextic PH curves. Following the geometric characteristics, we propose an efficient method for the identification of sextic PH curves. The key technique for making our method both intuitive and flexible is the introduction of auxiliary control points. In particular, we geometrically compute auxiliary control points for given sextic curves and obtain auxiliary control polygons. The geometric characteristics are then described using the leg lengths of the auxiliary control polygons and angles formed by the legs. Our method is suitable for general sextic polynomial curves, and we only resort to solve lower-degree (quadratic) systems.
The remainder of this paper is organized as follows. Section 2 introduces some fundamental aspects of sextic PH curves in the Bézier form. Section 3 presents some geometric properties of sextic PH curves. In Section 4, sextic PH curves are classified into two classes and their geometric characteristics are presented, respectively. Finally, in Section 5, we conclude the paper.

2. Bézier Representation for Sextic PH Curves

In this section, we recall the Bézier representation for planar parametric sextic PH curves. We first introduce some notations and, then, use the complex representation of R 2 to facilitate the computation in the later analysis of planar PH curves [4,11].
Let R and C denote the sets of all real numbers and complex numbers, respectively. A complex number is denoted by a single bold character, e.g., z . The complex norm, the real part, the imaginary part, and the argument of z in C are denoted by z , Re z , Im z , and arg z [ 0 , 2 π ) , respectively. A planar point ( x , y ) is interchangeably represented as a complex number x + i y . Similarly, a planar parametric curve P ( t ) = ( x ( t ) , y ( t ) ) , t [ 0 , 1 ] , can be identified with a complex-valued function P ( t ) = x ( t ) + i y ( t ) , and vice versa.
The Bézier representation for a planar sextic polynomial curve is defined by seven control points in the plane, { P i } i = 0 6 , as
P ( t ) = ( x ( t ) , y ( t ) ) = i = 0 6 B i 6 ( t ) P i , 0 t 1 ,
where B i n ( t ) = n i ( 1 t ) n i t i , i = 0 , , n , are the Bernstein polynomials. The polygon formed by consecutive control points { P i } i = 0 6 is called the control polygon. Let Δ P i denote the first forward-difference of the i-th control point, i.e., Δ P i = P i + 1 P i ; thus, the derivative of the curve (1) can be represented as
P ( t ) = ( x ( t ) , y ( t ) ) = 6 i = 0 5 B i 5 ( t ) Δ P i .
A planar sextic Bézier curve (1) is called a Pythagorean hodograph (PH) curve if and only if its hodograph P ( t ) = x ( t ) + i y ( t ) satisfies the Pythagorean condition:
x 2 ( t ) + y 2 ( t ) = σ 2 ( t ) ,
for some real polynomial σ ( t ) [1]. Equivalently, a planar sextic Bézier curve (1) is a PH curve if and only if
P ( t ) = w ( t ) [ u ( t ) + i v ( t ) ] 2 ,
for some real polynomials u ( t ) , v ( t ) , and w ( t ) , where u ( t ) and v ( t ) are relatively prime [12].

3. Properties of Sextic PH Curves

In this section, we investigate the presence of the cusps, inflection points, and curvature of a sextic PH curve. According to the sufficient and necessary condition (3) for a sextic curve to be PH, the degrees of u ( t ) + i v ( t ) and w ( t ) are either 2 and 1 or 1 and 3, which are mutually exclusive. The sextic PH curves corresponding to these two cases are referred to as Class I and Class II sextic PH curves, respectively, and will be discussed in detail in this section.
Class I
Polynomials u ( t ) + i v ( t ) and w ( t ) are quadratic and linear, respectively. Regardless of the symmetry of control points, we rewrite the factors of P ( t ) in Bernstein form as
w ( t ) = a 0 ( 1 t ) + t , u ( t ) + i v ( t ) = i = 0 2 z i B i 2 ( t ) ,
where z 0 , z 2 0 . Remark: if z 0 = 0 or z 2 = 0 , the curve is indeed a Class II sextic PH curve, which will be discussed later. If a 0 1 , then w ( t ) is a linear polynomial with real coefficients, which has exactly one real root. Denoting by ξ 0 the root of w ( t ) , we have
a 0 = ξ 0 ξ 0 1 .
That is to say, a sextic PH curve with the hodograph of the form (4) has a cusp occurring at ξ 0 . Otherwise, if a 0 = 1 (i.e., w ( t ) 1 ), the corresponding sextic curve degenerates to a Class I quintic PH curve [10]. On the contrary, a quintic Class I PH curve can also be considered as a sextic PH curve of this class using one step of degree elevation.
Class II
Polynomials u ( t ) + i v ( t ) and w ( t ) are linear and cubic, respectively. We rewrite the factors of P ( t ) in Bernstein form as
w ( t ) = i = 0 3 a i B i 3 ( t ) , u ( t ) + i v ( t ) = z 0 ( 1 t ) + z 1 t ,
for some a i R , i = 0 , , 3 , and z j C , j = 0 , 1 . Let ξ i , i = 0 , 1 , 2 , be roots of the cubic polynomial w ( t ) in Equation (6). According to the fundamental theorem of algebra, there is at least one real root. We assume ξ 0 R , then ξ 1 and ξ 2 are either real numbers or a pair of conjugate complex numbers. Hence, a Class II sextic PH curve presents one or three cusps, depending on the real root number of w ( t ) . Note that, when the cubic real polynomial w ( t ) degenerates to a quadratic, linear, or constant polynomial, the sextic PH curve degenerates to a Class II quintic, quartic, or cubic PH curve, respectively. Equally, cubic, quartic, and Class II quintic PH curves can be converted to Class II sextic PH curves using three steps to one step of degree elevation. A Class II sextic PH curve degenerates to a line segment when z 0 = 0 or z 1 = 0 , which is not discussed in this paper.
Let us now consider the curvature of sextic PH curves. In the following, the variables of functions in the notations are omitted for simplicity, e.g., u means u ( t ) . Following (3), the second derivative P = ( x , y ) is
x = w ( u 2 v 2 ) + 2 w ( u u v v ) , y = 2 ( u v w + u v w + u v w ) .
Therefore,
P × P = 2 w ( u 2 v 2 ) ( u v w + u v w + u v w ) 2 u v w [ w ( u 2 v 2 ) + 2 w ( u u v v ) ] = 2 u 2 u v w 2 2 u v 3 w 2 + 2 u 3 v w 2 2 u v 2 v w 2 2 u 2 u v w 2 + 2 u v 2 v w 2 = 2 v 2 w 2 ( u v u v ) + 2 u 2 w 2 ( u v u v ) = 2 w 2 ( u 2 + v 2 ) ( u v u v ) .
Thus, its curvature k can be computed by
k = | P × P | | P | 3 = u v u v w u 2 + v 2 .
Note that linear polynomials u, v lead to a constant numerator of curvature u v u v , which means that the curvature does not vanish, i.e., no inflections occur on Class II sextic PH curves. Furthermore, we may claim that there are no inflection points on cubic [1], quartic [12], Class II quintic [10], and Class III septic [18] PH curves, i.e., they either are convex curves or have cusps.

4. Geometric Characteristics of Sextic PH Curves

In this section, we describe the geometric characteristics of sextic PH curves using auxiliary control points, based on which we can identify whether a given sextic polynomial curve is a PH curve or not. We introduce auxiliary control points, which are linear combinations of the original control points. Following that, we discuss the sufficient and necessary conditions for a planar sextic Bézier curve to possess a PH based on the auxiliary control points.

4.1. Class I Sextic PH Curve

For a Class I sextic curve in Bézier representation, we introduce eight auxiliary control points Q i for i = 0 , , 7 , by the following four steps:
  • Let α = arg Δ P 0 + arg Δ P 5 2 , and
    Q 0 = P 1 + k 0 Δ P 0 , Q 7 = P 5 k 1 Δ P 5 ,
    where
    k 0 = tan α · Re Δ P 1 Im Δ P 1 tan α · Re Δ P 0 Im Δ P 0 , k 1 = tan α · Re Δ P 4 Im Δ P 4 tan α · Re Δ P 5 Im Δ P 5 .
  • Let
    c = P 4 Q 7 2 Δ P 0 · Δ P 5 , d = P 2 Q 0 2 Δ P 0 · Δ P 5 ;
    thus, we obtain
    Q 1 = P 1 + Δ P 0 5 x , Q 6 = P 5 x 5 Δ P 5 ,
    where x is a solution of the system of cubic equations as follows:
    x 3 10 k 1 x 2 + 25 k 1 2 x 25 d = 0 , 25 c x 3 25 k 0 2 x 2 + 10 k 0 x 1 = 0 , 5 k 0 x 2 + ( ± 25 c d 25 k 0 k 1 1 ) x + 5 k 1 = 0 .
    We defer deriving this system until the proof of Theorem 1. Note that the roots of the third equation of System (13) come in pairs:
    x 1 = ( ± 25 c d 25 k 0 k 1 1 ) + ( ± 25 c d 25 k 0 k 1 1 ) 2 100 k 0 k 1 10 k 0 , x 2 = ( ± 25 c d 25 k 0 k 1 1 ) ( ± 25 c d 25 k 0 k 1 1 ) 2 100 k 0 k 1 10 k 0 .
    Hence, we solve the system (13) by verifying if the roots of the third equation of Equation (13) satisfy the first two equations. Although there are four candidate solutions due to different signs, the number of roots is no more than three for the cubic system. For each real solution, Steps 3–4 are further performed, respectively. If there are no real roots for the above system, the curve is not a PH curve, and the procedure is terminated immediately.
  • Let Q 2 and Q 5 be points on lines Q 1 P 2 and Q 6 P 4 , respectively, such that
    Q 2 = P 2 + P 2 Q 1 2 x , Q 5 = P 4 x 2 ( Q 6 P 4 ) .
  • Finally, let β = arg ( P 2 Q 1 ) + arg ( Q 6 P 4 ) 2 , and
    Q 3 = P 3 + Im ( Q 2 P 3 ) · cos β Re ( Q 2 P 3 ) · sin β sin ( α β ) · e i α = Q 2 + Im ( Q 2 P 3 ) · cos α Re ( Q 2 P 3 ) · sin α sin ( α β ) · e i β , Q 4 = P 3 + Im ( Q 5 P 3 ) · cos β Re ( Q 5 P 3 ) · sin β sin ( α β ) · e i α = Q 5 + Im ( Q 5 P 3 ) · cos α Re ( Q 5 P 3 ) · sin α sin ( α β ) · e i β .
Now, we give a necessary and sufficient condition for Class I sextic PH curves.
Theorem 1.
Given a planar sextic curve with Bézier control points P i , i = 0 , , 6 , it is a Class I sextic PH curve if and only if there are auxiliary points Q j , j = 1 , , 7 , constructed following (12)(16), such that Q 2 Q 3 Q 4 Q 5 , and they further satisfy
5 P 2 Q 1 2 = 8 Δ P 0 · Δ Q 2 , 5 Q 6 P 4 2 = 8 Δ Q 4 · Δ P 5 , Q 2 P 2 2 = 2 Q 1 P 1 · Δ Q 4 , P 4 Q 5 2 = 2 Δ Q 2 · P 5 Q 6 , 5 P 3 Q 3 2 = Δ P 0 · P 5 Q 6 , 5 Q 4 P 3 2 = Q 1 P 1 · Δ P 5 .
Proof. 
A planar sextic curve with Bézier control points P i , i = 0 , , 6 , is a Class I sextic PH curve if and only if its derivative with factors in Equation (4) satisfies Equation (3), in other words, if and only if there are a 0 , z 0 , z 1 , z 2 satisfying:
P ( t ) = a 0 z 0 2 ( 1 t ) 5 + ( z 0 2 + 4 a 0 z 0 z 1 ) ( 1 t ) 4 t + ( 4 z 0 z 1 + 4 a 0 z 1 2 + 2 a 0 z 0 z 2 ) ( 1 t ) 3 t 2 + ( 4 z 1 2 + 2 z 0 z 2 + 4 a 0 z 1 z 2 ) ( 1 t ) 2 t 3 + ( 4 z 1 z 2 + a 0 z 2 2 ) ( 1 t ) t 4 . + z 2 2 t 5 .
By matching the coefficients of Bernstein polynomials in the above equations with Equation (2), we obtain the decompositions of the control polygon legs:
1 6 a 0 z 0 2 = Δ P 0 , 1 30 ( z 0 2 + 4 a 0 z 0 z 1 ) = Δ P 1 , 1 60 ( 4 z 0 z 1 + 4 a 0 z 1 2 + 2 a 0 z 0 z 2 ) = Δ P 2 , 1 60 ( 4 z 1 2 + 2 z 0 z 2 + 4 a 0 z 1 z 2 ) = Δ P 3 , 1 30 ( 4 z 1 z 2 + a 0 z 2 2 ) = Δ P 4 , 1 6 z 2 2 = Δ P 5 .
We define six auxiliary control points Q ˜ i , i = 1 , , 6 :
Q ˜ 1 = P 1 + 1 30 z 0 2 = P 2 2 a 0 15 z 0 z 1 , Q ˜ 2 = P 2 + 1 15 z 0 z 1 , Q ˜ 3 = Q ˜ 2 + a 0 15 z 1 2 = P 3 a 0 30 z 0 z 2 , Q ˜ 4 = P 3 + 1 30 z 0 z 2 , Q ˜ 5 = Q ˜ 4 + 1 15 z 1 2 = P 4 a 0 15 z 1 z 2 , Q ˜ 6 = P 4 + 2 15 z 1 z 2 = P 5 a 0 30 z 2 2 .
Let Q ˜ 0 , Q ˜ 7 be points on lines P 0 P 1 and P 5 P 6 , respectively, such that
arg P 2 Q ˜ 0 Δ P 0 = arg Δ P 5 Q ˜ 7 P 4 = 1 2 arg Δ P 5 Δ P 0 , Re ( Q ˜ 0 P 2 ) Im ( P 4 Q ˜ 7 ) = Im ( Q ˜ 0 P 2 ) Re ( P 4 Q ˜ 7 ) ,
that is Q ˜ 0 P 2 P 4 Q ˜ 7 . By computing the intersecting point of lines P 0 P 1 and Q ˜ 0 P 2 , we obtain Q ˜ 0 = Q 0 ; similarly, we have Q ˜ 7 = Q 7 ; see Figure 1a. Note that triangles P 2 Q 0 Q ˜ 1 and Q ˜ 6 Q 7 P 4 are similar triangles, as their edges are parallel with z 2 2 , z 0 z 2 , z 1 z 2 , respectively. Hence, the ratios of the lengths of their corresponding edges are equal:
Q ˜ 1 Q 0 P 4 Q 7 = P 2 Q 0 Q 7 Q ˜ 6 = P 2 Q ˜ 1 P 4 Q ˜ 6 .
Following Equations (18) and (19), we have
Δ P 0 5 ( Q ˜ 1 P 1 ) = 5 ( P 5 Q ˜ 6 ) Δ P 5 = a 0 , P 2 Q ˜ 1 2 P 4 Q ˜ 6 2 = a 0 2 z 0 2 z 2 2 = | a 0 | · Δ P 0 Δ P 5 .
Using the notations from Equation (10), we have
Q ˜ 1 Q ˜ 0 = Q ˜ 1 P 1 + P 1 Q ˜ 0 = | k 0 1 5 a 0 | · Δ P 0 , Q ˜ 7 Q ˜ 6 = Q ˜ 7 P 5 + P 5 Q ˜ 6 = | k 1 a 0 5 | · Δ P 5 .
Substituting Equation (11) into Equation (21) and, then, squaring both sides of Equation (21), we obtain
| k 0 1 5 a 0 | 2 c = d | k 1 a 0 5 | 2 = | a 0 | .
Thus, the cubic equations in Equation (13) can be derived by considering every two items to be equal. In other words, there is a real root of a 0 such that Q ˜ i = Q i , i = 1 , 2 , 5 , 6 , as shown in Figure 1b,c.
From Equation (19), we have arg ( P 3 Q ˜ 3 ) = arg ( Q ˜ 4 P 3 ) = arg z 0 z 2 , and
arg P 3 Q ˜ 3 Δ P 0 = arg Δ P 5 Q ˜ 4 P 3 = arg z 2 z 0 , arg Δ Q ˜ 2 P 2 Q 1 = arg P 2 Q 1 Δ P 0 = arg z 1 z 0 , arg Δ P 5 Q 6 P 4 = arg Q 6 P 4 Δ Q ˜ 4 = arg z 2 z 1 ;
thus, we know the points Q ˜ 3 , P 3 , Q ˜ 4 are co-linear, and we have Q ˜ 3 Q ˜ 4 Q 0 P 2 P 4 Q 7 . Since non-parallel lines intersect at a unique point, we have Q ˜ i = Q i , i = 3 , 4 , where Q i are defined in Equation (16); see Figure 1d. Thus far, we have proven that Q ˜ i , i = 1 , , 6 , coincident with the auxiliary control points Q i constructed in Equations (12)–(16); see Figure 1d. Moreover, since arg Δ Q 2 = arg Δ Q 4 = arg z 1 2 , we have Q 2 Q 3 Q 4 Q 5 . Following Equation (19), Equation (17) can be computed as follows:
5 P 2 Q 1 2 = 8 Δ P 0 · Δ Q 2 = 4 45 a 0 2 z 0 2 z 1 2 , 5 Q 6 P 4 2 = 8 Δ Q 4 · Δ P 5 = 4 45 z 1 2 z 2 2 , Q 2 P 2 2 = 2 Q 1 P 1 · Δ Q 4 = 1 225 z 0 2 z 1 2 , P 4 Q 5 2 = 2 Δ Q 2 · P 5 Q 6 = 1 225 a 0 2 z 1 2 z 2 2 , 5 P 3 Q 3 2 = Δ P 0 · P 5 Q 6 = 1 180 a 0 2 z 0 2 z 2 2 , 5 Q 4 P 3 2 = Q 1 P 1 · Δ P 5 = 1 180 z 0 2 z 2 2 .
On the contrary, given a planar sextic curve with Bézier control points P i , i = 0 , , 6 , if Equation (17) holds for auxiliary control points Q i , i = 0 , , 7 , we have
a 0 = Δ P 0 5 ( Q 1 P 1 ) = P 2 Q 1 2 ( Q 2 P 2 ) = Δ Q 2 Δ Q 4 = P 3 Q 3 Q 4 P 3 = 2 ( P 4 Q 5 ) Q 6 P 4 = 5 ( P 5 Q 6 ) Δ P 5 .
Note that, although the denominators and numerators are all complex numbers, the quotient a 0 is well-defined because the denominator and numerator in each item have the same arguments. Moreover, the complex coefficients z i , i = 0 , 1 , 2 , can be solved by Equation (19), that is
z 0 = ± 30 ( Q 1 P 1 ) , z 1 = ± 15 ( Q 2 P 2 ) 30 ( Q 1 P 1 ) , z 2 = ± 30 ( Q 4 P 3 ) 30 ( Q 1 P 1 ) .
By substituting Δ P i , i = 0 , , 5 , into Equation (2), we can verify that the curve is a Class I sextic PH curve because its hodograph is actually of the form (3), which proves the sufficiency. □
The construction of auxiliary control points depends only on the leg lengths and angles of the Bézier control polygon, which can be obtained by measurement; hence, our method is geometric-invariant.

An Example of a Class I Sextic PH Curve with a Concave Control Polygon

Now, we give an example of a Class I sextic PH curve with a concave control polygon, as shown in Figure 2. In this example, we assign the control points as follows:
P 0 = 0 + 0 i , P 1 = 1.228732694000000 + 2.054457192000000 i , P 2 = 1.024962887533333 + 5.016989645599999 i , P 3 = 0.519266461400000 + 7.371566599600000 i , P 4 = 2.371552422866667 + 10.372970071533333 i , P 5 = 7.268585548066668 + 11.789158331266666 i , P 6 = 12.868452536400000 + 4.155075291266666 i .
Following Equation (10), we compute
k 0 = 1.488475973647499 , k 1 = 0.150203163391209 ,
which immediately gives
Q 0 = 3.057671787054164 + 5.112467361179306 i Q 7 = 1.433521476333333 + 2.396866724000000 i .
By solving Equation (13), we obtain the unique solution:
a 0 = x = 1.200000000000002 ;
thus, we have
Q 1 = 1.433521476333333 2.396866724000000 i , Q 2 = 0.854730142200000 6.108707529599998 i , Q 3 = 0.522605614600000 7.322629161599998 i , Q 4 = 1.387493191400000 7.412347797933333 i , Q 5 = 0.239713394066667 8.423949157933333 i , Q 6 = 5.924617470866667 13.621338260866665 i .
We numerically verify that the leg lengths and angles of control polygon with auxiliary control points defined above satisfy the conditions in Theorem 1 within tolerance 10 13 [9]; thus, the corresponding curve is a Class I sextic PH curve. Moreover, we can further compute the complex coefficients of Equation (4):
z 0 = ± ( 3.009399999999999 1.706700000000001 i ) , z 1 = ± ( 1.693000000000000 4.481399999999998 i ) , z 2 = ± ( 6.723300000000002 + 3.406399999999999 i ) ,
which can facilitate algebraic computations of the curves, e.g., arc lengths, curvatures (8), etc.

4.2. Class II Sextic PH Curve

In this section, we propose the geometric characteristics of Class II sextic PH curves in the same fashion as we did in the previous section for Class I sextic PH curves. The auxiliary control points Q i , i = 1 , , 6 , for Class II sextic PH curves are constructed according to the following steps:
  • Let α = arg Δ P 0 + arg Δ P 5 2 ; we obtain
    Q 1 = P 1 + k 0 Δ P 0 , Q 6 = P 5 k 1 Δ P 5 ,
    where
    k 0 = Re Δ P 1 · tan α Im Δ P 1 Re Δ P 0 · tan α Im Δ P 0 , k 1 = Re Δ P 4 · tan α Im Δ P 4 Re Δ P 5 · tan α Im Δ P 5 .
  • Compute
    D Δ Q 2 = 5 ( Q 1 P 1 ) 2 Δ P 0 · ( P 2 Q 1 ) , D Δ Q 4 = 5 ( P 5 Q 6 ) 2 Δ P 5 · ( Q 6 P 4 ) , D P 2 Q 2 = D Δ Q 4 5 ( P 2 Q 1 ) · Δ P 0 , D P 3 Q 4 = ( Q 6 P 4 ) 10 ( P 2 Q 1 ) · Δ P 0 , D Q 3 P 3 = ( P 2 Q 1 ) 10 ( Q 6 P 4 ) · Δ P 5 , D Q 5 P 4 = D Δ Q 2 5 ( Q 6 P 4 ) · Δ P 5 .
  • Let Q i , i = 2 , , 5 , be points such that
    Q 2 = P 2 + D P 2 Q 2 , Q 3 = Q 2 + D Δ Q 2 , Q 4 = P 3 + D P 3 Q 4 , Q 5 = Q 4 + D Δ Q 4 ;
    thus, we have P 0 P 1 P 2 Q 2 P 3 Q 4 , Q 1 P 2 Q 2 Q 3 Q 4 Q 5 P 4 Q 6 , Q 3 P 3 Q 5 Q 4 P 5 P 6 , as shown in Figure 3.
The following theorem describes the necessary and sufficient conditions for a sextic curves to be a Class II PH curve.
Theorem 2.
Let P i , i = 0 , , 6 , be Bézier control points of a planar sextic curve and Q i , i = 1 , , 6 , be auxiliary control points constructed following Equations (23)(26), then it is a Class II sextic PH curve if and only if
5 P 2 Q 1 2 = 8 Δ P 0 · P 3 Q 3 , Δ Q 2 2 = 2 Q 1 P 1 · P 4 Q 5 , Δ Q 4 2 = 2 Q 2 P 2 · P 5 Q 6 , 5 Q 6 P 4 2 = 8 Q 4 P 3 · Δ P 5 .
Proof. 
For a Class II sextic PH curve, there are a i R , i = 0 , , 3 , z 0 , and z 1 such that its hodograph has the form (3), where u ( t ) + i v ( t ) , w ( t ) are polynomials of the form (6). Since all control points are distinct, we assume that a 3 = 1 ; thus, Equation (3) can be rewritten as
P ( t ) = a 0 z 0 2 ( 1 t ) 5 + ( 2 a 0 z 0 z 1 + 3 a 1 z 0 2 ) ( 1 t ) 4 t + ( a 0 z 1 2 + 6 a 1 z 0 z 1 + 3 a 2 z 0 2 ) ( 1 t ) 3 t 2 + ( 3 a 1 z 1 2 + 6 a 2 z 0 z 1 + z 0 2 ) ( 1 t ) 2 t 3 + ( 3 a 2 z 1 2 + 2 z 0 z 1 ) ( 1 t ) t 4 + z 1 2 t 5 .
By matching the coefficients of Bernstein polynomials in the above equation with Equation (2), we have
1 6 a 0 z 0 2 = Δ P 0 , 1 30 ( 2 a 0 z 0 z 1 + 3 a 1 z 0 2 ) = Δ P 1 , 1 60 ( a 0 z 1 2 + 6 a 1 z 0 z 1 + 3 a 2 z 0 2 ) = Δ P 2 , 1 60 ( 3 a 1 z 1 2 + 6 a 2 z 0 z 1 + z 0 2 ) = Δ P 3 , 1 30 ( 3 a 2 z 1 2 + 2 z 0 z 1 ) = Δ P 4 , 1 6 z 1 2 = Δ P 5 .
Let points Q ˜ i , i = 1 , , 6 , be defined by
Q ˜ 1 = P 1 + a 1 10 z 0 2 = P 2 a 0 15 z 0 z 1 , Q ˜ 2 = P 2 + a 2 20 z 0 2 , Q ˜ 3 = Q ˜ 2 + a 1 10 z 0 z 1 = P 3 a 0 60 z 1 2 , Q ˜ 4 = P 3 + 1 60 z 0 2 , Q ˜ 5 = Q ˜ 4 + a 2 10 z 0 z 1 = P 4 a 1 20 z 1 2 , Q ˜ 6 = P 4 + 1 15 z 0 z 1 = P 5 a 2 10 z 1 2 .
Thus, points Q ˜ 1 and Q ˜ 6 are located on lines P 0 P 1 and P 5 P 6 , respectively, and satisfy
arg P 2 Q ˜ 1 Δ P 0 = arg Δ P 5 Q ˜ 6 P 4 = arg z 1 z 0 = 1 2 arg Δ P 5 Δ P 0 , Re ( P 2 Q 1 ) Im ( P 4 Q 6 ) = Im ( P 2 Q 1 ) Re ( P 4 Q 6 ) ,
that is P 2 Q 1 P 4 Q 6 .
Note that points Q ˜ i and Q i are coincident with each other for i = 1 , 6 as the lines intersect at a unique point in a Euclidean plane. From Equation (29), we have P 0 P 1 P 2 Q ˜ 2 P 3 Q ˜ 4 , Q ˜ 1 P 2 Q ˜ 2 Q ˜ 3 Q ˜ 4 Q ˜ 5 P 4 Q ˜ 6 , Q ˜ 3 P 3 Q ˜ 5 P 4 P 5 P 6 , and
3 Δ P 0 : 5 ( Q ˜ 1 P 1 ) : 10 ( Q ˜ 2 P 2 ) : 30 ( Q ˜ 4 P 3 ) = 3 ( P 2 Q ˜ 1 ) : 2 Δ Q ˜ 2 : 2 Δ Q ˜ 4 : 3 ( Q ˜ 6 P 4 ) = 30 ( P 3 Q ˜ 3 ) : 10 ( P 4 Q ˜ 5 ) : 5 ( P 5 Q ˜ 6 ) : 3 Δ P 5 = a 0 : a 1 : a 2 : 1 .
We can verify that Q ˜ i coincides with Q i defined in Equation (26) for i = 2 , , 5 . In other words, points Q ˜ i , i = 1 , , 6 , are exactly the points constructed by Steps 1–3 in this section. Moreover, we have
5 P 2 Q 1 2 = 8 Δ P 0 · P 3 Q 3 = 1 45 a 0 2 z 0 2 z 1 2 , Δ Q 2 2 = 2 Q 1 P 1 · P 4 Q 5 = 1 100 a 1 2 z 0 2 z 1 2 , Δ Q 4 2 = 2 Q 2 P 2 · P 5 Q 6 = 1 100 a 2 2 z 0 2 z 1 2 , 5 Q 6 P 4 2 = 8 Q 4 P 3 · Δ P 5 = 1 45 z 0 2 z 1 2 .
Thus, we have proven the necessity.
On the contrary, if the conditions in Theorem 2 hold for a given planar sextic curve, we can obtain a group of solutions for a i , i = 0 , 1 , 2 , 3 , z 0 and z 1 from Equation (29), e.g.,
a 0 = 10 ( P 3 Q 3 ) Δ P 5 , a 1 = 10 ( P 4 Q 5 ) 3 Δ P 5 , a 2 = 5 ( P 5 Q 6 ) 3 Δ P 5 , a 3 = 1 , z 0 = ± 15 ( Q 6 P 4 ) 6 Δ P 5 , z 1 = ± 6 Δ P 5 ;
therefore, the curve is a Class II sextic PH curve. □

An Example of a Class II Sextic PH Curve with a Concave Control Polygon

In this section, we show an example of a Class II sextic PH curve with a concave control polygon given by (see Figure 4)
P 0 = 0 + 0 i , P 1 = 2.000009978666667 5.999907840000000 i , P 2 = 0.915235119476400 8.304442981248000 i , P 3 = 2.922077810388467 19.175025906393600 i , P 4 = 14.418872723033649 18.001271454000001 i , P 5 = 10.749764840044783 5.749685950089601 i , P 6 = 8.749765161711450 1.749702910089601 i .
Following Equation (24), we have
k 0 = 0.406575000000000 , k 1 = 3.020760000000000 .
Then, from Equations (23) and (26), we have
Q 1 = 2.813164035743067 8.439320370048000 i , Q 2 = 4.691203958949600 19.632118985472001 i , Q 3 = 2.762077836121800 19.495024549593602 i , Q 4 = 3.172079057721800 19.925014386393599 i , Q 5 = 14.744132670720999 18.651788695795201 i , Q 6 = 16.791283868366982 17.832674718000000 i .
We can numerically verify that Theorem 2 holds; thus, the curve shown in Figure 4 is a Class II sextic PH curve. Moreover, we obtain the coefficients of Equation (6) as follows:
a 0 = 0.800000000000004 , a 1 = 0.542100000000000 , a 2 = 5.034599999999999 , z 0 = ± ( 4.027000000000003 + 5.587199999999996 i ) , z 1 = ± ( 2.723300000000000 4.406400000000000 i ) .

5. Conclusions

In this paper, we developed a simple and intuitive geometric method for identifying sextic PH curves according to the geometric characteristics of sextic PH curves. The geometric characteristics were described by the leg lengths and angles of the Bézier control polygon with auxiliary control points. Different from previous work, which was limited to convex control polygons, our method applies to general sextic polynomial curves. While we focused on planar polynomial curves here, we believe that this work could be extended to Minkowski PH (MPH) curves [14,15], rational PH curves, and indirect PH curves [19]. These will be considered as research directions in our future work.

Author Contributions

Conceptualization, L.F.; Methodology, Y.L.; Software, Y.L. and L.F.; Validation, L.F.; Investigation, Z.Z.; Writing—original draft, Y.L.; Writing—review & editing, L.F.; Supervision, J.C.; Funding acquisition, J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China grant number 61872308, 61972327, 62272402.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors would like to thank the anonymous reviewers for their valuable comments and suggestions that helped to improve the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Farouki, R.T.; Sakkalis, T. Pythagorean hodographs. IBM J. Res. Dev. 1990, 34, 736–752. [Google Scholar] [CrossRef]
  2. Farouki, R.T.; Manjunathaiah, J.; Nichlas, D.; Yuan, G.F.; Jee, S. Variable-feedrate CNC interpolators for constant material removal rates along Pythagorean hodograph curves. Comput.-Aided Des. 1998, 30, 631–640. [Google Scholar] [CrossRef]
  3. Elber, G.; Lee, I.-K.; Kim, M.-S. Comparing offset curve approximation methods. IEEE Comput. Graph. Appl. 1997, 17, 62–71. [Google Scholar] [CrossRef] [Green Version]
  4. Farouki, R.T. Pythagorean-Hodograph Curves: Algebra and Geometry Inseparable; Springer: Berlin/Heidelberg, Germany, 2008. [Google Scholar]
  5. Farouki, R.T.; Giannelli, C.; Manni, C.; Sestini, A. Identification of spatial PH quintic Hermite interpolants with near-optimal shape measures. Comput. Aided Geom. Des. 2008, 25, 274–297. [Google Scholar] [CrossRef]
  6. Farouki, R.T.; Šír, Z. Rational Pythagorean hodograph space curves. Comput. Aided Geom. Des. 2011, 28, 75–88. [Google Scholar] [CrossRef]
  7. Kosinka, J.; Lávička, M. Pythagorean hodograph curves: A survey of recent advances. J. Geom. Graph. 2014, 18, 23–43. [Google Scholar]
  8. Sakkalis, T.; Farouki, R.T. Pythagorean hodograph curves in Euclidean spaces of dimension greater than 3. J. Comput. Appl. Math. 2012, 236, 4375–4382. [Google Scholar] [CrossRef] [Green Version]
  9. Farouki, R.T.; Giannelli, C.; Sestini, A. Identification and “reverse engineering” of Pythagorean-hodograph curves. Comput. Aided Geom. Des. 2015, 34, 21–36. [Google Scholar] [CrossRef] [Green Version]
  10. Fang, L.; Wang, G. Geometric characteristics of planar quintic Pythagorean hodograph curves. J. Comput. Appl. Math. 2018, 330, 117–127. [Google Scholar] [CrossRef]
  11. Farouki, R.T. The conformal map zz2 of the hodograph plane. Comput. Aided Geom. Des. 1994, 11, 363–390. [Google Scholar] [CrossRef]
  12. Wang, G.; Fang, L. On control polygons of quartic Pythagorean hodograph curves. Comput. Aided Geom. Des. 2009, 26, 1006–1015. [Google Scholar] [CrossRef]
  13. Wang, H.; Zhu, C.; Li, C. Identification of planar sextic Pythagorean-Hodograph curves. J. Math. Res. Appl. 2017, 37, 59–72. [Google Scholar]
  14. Šír, Z.; Kosinka, J. Low degree Euclidean and Minkowski Pythagorean hodograph curves. In Proceedings of the International Conference on Mathematical Methods for Curves and Surfaces, Tønsberg, Norway, 26 June–1 July 2008; Dæhlen, M., Floater, M., Lyche, T., Jean-Louis, M., Mørken, K., Schumaker, L.L., Eds.; Springer: Berlin/Heidelberg, Germany, 2010; pp. 394–418. [Google Scholar]
  15. Kosinka, J.; Lávička, M. On rational Minkowski Pythagorean hodograph curves. Comput. Aided Geom. Des. 2010, 27, 514–524. [Google Scholar] [CrossRef] [Green Version]
  16. Šír, Z. Classification of planar Pythagorean hodograph curves. Comput. Aided Geom Des. 2020, 80, 101866. [Google Scholar] [CrossRef]
  17. Fang, L.; Peng, Y.; Li, Y.; Cao, J. Classification of polynomial minimal surfaces. Comput. Aided Geom. Des. 2022, 96, 102106. [Google Scholar] [CrossRef]
  18. Li, Y.; Fang, L.; Cao, J.; Wang, G. Identification of two classes of planar septic Pythagorean hodograph curves. J. Comput. Appl. Math. 2019, 348, 383–400. [Google Scholar] [CrossRef]
  19. Lu, X.-J.; Zheng, J.; Cai, Y.; Zhao, G. Geometric characteristics of a class of cubic curves with rational offsets. Comput.-Aided Des. 2016, 70, 36–45. [Google Scholar] [CrossRef]
Figure 1. Construction of auxiliary points Q i , i = 0 , , 7 , for a Class I sextic PH curve with its Bézier control points P i , i = 0 , , 6 . (a) Points Q 0 and Q 7 are on lines P 0 P 1 and P 5 P 6 , respectively, such that arg P 2 Q 0 Δ P 0 = arg Δ P 5 Q 7 P 4 and Q 0 P 2 P 4 Q 7 . (b) Points Q 1 and Q 6 are on lines P 0 P 1 and P 5 P 6 , respectively, such that the triangles Q 0 Q 1 P 2 and Q 7 P 4 Q 6 are similar; thus, we have Δ Q 0 Q 7 P 4 = P 2 Q 0 Δ Q 6 = P 2 Q 1 Q 6 P 4 . (c) Points Q 2 and Q 5 are on lines Q 1 P 2 and Q 6 P 4 , respectively, such that P 2 Q 1 2 ( Q 2 P 2 ) = Δ P 0 5 ( Q 1 P 1 ) and 2 ( P 4 Q 5 ) Q 6 P 4 = 5 ( P 5 Q 6 ) Δ P 5 . (d) The line Q 3 Q 4 is passing through the point P 3 , such that Q 3 Q 4 Q 0 P 2 P 4 Q 7 , and they further satisfy arg Δ Q 2 P 2 Q 1 = arg P 2 Q 1 Δ P 0 and arg Δ P 5 Q 6 P 4 = arg Q 6 P 4 Δ Q 4 .
Figure 1. Construction of auxiliary points Q i , i = 0 , , 7 , for a Class I sextic PH curve with its Bézier control points P i , i = 0 , , 6 . (a) Points Q 0 and Q 7 are on lines P 0 P 1 and P 5 P 6 , respectively, such that arg P 2 Q 0 Δ P 0 = arg Δ P 5 Q 7 P 4 and Q 0 P 2 P 4 Q 7 . (b) Points Q 1 and Q 6 are on lines P 0 P 1 and P 5 P 6 , respectively, such that the triangles Q 0 Q 1 P 2 and Q 7 P 4 Q 6 are similar; thus, we have Δ Q 0 Q 7 P 4 = P 2 Q 0 Δ Q 6 = P 2 Q 1 Q 6 P 4 . (c) Points Q 2 and Q 5 are on lines Q 1 P 2 and Q 6 P 4 , respectively, such that P 2 Q 1 2 ( Q 2 P 2 ) = Δ P 0 5 ( Q 1 P 1 ) and 2 ( P 4 Q 5 ) Q 6 P 4 = 5 ( P 5 Q 6 ) Δ P 5 . (d) The line Q 3 Q 4 is passing through the point P 3 , such that Q 3 Q 4 Q 0 P 2 P 4 Q 7 , and they further satisfy arg Δ Q 2 P 2 Q 1 = arg P 2 Q 1 Δ P 0 and arg Δ P 5 Q 6 P 4 = arg Q 6 P 4 Δ Q 4 .
Mathematics 11 00383 g001
Figure 2. A Class I sextic PH curve with its Bézier control points P i , i = 0 , , 6 forming a concave control polygon. Auxiliary points Q i , i = 0 , , 7 are computed to verify the PH property.
Figure 2. A Class I sextic PH curve with its Bézier control points P i , i = 0 , , 6 forming a concave control polygon. Auxiliary points Q i , i = 0 , , 7 are computed to verify the PH property.
Mathematics 11 00383 g002
Figure 3. Geometric characteristics of a Class II sextic PH curve. (a) Points Q 1 and Q 6 are on lines P 0 P 1 and P 5 P 6 , respectively, such that arg P 2 Q 1 Δ P 0 = arg Δ P 5 Q 6 P 4 and P 2 Q 1 P 4 Q 6 . (b) Points Q i , i = 2 , , 5 , satisfy that P 0 P 1 P 2 Q 2 P 3 Q 4 , Q 1 P 2 Q 2 Q 3 Q 4 Q 5 P 4 Q 6 , Q 3 P 3 Q 5 Q 4 P 5 P 6 , and 3 Δ P 0 : 5 ( Q 1 P 1 ) : 10 ( Q 2 P 2 ) : 30 ( Q 4 P 3 ) = 3 ( P 2 Q 1 ) : 2 Δ Q 2 : 2 Δ Q 4 : 3 ( Q 6 P 4 ) = 30 ( P 3 Q 3 ) : 10 ( P 4 Q 5 ) : 5 ( P 5 Q 6 ) : 3 Δ P 5 .
Figure 3. Geometric characteristics of a Class II sextic PH curve. (a) Points Q 1 and Q 6 are on lines P 0 P 1 and P 5 P 6 , respectively, such that arg P 2 Q 1 Δ P 0 = arg Δ P 5 Q 6 P 4 and P 2 Q 1 P 4 Q 6 . (b) Points Q i , i = 2 , , 5 , satisfy that P 0 P 1 P 2 Q 2 P 3 Q 4 , Q 1 P 2 Q 2 Q 3 Q 4 Q 5 P 4 Q 6 , Q 3 P 3 Q 5 Q 4 P 5 P 6 , and 3 Δ P 0 : 5 ( Q 1 P 1 ) : 10 ( Q 2 P 2 ) : 30 ( Q 4 P 3 ) = 3 ( P 2 Q 1 ) : 2 Δ Q 2 : 2 Δ Q 4 : 3 ( Q 6 P 4 ) = 30 ( P 3 Q 3 ) : 10 ( P 4 Q 5 ) : 5 ( P 5 Q 6 ) : 3 Δ P 5 .
Mathematics 11 00383 g003
Figure 4. A Class II sextic PH curve with its Bézier control points P i , i = 0 , , 6 forming a concave control polygon. Auxiliary points Q i , i = 1 , , 6 are computed to verify the PH property.
Figure 4. A Class II sextic PH curve with its Bézier control points P i , i = 0 , , 6 forming a concave control polygon. Auxiliary points Q i , i = 1 , , 6 are computed to verify the PH property.
Mathematics 11 00383 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, Y.; Fang, L.; Zheng, Z.; Cao, J. On Control Polygons of Planar Sextic Pythagorean Hodograph Curves. Mathematics 2023, 11, 383. https://doi.org/10.3390/math11020383

AMA Style

Li Y, Fang L, Zheng Z, Cao J. On Control Polygons of Planar Sextic Pythagorean Hodograph Curves. Mathematics. 2023; 11(2):383. https://doi.org/10.3390/math11020383

Chicago/Turabian Style

Li, Yujun, Lincong Fang, Zhihao Zheng, and Juan Cao. 2023. "On Control Polygons of Planar Sextic Pythagorean Hodograph Curves" Mathematics 11, no. 2: 383. https://doi.org/10.3390/math11020383

APA Style

Li, Y., Fang, L., Zheng, Z., & Cao, J. (2023). On Control Polygons of Planar Sextic Pythagorean Hodograph Curves. Mathematics, 11(2), 383. https://doi.org/10.3390/math11020383

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop