ϵκ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon \kappa $$\end{document}-Curves: controlled local curvature extrema

The κ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa $$\end{document}-curve is a recently published interpolating spline which consists of quadratic Bézier segments passing through input points at the loci of local curvature extrema. We extend this representation to control the magnitudes of local maximum curvature in a new scheme called extended- or ϵκ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon \kappa $$\end{document}-curves.κ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa $$\end{document}-curves have been implemented as the curvature tool in Adobe Illustrator® and Photoshop® and are highly valued by professional designers. However, because of the limited degrees of freedom of quadratic Bézier curves, it provides no control over the curvature distribution. We propose new methods that enable the modification of local curvature at the interpolation points by degree elevation of the Bernstein basis as well as application of generalized trigonometric basis functions. By using ϵκ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon \kappa $$\end{document}-curves, designers acquire much more ability to produce a variety of expressions, as illustrated by our examples.


Introduction
The κ-curve, proposed recently by [28], is an interpolating spline which is curvature-continuous almost everywhere and passes through input points at the local curvature extrema.It has been implemented as the curvature tool in Adobe Illustrator ® and Photoshop ® and is accepted as a favored curve design tool by many designers (see, e.g., [4,6]).
We consider the reasons for the success of κ-curve to be: 1. Information along contours is concentrated at local maxima of curvature.2. Curves of low degree have smooth distribution of curvature.3. G 2 -continuous curves tend to look fairer than only G 1continuous ones.
Attneave [1] suggested, based on his empirical study, that information along contours is concentrated in regions of high magnitude of curvature, as opposed to being distributed uniformly along the contour, and it is further concentrated at local maxima of curvature (see also [33]).Although Attneave never published the details of his methods, [20] conducted a similar experiment and obtained the same results.Levien and Séquin [12] argue similarly and assert that points of maximal curvature are salient features.
The curvature of a polynomial curve is given by a relatively complicated rational function [7], and its distribution might not be globally smooth.However, if the curve is of a low degree, the curvature distribution is more uniform and the curve is fairer, thus more suitable for illustration.The quadratic polynomial curve has the nice property that its curvature has only one local maximum, and its location is easily computable [28], which makes the handling of curvature extrema much easier.
Graphic designers often accept G 1 continuity as good enough for illustration.However, discontinuity remains; for example, if you join a straight line and a circular arc with G 1 continuity, the rhythm of the curve will be broken at the joint.For this reason, we give preference to G 2 -continuous curves.
Nonetheless, κ-curves are not perfect, and their further investigation is necessary [29].The following are two important shortcomings of κ-curves: 1.They are not curvature-continuous everywhere: at inflection points only G 1 continuity is guaranteed.2. Since the degree of freedom (DoF) of the quadratic segments is limited, it is impossible to control the magnitudes of local maximum curvature at the input points.
For the first shortcoming, Wang et al. [26] provided a solution through the use of log-aesthetic curves [15,17] instead of polynomial curves.Log-aesthetic curves have a shape parameter (α), which can be utilized to control the curvature distribution as shown in Fig. 1.Their method guarantees G 2 continuity everywhere, including inflection points, since log-aesthetic curves with negative α values can represent Sshaped curves with G 2 continuity.(Note that these cannot be represented by quadratic Bézier curves.) These curves, however, are defined by a Cesàro equation and thus take extra time to evaluate, making the interpolation method impractical for real-time design purposes.
Because of the second shortcoming, if the designer wants to increase or decrease the magnitude of the curvature extremum, she needs to add extra input points, as shown in Fig. 2. Yan et al. [27] proposed a piecewise rational, quadratic, interpolatory curve that is able to reproduce circles and other elliptical or hyperbolic shapes.Although their main intention was to reproduce circles, their method could also control the magnitude of local maximum curvature.However, only rational quadratic curves are applicable, and it is not possible to extend their method for other types.
In this paper, we propose a new method to solve the second shortcoming by degree elevation of the Bernstein basis functions, giving an extra DoF to each quadratic curve segment, Fig. 1 Log-aesthetic curves with various α values [26].Input points are depicted by black boxes, and green points correspond to positions of local maximum curvature.The blue curves show the normal curvature.These curves are G 2 -continuous everywhere Fig. 2 Addition of extra input points to control the magnitude of curvature extrema: b is the original κ-curve, a, c shows deformed curves with large and small curvature extrema by adding two extra input points providing control over the magnitudes of local maximum curvature at input points as shown in Fig. 3.In order to increase the designers' possible choices, we also introduce a new trigonometric basis for which we can perform degree elevation.In addition, we propose a general method for bases with extra shape parameters.By adding one more parameter to each of the curve segments, the designers obtain more expressive power for their illustration.The family of this new curve is denoted as κ-curves.κ-curves preserve all of the appealing properties of κcurves, i.e., point interpolation, G 2 continuity (except at inflection points), continuous modification (changes smoothly when the input points move), local influence, and real-time generation.With a small processing overhead, κ-curves offer the ability to control the magnitude of local maximum curvature.
We have implemented κ-curves in MATLAB ® and Julia [11].The source of the Julia code is available online [22].
The rest of this paper is organized as follows.Section 2 reviews the related work.Section 3 presents a method to control the magnitudes of local maximum curvature by degree elevation of the Bernstein basis functions.Section 4 introduces a new trigonometric basis and proposes a method with degree elevation similar to the one proposed in the previous section.Finally, we end with conclusions and discussion of future work.
(a=2/3) a=0.7 a=0.8 (a=2/3) a=0.95 Fig. 3 The left-and rightmost curves are κ-curves; the others are κ-curves with gradual changes in the global shape parameter a.The face part (for a = 0.95 and a = 2/3) is zoomed in for better comparison.When a=2/3, κand κ-curves are identical

