Bézier curve

A Bézier curve ( ) is a parametric curve used in computer graphics and related fields. A set of discrete "control points" defines a smooth, continuous curve by means of a formula. Usually the curve is intended to approximate a real-world shape that otherwise has no mathematical representation or whose representation is unknown or too complicated. The Bézier curve is named after French engineer Pierre Bézier (1910–1999), who used it in the 1960s for designing curves for the bodywork of Renault cars. Other uses include the design of computer fonts and animation. Bézier curves can be combined to form a Bézier spline, or generalized to higher dimensions to form Bézier surfaces. The Bézier triangle is a special case of the latter.

In vector graphics, Bézier curves are used to model smooth curves that can be scaled indefinitely. "Paths", as they are commonly referred to in image manipulation programs, are combinations of linked Bézier curves. Paths are not bound by the limits of rasterized images and are intuitive to modify.

Bézier curves are also used in the time domain, particularly in animation, user interface design and smoothing cursor trajectory in eye gaze controlled interfaces. For example, a Bézier curve can be used to specify the velocity over time of an object such as an icon moving from A to B, rather than simply moving at a fixed number of pixels per step. When animators or interface designers talk about the "physics" or "feel" of an operation, they may be referring to the particular Bézier curve used to control the velocity over time of the move in question.

This also applies to robotics where the motion of a welding arm, for example, should be smooth to avoid unnecessary wear.

Invention

The mathematical basis for Bézier curves—the Bernstein polynomials—was established in 1912, but the polynomials were not applied to graphics until some 50 years later when mathematician Paul de Casteljau in 1959 developed de Casteljau's algorithm, a numerically stable method for evaluating the curves, and became the first to apply them to computer-aided design at French automaker Citroën. Yet, de Casteljau's method was patented in France but not published until the 1980s while the Bézier polynomials were widely publicised in the 1960s by the French engineer Pierre Bézier, who discovered them independently and used them to design automobile bodies at Renault.

Specific cases

A Bézier curve is defined by a set of '' control points'' P₀ through P_''n'', where ''n'' is called the order of the curve (''n'' = 1 for linear, 2 for quadratic, 3 for cubic, etc.). The first and last control points are always the endpoints of the curve; however, the intermediate control points (if any) generally do not lie on the curve. The sums in the following sections are to be understood as affine combinations – that is, the coefficients sum to 1.

Linear Bézier curves

Given distinct points P₀ and P₁, a linear Bézier curve is simply a line between those two points. The curve is given by

:$\mathbf(t)=\mathbf_0 + t(\mathbf_1-\mathbf_0)=(1-t)\mathbf_0 + t\mathbf_1,\ 0 \le t \le 1$
and is equivalent to linear interpolation.

Quadratic Bézier curves

A quadratic Bézier curve is the path traced by the function B(''t''), given points P₀, P₁, and P₂,

: \mathbf(t) = (1 - t) 1 - t) \mathbf P_0 + t \mathbf P_1+ t 1 - t) \mathbf P_1 + t \mathbf P_2\ 0 \le t \le 1,

which can be interpreted as the linear interpolant of corresponding points on the linear Bézier curves from P₀ to P₁ and from P₁ to P₂ respectively. Rearranging the preceding equation yields:

: $\mathbf(t) = (1 - t)^\mathbf_0 + 2(1 - t)t\mathbf_1 + t^\mathbf_2,\ 0 \le t \le 1.$

This can be written in a way that highlights the symmetry with respect to P₁:

: $\mathbf(t) = \mathbf_1+(1 - t)^(\mathbf_0 - \mathbf_1) + t^(\mathbf_2-\mathbf_1),\ 0 \le t \le 1.$

Which immediately gives the derivative of the Bézier curve with respect to ''t'':

: $\mathbf'(t) = 2 (1 - t) (\mathbf_1 - \mathbf_0) + 2 t (\mathbf_2 - \mathbf_1),$

from which it can be concluded that the tangents to the curve at P₀ and P₂ intersect at P₁. As ''t'' increases from 0 to 1, the curve departs from P₀ in the direction of P₁, then bends to arrive at P₂ from the direction of P₁.

The second derivative of the Bézier curve with respect to ''t'' is

