Bézierkurve

Bézier-Kurve 3. Grades und ihr Kontrollpolygon
Bézierkurven der Grade 1, 2 und 3 (rot) und die zugehörigen Kontrollpolygone (grau). Von links nach rechts wurde jeweils ein weiterer Kontrollpunkt (blau) hinzu­gefügt. Man erkennt, wie die Kurve bei Einfügen/Verändern eines Kontrollpunkts ihre Richtung und/oder Krümmung variiert.
Bézierkurven 1. 2. 3. Ordnung in Geogebra – siehe auch interaktives Geogebra-Applet.

Die Bézierkurve [be'zje…] ist eine parametrisch modellierte Kurve, die ein wichtiges Werkzeug bei der Beschreibung von Freiformkurven und -flächen darstellt.

In der Computergrafik finden Bézierkurven wegen ihrer optischen Eleganz und der verhältnismäßig leichten mathematischen Handhabbarkeit häufig Anwendung. Sie werden zur Definition von Kurven und Flächen in Vektorgrafiken genutzt. Mögliche Anwendungsfälle finden sich z. B. im Computer Aided Design, bei der Erstellung von Illustrationen (siehe z. B. SVG) oder der Beschreibung von Schrifttypen (z. B. Postscript, Type1, TrueType und CFF-OpenType).

Die Bézierkurve wurde Anfang der 1960er Jahre unabhängig voneinander von Pierre Bézier bei Renault und Paul de Casteljau bei Citroën für Computer-Aided Design (computerunterstützte Konstruktion) entwickelt. Paul de Casteljau gelang zwar die Entdeckung früher, Citroën hielt seine Forschungen jedoch bis zum Ende der 1960er Jahre als Betriebsgeheimnis zurück.

Verallgemeinerungen des Konzepts der Bézierkurven führen zu den Bézierflächen.

Motivation, Definition

Numerisch einfache Kurven in der Ebene sind solche, die mit Hilfe einer Parameterdarstellung x ( t ) = ( x ( t ) , y ( t ) ) T , t 1 t t 2 , {\displaystyle {\bf {x}}(t)=(x(t),y(t))^{T},\;t_{1}\leq t\leq t_{2},} beschrieben werden, wobei x ( t ) {\displaystyle x(t)} und y ( t ) {\displaystyle y(t)} Polynome in t {\displaystyle t} sind. Ist

x ( t ) = a 0 + a 1 t + a 2 t 2 + + a n t n {\displaystyle x(t)=a_{0}+a_{1}t+a_{2}t^{2}+\cdots +a_{n}t^{n}} ,
y ( t ) = b 0 + b 1 t + b 2 t 2 + + b n t n {\displaystyle y(t)=b_{0}+b_{1}t+b_{2}t^{2}+\cdots +b_{n}t^{n}}

und setzt man a i = ( a i , b i ) T {\displaystyle {\bf {a}}_{i}=(a_{i},b_{i})^{T}} , so lässt sich die Kurve übersichtlicher durch

x ( t ) = a 0 + a 1 t + a 2 t 2 + + a n t n   {\displaystyle {\bf {x}}(t)={\bf {a}}_{0}+{\bf {a}}_{1}t+{\bf {a}}_{2}t^{2}+\cdots +{\bf {a}}_{n}t^{n}\ }

beschreiben.

Im Allgemeinen lässt sich direkt aus den Koeffizienten-Punkten a i {\displaystyle {\bf {a}}_{i}} nur wenig über den Kurvenverlauf aussagen. Lediglich a 0 {\displaystyle {\bf {a}}_{0}} (Anfangspunkt der Kurve) und a 1 {\displaystyle {\bf {a}}_{1}} (Tangentialvektor der Kurve an a 0 {\displaystyle {\bf {a}}_{0}} ) haben konkrete geometrische Bedeutungen. Dies ändert sich, wenn man die Polynome x ( t ) , y ( t ) {\displaystyle x(t),y(t)} nicht in der Monom-Basis { 1 , t , t 2 , , t n } {\displaystyle \{1,t,t^{2},\dotsc ,t^{n}\}} , sondern in der folgenden Bernsteinbasis { B 0 n ( t ) , B 1 n ( t ) , , B n n ( t ) } {\displaystyle \{B_{0}^{n}(t),B_{1}^{n}(t),\dotsc ,B_{n}^{n}(t)\}} darstellt: B i n ( t ) := ( n i ) t i ( 1 t ) n i , 0 i n . {\displaystyle \quad B_{i}^{n}(t):={n \choose i}t^{i}(1-t)^{n-i},\quad 0\leq i\leq n.}

Es sei nun n > 0 {\displaystyle n>0} festgewählt und die Vektoren b 0 , b 1 , , b n {\displaystyle {\bf {b}}_{0},{\bf {b}}_{1},\dotsc ,{\bf {b}}_{n}} beschreiben ein ebenes oder räumliches Polygon. Dann heißt die Darstellung x = b ( t )   =   i = 0 n ( n i ) t i ( 1 t ) n i b i   =   b 0 B 0 n ( t ) + b 1 B 1 n ( t ) + + b n B n n ( t ) , 0 t 1 , {\displaystyle \quad {\begin{aligned}{\bf {x}}={\bf {b}}(t)\ &=\ \sum _{i=0}^{n}{\binom {n}{i}}t^{i}(1-t)^{n-i}{\bf {b}}_{i}\\\ &=\ {\bf {b}}_{0}B_{0}^{n}(t)+{\bf {b}}_{1}B_{1}^{n}(t)+\cdots +{\bf {b}}_{n}B_{n}^{n}(t),\quad 0\leq t\leq 1,\end{aligned}}} eine Bézierkurve[1][2] vom (maximalen) Grad n {\displaystyle n} . Die Punkte b 0 , , b n {\displaystyle {\bf {b}}_{0},\dotsc ,{\bf {b}}_{n}} nennt man Kontrollpunkte der Bézierkurve.

Eigenschaften der Bernsteinpolynome:

Bernsteinpolynome
  1. B 0 n ( t ) + B 1 n ( t ) + + B n n ( t ) = 1 , {\displaystyle B_{0}^{n}(t)+B_{1}^{n}(t)+\cdots +B_{n}^{n}(t)=1,}
  2. B 0 n ( 0 ) = 1 , B i n ( 0 ) = 0 {\displaystyle B_{0}^{n}(0)=1,\;B_{i}^{n}(0)=0} für i > 0 , B n n ( 1 ) = 1 , B i n ( 1 ) = 0 {\displaystyle i>0,\quad B_{n}^{n}(1)=1,\;B_{i}^{n}(1)=0} für i < n . {\displaystyle i<n.}
  3. Das Bernsteinpolynom B i n {\displaystyle B_{i}^{n}} hat genau ein Maximum und zwar an der Stelle t = i / n . {\displaystyle t=i/n.} D. h. eine leichte Veränderung des Punktes b i {\displaystyle {\bf {b}}_{i}} hat nur in der Umgebung von x ( i / n ) {\displaystyle {\bf {x}}(i/n)} eine wesentliche Veränderung der Kurve zur Folge.

Eigenschaften einer Bézierkurve:

  1. b 0 {\displaystyle {\bf {b}}_{0}} ist der Anfangs-, b n {\displaystyle {\bf {b}}_{n}} der Endpunkt
  2. b 1 b 0 {\displaystyle {\bf {b}}_{1}-{\bf {b}}_{0}} ist die Richtung der Tangente im Punkt b 0 = x ( 0 ) , b n 1 b n {\displaystyle {\bf {b}}_{0}={\bf {x}}(0),\;{\bf {b}}_{n-1}-{\bf {b}}_{n}} ist die Richtung der Tangente im Punkt b n = x ( 1 ) . {\displaystyle {\bf {b}}_{n}={\bf {x}}(1).}
  3. Das Polygon b 0 , b 1 , , b n {\displaystyle {\bf {b}}_{0},{\bf {b}}_{1},\dotsc ,{\bf {b}}_{n}} gibt einen ungefähren Verlauf der Kurve an.

Weitere Eigenschaften der Bernsteinbasis

Für Untersuchungen von Bézierkurven sind die folgenden Eigenschaften[3] nützlich:

Beziehung zwischen der Bernstein- und der Monom-Basis
(MB) t i = j = i n ( j i ) ( n i ) B j n ( t ) , B i n ( t ) = j = i n ( 1 ) j i ( n j ) ( j i ) t j , {\displaystyle \qquad t^{i}=\sum _{j=i}^{n}{\frac {j \choose i}{n \choose i}}B_{j}^{n}(t),\qquad B_{i}^{n}(t)=\sum _{j=i}^{n}(-1)^{j-i}{n \choose j}{j \choose i}t^{j},}
Rekursion
(R) B i n ( t ) = ( 1 t ) B i n 1 ( t ) + t B i 1 n 1 ( t ) , {\displaystyle \qquad B_{i}^{n}(t)=(1-t)B_{i}^{n-1}(t)+tB_{i-1}^{n-1}(t),}
Skalierung
(S) B i n ( c t ) = j = 0 n B i j ( c ) B j n ( t ) , {\displaystyle \qquad B_{i}^{n}(ct)=\sum _{j=0}^{n}B_{i}^{j}(c)B_{j}^{n}(t),}
Ableitung
(A) d d t B i n ( t ) = n ( B i 1 n 1 ( t ) B i n 1 ( t ) ) , {\displaystyle \qquad {\frac {d}{dt}}B_{i}^{n}(t)=n\left(B_{i-1}^{n-1}(t)-B_{i}^{n-1}(t)\right),}
(Man beachte, dass B 1 ( t ) = 0 , B n n 1 ( t ) = 0 {\displaystyle B_{-1}^{\cdots }(t)=0,B_{n}^{n-1}(t)=0} ist.)
Produkt
(P) B i m ( t ) B j n ( t ) = ( m i ) ( n j ) ( m + n i + j ) B i + j m + n ( t ) . {\displaystyle \qquad B_{i}^{m}(t)B_{j}^{n}(t)={\frac {{m \choose i}{n \choose j}}{m+n \choose i+j}}B_{i+j}^{m+n}(t).}

Weitere Eigenschaften einer Bézierkurve

In der Literatur[4] werden noch weitere Eigenschaften einer Bézierkurve aufgelistet:

  • Die Kurve liegt innerhalb der konvexen Hülle des Kontrollpolygons. Dies folgt daraus, dass die Bernsteinpolynome vom Grad n {\displaystyle n} eine Zerlegung der Eins sind:
1 = i = 0 n B i n ( t ) t [ 0 , 1 ] {\displaystyle 1=\sum _{i=0}^{n}B_{i}^{n}(t)\quad t\in [0,1]}
  • Die ersten Summanden des Taylorpolynoms bei t = 0 {\displaystyle t=0} bzw. bei t = 1 {\displaystyle t=1} lauten für n 2 {\displaystyle n\geq 2} :
b ( t ) = b 0 + n ( b 1 b 0 ) t + n ( n 1 ) 2 ( b 0 2 b 1 + b 2 ) t 2 + O ( t 3 ) {\displaystyle {\bf {b}}(t)={\bf {b}}_{0}+n({\bf {b}}_{1}-{\bf {b}}_{0})t+{\frac {n(n-1)}{2}}({\bf {b}}_{0}-2{\bf {b}}_{1}+{\bf {b}}_{2})t^{2}+{\mathcal {O}}(t^{3})}
b ( t ) = b n + n ( b n 1 b n ) ( 1 t ) + n ( n 1 ) 2 ( b n 2 2 b n 1 + b n ) ( 1 t ) 2 + O ( ( 1 t ) 3 ) {\displaystyle {\bf {b}}(t)={\bf {b}}_{n}+n({\bf {b}}_{n-1}-{\bf {b}}_{n})(1-t)+{\frac {n(n-1)}{2}}({\bf {b}}_{n-2}-2{\bf {b}}_{n-1}+{\bf {b}}_{n})(1-t)^{2}+{\mathcal {O}}((1-t)^{3})}
  • Eine Gerade schneidet eine Bézierkurve höchstens so oft, wie sie ihr Kontrollpolygon schneidet (die Kurve ist variationsreduzierend, bzw. hat eine beschränkte Schwankung).
  • Eine affine Transformation (Verschiebung, Skalierung, Rotation, Scherung) kann auf die Bézierkurve durch Transformation des Kontrollpolygons angewendet werden („affine Invarianz“).
  • Liegen alle Kontrollpunkte auf einer Geraden, so wird die Bézierkurve zu einer Strecke (Vorteil gegenüber der Polynominterpolation).
  • Der Einfluss eines Kontrollpunktes auf die Kurve ist global. Das heißt: Verschiebt man einen Punkt, verändert sich die gesamte Kurve. Daher verwendet man in der Praxis meist Splines, zusammengesetzte Kurven festen Grades, die stetig ineinander übergehen.
  • Eine Bézierkurve kann immer in zwei Bézierkurven gleicher Ordnung geteilt werden, wobei sich die neuen Kontrollpunkte aus den alten mit Hilfe des De-Casteljau-Algorithmus ergeben (s. Abschnitt Teilung einer Bézierkurve).

Der De-Casteljau-Algorithmus

Hauptartikel: De-Casteljau-Algorithmus

Höhere Potenzen von t {\displaystyle t} auszurechnen ist numerisch instabil. Der folgende Algorithmus führt deshalb die Berechnung eines Kurvenpunktes auf wiederholte lineare Interpolation zurück. In jedem Schritt wird mittels linearer Interpolation ein neues um 1 kürzeres Polygon berechnet (s. Bild). Bei der letzten Interpolation entsteht schließlich der Kurvenpunkt:

De-Casteljau-Algorithmus für eine Bézier-Kurve 3. Grades

Für das Polygon b 0 , b 1 , , b n {\displaystyle \;{\bf {b}}_{0},{\bf {b}}_{1},\dots ,{\bf {b}}_{n}\;} im R 2 {\displaystyle \mathbb {R} ^{2}} (oder R 3 {\displaystyle \mathbb {R} ^{3}} ) und einem t R {\displaystyle t\in \mathbb {R} } definiert man rekursiv für jedes r = 1 , , n {\displaystyle r=1,\dots ,n} das Polygon

  • b i r ( t ) = ( 1 t ) b i r 1 ( t ) + t b i + 1 r 1 ( t ) , i = 0 , , n r   {\displaystyle \quad {\bf {b}}_{i}^{r}(t)=(1-t){\bf {b}}_{i}^{r-1}(t)+t{\bf {b}}_{i+1}^{r-1}(t),\quad i=0,\dots ,n-r\ } erzeugt. Dabei sei   b i 0 ( t ) = b i   {\displaystyle \ {\bf {b}}_{i}^{0}(t)={\bf {b}}_{i}\ } .

Das Polygon der Stufe r = 0 {\displaystyle r=0} ist identisch mit dem Ausgangspolygon, das Polygon der Stufe r = n {\displaystyle r=n} ist ein Punkt, der Kurvenpunkt.

Aus der Rekursionseigenschaft (R) der Bernsteinpolynome folgt

b i r ( t ) = j = 0 r b i + j B j r ( t ) ,   {\displaystyle {\bf {b}}_{i}^{r}(t)=\sum _{j=0}^{r}{\bf {b}}_{i+j}B_{j}^{r}(t),\ }
für r = 0 , , n ,   i = 0 , , n r   . {\displaystyle r=0,\dots ,n,\ i=0,\dots ,n-r\ .}

(Beweis mit Hilfe vollständiger Induktion über r.) Also ist

  • b 0 n ( t ) = j = 0 n b j B j n ( t ) {\displaystyle \quad {\bf {b}}_{0}^{n}(t)=\sum _{j=0}^{n}{\bf {b}}_{j}B_{j}^{n}(t)}

die Bézierkurve mit dem Kontrollpolygon b 0 , b 1 , , b n {\displaystyle {\bf {b}}_{0},{\bf {b}}_{1},\dots ,{\bf {b}}_{n}} . Diese Methode, einen Punkt der Bézierkurve durch lineare Interpolationen zu bestimmen, heißt De-Casteljau-Algorithmus.

Wie für ein t R {\displaystyle t\in \mathbb {R} } mit Hilfe des Casteljau-Algorithmus aus dem Kontrollpolygon die Zwischenpolygone und schließlich der Punkt der Bézierkurve entsteht, zeigt die Abbildung für n = 3 {\displaystyle n=3} . Die neuen Punkte teilen immer die alten Strecken, auf denen sie liegen, im gleichen Verhältnis ( 1 t ) : t {\displaystyle (1-t):t} .

Bézierkurven bis zum dritten Grad

Lineare Bézierkurve

Lineare Bézierkurven ( n = 1 ) {\displaystyle (n=1)} :

Zwei Kontrollpunkte P 0 : b 0 {\displaystyle P_{0}:{\bf {b}}_{0}} und P 1 : b 1 {\displaystyle P_{1}:{\bf {b}}_{1}} bestimmen eine gerade Strecke zwischen diesen beiden Punkten. Der Verlauf dieser linearen Bézier„kurve“ wird beschrieben durch

b ( t ) =   ( 1 t ) b 0 + t b 1  ,    t [ 0 , 1 ] {\displaystyle {\bf {b}}(t)=\ (1-t){\bf {b_{0}}}+t{\bf {b}}_{1}{\mbox{ , }}\ t\in [0,1]}

Quadratische Bézierkurven ( n = 2 ) {\displaystyle (n=2)} :

Quadratische Bézierkurve mit Casteljau-Algorithmus

Eine quadratische Bézierkurve ist der Pfad, der durch die Funktion b ( t ) {\displaystyle {\bf {b}}(t)} für die Punkte P 0 : b 0 {\displaystyle P_{0}:{\bf {b}}_{0}} , P 1 : b 1 {\displaystyle P_{1}:{\bf {b}}_{1}} und P 2 : b 2 {\displaystyle P_{2}:{\bf {b}}_{2}} beschrieben wird:

b ( t )   =   i = 0 2 ( 2 i ) t i ( 1 t ) 2 i b i   =   ( 1 t ) 2 b 0 + 2 t ( 1 t ) b 1 + t 2 b 2   =   ( b 0 2 b 1 + b 2 ) t 2 + ( 2 b 0 + 2 b 1 ) t + b 0  ,    t [ 0 , 1 ]   . {\displaystyle {\begin{aligned}{\bf {b}}(t)\ &=\ \sum _{i=0}^{2}{\binom {2}{i}}t^{i}(1-t)^{2-i}{\bf {b}}_{i}\\\ &=\ (1-t)^{2}{\bf {b}}_{0}+2t(1-t){\bf {b}}_{1}+t^{2}{\bf {b}}_{2}\\\ &=\ ({\bf {b}}_{0}-2{\bf {b}}_{1}+{\bf {b}}_{2})t^{2}+(-2{\bf {b}}_{0}+2{\bf {b}}_{1})t+{\bf {b}}_{0}{\mbox{ , }}\ t\in [0,1]\ .\end{aligned}}}

Die letzte Zeile zeigt: Eine quadratische Bézierkurve ist eine Parabel.

Mit Hilfe des De-Casteljau-Algorithmus ausgedrückt: b ( t ) = b 0 2 ( t ) {\displaystyle {\bf {b}}(t)={\bf {b}}_{0}^{2}(t)}

b 0 2 ( t ) = ( 1 t ) [ ( 1 t ) b 0 b 0 0 + t b 1 b 1 0 ] b 0 1 ( t ) + t [ ( 1 t ) b 1 b 1 0 + t b 2 b 2 0 ] b 1 1 ( t ) {\displaystyle {\bf {b}}_{0}^{2}(t)=(1-t)\underbrace {{\Bigl [}(1-t)\overbrace {{\bf {b}}_{0}} ^{{\bf {b}}_{0}^{0}}+t\overbrace {{\bf {b}}_{1}} ^{{\bf {b}}_{1}^{0}}{\Bigr ]}} _{{\bf {b}}_{0}^{1}(t)}+t\underbrace {{\Bigl [}(1-t)\overbrace {{\bf {b}}_{1}} ^{{\bf {b}}_{1}^{0}}+t\overbrace {{\bf {b}}_{2}} ^{{\bf {b}}_{2}^{0}}{\Bigr ]}} _{{\bf {b}}_{1}^{1}(t)}}


Kubische Bézierkurven ( n = 3 ) {\displaystyle (n=3)} :

Kubische Bézierkurve mit Casteljau-Algorithmus

Kubische Bézierkurven sind in der Praxis von großer Bedeutung, da sowohl B-Spline-Kurven als auch NURBS stückweise in kubische Bézierkurven umgewandelt werden, um dann effizient mit dem De-Casteljau-Algorithmus gezeichnet zu werden. Selbiges gilt für hermitesche Splines, die in ihrer kubischen Form vor allem in der Computeranimation zur Interpolation zwischen Keyframes verwendet werden.

b ( t )   =   i = 0 3 ( 3 i ) t i ( 1 t ) 3 i b i =   ( 1 t ) 3 b 0 + 3 t ( 1 t ) 2 b 1 + 3 t 2 ( 1 t ) b 2 + t 3 b 3 =   ( b 0 + 3 b 1 3 b 2 + b 3 ) t 3 + ( 3 b 0 6 b 1 + 3 b 2 ) t 2 + ( 3 b 0 + 3 b 1 ) t + b 0   ,   t [ 0 , 1 ]   . {\displaystyle {\begin{aligned}{\bf {b}}(t)\ &=\ \sum _{i=0}^{3}{\binom {3}{i}}t^{i}(1-t)^{3-i}{\bf {b}}_{i}\\&=\ (1-t)^{3}{\bf {b}}_{0}+3t(1-t)^{2}{\bf {b}}_{1}+3t^{2}(1-t){\bf {b}}_{2}+t^{3}{\bf {b}}_{3}\\&=\ (-{\bf {b}}_{0}+3{\bf {b}}_{1}-3{\bf {b}}_{2}+{\bf {b}}_{3})t^{3}+(3{\bf {b}}_{0}-6{\bf {b}}_{1}+3{\bf {b}}_{2})t^{2}+(-3{\bf {b}}_{0}+3{\bf {b}}_{1})t+{\bf {b}}_{0}\ ,\ t\in [0,1]\ .\end{aligned}}}

Mit Hilfe des De-Casteljau-Algorithmus ausgedrückt: b ( t ) = b 0 3 ( t ) {\displaystyle {\bf {b}}(t)={\bf {b}}_{0}^{3}(t)}

b 0 3 ( t ) = ( 1 t ) { ( 1 t ) [ ( 1 t ) b 0 + t b 1 ] b 0 1 ( t ) + t [ ( 1 t ) b 1 + t b 2 ] b 1 1 ( t ) } b 0 2 ( t ) + t { ( 1 t ) [ ( 1 t ) b 1 + t b 2 ] b 1 1 ( t ) + t [ ( 1 t ) b 2 + t b 3 ] b 2 1 ( t ) } b 1 2 ( t ) {\displaystyle {\begin{aligned}{\bf {b}}_{0}^{3}(t)=(1-t)&\underbrace {{\Bigl \{}(1-t)\overbrace {{\bigl [}(1-t){\bf {b}}_{0}+t{\bf {b}}_{1}{\bigr ]}} ^{{\bf {b}}_{0}^{1}(t)}+t\overbrace {{\bigl [}(1-t){\bf {b}}_{1}+t{\bf {b}}_{2}{\bigr ]}} ^{{\bf {b}}_{1}^{1}(t)}{\Bigr \}}} _{{\bf {b}}_{0}^{2}(t)}+\\t&\underbrace {{\Bigl \{}(1-t)\overbrace {{\bigl [}(1-t){\bf {b}}_{1}+t{\bf {b}}_{2}{\bigr ]}} ^{{\bf {b}}_{1}^{1}(t)}+t\overbrace {{\bigl [}(1-t){\bf {b}}_{2}+t{\bf {b}}_{3}{\bigr ]}} ^{{\bf {b}}_{2}^{1}(t)}{\Bigr \}}} _{{\bf {b}}_{1}^{2}(t)}\end{aligned}}}

Vergleich der Kurvendarstellungen

Vergleich zwischen Kurvendarstellungen

Polynomiale Kurven (d. h. die Koordinaten sind Polynome bzgl. t {\displaystyle t} ) lassen sich in der Monom-Basis, der Bernsteinbasis und mit Hilfe des De-Casteljau-Algorithmus (fortgesetzte lineare Interpolation) beschreiben. Da die Koeffizientenpunkte der Monom-Basis nicht viel über den Kurvenverlauf aussagen, entstanden die Darstellungen mit der Bernstein-Basis (Bézier-Kurven) und mit dem De-Casteljau-Algorithmus. Die letzten beiden haben allerdings auch Vor- und Nachteile. Der De-Casteljau-Algorithmus hat gegenüber der Bézierdarstellung bei der Berechnung der Punkte (Numerik) Vorteile. Bézierkurven lassen sich dafür durch die vielen formalen Eigenschaften (s. o.) der Bernstein-Polynome leichter theoretisch (z. B. Krümmung) untersuchen. Der numerische Nachteil der Bézier-Kurven (Auswertung der Bernstein-Polynome) lässt sich durch eine dem Horner-Schema ähnlichen Methode ausgleichen:[5]

function bezier_comp(degree: integer; coeff: r_array; t: real): real;
{Berechnet eine Komponente einer Bezier-Kurve. (Aus FARIN: Curves and Surfaces...)}
var
  i, n_choose_i: integer;  fact, t1, aux: real;
begin
  t1 := 1 - t;  fact := 1;  n_choose_i := 1;
  aux := coeff[0] * t1;
  for i := 1 to degree - 1 do
  begin
    fact := fact * t;
    n_choose_i := n_choose_i * (degree - i + 1) div i;
    aux := (aux + fact * n_choose_i * coeff[i]) * t1;
  end;
  aux := aux + fact * t * coeff[degree];
  bezier_comp := aux;
end;  {bezier_comp}

Ableitungen einer Bézier-Kurve

Mit Hilfe der Ableitungen der Bernsteinpolynome ergibt sich für die 1. Ableitung der Bézierkurve   b ( t ) = i = 0 n B i n ( t ) b i   {\displaystyle \ {\bf {b}}(t)=\sum _{i=0}^{n}B_{i}^{n}(t){\bf {b}}_{i}\ } :

d d t b ( t ) = n i = 0 n 1 B i n 1 ( t ) ( b i + 1 b i )   . {\displaystyle {\frac {d}{dt}}{\bf {b}}(t)=n\sum _{i=0}^{n-1}B_{i}^{n-1}(t)({\bf {b}}_{i+1}-{\bf {b}}_{i})\ .}

Lässt man die Tangentenvektoren alle im Nullpunkt des Koordinatensystems beginnen, so beschreiben sie eine weitere Bézierkurve mit den Kontrollpunkten   Δ b i := b i + 1 b i   {\displaystyle \ \Delta {\bf {b}}_{i}:={\bf {b}}_{i+1}-{\bf {b}}_{i}\ } .

Speziell gilt:

d d t b ( 0 ) = n ( b 1 b 0 ) {\displaystyle {\frac {d}{dt}}{\bf {b}}(0)=n({\bf {b}}_{1}-{\bf {b}}_{0})\quad } und d d t b ( 1 ) = n ( b n b n 1 ) . {\displaystyle \quad {\frac {d}{dt}}{\bf {b}}(1)=n({\bf {b}}_{n}-{\bf {b}}_{n-1}).}

Um höhere Ableitungen übersichtlich schreiben zu können, führt man folgenden Differenzenoperator ein:[6]

Δ r b i := Δ r 1 b i + 1 Δ r 1 b i , {\displaystyle \Delta ^{r}{\bf {b}}_{i}:=\Delta ^{r-1}{\bf {b}}_{i+1}-\Delta ^{r-1}{\bf {b}}_{i},}

Es ist

Δ 0 b i = b i Δ 1 b i = b i + 1 b i Δ 2 b i = b i + 2 2 b i + 1 + b i Δ 3 b i = b i + 3 3 b i + 2 + 3 b i + 1 b i Δ r b i = j = 0 r ( r j ) ( 1 ) r j b i + j {\displaystyle {\begin{array}{rcl}\Delta ^{0}{\bf {b}}_{i}&=&{\bf {b}}_{i}\\\Delta ^{1}{\bf {b}}_{i}&=&{\bf {b}}_{i+1}-{\bf {b}}_{i}\\\Delta ^{2}{\bf {b}}_{i}&=&{\bf {b}}_{i+2}-2{\bf {b}}_{i+1}+{\bf {b}}_{i}\\\Delta ^{3}{\bf {b}}_{i}&=&{\bf {b}}_{i+3}-3{\bf {b}}_{i+2}+3{\bf {b}}_{i+1}-{\bf {b}}_{i}\\&\vdots \\\Delta ^{r}{\bf {b}}_{i}&=&\sum _{j=0}^{r}{r \choose j}(-1)^{r-j}{\bf {b}}_{i+j}\end{array}}}

Die r {\displaystyle r} -te Ableitung der Bézierkurve   b ( t ) = i = 0 n b i B i n ( t )   {\displaystyle \ {\bf {b}}(t)=\sum _{i=0}^{n}{\bf {b}}_{i}B_{i}^{n}(t)\ } lässt sich jetzt wie folgt schreiben:

d r d t r b ( t ) = n ! ( n r ) ! i = 0 n r Δ r b i B i n r ( t ) {\displaystyle {\frac {d^{r}}{dt^{r}}}{\bf {b}}(t)={\frac {n!}{(n-r)!}}\cdot \sum _{i=0}^{n-r}\Delta ^{r}{\bf {b}}_{i}B_{i}^{n-r}(t)}

Speziell für t = 0 {\displaystyle t=0} und t = 1 {\displaystyle t=1} erhält man

d r d t r b ( 0 ) = n ! ( n r ) ! Δ r b 0 {\displaystyle {\frac {d^{r}}{dt^{r}}}{\bf {b}}(0)={\frac {n!}{(n-r)!}}\Delta ^{r}{\bf {b}}_{0}\quad } und d r d t r b ( 1 ) = n ! ( n r ) ! Δ r b n r {\displaystyle \quad {\frac {d^{r}}{dt^{r}}}{\bf {b}}(1)={\frac {n!}{(n-r)!}}\Delta ^{r}{\bf {b}}_{n-r}}

Graderhöhung einer Bézierkurve

Eine wichtige Manipulation der Darstellung einer vorgegebenen Bézierkurve ist die sog. Graderhöhung. Sie ist vergleichbar mit dem Anfügen von Termen 0 t n + 1 , 0 t n + 2 , {\displaystyle 0t^{n+1},0t^{n+2},\dots } an ein Polynom   a 0 + a 1 t + a n t n   {\displaystyle \ a_{0}+a_{1}t+\cdots a_{n}t^{n}\ } . Dabei ändert sich das Polynom nicht und der (scheinbare) Grad wird erhöht. Analog stellt man eine fest vorgegebene Bézierkurve   b ( t ) = i = 0 n b i B i n ( t )   {\displaystyle \ {\bf {b}}(t)=\sum _{i=0}^{n}{\bf {b}}_{i}B_{i}^{n}(t)\ } in der Form   b ( t ) = i = 0 n + 1 b i ( 1 ) B i n + 1 ( t )   {\displaystyle \ {\bf {b}}(t)=\sum _{i=0}^{n+1}{\bf {b}}_{i}^{(1)}B_{i}^{n+1}(t)\ } mit geeigneten neuen Kontrollpunkten   b 0 ( 1 ) , , b n + 1 ( 1 )   {\displaystyle \ {\bf {b}}_{0}^{(1)},\dots ,{\bf {b}}_{n+1}^{(1)}\ } dar. Um die neuen Kontrollpunkte zu bestimmen, multipliziert man die ursprüngliche Darstellung mit dem Faktor   t + ( 1 t )   {\displaystyle \ t+(1-t)\ } :

b ( t ) = i = 0 n B i n ( t ) ( t + ( 1 t ) ) b i = i = 0 n ( n i ) t i ( 1 t ) n i ( t + ( 1 t ) ) b i = i = 0 n ( n i ) t i + 1 ( 1 t ) n i b i + i = 0 n ( n i ) t i ( 1 t ) n i + 1 b i = i = 1 n + 1 ( n i 1 ) t i ( 1 t ) n + 1 i b i 1 + i = 0 n ( n i ) t i ( 1 t ) n i + 1 b i = i = 0 n + 1 t i ( 1 t ) n + 1 i ( ( n i 1 ) b i 1 + ( n i ) b i ) = i = 0 n + 1 ( n + 1 i ) t i ( 1 t ) n + 1 i b i ( 1 ) = i = 0 n + 1 B i n + 1 ( t ) b i ( 1 ) , wobei b i ( 1 ) := i n + 1 b i 1 + ( 1 i n + 1 ) b i ,   i = 0 , , n + 1 = ( i n + 1 ) b i 1 für  1 i n  (s. Casteljau-Algor.) {\displaystyle {\begin{array}{rcl}{\bf {b}}(t)&=&\sum _{i=0}^{n}B_{i}^{n}(t)(t+(1-t)){\bf {b}}_{i}\\&=&\sum _{i=0}^{n}{n \choose i}t^{i}(1-t)^{n-i}(t+(1-t)){\bf {b}}_{i}\\&=&\sum _{i=0}^{n}{n \choose i}t^{i+1}(1-t)^{n-i}{\bf {b}}_{i}+\sum _{i=0}^{n}{n \choose i}t^{i}(1-t)^{n-i+1}{\bf {b}}_{i}\\&=&\sum _{i=1}^{n+1}{n \choose i-1}t^{i}(1-t)^{n+1-i}{\bf {b}}_{i-1}+\sum _{i=0}^{n}{n \choose i}t^{i}(1-t)^{n-i+1}{\bf {b}}_{i}\\&=&\sum _{i=0}^{n+1}t^{i}(1-t)^{n+1-i}\cdot \left({n \choose i-1}{\bf {b}}_{i-1}+{n \choose i}{\bf {b}}_{i}\right)\\&=&\sum _{i=0}^{n+1}{n+1 \choose i}t^{i}(1-t)^{n+1-i}{\bf {b}}_{i}^{(1)}\\&=&\sum _{i=0}^{n+1}B_{i}^{n+1}(t){\bf {b}}_{i}^{(1)},\quad {\text{wobei}}\\{\bf {b}}_{i}^{(1)}&:=&{\frac {i}{n+1}}{\bf {b}}_{i-1}+\left(1-{\frac {i}{n+1}}\right){\bf {b}}_{i},\ i=0,\dots ,n+1\\&=&\left({\frac {i}{n+1}}\right){\bf {b}}_{i}^{1}\quad {\text{für }}1\leq i\leq n{\text{ (s. Casteljau-Algor.)}}\end{array}}}

Die neuen Kontrollpunkte sind also:[7]

b 0 , ( 1 n + 1 ) b 1 1 , ( 2 n + 1 ) b 2 1 , , ( n n + 1 ) b n 1 , b n . {\displaystyle {\bf {b}}_{0},\quad \left({\frac {1}{n+1}}\right){\bf {b}}_{1}^{1},\quad \left({\frac {2}{n+1}}\right){\bf {b}}_{2}^{1}\quad ,\ldots ,\quad \left({\frac {n}{n+1}}\right){\bf {b}}_{n}^{1},\quad {\bf {b}}_{n}.}
Graderhöhungen einer Bézier-Kurve:
Kontrollpolygone bei ein-, zwei- und 10-maliger Graderhöhung}