Related work
In this section, we first review [28] and their underlying strategy developing κ-curves.Next, we discuss related researches on various kinds of basis function formulations for curve design.

Ä-Curve
The basic framework of our method is adopted from [28], generating curves controlled by interpolation points.κcurves have stimulated the field of interpolatory curve generation, resulting in works such as [5], which proposes a method for good control over the location and type of geometric feature points (e.g., cusps and loops).
Yan et al. [28] create a sequence of quadratic curves with G 2 continuity almost everywhere.They derived explicit formulae for the point between two quadratic segments to guarantee G 2 continuity and for the additional condition that input points should be interpolated at maximum curvature magnitude positions.κ-curves are determined by the locations of the middle control points of quadratic Bézier segments-the rest is easily derived from the continuity constraints.Hence, the variables are the locations of these middle control points.
Their basic strategy to determine these locations is to adopt a local/global approach [13,23]: G 2 -continuous connection is performed locally, while the interpolation at maximum curvature magnitude positions is done globally.
Regarding control of the magnitude of local curvature, the most common technique is to change the weight of a control point of a rational curve [7].A larger weight attracts the curve to its control point, which makes local curvature larger.However, this technique is not applicable for interpolatory curves.Another technique is to introduce extra parameters called bias and tension to the B-spline formulation for controlling local curvature [2], but this is also not applicable to interpolatory curves.

