CN109613547B - 一种基于反光板的占据栅格地图构建方法 - Google Patents

Abstract

本发明适用于自动控制技术领域，提供了一种基于反光板的占据栅格地图构建方法，步骤如下：S1、检测激光雷达在上一帧中是否定位成功，若检测结果为是，则基于动态匹配来获取当前帧中的已知反光板和未知反光板，否则基于静态匹配来获取当前帧中的已知反光板和未知反光板；S2、基于已知反光板来计算激光雷达在路标地图中的位姿，检测是否存在未知反光板，检测结果为是，计算未知反光板在路标地图中的坐标，执行步骤S3，检测结果为否，执行步骤S3；S3、更新占据栅格地图。反光板的高精度定位结果来降低栅格地图构建的累积误差，无需要其它传感器的辅助测量，在地图构建完毕之后，可以撤掉全部反光板，或在精确定位区域留下反光板，场景的干涉较低。

Description

一种基于反光板的占据栅格地图构建方法

技术领域

本发明属于自动控制技术领域，提供了一种基于反光板的占据栅格地图构建方法。

背景技术

随着社会的发展和技术的进步，移动机器人越来越深的介入到人类的日常生活中，例如家庭中的清洁机器人、工厂中的搬运机器人以及餐馆中的送餐机器人等。移动机器人想要实现上述功能必须准确的知道自身所在的位置，即实时定位。而移动机器人实时定位的先决条件是建立地图，这是机器人导航和其它智能行为的关键。激光雷达传感器因为测量精度高、不受光照条件影响等优点，被广泛的应用于移动机器人定位导航中。基于激光雷达的定位导航方法可以分为基于反光标识的定位导航方法和基于轮廓的定位导航方法，其中基于反光标识的激光雷达定位导航方法有着精度高、稳定性好等优点，但是有些场景可能会不方便布设反光标识，此外大量布设反光标识会增加成本。基于轮廓的定位导航方法不需要额外的布设反光标识，应用场景更广泛，物料成本更低。基于轮廓的定位导航方法首先要构建栅格地图，现有的基于激光雷达栅格地图构建方法累积误差大、需要里程计或IMU等其它传感器的辅助。

发明内容

本发明实施例提供一种基于反光板的占据栅格地图构建方法，基于反光板来构建精准的占据栅格地图，不依赖于里程计和IMU等辅助测量设备。

为了实现上述目的，本发明提供了一种基于反光板的占据栅格地图构建方法，所述方法包括如下步骤：

S1、检测激光雷达在上一帧中是否定位成功，若检测结果为是，则基于动态匹配来获取当前帧中的已知反光板和未知反光板，否则基于静态匹配来获取当前帧中的已知反光板和未知反光板；

S2、基于三个已知反光板来计算激光雷达在路标地图中的位姿，并检测是否存在未知反光板，若检测结果为是，则计算未知反光板在路标地图中的坐标，再执行步骤S3，若检测结果为否，则直接执行步骤S3；

S3、构建并更新栅格地图。

进一步的，在步骤S1之前还包括：

S4、检测路标地图中的路标数量是否零；

S5、若检测结果为是，则以激光雷达中心为路标地图的原点，计算当前帧中的反光板在路标地图中的坐标；若检测结果为否，则执行步骤S1。

进一步的，在步骤S4之前还包括：

S6、检测当前帧是否提取到了至少3个反光板；

S7、若检测结果为是，则执行步骤S4，若检测结果为否，则返回步骤S6。

进一步的，步骤S3具体包括如下步骤：

S31、计算光束起点及终点在栅格地图中的坐标；

S32、对光束终点对

所在的栅格进行占据更新，对从

到

所经过的栅格进行自由更新；

S33、计算占据更新栅格及自由更新栅格的栅格占据概率，并更新占据栅格地图。

本发明提供的基于反光板的占据栅格地图构建方法具有如下有益效果：反光板的高精度定位结果来降低栅格地图构建的累积误差，无需要其它传感器的辅助测量，在地图构建完毕之后，可以撤掉全部反光板，或在需要精确定位的区域留下反光板，对场景的干涉较低。

附图说明

