Garland-Heckbert-Algorithmus

Der Garland-Heckbert-Algorithmus ist ein Verfahren der Computergrafik, mit dem der Detailgrad eines als Dreiecksnetz modellierten 3D-Körpers kontrolliert reduziert werden kann, ohne dabei das grundsätzliche Aussehen des Körpers zu verändern. Dies geschieht inkrementell, wobei immer eine Kante des Dreiecksnetzes gelöscht und die beiden adjazenten Knoten zu einem einzelnen Knoten zusammengefasst werden. Es wird immer diejenige Kante gelöscht, deren Entfernen die geringste optische Veränderung des Körpers bewirkt.

Der Algorithmus wurde 1997 in einer Veröffentlichung von Michael Garland und Paul Heckbert vorgestellt[1].

Algorithmus

Kantenkontraktion

In einem Dreiecksnetz wird das Löschen einer Kante und Zusammenfassen der übrig gebliebenen Knoten als Kantenkontraktion oder -kollaps bezeichnet. Der Algorithmus gibt Auskunft darüber, welche Kante gelöscht werden muss, und an welcher Stelle der neu entstandene Knoten platziert wird.

Dazu wird eine Fehlerquadrik verwendet, die jedem Knoten und jeder Kante zugeordnet wird. Diese gibt Auskunft darüber, wie groß der entstehende optische Fehler wäre, wenn die entsprechende Kante kontrahiert werden würde. Gleichzeitig kann mit ihr berechnet werden, wo der neu entstehende Knoten platziert werden muss, damit genau dieser Fehler (und kein größerer) entsteht.

Pseudocode

Berechne zunächst für jeden Knoten b_i die zugehörige initiale Fehlerquadrik Q_bi. Siehe dazu weiter unten.

Danach iteriere solange, bis der gewünschte Detailgrad erreicht ist:

1. Jede Kante erhält als Quadrik Q_Ej die Summe der Quadriken ihrer beiden Endpunkte.
2. Berechne anhand der Quadrik für jede Kante, an welcher Stelle der zusammengezogene Knoten erstellt werden müsste, um minimalen Fehler beim Kollabieren der Kante zu verursachen.
3. Ordne alle Kanten anhand ihrer Fehler. Kollabiere die Kante mit dem kleinsten Fehler.
4. Ordne dem neuen Knoten als Quadrik die Summer der Quadriken der beiden Endpunkte der Kante zu.
5. Update die Quadriken aller Kanten, die durch den Kollaps der Kante verändert wurden.

Detaillierte Beschreibung

Initiale Fehlerquadriken

Den Fehler, der beim Zusammenfassen zweier Knoten v₁ und v₂ zu einem neuen v_neu entsteht, definiert man als: Quadrat der Summe der Abstände, die v_neu zu allen Flächen hat, die durch irgendeins der Dreiecke definiert werden, welche an v₁ oder v₂ grenzen.

Um zunächst den Abstand eines Punktes v_neu zu einer Fläche F zu berechnen, kann folgende Formel verwendet werden:

${\displaystyle dist(v_{\text{neu}},F)=n^{T}\cdot (v_{\text{neu}}-b)}$

Dabei ist b ein beliebiger Punkt in der Fläche (z. B. einer der Eckpunkte des Dreiecks, das in der Fläche liegt), und n^T ist der Normalenvektor der Fläche (zum Zeilenvektor transponiert).

Bildet man nun wie vorgesehen das Quadrat dieses Abstands, erhält man:

${\displaystyle {\begin{aligned}dist(v_{\text{neu}},F)^{2}&=(n^{T}\cdot (v_{\text{neu}}-b))^{2}\\&=(n^{T}\cdot (v_{\text{neu}}-b))\cdot (n^{T}\cdot (v_{\text{neu}}-b))\\&=(n^{T}v_{\text{neu}}-n^{T}b)\cdot (n^{T}v_{\text{neu}}-n^{T}b)\\&=(n^{T}v_{\text{neu}})^{2}-n^{T}v_{\text{neu}}\cdot n^{T}b-n^{T}b\cdot n^{T}v_{\text{neu}}+(n^{T}b)^{2}\\&={v_{\text{neu}}}^{T}n\cdot n^{T}v_{\text{neu}}-b^{T}n\cdot n^{T}v_{\text{neu}}-{v_{\text{neu}}}^{T}n\cdot n^{T}b+b^{T}n\cdot n^{T}b\end{aligned}}}$

Verwendet man homogene Koordinaten, so ist diese Rechnung gleichbedeutend mit:

${\displaystyle dist(v_{\text{neu}},F)^{2}=(v_{{\text{neu}}_{x}}v_{{\text{neu}}_{y}}v_{{\text{neu}}_{z}}1)\cdot {\begin{pmatrix}nn^{T}&-nn^{T}b\\-b^{T}nn^{T}&b^{T}nn^{T}b\end{pmatrix}}\cdot {\begin{pmatrix}v_{{\text{neu}}_{x}}\\v_{{\text{neu}}_{y}}\\v_{{\text{neu}}_{z}}\\1\end{pmatrix}}}$