Basis functions
As mentioned in [28], curve modeling has a long history, especially in computer-aided geometric design, as well as computer graphics.In CAGD and applied mathematics, to extend the expressive power of curves, many researchers have been trying to develop new bases with extra shape parameters.The following is a (nonexhaustive) list of such bases: 1. C-Bézier curve [31] 2. Cubic alternative curve [10] 3. Cubic trigonometric Bézier curve (T-Bézier basis) [8] 4. αβ-Bernstein-like basis [32] 5. Quasi-cubic trigonometric Bernstein basis [30] 6. Trigonometric cubic Bernstein-like basis [24] Our method proposed in the next section can be applied for curves based not only on polynomials, but also other bases such as a trigonometric basis with degree elevation property which we will introduce in Sect. 4. All of the representations listed above have extra parameters for shape control, and we can utilize these parameters to control the magnitudes of local maximum curvature.(Not all curve types are applicable, however, as explained in Appendix C.) It is interesting to note that most researchers to date have attempted to develop new bases using four control points, based on cubic polynomials or quadratic trigonometric functions.They prefer using four control points out of concern for connections at both ends of the curve.In order to control the magnitude of curvature at the two ends independently, at least two control points are necessary at each end.This differs fundamentally from Yan et al.'s (and our) approach, which uses only three control points.
The importance of [28] is their proposed paradigm shift for curve generation, by considering the local maximum curvature in the middle, instead of focusing on the endpoints.If we can assume that the curvature has just one local maximum in each segment, then only one extra parameter per segment is adequate to control the magnitude of its local maximum curvature.
To our best knowledge, no trigonometric basis family for arbitrary degree has been published yet.Our novel generalized trigonometric basis functions range from linear, using three control points, to any higher degree n, using 2n + 1 control points.The curve can be evaluated by a recursive method , similar to de Casteljau's algorithm, as explained in Appendix B. Since the curve uses trigonometric functions as blending functions, it can represent a circular arc exactly, without using a rational form.

Cubic Bernstein polynomials
In this section, we extend κ-curves in a direct manner, by elevating the degree of quadratic Bézier segments to cubic.
2. All local maximum curvature points are the same as the input points.3. G 2 continuity is guaranteed almost everywhere (except for inflection points).
In the following, we discuss only closed curves, but it is straightforward to extend our methods to open curves as has been demonstrated for κ-curves.
If we elevate the degree of a planar Bézier curve, we obtain an additional control point, which has two DoFs (the x and y coordinates).To reduce these to one, we add a geometric constraint on the location of the second and third control points of the cubic Bézier curve, as shown in Fig. 4. Here, a is an internal division ratio, where the larger a is, the closer the control points P 1 and P 2 are to the control point Q 1 .We make the restriction 2/3 ≤ a < 1 because the curve should not have a complicated curvature distribution.Using Q i , the curve C(t; a) is expressed by Note that if a = 2/3, the curve degenerates to quadratic.
We have proved that by constraining the construction of the cubic polynomial curve as in Eq. (1), using only three control points instead of four, the curvature in one curve segment has at most one local maximum for 2/3 ≤ a < 1; see details in [16], as well as Appendix A for a general discussion on the curvature extrema of cubic polynomial curves, and a high-level summary of the proof.In addition, we have made MATHEMATICA simulation available in [9] to compute the number of curvature extrema for this curve using a : 1 -a P 0 =Q 0 Fig. 4 Constrained cubic Bézier curve.If a = 2/3, the curve becomes quadratic both classical approach and Sturm's theorem.Hence, we can safely assume that the curvature in one curve segment has at most one local maximum, and a single extra parameter for each segment is enough to control the magnitudes of local maximum curvature.

