Non-Uniform Rational B-Spline
# Non-Uniform Rational B-Spline

Non-uniform rational B-Splines (deutsch: nicht-uniforme rationale B-Splines, kurz NURBS) sind mathematisch definierte Kurven oder Flächen, die im Computergrafik-Bereich, beispielsweise im CGI oder CAD, zur Modellierung beliebiger Formen verwendet werden. Ein NURBS kann jeden beliebigen nicht-verzweigenden stetigen Linienzug darstellen. Im Computergrafik-Bereich wird die Geometrieinformation jedoch meist über stückweise funktional definierte Geometrie-Elemente dargestellt, statt die komplette Formgebung über einen einzigen NURBS abzubilden.

## Entwicklung

In den 1950er Jahren wurden besonders im Automobil- und Schiffbau für die fehlerfreie Reproduzierbarkeit technischer Bauteile mathematisch exakte Beschreibungen von Freiformflächen benötigt. Vor dieser Zeit wurden derartige Oberflächen durch einzelne von einem Konstrukteur hergestellte physikalische Modelle beschrieben.

So begann unter anderem der Ingenieur Pierre Étienne Bézier, zu dieser Zeit bei Renault in Frankreich, mit der Entwicklung der nach ihm benannten Bézierkurve. Unabhängig von Bézier arbeitete Paul de Casteljau, angestellt bei Citroën, zur gleichen Zeit auch an diesem mathematischen Problem. Weil Bézier die Ergebnisse seiner Arbeit veröffentlichte, werden heutzutage in der graphischen Datenverarbeitung Splines, deren Kontrollpunkte nicht auf der Kurve selbst liegen, als „Bézier-Spline“ bezeichnet, während der Name von de Casteljau in dem nach ihm benannten Algorithmus fortlebt, der für die numerische Verarbeitung parametrischer Flächen eingesetzt wird. In den 1960er Jahren wurde klar, dass non-uniform rational B-Splines (NURBS) eine Generalisierung von Bézier-Splines sind, die als uniform non-rational B-Splines angesehen werden können.

Zunächst wurden NURBS nur in proprietären CAD-Werkzeugen von Automobilunternehmen verwendet. Später hielten sie Einzug in weiter verbreitete Computergrafik-Anwendungen, beispielsweise die Open Graphics Library (OpenGL).

Auf Systemen von Silicon Graphics wurde erstmals 1989 die Echtzeitverarbeitung und interaktives Rendering von NURBS-Kurven und -Flächen möglich. Auf der Cebit 1994 stellten die Firmen CAS Berlin und CAA in Zusammenarbeit mit der TU Berlin den ersten, auf den Namen NöRBS getauften, interaktiven NURBS-Modellierer für PCs vor. Heutzutage enthalten die meisten professionellen Computergrafik-Anwendungen NURBS-Technologie, die in den meisten Fällen durch die Integration einer NURBS-Engine realisiert ist, die von einer darauf spezialisierten Firma bereitgestellt wird.

## Anwendung

NURBS sind bei der computergestützten Konstruktion (CAD) und Fertigung (CAM) beinahe unersetzlich und Teil zahlreicher Industriestandards, wie IGES (Initial Graphics Exchange Specification), STEP (STandard for the Exchange of Product model data) und PHIGS (Programmer’s Hierarchical Interactive Graphics System). Im Allgemeinen ist die interaktive Bearbeitung von NURBS-Kurven und -Flächen sehr intuitiv und vorhersagbar. Kontrollpunkte sind stets entweder direkt mit der Kurve oder Fläche verbunden oder wirken, als seien sie mit einem Gummiband verbunden. Eine Manipulation der Geometrieelemente kann – am offensichtlichsten bei Bézierkurven – direkt an den Kontrollpunkten durchgeführt oder durch übergeordnete Werkzeuge realisiert werden. Derartige Werkzeuge basieren zum Teil auf der Eigenschaft von NURBS, Kurven und Flächen unterschiedlicher Stetigkeit darstellen zu können.

Hierbei werden unterschiedlich starke Stetigkeitsforderungen gestellt. Neben der üblichen ${\mathcal {C}}$-Stetigkeit wird hier auch die geometrische Stetigkeit betrachtet:

* ${\mathcal {G}}^{0}:$ Positions-Stetigkeit gilt, wenn die Endpunkte zweier Kurven oder Flächen zusammentreffen. Trotzdem können Kurven oder Flächen sich in einem Winkel berühren, der zu einer scharfen Kante oder Ecke an dieser Stelle führt und Störungen bei Lichteffekten verursacht.
* ${\mathcal {G}}^{1}:$ Bei tangentialer Stetigkeit verhindert die Parallelität der Endvektoren von Kurven oder Flächen das Auftreten scharfer Kanten. Tangentiale Stetigkeit ist oftmals ausreichend, da die Beleuchtung derartiger Geometrien kontinuierlich und damit natürlich erscheint.
* ${\mathcal {G}}^{2}:$ Krümmungsstetigkeit herrscht bei gleichen Krümmungswerten im gemeinsamen Punkt von aneinander grenzender Kurven, respektive gemeinsamen Linien aneinander grenzender Flächen. Diese Stetigkeit ist beispielsweise notwendig für überströmte Flächen, um ein Ablösen der Strömung ausschließen zu können.

Höhere Stetigkeitsgrade sind mit NURBS ebenfalls möglich.

Die Realisierung glatter Oberflächen erfordert NURBS-Flächen, die wenigstens ${\mathcal {G}}^{1}$-Stetigkeit erreichen. Fehler in der Oberfläche können durch Beleuchtung und Reflexion detektiert werden. Eine Methode zur Bewertung einer Fläche basiert auf der Auswertung eines durch Ray Tracing oder Reflexion und Abbildung aufgenommenen Bilds einer Fläche mit weißen auf sie reflektierten Streifen, womit selbst die kleinsten Abweichungen der Fläche aufgedeckt werden können. Diese Methode stammt aus dem Prototypenbau im Automobilbereich, wo die Oberflächenqualität durch Überprüfung der Reflexionen von einem Neonlichthimmel auf der Karosserie sichergestellt wird.

## Mathematische Beschreibung

NURBS-Kurven und -Flächen haben eine Reihe interessanter Eigenschaften:

* Sie sind invariant für projektive Transformationen.
* Sie bieten eine gemeinsame mathematische Darstellung für sowohl analytische Standardformen (z. B. Kegelschnitte) als auch Freiformflächen.
* Sie reduzieren den Speicheraufwand für geometrische Objekte (im Vergleich zu einfacheren Methoden).
* Sie können durch numerisch stabile und präzise Algorithmen verhältnismäßig schnell ausgewertet werden.
* Sie sind Verallgemeinerungen von nicht-rationalen B-Splines und nicht-rationalen und rationalen Bézier-Kurven und -Flächen. Die Aussage, dass NURBS-Kurven eine Verallgemeinerung von Bézierkurven sind, bedeutet, dass alle Bézierkurven NURBS-Kurven, aber nicht alle NURBS-Kurven Bézierkurven sind.

Eine NURBS-Kurve $\textstyle C(u)$ ist definiert durch den Grad $\textstyle k$ ihrer Basispolynome (= Ordnung p der NURBS-Kurve – 1), eine Menge $\textstyle \mathbb {P}$ gewichteter ($\textstyle w_{i}$) Kontrollpunkte $\textstyle P_{i}$ und einen Knotenvektor $\textstyle U$. NURBS-Kurven und -Flächen sind Verallgemeinerungen von sowohl B-Splines als auch Bézierkurven und -flächen. Der hauptsächliche Unterschied zu diesen beiden Splinearten ist die Gewichtung der Kontrollpunkte mit den Gewichten $\textstyle w_{i}$. Durch die $\textstyle w_{i}$ werden NURBS-Kurven rational.

Eine NURBS-Kurve ist über die Summe der mit rationalen B-Spline-Basisfunktionen $\textstyle R_{i,k}$ gewichteten Kontrollpunkte $\textstyle P_{i}$ , also über die Formel

    $C(u)=\sum _{i=0}^{n}R_{i,k}(u)P_{i}$

vollständig definiert. Die rationale B-Spline-Basisfunktion errechnet sich aus B-Spline-Basisfunktionen $\textstyle N_{i,k}$ des Grads der Basispolynome $\textstyle k$ und den zu den Kontrollpunkten zugehörigen Gewichten $\textstyle w_{i}$ zu

    $R_{i,k}(u)={\frac {N_{i,k}(u)w_{i}}{\sum _{j=0}^{n}N_{j,k}(u)w_{j}}}$.

