#!/usr/bin/env python3 # -*- coding: utf-8 -*- """ @author: huiming zhou """ import queue import environment import tools class Astar: def __init__(self, Start_State, Goal_State, n, m, heuristic_type): self.xI = Start_State self.xG = Goal_State self.u_set = environment.motions # feasible input set self.obs_map = environment.map_obs() # position of obstacles self.n = n self.m = m self.heuristic_type = heuristic_type def searching(self): """ Searching using A_star. :return: planning path, action in each node, visited nodes in the planning process """ q_astar = queue.QueuePrior() # priority queue q_astar.put(self.xI, 0) parent = {self.xI: self.xI} # record parents of nodes actions = {self.xI: (0, 0)} # record actions of nodes cost = {self.xI: 0} visited = [] while not q_astar.empty(): x_current = q_astar.get() visited.append(x_current) # record visited nodes if x_current == self.xG: # stop condition break for u_next in self.u_set: # explore neighborhoods of current node x_next = tuple([x_current[i] + u_next[i] for i in range(len(x_current))]) # if neighbor node is not in obstacles -> ... if 0 <= x_next[0] < self.n and 0 <= x_next[1] < self.m \ and not tools.obs_detect(x_current, u_next, self.obs_map): new_cost = cost[x_current] + int(self.get_cost(x_current, u_next)) if x_next not in cost or new_cost < cost[x_next]: # conditions for updating cost cost[x_next] = new_cost priority = new_cost + self.Heuristic(x_next, self.xG, self.heuristic_type) q_astar.put(x_next, priority) # put node into queue using priority "f+h" parent[x_next] = x_current actions[x_next] = u_next [path_astar, actions_astar] = tools.extract_path(self.xI, self.xG, parent, actions) return path_astar, actions_astar, visited def get_cost(self, x, u): """ Calculate cost for this motion :param x: current node :param u: input :return: cost for this motion :note: cost function could be more complicate! """ return 1 def Heuristic(self, state, goal, heuristic_type): """ Calculate heuristic. :param state: current node (state) :param goal: goal node (state) :param heuristic_type: choosing different heuristic functions :return: heuristic """ if heuristic_type == "manhattan": return abs(goal[0] - state[0]) + abs(goal[1] - state[1]) elif heuristic_type == "euclidean": return ((goal[0] - state[0]) ** 2 + (goal[1] - state[1]) ** 2) ** (1 / 2) else: print("Please choose right heuristic type!") if __name__ == '__main__': x_Start = (15, 10) # Starting node x_Goal = (48, 15) # Goal node astar = Astar(x_Start, x_Goal, environment.col, environment.row, "manhattan") [path_astar, actions_astar, visited_astar] = astar.searching() tools.showPath(x_Start, x_Goal, path_astar, visited_astar, 'Astar_searching') # Plot path and visited nodes