图1为本发明实施例提供的基于反光板的占据栅格地图构建方法流程图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

本发明利用反光板的高精度定位结果来降低栅格地图构建的累积误差，不需要其它传感器的辅助测量，在地图构建完毕之后，可以撤掉全部反光板，或在需要精确定位的区域留下反光板，对场景的干涉较低。

图1为本发明实施例提供的基于反光板的占据栅格地图构建方法流程图，该方法包括如下步骤：

S1、检测激光雷达在上一帧中是否定位成功，若检测结果为是，则基于动态匹配来匹配出当前帧中的已知反光板和未知反光板；若检测结果为否，则基于静态匹配来匹配出当前帧中的已知反光板和未知反光板。

反光板即为路标的地图中的路标，已知反光板是指已经纳入路标地图中的反光板，未知反光板是指未纳入路标地图中的反光板。

在本发明实施例中，在步骤S1之前还包括检测当前帧提取到的反光板是否首次构建路标地图，其判断方法处理方法具体如下：

S4、检测路标地图中的路标数量是否零；

S5、若检测结果为是，则以激光雷达为路标地图的原点，计算当前帧中的反光板在路标地图中的位置；若检测结果为否，则执行步骤S1。

在本发明实施例中，在步骤S4之前还包括：

检测当前帧是否至少提取到了3个反光板，提取到的反光板数量低于3，无法实现激光雷达在当前帧的定位、未知反光板在路标地图中的定位，及占据栅格地图的构建，因此，在步骤S4之前还包括：

S6、检测当前帧是否提取到了至少3个反光板；

S7、若检测结果为是，检测反光板中心到激光雷达中心的距离及夹角，计算当前帧反光板在激光雷达坐标系中的坐标，并执行步骤S4，若检测结果为否，则返回步骤S6。

根据激光雷达探测距离与反光板反射强度之间的关系筛选出反光板上的扫描点，利用单个反光板上的反射点集拟合得到反光板中心，检测反光板中心距激光雷达中心的距离d_i及角度θ_i，用r_i(d_i,θ_i)来表示反光板r_i的中心距激光雷达中心的距离d_i及角度θ_i，基于r_i(d_i,θ_i)计算反光板中心在激光雷达坐标系中的坐标。

S2、基于三个已知反光板来计算激光雷达在路标地图中的位姿，并检测是否存在未知反光板，若检测结果为是，则计算未知反光板在路标地图中的坐标，并执行步骤S3，若检测结果为否，则执行步骤S3；

基于三个已知反光板来计算激光雷达在路标地图中的位姿，激光雷达在路标地图中的位姿计算方法具体如下：

已知反光板在雷达坐标系中的局部坐标X_loc、在路标地图坐标系中的坐标X，从反光板的局部坐标到反光板地图坐标之间转换满足刚体变换，即满足公式X＝RX_loc+t，基于三个已知反光板可以计算出雷达坐标系与地图坐标系之间的旋转矩阵R、平移向量t，因为激光雷达的局部坐标为

则激光雷达在地图中的坐标X^scan＝t。

计算未知反光板在路标地图中的坐标，计算方法具体如下：

在当前位置上连续定位50次，计算得到激光雷达的当前位姿相对于反光板地图坐标系的旋转平移矩阵的平均值

以及未知的反光板在激光雷达坐标系中的平均坐标

则未知反光板在路标地图坐标系中的坐标为

S3、更新占据栅格地图。

在本发明实施例中，需要将激光雷达水平放置，反光板布置于需要构建栅格地图的区域，使激光雷达时刻都能扫描到至少3个反光板。已知栅格地图尺寸为W*H，单个栅格的边长所代表的实际尺寸为l，设置栅格地图的原点为

占据栅格地图的更新方法具体如下：

S31、计算光束起点及终点在栅格地图中的坐标；

已知激光雷达当前帧在路标地图中的位姿P(x,y,θ)，则激光雷达在栅格地图中的位置为

激光雷达坐标系中的扫描点到栅格地图坐标系的转换矩阵为M，M的表达式具体如下：

其中P^m即是激光雷达的光束起点，已知第i个激光扫描点在激光雷达坐标系中的齐次坐标为p_i(x_i,y_i,1)，则第i个光束的终点在栅格地图中的坐标为

