Tutorial Separating Axis Theorem: Unterschied zwischen den Versionen

Aus DGL Wiki
Wechseln zu: Navigation, Suche
K (Neue Datei Links)
K (Rechtschreibfehler ausgebessert)
 
(3 dazwischenliegende Versionen von 3 Benutzern werden nicht angezeigt)
Zeile 106: Zeile 106:
  
 
'''Und hier das Beispielprojekt:'''
 
'''Und hier das Beispielprojekt:'''
* [[Datei:tut_SAT_exe|Windows-Binary]]
+
* {{ArchivLink|file=tut_SAT_exe|text=Windows-Binary}}
* [[Datei:tut_SAT_src|Delphi-Quelltext]]
+
* {{ArchivLink|file=tut_SAT_src|text=Delphi-Quelltext}}
  
 
=== Beispiel 2 - Kreis <> Polygon ===  
 
=== Beispiel 2 - Kreis <> Polygon ===  
Zeile 118: Zeile 118:
  
 
'''Und hier das Beispielprojekt für die Kollision Kreis <> Polygon:'''
 
'''Und hier das Beispielprojekt für die Kollision Kreis <> Polygon:'''
* [[Datei:tut_SAT_Kreis_exe|Windows-Binary]]
+
* {{ArchivLink|file=tut_SAT_Kreis_exe|text=Windows-Binary}}
* [[Datei:tut_SAT_Kreis_src|Delphi-Quelltext]]
+
* {{ArchivLink|file=tut_SAT_Kreis_src|text=Delphi-Quelltext}}
  
 
=== Beispiel 3 - Trennung von Polygonen ===
 
=== Beispiel 3 - Trennung von Polygonen ===
Zeile 131: Zeile 131:
  
 
'''Hier gibt es die Exe und den Code für die Trennung von Polygonen:'''
 
'''Hier gibt es die Exe und den Code für die Trennung von Polygonen:'''
* [[Datei:tut_SAT_Trennung_exe|Windows-Binary]]
+
* {{ArchivLink|file=tut_SAT_Trennung_exe|text=Windows-Binary}}
* [[Datei:tut_SAT_Trennung_src|Delphi-Quelltext]]
+
* {{ArchivLink|file=tut_SAT_Trennung_src|text=Delphi-Quelltext}}
  
 
=== Beispiel 4 - Rotation ===
 
=== Beispiel 4 - Rotation ===
Zeile 149: Zeile 149:
  
 
'''Hier gibt es den Code für die Behandlung von Rotation:'''
 
'''Hier gibt es den Code für die Behandlung von Rotation:'''
* [[Datei:tut_SAT_Rotation_src|Delphi-Quelltext]]
+
* {{ArchivLink|file=tut_SAT_Rotation_src|text=Delphi-Quelltext}}
 +
 
 
== Links ==
 
== Links ==
  