Wesentliche Eigenschaften der Graderhöhung sind:

  1. Wiederholte Graderhöhung führt zu einer Approximation der Bézierkurve durch das Kontrollpolygon.
  2. Die größere Anzahl von Kontrollpunkten bietet mehr Freiheitsgrade, die Kurve zu verändern.
  3. Mehrere Bézierkurven lassen sich auf einen einheitlichen Grad bringen. Dies ist wichtig bei Tensorprodukt-Bézierflächen.
  4. Damit lassen sich auch dann quadratische Bézierkurven als kubische darstellen, falls ein Vektorzeichenprogramm (z. B. Inkscape) bzw. eine Grafikbibliothek (z. B. Cairo) nur kubische unterstützt.

Teilung einer Bézierkurve

Zerlegung einer Bézierkurve

Eine Bézierkurve Γ 0 : b ( t ) {\displaystyle \Gamma _{0}:{\bf {b}}(t)} ist normalerweise definiert für 0 t 1 {\displaystyle 0\leq t\leq 1} . Sei nun c ( 0 , 1 ) {\displaystyle c\in (0,1)} . Dann ist Γ 1 : b ( t ) {\displaystyle \Gamma _{1}:{\bf {b}}(t)} mit 0 t c {\displaystyle 0\leq t\leq c} ein Teil der gegebenen Bézierkurve. Es soll nun die Teilkurve Γ 1 {\displaystyle \Gamma _{1}} als Bézierkurve c ( u ) = i = 0 n c i B i n ( u ) {\displaystyle {\bf {c}}(u)=\sum _{i=0}^{n}{\bf {c}}_{i}B_{i}^{n}(u)} mit u [ 0 , 1 ] {\displaystyle u\in [0,1]} vom (selben) Grad n {\displaystyle n} mit geeigneten Kontrollpunkten c 0 , c 1 , , c n {\displaystyle {\bf {c}}_{0},{\bf {c}}_{1},\dots ,{\bf {c}}_{n}} dargestellt werden. Setzt man u := t c {\displaystyle u:={\frac {t}{c}}} , so muss die folgende Gleichung erfüllt sein:

r = 0 n c r B r n ( t c ) = i = 0 n b i B i n ( t ) , {\displaystyle \sum _{r=0}^{n}{\bf {c}}_{r}B_{r}^{n}\left({\frac {t}{c}}\right)=\sum _{i=0}^{n}{\bf {b}}_{i}B_{i}^{n}(t),\quad } für t [ 0 , c ] {\displaystyle \quad t\in [0,c]}

Dies gilt für[8][9]

  • c r = b 0 r ( c ) {\displaystyle \quad {\bf {c}}_{r}={\bf {b}}_{0}^{r}(c)\quad } (s. Casteljau-Alg. und Abbildung)

Denn

r = 0 n c r B r n ( t c ) = r = 0 n i = 0 r b i B i r ( c ) B r n ( t c ) {\displaystyle \sum _{r=0}^{n}{\bf {c}}_{r}B_{r}^{n}\left({\frac {t}{c}}\right)=\sum _{r=0}^{n}\sum _{i=0}^{r}{\bf {b}}_{i}B_{i}^{r}(c)B_{r}^{n}\left({\frac {t}{c}}\right)}
= r = 0 n i = 0 n b i B i r ( c ) B r n ( t c ) {\displaystyle =\sum _{r=0}^{n}\sum _{i=0}^{n}{\bf {b}}_{i}B_{i}^{r}(c)B_{r}^{n}\left({\frac {t}{c}}\right)\quad } wegen B i r ( ) = 0 für i > r {\displaystyle \quad B_{i}^{r}(\cdots )=0\quad {\text{für}}\quad i>r}
= i = 0 n b i r = 0 n B i r ( c ) B r n ( t c ) {\displaystyle =\sum _{i=0}^{n}{\bf {b}}_{i}\sum _{r=0}^{n}B_{i}^{r}(c)B_{r}^{n}\left({\frac {t}{c}}\right)}
= i = 0 n b i B i n ( t ) {\displaystyle =\sum _{i=0}^{n}{\bf {b}}_{i}B_{i}^{n}(t)\quad } (s. Eigensch. (S) der Bernst.-Pol.)