Geometric constraints
We assume that κ-curves consist of a sequence of constrained cubic polynomial curves parameterized by t and also a i , which is an extra shape parameter.The control points are given by c i,0 , c i,1 and c i,2 ∈ R 2 , corresponding to Q i , i = 0, 1, 2 in Fig. 4. The a i 's are reserved for designers and can be manipulated independently.
The curvature of this curve c i (t; a i ) is given by where indicates the area of the triangle specified by its arguments, and In the quadratic case (i.e., κ-curves), as the curvature has such a simple formula, we can express the parameter t i at the point of maximal curvature explicitly, in terms of the Bézier coefficients of the i th quadratic Bézier curve as where a, b means the scalar product of vectors a and b.Then, we add the condition where p i is the i th input point.Solving for c i,1 and substituting into Eq.( 4), we get a cubic equation in t i that depends only on the endpoints c i,0 and c i,2 , and the input points p i .Unfortunately, we cannot obtain an explicit formula like Eq. (4) for the parameter t i in the cubic case, because of its high degree (see details in Appendix A.1), but this is not a problem.Solving Eq. ( 5) for c i,1 , we arrive at Substituting this into the derivative of Eq. ( 3), and letting it equal 0, we obtain (after some simplification) a polynomial equation of degree 9 in t i .This equation can be derived by the Maxima [14] code in Fig. 13 (Appendix A).We solve this equation and select a real root in [0, 1].Note that we have proved that there is one and only one solution for the polynomial equation of degree 9 in t i ∈ [0, 1] as in the case of κ-curves; see details in [16], as well as Appendix A.2.However, when we use other types of curves with more complicated representations (see examples in Appendix C), this kind of formula may be hard to derive.In these cases, we can use the relaxed Newton's method to compute the maximum curvature.For this, we need to be able to compute the curvature and its derivative; we do this using a quadratic Taylor series approximation around the last value of t i .
We introduce the constant λ i (0 < λ i < 1) according to the construction method of κ-curves and set Let the curvatures at the endpoints of the curve segment be denoted by κ i (0; a i ) and κ i (1; a i ), then from Eq. ( 3) introducing the notations and By adopting the local/global approach, we treat c i,0 as fixed for the computation of κ i (1; a i ), although it depends on λ i−1 (similarly for c i+1,2 ).
In order to guarantee G 2 continuity at the joint of two consecutive segments, the following equations should be satisfied: Hence, Since 0 < a i , a i+1 < 1, λ i is real and 0 < λ i < 1.

Optimization
In the global phase, we calculate the positions of the middle control points c i,1 by solving a linear system of equations.We treat the current values of λ i (internal division ratios of c i,1 and c i+1,1 ) and t i (parameters of local maximum curvature) as fixed.Substituting Eq. ( 7) into Eq.( 5), we get which can be solved for c i,1 .
The optimization process is summarized in Algorithm 1.

Algorithm 1:
The optimization process.
Result: control points c i,k and parameters of maximal curvature t i Set all λ i to 0.5; Compute all c i,0 and c i,2 by Eq. ( 7); while not convergent do Compute all λ i by Eq. ( 10); Compute all c i,0 and c i,2 by Eq. ( 7); Compute all t i by polynomial root finding; Compute all c i,1 by Eq. ( 11); end Compute all c i,0 and c i,2 by Eq. ( 7);

Results
Figure 5 shows examples of closed κ-curves along with the original κ-curve.The input points are located at the same positions.The a i values of these curves are equal to 2/3, except for one, two and three input points, respectively, where a i is set to 0.85.If the a i of all input points are 2/3, the κcurve on the left is generated.Since we specify a larger value for some input points, the magnitudes of the corresponding local maximum curvature increase, as we expected.Fig. 7 The leftmost curve is a κ-curve; the other curves are κ-curves: a is equal to 0.75, 0.85 and 0.95, respectively.Note the curvature at the inflection points Figure 6 shows another example of local curvature control.From the left to right, the first drawing shows a bird using κ-curves.In the second, we set a at the wing tips to 0.95, while leaving all others at the default 2/3.This has the effect of sharpening the wing tips.In the third figure, we reversed the role of the input points, giving a = 2/3 to those at the wing tips and 0.95 to all other points.Here, the wings are rounded, while other parts of the bird get sharper.
Notice that the bird's beak resembles a cusp, but is actually the start and end points of an open curve located at the same position.In our implementation, we limit 2/3 ≤ a ≤ 1 to make a curve with smooth curvature distribution, which disallows the generation of a cusp even at a = 1.In cases where the designer wants to use a cusp, the curve should be cut in two, or the input points should be relocated to form a cusp as explained in [28].
Figure 7 shows examples of global curvature control.There are three κ-curves, with a set to 0.75, 0.85 and 0.95, respectively, along with the original κ-curve (a = 2/3) for comparison.By increasing a, the magnitudes of local maximum curvature increase.As the close-up windows indicate, at the inflection point G 2 continuity is violated for κ-curves, and only G 1 continuity is guaranteed.However, for κ-curves with a larger a, the magnitude of curvature at inflection points, and consequently the G 2 error, becomes smaller.Note that although these curves are almost G 2 -continuous everywhere, they are quite different from those in Fig. 1 (generated using the same input points).
Figures 3, 8 and 9 show the effect of changes of the global shape parameter a on various designs.As discussed above, larger a values generally induce larger local curvature extrema and steeper curvature variation.The resulting curves look more sharp at the input points and more flat between them.

