Bézierkurve

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.

Kubische Bézierkurve

Definition

Konstruktion einer kubischen Bézierkurve

Eine Bézierkurve $n$$n$-ten Grades zu gegebenen $n+1$$n+1$ Kontroll- oder Bézierpunkten $P_{0},P_{1},\dotsc ,P_{n}$$P_{0},P_{1},\dotsc ,P_{n}$, die das so genannte Kontrollpolygon bilden, ist für $t\in [0,1]$$t\in [0,1]$ definiert als

$C(t)=\sum _{i=0}^{n}B_{i,n}(t)P_{i}\,$$C(t)=\sum _{{i=0}}^{n}B_{{i,n}}(t)P_{i}\,$

wobei

$B_{i,n}(t)={\binom {n}{i}}t^{i}(1-t)^{n-i}$$B_{{i,n}}(t)={\binom ni}t^{i}(1-t)^{{n-i}}$

das i-te Bernsteinpolynom n-ten Grades ist. Diese bilden eine Basis des Vektorraums der Polynome und erfüllen die Rekursionsformel

$B_{i,n}(t)=(1-t)\cdot B_{i,n-1}(t)\,+\,t\cdot B_{i-1,n-1}(t)$$B_{{i,n}}(t)=(1-t)\cdot B_{{i,n-1}}(t)\,+\,t\cdot B_{{i-1,n-1}}(t)$

mit $B_{i,n}:=0$$B_{{i,n}}:=0$ für $i<0$$i<0$ oder $i>n$$i>n$, sowie $B_{0,0}:=1$$B_{{0,0}}:=1$. Dies erlaubt eine numerisch stabile rekursive Berechnung der Werte einer Bézierkurve mit Hilfe des De-Casteljau-Algorithmus:

${\begin{aligned}C_{i}^{0}(t)&\;:=\;P_{i}\\C_{i}^{j}(t)&\;:=\;(1-t)C_{i}^{j-1}(t)\,+\,tC_{i+1}^{j-1}(t)\quad {\text{mit}}\quad j=1,\ldots ,n\quad {\text{und}}\quad i=0,\ldots ,n-j\end{aligned}}$${\begin{aligned}C_{i}^{0}(t)&\;:=\;P_{i}\\C_{i}^{j}(t)&\;:=\;(1-t)C_{i}^{{j-1}}(t)\,+\,tC_{{i+1}}^{{j-1}}(t)\quad {\text{mit}}\quad j=1,\ldots ,n\quad {\text{und}}\quad i=0,\ldots ,n-j\end{aligned}}$

Eigenschaften

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 „oszilliert“, d. h., ihre Richtung und Krümmung verändert.

• Die Kurve liegt innerhalb der konvexen Hülle des Kontrollpolygons. Dies folgt daraus, dass die Bernsteinpolynome vom Grad n eine Zerlegung der Eins sind:

$1=\sum _{i=0}^{n}B_{i,n}(t)\quad t\in [0,1]$$1=\sum _{{i=0}}^{n}B_{{i,n}}(t)\quad t\in [0,1]$

• Die Kurve geht genau durch die Endpunkte $P_{0}$$P_{0}$ und $P_{n}$$P_{n}$:

${\begin{aligned}C(0)\;&=\;\sum _{i=0}^{n}{\binom {n}{i}}0^{i}(1-0)^{n-i}P_{i}\\&=\;P_{0}\\C(1)\;&=\;P_{n}\end{aligned}}$${\begin{aligned}C(0)\;&=\;\sum _{{i=0}}^{n}{\binom ni}0^{i}(1-0)^{{n-i}}P_{i}\\&=\;P_{0}\\C(1)\;&=\;P_{n}\end{aligned}}$

• Die Tangenten in den Endpunkten verlaufen entlang der Kanten $P_{0}P_{1}$$P_{0}P_{1}$ bzw. $P_{n-1}P_{n}$$P_{{n-1}}P_{n}$ des Kontrollpolygons:

${\begin{aligned}C'(0)\;&=\;n\,(P_{1}-P_{0})\\C'(1)\;&=\;n\,(P_{n}-P_{n-1})\end{aligned}}$${\begin{aligned}C'(0)\;&=\;n\,(P_{1}-P_{0})\\C'(1)\;&=\;n\,(P_{n}-P_{{n-1}})\end{aligned}}$

