Skip to content

pedrohm/Computational-geometry

 
 

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

73 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

Computational Geometry Unity Library

This library consists of two folders. The idea is that one is for testing purposes and the other folder is the folder you drag into your project.

Make sure all input coordinates are normalized to range 0-1 to avoid floating point precision issues! Normalizing methods exists in HelpMethods. This is not always needed but if you notice that an algorithm doesn't work, try to normalize the input coordinates.

Some of these algorithms are available in tutorial form here: https://www.habrador.com/tutorials/math/

Finished

1. Intersection

1.1 2d-space

Point-triangle

Intersection point-triangle

Point-polygon. Suffers from floating point precision issues

Intersection point-polygon

Triangle-triangle

Intersection triangle-triangle

AABB-AABB

Intersection point-triangle

Line-line

Intersection line-line

Ray-plane

Intersection ray-plane

Line-plane

Intersection line-plane

Plane-plane

Intersection plane-plane

Point-circle

Intersection point-circle

2. Generate mesh

Grid mesh

Mesh grid

Mesh shapes: Arrow, circles, lines

Mesh shapes

3. Convex Hull

*Jarvis March. Is also known as "Gift wrapping"

Convex hull jarvis march

4. Triangulation

Triangulate convex polygon.

You have points on a convex hull you want to triangulate. You have four options here if you have colinear points (points on the same line):

  1. Triangulate the convex hull while ignoring the colinear points. The area covered will be the same anyway.
  2. Triangulate the convex hull and add the colinear points by splitting triangle edges.
  3. Add a point inside of the convex hull.
  4. Use the algorithm below called "Triangulate points with 'visible edge' algorithm."

Triangulation convex polygon

Triangulate points with "visible edge" algorithm.

You have some points you want to triangulate, you follow the steps:

  1. Sort all points in x and then y direction
  2. Find the first triangle
  3. Add the rest of the sorted points one-by-one and build triangles to visible edges on the existing triangulation. To determine if an edge is visible from the point you build the convex shape from the existing triangles. Then for each edge in the convex hull, you build a triangle with the point. If this triangle is oriented clockwise, the edge is visible and you can add a new triangle.

A visualization of this algorithm can be found here: https://www.youtube.com/watch?v=MkMXKu1m6A4

Triangulation visible edges

Triangulate points with "point-by-point" algorithm.

You have some points you want to triangulate, you follow the steps:

  1. Generate the convex hull of all points.
  2. Triangulate the convex hull with one of several algorithms mentioned above.
  3. Add the rest of the points one-by-one by splitting the triangles they end up in into three new triangles.

Triangulation point-by-point

4.1 Delaunay triangulation

"point-by-point" method

You generate a big triangle around all points you want to triangulate. Then you add each point one after the other. The triangle the point ends up in is split into three new triangles. After the split you restore the Delaunay triangulation by flipping edges. When all points have been added you remove the remains of the first big triangle. A visualization of this algorithm can be found here: https://www.youtube.com/watch?v=YNQR5tH-s40

Triangulation Delaunay point-by-point

"flip edges" method

You triangulate the points by using a "bad" triangulation method (which is in this case either "visible edge" or "point-by-point" from above). Then you go through all edges and check if the edge should be flipped to make a better triangle. When no more edges can be flipped you are done! A visualization of this algorithm can be found: https://www.youtube.com/watch?v=-d7Nb4fxL5s and https://www.youtube.com/watch?v=lR_SzgEkDwk

Triangulation Delaunay flip edges

Constrained triangulation

You add the constraints to the points and generate a Delaunay triangulation by using one of the above methods. Use this triangulation to find which edges interesect with the constraints. Then you flip these edges until they no longer interesect with the constraint. You finally remove the triangles that are "inside" of the constraint.

Triangulation Delaunay constrained

5. Voronoi diagram

From a Delaunay triangulation

You first generate a Delaunay triangulation by using some method. Then you use the fact that you can get the Voronoi diagram from the Delaunay triangulation.

Voronoi from delaunay

6. Polygon clipping

Greiner-Hormann method

Polygon clipping greiner

Sutherland-Hodgman method

Polygon clipping sutherland

7. Other

Is a triangle oriented clockwise?

Triangle orientation

Is a point left, on, or right of vector?

Point vector orientation

Is a point left, on, or right of a plane? Which is the same as the distance to the plane.

Point plane orientation

Is a quadrilateral convex?

Quadrilateral convex or concave

Is a point between two other points on the same line?

Is point between points on line segment

Closest point on a line-segment?

Closest point on line segment

TODO

Algorithms to implement

  • Dynamic constrained delaunay triangulation
  • Convex hull: Quickhull from the Valve paper
  • Convex hull: Graham scan
  • Marching cubes
  • Cut 3d mesh with plane
  • Metaballs
  • Voronoi with Fortune's algorithm
  • Voronoi point-by-point
  • Triangulation concave polygon by ear clipping
  • Rectangle-rectangle with SAT
  • Triangulate with marching squares

Stuff to fix

  • Optimize Constrained Delaunay - there's a faster method to find edges that intersects with the constrained edge. I also think the method where triangles within the constrain is removed can be faster.
  • Make a test scene to test that the "find which triangle a point is in by triangulation walk" is working

Socials

Follow me on Twitter for more Unity stuff: https://twitter.com/eriknordeus

About

Computational Geometry Unity library with implementations of intersection algorithms, triangulations like delaunay, voronoi diagrams, polygon clipping, etc

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages

  • C# 100.0%