S32、对光束终点对

所在的栅格进行占据更新，对从

到

所经过的栅格进行自由更新；

已知一个光束的起点

和终点

k是光束经过的栅格的数量，可以利用优势对数积分法对光束经过栅格的概率进行更新。对

所在的栅格进行占据更新，对从

到

所经过的栅格进行自由更新。设置自由更新的值oddsFree＝0.4，占据更新的值oddsOcc＝0.6。

S33、计算占据更新栅格及自由更栅格的栅格占据概率，并更新栅格地图。

假定更新前的栅格占据概率为p，则更新前栅格的odds值为odds＝log(p/(1-p))，若当前栅格是进行自由更新，则令odds＝odds+oddsFree，若当前栅格是进行占据更新，令odds＝odds+oddsOcc，更新后，栅格的占据概率为p＝exp(odds)/(exp(odds)+1)。

本发明提供的基于反光板的占据栅格地图构建方法具有如下有益效果：反光板的高精度定位结果来降低栅格地图构建的累积误差，无需要其它传感器的辅助测量，在地图构建完毕之后，可以撤掉全部反光板，或在需要精确定位的区域留下反光板，对场景的干涉较低。

在本发明涉及的静态匹配及动态匹配可以采用现有的方法，或者是基于如下如下方法来实现静态匹配或动态匹配。

在本发明实施例中，设置反光板匹配允许的误差为matchW，定义全局最优匹配的误差best_err等于无穷大，本发明中的静态匹配具体包括如下步骤：

S11、根据路标地图中m个路标点R₁(X₁,Y₁),...,R_m(X_m,Y_m)，计算两两路标点之间的距离，得到全局距离列表D_1,2,...,D_1,m；...；D_m,1,...,D_m,m-1；

其中，D_m,m-1表示路标点R_m与路标点R_m-1之间的距离，(X_m,Y_m)为第m个路标点R_m在路标地图中的坐标；

S12、根据当前帧中的n个反光板r₁(x₁,y₁),...,r_n(x_n,y_n)，计算两两反光板之间的距离，得到局部距离列表d_1,2,....,d_1,n；....；d_n,1,...,d_n,n-1；

d_n,n-1表示反光板r_n与反光板r_n-1之间的距离，(x_n,y_n)为第n个反光板r_n在激光雷达坐标系中的坐标，根据激光雷达探测距离与反光板反射强度之间的关系筛选出反光板上的扫描点，利用单个反光板上的反射点集拟合得到反光板中心，检测反光板中心距激光雷达中心的距离d_i及角度θ_i，用r_i(d_i,θ_i)来表示反光板r_i的中心距激光雷达中心的距离d_i及夹角θ_i，基于r_i(d_i,θ_i)计算反光板中心在激光雷达坐标系中的坐标。

S13、选择全局距离列表中的第i行D_i与局部距离列表中的第j行d_j，假设路标地图中第i个路标点与当前帧中的第j个反光板匹配成功，将D_i行中的每一个元素与d_j行中的每一个元素进行对比，若D_i行中的第k个元素，与d_j行中的第l个元素满足abs(D_i,k-d_j,l)＜matchW，则将第k个路标点加入到第l个反光板的候选匹配中，对比结束得到D_i与d_j的匹配表matchList_i,j＝{r₁:R_i，...；...；r_n:R_j，...}；

S14、对匹配表matchList_i,j进行最佳匹配查找，得到当前最佳的匹配组合、匹配数目macth_count和匹配误差macth_err。若macth_count≥3且macth_err小于全局最优匹配的误差best_err，用当前匹配组合替换全局最优的匹配组合，令best_err＝macth_err；

S15、检测当前j的取值是否等于n，若j≠n，则令j＝j+1，并执行步骤S13，若j＝n，则检测当前i的取值是否为m，若i≠m，则令i＝i+1，j＝1，执行步骤S13，若i＝m，则输出全局最优匹配，若全局最优匹配的best_err＜matchW，则反光板匹配成功，返回全局最优匹配组合，若best_err≥matchW，则输出匹配失败。

在本发明实施例中，反光板的动态匹配具体包括如下步骤：

