SLAM从入门到精通(bresenham绘制算法)

2023-10-15 13:12

本文主要是介绍SLAM从入门到精通(bresenham绘制算法),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

【 声明:版权所有,欢迎转载,请勿用于商业用途。 联系信箱:feixiaoxing @163.com】

        前面我们说过,学术界和工业界对于slam的要求是不一样的。前者要求robot在运动的过程中,同步实现定位和制图的操作。但是工业界中,一般两者是分开来的。首先进行制图,一旦制图完成之后就专注于定位操作。但是这两者之间没有明显的差别。

1、学术slam和工业slam的关系

        整理来说,完全可以把工业界的应用看成是slam的一个特例。真正的slam一般是这样的,首先在robot开始的位置设置为原点,lidar扫描的地方,绘制出一个地图。接着,通过遥控按钮操作robot,这个时候robot就走到了一个新的位置。借助于odom数据、lidar数据和之前的地图数据,我们可以求解出robot新位置的坐标。这样,我们通过robot的新位置、lidar扫描数据,又可以绘制出新的地图了。循环往复,不断重复上面的过程,最终就会可以得到一个完整的地图。这就是slam的基本原理。

2、slam和定位的转换

        有了上面的分析,大家也就明白了为什么很多论文都是讲slam算法,但是不谈具体的定位算法。这主要就是因为如果把slam算法中的制图删除、定位算法中的地图直接换成全局地图,这本身就是一个相对完整的定位算法方案了。

3、制图的算法

        有了准确的位姿之后,制图就不复杂了。前面我们说过,制图其实就是绘制很多的射线。雷达无法扫描到的地方一般被认为是灰色。没有达到最远测量距离,提前返回来的地方,那就是黑色,也就是有障碍物的地方。而lidar和障碍物之间的区域,这就是白色,可以被认为是能够自由活动的区域。

        上面这些描述,都是建立在lidar和base_link基本一致的基础上得出来的。如果两者不一致,那么还需要做一个tf转换。所以大家看到得是地图文件,实际上是通过一条、一条得射线绘制出来的。

4、bresenham算法

        这个算法本身是属于计算机图形学的一个基本算法。我们平时在绘制直线的时候,很容易遇到浮点数。但是计算机屏幕上面是没有浮点数的,每一个x、y坐标都是一个整数。那么怎么把这些直线上的数据在屏幕上面显示出来,这就是bresenham算法所要解决的事情。

        希望了解算法详细内容的同学,可以看下这个wiki连接,讲的还是很仔细的,

https://en.wikipedia.org/wiki/Bresenham%27s_line_algorithm

5、编译和测试

        不喜欢看算法,想直接看代码的同学也可以直接github上查找对应的代码。网上有很多bresenham算法的实现代码,这里给出一个链接,

https://gist.github.com/nowke/965fed0d5191bf373f1262be584207bb

        如果链接也不想看,这里给出一个可以编译版本的代码方案。主要是头文件需要修改下,

#include <GL/glut.h>
#include <stdio.h>int x1, y1, x2, y2;void myInit() {glClear(GL_COLOR_BUFFER_BIT);glClearColor(0.0, 0.0, 0.0, 1.0);glMatrixMode(GL_PROJECTION);gluOrtho2D(0, 500, 0, 500);
}void draw_pixel(int x, int y) {glBegin(GL_POINTS);glVertex2i(x, y);glEnd();
}void draw_line(int x1, int x2, int y1, int y2) {int dx, dy, i, e;int incx, incy, inc1, inc2;int x,y;dx = x2-x1;dy = y2-y1;if (dx < 0) dx = -dx;if (dy < 0) dy = -dy;incx = 1;if (x2 < x1) incx = -1;incy = 1;if (y2 < y1) incy = -1;x = x1; y = y1;if (dx > dy) {draw_pixel(x, y);e = 2 * dy-dx;inc1 = 2*(dy-dx);inc2 = 2*dy;for (i=0; i<dx; i++) {if (e >= 0) {y += incy;e += inc1;}elsee += inc2;x += incx;draw_pixel(x, y);}} else {draw_pixel(x, y);e = 2*dx-dy;inc1 = 2*(dx-dy);inc2 = 2*dx;for (i=0; i<dy; i++) {if (e >= 0) {x += incx;e += inc1;}elsee += inc2;y += incy;draw_pixel(x, y);}}
}void myDisplay() {draw_line(x1, x2, y1, y2);glFlush();
}void main(int argc, char **argv) {printf( "Enter (x1, y1, x2, y2)\n");scanf("%d %d %d %d", &x1, &y1, &x2, &y2);glutInit(&argc, argv);glutInitDisplayMode(GLUT_SINGLE|GLUT_RGB);glutInitWindowSize(500, 500);glutInitWindowPosition(0, 0);glutCreateWindow("Bresenham's Line Drawing");myInit();glutDisplayFunc(myDisplay);glutMainLoop();
}

        这里的编译命令是在ubuntu 20.04上实现的,如果有glut.h文件无法找到的错误,可以先安装一下这个库,