Der Parameter $\textstyle u\in [a,b]$ schaltet im Bereich des Knotenvektors

    $\displaystyle U=\{\underbrace {a_{(0)},\ldots ,a_{(p-1)}} _{p},u_{p},\ldots ,u_{n},\underbrace {b_{(n+1)},\ldots ,b_{(n+p)}} _{p}\}$

die einzelnen Segmente der Spline-Kurve aktiv. Die Elemente des Knotenvektors sind monoton steigend, wobei alle $\textstyle a_{(i)}=a_{(j)}$ sowie alle $\textstyle b_{(i)}=b_{(j)}$ sind.

### Randbedingungen und Definitionen

Kontrollpunkte

    $\mathbb {P} =\{P_{i}\mid P_{i}=(x,y){\text{ mit }}x,y\in \mathbb {R} {\text{ und }}i=0..n\in \mathbb {N} \}$
    $|\mathbb {P} |=m=(n+1)$
    $\mathbb {W} =\{w_{i}\in \mathbb {R} ,i=0..n\in \mathbb {N} \}$
    $|\mathbb {W} |=m=(n+1)$

Knotenvektor

Der Knotenvektor $\textstyle U$ besteht aus Parameterwerten $\textstyle u_{i}$, die den Einfluss der Kontrollpunkte $\textstyle P_{i}$ auf die NURBS-Kurve festlegen. Die Anzahl der Knoten ist immer gleich der Anzahl der Kontrollpunkte plus dem Grad der Kurve plus eins ($\textstyle |U|=r=m+k+1=n+p+1$). Beispielsweise hat eine Kurve dritten Grades mit vier Kontrollpunkten acht Knoten (4 + 3 + 1 = 8). Für den Graphik-Designer sind Knoten gewöhnlich nicht hilfreich; sie werden nur für interne Berechnungen benötigt. Aus diesem Grund sind Knoten in vielen Graphik-Design-Programmen nicht änderbar oder überhaupt sichtbar. Neuere Versionen von NURBS-Software (z. B. Alias Maya) erlauben interaktive Änderungen an Knotenpositionen, was allerdings bedeutend weniger intuitiv ist als Änderungen an Kontrollpunkten.

Die Werte des Knotenvektors müssen in aufsteigender Reihenfolge vorliegen. Damit ist (0, 0, 1, 2, 3) gültig, (0, 0, 2, 1, 3) dagegen nicht. Die einzelnen Knotenwerte haben keine Aussage für sich selbst; einzig und allein die Verhältnisse der Differenzen zwischen den Knotenwerten haben eine Bedeutung. Demzufolge ergeben die Knotenvektoren (0, 0, 1, 2, 3), (0, 0, 2, 4, 6), und (1, 1, 2, 3, 4) alle die gleiche Kurve. Weiterhin muss die Vielfachheit eines Knotens kleiner-gleich sein als der Grad der Kurve (kein Knoten darf öfter als der Grad der Kurve auftauchen). Für NURBS ersten Grades ist jeder Knoten gepaart mit einem Kontrollpunkt.

Ordnung einer NURBS-Kurve

Eine NURBS-Kurve ist von der Ordnung $\textstyle p=k+1$, sie besteht aus Polynomen vom Grad $\textstyle k$. Die Ordnung $\textstyle p$ einer NURBS-Kurve ist festgelegt durch die Anzahl benachbarter Kontrollpunkte, die die Kurve beeinflussen. Die Kurve setzt sich mathematisch aus Polynomen zusammen, deren Grad $\textstyle k$ eins kleiner als die Ordnung der Kurve ist ($\textstyle k=p-1$). Also werden Kurven zweiter Ordnung (die durch lineare Polynome dargestellt werden) lineare Kurven, Kurven dritter Ordnung quadratische Kurven und Kurven vierter Ordnung kubische Kurven genannt. Die Anzahl $\textstyle m$ der Kontrollpunkte muss größer oder gleich der Ordnung $\textstyle p$ der Kurve sein.

In der Praxis werden am häufigsten kubische Kurven gebraucht. Kurven fünften oder sechsten Grades sind manchmal nützlich, gerade für Ableitungen, Kurven höheren Grades werden in der Praxis aber nie benutzt, da sie zu internen numerischen Problemen führen und ihre Berechnung tendenziell unverhältnismäßig viel Berechnungszeit erfordert.

### Einheitskreis

