using astar.PathingHelper; using Graph; using Microsoft.Extensions.Logging; using Graph.Utils; using OSM_Graph.Enums; using OSM_Regions; namespace astar { public static class Astar { public static Route FindPath(float startLat, float startLon, float endLat, float endLon, float regionSize, bool car = true, string? importFolderPath = null, ILogger? logger = null) { RegionLoader rl = new(regionSize, importFolderPath, logger: logger); Graph graph = Spiral(rl, startLat, startLon, regionSize); Graph endRegion = Spiral(rl, endLat, endLon, regionSize); graph.ConcatGraph(endRegion); KeyValuePair startNode = graph.ClosestNodeToCoordinates(startLat, startLon, car); startNode.Value.PreviousIsFromStart = true; startNode.Value.PreviousNodeId = startNode.Key; startNode.Value.Distance = 0f; KeyValuePair endNode = graph.ClosestNodeToCoordinates(endLat, endLon, car); endNode.Value.PreviousIsFromStart = false; endNode.Value.PreviousNodeId = endNode.Key; endNode.Value.Distance = 0f; double totalDistance = NodeUtils.DistanceBetween(startNode.Value, endNode.Value); PriorityHelper priorityHelper = new(totalDistance, SpeedHelper.GetMaxSpeed(car)); logger?.Log(LogLevel.Information, "From {0:00.00000}#{1:000.00000} to {2:00.00000}#{3:000.00000} Great-Circle {4:00000.00}km", startNode.Value.Lat, startNode.Value.Lon, endNode.Value.Lat, endNode.Value.Lon, totalDistance / 1000); PriorityQueue toVisitStart = new(); toVisitStart.Enqueue(startNode.Key, 0); PriorityQueue toVisitEnd = new(); toVisitEnd.Enqueue(endNode.Key, 0); while (toVisitStart.Count > 0 && toVisitEnd.Count > 0) { ulong currentNodeStartId = toVisitStart.Dequeue(); Node currentNodeStart = graph.Nodes[currentNodeStartId]; foreach ((ulong neighborId, KeyValuePair wayId) in currentNodeStart.Neighbors) { if (!graph.ContainsNode(neighborId)) graph.ConcatGraph(Graph.FromGraph(rl.LoadRegionFromNodeId(neighborId))); if (!graph.ContainsWay(wayId.Key)) { foreach (global::Graph.Graph? g in rl.LoadRegionsFromWayId(wayId.Key)) graph.ConcatGraph(Graph.FromGraph(g)); } OSM_Graph.Way way = graph.Ways[wayId.Key]; byte speed = SpeedHelper.GetSpeed(way, car); if(speed < 1) continue; if(wayId.Value && way.GetDirection() == WayDirection.Forwards && car) continue; if(!wayId.Value && way.GetDirection() == WayDirection.Backwards && car) continue; Node neighborNode = graph.Nodes[neighborId]; if (neighborNode.PreviousIsFromStart is false)//Check if we found the opposite End return PathFound(graph, currentNodeStart, neighborNode, logger); float distance = (currentNodeStart.Distance??float.MaxValue) + (float)currentNodeStart.DistanceTo(neighborNode); if (neighborNode.PreviousNodeId is null || neighborNode.Distance > distance) { neighborNode.PreviousNodeId = currentNodeStartId; neighborNode.Distance = distance; neighborNode.PreviousIsFromStart = true; toVisitStart.Enqueue(neighborId, priorityHelper.CalculatePriority(currentNodeStart, neighborNode, endNode.Value, speed)); } logger?.LogTrace($"Neighbor {neighborId} {neighborNode}"); } ulong currentNodeEndId = toVisitEnd.Dequeue(); Node currentNodeEnd = graph.Nodes[currentNodeEndId]; foreach ((ulong neighborId, KeyValuePair wayId) in currentNodeEnd.Neighbors) { if (!graph.ContainsNode(neighborId)) graph.ConcatGraph(Graph.FromGraph(rl.LoadRegionFromNodeId(neighborId))); if (!graph.ContainsWay(wayId.Key)) { foreach (global::Graph.Graph? g in rl.LoadRegionsFromWayId(wayId.Key)) graph.ConcatGraph(Graph.FromGraph(g)); } OSM_Graph.Way way = graph.Ways[wayId.Key]; byte speed = SpeedHelper.GetSpeed(way, car); if(speed < 1) continue; if(wayId.Value && way.GetDirection() == WayDirection.Backwards && car) continue; if(!wayId.Value && way.GetDirection() == WayDirection.Forwards && car) continue; Node neighborNode = graph.Nodes[neighborId]; if (neighborNode.PreviousIsFromStart is true)//Check if we found the opposite End return PathFound(graph, neighborNode, currentNodeEnd, logger); float distance = (currentNodeEnd.Distance??float.MaxValue) + (float)currentNodeEnd.DistanceTo(neighborNode); if (neighborNode.PreviousNodeId is null || neighborNode.Distance > distance) { neighborNode.PreviousNodeId = currentNodeEndId; neighborNode.Distance = distance; neighborNode.PreviousIsFromStart = false; toVisitEnd.Enqueue(neighborId, priorityHelper.CalculatePriority(currentNodeEnd, neighborNode, startNode.Value, speed)); } logger?.LogTrace($"Neighbor {neighborId} {neighborNode}"); } logger?.LogDebug($"Distance {currentNodeStart.DistanceTo(currentNodeEnd):000000.00}m toVisit-Queues: {toVisitStart.Count} {toVisitStart.UnorderedItems.MinBy(i => i.Priority).Priority} {toVisitEnd.Count} {toVisitEnd.UnorderedItems.MinBy(i => i.Priority).Priority}"); } return new Route(graph, Array.Empty().ToList(), false); } private static Route PathFound(Graph graph, Node fromStart, Node fromEnd, ILogger? logger = null) { logger?.LogInformation("Path found!"); List path = new(); path.Add(new Step((float)NodeUtils.DistanceBetween(fromStart, fromEnd), fromStart, fromEnd)); Node current = fromStart; while (current.Distance != 0f) { Step step = new((float)NodeUtils.DistanceBetween(graph.Nodes[(ulong)current.PreviousNodeId!], current), graph.Nodes[(ulong)current.PreviousNodeId!], current); path.Add(step); current = graph.Nodes[(ulong)current.PreviousNodeId!]; } path.Reverse();//Since we go from the middle backwards until here current = fromEnd; while (current.Distance != 0f) { Step step = new((float)NodeUtils.DistanceBetween(graph.Nodes[(ulong)current.PreviousNodeId!], current), current, graph.Nodes[(ulong)current.PreviousNodeId!]); path.Add(step); current = graph.Nodes[(ulong)current.PreviousNodeId!]; } Route r = new (graph, path, true); logger?.LogInformation(r.ToString()); return r; } private static Graph Spiral(RegionLoader loader, float lat, float lon, float regionSize) { Graph? ret = Graph.FromGraph(loader.LoadRegionFromCoordinates(lat, lon)); int iteration = 1; while (ret is null) { for (int x = -iteration; x <= iteration; x++) { Graph? g1 = Graph.FromGraph(loader.LoadRegionFromCoordinates(lat + x * regionSize, lon - iteration * regionSize)); Graph? g2 = Graph.FromGraph(loader.LoadRegionFromCoordinates(lat + x * regionSize, lon + iteration * regionSize)); if (ret is not null) { ret.ConcatGraph(g1); ret.ConcatGraph(g2); } else if (ret is null && g1 is not null) { ret = g1; ret.ConcatGraph(g2); }else if (ret is null && g2 is not null) ret = g2; } for (int y = -iteration + 1; y < iteration; y++) { Graph? g1 = Graph.FromGraph(loader.LoadRegionFromCoordinates(lat - iteration * regionSize, lon + y * regionSize)); Graph? g2 = Graph.FromGraph(loader.LoadRegionFromCoordinates(lat + iteration * regionSize, lon + y * regionSize)); if (ret is not null) { ret.ConcatGraph(g1); ret.ConcatGraph(g2); } else if (ret is null && g1 is not null) { ret = g1; ret.ConcatGraph(g2); }else if (ret is null && g2 is not null) ret = g2; } iteration++; } return ret; } } }