A星算法(A* A Star algorithm)原理以及代码实例,超详细,超简单,大白话谁都能看懂

news2024/11/23 20:29:23

本文以这篇博主的文章为基础【精选】A*算法(超级详细讲解,附有举例的详细手写步骤)-CSDN博客

这篇文章的博主做了一个UI界面,但我感觉,这样对新手关注算法和代码本身反而不利,会被界面的代码所干扰。所以笔者在他的基础上舍去了界面的内容,只关注算法本身。

A星算法的作用就是:已知起点和终点坐标,将地图离散化为网格,可以使用A star算法寻路。

A star算法简单来说三个步骤:一是准备两个列表,开列表和闭列表,开列表将节点移进移出,闭列表只进不出。二是在每一步探索的时候,计算三个“花费”,起点走到当前点的实际花费 G,从当前点到目标点的预估花费 H, 总花费F = G + H。三是计算父节点,父节点的作用是推算花费以及到达终点后倒推路径。

具体来说,我们称现在所处的网格为 checking point k检查checking point k周围所有的下一步可到达点并『临时』将这些可到达点的父节点记为checking point k。可到达点是没有障碍物且不在闭列表中的网格(near point 1, near point 2, ......, near point i, ......, near point n),对于near point i,计算起点到near point i的实际花费:

起点到near point i的实际花费 = 起点到checking point k的实际花费 + checking point k到near point i的实际花费。

计算near point i到终点的预估花费:

near point i到终点的预估花费 = near point i到终点的曼哈顿距离。

near point i的总花费 = 实际花费 + 预估花费。这里只是“花费“的一种定义形式,也可以用其他的定义形式。

每个点都是我们建立的点类的类实例,在点类中储存这三个花费,以及储存父节点

如果near point i不在开列表中,将其加进去。注意我们前面『临时』标记的父节点,这里正式标记checking point k为父节点

如果near point i已经在开列表中,则说明,near point i在我们到达checking point k之前,处于另一个checking point j时,作为checking point j的可到达点被计算过了。这里我们比较一下near point i旧的总花费和新的总花费,如果新的总花费小,则更新花费,并把near point i的父节点更新为checking point k;如果旧的总花费小,则保持旧的花费,以及保持旧的父节点。

当checking point k的所有可到达点都被检查完后,将其移进闭列表

之后,从开列表中选一个总花费最小的点作为新的checking point,重复上述的检查可达点操作,直到找到终点,起始时,将起点加入开列表,如果在途中开列表空了,则不存在可达路径。

完整代码如下:

import numpy as np
import json
import matplotlib.pyplot as plt


class Map:
    def __init__(self):
        # 设置起点,终点所在的行列数,左上角为0,0
        start_row = 17
        start_col = 2
        end_row = 9
        end_col = 16
        with open('map.txt', 'r') as f:
            my_map = json.loads(f.read())

        my_map = np.array(my_map)
        self.map = np.where(my_map == 0, 1, np.where(my_map == 1, 0, my_map))
        self.map[start_row, start_col] = -100   # 起点
        self.map[end_row, end_col] = 100        # 终点
        # self.map = np.array([[1, 1, 1, 1, 0, 100],
        #                      [1, 0, 0, 1, 0, 1],
        #                      [1, -100, 0, 1, 0, 1],
        #                      [1, 1, 0, 1, 0, 1],
        #                      [1, 1, 1, 1, 1, 1]])

    def get_start_point(self):
        indices = np.where(self.map == -100)
        return indices[0][0], indices[1][0]

    def get_end_point(self):
        indices = np.where(self.map == 100)
        return indices[0][0], indices[1][0]

    def check_grid(self, point):
        return self.map[point.x, point.y]


class Point:
    def __init__(self, x_, y_):
        self.x = x_
        self.y = y_
        self.father = None
        self.G = 0  # 起点到当前节点所花费的消耗
        self.H = 0  # 到终点的预估消耗
        self.F = 0

    def get_x_y(self):
        return self.x, self.y

    def set_GHF(self, G, H, F):
        self.H = H
        self.G = G
        self.F = F


