Approximate multi-degree reduction of Q-Bézier curves via generalized Bernstein polynomial functions

Quasi Bézier curves (or QB-curves, for short) possess the excellent geometric features of classical Bézier curves and also have good shape adjustability. In this paper, an algorithm for a multi-degree reduction of QB-curves based on L2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$L_{2}$\end{document} norm and by the analysis of geometric characteristics of QB-curves is constructed. The approximating approach for QB-curves of degree n+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$n+1$\end{document} by degree m (m≤n\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$m\leq n$\end{document}) is also given. Secondly, by solving the linear equations under the constraints of C0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$C^{0}$\end{document} and C1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$C^{1}$\end{document} and without constraints, the explicit expression of the points of the approximating curve is obtained, which minimizes the error between the original curve and the approximating curve using the least square method. Some numerical examples of degree reduction under different constraints are given, and the corresponding errors are calculated as well. The results show that this method can be easily implemented, is highly precise and very effective.


Introduction
A parametric curve, such as Bézier curve, is a fundamental tool and research content in the field of Computer Aided Geometric Design (CAGD)/Computer Aided Manufacturing (CAM). It plays a major role in geometric modeling of various products, shape designing, sketching, etc. It has also become one of the most significant schemes for representing curves in Computer Aided Design (CAD)/CAM systems due to straightforward, instinctive construction and several significant characteristics of a Bézier curve. In engineering applications and modeling the shape of a Bézier curve is not sufficient because it is uniquely determined by the control points. Researchers introduced a rational Bézier curve for modification or adjustment of the shape of curves without changing the control points with the help of weight factor to overcome this shortcoming. However, many problems occur when we solve these rational factors such as computational complexity, repeated derivation, inconvenient integration, and so on [1][2][3][4][5][6][7][8][9].
To defeat these problems and to get better approximation to real curves, researchers have introduced many Bézier curves in non-rational form (including trigonometric, poly-nomial, and hybrid trigonometric Bézier curves) with shape parameters [10][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29]. In [25], generalized Bernstein basis functions were used to construct a QB-curve with multiple shape parameters. The curve of this kind also possesses several significant geometric characteristics of classical Bézier curves, and it also has bendable shape modification, i.e., the shape of a curve can be adjusted by altering the values of shape parameters to construct more complicate curves of the same degree. This makes QB-curves extensively used in different modelings and CAD/CAM systems and has certain application value in describing curves and surfaces. Hu et al. presented the continuity conditions of the smooth splicing of two adjacent QB-curves G1 and G2 and analyzed the result of shape parameters on the shape of the combined curves after splicing [30]. In [31], the continuity constraints for QB surfaces of degree (m, n) are studied by using the end point properties of generalized QB basis functions. As far as we are aware, the research on a QB curve regarding its degree reduction has not been reported before.
For the degree of curves, different CAD/CAM systems have different requirements. It is necessary to reduce the degree of curves to realize the data conversion and transmission between curves of different degree. There are two types of methods for degree reduction of curves: (a) first is the geometric method which is based on control points. The inverse process of degree elevation is used to solve the control points of degree reduction curves in [32]. Young et al. used the geometric properties of a Bézier curve and combined them with generalized inverse matrix and least squares theory to achieve degree reduction approximation [33,34]. (b) The second type is the algebraic method which is based on basis transformation. The degree reduction of a Bézier curve is achieved by Chebyshev-basis transformation [35,36], and the best approximation degree reduction problem is studied by using Legendre polynomial theory [37]. Xu et al. presented a method for degree reduction of a Bézier curve based on constrained Jacobi polynomials [38]. Li et al. [39] constructed some geometric continuity conditions for the generalized cubic H-Bézier model for the purpose of constructing shape-controlled complex curves and surfaces in engineering. The authors in [40,41] constructed generalized trigonometric Bézier curves with shape parameters for the purpose of constructing some complex curves and surface applications in engineering. The basis functions proposed in [39][40][41] are different from the basis functions utilized in this study.
In this paper, based on L 2 -norm, the least square of mth-degree QB-curves is used to approximate n + 1th-degree QB-curves. The degree reduction without constraints and under C 0 , and C 1 constraints is considered. The specific expression for calculating the degree reduction of curves control points is given.

The definition of QB-curves
Let P i ∈ R d , i = 0, 1, . . . , n and d = 2, 3, be the set of control points, and its corresponding polynomial curve with shape parameters {λ i } n i=1 of degree n is called QB curve, which can be defined as follows: where nth-degree basis function {b 0,n (θ ), b 1,n (θ ), . . . , b n,n (θ )} is the nth-degree QB basis functions with shape parameters {λ i } n i=1 , the specific form is [25] Here

Problem description
Problem 1 Let us consider that the n + 1th-degree QB curve determined by control points {P * i } n+1 i=0 ∈ R d , d = 2, 3, has the following form: where {b i,n+1 (θ )} n+1 i=0 are n + 1th degree QB basis functions. The degree reduction of the so-called QB curve refers to finding the low mth-degree QB curve (m ≤ n) whose control such that the distance between the two curves is minimized under a certain distance function d(r * (θ ), r(θ )).
In order to obtain the explicit expression of an approximating QB curve, we choose to use the L 2 -norm to measure the degree of the approximating curve before and after reduction as a whole, and define its "distance" as follows: The above equation is a vector function, and the labels r * (θ ) = (r * 1 (θ ), r * 2 (θ ), . . . , r * s (θ )) and r(θ ) = (r 1 (θ ), r 2 (θ ), . . . , r s (θ )) are introduced to transform the vector function into the minimized component function Thus, Eq. (5) can be determined by Eq. (7) Therefore, when each component distance function d(r * j (θ ), r j (θ )) gets the minimum value, d(r * (θ ), r(θ )) reaches the minimum value. In this paper, we only discuss the minimum problem in the form of component function and introduce the problem of degree reduction from the solution problem.
i=0 be the sequence of real numbers, and its corresponding n + 1thdegree QB function can be defined as follows: Then we will seek real numbers {P i } m i=0 corresponding to the mth-degree QB function (m ≤ n) such that d 2 (r * (θ ), r(θ )) = 1 0 [r * (θ ) -r(θ )] 2 dθ minimizes by least square distance.
In order to determine the coefficients of the approximating function r(θ ), the next main purpose is to solve the unknowns {P i } m i=0 .

Degree reduction of QB-curves by C 0 continuity conditions
To reduce the degree of the curves, if the C 0 continuity is satisfied, that is, the first and last points of the two curves coincide, then there are The remaining m -1 points are calculated according to Theorem 2.

Degree reduction of QB-curves by C 1 continuity conditions
To reduce the degree of the curves, if the C 1 continuity is required, that is, the first and last points of two curves have C 1 continuity interpolation, then there are P 0 = P * 0 , P m = P * n+1 , P 1 = P * 0 + n + 1 + λ 1 m + λ 1 P * 1 -P * 0 , P m-1 = P * n+1 - The remaining m -3 points are calculated according to Theorem 3.

Examples of degree reduction curves
In this paper, a wide numerical study has been carried out for the verification of correctness of the algorithm. The following is a numerical example of the application of the algorithm to QB-curves, in which the square distance formula is used to determine the error between the curves before and after reduction.

Constraint condition Control points Error
Under unrestricted condition to construct eight QB-curves (blue solid lines) without constraints and under the constraint of C 0 and C 1 , the curve is reduced to quintic QB-curves (red dashed lines).
Here we give two different shape parameters to the reduced quintic QB curve. That is, λ 1 = 1, λ 2 = 0, λ 3 = 0, λ 4 = 1, λ 5 = 0; The curves before and after reduction are shown in Fig. 2, and the control points and errors after reduction are shown in Tables 3 and 4, respectively.     to construct eight QB-curves (blue solid lines) without constraints and under the constraint of C 0 and C 1 , the curve is reduced to quartic QB-curves (red dashed lines). Here we give two different shape parameters to the reduced quartic QB curve. That is,

Example 5.3 Given the shape parameters
The curves before and after reduction are shown in Fig. 3, and the control points and errors after reduction are shown in Tables 5 and 6, respectively.

Conclusions
In this paper, the least square degree reduction approximation problem for QB-curves based on L 2 -norm without constrains and under the C 0 and C 1 constraints is studied. An algorithm for control points of approximating curves is also given. Three practical examples and their specific errors under three conditions reveal that the method achieves one-time reduction and multi-degree least square approximation of QB curve under various constraints. That is to say, this method is applicable for the system of CAD/CAM modeling. The degree reduction for QB surfaces will be studied in future work.