Im Gegensatz zu nicht-rationalen Kurven, die nicht geeignet sind, Kreise darzustellen, klappt dies mit NURBS problemlos. Der Einheitskreis in der xy-Ebene kann beispielsweise als NURBS-Kurve vom Grad 2 mit Knotenvektor (0,0,0,1,1,2,2,3,3,4,4,4) und folgenden Kontrollpunkten und Gewichten konstruiert werden.

| x | y | z | Gewicht |
| - | - | - | - |
| 1 | 0 | 0 | 1 |
| 1 | 1 | 0 | √2 / 2 |
| 0 | 1 | 0 | 1 |
| −1 | 1 | 0 | √2 / 2 |
| −1 | 0 | 0 | 1 |
| −1 | −1 | 0 | √2 / 2 |
| 0 | −1 | 0 | 1 |
| 1 | −1 | 0 | √2 / 2 |
| 1 | 0 | 0 | 1 |

### NURBS-Flächen

Während eine NURBS-Kurve ausschließlich in eine parametrische Richtung $u$ aufgespannt ist, wird eine NURBS-Fläche durch zwei Parameter, genannt $u$ und $v$, aufgespannt. Die Kurve kann durch Auswertung an unterschiedlichen Parametern im kartesischen zwei- oder dreidimensionalen Raum abgebildet werden. Analog erfolgt die Abbildung einer NURBS-Fläche im kartesischen Raum durch Auswertung mit verschiedenen Werten für zwei Parameter.

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

sind definiert durch ein Kontrollgitter $P_{i,j}$ und die rationale Basisfunktion

    $R_{i,j}(u,v)={\frac {N_{i,p}(u)N_{j,q}(v)w_{i,j}}{\sum _{k=0}^{n}\sum _{l=0}^{m}N_{k,p}(u)N_{l,q}(v)w_{k,l}}}$

mit einer zweidimensionalen Gewichtematrix $w_{i,j}$. Die zweite Dimension der Fläche $S(u,v)$ mit dem Parameter $v$ wird geschaltet durch den analog zu $U$ aufgebauten Knotenvektor

    $V=\{\underbrace {c,\ldots ,c} _{q},v_{q+1},\ldots ,v_{s-q-1},\underbrace {d,\ldots ,d} _{q}\}$

mit der Ordnung $q$ und der Länge $s$.

## Einfache Erklärung eines komplexen Sachverhaltes

In den 1950er-Jahren hat Bézier eine mathematische Formel entdeckt, mit der sich eine geschwungene Kurve beschreiben lässt. Die Kurve wird dabei durch wenige sogenannte Kontrollpunkte definiert. In der Folge hat sich herausgestellt, dass diese Bézierkurven noch nicht ausreichend sind, um beispielsweise einen Kreis exakt darzustellen. Daraufhin wurden mathematische Formeln entwickelt, bei denen die einzelnen Kontrollpunkte unterschiedlich stark (rational) auf den Kurvenverlauf einwirken können, die sogenannten B-Splines. Der Begriff Spline ist historisch und stammt aus dem Schiffbau. Linealartige biegsame Metallstreifen, sogenannte Straklatten (engl. splines), wurden seitlich mit Gewichten beschwert und konnten somit eine definierte Kurve darstellen. Mehrere dieser Kurven hintereinander stellten die Spanten dar, die wiederum den Rumpf des Schiffes ergaben.

Später wurden dann mathematische Formeln entwickelt, die sowohl Béziers Formel als auch die B-Splines verallgemeinerten, die sogenannten „Non uniform rational B-Splines“. Damit kann man also beliebige Kurven darstellen und wenn man nicht nur eine Richtung („Länge“) wie bei einer Kurve, sondern eine zweite („Breite“) hinzufügt, können damit beliebige Flächen dargestellt werden. Theoretisch können diese Flächen beliebig komplex und groß sein. Da der Aufwand zur Berechnung der Formeln allerdings mit zunehmender Komplexität stark ansteigt und auch mit einem leistungsfähigen Computer nicht mehr in angemessener Zeit berechnet werden kann, werden kompliziertere Flächen durch mehrere aneinandergesetzte sogenannte Patches dargestellt. In der industriellen Produktion werden mittlerweile fast ausschließlich NURBS eingesetzt, um jede Art von herzustellendem Objekt (Joghurtbecher, Motorgehäuse, Sonnenbrillen usw.) mathematisch exakt abbilden zu können.