Der restliche Bogen ist die Bézierkurve Γ 2 : d ( v ) {\displaystyle \Gamma _{2}:{\bf {d}}(v)} mit den Kontrollpunkten

  • d r = b r n r ( c ) {\displaystyle \quad {\bf {d}}_{r}={\bf {b}}_{r}^{n-r}(c)\quad } (s. Abbildung)

Rationale Bézierkurven

Rationale Kurven und projektive Kurven

Bézierkurven sind parametrisierte Kurven, deren Parameterdarstellungen nur Polynome verwenden. Leider lassen sich so wichtige und geometrisch einfache Kurven wie Kreise nicht durch polynomiale Parameterdarstellungen beschreiben. Dieser Nachteil ist u. a. das Motiv für die Erweiterung der als Parameterfunktionen zulässigen Funktionen auf rationale Funktionen. Denn jeder Kegelschnitt hat eine rationale Darstellung. Da eine Kurve mit einer rationalen Darstellung

( p 1 ( t ) q 1 ( t ) , p 2 ( t ) q 2 ( t ) ) T {\displaystyle \left({\frac {p_{1}(t)}{q_{1}(t)}},{\frac {p_{2}(t)}{q_{2}(t)}}\right)^{T}}

wobei die Funktionen p i {\displaystyle p_{i}} und q i {\displaystyle q_{i}} Polynome sind, in homogenen Koordinaten die polynomiale Darstellung