在激光雷达在上一帧中的位姿获取成功，在上一帧中的位姿为pose(x,y,θ)，基于动态匹配来寻找反光板对应的路标点，动态匹配的步骤如下：

S21、计算路标地图中所有路标点R₁(X₁,Y₁),....,R_m(X_m,Y_m)距激光雷达中心的距离，提取距离大于最小探测距离dist_min，且小于最大探测距离dist_max的路标点R₁(X₁,Y₁),....,R_z(X_z,Y_z)

S22、基于激光雷达在上一帧中的位姿pose(x,y,θ)，计算n个反光板r₁(x₁,y₁),...,r_n(x_n,y_n)在路标地图中的坐标P(x,y)，

S23、计算反光板r_k与路标点R_l的距离

若

将路标点R_l加入到反光板r_k的匹配列表中；

S24、检测当前l的取值是否等于z，若l≠z，则令l＝l+1，执行步骤S23，若l＝z，则检测当前k的取值是否为n，若k≠n，则令k＝k+1，l＝1，执行步骤S23，若k＝n，则输出匹配表matchList r₁:R_i,....；....；r_n:R_j,...；

S25、查找匹配表matchList中的最优匹配组合，若最优匹配组合的匹配数目大于2，且macth_err最小，若最优匹配组合的匹配误差小于matchW，则认为反光板匹配成功，返回最优匹配组合，否则，匹配失败。

在本发明实施例中，静态匹配及动态匹配中匹配表中的最佳匹配查找方法具体如下：

S31、按照设定的顺序从每个反光板的可能匹配路标点中选择一个，得到一个匹配组合

已知当前扫描到n个反光板，每个反光板所有的可能匹配为r₁:R_i,....；....；r_n:R_j,...，设置当前最佳匹配的误差为minErr。

S32、删除匹配组合

中错误匹配，获得匹配组合

计算

的匹配误差matchErr，若matchErr小于当前最佳匹配组合的匹配误差minErr，则将所述匹配组合

设为当前最佳匹配，令minErr＝matchErr，返回步骤S31，直至遍历匹配表中所有的匹配组合，输出最佳匹配、最佳匹配的匹配数目及匹配误差。

在本发明实施例中，错误匹配的删除方法具体包括如下步骤：

S321、从匹配组合

中顺序提取n个反光板及其对应的路标点；

S322、计算第i个反光板r_i与剩余反光板r_j(j≠i)间距离

以及路标点R_i与剩余路标点R_j(j≠i)间的距离

路标点R_i为反光板r_i的可能匹配路标点，路标点R_j为反光板r_i的可能匹配路标点；

S323、计算

与

之差deta_dist，统计deta_dist小于matchW的个数macth_count_i，若macth_count_i≥3，则匹配组合

的匹配数目macth_count加1，否则，将反光板r_i的匹配路标点设为空；

S324、检测i是否等于n，若i≠n，则令i＝i+1，并执行步骤S322，若i＝n，则输出匹配组合

即将

更新为

及匹配组合

的数目macth_count，若macth_count＜3，则返回步骤S321，若macth_count≥3，执行步骤S325；

S325、将匹配组合

中的反光板r_i与反光板r_j进行比较，检测反光板r_i与反光板r_j是否匹配到同一路标点；

S326、若反光板r_i与反光板r_j未匹配到同一路标点，则检测j是否等于n，若j≠n，令j＝j+1，执行步骤S325，若j＝n，则检测i是否等于n，若i≠n，令i＝i+1，j＝i+1执行步骤S325，若i＝n，则输出匹配组合

若反光板r_i与反光板r_j匹配到同一路标点，则匹配组合

的macth_count减去1，检测macth_count是否小于3，若macth_count＜3，则返回步骤S321，若macth_count≥3，则执行步骤S327；

S327、提取反光板r_i、r_j及两反光板r_i、r_j的剩余反光板r_m(m≠i,j)，剩余反光板中不包括匹配路标点为空的反光板；

S328、计算反光板r_i与剩余反光板r_m之间的距离

及路标点R_i与剩余路标点R_m间的距离

计算距离差

将距离差进行累积

计算反光板r_j与剩余反光板r_m之间的距离