class Astar:
    def __init__(self):
        self.openlist = []
        self.closelist = []
        self.map = Map()
        self.start_x, self.start_y = self.map.get_start_point()
        self.start_position = None
        self.end_x, self.end_y = self.map.get_end_point()
        self.find_path = False
        self.path = []

    def cal_GHF(self, checkpoint):
        if checkpoint.father is not None:
            G = checkpoint.father.G + 1   # 起点到父节点的花费加上父节点到本节点的花费
        else:
            G = 0
        H = abs(checkpoint.x - self.end_x) + abs(checkpoint.y - self.end_y)
        F = G + H
        return G, H, F

    def add_near_point(self, check_point):
        x, y = check_point.get_x_y()
        tmp_list = [Point(x-1, y-1), Point(x-1, y), Point(x-1, y+1),
                    Point(x, y-1), Point(x, y+1),
                    Point(x+1, y-1), Point(x+1, y), Point(x+1, y+1)]
        near_list = []
        for pi in tmp_list:
            if self.map.map.shape[0] > pi.x >= 0 and self.map.map.shape[1] > pi.y >= 0:     # 在地图范围内
                if self.map.check_grid(pi) == 100:
                    return [pi]
                elif self.map.check_grid(pi) == 1 and self.not_in_closelist(pi):
                    near_list.append(pi)

        return near_list

    def choose_min_F_point(self):
        minF = 1e10
        choosed_point = None
        for pi in self.openlist:
            if pi.F < minF:
                minF = pi.F
                choosed_point = pi
        return choosed_point

    def not_in_openlist(self, pi):
        not_in = True
        for pii in self.openlist:
            if pii.x == pi.x and pii.y == pi.y:
                not_in = False
        return not_in

    def not_in_closelist(self, pi):
        not_in = True
        for pii in self.closelist:
            if pii.x == pi.x and pii.y == pi.y:
                not_in = False
        return not_in

    def run(self):
        self.start_position = Point(self.start_x, self.start_y)
        G, H, F = self.cal_GHF(self.start_position)
        self.start_position.set_GHF(G, H, F)
        self.openlist.append(self.start_position)

        while True:
            checking_point = self.choose_min_F_point()
            if checking_point is None or self.find_path:
                print("End!")
                break
            self.openlist.remove(checking_point)
            self.closelist.append(checking_point)
            near_list = self.add_near_point(checking_point)
            for pi in near_list:
                if self.map.check_grid(pi) == 100:
                    self.find_path = True
                    # print("find path:\n{}".format(checking_point.get_x_y()))
                    self.path.append([checking_point.get_x_y()[0], checking_point.get_x_y()[1]])
                    reverse_point_father = checking_point.father
                    while reverse_point_father.father is not None:
                        # print(reverse_point_father.get_x_y())
                        self.path.append([reverse_point_father.get_x_y()[0], reverse_point_father.get_x_y()[1]])
                        reverse_point_father = reverse_point_father.father
                    break

                if self.not_in_openlist(pi):
                    pi.father = checking_point
                    G, H, F = self.cal_GHF(pi)
                    pi.set_GHF(G, H, F)
                    self.openlist.append(pi)
                else:
                    G_old = pi.G + 1
                    G_new = checking_point.G + 1
                    if G_new < G_old:
                        pi.father = checking_point
                        G, H, F = self.cal_GHF(pi)
                        pi.set_GHF(G, H, F)

        # 打印路性
        print("path: ")
        print(self.path)
        self.path = self.path[::-1]
        with open('my_path.txt', 'w') as file:
            for item in self.path:
                file.write(str(item) + '\n')

    def check(self):
        for pi in self.path:
            self.map.map[pi[0], pi[1]] = 2

        fig = plt.figure("A Star Algorithm")

        cmap = plt.cm.colors.ListedColormap(['yellow', 'black', 'white', 'blue', 'green'])

        # 创建一个离散的归一化器,根据不同数值映射到不同颜色
        bounds = [-101, -99, 0, 1, 2, 99, 101]
        norm = plt.cm.colors.BoundaryNorm(bounds, cmap.N)

        # 显示二维数组
        plt.imshow(self.map.map, cmap=cmap, norm=norm)

        # 添加颜色条,以便查看数值与颜色的对应关系
        cb = plt.colorbar()

        # 显示图
        plt.show()