: $\mathbf''(t) = 2(\mathbf_2 - 2 \mathbf_1 + \mathbf_0).$

Cubic Bézier curves

Four points P₀, P₁, P₂ and P₃ in the plane or in higher-dimensional space define a cubic Bézier curve.
The curve starts at P₀ going toward P₁ and arrives at P₃ coming from the direction of P₂. Usually, it will not pass through P₁ or P₂; these points are only there to provide directional information. The distance between P₁ and P₂ determines "how far" and "how fast" the curve moves towards P₁ before turning towards P₂.

Writing B_{P''i'',P''j'',P''k''}(''t'') for the quadratic Bézier curve defined by points P_''i'', P_''j'', and P_''k'', the cubic Bézier curve can be defined as an affine combination of two quadratic Bézier curves:

:$\mathbf(t) = (1-t)\mathbf_(t) + t \mathbf_(t),\ 0 \le t \le 1.$

The explicit form of the curve is:

:$\mathbf(t) = (1-t)^3\mathbf_0+3(1-t)^2t\mathbf_1+3(1-t)t^2\mathbf_2+t^3\mathbf_3,\ 0 \le t \le 1.$

For some choices of P₁ and P₂ the curve may intersect itself, or contain a cusp.

Any series of 4 distinct points can be converted to a cubic Bézier curve that goes through all 4 points in order.
Given the starting and ending point of some cubic Bézier curve, and the points along the curve corresponding to ''t'' = 1/3 and ''t'' = 2/3, the control points for the original Bézier curve can be recovered.

The derivative of the cubic Bézier curve with respect to ''t'' is
: $\mathbf'(t) = 3(1-t)^2(\mathbf_1 - \mathbf_0) + 6(1-t)t(\mathbf_2 - \mathbf_1) + 3t^2(\mathbf_3 - \mathbf_2) \,.$

The second derivative of the Bézier curve with respect to ''t'' is
: $\mathbf''(t) = 6(1-t)(\mathbf_2 - 2 \mathbf_1 + \mathbf_0) + 6t(\mathbf_3 - 2 \mathbf_2 + \mathbf_1) \,.$

General definition

Bézier curves can be defined for any degree ''n''.

Recursive definition

A recursive definition for the Bézier curve of degree ''n'' expresses it as a point-to-point linear combination (linear interpolation) of a pair of corresponding points in two Bézier curves of degree ''n'' − 1.

Let $\mathbf_$ denote the Bézier curve determined by any selection of points P₀, P₁, ..., P_''k''. Then to start,

:$\mathbf_(t) = \mathbf_0 \text$

:$\mathbf(t) = \mathbf_(t) = (1-t)\mathbf_(t) + t\mathbf_(t)$

This recursion is elucidated in the animations below.

Explicit definition