计算路标点R_j与剩余路标点R_m间的距离

计算距离差

将距离差进行累积

其中，路标点R_i为反光板r_i在匹配组合

中的匹配路标点，路标点R_m为反光板r_m在该匹配组合

中的匹配路标点；

S329、比较sum_i与sum_j，若sum_i＜sum_j，则保留反光板r_i匹配的路标点，将反光板r_j的匹配路标点设为空，若sum_i＞sum_j，保留反光板r_j匹配的路标点，将反光板r_i的匹配路标点设为空；

S330、在匹配组合

中提取匹配路标点不为空的反光板r₁,...,r_s；

S331、计算反光板r_i与反光板r_j之间的距离

及路标点R_i与路标点R_j之间的距离

路标点R_i与路标点R_j分别为反光板r_i与反光板r_j在匹配组合

中的匹配路标点；

S332、计算

与

的差值，若

与

的差值大于matchW，则执行步骤S333；若

与

的差值小于或等于matchW，检测j是否等于s，若j≠s，令j＝j+1，执行步骤S331，若j＝s，则检测i是否等于s，若i≠s，则令i＝i+1，j＝i+1，执行步骤S331，若i＝s，则输出

S333、反光板r_i与反光板r_j之间至少有一个匹配错误，令匹配组合

的macth_count减1，若macth_count＜3，则返回步骤S321，若macth_count≥3，则执行步骤S334；

S334、提取反光板r_i、r_j及两反光板r_i、r_j的剩余反光板r_m(m≠i,j)，且剩余反光板中不包括匹配路标点为空的反光板；

S335、计算反光板r_i与剩余反光板r_m间的距离

及路标点R_i与剩余路标点R_m间的距离

计算距离差

将距离差进行累积

计算反光板r_j与剩余反光板r_m间的距离

及路标点R_j与剩余路标点R_m间的距离

计算距离差

将距离差进行累积

其中，路标点R_i为反光板r_i在匹配组合

中的匹配路标点，路标点R_m为反光板r_m在匹配组合

中的匹配路标点；

S336、将sum_i与sum_j进行比较，若sum_i＜sum_j，则保留反光板r_i匹配的路标点，将反光板r_j的匹配路标点设为空，sum_i≥sum_j，保留反光板r_j匹配的路标点，将反光板r_i的匹配路标点设为空。

该反光板的匹配方法具有如下有益效果：分层进行匹配，提高了匹配的速度，提高匹配的可靠性。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (3)

1.一种基于反光板的占据栅格地图构建方法，其特征在于，所述方法包括如下步骤：

S1、检测激光雷达在上一帧中是否定位成功，若检测结果为是，则基于动态匹配来获取当前帧中的已知反光板和未知反光板，若检测结果为否，基于静态匹配来获取当前帧中的已知反光板和未知反光板；

S2、基于三个已知反光板来计算激光雷达在路标地图中的位姿，并检测是否存在未知反光板，若检测结果为是，则计算未知反光板在路标地图中的坐标，再执行步骤S3，若检测结果为否，则直接执行步骤S3；

S3、构建并更新栅格地图，步骤S3具体包括如下步骤：

S31、计算光束起点及终点在栅格地图中的坐标；

S32、对光束终点

所在的栅格进行占据更新，对从

到

所经过的栅格进行自由更新；

S33、计算占据更新栅格及自由更新栅格的栅格占据概率，并更新占据栅格地图；

其中，

分别表示光束的起点和终点所在栅格，

表示光束终点所在栅格

的前一个栅格。

2.如权利要求1所述基于反光板的占据栅格地图构建方法，其特征在于，在步骤S1之前还包括：

S4、检测路标地图中的路标数量是否零；

S5、若检测结果为是，则以激光雷达中心为路标地图的原点，计算当前帧中的反光板在路标地图中的坐标；若检测结果为否，则执行步骤S1。

3.如权利要求1所述基于反光板的占据栅格地图构建方法，其特征在于，在步骤S4之前还包括：

S6、检测当前帧是否提取到了至少3个反光板；

S7、若检测结果为是，则执行步骤S4，若检测结果为否，则返回步骤S6。

Images