• Die ersten Summanden des Taylorpolynoms bei $t=0$$t=0$ bzw. bei $t=1$$t=1$ lauten für $n\geq 2$$n\geq 2$:

$C(t)=P_{0}+n(P_{1}-P_{0})t+{\frac {n(n-1)}{2}}(P_{0}-2P_{1}+P_{2})t^{2}+{\mathcal {O}}(t^{3})$$C(t)=P_{{0}}+n(P_{{1}}-P_{{0}})t+{\frac {n(n-1)}{2}}(P_{{0}}-2P_{{1}}+P_{{2}})t^{{2}}+{\mathcal {O}}(t^{{3}})$

$C(t)=P_{n}+n(P_{n-1}-P_{n})(1-t)+{\frac {n(n-1)}{2}}(P_{n-2}-2P_{n-1}+P_{n})(1-t)^{2}+{\mathcal {O}}((1-t)^{3})$$C(t)=P_{{n}}+n(P_{{n-1}}-P_{{n}})(1-t)+{\frac {n(n-1)}{2}}(P_{{n-2}}-2P_{{n-1}}+P_{{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 vorherigen Stützpunkten ergeben. Dabei ist der Trennungspunkt vom Parameter $t$$t$ abhängig. Aus der Abbildung für die Konstruktion einer Bézierkurve ist ersichtlich, dass sich die neue erste Kurve aus den Kontrollpunkten $P_{0},Q_{0},R_{0},B$$P_{0},Q_{0},R_{0},B$ zusammensetzt, während die zweite Kurve aus den Punkten $B,R_{1},Q_{2},P_{3}$$B,R_{1},Q_{2},P_{3}$ besteht. Diese Eigenschaft kann genutzt werden, um eine Kurve rekursiv mit Hilfe des De-Casteljau-Algorithmus durch Geraden anzunähern.

Als verallgemeinerte Form der Bézierkurve kann die Bézierfläche gesehen werden. Eine Bézierfläche $(n,m)$$(n,m)$-ter Ordnung ist eine Fläche der Form

$C(u,v)=\sum _{i=0}^{n}\sum _{j=0}^{m}P_{i,j}B_{i,n}(u)B_{j,m}(v)$$C(u,v)=\sum _{{i=0}}^{n}\sum _{{j=0}}^{m}P_{{i,j}}B_{{i,n}}(u)B_{{j,m}}(v)$

mit den Kontrollpunkten $P_{i,j}$$P_{{i,j}}$ und den Bernsteinpolynomen $B_{i,n}(u)$$B_{{i,n}}(u)$ und $B_{j,m}(v)$$B_{{j,m}}(v)$.

Eine Bézierfläche kann also durch zwei zueinander orthogonale Bézierkurven beschrieben werden.

Bézierkurven bis zum dritten Grad

Lineare Bézierkurven (n=1)

Zwei Kontrollpunkte $P_{0}$$P_{0}$ und $P_{1}$$P_{1}$ bestimmen eine gerade Strecke zwischen diesen beiden Punkten. Der Verlauf dieser linearen Bézier„kurve“ wird angegeben durch

$C(t)=\sum _{i=0}^{1}t^{i}(1-t)^{1-i}P_{i}=\ (1-t)P_{0}+tP_{1}{\mbox{ , }}t\in [0,1]$$C(t)=\sum _{{i=0}}^{1}t^{i}(1-t)^{{1-i}}P_{i}=\ (1-t)P_{0}+tP_{1}{\mbox{ , }}t\in [0,1]$

Quadratische Bézierkurven (n=2)

Eine quadratische Bézierkurve ist der Pfad, der durch die Funktion $C(t)$$C(t)$ für die Punkte $P_{0}$$P_{0}$, $P_{1}$$P_{1}$ und $P_{2}$$P_{2}$ verfolgt wird:

${\begin{aligned}C(t)\ &=\ \sum _{i=0}^{2}{\binom {2}{i}}t^{i}(1-t)^{2-i}P_{i}\\\ &=\ (1-t)^{2}P_{0}+2t(1-t)P_{1}+t^{2}P_{2}\\\ &=\ (P_{0}-2P_{1}+P_{2})t^{2}+(-2P_{0}+2P_{1})t+P_{0}{\mbox{ , }}t\in [0,1]\end{aligned}}$${\begin{aligned}C(t)\ &=\ \sum _{{i=0}}^{2}{\binom 2i}t^{i}(1-t)^{{2-i}}P_{i}\\\ &=\ (1-t)^{{2}}P_{0}+2t(1-t)P_{1}+t^{{2}}P_{2}\\\ &=\ (P_{0}-2P_{1}+P_{2})t^{{2}}+(-2P_{0}+2P_{1})t+P_{0}{\mbox{ , }}t\in [0,1]\end{aligned}}$

Mit Hilfe des De-Casteljau-Algorithmus ausgedrückt: $C(t)=C_{0}^{2}(t)$$C(t)=C_{{0}}^{{2}}(t)$

$C_{0}^{2}(t)=(1-t)\underbrace {{\Bigl [}(1-t)\overbrace {P_{0}} ^{C_{0}^{0}}+t\overbrace {P_{1}} ^{C_{1}^{0}}{\Bigr ]}} _{C_{0}^{1}(t)}+t\underbrace {{\Bigl [}(1-t)\overbrace {P_{1}} ^{C_{1}^{0}}+t\overbrace {P_{2}} ^{C_{2}^{0}}{\Bigr ]}} _{C_{1}^{1}(t)}$$C_{{0}}^{{2}}(t)=(1-t)\underbrace {{\Bigl [}(1-t)\overbrace {P_{{0}}}^{{C_{{0}}^{{0}}}}+t\overbrace {P_{{1}}}^{{C_{{1}}^{{0}}}}{\Bigr ]}}_{{C_{{0}}^{{1}}(t)}}+t\underbrace {{\Bigl [}(1-t)\overbrace {P_{{1}}}^{{C_{{1}}^{{0}}}}+t\overbrace {P_{{2}}}^{{C_{{2}}^{{0}}}}{\Bigr ]}}_{{C_{{1}}^{{1}}(t)}}$

Die Strecken $C_{0}^{1}(t)$$C_{{0}}^{{1}}(t)$ und $C_{1}^{1}(t)$$C_{{1}}^{{1}}(t)$ sind die Kanten des Kontrollpolygons (graue Linien in nebenstehender Animation). Die Kurvenschar ${\tilde {C}}_{0}^{2}({\tilde {t}},t)=(1-{\tilde {t}})C_{0}^{1}(t)+{\tilde {t}}C_{1}^{1}(t)$${\tilde {C}}_{{0}}^{{2}}({\tilde {t}},t)=(1-{\tilde {t}})C_{{0}}^{{1}}(t)+{\tilde {t}}C_{{1}}^{{1}}(t)$ mit ${\tilde {t}}\in [0,1]$${\tilde {t}}\in [0,1]$ entspricht den grünen Linien in der Animation. Die Auswertung an den Stellen mit ${\tilde {t}}=t$${\tilde {t}}=t$ gibt den Verlauf der Bézierkurve an: $C_{0}^{2}(t)={\tilde {C}}_{0}^{2}(t,t)$$C_{{0}}^{{2}}(t)={\tilde {C}}_{{0}}^{{2}}(t,t)$.

Kubische Bézierkurven (n=3)

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.

Vier Punkte ($P_{0}$$P_{0}$, $P_{1}$$P_{1}$, $P_{2}$$P_{2}$ und $P_{3}$$P_{3}$) bestimmen eine kubische Bézierkurve. Die Kurve beginnt bei $P_{0}$$P_{0}$ und geht in Richtung $P_{1}$$P_{1}$ und dann aus Richtung $P_{2}$$P_{2}$ zu $P_{3}$$P_{3}$. Im Allgemeinen geht die Kurve nicht durch $P_{1}$$P_{1}$ und $P_{2}$$P_{2}$ – diese Punkte dienen nur der Richtungsangabe, wobei $P_{1}$$P_{1}$ die Richtung bestimmt, in welche die Kurve in $P_{0}$$P_{0}$ geht. $P_{2}$$P_{2}$ legt die Richtung fest, aus welcher die Kurve zu $P_{3}$$P_{3}$ geht. Der Abstand zwischen $P_{0}$$P_{0}$ und $P_{1}$$P_{1}$ und der Abstand von $P_{2}$$P_{2}$ und $P_{3}$$P_{3}$ bestimmen, „wie weit“ sich die Kurve in Richtung der Kontrollpunkte $P_{1}$$P_{1}$ und $P_{2}$$P_{2}$ bewegt, bevor sie in Richtung $P_{3}$$P_{3}$ läuft.

${\begin{aligned}C(t)\ &=\ \sum _{i=0}^{3}{\binom {3}{i}}t^{i}(1-t)^{3-i}P_{i}\\&=\ (1-t)^{3}P_{0}+3t(1-t)^{2}P_{1}+3t^{2}(1-t)P_{2}+t^{3}P_{3}\\&=\ (-P_{0}+3P_{1}-3P_{2}+P_{3})t^{3}+(3P_{0}-6P_{1}+3P_{2})t^{2}+(-3P_{0}+3P_{1})t+P_{0},t\in [0,1]\end{aligned}}$${\begin{aligned}C(t)\ &=\ \sum _{{i=0}}^{3}{\binom 3i}t^{i}(1-t)^{{3-i}}P_{i}\\&=\ (1-t)^{3}P_{0}+3t(1-t)^{2}P_{1}+3t^{2}(1-t)P_{2}+t^{3}P_{3}\\&=\ (-P_{0}+3P_{1}-3P_{2}+P_{3})t^{3}+(3P_{0}-6P_{1}+3P_{2})t^{2}+(-3P_{0}+3P_{1})t+P_{0},t\in [0,1]\end{aligned}}$

Mit Hilfe des De-Casteljau-Algorithmus ausgedrückt: $C(t)=C_{0}^{3}(t)$$C(t)=C_{{0}}^{{3}}(t)$

${\begin{aligned}C_{0}^{3}(t)=(1-t)&\underbrace {{\Bigl \{}(1-t)\overbrace {{\bigl [}(1-t)P_{0}+tP_{1}{\bigr ]}} ^{C_{0}^{1}(t)}+t\overbrace {{\bigl [}(1-t)P_{1}+tP_{2}{\bigr ]}} ^{C_{1}^{1}(t)}{\Bigr \}}} _{C_{0}^{2}(t)}+\\+t&\underbrace {{\Bigl \{}(1-t)\overbrace {{\bigl [}(1-t)P_{1}+tP_{2}{\bigr ]}} ^{C_{1}^{1}(t)}+t\overbrace {{\bigl [}(1-t)P_{2}+tP_{3}{\bigr ]}} ^{C_{2}^{1}(t)}{\Bigr \}}} _{C_{1}^{2}(t)}\end{aligned}}$${\begin{aligned}C_{{0}}^{{3}}(t)=(1-t)&\underbrace {{\Bigl \{}(1-t)\overbrace {{\bigl [}(1-t)P_{{0}}+tP_{{1}}{\bigr ]}}^{{C_{{0}}^{{1}}(t)}}+t\overbrace {{\bigl [}(1-t)P_{{1}}+tP_{{2}}{\bigr ]}}^{{C_{{1}}^{{1}}(t)}}{\Bigr \}}}_{{C_{{0}}^{{2}}(t)}}+\\+t&\underbrace {{\Bigl \{}(1-t)\overbrace {{\bigl [}(1-t)P_{{1}}+tP_{{2}}{\bigr ]}}^{{C_{{1}}^{{1}}(t)}}+t\overbrace {{\bigl [}(1-t)P_{{2}}+tP_{{3}}{\bigr ]}}^{{C_{{2}}^{{1}}(t)}}{\Bigr \}}}_{{C_{{1}}^{{2}}(t)}}\end{aligned}}$

Die Strecken $C_{0}^{1}(t)$$C_{{0}}^{{1}}(t)$ und $C_{1}^{1}(t)$$C_{{1}}^{{1}}(t)$ sowie $C_{2}^{1}(t)$$C_{{2}}^{{1}}(t)$ sind die Kanten des Kontrollpolygons (graue Linien in nebenstehender Animation). Die beiden Kurvenscharen ${\tilde {C}}_{0}^{2}({\tilde {t}},t)=(1-{\tilde {t}})C_{0}^{1}(t)+{\tilde {t}}C_{1}^{1}(t)$${\tilde {C}}_{{0}}^{{2}}({\tilde {t}},t)=(1-{\tilde {t}})C_{{0}}^{{1}}(t)+{\tilde {t}}C_{{1}}^{{1}}(t)$ und ${\tilde {C}}_{1}^{2}({\tilde {t}},t)=(1-{\tilde {t}})C_{1}^{1}(t)+{\tilde {t}}C_{2}^{1}(t)$${\tilde {C}}_{{1}}^{{2}}({\tilde {t}},t)=(1-{\tilde {t}})C_{{1}}^{{1}}(t)+{\tilde {t}}C_{{2}}^{{1}}(t)$ mit ${\tilde {t}}\in [0,1]$${\tilde {t}}\in [0,1]$ entsprechen den grünen Linien in der Animation.

Die Kurvenschar ${\tilde {C}}_{0}^{3}({\tilde {t}},t)=(1-{\tilde {t}})C_{0}^{2}(t)+{\tilde {t}}C_{1}^{2}(t)$${\tilde {C}}_{{0}}^{{3}}({\tilde {t}},t)=(1-{\tilde {t}})C_{{0}}^{{2}}(t)+{\tilde {t}}C_{{1}}^{{2}}(t)$ mit ${\tilde {t}}\in [0,1]$${\tilde {t}}\in [0,1]$ entspricht den blauen Linien in der Animation. Die Auswertung an den Stellen mit ${\tilde {t}}=t$${\tilde {t}}=t$ gibt den Verlauf der Bézierkurve an: $C_{0}^{3}(t)={\tilde {C}}_{0}^{3}(t,t)$$C_{{0}}^{{3}}(t)={\tilde {C}}_{{0}}^{{3}}(t,t)$.

Kubische Darstellung quadratischer Bézierkurven

Wählt man die mittleren Bézierpunkte $P_{1}$$P_{1}$ und $P_{2}$$P_{2}$ einer kubischen Bézierkurve wie folgt

$P_{1}:=P_{0}+{\frac {2}{3}}(P_{12}-P_{0})\,,\ P_{2}:=P_{3}+{\frac {2}{3}}(P_{12}-P_{3})$$P_{{1}}:=P_{{0}}+{\frac {2}{3}}(P_{{12}}-P_{{0}})\,,\ P_{{2}}:=P_{{3}}+{\frac {2}{3}}(P_{{12}}-P_{{3}})$

erhält man eine Kurve, die exakt der quadratischen Bézierkurve mit den Punkten $P_{0}$$P_{0}$, $P_{12}$$P_{{12}}$ und $P_{3}$$P_{3}$ entspricht:

${\begin{aligned}C(t)&=(1-t)^{3}P_{0}+3t(1-t)^{2}\left({\frac {2}{3}}P_{12}+{\frac {1}{3}}P_{0}\right)+3t^{2}(1-t)\left({\frac {2}{3}}P_{12}+{\frac {1}{3}}P_{3}\right)+t^{3}P_{3}\\&=(1-t)^{2}P_{0}+2t(1-t)P_{12}+t^{2}P_{3}\end{aligned}}$${\begin{aligned}C(t)&=(1-t)^{{3}}P_{{0}}+3t(1-t)^{{2}}\left({\frac {2}{3}}P_{{12}}+{\frac {1}{3}}P_{{0}}\right)+3t^{{2}}(1-t)\left({\frac {2}{3}}P_{{12}}+{\frac {1}{3}}P_{{3}}\right)+t^{{3}}P_{{3}}\\&=(1-t)^{{2}}P_{{0}}+2t(1-t)P_{{12}}+t^{{2}}P_{{3}}\end{aligned}}$

Damit lassen sich auch dann quadratische Bézierkurven darstellen, falls ein Vektorzeichenprogramm (z. B. Inkscape) bzw. eine Grafikbibliothek (z. B. Cairo) nur kubische unterstützt.

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$$n$ keine Bézierkurve $t\mapsto (x(t)|y(t))$$t\mapsto (x(t)|y(t))$ $n$$n$-ten Grades in einem noch so kleinen Kreisbogen mit Radius $r$$r$ zum Mittelpunkt $(z_{1}|z_{2})$$(z_{1}|z_{2})$, denn $(x(t_{0})|y(t_{0}))$$(x(t_{0})|y(t_{0}))$ liegt genau dann auf dem Kreisbogen, wenn $t_{0}$$t_{0}$ Nullstelle der Polynomfunktion $t\mapsto [x(t)-z_{1}]^{2}+[y(t)-z_{2}]^{2}-r^{2}$$t\mapsto [x(t)-z_{1}]^{2}+[y(t)-z_{2}]^{2}-r^{2}$ vom Grad $2n$$2n$ ist, was höchstens $2n$$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$$Q$ mit folgenden Parametern:

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

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

Beispielkoordinaten Viertelkreis

Als einfaches Beispiel einer Viertelkreisapproximation wählen wir:

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

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

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

$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$$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 $\kappa =0{,}552$$\kappa =0{,}552$, wie die nachfolgende Betrachtung zeigt.

Fehleranalyse

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

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

$f_{\kappa }(t)={\frac {x_{\kappa }(t)^{2}+y_{\kappa }(t)^{2}-r^{2}}{r^{2}}};\quad 0\leq t\leq 1,$$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_{\kappa }(t_{0})=0$$f_{\kappa }(t_{0})=0$. $f_{\kappa }(t)$$f_{\kappa }(t)$ ist ein Maß für die Abweichung der Approximation $(x_{\kappa }|y_{\kappa })$$(x_{\kappa }|y_{\kappa })$ von der Kreisgestalt.

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

$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$$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\in \left\lbrace 0;0{,}5;1\right\rbrace$$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_{\kappa }(t)=5{,}45\cdot 10^{-4}$$f_{\kappa }(t)=5{,}45\cdot 10^{{-4}}$ bei $t=0{,}211$$t=0{,}211$ und bei $t=0{,}789$$t=0{,}789$.

Fordert man, dass die aufsummierten Fehler über die gesamte Kurve verschwinden ($f_{\kappa }(t)$$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=\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$$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_{\kappa }(0{,}5)=-5{,}3\cdot 10^{-4}$$f_{\kappa }(0{,}5)=-5{,}3\cdot 10^{{-4}}$ und bei $f_{\kappa }(0{,}173)=f_{\kappa }(0{,}827)=3{,}5\cdot 10^{-4}$$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)$$M=(0|0)$, $P_{0}=(0|-r)$$P_{0}=(0|-r)$, $P_{3}=(0|r)$$P_{3}=(0|r)$, und $P_{1}=(\kappa r|-r)$$P_{1}=(\kappa r|-r)$, $P_{2}=(\kappa r|r)$$P_{2}=(\kappa r|r)$ mit $\kappa =1{,}3156$$\kappa =1{,}3156$ beträgt die maximale Abweichung $f_{\kappa }(t)=2{,}65\%$$f_{\kappa }(t)=2{,}65\%$. Dies ist bzgl. der maximalen Abweichung etwa 50 mal schlechter als die Viertelkreisapproximation.

Gebrochenrationale Bézierkurve

Gebrochenrationale Bezierkurve dritten Grades mit unterschiedlicher Gewichtung der Kontrollpunkte

Gebrochenrationale Bézierkurven können vereinfacht als Bézierkurven betrachtet werden, deren Kontrollpunkte $P_{i}$$P_{i}$ mit Gewichten/Anziehungskräften $w_{i}$$w_{i}$ versehen sind und die damit den Kurvenverlauf in ihre Richtung hin beeinflussen.

Dazu stellt man die Bézierkurve mittels homogener Koordinaten ${\overline {P}}_{i}=w_{i}{\bigl (}{\begin{smallmatrix}P_{i}\\1\end{smallmatrix}}{\bigr )}$$\overline {P}_{i}=w_{i}{\bigl (}{\begin{smallmatrix}P_{i}\\1\end{smallmatrix}}{\bigr )}$ im projektiven Raum dar.

${\overline {C}}(t)\ =\ \sum _{i=0}^{n}B_{i,n}(t){\overline {P}}_{i}\ =\ {\begin{pmatrix}\sum _{i=0}^{n}B_{i,n}(t)w_{i}P_{i}\\\sum _{i=0}^{n}B_{i,n}(t)w_{i}\end{pmatrix}}$$\overline {C}(t)\ =\ \sum _{{i=0}}^{n}B_{{i,n}}(t)\overline {P}_{i}\ =\ {\begin{pmatrix}\sum _{{i=0}}^{n}B_{{i,n}}(t)w_{i}P_{i}\\\sum _{{i=0}}^{n}B_{{i,n}}(t)w_{i}\end{pmatrix}}$

Die Bezeichnung als gebrochenrationale Bézierkurven ergibt sich aus der Rücktransformation (Dehomogenisierung oder Projektion) in den Ursprungsraum, was praktisch durch das Teilen der zusätzlichen Komponente erfolgt.

$C(t)\ =\ {\frac {\sum _{i=0}^{n}B_{i,n}(t)w_{i}P_{i}}{\sum _{i=0}^{n}B_{i,n}(t)w_{i}}}$$C(t)\ =\ {\frac {\sum _{{i=0}}^{n}B_{{i,n}}(t)w_{i}P_{i}}{\sum _{{i=0}}^{n}B_{{i,n}}(t)w_{i}}}$

Genese

Nach der Transformation in den projektiven Raum wird die Kurve nach dem normalen Bildungsgesetz erzeugt und anschließend wieder in den Ursprungsraum zurück transformiert.

Projektiver Raum und homogene Koordinaten

Die Transformation der Kontrollpunkte in den projektiven Raum mittels ihrer Gewichte verändert im Allgemeinen (also wenn nicht alle Gewichte gleich 1 sind) die Lage der Kontrollpunkte zueinander, verzerrt also das Kontrollpolygon. Diese Verzerrung wirkt sich nun so aus, dass die Kurve sich in der zurücktransformierten Darstellung den Punkten mit höherer Gewichtung stärker nähert.

Veranschaulichen kann man sich dies, indem man zunächst eine 2-dimensionale Bézierkurve im 3-dimensionalen Raum (projektiver Raum) in der $z=1$$z=1$ Ebene (Ursprungsraum) zeichnet. Wobei die $z$$z$-Komponente die zusätzliche homogene Komponente darstellen soll.

${\overline {P}}_{i}=w_{i}{\bigl (}{\begin{smallmatrix}{\bigl (}{\begin{smallmatrix}x\\y\end{smallmatrix}}{\bigr )}=P_{i}\\z=1\end{smallmatrix}}{\bigr )}$$\overline {P}_{i}=w_{i}{\bigl (}{\begin{smallmatrix}{\bigl (}{\begin{smallmatrix}x\\y\end{smallmatrix}}{\bigr )}=P_{i}\\z=1\end{smallmatrix}}{\bigr )}$

Demnach sind die Gewichte $w_{i}=1$$w_{i}=1$ zu setzen. Fasst man die homogenen Kontrollpunkte ${\overline {P}}_{i}$$\overline {P}_{i}$ als Ortsvektoren des projektiven Raums auf, so entfernen sich diese mit zunehmenden Gewichten $w_{i}>1$$w_{i}>1$ vom Ursprung. Die Bézierkurve verlässt also die $z=1$$z=1$ Ebene und wird 3-dimensional. Eine Rücktransformation auf die Projektionsebene $z=1$$z=1$ führt zu einer, zu den stärker gewichteten Kontrollpunkten hingezogenen, rationalen Bézierkurve.

Weblinks

• Eric W. Weisstein: Bézier Curve. In: MathWorld (englisch).

Vier Stützpunkte

• Mehrere Java-Applets mit Erläuterungen
• Java-Applet

Mehr als vier Stützpunkte

• TinySpline: Open Source C-Programmbibliothek für NURBS, B-Splines und Bezier Splines mit Bindings für verschiedene Sprachen

Literatur

• Gerald Farin: Curves and Surfaces for CAGD. A practical guide. 5. Aufl. Academic Press, San Diego 2002, ISBN 1-55860-737-4
• David Salomon: Curves and Surfaces for Computer Graphics. Springer Science+Business Media, Inc., 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, Comput. Aided Geom. Des. 1, S. 1-60, 1984