if __name__ == "__main__":
    astar = Astar()
    astar.run()
    astar.check()

使用我的代码前,可以先运行文章开头提到的那个博主的代码,生成一个地图保存,我的代码加载他的地图(也可以使用我注释掉的那个地图),我的地图中,1表示可行网格,0表示障碍物网格。如果我们足够幸运的话,我们两篇文章的结果应该是一致的。

 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1144548.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

Win11 安装wsl遇到的问题解决

Win11 安装wsl遇到的问题解决 Win11 安装wsl遇到的问题解决WslRegisterDistribution failed:0x8007019eWslRegisterDistribution failed:0x800701bcUbuntu换源WSL通过网络访问Windows Win11 安装wsl遇到的问题解决 WslRegisterDistribution failed:0x8007019e 参考Link WslR…

关于 @Transactional 注解的类中使用 this 调用问题

在一个类中打断点的时候报了一个异常&#xff0c;这个异常&#xff0c;但是一直找不到&#xff0c;直到我在类中调用另外一个方法的时候&#xff0c;看到该方法里面用了 this 调用了 mybatis-plus 的 Api 去操作数据库&#xff0c;而最外层的方法却添加了 Transactional(rollba…

【python debug】python常见编译问题解决方法_1

1. ValueError: could not convert string to float: ‘File’ 问题原因&#xff1a;这里的报错是要把’File’这个转成float失败&#xff0c;非数字字符转换为float失败解决方法&#xff1a;然后去检查哪里出现了这个’File’&#xff0c;data load读入多个文件的数据&#x…

C++-类与对象(上)

一、 auto关键字 1.自动识别数据类型 2.auto的初始化 3.auto简化for循环 nullptr的使用 二、类与对象 1.c中类的定义 2.c语言与c的比较 3.类的访问限定符以及封装 3.1访问限定符 3.2封装 3.3类的作用域 3.4类的声明与定义分离 &#x1f5e1;CSDN主页&#xff1a;d1ff1cult.&…

转变命运!揭秘反转链表的神奇算法!

目录 使用虚拟头节点来辅助实现链表反转直接操作链表实现反转使用递归来实现链表反转 链表是计算机科学中常用的数据结构之一&#xff0c;它由一系列节点构成&#xff0c;每个节点包含一个值和指向下一个节点的指针。链表的灵活性使其在许多场景下被广泛应用&#xff0c;但其中…

S32K144芯片焊接完成后使用S32DS初次下载无法下载解决方法

一、问题现象如下&#xff0c;S32DS Debug下报错 二、原因&#xff0c;原厂芯片出厂时的FLASH Memory的安全机制是激活的&#xff0c;仿真器是可以连上&#xff0c;但是没法读取Flash Memory的内容 三、解决方法 参考图示&#xff0c;解锁后即可正常Debug

2023mathorcup大数据竞赛选题建议及思路

大家好呀&#xff0c;昨天6点2023年第四届MathorCup高校数学建模挑战赛——大数据竞赛开赛&#xff0c;在这里给大家带来初步的选题建议及思路。 注意&#xff0c;本文章只是比较简略的图文讲解&#xff0c;更加详细完整的视频讲解请移步&#xff1a; 2023mathorcup大数据数学…

【Opencv4快速入门】轮廓检测findContours

7.2 轮廓检测findContours 7.2.1 轮廓查找findContours7.2.2 轮廓绘制drawContours图像轮廓是指图像中对象的边界,是图像目标的外部特征,这个特征对于图像分析、目标识别和理解更深层次的含义具有重要的作用。 7.2.1 轮廓查找findContours 图像的轮廓补单能够提供物体的边缘,…

Crypto(6)攻防世界-babyrsa

参考文章&#xff1a; [攻防世界adworld] Crypto - babyrsa - 知乎 (zhihu.com) 涉及到的RSA知识点&#xff1a;