Generalized trigonometric basis
In this section, we describe our new generalized trigonometric basis.This is based on the trigonometric cubic Bernstein-like basis [24], which we are going to review first.The trigonometric cubic Bernstein-like basis functions have an extra shape parameter α and are defined by (12) where S = sin π t 2 , C = cos π t 2 , for α ∈ (0, 2), t ∈ [0, 1].Note that these functions satisfy partition of unity, i.e., 3 i=0 f i (t) = 1 for any α.When α = 1, the above functions are simplified to If we add the second and third functions together and rename them to u, v and w, we obtain blending functions {u, v, w} as follows: It is straightforward to define a curve by these blending functions with three control points, which we can regard as a "linear" trigonometric curve since the highest degree the trigonometric functions are in is one.
One interesting relationship among these functions is which enables and yields the five blending functions {u 2 , 2uv, 4uw, 2vw, w 2 }, associated with five control points.We can define a curve using these blending functions and regard it as a "quadratic" trigonometric curve since the highest power of each blending function is now degree two.In a similar way, we can extend blending functions of "degree" n with 2n + 1 control points.As explained in Appendix B, we can perform a recursive procedure to evaluate a curve of any degree similar to de Casteljau's algorithm avoiding the overhead of trigonometric function evaluation.This means that it is not necessary to calculate the coefficients of blending functions, or keep a coefficient table.
We formulate the κ-curve in this basis using a strategy similar to that in the previous section, i.e., using a sequence of quadratic trigonometric curves with a constraint on the positions of their control points, as shown in Fig. 10.Note the location of the control point 123 a : 1 -a 1a :a 1:1 Fig. 10 Constrained quadratic trigonometric curve.When a = 1/2 (bottom), the curve becomes linear When a is equal to 1/2, the curve degenerates to a linear trigonometric curve.

Geometric constraints and optimization
First, we analyze the linear trigonometric curve since it corresponds to the original κ-curve.Let c i (t) be a linear trigonometric curve with control points c i,0 , c i,1 and c i,2 and defined by where S = sin π t 2 , C = cos π t 2 and t ∈ [0, 1].Its curvature is given by where The numerator of the above formula does not depend on t, so the extrema of the following f i (t) corresponds to those of κ i (t): and its derivative with respect to t is given by By assuming d f i (t)/dt = 0 with S, C = 0, we obtain where γ = ( s i 2 − r i 2 )/ r i , s i .We can solve the above equation and obtain Since 0 ≤ S, C ≤ 1, we have the unique solution Hence, Note that when r i and s i are perpendicular to each other, if r i = s i , then the curve becomes a circular arc, and no local maximum curvature exists.If r i > s i , then the curvature at t = 1 will be maximum, and if r i < s i , the curvature at t = 0 will be maximum in this curve segment.
For a quadratic trigonometric curve, the curvatures κ i (1; a i ) and κ i+1 (0; a i ) at the endpoints of the constrained quadratic trigonometric curve c i are given by We can calculate λ i by guaranteeing G 2 continuity at the joint of c i (1; a i ) and c i+1 (0; a i+1 ): As before, we get a linear system of equations for c i,1 : a=0.75 a=0.9 a=0.75 a=0.9 a=0.55 a=0.6 a=2/3 Fig. 11 The two curves on the top left (red) are κ-curves; the two on the top right (brown) are κ-curves using cubic Bernstein basis functions with a = 0.75 and 0.9.The bottom row shows κ-curves using quadratic trigonometric basis functions with a = 0.55, 0.6, 0.75 and 0.9 a=0.7 a=0.55 Fig. 12 A Christmas tree drawn with κ-curves using the cubic Bernstein basis functions (left) and the quadratic trigonometric basis functions (right).Note that the latter has more rounded forms and-in this case-preferable where