sudo apt-get install libglfw3-dev libgl1-mesa-dev libglu1-mesa-dev

        安装好库之后,就可以用命令gcc来编译了,

gcc draw.c -lglut -lGL -lGLEW -lGLU -o draw

        一切顺利的话,就可以生成./draw文件。执行后,会看到一个提示,要求我们输入起点和终点的坐标,

feixiaoxing@feixiaoxing-VirtualBox:~/Desktop$ ./draw
Enter (x1, y1, x2, y2)
100 3 200 500

        老老实实输入坐标后,按回车就可以看到对应的直线显示了,

        整个算法借助于opengl实现了显示过程。有兴趣的同学可以单步调试一下。

这篇关于SLAM从入门到精通(bresenham绘制算法)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/217949

相关文章

Python中的随机森林算法与实战

《Python中的随机森林算法与实战》本文详细介绍了随机森林算法,包括其原理、实现步骤、分类和回归案例,并讨论了其优点和缺点,通过面向对象编程实现了一个简单的随机森林模型,并应用于鸢尾花分类和波士顿房... 目录1、随机森林算法概述2、随机森林的原理3、实现步骤4、分类案例:使用随机森林预测鸢尾花品种4.1

使用Python绘制蛇年春节祝福艺术图

《使用Python绘制蛇年春节祝福艺术图》:本文主要介绍如何使用Python的Matplotlib库绘制一幅富有创意的“蛇年有福”艺术图,这幅图结合了数字,蛇形,花朵等装饰,需要的可以参考下... 目录1. 绘图的基本概念2. 准备工作3. 实现代码解析3.1 设置绘图画布3.2 绘制数字“2025”3.3

使用Python绘制可爱的招财猫

《使用Python绘制可爱的招财猫》招财猫,也被称为“幸运猫”,是一种象征财富和好运的吉祥物,经常出现在亚洲文化的商店、餐厅和家庭中,今天,我将带你用Python和matplotlib库从零开始绘制一... 目录1. 为什么选择用 python 绘制?2. 绘图的基本概念3. 实现代码解析3.1 设置绘图画

Python绘制土地利用和土地覆盖类型图示例详解

《Python绘制土地利用和土地覆盖类型图示例详解》本文介绍了如何使用Python绘制土地利用和土地覆盖类型图,并提供了详细的代码示例,通过安装所需的库,准备地理数据,使用geopandas和matp... 目录一、所需库的安装二、数据准备三、绘制土地利用和土地覆盖类型图四、代码解释五、其他可视化形式1.

如何用Python绘制简易动态圣诞树

《如何用Python绘制简易动态圣诞树》这篇文章主要给大家介绍了关于如何用Python绘制简易动态圣诞树,文中讲解了如何通过编写代码来实现特定的效果,包括代码的编写技巧和效果的展示,需要的朋友可以参考... 目录代码:效果:总结 代码:import randomimport timefrom math

Spring Security 从入门到进阶系列教程

Spring Security 入门系列 《保护 Web 应用的安全》 《Spring-Security-入门(一):登录与退出》 《Spring-Security-入门(二):基于数据库验证》 《Spring-Security-入门(三):密码加密》 《Spring-Security-入门(四):自定义-Filter》 《Spring-Security-入门(五):在 Sprin

不懂推荐算法也能设计推荐系统

本文以商业化应用推荐为例,告诉我们不懂推荐算法的产品,也能从产品侧出发, 设计出一款不错的推荐系统。 相信很多新手产品,看到算法二字,多是懵圈的。 什么排序算法、最短路径等都是相对传统的算法(注:传统是指科班出身的产品都会接触过)。但对于推荐算法,多数产品对着网上搜到的资源,都会无从下手。特别当某些推荐算法 和 “AI”扯上关系后,更是加大了理解的难度。 但,不了解推荐算法,就无法做推荐系

无人叉车3d激光slam多房间建图定位异常处理方案-墙体画线地图切分方案

墙体画线地图切分方案 针对问题:墙体两侧特征混淆误匹配,导致建图和定位偏差,表现为过门跳变、外月台走歪等 ·解决思路:预期的根治方案IGICP需要较长时间完成上线,先使用切分地图的工程化方案,即墙体两侧切分为不同地图,在某一侧只使用该侧地图进行定位 方案思路 切分原理:切分地图基于关键帧位置,而非点云。 理论基础:光照是直线的,一帧点云必定只能照射到墙的一侧,无法同时照到两侧实践考虑:关

康拓展开(hash算法中会用到)

康拓展开是一个全排列到一个自然数的双射(也就是某个全排列与某个自然数一一对应) 公式: X=a[n]*(n-1)!+a[n-1]*(n-2)!+...+a[i]*(i-1)!+...+a[1]*0! 其中,a[i]为整数,并且0<=a[i]<i,1<=i<=n。(a[i]在不同应用中的含义不同); 典型应用: 计算当前排列在所有由小到大全排列中的顺序,也就是说求当前排列是第

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个