Die Matrix ${\displaystyle Q_{F}={\begin{pmatrix}nn^{T}&-nn^{T}b\\-b^{T}nn^{T}&b^{T}nn^{T}b\end{pmatrix}}}$ wird als Fehlerquadrik der Fläche F bezeichnet.

Betrachtet man nun einen Knoten b des Dreiecksnetzes, so kann ihm als Fehlerquadrik Q_b die Summe der Quadriken seiner angrenzenden Dreiecksflächen F_i zugewiesen werden.

${\displaystyle Q_{b}=\sum _{i}Q_{F_{i}}}$

Die Quadrik Q_b wird zur Initialisierung des Algorithmus für jeden Knoten des Dreiecksnetzes berechnet.

Iterationsschritt

Berechnung der zu kontrahierenden Kante

Jede Kante e_j des Dreiecksnetzes erhält als Quadrik Q_ej die Summe der Quadriken ihrer Endpunkte b₁ und b₂:

${\displaystyle Q_{e_{j}}=Q_{b_{1}}+Q_{b_{2}}}$

Falls der Kante e_j vorher noch keine Quadrik zugeordnet war (nur im ersten Iterationsschritt), oder falls sich die Quadrik im Vergleich zum letzten Iterationsschritt verändert hat (weil einer der Endpunkte durch Kantenkontraktion verschoben wurde), muss nun berechnet werden:

• Welcher Punkt v_neu beim Kontrahieren von e_j einen möglichst kleinen Fehler verursachen würde
• Wie groß der verursachte Fehler wäre.

Um den optimalen Punkt v_neu zu bestimmen, muss also der Minimalwert der Fehlerfunktion gefunden werden, d. h. alle partiellen Ableitungen von ${\displaystyle ({v_{\text{neu}}}^{T}\quad 1)\cdot Q_{e_{j}}\cdot {\begin{pmatrix}v_{\text{neu}}\\1\end{pmatrix}}}$ müssen gleich 0 sein.

Berechnet man diese Ableitungen, erhält man folgendes zu lösende Gleichungssystem:

${\displaystyle {\begin{pmatrix}&nn^{T}&&-nn^{T}b\\0&0&0&1\end{pmatrix}}\cdot {\begin{pmatrix}v_{{\text{neu}}_{x}}\\v_{{\text{neu}}_{y}}\\v_{{\text{neu}}_{z}}\\1\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\1\end{pmatrix}}}$

Dabei können sich mehrere Fälle ergeben:

• Die Matrix ${\displaystyle {\begin{pmatrix}&nn^{T}&&-nn^{T}b\\0&0&0&1\end{pmatrix}}}$ ist invertierbar: Der Punkt v_neu kann einfach durch das Gleichungssystem berechnet werden.
• Die Matrix ${\displaystyle {\begin{pmatrix}&nn^{T}&&-nn^{T}b\\0&0&0&1\end{pmatrix}}}$ ist nicht invertierbar: Der optimale Wert für v_neu kann auf der Verbindungsstrecke der beiden Kantenendpunkte b₁ und b₂ gefunden werden. Man beschreibt dann v_neu als Linearkombination der beiden Endpunkte: ${\displaystyle {v_{\text{neu}}}^{*}(\lambda )=(1-\lambda )\cdot b_{1}+\lambda \cdot b_{2}}$ und sucht stattdessen das Minimum von ${\displaystyle dist(v_{\text{neu}},F)^{2}=(({v_{\text{neu}}}^{*})^{T}\quad 1)\cdot Q_{e_{j}}\cdot {\begin{pmatrix}{v_{\text{neu}}}^{*}\\1\end{pmatrix}}}$ bezüglich λ. Dies geschieht durch einfache Ableitung der Funktion nach λ und setzen dieser Ableitung = 0.
• Sollte auch die zweite Methode kein Ergebnis liefern, so wird einfach der Mittelpunkt zwischen b₁ und b₂ verwendet, oder b₁ oder b₂ selbst.

Nachdem nun v_neu bekannt ist, wird der entstandene Fehler ${\displaystyle ({v_{\text{neu}}}^{T}\quad 1)\cdot Q_{e_{j}}\cdot {\begin{pmatrix}v_{\text{neu}}\\1\end{pmatrix}}}$ berechnet. Dieser Fehler wird in eine Prioritätswarteschlange eingetragen, wobei der niedrigste Fehler das Top-Element der Warteschlange ist.

Kontrahieren der Kante

Aus der Warteschlange wird das Top-Element entfernt. Dieses Element entspricht derjenigen Kante, deren Kontraktion den geringsten Fehler verursacht. Die Kante wird nun kontrahiert, und dem dabei neu entstehenden Punkt v_neu wird die Fehlerquadrik der dabei zerstörten Kante zugewiesen. Dies bewirkt, dass der Fehler in den folgenden Iterationsschritten nicht nur von dem jeweils aktuellen Dreiecksnetz abhängt, sondern auch die vorherigen Veränderungen mit einbezogen werden.

Referenzen

1. M. Garland and P. Heckbert. Surface Simplification Using Quadric Error Metrics. In Proceedings of SIGGRAPH 97. (PDF)