Results
Figures 11 and 12 show examples of κ-curves using the quadratic trigonometric basis functions explained in this section.In the first figure, the top left two curves (red) are κ-curves.The top right two curves (brown) are κ-curves, using cubic Bernstein basis functions with a = 0.75 and 0.9.The bottom row (green) shows κ-curves using quadratic trigonometric basis functions with a = 0.55, 0.6, 0.75 and 0.9.The curves in the bottom row are more rounded than those of the Bernstein basis.By increasing a, the differences between the two types of curves become smaller as they approach a polyline generated by connecting the input points.
The second figure shows a case where these more rounded forms are clearly preferable.

Use of built-in shape parameters
So far our strategy made use of degree elevation-from quadratic to cubic in the polynomial case, and from linear to quadratic in the trigonometric case.Yet another strategy is to use extra shape parameters built into the basis functions.
As an illustrative example, take the trigonometric cubic Bernstein-like basis functions reviewed in the previous section.In the framework of κ-curves, we have to degenerate the curve by relocating the positions of its control points, essentially reducing its degree.
The trigonometric cubic Bernstein-like basis functions need four control points to define a curve.To construct "quadratic" curves corresponding to the quadratic Bézier segments of κ-curves, we make the second and third control points collocate.Hence, the blending functions become where S = sin π t 2 , C = cos π t 2 , for α ∈ (0, 2), t ∈ [0, 1].Incidentally, this will result in the same curve as Eq. ( 17), if we substitute 2a for α.
Note that several curve types-such as the cubic alternative curve and the αβ-Bernstein-like basis functions, for specific extra parameters-have zero curvature at the endpoints, so we cannot obtain the internal division ratio λ i .For these curves, the method shown in this section cannot be applied.See details in Appendix C.

Conclusions and future work
We have proposed two types of κ-curves as extensions of κ-curves, for controlling the magnitudes of local maximum curvature.Our methods use degree elevation of the Bernstein basis functions and a new family of trigonometric basis functions.
In line with Yan et al.'s paradigm shift for curve generation, we consider the local maximum curvature in the middle, instead of focusing on the endpoints.The new curves preserve all of the nice properties of κ-curves, i.e., point interpolation, G 2 continuity (except at inflection points), continuous modification (changes smoothly when the input points move) and local influence.Computing κ-curves is quite fast, and designers can manipulate them interactively, acquiring much more expressive power for curve design, as illustrated by our examples.The processing time for generating κ-curves is similar to that of κ-curves, especially for κ-curves with cubic Bernstein basis.For example, κ-curves consume 0.07 sec to draw Fig. 8 (bear).Under similar conditions, cubic Bézier κ-curves consume 0.08 sec and generalized trigonometric κ-curves takes 0.37 sec on average.
Future work includes the development of κ-curve plugins for Adobe Illustrator ® and Photoshop ® .Another possible research direction is to apply the proposed method to different types of aesthetic curves, e.g., log-aesthetic curves [26], σcurves [18] or τ -curves [19].
where x y − x y = 0 means the curvature is equal to 0, and it corresponds to inflection points.Consequently, (x y − x y )(x 2 + y 2 ) − 3(x x + y y )(x y − x y ) = 0 (33) corresponds to curvature extrema.The above equation is regarded as (quadratic × quartic) − (cubic × cubic), but the coefficient of sextic terms vanishes and it becomes quintic.Therefore, a cubic polynomial curve can have at most 5 curvature extrema.Please refer to [25] for more details.Since Eq. ( 33) is quintic, it is generally not possible to have analytical solutions, so we must use numerical approach to obtain the points on the curve where curvature extrema occurs.