( p 1 ( t ) q 2 ( t ) , p 2 ( t ) q 1 ( t ) , q 1 ( t ) q 2 ( t ) ) T {\displaystyle \langle \left(p_{1}(t)q_{2}(t)\;,\;p_{2}(t)q_{1}(t)\;,\;q_{1}(t)q_{2}(t)\right)^{T}\rangle }

besitzt, lassen sich ebene Kurven mit rationalen Koeffizientenfunktionen als Zentralprojektion einer Bézierkurve im R 3 {\displaystyle \mathbb {R} ^{3}} auf die Einbettungsebene x 3 = 1 {\displaystyle x_{3}=1} auffassen.

Die analoge Aussage gilt für Kurven im R 3 {\displaystyle \mathbb {R} ^{3}} . Sie lassen sich als Zentralprojektion einer Bézierkurve im R 4 {\displaystyle \mathbb {R} ^{4}} auf den Einbettungsraum x 4 = 1 {\displaystyle x_{4}=1} auffassen. Damit lassen sich die Vorteile der Bézier-Darstellung einer polynomialen Kurve auch für rationale Kurven nutzen.

Rationale Bézierkurve (rot) als Projektion einer gewöhnlichen räumlichen Bézierkurve (blau). Die Kontrollpunkte a i {\displaystyle {\bf {a}}_{i}} sind alle eigentlich, da keiner der Kontrollpunkte b i {\displaystyle {\bf {b}}_{i}} in der x 1 {\displaystyle x_{1}} - x 2 {\displaystyle x_{2}} -Ebene liegt.

Ebene rationale Bézierkurven

Es sei nun n > 0 {\displaystyle n>0} festgewählt und die Vektoren b 0 , b 1 , , b n ,   b i 0 , {\displaystyle {\bf {b}}_{0},{\bf {b}}_{1},\cdots ,{\bf {b}}_{n},\ {\bf {b}}_{i}\neq {\bf {0}},} beschreiben ein Polygon im R 3 {\displaystyle \mathbb {R} ^{3}} . Dann ist

x ( t ) := b 0 B 0 n ( t ) + b 1 B 1 n ( t ) + + b n B n n ( t ) , 0 t 1 , {\displaystyle {\bf {x}}(t):={\bf {b}}_{0}B_{0}^{n}(t)+{\bf {b}}_{1}B_{1}^{n}(t)+\cdots +{\bf {b}}_{n}B_{n}^{n}(t),\quad 0\leq t\leq 1,}

eine (räumliche) Bézier-Kurve vom Grad n {\displaystyle n} . Die Punkte b 0 , , b n {\displaystyle {\bf {b}}_{0},\cdots ,{\bf {b}}_{n}} sind die Kontrollpunkte der (räumlichen) Bézierkurve. Fasst man die 1-dimensionalen Unterräume

b 0 B 0 n ( t ) + b 1 B 1 n ( t ) + + b n B n n ( t ) {\displaystyle \langle {\bf {b}}_{0}B_{0}^{n}(t)+{\bf {b}}_{1}B_{1}^{n}(t)+\cdots +{\bf {b}}_{n}B_{n}^{n}(t)\rangle }

als Punkte der reellen projektiven Ebene mit der Ferngerade x 3 = 0 {\displaystyle x_{3}=0} auf, so bezeichnet man den affinen Anteil (Projektion vom Nullpunkt aus auf die Ebene x 3 = 1 {\displaystyle x_{3}=1} ) dieser projektiven Kurve als rationale Bézierkurve.

Die Kontrollpunkte der Bézierkurve im R 3 {\displaystyle \mathbb {R} ^{3}} lassen sich folgendermaßen beschreiben:

b i = w i ( x i , y i , 1 ) T ,   w i > 0 ,   {\displaystyle {\bf {b}}_{i}=w_{i}(x_{i},y_{i},1)^{T},\ w_{i}>0,\ } falls b i {\displaystyle {\bf {b}}_{i}} nicht auf der Ferngerade x 3 = 0 {\displaystyle x_{3}=0} und
b i = w i ( x i , y i , 0 ) T ,   w i > 0 ,   {\displaystyle {\bf {b}}_{i}=w_{i}(x_{i},y_{i},0)^{T},\ w_{i}>0,\ } falls b i {\displaystyle {\bf {b}}_{i}} auf der Ferngerade liegt.

Beim Übergang zu inhomogenen Koordinaten wird ein Kontrollpunkt entweder auf den affinen Punkt a i := ( x i , y i ) T {\displaystyle {\bf {a}}_{i}:=(x_{i},y_{i})^{T}} oder auf den Fernpunkt in Richtung a i := ( x i , y i ) T {\displaystyle {\bf {a}}_{i}:=(x_{i},y_{i})^{T}} abgebildet. Der Punkt a i {\displaystyle {\bf {a}}_{i}} heißt eigentlicher bzw. uneigentlicher Kontrollpunkt der rationalen Bézierkurve und die Zahl w i {\displaystyle w_{i}} heißt das Gewicht des Kontrollpunktes a i {\displaystyle {\bf {a}}_{i}} . Eine rationale Bézierkurve hat folgende affine Beschreibung:[10][11]

  • i = 0 n w i a i B i n ( t ) i = 0 n δ i w i B i n ( t ) , {\displaystyle \qquad {\frac {\sum _{i=0}^{n}w_{i}{\bf {a}}_{i}B_{i}^{n}(t)}{\sum _{i=0}^{n}\delta _{i}w_{i}B_{i}^{n}(t)}},}

