CN116295426A - UAV autonomous exploration method and device based on 3D reconstruction quality feedback - Google Patents
UAV autonomous exploration method and device based on 3D reconstruction quality feedback Download PDFInfo
- Publication number
- CN116295426A CN116295426A CN202310294785.7A CN202310294785A CN116295426A CN 116295426 A CN116295426 A CN 116295426A CN 202310294785 A CN202310294785 A CN 202310294785A CN 116295426 A CN116295426 A CN 116295426A
- Authority
- CN
- China
- Prior art keywords
- voxel
- point
- exploration
- candidate target
- gain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000011156 evaluation Methods 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000012216 screening Methods 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 4
- 230000007613 environmental effect Effects 0.000 description 3
- AYFVYJQAPQTCCC-GBXIJSLDSA-N L-threonine Chemical compound C[C@@H](O)[C@H](N)C(O)=O AYFVYJQAPQTCCC-GBXIJSLDSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000008447 perception Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/89—Lidar systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/88—Lidar systems specially adapted for specific applications
- G01S17/93—Lidar systems specially adapted for specific applications for anti-collision purposes
- G01S17/933—Lidar systems specially adapted for specific applications for anti-collision purposes of aircraft or spacecraft
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Electromagnetism (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
Description
技术领域technical field
本发明涉及测绘技术领域,特别是涉及基于三维重建质量反馈的无人机自主探索方法及装置。The invention relates to the technical field of surveying and mapping, in particular to an autonomous exploration method and device for unmanned aerial vehicles based on three-dimensional reconstruction quality feedback.
背景技术Background technique
大型工业建筑结构定期巡检等应用领域中,对三维建模准确性和系统自主性的要求越来越高。传统三维重建任务由操作人员手持测量设备执行,近年来有部分研究人员尝试利用无人载体来实现对小型物体或者小区域环境的自主探索和三维重建。然而,目前自主环境覆盖和探索规划领域中的大部分研究工作都侧重于,如何通过分析环境地图搜索出能以更小的时间和运动代价覆盖更多未知区域的最优运动轨迹,驱使载体沿着该路径运动的同时,收集传感器感知数据,最终获得完整覆盖环境的地图表达。该类探索算法有效性只根据空间中已知区域体积的比重度量,并未考虑到重建模型的不确定性。传统算法在环境特征非结构化、高度变化范围大的场景中,容易丢失局部区域的细节信息,降低环境模型置信度和完整性。为了补全局部细节,往往需要多次重复作业,导致采集的信息冗余、任务完成效率低下。In applications such as regular inspections of large industrial building structures, the requirements for 3D modeling accuracy and system autonomy are getting higher and higher. Traditional 3D reconstruction tasks are performed by operators with handheld measurement equipment. In recent years, some researchers have tried to use unmanned vehicles to realize autonomous exploration and 3D reconstruction of small objects or small area environments. However, most of the current research work in the field of autonomous environment coverage and exploration planning focuses on how to find the optimal trajectory that can cover more unknown areas with less time and movement costs by analyzing the environment map, so as to drive the carrier along the While moving along the path, sensor perception data is collected, and finally a map expression that completely covers the environment is obtained. The effectiveness of this type of exploration algorithm is only measured according to the specific gravity of the volume of the known region in space, and does not take into account the uncertainty of the reconstruction model. In scenes with unstructured environmental features and a large range of height changes, traditional algorithms tend to lose detailed information in local areas, reducing the confidence and integrity of the environmental model. In order to complete the local details, it is often necessary to repeat the operation many times, resulting in redundant collected information and low efficiency of task completion.
发明内容Contents of the invention
本发明的目的是针对传统大范围环境探索算法未考虑场景重建地图不确定性这一问题进行改进,在选择备选目标点时,不仅计算该点的未知信息增益,同时加入对表面重建效果的考虑,最终建立地下未知环境的高质量场景重建模型。The purpose of the present invention is to improve the problem that the traditional large-scale environment exploration algorithm does not consider the uncertainty of the scene reconstruction map. When selecting a candidate target point, it not only calculates the unknown information gain of the point, but also adds the effect of surface reconstruction. Consider, and finally build a high-quality scene reconstruction model of the underground unknown environment.
为实现上述目的,本发明提供了如下方案:To achieve the above object, the present invention provides the following scheme:
基于三维重建质量反馈的无人机自主探索方法,其特征在于,包括:The UAV autonomous exploration method based on 3D reconstruction quality feedback is characterized in that it includes:
采集激光雷达点云数据,基于所述激光雷达点云数据构建截断式有符号距离场地图;Gather lidar point cloud data, build a truncated signed distance field map based on the lidar point cloud data;
在探索空间内构建随机图,搜索候选目标点,获取各个所述候选目标点到无人机当前位置的最短路径;Construct a random graph in the exploration space, search for candidate target points, and obtain the shortest path from each candidate target point to the current position of the drone;
通过所述截断式有符号距离场地图判断体素状态,计算所述候选目标点的探索增益;Judging the voxel state by the truncated signed distance field map, and calculating the exploration gain of the candidate target point;
结合所述探索增益对所述最短路径进行评估,获取最佳路径。Evaluate the shortest path in combination with the exploration gain to obtain the best path.
可选地,根据所述激光雷达点云数据构建截断式有符号距离场地图包括:Optionally, constructing a truncated signed distance field map according to the lidar point cloud data includes:
将新点云中的每个激光反射点投影到全局地图的体素栅格中,合并所述激光反射点落入同一个体素中的所有激光射线,获得激光射线组;Projecting each laser reflection point in the new point cloud to the voxel grid of the global map, merging all laser rays whose laser reflection point falls into the same voxel, and obtaining a laser ray group;
在所述体素栅格处基于所述激光射线组执行光线追踪,遍历追踪光线末端点到传感器质心的每个体素,加权更新所述体素的投影距离和权重,构建所述截断式有符号距离场地图。Perform ray tracing based on the laser ray group at the voxel grid, traverse each voxel from the end point of the traced ray to the centroid of the sensor, update the projection distance and weight of the voxel by weighting, and construct the truncated signed Distance field map.
可选地,加权更新所述体素的投影距离和权重的方法为:Optionally, the weighted method for updating the projection distance and weight of the voxel is:
Wi+1(x,p)=min(Wi(x)+w(x,p),Wmax)W i+1 (x,p)=min(W i (x)+w(x,p),W max )
其中,x为该体素中心在导航系下的三维位置,p为新一帧点云中的某一个激光反射点在导航系下的三维位置,s为传感器在导航系下的三维位置,Di为地图中该体素当前存储的投影距离,Wi为存储的权重,Wmax为体素的权重上限,d为根据新一帧点云数据得到的沿着激光射线到物体表面的距离,w为根据新一帧点云数据得到的权重,Di+1为更新后的投影距离,Wi+1为更新后的权重,i为体素数据当前已更新的次数。Among them, x is the three-dimensional position of the voxel center in the navigation system, p is the three-dimensional position of a certain laser reflection point in the new point cloud in the navigation system, s is the three-dimensional position of the sensor in the navigation system, D i is the currently stored projection distance of the voxel in the map, W i is the stored weight, W max is the upper limit of the weight of the voxel, d is the distance from the laser ray to the surface of the object obtained according to the new frame of point cloud data, w is the weight obtained from the new frame of point cloud data, D i+1 is the updated projection distance, W i+1 is the updated weight, and i is the number of voxel data currently updated.
可选地,在所述探索空间内构建随机图包括:Optionally, constructing a random graph in the exploration space includes:
将所述无人机的当前位置作为根顶点加入到随机图中,在空间中通过Sample()函数采样生成随机点,并在所述随机图中搜索距离所述随机点最近的顶点;The current position of the drone is added to the random graph as the root vertex, a random point is generated by sampling the Sample() function in space, and the nearest vertex to the random point is searched in the random graph;
根据环境地图判断所述随机点指向所述顶点的连接是否可通行,若连线无碰撞,则将所述随机点作为新的顶点,即最新点,所述最新点指向所述随机点的连接作为新连接加入到所述随机图中;According to the environment map, it is judged whether the connection from the random point to the vertex is passable. If the connection has no collision, the random point is used as a new vertex, that is, the latest point, and the latest point points to the connection of the random point. join the random graph as a new connection;
所述随机图中所有位于所述最新点邻域内的顶点与所述最新点连接,若所述连线无碰撞,则将所述最新点邻域内的顶点与所述最新点连接加入到所述随机图中,获取最终的所述随机图。All vertices in the neighborhood of the latest point in the random graph are connected to the latest point, and if the connection has no collision, the vertices in the neighborhood of the latest point are connected to the latest point and added to the In the random graph, the final random graph is obtained.
可选地,搜索候选目标点,获取各个所述候选目标点到无人机当前位置的最短路径包括:Optionally, searching for candidate target points and obtaining the shortest path from each of the candidate target points to the current position of the drone includes:
随机图中的各顶点均为候选目标点,基于加入到所述随机图中各点的连接关系,通过Dijkstra算法搜索各个所述候选目标点到无人机当前位置的最短路径。Each vertex in the random graph is a candidate target point, and based on the connection relationship added to each point in the random graph, the Dijkstra algorithm is used to search for the shortest path from each of the candidate target points to the current position of the drone.
可选地,通过所述截断式有符号距离场地图判断体素状态包括:Optionally, judging the voxel state through the truncated signed distance field map includes:
所述体素状态分为未知、已知有障碍物和已知无障碍物;The voxel state is divided into unknown, known with obstacles and known without obstacles;
如果所述体素状态判断为未知,则判断以指定体素为中心的预设区域中周围预设个数的体素是否存在已知有障碍物体素。If the voxel state is determined to be unknown, it is determined whether there are known obstacle pixels in a preset number of voxels around the specified voxel in the preset area.
可选地,计算所述候选目标点的探索增益包括:Optionally, calculating the exploration gain of the candidate target point includes:
如果所述体素状态为未知且周围区域中不存在已知有障碍物体素,则累计该体素数量作为所述候选目标点的未知体素增益;If the voxel state is unknown and there is no known obstacle pixel in the surrounding area, then accumulating the number of voxels as the unknown voxel gain of the candidate target point;
如果所述体素状态为未知且周围区域中存在已知有障碍物,则所述体素为边界体素,计入所述候选目标点的重建表面增益,如果所述体素状态为已知有障碍物,当该体素的权重低于一定阈值时,同样计入所述候选目标点的重建表面增益;If the state of the voxel is unknown and there are known obstacles in the surrounding area, the voxel is a boundary voxel, which is included in the reconstruction surface gain of the candidate target point, if the state of the voxel is known There are obstacles, when the weight of the voxel is lower than a certain threshold, it is also included in the reconstruction surface gain of the candidate target point;
根据所述未知体素增益和所述重建表面增益进行加权计算,获取所述候选目标点的探索增益。Perform weighted calculation according to the unknown voxel gain and the reconstructed surface gain to obtain the exploration gain of the candidate target point.
可选地,获取所述最佳路径包括:Optionally, obtaining the best path includes:
根据路径长度、角度偏差和所述探索增益之和评估所有最短路径,获取路径探索增益最高的最短路径作为所述最佳路径。Evaluate all the shortest paths according to the sum of the path length, angle deviation and the exploration gain, and obtain the shortest path with the highest path exploration gain as the optimal path.
为实现上述目的,本发明提供了基于三维重建质量反馈的无人机自主探索装置,包括:In order to achieve the above object, the present invention provides a UAV autonomous exploration device based on 3D reconstruction quality feedback, including:
构图模块,用于根据激光雷达点云数据构建截断式有符号距离场地图;A composition module for constructing a truncated signed distance field map based on lidar point cloud data;
筛选点模块,用于在载体附近搜索候选目标点;Screening point module for searching candidate target points near the carrier;
路径优化模块,用于获取各个所述候选目标点到无人机的当前位置的的最短路径;A path optimization module, configured to obtain the shortest path from each of the candidate target points to the current position of the drone;
体素状态判断模块,用于通过所述截断式有符号距离场地图判断体素状态,计算所述候选目标点的探索增益;A voxel state judging module, configured to judge the voxel state through the truncated signed distance field map, and calculate the exploration gain of the candidate target point;
路径评估模块,用于计算所有最短路径的路径探索增益,获取最短路径和探索增益最高的最佳路径。The path evaluation module is used to calculate the path exploration gains of all the shortest paths, and obtain the shortest path and the best path with the highest exploration gain.
本发明的有益效果为:The beneficial effects of the present invention are:
本发明首先根据激光雷达点云数据构建截断式有符号距离场地图,然后,在探索空间内通过RRG算法搜索候选目标点。其次,通过Dijkstra算法搜索各个候选目标点到当前位置的最短路径,再基于地图数据判断体素状态,计算候选目标点的探索增益。最后,评估所有最短路径,执行路径探索增益最高的最佳路径。采用本发明的技术方案,无人机在执行自主探索任务时可以判断三维重建地图的不确定性,最终能提升场景重建的质量和置信度。The present invention first constructs a truncated signed distance field map according to the laser radar point cloud data, and then searches for candidate target points through the RRG algorithm in the exploration space. Secondly, the Dijkstra algorithm is used to search the shortest path from each candidate target point to the current location, and then the voxel state is judged based on the map data, and the exploration gain of the candidate target point is calculated. Finally, all shortest paths are evaluated, and path exploration is performed to find the best path with the highest gain. By adopting the technical solution of the present invention, the UAV can judge the uncertainty of the three-dimensional reconstruction map when performing autonomous exploration tasks, and finally can improve the quality and confidence of scene reconstruction.
附图说明Description of drawings
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the accompanying drawings required in the embodiments. Obviously, the accompanying drawings in the following description are only some of the present invention. Embodiments, for those of ordinary skill in the art, other drawings can also be obtained according to these drawings without paying creative labor.
图1为本发明实施例的基于三维重建质量反馈的无人机自主探索方法的流程图;Fig. 1 is the flow chart of the UAV autonomous exploration method based on the three-dimensional reconstruction quality feedback of the embodiment of the present invention;
图2为本发明实施例的基于三维重建质量反馈的无人机自主探索装置的结构图。FIG. 2 is a structural diagram of a UAV autonomous exploration device based on 3D reconstruction quality feedback according to an embodiment of the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some, not all, embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts belong to the protection scope of the present invention.
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。In order to make the above objects, features and advantages of the present invention more comprehensible, the present invention will be further described in detail below in conjunction with the accompanying drawings and specific embodiments.
基于三维重建质量反馈的无人机自主探索方法,包括:采集激光雷达点云数据,基于激光雷达点云数据构建截断式有符号距离场地图;在探索空间内构建随机图,搜索候选目标点,获取各个候选目标点到无人机当前位置的最短路径;通过截断式有符号距离场地图判断体素状态,计算候选目标点的探索增益;结合探索增益对最短路径进行评估,获取最佳路径。UAV autonomous exploration method based on 3D reconstruction quality feedback, including: collecting lidar point cloud data, constructing a truncated signed distance field map based on lidar point cloud data; constructing a random graph in the exploration space, searching for candidate target points, Obtain the shortest path from each candidate target point to the current position of the UAV; judge the voxel state through the truncated signed distance field map, and calculate the exploration gain of the candidate target point; combine the exploration gain to evaluate the shortest path to obtain the best path.
根据激光雷达点云数据构建截断式有符号距离场地图包括:将新点云中的每个激光反射点投影到全局地图的体素栅格中,合并激光反射点落入同一个体素中的所有激光射线,获得激光射线组;在体素栅格处基于激光射线组执行光线追踪,遍历追踪光线末端点到传感器质心的每个体素,加权更新体素的投影距离和权重,构建截断式有符号距离场地图。Constructing a truncated signed distance field map from lidar point cloud data includes: projecting each laser reflection point in the new point cloud into the voxel grid of the global map, and merging all laser reflection points falling into the same voxel Laser ray, get the laser ray group; perform ray tracing based on the laser ray group at the voxel grid, traverse each voxel from the end point of the traced ray to the sensor centroid, update the projection distance and weight of the voxel weightedly, and construct a truncated signed Distance field map.
在探索空间内构建随机图包括:将无人机的当前位置作为根顶点加入到随机图中,在空间中通过Sample()函数采样生成随机点,并在随机图中搜索距离随机点最近的顶点;根据环境地图判断随机点指向顶点的连接是否可通行,若连线无碰撞,则将随机点作为新的顶点,即最新点,最新点指向随机点的连接作为新连接加入到随机图中;随机图中所有位于最新点邻域内的顶点与最新点连接,若连线无碰撞,则将最新点邻域内的顶点与最新点连接加入到随机图中,获取最终的随机图。Constructing a random graph in the exploration space includes: adding the current position of the drone to the random graph as the root vertex, generating random points by sampling in the space through the Sample() function, and searching for the closest vertex to the random point in the random graph ;According to the environment map, it is judged whether the connection from the random point to the vertex is passable. If the connection has no collision, the random point is used as the new vertex, that is, the latest point, and the connection from the latest point to the random point is added to the random graph as a new connection; All the vertices in the neighborhood of the latest point in the random graph are connected to the latest point. If the connection has no collision, the vertices in the neighborhood of the latest point and the latest point are connected to the random graph to obtain the final random graph.
搜索候选目标点,获取各个候选目标点到无人机当前位置的最短路径包括:随机图中的各顶点均为候选目标点,基于加入到最终的随机图中各点的连接关系,通过Dijkstra算法搜索各个候选目标点到无人机当前位置的最短路径。Searching for candidate target points and obtaining the shortest path from each candidate target point to the current position of the UAV includes: each vertex in the random graph is a candidate target point, based on the connection relationship of each point added to the final random graph, through the Dijkstra algorithm Search the shortest path from each candidate target point to the current position of the UAV.
通过截断式有符号距离场地图判断体素状态包括:体素状态分为未知、已知有障碍物和已知无障碍物;如果体素状态判断为未知,则判断以指定体素为中心的预设区域中周围预设个数的体素是否存在已知有障碍物体素。Judging the voxel state through the truncated signed distance field map includes: the voxel state is divided into unknown, known obstacles and known obstacles; if the voxel state is judged as unknown, then judge the specified voxel as the center Whether there are known obstacle pixels in the preset number of surrounding voxels in the preset area.
计算候选目标点的探索增益包括:Calculating the exploration gain of candidate target points includes:
如果所述体素状态为未知且周围区域中不存在已知有障碍物体素,则累计该种体素数量作为所述候选目标点的未知体素增益;If the voxel state is unknown and there is no known obstacle pixel in the surrounding area, accumulating the number of voxels of this type as the unknown voxel gain of the candidate target point;
如果所述体素状态为未知且周围区域中存在已知有障碍物,则所述体素为边界体素,计入所述候选目标点的重建表面增益,如果所述体素状态为已知有障碍物,当该体素的权重低于一定阈值时,同样计入所述候选目标点的重建表面增益;If the state of the voxel is unknown and there are known obstacles in the surrounding area, the voxel is a boundary voxel, which is included in the reconstruction surface gain of the candidate target point, if the state of the voxel is known There are obstacles, when the weight of the voxel is lower than a certain threshold, it is also included in the reconstruction surface gain of the candidate target point;
根据未知体素增益和重建表面增益进行加权计算,获取候选目标点的探索增益。The weighted calculation is performed according to the unknown voxel gain and the reconstructed surface gain to obtain the exploration gain of the candidate target point.
获取最佳路径包括:根据路径长度、角度偏差和探索增益之和评估所有最短路径,获取路径探索增益最高的最短路径作为最佳路径。Obtaining the best path includes: evaluating all shortest paths according to the sum of path length, angle deviation and exploration gain, and obtaining the shortest path with the highest path exploration gain as the best path.
为实现上述目的,本实施例还提供了基于三维重建质量反馈的无人机自主探索装置,包括:构图模块,用于根据激光雷达点云数据构建截断式有符号距离场地图;筛选点模块,用于在载体附近搜索候选目标点;路径优化模块,用于获取各个候选目标点到无人机的当前位置的的最短路径;体素状态判断模块,用于通过截断式有符号距离场地图判断体素状态,计算候选目标点的探索增益;路径评估模块,用于计算所有最短路径的路径探索增益,获取最短路径和探索增益最高的最佳路径。In order to achieve the above purpose, this embodiment also provides an autonomous exploration device for UAVs based on 3D reconstruction quality feedback, including: a composition module, which is used to construct a truncated signed distance field map based on lidar point cloud data; a screening point module, It is used to search for candidate target points near the carrier; the path optimization module is used to obtain the shortest path from each candidate target point to the current position of the UAV; the voxel state judgment module is used to judge by truncated signed distance field map The voxel state calculates the exploration gain of the candidate target point; the path evaluation module is used to calculate the path exploration gain of all the shortest paths, and obtains the shortest path and the best path with the highest exploration gain.
如图1所示,本发明提供基于三维重建质量反馈的无人机自主探索方法,包括:As shown in Figure 1, the present invention provides a UAV autonomous exploration method based on 3D reconstruction quality feedback, including:
步骤1、点云数据是关于探索任务环境空间的数据,启动激光雷达,传感器通过网线向处理器传输点云数据即可采集此数据,根据激光雷达点云数据构建截断式有符号距离场地图(该地图由一个个正方体体素构成,每个体素中存储了截断式有符号距离和权重数据。有符号是指:如果体素在靠近传感器一侧,距离即为正,反之,则为负。截断式是指,如果距离太远或者太近,则将距离截断取为定值);Step 1. The point cloud data is the data about the environment space of the exploration mission. Start the lidar, and the sensor can collect the data by transmitting the point cloud data to the processor through the network cable. According to the lidar point cloud data, construct a truncated signed distance field map ( The map is composed of square voxels, and each voxel stores truncated signed distance and weight data. Signed means: if the voxel is on the side close to the sensor, the distance is positive, otherwise, it is negative. The truncated type means that if the distance is too far or too close, the distance is truncated as a fixed value);
步骤2、在探索任务区域(探索空间)内通过RRG算法搜索候选目标点;Step 2, search for candidate target points by RRG algorithm in the exploration task area (exploration space);
步骤3、通过Dijkstra算法搜索各个候选目标点到当前位置的最短路径;Step 3, search for the shortest path from each candidate target point to the current position through the Dijkstra algorithm;
步骤4、基于地图数据判断体素状态,计算候选目标点的探索增益;Step 4. Determine the voxel state based on the map data, and calculate the exploration gain of the candidate target point;
步骤5、评估所有最短路径,执行路径探索增益最高的最佳路径。Step 5. Evaluate all the shortest paths, and perform path search for the best path with the highest gain.
进一步,步骤1具体包括:Further, step 1 specifically includes:
步骤1.1、将新一帧点云中的每个激光反射点投影到全局地图的体素栅格中,再合并反射点落入同一个体素中的激光射线为一组激光射线;Step 1.1. Project each laser reflection point in a new frame of point cloud to the voxel grid of the global map, and then merge the laser rays whose reflection points fall into the same voxel into a set of laser rays;
步骤1.2、在该体素栅格处基于这一组激光射线只执行一次光线追踪,遍历追踪光线末端点到传感器质心的每个体素,加权更新体素的投影距离Di+1和权重Wi+1,Step 1.2: Perform ray tracing only once based on this group of laser rays at the voxel grid, traverse each voxel from the end point of the traced ray to the sensor centroid, and weight update the projection distance D i+1 and weight W i of the voxel +1 ,
Wi+1(x,p)=min(Wi(x)+w(x,p),Wmax)W i+1 (x,p)=min(W i (x)+w(x,p),W max )
其中,x为该体素中心在导航系下的三维位置,p为新一帧点云中的某一个激光反射点在导航系下的三维位置,s为传感器在导航系下的三维位置,Di+1为更新后的投影距离,Wi+1为更新后的权重,i为体素数据当前已更新的次数。Di为地图中该体素当前存储的投影距离,Wi为存储的权重,Wmax为体素的权重上限。d为根据新一帧点云数据得到的沿着激光射线到物体表面的距离,w为根据新一帧点云数据得到的权重,计算公式如下:Among them, x is the three-dimensional position of the voxel center in the navigation system, p is the three-dimensional position of a certain laser reflection point in the new point cloud in the navigation system, s is the three-dimensional position of the sensor in the navigation system, D i+1 is the updated projection distance, W i+1 is the updated weight, and i is the number of times the voxel data has been updated currently. D i is the currently stored projection distance of the voxel in the map, W i is the stored weight, and W max is the upper limit of the weight of the voxel. d is the distance along the laser ray to the surface of the object obtained according to the new frame of point cloud data, w is the weight obtained according to the new frame of point cloud data, and the calculation formula is as follows:
d(x,p,s)=‖p-x‖sign((p-x)·(p-s))d(x,p,s)=‖p-x‖sign((p-x)·(p-s))
其中,sign()表示符号函数,为向量的点乘运算符号,wt是投影距离大于截断距离且靠近传感器一侧的体素更新权重,0<wt<1,distrun为截断距离,size为地图体素大小。Among them, sign() represents the sign function, which is the point multiplication operation symbol of the vector, w t is the update weight of the voxel whose projection distance is greater than the truncation distance and close to the sensor side, 0<w t <1, dis trun is the truncation distance, size is the map voxel size.
进一步,步骤2具体包括:Further, step 2 specifically includes:
步骤2.1、将无人机当前位置作为根顶点加入到随机图中;Step 2.1, add the current position of the drone to the random graph as the root vertex;
步骤2.2、在探索空间中通过Sample()函数采样生成随机点Xrand;Step 2.2, generate a random point X rand by sampling the Sample() function in the exploration space;
步骤2.3、在随机图中搜索距离该随机点最近的顶点Xnearest;Step 2.3, searching for the vertex X nearest closest to the random point in the random graph;
步骤2.4、根据当前构建的环境地图,判断随机点指向最近点的连接是否可通行。如果连接是无碰撞的,那么就将该随机点作为新顶点,随机点指向最近点的连接和最近点指向随机点的连接作为新连接加入到随机图中(该连接为矢量,带有方向),随机图中的所有顶点均为候选目标点;Step 2.4. According to the currently constructed environment map, it is judged whether the connection from the random point to the nearest point is passable. If the connection is collision-free, then the random point is used as a new vertex, the connection from the random point to the nearest point and the connection from the nearest point to the random point are added to the random graph as a new connection (the connection is a vector with direction) , all vertices in the random graph are candidate target points;
步骤2.5、随机图中所有位于新顶点邻域内的顶点,也都形成与新顶点的连接,如果连接是安全的,就将这些连接也加入到随机图中。Step 2.5. All vertices in the random graph within the neighborhood of the new vertex also form connections with the new vertex. If the connections are safe, these connections are also added to the random graph.
进一步,步骤4具体包括:Further, step 4 specifically includes:
步骤4.1、基于构建的截断式有符号距离场地图判断体素m状态Status(m),状态分为未知、已知有障碍物和已知无障碍物三种:Step 4.1. Based on the constructed truncated signed distance field map, determine the Status (m) of the voxel m. The status is divided into three types: unknown, known with obstacles and known without obstacles:
其中,W(m)和D(m)分别为该体素的权重和投影距离,surthre为表面体素投影距离上限;Among them, W(m) and D(m) are the weight and projection distance of the voxel, respectively, and sur thre is the upper limit of the projection distance of the surface voxel;
步骤4.2、如果该体素的状态为未知,判断以某个体素为中心的正方体区域中的周围26个体素是否存在体素为已知有障碍物体素;Step 4.2, if the state of the voxel is unknown, determine whether there are voxels in the surrounding 26 voxels in the cube area centered on a certain voxel as known obstacle pixels;
步骤4.2.1、如果周围体素不存在已知有障碍物体素,则该体素需计入候选目标点ξ的未知体素增益Gainukn(ξ):Step 4.2.1. If there is no known obstacle pixel in the surrounding voxel, this voxel needs to be included in the unknown voxel gain Gain ukn (ξ) of the candidate target point ξ:
其中,Vis(ξ)是在候选点ξ处三维激光雷达视场范围内可以感知到的所有体素,Ukn为未知体素集合,Occup为已知有障碍物体素集合,Free为已知无障碍物体素集合;Among them, Vis(ξ) is all the voxels that can be perceived within the 3D lidar field of view at the candidate point ξ, Ukn is the unknown voxel set, Occup is the known obstacle pixel set, Free is the known barrier-free collection of object voxels;
步骤4.2.2、如果周围体素存在已知有障碍物体素,则该体素为边界体素,需计入候选目标点ξ的重建表面增益Gainrec(ξ);Step 4.2.2. If there are known obstacle pixels in the surrounding voxels, the voxel is a boundary voxel, which needs to be included in the reconstruction surface gain Gain rec (ξ) of the candidate target point ξ;
步骤4.3、如果该体素的状态为已知有障碍物,当置信度较低时,需同边界体素计入候选目标点ξ的重建表面增益Gainrec(ξ):Step 4.3. If the state of the voxel is known to have obstacles, when the confidence level is low, the reconstruction surface gain Gain rec (ξ) of the candidate target point ξ needs to be included in the same boundary voxel:
其中,U为边界体素集合,S为表面体素集合,ηthre为影响因子阈值,η(m)为体素影响因子,用于衡量表面体素对重建效果的影响大小,其计算公式为:Among them, U is the boundary voxel set, S is the surface voxel set, η thre is the threshold value of the influence factor, and η(m) is the voxel influence factor, which is used to measure the impact of surface voxels on the reconstruction effect, and its calculation formula is :
angspan=2tan(size,2z)ang span = 2tan(size,2z)
其中,z为传感器到体素中心的距离信息,Nrays(m)为能穿过该体素的激光射线数量,angspan为该体素在X和Y方向上射线可横跨的角度范围,angresx和angresy分别为三维激光雷达在X和Y方向上的角分辨率。Among them, z is the distance information from the sensor to the voxel center, N rays (m) is the number of laser rays that can pass through the voxel, and ang span is the angle range that the voxel can span in the X and Y directions, angres x and angres y are the angular resolutions of the 3D lidar in the X and Y directions, respectively.
步骤4.4、根据未知体素增益Gainukn(ξ)和重建表面增益Gainrec(ξ)两种增益,加权计算得到候选目标点的探索增益ExpGain(ξ)。Step 4.4. Based on the unknown voxel gain Gain ukn (ξ) and the reconstructed surface gain Gain rec (ξ), the weighted calculation obtains the exploration gain ExpGain (ξ) of the candidate target point.
进一步,步骤5具体包括:Further, step 5 specifically includes:
步骤5.1、根据路径长度、角度偏差和路径上各点的探索增益之和评估所有最短路径,第j条最短路径Pathj的路径探索增益计算公式如下:Step 5.1. Evaluate all the shortest paths according to the path length, angle deviation and the sum of the exploration gains of each point on the path. The formula for calculating the path exploration gain of the jth shortest path Path j is as follows:
其中,γS和γD为大于0的可调系数,Pathexp是一条长度与Pathj相同但是方向与无人机当前航向保持一致的路径,S(Pathj,Pathexp)是Pathj与Pathexp之间的距离矩阵。mj是该路径上的顶点总数,是路径Pathj上的第k个顶点,/>是沿着Pathj从/>到根顶点/>的距离。Among them, γ S and γ D are adjustable coefficients greater than 0, Path exp is a path whose length is the same as Path j but whose direction is consistent with the current heading of the UAV, S(Path j , Path exp ) is Path j and Path The distance matrix between exp . m j is the total number of vertices on the path, is the kth vertex on the path Path j , /> is along Path j from /> to the root vertex /> distance.
步骤5.2、路径探索增益最高的最短路径为最佳路径,并执行该最佳路径。Step 5.2: The shortest path with the highest path search gain is the best path, and the best path is executed.
本发明首先根据激光雷达点云数据构建截断式有符号距离场地图,然后,在探索空间内通过RRG算法搜索候选目标点。其次,通过Dijkstra算法搜索各个候选目标点到当前位置的最短路径,再基于地图数据判断体素状态,计算候选目标点的探索增益。最后,评估所有最短路径,执行路径探索增益最高的最佳路径。采用本发明的技术方案,无人机在执行自主探索任务时可以判断三维重建地图的不确定性,最终能提升场景重建的质量和置信度。The present invention first constructs a truncated signed distance field map according to the laser radar point cloud data, and then searches for candidate target points through the RRG algorithm in the exploration space. Secondly, the Dijkstra algorithm is used to search the shortest path from each candidate target point to the current location, and then the voxel state is judged based on the map data, and the exploration gain of the candidate target point is calculated. Finally, all shortest paths are evaluated, and path exploration is performed to find the best path with the highest gain. By adopting the technical solution of the present invention, the UAV can judge the uncertainty of the three-dimensional reconstruction map when performing autonomous exploration tasks, and finally can improve the quality and confidence of scene reconstruction.
如图2所示,本发明还提供基于三维重建质量反馈的无人机自主探索装置,包括:As shown in Figure 2, the present invention also provides a UAV autonomous exploration device based on 3D reconstruction quality feedback, including:
构图模块,用于根据激光雷达点云数据构建环境地图;A composition module for constructing an environmental map based on lidar point cloud data;
筛选点模块,用于在载体附近搜索候选目标点;Screening point module for searching candidate target points near the carrier;
路径优化模块,用于搜索各候选目标点到当前位置的最短路径;Path optimization module, for searching the shortest path from each candidate target point to the current position;
路径评估模块,用于计算所有最短路径的路径探索增益,最终执行路径探索增益最高的最佳路径。The path evaluation module is used to calculate the path exploration gains of all the shortest paths, and finally execute the best path with the highest path exploration gain.
采用本发明的技术方案,无人机在执行自主探索任务时可以判断三维重建地图的不确定性,最终能提升场景重建的质量和置信度。By adopting the technical solution of the present invention, the UAV can judge the uncertainty of the three-dimensional reconstruction map when performing autonomous exploration tasks, and finally can improve the quality and confidence of scene reconstruction.
以上所述的实施例仅是对本发明优选方式进行的描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。The above-mentioned embodiments are only descriptions of the preferred modes of the present invention, and do not limit the scope of the present invention. Variations and improvements should fall within the scope of protection defined by the claims of the present invention.
Claims (9)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310294785.7A CN116295426A (en) | 2023-03-24 | 2023-03-24 | UAV autonomous exploration method and device based on 3D reconstruction quality feedback |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310294785.7A CN116295426A (en) | 2023-03-24 | 2023-03-24 | UAV autonomous exploration method and device based on 3D reconstruction quality feedback |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116295426A true CN116295426A (en) | 2023-06-23 |
Family
ID=86835814
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310294785.7A Pending CN116295426A (en) | 2023-03-24 | 2023-03-24 | UAV autonomous exploration method and device based on 3D reconstruction quality feedback |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116295426A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117519150A (en) * | 2023-11-02 | 2024-02-06 | 浙江大学 | Autonomous implicit indoor scene reconstruction method combined with boundary exploration |
CN118334263A (en) * | 2024-06-11 | 2024-07-12 | 中国科学技术大学 | High-precision modeling method of fused laser point cloud based on truncated signed distance function |
CN118605575A (en) * | 2024-07-26 | 2024-09-06 | 中兴海陆工程有限公司 | A UAV autonomous exploration method based on boundary guidance |
-
2023
- 2023-03-24 CN CN202310294785.7A patent/CN116295426A/en active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117519150A (en) * | 2023-11-02 | 2024-02-06 | 浙江大学 | Autonomous implicit indoor scene reconstruction method combined with boundary exploration |
CN117519150B (en) * | 2023-11-02 | 2024-06-04 | 浙江大学 | Autonomous implicit indoor scene reconstruction method combined with boundary exploration |
CN118334263A (en) * | 2024-06-11 | 2024-07-12 | 中国科学技术大学 | High-precision modeling method of fused laser point cloud based on truncated signed distance function |
CN118605575A (en) * | 2024-07-26 | 2024-09-06 | 中兴海陆工程有限公司 | A UAV autonomous exploration method based on boundary guidance |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Khan et al. | A comparative survey of lidar-slam and lidar based sensor technologies | |
CN106802668B (en) | Unmanned aerial vehicle three-dimensional collision avoidance method and system based on binocular and ultrasonic fusion | |
Hrabar | 3D path planning and stereo-based obstacle avoidance for rotorcraft UAVs | |
CN116295426A (en) | UAV autonomous exploration method and device based on 3D reconstruction quality feedback | |
CN111982114B (en) | Rescue robot for estimating three-dimensional pose by adopting IMU data fusion | |
US20160004795A1 (en) | Device, method, apparatus, and computer-readable medium for solar site assessment | |
Ibrahim et al. | BIM-driven mission planning and navigation for automatic indoor construction progress detection using robotic ground platform | |
KR20180079428A (en) | Apparatus and method for automatic localization | |
CN107966989A (en) | A kind of robot autonomous navigation system | |
CN113867340B (en) | Beyond-the-earth unknown environment beyond-view-range global path planning system and method | |
JP2019504418A (en) | Method and system for determining the position of a moving object | |
CN207965645U (en) | A kind of robot autonomous navigation system | |
CN115855062A (en) | A Method for Autonomous Mapping and Path Planning of Indoor Mobile Robots | |
CN113741503A (en) | Autonomous positioning type unmanned aerial vehicle and indoor path autonomous planning method thereof | |
CN118857278A (en) | A method and system for surveying and mapping geographic information | |
CN111197986B (en) | Real-time early warning and obstacle avoidance method for three-dimensional path of unmanned aerial vehicle | |
Rekleitis et al. | Autonomous planetary exploration using LIDAR data | |
Singh et al. | High-resolution mapping of forested hills using real-time UAV terrain following | |
Bailey | Unmanned aerial vehicle path planning and image processing for orthoimagery and digital surface model generation | |
Gómez Arnaldo et al. | Path Planning for Unmanned Aerial Vehicles in Complex Environments | |
Kamat et al. | A survey on autonomous navigation techniques | |
Wettergreen et al. | Developing a framework for reliable autonomous surface mobility | |
CN118362112A (en) | Collaborative positioning mapping method based on multiple robots | |
Orlov et al. | Machine vision system for autonomous agricultural vehicle | |
Priyasad et al. | Point cloud based autonomous area exploration algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |