Composite Bézier curve

File:Beziergon.svg

Béziergon – The red béziergon passes through the blue vertices, the green points are control points that determine the shape of the connecting Bézier curves

In geometric modelling and in computer graphics, a composite Bézier curve or Bézier spline is a spline made out of Bézier curves that is at least $C^0$ continuous. In other words, a composite Bézier curve is a series of Bézier curves joined end to end where the last point of one curve coincides with the starting point of the next curve. Depending on the application, additional smoothness requirements (such as $C^1$ or $C^2$ continuity) may be added.

Perhaps the most common use of composite Béziers is to describe the outline of each letter in a PostScript or PDF file. Such outlines are composed of one béziergon for open letters, or multiple béziergons for closed letters. Modern vector graphics and computer font systems like PostScript, Asymptote, Metafont, OpenType, and SVG use composite Bézier curves composed of cubic Bézier curves (3rd order curves) for drawing curved shapes.

Terminology

File:Sinc Function Approximation with Bezier Splines.svg

Sinc function approximated using a smooth Bézier spline, i.e., a series of smoothly-joined Bézier curves

A continuous composite Bézier is also called a polybézier, by similarity to polyline, but whereas in polylines the points are connected by straight lines, in a polybézier the points are connected by Bézier curves. A béziergon (also called bézigon) is a closed path composed of Bézier curves. It is similar to a polygon in that it connects a set of vertices by lines, but whereas in polygons the vertices are connected by straight lines, in a béziergon the vertices are connected by Bézier curves. Some authors even call a $C^0$ composite Bézier curve a "Bézier spline"; the latter term is however used by other authors as a synonym for the (non-composite) Bézier curve, and they add "composite" in front of "Bézier spline" to denote the composite case.

Smooth joining

A commonly desired property of splines is for them to join their individual curves together with a specified level of parametric or geometric continuity. While individual curves in the spline are fully $C^{\mathrm{\infty }}$ continuous within their own interval, there is always some amount of discontinuity where different curves meet.

The Bézier spline is fairly unique in that it's one of the few splines that doesn't guarantee any higher degree of continuity than $C^0$. It is, however, possible to arrange control points to guarantee various levels of continuity across joins, though this can come at a loss of local control if the constraint is too strict for the given degree of the Bézier spline.

Smoothly joining cubic Béziers

Given two cubic Bézier curves with control points $[{\mathbf{𝐏}}_0,{\mathbf{𝐏}}_1,{\mathbf{𝐏}}_2,{\mathbf{𝐏}}_3]$ and $[{\mathbf{𝐏}}_3,{\mathbf{𝐏}}_4,{\mathbf{𝐏}}_5,{\mathbf{𝐏}}_6]$ respectively, the constraints for ensuring continuity at ${\mathbf{𝐏}}_3$ can be defined as follows:

• $C^0/G^0$ (positional continuity) requires that they meet at the same point, which all polybéziers do by definition. In this example, the shared point is ${\mathbf{𝐏}}_3$
• $C^1$ (velocity continuity) requires the neighboring control points around the join to be mirrors of each other. In other words, they must follow the constraint of ${\mathbf{𝐏}}_4=2{\mathbf{𝐏}}_3-{\mathbf{𝐏}}_2$
• $G^1$ (tangent continuity) requires the neighboring control points to be collinear with the join. This is less strict than $C^1$ continuity, leaving an extra degree of freedom which can be parameterized using a scalar ${\beta }_1$. The constraint can then be expressed by ${\mathbf{𝐏}}_4={\mathbf{𝐏}}_3+({\mathbf{𝐏}}_3-{\mathbf{𝐏}}_2){\beta }_1$

While the following continuity constraints are possible, they are rarely used with cubic polybéziers, as other splines like the B-spline or the β-spline will naturally handle higher constraints without loss of local control.

• $C^2$ (acceleration continuity) is constrained by ${\mathbf{𝐏}}_5={\mathbf{𝐏}}_1+4({\mathbf{𝐏}}_3-{\mathbf{𝐏}}_2)$. However, applying this constraint across an entire cubic Bézier spline will cause a cascading loss of local control over the tangent points. The curve will still pass through every third point in the spline, but control over its shape will be lost. In order to achieve $C^2$ continuity using cubic curves, it's recommended to use a cubic uniform B-spline instead, as it ensures $C^2$ continuity without loss of local control, at the expense of no longer being guaranteed to pass through specific points
• $G^2$ (curvature continuity) is constrained by ${\mathbf{𝐏}}_5={\mathbf{𝐏}}_3+({\mathbf{𝐏}}_3-{\mathbf{𝐏}}_2)(2{\beta }_1+{\beta }_1^2+{\beta }_2/2)+({\mathbf{𝐏}}_1-{\mathbf{𝐏}}_2){\beta }_1^2$, leaving two degrees of freedom compared to $C^2$, in the form of two scalars ${\beta }_1$ and ${\beta }_2$. Higher degrees of geometric continuity is possible, though they get increasingly complex
• $C^3$ (jolt continuity) is constrained by ${\mathbf{𝐏}}_6={\mathbf{𝐏}}_3+({\mathbf{𝐏}}_3-{\mathbf{𝐏}}_0)+6({\mathbf{𝐏}}_1-{\mathbf{𝐏}}_2+{\mathbf{𝐏}}_3-{\mathbf{𝐏}}_2)$. Applying this constraint to the cubic Bézier spline will cause a complete loss of local control, as the entire spline is now fully constrained and defined by the first curve's control points. In fact, it is arguably no longer a spline, as its shape is now equivalent to extrapolating the first curve indefinitely, making it not only $C^3$ continuous, but $C^{\mathrm{\infty }}$, as joins between separate curves no longer exist

Approximating circular arcs

In case circular arc primitives are not supported in a particular environment, they may be approximated by Bézier curves. Commonly, eight quadratic segments or four cubic segments are used to approximate a circle. It is desirable to find the length $\mathbf{𝐤}$ of control points which result in the least approximation error for a given number of cubic segments.

Using four curves

Considering only the 90-degree unit-circular arc in the first quadrant, we define the endpoints $\mathbf{𝐀}$ and $\mathbf{𝐁}$ with control points ${\mathbf{𝐀}}^{\prime }$ and ${\mathbf{𝐁}}^{\prime }$, respectively, as:

$\begin{array}{rl}\mathbf{𝐀}& =[0,1]\\ {\mathbf{𝐀}}^{\prime }& =[\mathbf{𝐤},1]\\ {\mathbf{𝐁}}^{\prime }& =[1,\mathbf{𝐤}]\\ \mathbf{𝐁}& =[1,0]\\ \end{array}$

From the definition of the cubic Bézier curve, we have:

$\mathbf{𝐂}(t)=(1-t{)}^3\mathbf{𝐀}+3(1-t{)}^{2}t{\mathbf{𝐀}}^{\prime }+3(1-t)t^2{\mathbf{𝐁}}^{\prime }+t^3\mathbf{𝐁}$

With the point $\mathbf{𝐂}(t=0.5)$ as the midpoint of the arc, we may write the following two equations:

$\begin{array}{rlr}\mathbf{𝐂}& =\frac{1}{8}\mathbf{𝐀}+\frac{3}{8}{\mathbf{𝐀}}^{\prime }+\frac{3}{8}{\mathbf{𝐁}}^{\prime }+\frac{1}{8}\mathbf{𝐁}& \text{(by substitution into cubic Bézier definition)}\\ \mathbf{𝐂}& =[\sqrt{2}/2,\sqrt{2}/2]& \text{(by the definition of a circle)}\end{array}$

Solving these equations for the x-coordinate (and identically for the y-coordinate) yields:

$\frac{0}{8}\mathbf{+}\frac{3}{8}\mathbf{𝐤}+\frac{3}{8}+\frac{1}{8}=\sqrt{2}/2$

$\mathbf{𝐤}=\frac{4}{3}(\sqrt{2}-1)\approx 0.5522847498$

Error analysis

Note however that the resulting Bézier curve is entirely outside the circle, with a maximum deviation of the radius of about 0.00027. For the rest of the discussion, the radial drift for a given Template:Math is $\sqrt{x(t{)}^2+y(t{)}^2}-1$.

Unattributed correction: Adding a small correction to intermediate points such as

$\begin{array}{rl}{\mathbf{𝐀}}^{\prime }& =[\mathbf{𝐤}+0.0009,1-0.00103]\\ {\mathbf{𝐁}}^{\prime }& =[1-0.00103,\mathbf{𝐤}+0.0009],\end{array}$

the magnitude of the radius deviation to 1 is reduced by a factor of about 3, to 0.000068 (at the expense of the derivability of the approximated circle curve at endpoints).

Ellsworth's three-parameter correction: Instead of solving for a single k, solve for $\mathbf{𝐀}=[0,a],$ ${\mathbf{𝐀}}^{\prime }=[b,c],$ ${\mathbf{𝐁}}^{\prime }=[c,b],$ $\mathbf{𝐁}=[a,0].$ Minimize the maximum radial drift, which is $f(a,b,c)=\max (|d(t)|,0\le t\le 1)$. Solving numerically gives Template:Math Template:Math and Template:Math. The maximum radial drift is 0.000055.

Ellsworth's least-square correction: Solve for the k that minimizes ${\int }_0^1{(d(t))}^{2}dt$. Numerically $k\approx 0.55197038140$. The maximal radial drift is 0.0002, for a 1.9% reduction in the total squared radial drift compared to the three-parameter case.

General case

We may approximate a unit circle of radius from an arbitrary number of cubic Bézier curves. Let the arc start at point $\mathbf{𝐀}$ and end at point $\mathbf{𝐁}$, placed at equal distances above and below the x-axis, spanning an arc of angle $\theta =2\varphi \le =\frac{\pi }{2}$:

$\begin{array}{rl}\mathbf{𝐀}& =[\cos (\varphi ),\sin (\varphi )]\\ \mathbf{𝐁}& =[\cos (\varphi ),-\sin (\varphi )]=[{\mathbf{𝐀}}_x,-{\mathbf{𝐀}}_y]\end{array}$

After solving for the curve starting and ending at the given points, the arc's midpoint occurring at $t=0.5$, and the curve's end point tangents being the same as the circles, the control points may be written as:

$\begin{array}{rlr}{\mathbf{𝐀}}^{\prime }& =[\frac{4-{\mathbf{𝐀}}_x}{3},\frac{(1-{\mathbf{𝐀}}_x)(3-{\mathbf{𝐀}}_x)}{3{\mathbf{𝐀}}_y}]& =[\frac{4-\cos (\varphi )}{3},\frac{(1-\cos (\varphi ))(3-\cos (\varphi ))}{3\sin (\varphi )}]\\ {\mathbf{𝐁}}^{\prime }& =[{{\mathbf{𝐀}}^{\prime }}_x,-{{\mathbf{𝐀}}^{\prime }}_y]& =[\frac{4-\cos (\varphi )}{3},\frac{(1-\cos (\varphi ))(3-\cos (\varphi ))}{-3\sin (\varphi )}]\end{array}$

Multiplying the values provide general solutions for arcs of any radius. Using rotation matrices allows the drawing of any angle.

The constraints of Ellsworth's solutions can be applied to the general case and numerically solved as well.

Examples

• Eight-segment quadratic polybézier (red) approximating a circle (black) with control points
  Eight-segment quadratic polybézier (red) approximating a circle (black) with control points
• Four-segment cubic polybézier (red) approximating a circle (black) with control points
  Four-segment cubic polybézier (red) approximating a circle (black) with control points

Other approaches

For Bézier curves of degree n > 1, Ahn (2019) provide an $G^2$ approximation for any number of segments. The Hausdorff distance as well as the aforementioned "radial drift" show 2n-1 extreme points and Ahn's approximation is, as of the time of publication, the one with the smallest maximum Hausdorff distance. Explicit, analytical solutions are provided for the quadratic, cubic, and quartic cases.

Applying affine transformation to the constraints allows the drawing of ellipses and elliptical arcs.

Fonts

TrueType fonts use composite Béziers composed of quadratic Bézier curves (2nd order curves) while most PostScript Type 1 and OpenType CFF use cubic (3rd order) curves. To describe a typical type design as a computer font to any given accuracy, 3rd order Béziers require less data than 2nd order Béziers; and these in turn require less data than a series of straight lines. This is true even though any one straight line segment requires less data than any one segment of a parabola; and that parabolic segment in turn requires less data than any one segment of a 3rd order curve.

See also

• B-spline
• Confluent hypergeometric function