public void TestVoronoiDiagramBuilder()
{
const string wkt =
@"POLYGON ((561758.3399999999 4815264.26, 561758.2199999997 4815264.24, 561758.0899999999 4815264.25, 561757.9699999997 4815264.279999999, 561757.8499999996 4815264.34, 561757.7599999999 4815264.42, 561757.6799999997 4815264.52, 561757.6299999999 4815264.640000001, 561757.5999999996 4815264.76, 561757.5999999996 4815264.890000001, 561757.6299999999 4815265.02, 561757.6799999997 4815265.130000001, 561757.7599999999 4815265.230000001, 561757.8600000005 4815265.310000001, 561758.6100000005 4815265.724000001, 561759.3600000005 4815266.138000001, 561760.1100000005 4815266.552000001, 561760.8600000005 4815266.966000001, 561761.6100000005 4815267.380000001, 561761.75 4815267.41, 561761.9000000004 4815267.4, 561762.04 4815267.359999999, 561762.1699999999 4815267.289999999, 561762.2699999996 4815267.1899999995, 561762.3600000005 4815267.07, 561762.4100000003 4815266.930000001, 561762.4299999997 4815266.789999999, 561762.4199999999 4815266.65, 561762.3799999999 4815266.51, 561762.2999999998 4815266.380000001, 561761.5319999999 4815265.966000001, 561760.764 4815265.552000001, 561759.9959999999 4815265.138000001, 561759.228 4815264.724000001, 561758.4600000001 4815264.310000001, 561758.3399999999 4815264.26))";
var reader = new WKTReader(GeometryFactory.Default);
var geom = reader.Read(wkt);
Assert.IsTrue(geom.IsValid);
Assert.IsTrue(geom.IsSimple);
//var dtb = new DelaunayTriangulationBuilder();
//dtb.SetSites(geom);
QuadEdgeSubdivision qes = null;
//Assert.DoesNotThrow( () => qes = dtb.GetSubdivision());
//Assert.IsNotNull(qes);
//foreach (var t in qes.GetVoronoiCellPolygons(GeometryFactory.Default))
// Console.WriteLine(t.AsText());
//Console.WriteLine();
IGeometryCollection diag = null;
//Assert.DoesNotThrow(() => diag = qes.GetVoronoiDiagram(GeometryFactory.Default));
//Assert.IsNotNull(diag);
var vdb = new VoronoiDiagramBuilder();
vdb.SetSites(geom);
qes = null;
Assert.DoesNotThrow(() => qes = vdb.GetSubdivision());
Assert.IsNotNull(qes);
foreach (var t in qes.GetVoronoiCellPolygons(GeometryFactory.Default))
Console.WriteLine(t.AsText());
Console.WriteLine();
diag = null;
Assert.DoesNotThrow(() => diag = vdb.GetDiagram(GeometryFactory.Default));
Assert.IsNotNull(diag);
}
public static IGeometry VoronoiDiagramWithData(IGeometry geom, IGeometry g2)
{
GeometryDataUtil.SetComponentDataToIndex(geom);
VertexTaggedGeometryDataMapper mapper = new VertexTaggedGeometryDataMapper();
mapper.LoadSourceGeometries(geom);
VoronoiDiagramBuilder builder = new VoronoiDiagramBuilder();
builder.SetSites(mapper.Coordinates);
if (g2 != null)
builder.ClipEnvelope = g2.EnvelopeInternal;
builder.Tolerance = TriangulationTolerance;
IGeometry diagram = builder.GetDiagram(geom.Factory);
mapper.TransferData(diagram);
return diagram;
}
public void Generate(int pointCount)
{
//Generate random points
var hashSet = new HashSet<Coordinate>();
{
var rand = new Random(seed);
while (hashSet.Count < pointCount)
{
var x = rand.NextDouble()*2 - 1;
var y = rand.NextDouble()*2 - 1;
if (x < -0.99 || y < -0.99 || x > 0.99 || y > 0.99) continue;
hashSet.Add(new Coordinate(x, y));
}
}
//Optimize points
{
var points = hashSet.ToArray();
for (var i = 0; i < 2; i++)
{
var builder = new VoronoiDiagramBuilder();
builder.SetSites(points);
VoronoiDiagram = builder.GetDiagram(new GeometryFactory());
for (var j = 0; j < points.Length; j++)
{
var poly = VoronoiDiagram[j] as Polygon;
points[j] = new Coordinate(poly.Centroid.X, poly.Centroid.Y);
}
}
}
//Build graph
PlanarGraph graph;
{
VoronoiDiagram = ClipGeometryCollection(VoronoiDiagram, new Envelope(-1, 1, -1, 1));
graph = new PlanarGraph(new OverlayNodeFactory());
var edges = new List<Edge>();
for (var i = 0; i < VoronoiDiagram.Count; i++)
{
var poly = VoronoiDiagram[i] as Polygon;
var coords = poly.Coordinates;
for (var j = 1; j < coords.Length; j++)
{
edges.Add(new Edge(new[] {coords[j - 1], coords[j]}, new Label(Location.Boundary)));
}
}
graph.AddEdges(edges);
}
//Convert graph
Dictionary<Node, MapNode> nodeDict;
{
var polys = new Dictionary<MapPolygon, HashSet<MapPolygon>>();
nodeDict = new Dictionary<Node, MapNode>();
var dats = new Dictionary<MapNode, Tuple<HashSet<MapPolygon>, HashSet<MapEdge>>>();
for (var i = 0; i < VoronoiDiagram.Count; i++)
{
var nodes = new List<MapNode>();
var poly = new MapPolygon
{
CentroidX = VoronoiDiagram[i].Centroid.X,
CentroidY = VoronoiDiagram[i].Centroid.Y,
Polygon = VoronoiDiagram[i] as Polygon
};
foreach (var j in VoronoiDiagram[i].Coordinates.Skip(1))
{
var n = graph.Find(j);
MapNode mapNode;
if (!nodeDict.TryGetValue(n, out mapNode))
{
mapNode = new MapNode {X = j.X, Y = j.Y};
dats[mapNode] =
new Tuple<HashSet<MapPolygon>, HashSet<MapEdge>>(new HashSet<MapPolygon> {poly},
new HashSet<MapEdge>());
}
else
dats[mapNode].Item1.Add(poly);
nodes.Add(nodeDict[n] = mapNode);
}
poly.Nodes = nodes.ToArray();
polys.Add(poly, new HashSet<MapPolygon>());
}
foreach (var i in nodeDict)
{
foreach (var j in dats[i.Value].Item1)
foreach (var k in dats[i.Value].Item1)
if (j != k)
{
polys[j].Add(k);
polys[k].Add(j);
}
foreach (var j in i.Key.Edges)
{
var from = nodeDict[graph.Find(j.Coordinate)];
var to = nodeDict[graph.Find(j.DirectedCoordinate)];
dats[from].Item2.Add(new MapEdge {From = from, To = to});
}
}
var ftrh = dats.Count(_ => _.Value.Item2.Count == 0);
foreach (var i in dats)
i.Key.Edges = i.Value.Item2.ToArray();
var x = polys.ToArray();
for (var i = 0; i < x.Length; i++)
{
x[i].Key.Neighbour = x[i].Value.ToArray();
x[i].Key.Id = i;
}
Polygons = x.Select(_ => _.Key).ToArray();
}
//Generate map
DetermineLandmass();
FindOceans();
ComputeDistances();
RedistributeElevation(nodeDict.Values);
FindLakesAndCoasts();
}
public static IGeometry VoronoiDiagram(IGeometry geom, IGeometry g2)
{
VoronoiDiagramBuilder builder = new VoronoiDiagramBuilder();
builder.SetSites(geom);
if (g2 != null)
builder.ClipEnvelope = g2.EnvelopeInternal;
builder.Tolerance = TriangulationTolerance;
IGeometry diagram = builder.GetDiagram(geom.Factory);
return diagram;
}
public static IGeometry voronoiDiagram(IGeometry geom, IGeometry g2)
{
var builder = new VoronoiDiagramBuilder();
builder.SetSites(geom);
if (g2 != null)
builder.ClipEnvelope = g2.EnvelopeInternal;
builder.Tolerance = TRIANGULATION_TOLERANCE;
IGeometry diagram = builder.GetDiagram(geom.Factory);
return diagram;
}