SAT-Tutorial(Eng) [http://www.harveycartel.org/metanet/tutorials/tutorialA.html]
+
SAT-Tutorial(Eng) [http://www.metanetsoftware.com/technique/tutorialA.html]
  
 
SAT-Tutorial(Eng / VB) [http://gpwiki.org/index.php/VB:Tutorials:Building_A_Physics_Engine:Basic_Intersection_Detection]
 
SAT-Tutorial(Eng / VB) [http://gpwiki.org/index.php/VB:Tutorials:Building_A_Physics_Engine:Basic_Intersection_Detection]
Zeile 158: Zeile 159:
 
== Nachwort ==
 
== Nachwort ==
  
Ich hoffe das Tutorial war nicht zu trocken und hat vielleicht auch ein wenig Spass gemacht und weitergeholfen. Für Fragen, Vorschläge, Ergänzungen, etc. bin ich selbstverständlich offen.
+
Ich hoffe das Tutorial war nicht zu trocken und hat vielleicht auch ein wenig Spaß gemacht und weitergeholfen. Für Fragen, Vorschläge, Ergänzungen, etc. bin ich selbstverständlich offen.
  
 
mfg
 
mfg

Aktuelle Version vom 25. September 2013, 15:32 Uhr

Kollisionserkennung

von Polygonen mit dem Separating Axis Theorem

Vorwort

In diesem Tutorial möchte ich eine schnelle Variante zur Kollision zweier konvexer Polygone erläutern. Diese kann nachträglich für konkave Polygone und andere Objekte wie Kreise und abgerundete Objekte verwendet werden. Um den mathematischen Hintergrund zu verstehen, ist es sinnvoll ein wenig Kenntnis in Vektorrechnung zu haben.

Kollision zweier Polygone

Die Theorie

Normale der Seite des Quadrats

Das Separating Axis Theorem (kurz: SAT) besagt, dass zwei Polygone sich nicht schneiden, wenn es möglich ist, eine Gerade zu finden, die zwischen den beiden liegt, bzw. die beiden trennt. Der Begriff "Gerade" ist hier allerdings nicht ganz korrekt. Eine Gerade hat eine räumliche Lage, diese ist für unser Vorhaben jedoch nicht von Belangen, deshalb werde ich im Folgenden das Wort Achse benutzen, daher auch der Name Separating Axis (Trennungsachse). Nun gibt es unendlich viele Achsen die man testen könnte... Glücklicherweise kann man sich hier auf eine überschaubare Zahl beschränken, denn man braucht nur die Anzahl der Seiten beider Polygone. Bei einem Viereck wären das vier, bei einem Dreieck drei, etc. Hat man die Eckpunkte des Polygons als Vektoren (Ortsvektoren) gegeben, kann man durch Subtraktion zweier Ortsvektoren den Vektor bestimmen der zu der Seite gehört, die von den beiden Vektoren aufgespannt wird.

Auf dem Bild rechts ist ein Beispiel zu sehen. Dort ist ein Quadrat, die grünen Striche bezeichnen die Ortsvektoren der Eckpunkte, der rote Strich ist die berechnete Seite. Was wir aber brauchen ist der blaue Strich, das ist die Normale der Seite. Die Normale berechnet sich folgendermaßen:

SAT Normale Formel.jpg

die Koordinaten werden vertauscht und eine von beiden wird negiert, so erhält man einen Vektor der senkrecht zu dem original Vektor ist. Da es zwei Möglichkeiten gibt, spricht man von der linken oder rechten Normale (aus der Sicht des Vektors). Der Vektor muss dann noch normiert werden, sodass er die Länge 1 erhält,. Jetzt müssen beide Polygone auf diesen Vektor projiziert werden, denn dadurch erhalten wir ein eindimensionales Abbild unserer Polygone und können mittels eines einfachen Vergleichs überprüfen, ob sich die beiden eindimensionalen Strecken schneiden. Sollte ein Fall eintreffen bei dem kein Schnitt stattfindet, dann kollidieren die beiden Polygone nicht und die Prozedur kann abgebrochen werden. Dies ist auch der Grund, warum dieses Verfahren recht schnell ist, da im besten Fall schon im ersten Durchlauf abgebrochen werden kann.

Nun zur Projektion:

SAT Kollision.jpgSAT Keine Kollision.jpg

Auf dem linken Bild sieht man, wie beide Polygone auf die Achse projiziert werden, der pinke Bereich zeigt die Schnittmenge an. In diesem Fall ist die Achse die Normale der linken oder rechten Seite des Quadrats. Auf dem Bild rechts ist der Fall dargestellt, dass keine Kollision stattfindet, demzufolge gibt es auch keine Schnittmenge auf der Achse. Somit gibt es eine Gerade, bzw. eine Trennungsachse senkrecht zur unserer Projektionsachse, welche die beiden Polygone separiert.

Für die Projektion eines Vektors auf einen anderen verwenden wir das Skalarprodukt, bei diesem kommt ein Zahlenwert heraus, der die Position des Eckpunktes auf unserer eindimensionalen Achse darstellt. Haben wir sämtliche Punkte projiziert, so müssen wir für die jeweiligen Polygone noch jeweils den kleinsten und größten Wert heraussuchen, damit wir zwei Strecken erhalten. Diese werden dann mit einer einfachen Abfrage auf Schnitt geprüft und das wars.

Zusammenfassung

  • Jedes der beiden Polygone durchgehen und alle nötigen Achsen aus den Normalen der Seitenflächen bestimmen
    • Jeden Eckpunkt jedes Polygons auf diese Achsen projizieren
    • Die kleinsten und größten Werte ermitteln und auf Schnitt prüfen
  • Tritt der Fall auf, dass kein Schnitt statt findet, so kann sofort abgebrochen werden, es findet keine Kollision statt.

Kollision eines Kreises und eines Polygons

Theorie

Das Prinzip für die Kollision zweier Polygone ist denke ich jetzt klar geworden, doch was ist, wenn wir einen Kreis haben, der mit einem Polygon kollidiert? Dieser Fall ist leicht abzuhandeln, ein Kreis hat unendlich viele Normalen, die man testen könnte, uns reichen aber die, die die Vertices des Polygons schneiden würden, sprich: die Achsen, die den Kreismittelpunkt und die Ecken unseres Polygons verbinden würden.

Projektionsachsen eines Kreises

Die blauen Linien auf dem Bild sind wieder einmal die Achsen auf die wir projizieren, diese kommen zu den, die wir aus dem Polygon berechnen hinzu.

Zusammenfassung

  • Wir berechnen also den Vektor vom Kreis-Mittelpunkt zum Vertex, dieser wird normiert.
  • Dann projizieren wir das Polygon wie gehabt
  • Der Kreis wird projiziert, indem der Vektor, auf den wir projizieren, mit dem Radius des Kreises skaliert (also skalar multipliziert) wird. Mit diesem wird dann genauso weiter verfahren. Dies ist dann der Max-Wert für unseren Kreis, der Min-Wert ist einfach der Max-Wert * -1, also -max.

Hierzu ist zu sagen, dass dies noch dahingehend optimieren kann, dass man feststellt, in welcher Voronoi-Region des Polygons sich der Kreis befindet. Wenn man diese Region hat, reicht es diese eine Achse zu prüfen.

Kollision eines Punktes und eines Polygons

Ganz am Rande möchte ich noch erwähnen, dass es ebenso möglich ist, zu prüfen, ob sich ein Punkt in einem Polygon befindet. Hierzu wird einfach der Punkt auf die Achsen projiziert und geprüft ob er größer als das Minimum und kleiner als das Maximum des Polygons auf dieser Achse ist.

Polygone trennen

Damit unsere Kollisionserkennung praxistauglich wird, müssen die Polygone, wenn sie kollidieren auch wieder getrennt werden können. Hierzu benötigen wir einen Vektor, der unsere beiden Polygone wieder auseinander "schiebt". Selbstverständlich könnte man einen beliebigen Vektor nehmen, aber das Ergebnis wäre eher realtitätsfern, deshalb brauchen wir den Vektor, der den kürzesten Weg beschreibt, um die beiden Polygone zu trennen. Dieser Vektor wird auch die Normale der Kollision oder MTD-Vektor genannt (für Minimum Translation Distance). Hier kommt ein weiterer Vorteil des Separating Axis Theorems zum tragen, denn den Vektor den wir suchen, haben wir schon so gut wie berechnet. Wir multiplizieren alle unsere Projektionsachsen mit den Differenzen, die wir bei den eindimensionalen Kollisionen herausbekommen haben, somit erhalten wir Vektoren, die für das Außeinander-Schieben unserer Polygonein in Frage kommen. Nun müssen wir nur noch den kürzesten davon finden, dies ist unser MTD. Allerdings sollte man vor der Verwendung sicherstellen, dass dieser Vektor vom Polygon wegzeigt. Hierzu können wir ebenfalls das Skalarprodukt nutzen.

Rotation

Um Polygone behandeln zu können, die rotieren, sind lediglich einige Änderungen an der Polygonklasse erforderlich. Hierbei gibt es mehrere Optionen, zum einen können wir unsere Vertices in jedem Schritt drehen und speichern, allerdings kann dies auf Dauer zu Ungenauigkeiten führen. Eine zweite Möglichkeit wäre, bei jedem Zugriff auf die Vertices, einen gedrehten Vektor auszugeben. Dies erhöht zwar die benötigte Rechenleistung, ist aber deutlich genauer. Das Seperation Axis Theorem kann dann mit den gedrehten Vertices weiter verfahren.

Das Beispielprojekt

Beispiel 1 - Polygon <> Polygon

Was wäre ein Tutorial doch ohne Beispiel ;) Ich werde hier nur kurz die Verwendung des Codes erläutern und ein kleines Beispielprogramm anhängen. Ein Polygon muss natürlich erst einmal erzeugt werden:

A := TPolygon.Create;

Die einfachste Möglichkeit es zu gestalten funktioniert so:

  with A do
  begin
    position := to_v2f(200, 200);
    AddVertex(to_v2f(50, -50));
    AddVertex(to_v2f(50, 50));
    AddVertex(to_v2f(-50, 50));
    AddVertex(to_v2f(-50, -50));
  end;

Dies liefert ein Quadrat mit den Maßen 100*100 an der Position (200|200). Die Koordinaten der Vertices werden absolut zur Position und entgegen des Uhrzeigersinns angegeben.

Und hier das Beispielprojekt:

Beispiel 2 - Kreis <> Polygon

Für den Kreis ist die Verwendung im Prinzip die selbe, mit:

  C := TCircle.Create;

wird zunächst ein Kreis erzeugt und mit Position und Padius bekommt er seine Werte zugewiesen.

Und hier das Beispielprojekt für die Kollision Kreis <> Polygon:

Beispiel 3 - Trennung von Polygonen

In diesem Beispiel wird ein Polygon mittels:

  B.position := v2f_sub(B.position, MTD);

verschoben, sodass die beiden Polygone sich nicht schneiden. Es ließen sich auch beide um die Hälfte des Vektors verschieben, dies hängt von der Anwendung ab.

Hier gibt es die Exe und den Code für die Trennung von Polygonen:

Beispiel 4 - Rotation

Die Rotation des Polygons erfolgt beim Auslesen der Vertex Koordinaten mittels einer Rotationsmatrix:

function TPolygon.GetVertex(n: integer): TVector2f;
var
  acos, asin: extended;
begin
  sincos(degtorad(fangle),asin,acos);
  result.x := fvertices[n].x * acos + fvertices[n].y * -asin;
  result.y := fvertices[n].x * asin + fvertices[n].y * acos;
end;

Hier gibt es den Code für die Behandlung von Rotation:

Links

SAT-Tutorial(Eng) [1]

SAT-Tutorial(Eng / VB) [2]

Nachwort

Ich hoffe das Tutorial war nicht zu trocken und hat vielleicht auch ein wenig Spaß gemacht und weitergeholfen. Für Fragen, Vorschläge, Ergänzungen, etc. bin ich selbstverständlich offen.

mfg

Seth



Nächstes Tutorial:
Tutorial Kollision1

Schreibt was ihr zu diesem Tutorial denkt ins Feedbackforum von DelphiGL.com.
Lob, Verbesserungsvorschläge, Hinweise und Tutorialwünsche sind stets willkommen.