wobei δ i = 1 {\displaystyle \delta _{i}=1} für eigentliche und δ i = 0 {\displaystyle \delta _{i}=0} für uneigentliche Kontrollpunkte zu setzen ist.

Die rationalen Bézierkurven haben (u. a.) die folgenden Eigenschaften:

Sind a i , w i , i = 0 , , n , {\displaystyle {\bf {a}}_{i},w_{i},i=0,\dots ,n,} eigentliche Kontrollpunkte bzw. die Gewichte einer rationalen Bézierkurve Γ {\displaystyle \Gamma } , so gilt

  1. Die Kurve Γ {\displaystyle \Gamma } enthält die Kontrollpunkte a 0 , a n {\displaystyle {\bf {a}}_{0},{\bf {a}}_{n}} (erster bzw. letzter Punkt des Kontrollpolygons).
  2. Die Tangente im Punkt a 0 {\displaystyle {\bf {a}}_{0}} bzw. a n {\displaystyle {\bf {a}}_{n}} hat die Richtung a 1 a 0 {\displaystyle {\bf {a}}_{1}-{\bf {a}}_{0}} bzw. a n a n 1 {\displaystyle {\bf {a}}_{n}-{\bf {a}}_{n-1}} .
  3. Eine Erhöhung des Gewichts w i {\displaystyle w_{i}} bewirkt eine Veränderung der Kurve auf den Kontrollpunkt a i {\displaystyle {\bf {a}}_{i}} zu. (s. Abbildung)
Rationale Bézierkurven mit Gewichten w 0 = w 1 = w 3 = 1 {\displaystyle w_{0}=w_{1}=w_{3}=1} und verschiedenen Gewichten w 2 {\displaystyle w_{2}}

Zusammenfassung:
Eine ebene rationale Bézierkurve besitzt neben dem Kontrollpolygon noch die Gewichte w i {\displaystyle w_{i}} als Designparameter. Will man eine Kurve erzeugen, legt man zunächst die Kontrollpunkte a i {\displaystyle {\bf {a}}_{i}} und die Gewichte w i {\displaystyle w_{i}} fest. Dadurch wird dann auch eine räumliche (gewöhnliche) Bézierkurve mit den Kontrollpunkten b 0 , b 1 , , b n {\displaystyle {\bf {b}}_{0},{\bf {b}}_{1},\cdots ,{\bf {b}}_{n}} definiert. Die Projektion dieser Kurve (vom Nullpunkt aus) auf die x-y-Ebene ( x 3 = 1 {\displaystyle x_{3}=1} ) liefert dann die ebene rationale Bézierkurve. Eine Variation der Gewichte ändert nicht die Kontrollpunkte a i {\displaystyle {\bf {a}}_{i}} , aber die (räumlichen) Kontrollpunkte b i {\displaystyle {\bf {b}}_{i}} und damit die zugehörige räumliche Bézierkurve und schließlich auch die (ebene) rationale Bézierkurve. Erhöht man ein Gewicht w i {\displaystyle w_{i}} , so entfernt sich der zugehörige Kontrollpunkt b i {\displaystyle {\bf {b}}_{i}} vom Nullpunkt und zieht die räumliche Bézierkurve mit. Der zugehörige Kontrollpunkt a i {\displaystyle {\bf {a}}_{i}} dagegen bleibt unverändert. Die rationale Bézierkurve bewegt sich auf ihn zu (s. Bild). Verringert man das Gewicht, bewegt sich die Kurve von dem Kontrollpunkt weg. Falls alle Gewichte 1 sind, ist die rationale Bézierkurve eine gewöhnliche Bézierkurve mit den Kontrollpunkten a i {\displaystyle {\bf {a}}_{i}} .

Kegelschnitte als rationale Bézierkurven

Parabel:
Eine Bézierkurve vom Grad zwei mit nicht kollinearen Kontrollpunkten b 0 , b 1 , b 2 {\displaystyle {\bf {b}}_{0},{\bf {b}}_{1},{\bf {b}}_{2}} im R 3 / R 2 {\displaystyle \mathbb {R} ^{3}/\mathbb {R} ^{2}} ist immer eine Parabel (s. oben). Um eine Parabel als (ganz-)rationale Bézierkurve darzustellen, wählt man drei nicht kollineare Kontrollpunkte a 0 , a 1 , a 2 {\displaystyle {\bf {a}}_{0},{\bf {a}}_{1},{\bf {a}}_{2}} und setzt w 0 = w 1 = w 2 = 1 {\displaystyle w_{0}=w_{1}=w_{2}=1} und δ 0 = δ 1 = δ 2 = 1 {\displaystyle \delta _{0}=\delta _{1}=\delta _{2}=1} . Letzteres bedeutet: Die Kontrollpunkte sind alle eigentlich.

Ellipsen und Hyperbeln lassen sich durch Zentralprojektion von Parabeln im R 3 {\displaystyle \mathbb {R} ^{3}} , deren Ebenen nicht den Nullpunkt enthalten, auf die Einbettungsebene x 3 = 1 {\displaystyle x_{3}=1} erzeugen.

Hyperbel (rot) als rationale Bézierkurve mit den Kontrollpunkten a 0 , a 1 , a 2 {\displaystyle {\bf {a}}_{0},{\bf {a}}_{1},{\bf {a}}_{2}} . Sie entsteht durch Projektion der blauen Parabel vom Nullpunkt aus auf die Ebene x 3 = 1 {\displaystyle x_{3}=1} . Die Parabelpunkte b 0 , b 2 {\displaystyle {\bf {b}}_{0},{\bf {b}}_{2}} gehen dabei auf die Fernpunkte der Hyperbelasymptoten über.

Hyperbel:
Für die Kontrollpunkte

  b 0 = ( 1 , 0 , 0 ) T   ,   b 1 = ( 0 , 0 , 1 / 2 ) T   ,   b 2 = ( 0 , 1 , 0 ) T   {\displaystyle \ {\bf {b}}_{0}=(1,0,0)^{T}\ ,\ {\bf {b}}_{1}=(0,0,1/2)^{T}\ ,\ {\bf {b}}_{2}=(0,1,0)^{T}\ }

beschreibt

b ( t ) = t 2 b 0 + 2 t ( 1 t ) b 1 + ( 1 t ) 2 b 2 {\displaystyle {\bf {b}}(t)=t^{2}{\bf {b}}_{0}+2t(1-t){\bf {b}}_{1}+(1-t)^{2}{\bf {b}}_{2}}
= ( t 2 , ( 1 t ) 2 , t ( 1 t ) ) T , t R , {\displaystyle =\left(t^{2},(1-t)^{2},t(1-t)\right)^{T},\quad t\in \mathbb {R} ,}

eine Parabel, die auf dem Kegel mit der Gleichung   x 1 x 2 = x 3 2   {\displaystyle \ x_{1}x_{2}=x_{3}^{2}\ } liegt (s. Bild). Die Kontrollpunkte und Gewichte der zugehörigen (ebenen) rationalen Bézierkurve sind:

  a 0 = ( 1 , 0 ) T   ,   a 1 = ( 0 , 0 ) T   ,   a 2 = ( 0 , 1 ) T {\displaystyle \ {\bf {a}}_{0}=(1,0)^{T}\ ,\ {\bf {a}}_{1}=(0,0)^{T}\ ,\ {\bf {a}}_{2}=(0,1)^{T}\quad } bzw. :   w 0 = 1 , w 1 = 1 / 2 , w 2 = 1   {\displaystyle \ w_{0}=1\;,\;w_{1}=1/2\;,\;w_{2}=1\ } .

a 0 , a 2 {\displaystyle {\bf {a}}_{0},{\bf {a}}_{2}} sind uneigentliche Kontrollpunkte. Damit ist   δ 0 = 0 , δ 1 = 1 , δ 2 = 0   {\displaystyle \ \delta _{0}=0\;,\;\delta _{1}=1\;,\;\delta _{2}=0\ } und der Nenner (s. o.) der rationalen Komponenten ist   t ( 1 t )   {\displaystyle \ t(1-t)\ } .

Also ist die zugehörige rationale Bézierkurve

a ( t ) = t 2 a 0 + 1 2 2 t ( 1 t ) a 1 + ( 1 t ) 2 a 2 t ( 1 t ) {\displaystyle {\bf {a}}(t)={\frac {t^{2}{\bf {a}}_{0}+{\tfrac {1}{2}}2t(1-t){\bf {a}}_{1}+(1-t)^{2}{\bf {a}}_{2}}{t(1-t)}}}
= ( t 2 , ( 1 t ) 2 ) T t ( 1 t ) = ( t 1 t , 1 t t ) T {\displaystyle ={\frac {\left(t^{2},(1-t)^{2}\right)^{T}}{t(1-t)}}=\left({\frac {t}{1-t}},{\frac {1-t}{t}}\right)^{T}}

Dies ist eine rationale Parameterdarstellung eines Astes der Hyperbel mit der Gleichung   y = 1 / x   {\displaystyle \ y=1/x\ } .

Die Änderung   a 0 = ( c , 0 ) T   {\displaystyle \ {\bf {a}}_{0}=(c,0)^{T}\ } liefert eine rationale Bézierdarstellung der Hyperbel y = c / x {\displaystyle y=c/x} .

