/// <summary>
/// If this is a Delaunay Triangulation of a pointset we need to fill so the triangle mesh gets a ConvexHull
/// </summary>
private static void FinalizationConvexHull(DTSweepContext tcx)
{
AdvancingFrontNode n1, n2, n3;
DelaunayTriangle t1;
TriangulationPoint first, p1;
n1 = tcx.Front.Head.Next;
n2 = n1.Next;
n3 = n2.Next;
first = n1.Point;
TurnAdvancingFrontConvex(tcx, n1, n2);
n1 = tcx.Front.Tail.Prev;
if (n1.Triangle.Contains(n1.Next.Point) && n1.Triangle.Contains(n1.Prev.Point))
{
t1 = n1.Triangle.NeighborAcrossFrom(n1.Point);
RotateTrianglePair(n1.Triangle, n1.Point, t1, t1.OppositePoint(n1.Triangle, n1.Point));
tcx.MapTriangleToNodes(n1.Triangle);
tcx.MapTriangleToNodes(t1);
}
n1 = tcx.Front.Head.Next;
if (n1.Triangle.Contains(n1.Prev.Point) && n1.Triangle.Contains(n1.Next.Point))
{
t1 = n1.Triangle.NeighborAcrossFrom(n1.Point);
RotateTrianglePair(n1.Triangle, n1.Point, t1, t1.OppositePoint(n1.Triangle, n1.Point));
tcx.MapTriangleToNodes(n1.Triangle);
tcx.MapTriangleToNodes(t1);
}
// TODO: implement ConvexHull for lower right and left boundary
// Lower right boundary
first = tcx.Front.Head.Point;
n2 = tcx.Front.Tail.Prev;
t1 = n2.Triangle;
p1 = n2.Point;
do
{
tcx.RemoveFromList(t1);
p1 = t1.PointCCWFrom(p1);
if (p1 == first)
{
break;
}
t1 = t1.NeighborCCWFrom(p1);
} while (true);
// Lower left boundary
first = tcx.Front.Head.Next.Point;
p1 = t1.PointCWFrom(tcx.Front.Head.Point);
t1 = t1.NeighborCWFrom(tcx.Front.Head.Point);
do
{
tcx.RemoveFromList(t1);
p1 = t1.PointCCWFrom(p1);
t1 = t1.NeighborCCWFrom(p1);
} while (p1 != first);
tcx.FinalizeTriangulation();
}
private static void FinalizationConvexHull(DTSweepContext tcx)
{
DelaunayTriangle triangle;
AdvancingFrontNode next = tcx.Front.Head.Next;
AdvancingFrontNode c = next.Next;
AdvancingFrontNode node3 = c.Next;
TriangulationPoint point = next.Point;
TurnAdvancingFrontConvex(tcx, next, c);
next = tcx.Front.Tail.Prev;
if (next.Triangle.Contains(next.Next.Point) && next.Triangle.Contains(next.Prev.Point))
{
triangle = next.Triangle.NeighborAcrossFrom(next.Point);
RotateTrianglePair(next.Triangle, next.Point, triangle, triangle.OppositePoint(next.Triangle, next.Point));
tcx.MapTriangleToNodes(next.Triangle);
tcx.MapTriangleToNodes(triangle);
}
next = tcx.Front.Head.Next;
if (next.Triangle.Contains(next.Prev.Point) && next.Triangle.Contains(next.Next.Point))
{
triangle = next.Triangle.NeighborAcrossFrom(next.Point);
RotateTrianglePair(next.Triangle, next.Point, triangle, triangle.OppositePoint(next.Triangle, next.Point));
tcx.MapTriangleToNodes(next.Triangle);
tcx.MapTriangleToNodes(triangle);
}
point = tcx.Front.Head.Point;
c = tcx.Front.Tail.Prev;
triangle = c.Triangle;
TriangulationPoint point2 = c.Point;
while (true)
{
tcx.RemoveFromList(triangle);
point2 = triangle.PointCCWFrom(point2);
if (point2 == point)
{
break;
}
triangle = triangle.NeighborCCWFrom(point2);
}
point = tcx.Front.Head.Next.Point;
point2 = triangle.PointCWFrom(tcx.Front.Head.Point);
triangle = triangle.NeighborCWFrom(tcx.Front.Head.Point);
do
{
tcx.RemoveFromList(triangle);
point2 = triangle.PointCCWFrom(point2);
triangle = triangle.NeighborCCWFrom(point2);
}while (point2 != point);
tcx.FinalizeTriangulation();
}
/// <summary>
/// If this is a Delaunay Triangulation of a pointset we need to fill so the triangle mesh gets a ConvexHull
/// </summary>
private static void FinalizationConvexHull( DTSweepContext tcx )
{
AdvancingFrontNode n1, n2;
DelaunayTriangle t1;
TriangulationPoint first, p1;
n1 = tcx.Front.Head.Next;
n2 = n1.Next;
first = n1.Point;
TurnAdvancingFrontConvex(tcx, n1, n2);
n1 = tcx.Front.Tail.Prev;
if (n1.Triangle.Contains(n1.Next.Point) && n1.Triangle.Contains(n1.Prev.Point)) {
t1 = n1.Triangle.NeighborAcrossFrom(n1.Point);
RotateTrianglePair(n1.Triangle, n1.Point, t1, t1.OppositePoint(n1.Triangle, n1.Point));
tcx.MapTriangleToNodes(n1.Triangle);
tcx.MapTriangleToNodes(t1);
}
n1 = tcx.Front.Head.Next;
if (n1.Triangle.Contains(n1.Prev.Point) && n1.Triangle.Contains(n1.Next.Point)) {
t1 = n1.Triangle.NeighborAcrossFrom(n1.Point);
RotateTrianglePair(n1.Triangle, n1.Point, t1, t1.OppositePoint(n1.Triangle, n1.Point));
tcx.MapTriangleToNodes(n1.Triangle);
tcx.MapTriangleToNodes(t1);
}
// TODO: implement ConvexHull for lower right and left boundary
// Lower right boundary
first = tcx.Front.Head.Point;
n2 = tcx.Front.Tail.Prev;
t1 = n2.Triangle;
p1 = n2.Point;
do {
tcx.RemoveFromList(t1);
p1 = t1.PointCCWFrom(p1);
if (p1 == first) break;
t1 = t1.NeighborCCWFrom(p1);
} while (true);
// Lower left boundary
first = tcx.Front.Head.Next.Point;
p1 = t1.PointCWFrom(tcx.Front.Head.Point);
t1 = t1.NeighborCWFrom(tcx.Front.Head.Point);
do {
tcx.RemoveFromList(t1);
p1 = t1.PointCCWFrom(p1);
t1 = t1.NeighborCCWFrom(p1);
} while (p1 != first);
tcx.FinalizeTriangulation();
}