-
Notifications
You must be signed in to change notification settings - Fork 27
/
GBFS.py
222 lines (167 loc) · 6.26 KB
/
GBFS.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
# -*- coding: utf-8 -*-
"""
Created on Thu Mar 30 16:45:58 2023
@author: HJ
"""
# 贪婪最佳优先搜索(Greedy Best First Search, GBFS)算法
# A*: F = G + H
# GBFS: F = H
# https://zhuanlan.zhihu.com/p/346666812
from functools import lru_cache
from common import *
Queue_Type = 0
"""
# OpenList 采用的 PriorityQueue 的结构
## 0 -> SetQueue
## 1 -> ListQueue
## 2 -> PriorityQueuePro
List/Set可以实现更新OpenList中节点的parent和cost, 找到的路径更优\n
PriorityQueuePro速度最快, 但无法更新信息, 路径较差\n
List速度最慢, Set速度接近PriorityQueuePro甚至更快\n
"""
# 地图读取
IMAGE_PATH = 'image.jpg' # 原图路径
THRESH = 172 # 图片二值化阈值, 大于阈值的部分被置为255, 小于部分被置为0
HIGHT = 350 # 地图高度
WIDTH = 600 # 地图宽度
MAP = GridMap(IMAGE_PATH, THRESH, HIGHT, WIDTH) # 栅格地图对象
# 起点终点
START = (290, 270) # 起点坐标 y轴向下为正
END = (298, 150) # 终点坐标 y轴向下为正
""" ---------------------------- Greedy Best First Search算法 ---------------------------- """
# F = H
# 设置OpenList使用的优先队列
if Queue_Type == 0:
NodeQueue = SetQueue
elif Queue_Type == 1:
NodeQueue = ListQueue
else:
NodeQueue = PriorityQueuePro
# 贪婪最佳优先搜索算法
class GBFS:
"""GBFS算法"""
def __init__(
self,
start_pos = START,
end_pos = END,
map_array = MAP.map_array,
move_step = 3,
move_direction = 8,
):
"""GBFS算法
Parameters
----------
start_pos : tuple/list
起点坐标
end_pos : tuple/list
终点坐标
map_array : ndarray
二值化地图, 0表示障碍物, 255表示空白, H*W维
move_step : int
移动步数, 默认3
move_direction : int (8 or 4)
移动方向, 默认8个方向
"""
# 网格化地图
self.map_array = map_array # H * W
self.width = self.map_array.shape[1]
self.high = self.map_array.shape[0]
# 起点终点
self.start = Node(*start_pos) # 初始位置
self.end = Node(*end_pos) # 结束位置
# Error Check
if not self._in_map(self.start) or not self._in_map(self.end):
raise ValueError(f"x坐标范围0~{self.width-1}, y坐标范围0~{self.height-1}")
if self._is_collided(self.start):
raise ValueError(f"起点x坐标或y坐标在障碍物上")
if self._is_collided(self.end):
raise ValueError(f"终点x坐标或y坐标在障碍物上")
# 算法初始化
self.reset(move_step, move_direction)
def reset(self, move_step=3, move_direction=8):
"""重置算法"""
self.__reset_flag = False
self.move_step = move_step # 移动步长(搜索后期会减小)
self.move_direction = move_direction # 移动方向 8 个
self.close_set = set() # 存储已经走过的位置及其G值
self.open_queue = NodeQueue() # 存储当前位置周围可行的位置及其F值
self.path_list = [] # 存储路径(CloseList里的数据无序)
def search(self):
"""搜索路径"""
return self.__call__()
def _in_map(self, node: Node):
"""点是否在网格地图中"""
return (0 <= node.x < self.width) and (0 <= node.y < self.high) # 右边不能取等!!!
def _is_collided(self, node: Node):
"""点是否和障碍物碰撞"""
return self.map_array[node.y, node.x] == 0
def _move(self):
"""移动点"""
@lru_cache(maxsize=3) # 避免参数相同时重复计算
def _move(move_step:int, move_direction:int):
move = (
[0, move_step], # 上
[0, -move_step], # 下
[-move_step, 0], # 左
[move_step, 0], # 右
[move_step, move_step], # 右上
[move_step, -move_step], # 右下
[-move_step, move_step], # 左上
[-move_step, -move_step], # 左下
)
return move[0:move_direction] # 坐标增量+代价
return _move(self.move_step, self.move_direction)
def _update_open_list(self, curr: Node):
"""open_list添加可行点"""
for add in self._move():
# 更新可行位置
next_ = curr + add
# 新位置是否在地图外边
if not self._in_map(next_):
continue
# 新位置是否碰到障碍物
if self._is_collided(next_):
continue
# 新位置是否在 CloseList 中
if next_ in self.close_set:
continue
# 计算所添加的结点的代价
H = next_ - self.end # 剩余距离估计
next_.cost = H # G = 0
# open-list添加/更新结点
self.open_queue.put(next_)
# 当剩余距离小时, 走慢一点
if H < 20:
self.move_step = 1
def __call__(self):
"""GBFS路径搜索"""
assert not self.__reset_flag, "call之前需要reset"
print("搜索中\n")
# 初始化列表
self.open_queue.put(self.start) # 初始化 OpenList
# 正向搜索节点(CloseList里的数据无序)
tic()
while not self.open_queue.empty():
# 弹出 OpenList 代价 H 最小的点
curr = self.open_queue.get()
# 更新 OpenList
self._update_open_list(curr)
# 更新 CloseList
self.close_set.add(curr) # G始终为0
# 结束迭代
if curr == self.end:
break
print("路径搜索完成\n")
toc()
# 节点组合成路径
while curr.parent is not None:
self.path_list.append(curr)
curr = curr.parent
self.path_list.reverse()
# 需要重置
self.__reset_flag = True
return self.path_list
# debug
if __name__ == '__main__':
p = GBFS()()
MAP.show_path(p)