systrace/perfetto如何看surfaceflinger的vsync信号方法-android framework实战车载手机系统开发

背景&#xff1a; hi&#xff0c;粉丝朋友们&#xff1a; 大家好&#xff01;近期分享了surfaceflinger相关的一些blog&#xff0c;有同学就对相关的一些内容产生了一些疑问。 比如&#xff1a;vsync查看问题&#xff0c;即怎么才可以说是vsync到来了。 比如perfetto中surfac…

保姆级教学安装Linux操作系统,以及Linux的语法入门

&#x1f3c5;我是默&#xff0c;一个在CSDN分享笔记的博主。&#x1f4da;&#x1f4da; &#x1f31f;在这里&#xff0c;我要推荐给大家我的专栏《Linux》。&#x1f3af;&#x1f3af; &#x1f680;无论你是编程小白&#xff0c;还是有一定基础的程序员&#xff0c;这个专…

【性能测试】初识 Jmeter 中的 BeanShell

初识 Jmeter 中的 BeanShell 1.简介1.1 应用场景1.2 BeanShell 类型 2.常用内置变量2.1 log 日志模块2.2 vars 模块2.3 props 模块2.4 prev 模块 3.常见应用场景3.1 Java 文件处理3.2 导入外部 jar 包 BeanShell 是一个小型嵌入式 Java 源代码解释器&#xff0c;完全兼容 Java …

大数据采集技术与预处理学习一:大数据概念、数据预处理、网络数据采集

目录 大数据概念&#xff1a; 1.数据采集过程中会采集哪些类型的数据&#xff1f; 2.非结构化数据采集的特点是什么&#xff1f; 3.请阐述传统的数据采集与大数据采集的区别&#xff1f; ​​​​​​​ ​​​​​​​4.大数据采集的数据源有哪些&#xff1f;针对不同的数…

二叉树的前序、中序和后序非递归

目录 一、前序 二、中序 三、后序 一、前序 力扣&#xff08;LeetCode&#xff09;官网 - 全球极客挚爱的技术成长平台 /*** Definition for a binary tree node.* struct TreeNode {* int val;* TreeNode *left;* TreeNode *right;* TreeNode() : val(0)…

C++ 类和对象 (查漏补缺)

Inline 内联函数 内联函数是为了替代宏函数而出来的。 下面用宏实现一个ADD宏函数&#xff1a; 为什么这个ADD宏函数要这么写&#xff0c;首先我们来看&#xff0c;假设这样写&#xff1a; #define ADD(x,y)(xy) 会有什么问题呢&#xff1f; 宏函数是直接替换了&#xf…

FileInputStream文件字节输入流

一.概念 以内存为基准&#xff0c;把磁盘文件中的数据以字节形式读入内存中 二.构造器 public FileInputStream(File file) public FileInputStream(String pathname) 这两个都是创建字节输入流管道与源文件接通 三.方法 public int read() :每次读取一个字节返回&#xff0c;如…

设计模式:原型模式(C#、JAVA、JavaScript、C++、Python、Go、PHP)

上一篇《访问者模式》 下一篇《享元模式》 简介&#xff1a; 原型模式&#xff0c;它是一种创建型设计模式&#xff0c;它允许通过复制原型对象来创建新的对象&#xff0c;而无需知道创建的细节。其工作原…

安装Spring Tool Suite 4的一些问题

1、安装lombok.jar &#xff08;1&#xff09;把lombok.jar放到Spring Tool Suite 4的根目录下&#xff0c;java -jar lombok.jar&#xff0c;选择Spring Tool Suite 4&#xff0c;安装。然后SpringToolSuite4.ini里面会自动写入一行 sts4需要额外配置一行参数【-javaagent:lo…

springboot+mybatis-plus实现读写分离

shigen坚持日更的博客写手&#xff0c;擅长Java、python、vue、shell等编程语言和各种应用程序、脚本的开发。坚持记录和分享从业两年以来的技术积累和思考&#xff0c;不断沉淀和成长。 最近shigen加班也比较严重&#xff0c;很多天文章没有更新了&#xff0c;各位读者和伙伴见…