The formula can be expressed explicitly as follows (where t⁰ and (1-t)⁰ are extended continuously to be 1 throughout ,1:

:$\begin \mathbf(t) &= \sum_^n (1 - t)^t^i\mathbf_i \\ &=(1 - t)^n\mathbf_0 + (1 - t)^t\mathbf_1 + \cdots + (1 - t)t^\mathbf_ + t^n\mathbf_n, && 0 \leqslant t \leqslant 1 \end$

where $\scriptstyle$ are the binomial coefficients.

For example, when ''n'' = 5:

:$\begin \mathbf(t) &= (1 - t)^5\mathbf_0 + 5t(1 - t)^4\mathbf_1 + 10t^2(1 - t)^3 \mathbf_2 + 10t^3 (1-t)^2 \mathbf_3 + 5t^4(1-t) \mathbf_4 + t^5 \mathbf_5, && 0 \leqslant t \leqslant 1. \end$

Terminology

Some terminology is associated with these parametric curves. We have

:$\mathbf(t) = \sum_^n b_(t)\mathbf_i,\ \ \ 0 \le t \le 1$

where the polynomials

:$b_(t) = t^i (1 - t)^,\ \ \ i = 0, \ldots, n$

are known as Bernstein basis polynomials of degree ''n''.

Note that ''t''⁰ = 1, (1 − ''t'')⁰ = 1, and that the binomial coefficient, $\scriptstyle$, is:

:$= \frac.$

The points P_''i'' are called ''control points'' for the Bézier curve. The polygon formed by connecting the Bézier points with lines, starting with P₀ and finishing with P_''n'', is called the ''Bézier polygon'' (or ''control polygon''). The convex hull of the Bézier polygon contains the Bézier curve.

Polynomial form

Sometimes it is desirable to express the Bézier curve as a polynomial instead of a sum of less straightforward Bernstein polynomials. Application of the binomial theorem to the definition of the curve followed by some rearrangement will yield

:$\mathbf(t) = \sum_^n$
where
:$\mathbf_j = \frac \sum_^j \frac = \prod_^ (n - m) \sum_^j \frac.$

This could be practical if $\mathbf_j$ can be computed prior to many evaluations of $\mathbf(t)$; however one should use caution as high order curves may lack numeric stability (de Casteljau's algorithm should be used if this occurs). Note that the empty product is 1.

Properties

* The curve begins at P₀ and ends at P_''n''; this is the so-called ''endpoint interpolation'' property.
* The curve is a line if and only if all the control points are collinear.
* The start and end of the curve is tangent to the first and last section of the Bézier polygon, respectively.
* A curve can be split at any point into two subcurves, or into arbitrarily many subcurves, each of which is also a Bézier curve.
* Some curves that seem simple, such as the circle, cannot be described exactly by a Bézier or piecewise Bézier curve; though a four-piece cubic Bézier curve can approximate a circle (see composite Bézier curve), with a maximum radial error of less than one part in a thousand, when each inner control point (or offline point) is the distance $\textstyle\frac$ horizontally or vertically from an outer control point on a unit circle. More generally, an ''n''-piece cubic Bézier curve can approximate a circle, when each inner control point is the distance $\textstyle\frac\tan(t/4)$ from an outer control point on a unit circle, where ''t'' is 360/''n'' degrees, and ''n'' > 2.
* Every quadratic Bézier curve is also a cubic Bézier curve, and more generally, every degree ''n'' Bézier curve is also a degree ''m'' curve for any ''m'' > ''n''. In detail, a degree ''n'' curve with control points P₀, ..., P_''n'' is equivalent (including the parametrization) to the degree ''n'' + 1 curve with control points P'₀, ..., P'_{''n'' + 1}, where $\mathbf P'_k=\tfrac\mathbf P_+\left(1-\tfrac\right)\mathbf P_k$.
* Bézier curves have the variation diminishing property. What this means in intuitive terms is that a Bézier curve does not "undulate" more than the polygon of its control points, and may actually "undulate" less than that.
* There is no local control in degree ''n'' Bézier curves—meaning that any change to a control point requires recalculation of and thus affects the aspect of the entire curve, "although the further that one is from the control point that was changed, the smaller is the change in the curve."
* A Bézier curve of order higher than two may intersect itself or have a cusp for certain choices of the control points.

Second-order curve is a parabolic segment

A quadratic Bézier curve is also a segment of a parabola. As a parabola is a conic section, some sources refer to quadratic Béziers as "conic arcs". With reference to the figure on the right, the important features of the parabola can be derived as follows:

# Tangents to the parabola at the endpoints of the curve (A and B) intersect at its control point (C).
# If D is the midpoint of AB, the tangent to the curve which is perpendicular to CD (dashed cyan line) defines its vertex (V). Its axis of symmetry (dash-dot cyan) passes through V and is perpendicular to the tangent.
# E is either point on the curve with a tangent at 45° to CD (dashed green). If G is the intersection of this tangent and the axis, the line passing through G and perpendicular to CD is the directrix (solid green).
# The focus (F) is at the intersection of the axis and a line passing through E and perpendicular to CD (dotted yellow). The latus rectum is the line segment within the curve (solid yellow).

Derivative

The derivative for a curve of order ''n'' is

:$\mathbf'(t) = n \sum_^ b_(t) (\mathbf_ - \mathbf_i).$

Constructing Bézier curves

Linear curves

The ''t'' in the function for a linear Bézier curve can be thought of as describing how far B(''t'') is from P₀ to P₁. For example, when ''t''=0.25, B(''t'') is one quarter of the way from point P₀ to P₁. As ''t'' varies from 0 to 1, B(''t'') describes a line from P₀ to P₁.

Quadratic curves

For quadratic Bézier curves one can construct intermediate points Q₀ and Q₁ such that as ''t'' varies from 0 to 1:
* Point Q₀(''t'') varies from P₀ to P₁ and describes a linear Bézier curve.
* Point Q₁(''t'') varies from P₁ to P₂ and describes a linear Bézier curve.
* Point B(''t'') is interpolated linearly between Q₀(''t'') to Q₁(''t'') and describes a quadratic Bézier curve.

Higher-order curves

For higher-order curves one needs correspondingly more intermediate points. For cubic curves one can construct intermediate points Q₀, Q₁, and Q₂ that describe linear Bézier curves, and points R₀ and R₁ that describe quadratic Bézier curves:

For fourth-order curves one can construct intermediate points Q₀, Q₁, Q₂ and Q₃ that describe linear Bézier curves, points R₀, R₁ and R₂ that describe quadratic Bézier curves, and points S₀ and S₁ that describe cubic Bézier curves:

For fifth-order curves, one can construct similar intermediate points.

These representations rest on the process used in De Casteljau's algorithm to calculate Bézier curves.

Offsets (or stroking) of Bézier curves

The curve at a fixed offset from a given Bézier curve, called an offset or parallel curve in mathematics (lying "parallel" to the original curve, like the offset between rails in a railroad track), cannot be exactly formed by a Bézier curve (except in some trivial cases). In general, the two-sided offset curve of a cubic Bézier is a 10th-order algebraic curve and more generally for a Bézier of degree ''n'' the two-sided offset curve is an algebraic curve of degree 4''n'' − 2. However, there are heuristic methods that usually give an adequate approximation for practical purposes.

In the field of vector graphics, painting two symmetrically distanced offset curves is called ''stroking'' (the Bézier curve or in general a path of several Bézier segments). The conversion from offset curves to filled Bézier contours is of practical importance in converting fonts defined in Metafont, which require stroking of Bézier curves, to the more widely used PostScript type 1 fonts, which only require (for efficiency purposes) the mathematically simpler operation of filling a contour defined by (non-self-intersecting) Bézier curves.

Degree elevation

A Bézier curve of degree ''n'' can be converted into a Bézier curve of degree ''n'' + 1 ''with the same shape''. This is useful if software supports Bézier curves only of specific degree. For example, systems that can only work with cubic Bézier curves can implicitly work with quadratic curves by using their equivalent cubic representation.

To do degree elevation, we use the equality $\mathbf(t) = (1-t)\mathbf(t) + t\mathbf(t).$ Each component $\mathbf_(t)\mathbf_i$ is multiplied by (1 − ''t'') and ''t'', thus increasing a degree by one, without changing the value. Here is the example of increasing degree from 2 to 3.

:$\begin (1 - t)^2 \mathbf_0 + 2(1 - t)t\mathbf_1 + t^2 \mathbf_2 &= (1 - t)^3 \mathbf_0 + 2(1 - t)^2 t\mathbf_1 + (1 - t)t^2 \mathbf_2 + (1 - t)^t\mathbf_0 + 2(1 - t)t^2 \mathbf_1 + t^3 \mathbf_2 \\ &= (1 - t)^3 \mathbf_0 + (1 - t)^2 t \left( \mathbf_0 + 2\mathbf_1\right) + (1 - t) t^2 \left(2\mathbf_1 + \mathbf_2\right) + t^\mathbf_2 \\ &= (1 - t)^3 \mathbf_0 + 3(1 - t)^2 t \tfrac \left( \mathbf_0 + 2\mathbf_1 \right) + 3(1 - t) t^2 \tfrac \left( 2\mathbf_1 + \mathbf_2 \right) + t^\mathbf_2 \end$

For arbitrary ''n'' we use equalities

:$\begin (1 - t)\mathbf_ = \mathbf_ \\ t\mathbf_ = \mathbf_ \end \quad \implies \quad \begin (1 - t)\mathbf_ = \frac \mathbf_ \\ t\mathbf_ = \frac \mathbf_ \end$

Therefore:

:$\begin \mathbf(t) &= (1 - t)\sum_^n \mathbf_(t)\mathbf_i + t\sum_^n \mathbf_(t)\mathbf_i \\ &= \sum_^n \frac\mathbf_(t)\mathbf_i + \sum_^n \frac\mathbf_(t)\mathbf_i \\ &= \sum_^ \left(\frac\mathbf_ + \frac\mathbf_i\right) \mathbf_(t) \\ &= \sum_^ \mathbf_(t)\mathbf_i \end$

introducing arbitrary $\mathbf_$ and $\mathbf_$.

Therefore, new control points are

: $\mathbf_i = \frac\mathbf_ + \frac\mathbf_i,\quad i=0, \ldots, n + 1.$

Repeated degree elevation

The concept of degree elevation can be repeated on a control polygon R to get a sequence of control polygons R, R₁, R₂, and so on. After ''r'' degree elevations, the polygon R_''r'' has the vertices P_0,''r'', P_1,''r'', P_2,''r'', ..., P_{''n''+''r'',''r''} given by

:$\mathbf_ = \sum_^n \mathbf_j \tbinom nj \frac.$

It can also be shown that for the underlying Bézier curve ''B'',

:$\mathbf = \mathbf.$

Rational Bézier curves

The rational Bézier curve adds adjustable weights to provide closer approximations to arbitrary shapes. The numerator is a weighted Bernstein-form Bézier curve and the denominator is a weighted sum of Bernstein polynomials. Rational Bézier curves can, among other uses, be used to represent segments of conic sections exactly, including circular arcs.

Given ''n'' + 1 control points P₀, ..., P_''n'', the rational Bézier curve can be described by
:$\mathbf(t) = \frac ,$
or simply
:$\mathbf(t) = \frac .$
The expression can be extended by using number systems besides reals for the weights. In the complex plane the points , , and with weights , , and generate a full circle with radius one. For curves with points and weights on a circle, the weights can be scaled without changing the curve's shape. Scaling the central weight of the above curve by 1.35508 gives a more uniform parameterization.

Applications

Computer graphics

Bézier curves are widely used in computer graphics to model smooth curves. As the curve is completely contained in the convex hull of its control points, the points can be graphically displayed and used to manipulate the curve intuitively. Affine transformations such as translation and rotation can be applied on the curve by applying the respective transform on the control points of the curve.

Quadratic and cubic Bézier curves are most common. Higher degree curves are more computationally expensive to evaluate. When more complex shapes are needed, low order Bézier curves are patched together, producing a composite Bézier curve. A composite Bézier curve is commonly referred to as a "path" in vector graphics languages (like PostScript), vector graphics standards (like SVG) and vector graphics programs (like Artline, Timeworks Publisher, Adobe Illustrator, CorelDraw, Inkscape, and Allegro). In order to join Bézier curves into a composite Bézier curve without kinks, a property called ''G1 continuous'', it suffices to force the control point at which two constituent Bézier curves meet to lie on the line defined by the two control points on either side.

The simplest method for scan converting ( rasterizing) a Bézier curve is to evaluate it at many closely spaced points and scan convert the approximating sequence of line segments. However, this does not guarantee that the rasterized output looks sufficiently smooth, because the points may be spaced too far apart. Conversely it may generate too many points in areas where the curve is close to linear. A common adaptive method is recursive subdivision, in which a curve's control points are checked to see if the curve approximates a line to within a small tolerance. If not, the curve is subdivided parametrically into two segments, 0 ≤ ''t'' ≤ 0.5 and 0.5 ≤ ''t'' ≤ 1, and the same procedure is applied recursively to each half. There are also forward differencing methods, but great care must be taken to analyse error propagation.

Analytical methods where a Bézier is intersected with each scan line involve finding roots of cubic polynomials (for cubic Béziers) and dealing with multiple roots, so they are not often used in practice.

The rasterisation algorithm used in Metafont is based on discretising the curve, so that it is approximated by a sequence of "rook moves" that are purely vertical or purely horizontal, along the pixel boundaries. To that end, the plane is first split into eight 45° sectors (by the coordinate axes and the two lines $y = \pm x$), then the curve is decomposed into smaller segments such that the ''direction'' of a curve segment stays within one sector; since the curve velocity is a second degree polynomial, finding the $t$ values where it is parallel to one of these lines can be done by solving quadratic equations. Within each segment, either horizontal or vertical movement dominates, and the total number of steps in either direction can be read off from the endpoint coordinates; in for example the 0–45° sector horizontal movement to the right dominates, so it only remains to decide between which steps to the right the curve should make a step up.

There is also a modified curve form of Bresenham's line drawing algorithm by Zingl that performs this rasterization by subdividing the curve into rational pieces and calculating the error at each pixel location such that it either travels at a 45° angle or straight depending on compounding error as it iterates through the curve. This reduces the next step calculation to a series of integer additions and subtractions.Zingl, Alois "A Rasterizing Algorithm for Drawing Curves" (2012) http://members.chello.at/~easyfilter/Bresenham.pdf

Animation

In animation applications, such as Adobe Flash and Synfig, Bézier curves are used to outline, for example, movement. Users outline the wanted path in Bézier curves, and the application creates the needed frames for the object to move along the path.

In 3D animation, Bézier curves are often used to define 3D paths as well as 2D curves for keyframe interpolation. Bézier curves are now very frequently used to control the animation easing in CSS, JavaScript, JavaFx and Flutter SDK.

Fonts

TrueType fonts use composite Bézier curves composed of quadratic Bézier curves. Other languages and imaging tools (such as PostScript, Asymptote, Metafont, and SVG) use composite Béziers composed of cubic Bézier curves for drawing curved shapes. OpenType fonts can use either kind of curve, depending on which font technology underlies the OpenType wrapper.

The internal rendering of all Bézier curves in font or vector graphics renderers will split them recursively up to the point where the curve is flat enough to be drawn as a series of linear or circular segments. The exact splitting algorithm is implementation dependent, only the flatness criteria must be respected to reach the necessary precision and to avoid non-monotonic local changes of curvature. The "smooth curve" feature of charts in Microsoft Excel also uses this algorithm.

Because arcs of circles and ellipses cannot be exactly represented by Bézier curves, they are first approximated by Bézier curves, which are in turn approximated by arcs of circles. This is inefficient as there exists also approximations of all Bézier curves using arcs of circles or ellipses, which can be rendered incrementally with arbitrary precision. Another approach, used by modern hardware graphics adapters with accelerated geometry, can convert exactly all Bézier and conic curves (or surfaces) into NURBS, that can be rendered incrementally without first splitting the curve recursively to reach the necessary flatness condition. This approach also preserves the curve definition under all linear or perspective 2D and 3D transforms and projections.

Font engines, like FreeType, draw the font's curves (and lines) on a pixellated surface using a process known as font rasterization.

Robotics

Bézier curves can be used in robotics to produce trajectories of an end-effector due to the virtue of the control polygon’s ability to give a clear indication of whether the path is colliding with any nearby obstacle or object. Furthermore, joint space trajectories can be accurately differentiated using Bézier curves. Consequently, the derivatives of joint space trajectories are used in the calculation of the dynamics and control effort (torque profiles) of the robotic manipulator.

See also

* Bézier surface
* B-spline
* GEM/4 & GEM/5
* Hermite curve
* NURBS
* String art – Bézier curves are also formed by many common forms of string art, where strings are looped across a frame of nails.
* Variation diminishing property of Bézier curves

Notes

References

Citations

Sources

*
*
* Excellent discussion of implementation details; available for free as part of the TeX distribution.
*
*
*

Further reading

A Primer on Bézier Curves
an open source online book explaining Bézier curves and associated graphics algorithms, with interactive graphics

Cubic Bezier Curves – Under the Hood (video)
video showing how computers render a cubic Bézier curve, by Peter Nowell


Feature Column from American Mathematical Society
*
*
* This book is out of print and freely available from the author.
*
*
* (60 pages)
*
*
*
* Hovey, Chad (2022). Formulation and Python Implementation of Bézier and B-Spline Geometry
SAND2022-7702C
(153 pages)

External links

; Computer code
TinySpline: Open source C-library for NURBS, B-splines and Bézier curves with bindings for various languages

C++ library to generate Bézier functions at compile timeSimple Bézier curve implementation via recursive method in Python

{{DEFAULTSORT:Bezier curve
Graphic design
Interpolation
Curves
Design