Halbkreis (rot) als rationale Bézierkurve. a 1 {\displaystyle {\bf {a}}_{1}} ist uneigentlicher Kontrollpunkt, d. h. hier: Der Vektor a 1 {\displaystyle {\bf {a}}_{1}} gibt die Richtung der Tangenten in den Kreispunkten a 0 , a 2 {\displaystyle {\bf {a}}_{0},{\bf {a}}_{2}} an.

Kreis:
In dem folgenden Beispiel sind die Kontrollpunkte der (räumlichen) Parabel:

b 0 = ( 1 , 0 , 1 ) T , b 1 = ( 0 , 1 , 0 ) T , b 2 ( 1 , 0 , 1 ) T   {\displaystyle {\bf {b}}_{0}=(1,0,1)^{T}\;,\;{\bf {b}}_{1}=(0,1,0)^{T}\;,\;{\bf {b}}_{2}(-1,0,1)^{T}\ } .

Die Bézierkurve

b ( t ) = t 2 b 0 + 2 t ( 1 t ) b 1 + ( 1 t ) 2 b 2 {\displaystyle {\bf {b}}(t)=t^{2}{\bf {b}}_{0}+2t(1-t){\bf {b}}_{1}+(1-t)^{2}{\bf {b}}_{2}}
= ( t 2 ( 1 t ) 2 , 2 t ( 1 t ) , t 2 + ( 1 t ) 2 ) T {\displaystyle =\left(t^{2}-(1-t)^{2},2t(1-t),t^{2}+(1-t)^{2}\right)^{T}}

liegt in diesem Fall auf dem Kegel mit der Gleichung   x 1 2 + x 2 2 = x 3 2   {\displaystyle \ x_{1}^{2}+x_{2}^{2}=x_{3}^{2}\ } (s. Abbildung). Die Kontrollpunkte und Gewichte der zu gehörigen rationalen Bézierkurve sind:

  a 0 = ( 1 , 0 ) T   ,   a 1 = ( 0 , 1 ) T   ,   a 2 = ( 1 , 0 ) T {\displaystyle \ {\bf {a}}_{0}=(1,0)^{T}\ ,\ {\bf {a}}_{1}=(0,1)^{T}\ ,\ {\bf {a}}_{2}=(-1,0)^{T}\quad } bzw. w 0 = w 1 = w 2 = 1   {\displaystyle \quad w_{0}=w_{1}=w_{2}=1\ } .

a 1 {\displaystyle {\bf {a}}_{1}} ist uneigentlicher Kontrollpunkt. Damit ist   δ 0 = 1 , δ 1 = 0 , δ 2 = 1   {\displaystyle \ \delta _{0}=1\;,\;\delta _{1}=0\;,\;\delta _{2}=1\ } und der Nenner (s. o.) der rationalen Komponenten ist   t 2 + ( 1 t ) 2   {\displaystyle \ t^{2}+(1-t)^{2}\ } .

Also ist die zugehörige rationale Bézierkurve

a ( t ) = t 2 a 0 + 2 t ( 1 t ) a 1 + ( 1 t ) 2 a 2 t 2 + ( 1 t ) 2 {\displaystyle {\bf {a}}(t)={\frac {t^{2}{\bf {a}}_{0}+2t(1-t){\bf {a}}_{1}+(1-t)^{2}{\bf {a}}_{2}}{t^{2}+(1-t)^{2}}}}
= ( t 2 ( 1 t ) 2 , 2 t ( 1 t ) ) T t 2 + ( 1 t ) 2 {\displaystyle ={\frac {\left(t^{2}-(1-t)^{2},2t(1-t)\right)^{T}}{t^{2}+(1-t)^{2}}}}
= ( t 2 ( 1 t ) 2 t 2 + ( 1 t ) 2 , 2 t ( 1 t ) t 2 + ( 1 t ) 2 ) T {\displaystyle =\left({\frac {t^{2}-(1-t)^{2}}{t^{2}+(1-t)^{2}}}\;,\;{\frac {2t(1-t)}{t^{2}+(1-t)^{2}}}\right)^{T}}

Für 0 t 1 {\displaystyle 0\leq t\leq 1} ist dies eine rationale Parameterdarstellung eines halben Einheitskreises.

Setzt man   a 0 = ( a , 0 ) T   ,   a 1 = ( 0 , b ) T   ,   a 2 = ( a , 0 ) T , a > 0 , b > 0 ,   {\displaystyle \ {\bf {a}}_{0}=(a,0)^{T}\ ,\ {\bf {a}}_{1}=(0,b)^{T}\ ,\ {\bf {a}}_{2}=(-a,0)^{T},\;a>0,b>0,\ } erhält man eine rationale Bézierdarstellung der Ellipse mit der Gleichung   x 2 a 2 + y 2 b 2 = 1   {\displaystyle \ {\frac {x^{2}}{a^{2}}}+{\frac {y^{2}}{b^{2}}}=1\ } .

Anwendung: Kreisapproximation durch kubische Bézierkurven

Kreise bzw. Kreisbögen lassen sich durch Bézierkurven nicht exakt, sondern nur genähert darstellen. Eine solche Näherung ist z. B. für die Gestaltung einer Typ-1-PostScript-Schrift nötig, da hier nur Strecken und Bézierkurven dritten Grades erlaubt sind. (Jedoch verläuft auch für größere n {\displaystyle n} keine Bézierkurve t ( x ( t ) | y ( t ) ) {\displaystyle t\mapsto (x(t)|y(t))} n {\displaystyle n} -ten Grades in einem noch so kleinen Kreisbogen mit Radius r {\displaystyle r} zum Mittelpunkt ( z 1 | z 2 ) {\displaystyle (z_{1}|z_{2})} , denn ( x ( t 0 ) | y ( t 0 ) ) {\displaystyle (x(t_{0})|y(t_{0}))} liegt genau dann auf dem Kreisbogen, wenn t 0 {\displaystyle t_{0}} Nullstelle der Polynomfunktion t [ x ( t ) z 1 ] 2 + [ y ( t ) z 2 ] 2 r 2 {\displaystyle t\mapsto [x(t)-z_{1}]^{2}+[y(t)-z_{2}]^{2}-r^{2}} vom Grad 2 n {\displaystyle 2n} ist, was höchstens 2 n {\displaystyle 2n} Male vorkommt – vgl. Fehleranalyse.)

Teilt man einen Kreis in nur zwei (gleich große) Segmente und nähert die Halbkreise durch kubische Bézierkurven, zeigen sich größere Abweichung von der Kreisgestalt. Durch eine feinere Unterteilung in mehr Segmente lässt sich ein Kreis besser nähern. Je geringer der überstrichene Winkelbereich des Kreissegments ist, desto genauer ist die Näherung durch die Bézierkurve. Eine oft verwendete, einfache Realisierung eines Kreises verwendet vier Viertelkreisbögen, die als kubische Bézierkurven dargestellt werden. Um die Verbesserung der Näherung durch Verfeinerung der Unterteilung zu demonstrieren, werden in der Folge die Fehler der Halbkreisapproximation und der Viertelkreisapproximation miteinander verglichen.

Notation: Wir untersuchen Approximationen eines Kreises Q {\displaystyle Q} mit folgenden Parametern:

  • r {\displaystyle r} ist der Radius von Q {\displaystyle Q}
  • M {\displaystyle M} ist der Mittelpunkt von Q {\displaystyle Q}
  • die Kontrollpunkte P 0 {\displaystyle P_{0}} und P 3 {\displaystyle P_{3}} liegen vom Mittelpunkt M {\displaystyle M} im Abstand r {\displaystyle r} entfernt (also auf der Kreislinie von Q {\displaystyle Q} )
  • κ {\displaystyle \kappa } ist eine reelle Zahl zwischen 0 und 1 ( κ = 1 {\displaystyle \kappa =1} entspräche einer quadratischen Bézierapproximation).

Die zusätzlichen Kontrollpunkte P 1 {\displaystyle P_{1}} und P 2 {\displaystyle P_{2}} werden so gewählt, dass P 1 {\displaystyle P_{1}} zu P 0 {\displaystyle P_{0}} und P 2 {\displaystyle P_{2}} zu P 3 {\displaystyle P_{3}} den Abstand κ r {\displaystyle \kappa r} hat.

Beispielkoordinaten Viertelkreis:

Als einfaches Beispiel einer Viertelkreisapproximation wählen wir:

  • den Mittelpunkt M {\displaystyle M} des Kreises Q {\displaystyle Q} als ( 0 | 0 ) {\displaystyle (0|0)} ,
  • den Kontrollpunkt P 0 {\displaystyle P_{0}} auf der Kreislinie als ( r | 0 ) {\displaystyle (r|0)} ,
  • den Kontrollpunkt P 3 {\displaystyle P_{3}} auf der Kreislinie als ( 0 | r ) {\displaystyle (0|r)} – die Strecke M P 3 {\displaystyle MP_{3}} steht also senkrecht auf M P 0 {\displaystyle MP_{0}} , so dass beide Strecken einen Viertelkreissektor bilden –,
  • den Kontrollpunkt P 1 {\displaystyle P_{1}} als P 0 + κ ( P 3 M ) = ( r | r κ ) {\displaystyle P_{0}+\kappa (P_{3}-M)=(r|r\kappa )} (auf der Strecke P 0 ( r | r ) {\displaystyle P_{0}(r|r)} ),
  • den Kontrollpunkt P 2 {\displaystyle P_{2}} als P 3 + κ ( P 0 M ) = ( r κ | r ) {\displaystyle P_{3}+\kappa (P_{0}-M)=(r\kappa |r)} (auf der Strecke P 3 ( r | r ) {\displaystyle P_{3}(r|r)} ).