A.2 Uniqueness
Here, we give a high-level summary of the proof [16] that the curvature of a cubic curve of the form shown in Eq. ( 1) has at most one local extremum in the (0, 1) parameter interval.As a consequence, the degree-9 polynomial in Sect.3.1, computed by the Maxima [14] program in Fig. 13, has at most one real solution in the (0, 1) interval.When it has none, this means that no local extremum is present in the curvature, so the extremum occurs at Without loss of generality, place the control points as follows: where b ≥ 0 and h > 0. (The special cases where Q 0 = Q 2 or h = 0 are also handled in [16].)Let N (t, a) denote the left side of Eq. ( 33) applied to this curve.When b ≤ 3−2/a, it is easy to see that ∂ N (t, a)/∂t < 0 for any t ∈ (0, 1) and a ∈ [2/3, 1], so N is decreasing.Since N (0, a) is always positive, this means N has at most one 0-crossing.
In the following, let us also assume that b > 3 − 2/a.Then, the following statements can also be proven: From the above" it is easy to prove that N (t, a) = 0 has exactly one solution-and thus, the curvature has at most one local extremum-in (0, 1).

B.1 Recursive evaluation
For our new trigonometric basis, we can derive a recursive algorithm similar to de Casteljau's algorithm .For simplicity we explain only the quadratic case, but it can be extended to a general degree n by induction.To shorten expressions, we use u = 1 − S(t), v = S(t) + C(t) − 1 and w = 1 − C(t), where S(t) = sin π t 2 and C(t) = cos π t 2 .Note that v 2 = 2uw, and For a quadratic curve with this basis, five control points P i (i = 0 . . .4) are used, and the curve point at t is evaluated as Hence, the algorithm repeats a simple blending of three points u P i−1 + v P i + w P i+1 to generate a point on the given curve.

B.2 Triangle method
We can also construct a triangle using the coefficients of trigonometric basis functions, similarly to Pascal's triangle.
Below is a table of degree elevation, from the first row representing degree 1 to the sixth row representing degree 6:

C Various basis functions
Here, we check the applicability of the bases listed in Sect.2, except for the trigonometric cubic Bernstein-like basis functions [24], since that was already discussed in the paper.

C.1 C-Bézier curve [31]
The basis of the C-Bézier curve is {sin t, cos t, t, 1}, and the curve is defined by the following formula: 123 Fig. 13 Maxima [14] Here, α is a built-in shape parameter satisfying 0 < α ≤ π , and S = sin α, C = cos α.The parameter of the curve is t ∈ [0, α], and We degenerate the curve by adding its second and third basis functions, i.e., placing the second and third control points at the same position.The curve is then defined by three control points.However, even if we vary α from 0 to π , the blending functions do not vary much, as shown in Fig. 14

C.2 Cubic alternative curve [10]
The cubic alternative curve is similar to the cubic Bézier curve, and is defined by When α = β = 2, the curve becomes the cubic Ball curve; for α = β = 3, the classical cubic Bézier curve; and for α = β = 4, the cubic Timmer curve.The basis functions are nonnegative when 0 ≤ α ≤ 3.
When α = 3, this degenerates to the classical cubic Bézier curve with collocated second and third control points, and at the endpoints we get κ i (0, 3) = κ i (1, 3) = 0, since the directions of the first and second derivatives are the same.Hence, the proposed method is not applicable in this case.
By assuming a local maximum curvature at t i , we can express the input point p i as

C.3 Cubic trigonometric Bézier curve [8]
For where S = sin π t 2 and C = cos π t 2 .We assume λ = μ = α, and the curve is degenerated by adding the second and third blending functions.Hence, We define the curve c(t; 123

Fig. 5 Fig. 6
Fig.5 The leftmost curve is a κ-curve.The other curves are κ-curves: the a i values are equal to 2/3 except for one, two and three input points, respectively, where a i = 0.85