Die vier Kontrollpunkte liegen also auf dem Rand des Quadrats mit den Eckpunkten M = ( 0 | 0 ) {\displaystyle M=(0|0)} , P 0 = ( r | 0 ) {\displaystyle P_{0}=(r|0)} , P 3 = ( 0 | r ) {\displaystyle P_{3}=(0|r)} und ( r | r ) {\displaystyle (r|r)} . Dies gewährleistet immerhin, dass die Näherungskurve und die Kreislinie in P 0 {\displaystyle P_{0}} und P 3 {\displaystyle P_{3}} dieselbe Tangente haben. So ist auch die aus den Viertelkreisapproximationen zusammengesetzte Kurve in den Knotenpunkten „glatt“.

Die kubische Bézierkurve t ( x κ ( t ) | y κ ( t ) ) {\displaystyle t\mapsto (x_{\kappa }(t)|y_{\kappa }(t))} ( t [ 0 ; 1 ] {\displaystyle t\in [0;1]} ) hat mit diesen Kontrollpunkten folgende Form:

x κ ( t ) = ( 1 t ) 3 r + 3 t ( 1 t ) 2 r + 3 t 2 ( 1 t ) κ r , y κ ( t ) = t 3 r + 3 t 2 ( 1 t ) r + 3 t ( 1 t ) 2 κ r {\displaystyle x_{\kappa }(t)=(1-t)^{3}r+3t(1-t)^{2}r+3t^{2}(1-t)\kappa r\,,\quad y_{\kappa }(t)=t^{3}r+3t^{2}(1-t)r+3t(1-t)^{2}\kappa r}

Eine recht gute Approximation des oberen rechten Viertelkreisbogens erhält man mit κ = 0,552 {\displaystyle \kappa =0{,}552} , wie die nachfolgende Betrachtung zeigt.

Fehleranalyse:

Die Abweichung der gerade angegebenen Bézierkurve t ( x κ ( t ) | y κ ( t ) ) {\displaystyle t\mapsto (x_{\kappa }(t)|y_{\kappa }(t))} vom darzustellenden Kreis Q {\displaystyle Q} lässt sich folgendermaßen quantifizieren:

Ein Punkt ( x κ ( t 0 ) | y κ ( t 0 ) ) {\displaystyle (x_{\kappa }(t_{0})|y_{\kappa }(t_{0}))} der Bézierkurve ( x κ | y κ ) {\displaystyle (x_{\kappa }|y_{\kappa })} liegt genau dann auf der vorgegebenen Kreislinie mit Radius r {\displaystyle r} um den Mittelpunkt M = ( 0 | 0 ) {\displaystyle M=(0|0)} , wenn x κ ( t 0 ) 2 + y κ ( t 0 ) 2 = r 2 {\displaystyle x_{\kappa }(t_{0})^{2}+y_{\kappa }(t_{0})^{2}=r^{2}} („Koordinatengleichung“) gilt. Definiert man

f κ ( t ) = x κ ( t ) 2 + y κ ( t ) 2 r 2 r 2 ; 0 t 1 , {\displaystyle f_{\kappa }(t)={\frac {x_{\kappa }(t)^{2}+y_{\kappa }(t)^{2}-r^{2}}{r^{2}}};\quad 0\leq t\leq 1,}

so ist das äquivalent zu f κ ( t 0 ) = 0 {\displaystyle f_{\kappa }(t_{0})=0} . f κ ( t ) {\displaystyle f_{\kappa }(t)} ist ein Maß für die Abweichung der Approximation ( x κ | y κ ) {\displaystyle (x_{\kappa }|y_{\kappa })} von der Kreisgestalt.

Fordert man dann die Übereinstimmung der Bézierkurve ( x κ | y κ ) {\displaystyle (x_{\kappa }|y_{\kappa })} mit dem Kreis Q {\displaystyle Q} bei der Winkelhalbierenden, erhält man

0 = f κ ( 0 , 5 ) = 9 κ 2 + 24 κ 16 32 κ = 4 3 ( 2 1 ) 0,552 28 {\displaystyle 0=f_{\kappa }(0{,}5)={\frac {9\kappa ^{2}+24\kappa -16}{32}}\quad \Longrightarrow \quad \kappa ={\frac {4}{3}}({\sqrt {2}}-1)\approx 0{,}55228}

Der Fehler ist Null bei t { 0 ; 0 , 5 ; 1 } {\displaystyle t\in \left\lbrace 0;0{,}5;1\right\rbrace } , sonst überall positiv, d. h. die Bézierkurve liegt stets auf oder außerhalb des Kreisbogens. Der maximale Fehler beträgt f κ ( t ) = 5 , 45 10 4 {\displaystyle f_{\kappa }(t)=5{,}45\cdot 10^{-4}} bei t = 0,211 {\displaystyle t=0{,}211} und bei t = 0,789 {\displaystyle t=0{,}789} .

Fordert man, dass die aufsummierten Fehler über die gesamte Kurve verschwinden ( f κ ( t ) {\displaystyle f_{\kappa }(t)} kann sowohl positiv als auch negativ sein – die Bézierkurve verläuft teils außerhalb, teils innerhalb der Kreislinie – und das Integral darüber kann Null ergeben), erhält man

0 = 0 1 f κ ( t ) d t = 6 κ 2 + 13 κ 9 35 κ = 385 13 12 0,551 78 {\displaystyle 0=\int _{0}^{1}f_{\kappa }(t)\,\mathrm {d} t={\frac {6\kappa ^{2}+13\kappa -9}{35}}\quad \Longrightarrow \quad \kappa ={\frac {{\sqrt {385}}-13}{12}}\approx 0{,}55178}

Die größten Abweichungen liegen bei etwa f κ ( 0 , 5 ) = 5 , 3 10 4 {\displaystyle f_{\kappa }(0{,}5)=-5{,}3\cdot 10^{-4}} und bei f κ ( 0,173 ) = f κ ( 0,827 ) = 3 , 5 10 4 {\displaystyle f_{\kappa }(0{,}173)=f_{\kappa }(0{,}827)=3{,}5\cdot 10^{-4}} . Beide Approximationen sind somit für viele Anwendungsbereiche ausreichend.

Beispielkoordinaten Halbkreis:

Bei einer Halbkreisnäherung mit M = ( 0 | 0 ) {\displaystyle M=(0|0)} , P 0 = ( 0 | r ) {\displaystyle P_{0}=(0|-r)} , P 3 = ( 0 | r ) {\displaystyle P_{3}=(0|r)} , und P 1 = ( κ r | r ) {\displaystyle P_{1}=(\kappa r|-r)} , P 2 = ( κ r | r ) {\displaystyle P_{2}=(\kappa r|r)} mit κ = 1,315 6 {\displaystyle \kappa =1{,}3156} beträgt die maximale Abweichung f κ ( t ) = 2 , 65 % {\displaystyle f_{\kappa }(t)=2{,}65\,\%} . Dies ist bzgl. der maximalen Abweichung etwa 50 mal schlechter als die Viertelkreisapproximation.

Literatur

  • Gerald Farin: Curves and Surfaces for CAGD. A practical guide. 5. Auflage. Academic Press, San Diego 2002, ISBN 1-55860-737-4
  • J. Hoschek, D. Lasser: Grundlagen der geometrischen Datenverarbeitung. Vieweg+Teubner Verlag, 1989, ISBN 978-3-519-02962-5
  • David Salomon: Curves and Surfaces for Computer Graphics. Springer Science+Business Media, 2006, ISBN 0-387-24196-5
  • Boaswan Dzung Wong: Bézierkurven: gezeichnet und gerechnet. Orell Füssli Verlag, Zürich 2003, ISBN 3-280-04021-3
  • Wolfgang Boehm, Gerald Farin, Jürgen Kahmann: A survey of curve and surface methods in CAGD. In: Comput. Aided Geom. Des., 1, S. 1–60, 1984
  • Eric W. Weisstein: Bézier Curve. In: MathWorld (englisch).
  • Computerunterstützte Darstellende und Konstruktive Geometrie. (CDKG) Uni Darmstadt (PDF; 3,4 MB), S. 74,227.
  • Dr. Hans Walser, Mathematisches Institut der Universität Basel: Die Modellierung des schönen Scheins, Mathematikinformation Nr. 55, Seite 10, Abschnitt 10.2


Vier Stützpunkte

  • Java-Applet

Mehr als vier Stützpunkte

  • TinySpline. Open Source C-Programmbibliothek für NURBS, B-Splines und Bézier Splines mit Bindings für verschiedene Sprachen

Einzelnachweise

  1. Farin, S. 37
  2. Hoschek&Lasser, S. 120
  3. Farin S. 38
  4. Farin S. 40
  5. Farin: S. 48
  6. Farin: S. 44
  7. Farin: S. 51
  8. Farin S. 87
  9. Hoschek&Lasser S. 128
  10. Farin S. 231
  11. Hoschek&lasser S. 143