CN111765882A - Laser radar positioning method and related device thereof - Google Patents

Laser radar positioning method and related device thereof Download PDF

Info

Publication number
CN111765882A
CN111765882A CN202010561444.8A CN202010561444A CN111765882A CN 111765882 A CN111765882 A CN 111765882A CN 202010561444 A CN202010561444 A CN 202010561444A CN 111765882 A CN111765882 A CN 111765882A
Authority
CN
China
Prior art keywords
laser
points
matching
point set
point
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
Application number
CN202010561444.8A
Other languages
Chinese (zh)
Inventor
石鹏
林辉
卢维
殷俊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang Huaray Technology Co Ltd
Original Assignee
Zhejiang Dahua Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang Dahua Technology Co Ltd filed Critical Zhejiang Dahua Technology Co Ltd
Priority to CN202010561444.8A priority Critical patent/CN111765882A/en
Publication of CN111765882A publication Critical patent/CN111765882A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88Lidar systems specially adapted for specific applications
    • G01S17/89Lidar systems specially adapted for specific applications for mapping or imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/497Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

The application provides a laser radar positioning method and a related device thereof. The laser radar positioning method comprises the following steps: estimating the predicted pose of the laser radar; determining first coordinates of each laser point in the current frame laser point cloud data based on the predicted pose and the obtained current frame laser point cloud data; screening out laser points with poor matching degree based on the first coordinates to obtain a first laser point set; matching the first laser point set with the grid map based on the predicted pose to obtain a second coordinate of the laser point in the first laser point set; screening out laser points with poor matching degree based on the second coordinate to obtain a second laser point set; and matching the second laser point set with the map to determine the final pose of the laser radar. The method can remove the interference points so as to improve the positioning accuracy of the laser radar.

Description

Laser radar positioning method and related device thereof
Technical Field
The present application relates to the field of positioning technologies, and in particular, to a laser radar positioning method and a related apparatus.
Background
At present, among positioning methods based on laser radar, a positioning method based on a prior map is most commonly applied. The method has the main advantage that errors cannot be accumulated in the positioning process. However, a major drawback of such methods is: when the environment in which the map is constructed in advance changes, for example, a moving object and a local environment structure change (a case where some objects are stacked, etc.), the positioning performance is greatly affected.
Disclosure of Invention
The application provides a laser radar positioning method and a related device thereof, which can remove interference points so as to improve the positioning accuracy of the laser radar.
In order to achieve the above object, the present application provides a laser radar positioning method, including:
estimating the predicted pose of the laser radar;
determining first coordinates of each laser point in the current frame laser point cloud data based on the predicted pose and the obtained current frame laser point cloud data;
screening out laser points with poor matching degree based on the first coordinates to obtain a first laser point set;
matching the first laser point set with the grid map based on the predicted pose to obtain a second coordinate of the laser point in the first laser point set;
screening out laser points with poor matching degree based on the second coordinate to obtain a second laser point set;
and matching the second laser point set with the map to determine the final pose of the laser radar.
Screening out laser points with poor matching degree based on the second coordinate to obtain a second laser point set;
clustering laser points with similar distances in the first laser point set based on the second coordinate to obtain a plurality of classes;
screening out from each class a set of points of the laser spot that are greater than a first threshold from a second matching occupancy grid, wherein the second matching occupancy grid for each laser spot is the occupancy grid closest to each laser spot determined based on the second coordinates;
confirming the matching degree of the point concentration laser points and the second matching occupied grids, and screening out laser points with poor matching degree; and/or the presence of a gas in the gas,
and when the total number of the laser points of the point set is smaller than a second threshold value and the ratio of the total number of the laser points of the point set to the total number of the laser points in the class to which the point set belongs is smaller than a third threshold value, screening out all the laser points in the point set smaller than the third threshold value.
Wherein, confirm that the matching degree that laser point and second matching occupy the grid in the point set screens out the laser point that the matching degree is poor, include:
confirming the number of the laser points in the point set corresponding to the same second matching occupation grid, and screening out the laser points corresponding to the same second matching occupation grid when the number of the laser points in the point set exceeds a fourth threshold value and the distance between two laser points which are farthest away in the laser points corresponding to the same second matching occupation grid exceeds a fifth threshold value; and/or the presence of a gas in the gas,
taking the laser points and a plurality of adjacent laser points of the laser points in the point set in the class to which the point set belongs as a feature calculation point set, and calculating the geometric features of the geometric figure to which the laser points in the point set belong based on the feature calculation point set; and calculating a reference geometric characteristic based on the second matching occupation grids of all the laser points in the characteristic calculation point set, and screening out the laser points with low similarity degree based on the similarity degree of the geometric characteristic and the reference geometric characteristic.
The method for calculating the geometric characteristics of the geometric figure to which the laser points in the point set belong based on the characteristic calculation point set comprises the following steps: the normal vector and the curvature radius of the laser point are approximately solved by utilizing the eigenvalue and the eigenvector of the covariance of the distribution of the characteristic calculation point set;
computing a reference geometric feature based on a second matching occupancy grid of all laser points in the set of feature computation points, comprising: approximately solving a reference normal vector and a reference curvature radius of the laser point based on the eigenvalue and eigenvector of the covariance distributed by the second matching occupation of the center points of the grids in the feature calculation point set;
and screening out laser points with low similarity degree based on the similarity degree of the geometric features and the reference geometric features, wherein the screening out laser points comprises the following steps: when the included angle between the normal vector of the laser point and the reference direction quantity is larger than a sixth threshold value and smaller than a seventh threshold value, screening out the laser points which are larger than the sixth threshold value and smaller than the seventh threshold value;
and when the absolute value of the difference value between the curvature radius of the laser point and the reference curvature radius is larger than an eighth threshold value, screening out the laser points larger than the eighth threshold value.
Wherein screening out from each class a set of points of the laser spot larger than a first threshold from the first occupancy grid previously comprises: calculating a residual vector of the laser spot according to the second coordinate of the laser spot and the coordinate of the first occupied grid center point of the laser spot;
screening out from each class a set of points of the laser spot larger than a first threshold from the first occupancy grid, comprising: and screening out a point set of the laser points with the residual vector modular length larger than a first threshold value from each class, wherein the laser point residual vector modular length is calculated by the residual vector of the laser points.
And matching the second laser point set with a map to determine the final pose of the laser radar, wherein the method comprises the following steps: calculating a residual vector of the laser spot according to the second coordinate of the laser spot and the coordinate of the first occupied grid center point of the laser spot;
matching the second set of laser points to the map to determine a final pose of the lidar, previously comprising: and determining the weight coefficient of each laser point in the second laser point set, and determining the final pose of the laser radar according to the weight coefficient and the residual vector of each laser point in the second laser point set.
Wherein the weight coefficient of each laser spot in the second laser spot set is initially 1, and determining the weight coefficient of each laser spot in the second laser spot set includes:
calculating the direction angle of the residual vector of the laser point according to the residual vector of the laser point in the second laser point set;
clustering laser points with similar residual vector directions based on the residual vector direction angles to obtain a plurality of subsets;
and matching the subsets in the subset pairs with the opposite residual vector directions with the raster map one by one to confirm the matching scores of the subsets, and updating the weight coefficients of all the laser points in the subsets based on the matching scores.
Wherein, matching the subsets in the subset pairs with opposite residual vector directions with the grid map one by one, comprises:
confirming the sum of the residual vector moduli of all laser points in each subset of the subset pair with the opposite residual vector directions, and calculating the sum of the residual vector moduli of all the laser points in the second laser point set;
confirming whether the ratio of the sum of the modulus lengths of the residual vectors of all the subsets in the subset pair with the opposite residual vector direction to the sum of the modulus lengths of the residual vectors of the second laser point set is larger than a ninth threshold value;
and when the ratio is greater than a ninth threshold value and the ratio of the number of the laser points in all the subsets in the subset pairs with opposite residual vector directions to the total number of the laser points in the second laser point set exceeds a tenth threshold value, matching the subsets in the subset pairs with opposite residual vector directions with the grid map one by one.
Wherein, sieve out the laser point that the degree of matching is poor based on first coordinate, include:
taking the occupancy grid determined based on the first coordinates closest to each laser spot as a first matching occupancy grid for each laser spot;
laser spots that are greater than an eleventh threshold away from the first matching occupancy grid are screened out.
Wherein, estimating the predicted pose of the lidar comprises:
the predicted pose of the lidar is calculated based on other sensor information, or,
and estimating the linear velocity and the yaw angular velocity of the laser radar based on the final pose of the laser radar determined by the previous frames of laser point cloud data, and determining the predicted pose of the laser radar based on the linear velocity and the yaw angular velocity and the time interval between the current frame of laser point cloud data and the previous frame of laser point cloud data.
In order to achieve the above object, the present application provides a lidar positioning apparatus comprising a memory and a processor; the memory has stored therein a computer program for execution by the processor to perform the steps of the above method.
To achieve the above object, the present application provides a computer readable storage medium having a computer program stored thereon, which when executed by a processor, performs the steps of the above method.
The method comprises the following steps: firstly, determining a first coordinate of each laser point based on the predicted pose of the laser radar, determining the matching degree of the laser points and the grid map based on the first coordinate, and screening out the laser points with poor matching degree, so that the laser points with poor matching degree are screened out before matching, the matching precision can be improved, and the second coordinate of the laser points in the first laser point set is closer to a real coordinate; and after matching, laser points with poor matching degree are screened out based on the second coordinate, so that interference points can be further removed, and the positioning accuracy of the laser radar is improved.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.
FIG. 1 is a schematic flow chart diagram illustrating an embodiment of a lidar positioning method according to the present disclosure;
FIG. 2 is a schematic flowchart of an embodiment of a step of screening out laser spots with poor matching degrees based on a second coordinate in the laser radar positioning method of the present application;
FIG. 3 is a schematic structural diagram of an embodiment of a lidar positioning apparatus according to the present application;
FIG. 4 is a schematic structural diagram of an embodiment of a computer-readable storage medium according to the present application.
Detailed Description
In order to make those skilled in the art better understand the technical solution of the present application, the lidar positioning method and the related apparatus provided by the present application will be described in further detail below with reference to the accompanying drawings and the detailed description.
Referring to fig. 1, fig. 1 is a schematic flow chart of an embodiment of a laser radar positioning method according to the present application. The laser radar positioning method of the embodiment includes the following steps.
S110: and estimating the predicted pose of the laser radar.
Firstly, in the laser radar positioning process, the predicted pose of the laser radar can be predicted first, so that some laser points with poor matching degree can be screened out based on the predicted pose.
The pose information of the laser radar includes position information and yaw angle information of the laser radar. The position information of the laser radar may refer to two-dimensional position information or three-dimensional position information of the laser radar. For example, the pose information of the lidar may be used
Figure BDA0002546242470000051
And (4) showing. The predicted pose, the accurate pose and the initial pose of the laser radar can be used as the predicted pose of the laser radar.
In an implementation, the predicted pose of the lidar may be calculated based on other sensor information. Other sensors may be IMU, odometer, etc.
In another implementation, a predicted pose of the lidar is determined based on an accurate pose of the lidar determined from the previous frames of laser point cloud data. Specifically, in the case where the pose information includes position information and yaw angle information of the laser radar, the linear velocity and yaw angle velocity of the laser radar may be estimated from the accurate pose of the laser radar determined from the previous frames of laser point cloud data, and the predicted pose of the laser radar may be determined based on the linear velocity and yaw angle velocity, and the time interval between the current frame of laser point cloud data and the previous frame of laser point cloud data.
Illustratively, the first two frames of laser point cloud data S are utilizedk-1And Sk-2Linear velocity of determined accurate pose estimation current time k of laser radar
Figure BDA0002546242470000061
And yaw rate
Figure BDA0002546242470000062
Using estimated linear velocity
Figure BDA0002546242470000063
And yaw rate
Figure BDA0002546242470000064
And calculating the time interval delta t to predict the predicted pose of the laser radar at the current moment:
Figure BDA0002546242470000065
wherein the content of the first and second substances,
Figure BDA0002546242470000066
and
Figure BDA0002546242470000067
respectively, linear velocity and yaw rate, (x)k-1,yk-1,ψk-1) Is the accurate pose of the laser radar determined by the previous frame of laser point cloud data
Figure BDA0002546242470000068
And the predicted pose of the laser radar of the current frame laser point cloud data.
Certainly, in order to obtain a more accurate pose of the laser radar through the laser radar positioning method, a frame of laser point cloud data is used before estimating the predicted pose of the laser radar, and the initial pose of the laser radar is solved by using a CSM (correlation matching) method and a branch and bound method according to a priori initial value range.
S120: and determining first coordinates of each laser point in the current frame laser point cloud data based on the predicted pose and the obtained current frame laser point cloud data.
Calculating current frame excitation under the condition of predicting poseLight spot cloud data SkEach laser point p ini(i ═ 1,2, …, N) coordinates in the map coordinate system, which were taken as the first coordinates for each laser spot.
S130: and screening out laser points with poor matching degree based on the first coordinate to obtain a first laser point set.
And determining the matching degree of each laser point and the grid map based on the first coordinate of each laser point, screening out the laser points with poor matching degree, and reserving the laser points with good matching degree to obtain a first laser point set.
It will be appreciated that there will typically be an occupancy grid on the grid map corresponding to the laser points in the laser point cloud data, so the first matching occupancy grid for each laser point may be determined based on the first coordinates of each laser point. Alternatively, a grid corresponding to the laser point on the grid map may be determined based on the first coordinates of the laser point, and then an occupied grid closest to the grid corresponding to the laser point may be used as the first matching occupied grid of the laser point. And the grid where the coordinates of the laser points on the grid map are located is the grid corresponding to the laser points on the grid map. The first coordinates of the laser spot and the second coordinates of the laser spot may be used as the coordinates of the laser spot.
The degree of matching between the laser spot and its first matching occupied grid may refer to the degree of positional similarity between the two and the degree of shape similarity between the two, and the like. Accordingly, a laser spot having a low degree of similarity in occupying grid position and/or shape with the first matching can be used as a laser spot having a poor degree of matching.
For example, the degree of match between a laser spot and its first matched occupied grid may be expressed in terms of the distance between the laser spot and its first matched occupied grid. Such that a laser spot occupying the grid with the first matching at a distance exceeding the eleventh threshold may be regarded as a laser spot having a poor degree of matching. It will be appreciated that the distance between the two grids may be calculated from the coordinates of the two grids. The coordinates of the grid may be expressed in coordinates of a fixed point of the grid. The fixed point of the grid can be any point in the grid, such as the center point of the grid, the top left vertex of the grid, the top right vertex of the grid and the like.
Of course, the distance between the laser spot and its first matching occupied grid can also be determined in other ways. For example, a distance map corresponding to the grid map is established before determining that the laser spot occupies the grid with its first match. In the distance map, the distance corresponding to the occupied grid is 0, and the distance corresponding to the unoccupied grid is the distance between the occupied grid closest to the unoccupied grid and the unoccupied grid. And then when the distance between the grids occupied by the laser spot and the first matching of the laser spot is determined, firstly determining the grids corresponding to the laser spot, and determining the distance corresponding to the grids corresponding to the laser spot based on the distance mapping map.
S140: and matching the first laser point set with the grid map based on the predicted pose to obtain a second coordinate of the laser point in the first laser point set.
And matching the first laser point set with the grid map by taking the predicted pose as a matching initial value to obtain a second coordinate of the laser point in the first laser point set. The matching method may be ICP, NDT, or a correlation matching method, etc.
S150: and screening out the laser points with poor matching degree based on the second coordinate to obtain a second laser point set.
And determining the matching degree of each laser point in the first laser point set and the grid map based on the second coordinate of each laser point, screening out the laser points with poor matching degree, and reserving the laser points with good matching degree in the first laser point set to obtain a second laser point set.
Alternatively, a second matching occupied grid of each laser spot may be determined based on the second coordinates of each laser spot, then the matching degree between each laser spot in the first set of laser spots and the second matching occupied grid thereof is determined, and laser spots with poor matching degree in the first set of laser spots are screened out to obtain a second set of laser spots.
The grid corresponding to the laser point on the grid map may be determined based on the second coordinate of the laser point, and then the occupied grid closest to the grid corresponding to the laser point may be used as the second matching occupied grid of the laser point.
In addition, the degree of matching between each laser spot and its second matching occupied grid may refer to the degree of positional similarity between the two and/or the degree of shape similarity between the two, and the like. Accordingly, a laser spot having a low degree of similarity in position and/or shape to the second occupied grid can be used as a laser spot having a poor degree of matching.
Further, a portion of the laser spot occupying the grid with a lower degree of positional similarity or shape similarity with the second matching may be regarded as a laser spot having a poor degree of matching. For example, a part of the laser points occupying the position similarity of the grid with the second matching is taken as the laser points having a poor matching degree. Based on this, as shown in fig. 2, step S150 may include the following steps.
S151: and clustering the laser points with similar distances in the first laser point set based on the second coordinate to obtain a plurality of classes.
Specifically, the distance between two laser points adjacent to each other in the position is calculated by using the second coordinates of the two laser points adjacent to each other in the position, and if the distance between the two laser points adjacent to each other in the position is smaller than a twelfth threshold, the two laser points adjacent to each other in the position belong to the same class; and if the distance between two laser points adjacent to each other is larger than the twelfth threshold, the two laser points adjacent to each other do not belong to the same class. And judging whether the distance between the two laser points adjacent to each other in the position is smaller than a twelfth threshold value or not in sequence, and classifying the two laser points adjacent to each other in the position smaller than the threshold value into the same class to obtain a plurality of classes.
Each laser point has a corresponding serial number, and two laser points adjacent to the serial number are generally adjacent in position, so that whether the distance between two laser points adjacent to the serial number is smaller than a twelfth threshold value or not can be directly judged, whether the distance between two laser points adjacent to the position is smaller than the twelfth threshold value or not is judged, and the calculation process is simplified.
S152: a set of points from each class is screened from the laser spot whose distance to the second match occupies a grid greater than a first threshold.
The distance between each laser point and the second matching occupied grid of each class can be calculated, and then laser points which are larger than a first threshold value from the second matching occupied grid are screened out from each class, so as to obtain a point set of the laser points with the distance larger than the first threshold value in each class.
It will be appreciated that the step of calculating the distance between the respective laser spot and its second matching occupied grid in each class may be performed before step S151.
In one implementation, the distance between the grid corresponding to each laser point in each class and the second matching occupied grid of each laser point can be directly calculated.
In another implementation, the residual vectors of the laser spots may be calculated from the second coordinates of the laser spots and the coordinates of the second matching occupied grid of laser spots, and then the modulo lengths of the residual vectors of the respective laser spots may be calculated. And the mode length of the residual vector of each laser point is the distance between each laser point and the second matching occupied grid of each laser point. It will be appreciated that the corresponding grid of laser points may be determined using the second coordinates of the laser points, and then the residual vector of each laser point may be calculated using the coordinates of the grid corresponding to each laser point and the coordinates of the grid corresponding to each laser point occupied by the second matching of each laser point in each class. The specific calculation formula of the residual vector of each laser spot is as follows:
Figure BDA0002546242470000091
wherein the content of the first and second substances,
Figure BDA0002546242470000092
coordinates that refer to a second match of individual laser points occupying the grid;
Figure BDA0002546242470000093
the coordinates of grids corresponding to each laser point are referred to;
Figure BDA0002546242470000094
refers to the residual vector of each laser spot.
S153: and when the total number of the laser points of the point set is smaller than a second threshold value and the ratio of the total number of the laser points of the point set to the total number of the laser points in the class to which the point set belongs is smaller than a third threshold value, screening out all the laser points in the point set smaller than the third threshold value.
Discrete interference points in the class are eliminated by step S153.
S154: and confirming the matching degree of the laser points in the point set and the second matching occupying grid, and screening out the laser points with poor matching degree.
And determining the matching degree of each laser point occupying the grid by the second matching of each laser point in the point set of the laser points with the distance of each class larger than the first threshold value, and screening out the laser points with poor matching degree.
In one implementation, if the number of laser points in the point set corresponding to the same second matching occupancy grid is greater than a fourth threshold, it may be said that at least part of the laser points have an occupancy grid matching error. Wherein the fourth threshold value can be adjusted according to the laser scanning precision. Alternatively, in the case where the fourth threshold is exceeded, the multiple laser points in the set of points matching the same matching occupancy grid may all be deleted. Alternatively, in the case where the fourth threshold is exceeded, the laser spot whose occupancy grid matches erroneously may be further selected from all the laser spots in the spot set that match the same occupancy grid, and then the laser spot whose occupancy grid matches erroneously may be screened out. Still alternatively, in the case where the fourth threshold is exceeded, it may be further confirmed whether the distance of two laser dots farthest from among all the laser dots in the dot set matching the same matching occupancy grid exceeds a fifth threshold, and if the distance exceeds the fifth threshold, the laser dots corresponding to the same second matching occupancy grid are screened out. It will be appreciated that the maximum and minimum values may be taken from the sequence numbers of all laser points in the set of points where the same match occupies a grid match; and taking the distance between the laser point corresponding to the maximum value of the serial number and the laser point corresponding to the minimum value of the serial number as the distance between the two laser points which are farthest away from all the laser points in the point set matched with the same matching occupation grid.
In another implementation, the shape similarity between the respective laser spot and its second matching occupancy grid in each spot set can be further confirmed, and laser spots with poor shape similarity can be screened out. In particular, the similarity of the shape of the laser spot and its second matching occupied grid can be confirmed by comparing the geometrical features of the laser spot and its second matching occupied grid: if the difference between the geometrical characteristics of the laser point and the geometrical characteristics of the grid occupied by the second matching of the laser point is larger, the shape similarity of the grid occupied by the laser point and the second matching of the laser point is poor; otherwise, the shape similarity is good. Where the geometric feature may be a radius of curvature, a normal vector, or other value.
Specifically, the step of comparing the geometrical characteristics of the laser spot and the second matching occupied grid, confirming the matching degree of the laser spot and the second matching occupied grid, and screening out the laser spots with poor matching degree may include: taking the laser points and a plurality of adjacent laser points of the laser points in the point set in the class to which the point set belongs as a feature calculation point set, and calculating the geometric features of the geometric figure to which the laser points in the point set belong based on the feature calculation point set; and calculating a reference geometric characteristic based on the second matching occupation grids of all the laser points in the characteristic calculation point set, and screening out the laser points with low similarity degree based on the similarity degree of the geometric characteristic and the reference geometric characteristic. Since the sequence number adjacency basically means that the laser point positions are adjacent, a plurality of points of the laser point adjacent to the sequence number in the class to which the point set belongs can be regarded as a plurality of laser points of the laser point adjacent in the class to which the point set belongs. For example, a class includes a number of points with serial numbers 2, 5, 6, 7, 8, 12, 15, 19, 22, and 26, respectively, where 7 is a laser point whose distance is greater than a first threshold, and when calculating the geometric feature of 7, 5, 6, 7, 8, and 12 may be used as a feature calculation point set, the geometric feature of 7 may be calculated using the feature calculation point set, and the reference geometric feature may be calculated using the coordinates of the second matching occupied grid of 5, 6, 7, 8, and 12, respectively.
Further, when the geometric features of the shape similarity are judged to be normal vectors and curvature vectors, the calculating of the geometric features of the geometric figure to which the laser point belongs in the point set based on the feature calculation point set may include: and (4) approximately solving the normal vector and the curvature radius of the laser point by utilizing the eigenvalue and the eigenvector of the covariance of the distribution of the characteristic calculation point set. Wherein computing the reference geometric feature based on the second matching occupancy grid of all laser points in the set of feature computation points may include: and approximately solving a reference normal vector and a reference curvature radius of the laser point based on the eigenvalue and eigenvector of the covariance of the distribution of the central points of the second matching occupation grids of all the laser points in the feature calculation point set. When the included angle between the normal vector of the laser point and the reference normal vector is larger than a sixth threshold and smaller than a seventh threshold, the laser points larger than the sixth threshold and smaller than the seventh threshold are screened out; and/or when the absolute value of the difference value between the curvature radius of the laser point and the reference curvature radius is larger than an eighth threshold value, screening out the laser points larger than the eighth threshold value. The sixth threshold, the seventh threshold and the eighth threshold can be selected according to actual conditions. Preferably, the sixth threshold is close to 0 °, and may take the value 5 °. Preferably, the seventh threshold is close to 180 °, and may take 175 °.
The execution sequence of step S154 and step S153 is not limited. For example, step S154 may be performed before step S153, or simultaneously with step S153. In addition, in other embodiments, only step S154 may be performed, or only step S153 may be performed.
S160: and matching the second laser point set with the map to determine the accurate pose of the laser radar.
In an implementation mode, the second laser point set obtained through twice screening can be directly matched with a map for calculation so as to obtain the accurate pose of the laser radar.
In another implementation manner, the weight coefficient of each laser point can be determined based on the matching score of each laser point in the second laser point set and the grid map, then the second laser point set and the grid map are subjected to weighted matching based on the weight coefficient of each laser point, and the accurate pose of the laser radar is obtained through matching. And the residual vector corresponding to each laser point in the target function of the weighted matching process needs to be multiplied by the weight coefficient of the laser point. In addition, the method of weighted matching may include matching methods such as ICP, NDT, CSM, and the like.
Further, all laser points may be directly matched with the grid map to determine matching scores of the respective laser points with the grid map.
Or all laser points of the second laser point set can be divided into a plurality of subsets according to the residual vector, the subsets of the weight coefficients needing to be confirmed are matched with the grid map one by one to obtain matching scores of the subsets of the weight coefficients needing to be confirmed, and then the weight coefficients of all the laser points in the subsets are determined based on the matching scores of the subsets.
It can be understood that, in order to avoid the mutual influence of the two groups of laser points with opposite residual vector directions when the laser points are matched with the raster map, the two groups of laser points with opposite residual vector directions can be respectively classified into different subsets, and then the two groups of laser points with opposite residual vector directions are respectively matched with the map, so as to reduce the matching error and obtain an accurate matching score.
Preferably, all laser points of the second set of laser points may be divided into a plurality of subsets according to the direction of the residual vector, and the directions of the residual vectors of all laser points in each subset are the same or similar, so that two groups of laser points with opposite directions of the residual vectors necessarily belong to different subsets.
Illustratively, the method of dividing all laser points of the second set of laser points into a plurality of subsets may comprise: calculating the direction angle of the residual vector of the laser point according to the residual vector of the laser point in the second laser point set; and clustering the laser points with similar residual vector directions based on the residual vector direction angles to obtain a plurality of subsets. Wherein, the direction angle of the residual vector
Figure BDA0002546242470000121
The calculation formula of (a) is as follows:
Figure BDA0002546242470000122
wherein the content of the first and second substances,
Figure BDA0002546242470000123
is the residual vector of the laser spot;
Figure BDA0002546242470000124
is the residual vector direction angle of the laser spot.
It can be understood that the laser points with the residual vector direction angles in the same angle interval may be laser points with similar residual vector directions, and the residual vector direction angles of all the laser points in the second set of laser points are in the combination of all the angle intervals, and the different angle intervals do not overlap with each other. The difference between the upper limit value of the angle interval and the lower limit value of the angle interval may be an angle window threshold. In addition, the angle window thresholds for different angle intervals may be the same or different.
It should be noted that not all subsets need to determine the weighting coefficients based on the matching scores, that is, the weighting coefficients of all laser points in some subsets may be directly fixed values, for example, 1, although the value of the fixed value is not limited thereto.
Here, the subset having a large influence on the matching result may be determined as the weight coefficient based on the matching score, and the weight coefficients of all laser points in the subset having a small influence on the matching result may be set as fixed values. It will be appreciated that there are subsets where the number of laser spots occupies a larger grid distance from its second matching, and the number of laser spots is large, which may be a subset that has a large impact on the matching result. The sum of the modular lengths of all the residual vectors of laser spots in the subset can be compared with the sum of the modular lengths of all the residual vectors of laser spots in the second set of laser spots to determine whether there are laser spots in the subset that are further away from their second matching occupancy grid. Illustratively, if the ratio of the sum of the residual vector modulo lengths of the subset to the sum of the residual vector modulo lengths of the second set of laser points is greater than the ninth threshold, this indicates that there are several laser points in the subset that are further away from the second matching occupancy grid, and if the total number of laser points in the subset exceeds the tenth threshold, this subset is the subset that has a large influence on the matching result.
In other implementations, it may be determined whether it is necessary to determine the weight coefficients based on the matching scores for all subsets within a subset pair in units of determination of the subset pair. Wherein a subset pair comprises two subsets with opposite residual vector directions. For example, when both subsets of a subset pair are subsets that have a large influence on the matching result, all subsets of the subset pair need to be validated for weight coefficients based on the matching scores; when at least one of the two subsets of the subset pair is a subset having a small influence on the matching result, the weight coefficients of all laser points of all subsets within the subset pair are set to a fixed value without confirming the weight coefficients based on the matching scores for all subsets of the subset pair.
To better illustrate the steps of determining the weighting factors of all laser points in the second set of laser points in the positioning method of the present application, the following specific embodiments are provided for illustrative purposes:
counting the residual vectors of each laser point in the second laser point set, clustering the laser points with similar residual vector directions, and extracting a laser point set pair with opposite residual vector directions, wherein the method comprises the following steps:
1) calculating the direction angle of the matched residual vector:
traversing the second laser point set, and calculating the angle and the modular length of the residual vector direction of each laser point; wherein the direction angle of the vector pointing to the positive direction of the X axis of the map coordinate system is 0, and the anticlockwise direction is positive; matching residual vector direction angle of ith laser point
Figure BDA0002546242470000141
Figure BDA0002546242470000142
2) Clustering of matching residual vectors:
setting an angle window threshold β to be 5 degrees × i (i is not less than 1 and not more than 4), starting from-180 degrees, clustering in a β angle range (such as [ -180 degrees, β -180 degrees), [ β -180 degrees, 2 β -180 degrees °) …), clustering according to the matching residual vector direction angles of the laser points to obtain a plurality of subsets which are recorded as
Figure BDA0002546242470000143
3) Screening out laser point subset pairs which need to confirm the weight coefficient based on the matching score and have opposite residual vector directions:
wherein, in
Figure BDA0002546242470000144
When the temperature of the water is higher than the set temperature,
Figure BDA0002546242470000145
and
Figure BDA0002546242470000146
pairs of subsets of laser points with opposite residual vector directions.
Statistics of
Figure BDA0002546242470000147
And the sum of the modulo lengths of the residual vectors in (1) is recorded as
Figure BDA0002546242470000148
The sum of the modulus lengths of the residual vectors of all laser points is recorded as
Figure BDA0002546242470000149
If it is not
Figure BDA00025462424700001410
The laser point subset pair (marked as matching score-based) requiring the confirmation of the weight coefficient is satisfied
Figure BDA00025462424700001411
And
Figure BDA00025462424700001412
):
A. setting a proportional threshold
Figure BDA00025462424700001413
B. Laser spot set
Figure BDA00025462424700001414
Modulo sum of residual vector of
Figure BDA00025462424700001415
C. Setting a proportional threshold
Figure BDA00025462424700001416
Total number of laser spots of
Figure BDA00025462424700001417
Laser spot set
Figure BDA00025462424700001418
And
Figure BDA00025462424700001419
number of points (0)<α<1) All amounts are greater than
Figure BDA00025462424700001420
4) Calculating the weight coefficient of each laser point
A. Initializing the weight coefficients for each laser spot
Figure BDA00025462424700001421
Is 1.
B. Will need to confirm the weight coefficients based on the matching scores
Figure BDA00025462424700001422
Of all pairs of subsets of laser spots
Figure BDA00025462424700001423
And
Figure BDA00025462424700001424
respectively matching with map by CSM (correlation matching) method, and matching with lower-grade point set laser point weight coefficient
Figure BDA00025462424700001425
Set to α (0)<α<1) Wherein α is positively correlated with the match score.
In the embodiment, the first coordinate of each laser point is determined based on the predicted pose of the laser radar, the matching degree of the laser points and the grid map is determined based on the first coordinate, and the laser points with poor matching degree are screened out, so that the laser points with poor matching degree are screened out before matching, the matching precision can be improved, and the second coordinate of the laser points in the first laser point set is closer to the real coordinate; and after matching, laser points with poor matching degree are screened out based on the second coordinate, so that interference points can be further removed, and the accuracy of laser radar positioning is improved.
It can be understood that the laser radar positioning method can be applied to robot positioning, and because the laser radar is generally installed on the robot when the robot is positioned by the laser radar, the accurate pose of the laser radar determined based on the laser radar positioning method can be used as the accurate pose of the robot.
Referring to fig. 3, fig. 3 is a schematic structural diagram of an embodiment of a lidar positioning apparatus according to the present application. The laser radar positioning device 10 comprises a processor 12 and a memory 11; memory 11 is configured to store program instructions for implementing the laser radar location methods described above, and processor 12 is configured to execute the program instructions stored by memory 11.
The logic processes of the above-mentioned lidar positioning method are presented as a computer program, which may be stored in a computer-readable storage medium if it is sold or used as a stand-alone software product, and thus the present application proposes a computer-readable storage medium. Referring to fig. 4, fig. 4 is a schematic structural diagram of an embodiment of a computer-readable storage medium 20 of the present application, in which a computer program 21 is stored, and when the computer program is executed by a processor, the steps in the laser radar positioning method are implemented.
The computer-readable storage medium 20 may be a medium that can store a computer program, such as a usb disk, a removable hard disk, a Read-only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, or may be a server that stores the computer program, and the server can send the stored computer program to another device for running or can run the stored computer program by itself. The computer readable storage medium 20 may be a combination of a plurality of entities from a physical point of view, for example, a plurality of servers, a server plus a memory, or a memory plus a removable hard disk.
The above embodiments are merely examples and are not intended to limit the scope of the present disclosure, and all modifications, equivalents, and flow charts using the contents of the specification and drawings of the present disclosure or those directly or indirectly applied to other related technical fields are intended to be included in the scope of the present disclosure.

Claims (12)

1. A lidar positioning method, the method comprising:
estimating the predicted pose of the laser radar;
determining first coordinates of each laser point in the current frame laser point cloud data based on the predicted pose and the obtained current frame laser point cloud data;
screening out laser points with poor matching degree based on the first coordinates to obtain a first laser point set;
matching the first laser point set with a grid map based on the predicted pose to obtain a second coordinate of the laser point in the first laser point set;
screening out laser points with poor matching degree based on the second coordinate to obtain a second laser point set;
and matching the second laser point set with the map to determine the final pose of the laser radar.
2. The lidar positioning method according to claim 1, wherein the laser points with poor matching degree are screened out based on the second coordinate to obtain a second laser point set;
clustering laser points with similar distances in the first laser point set based on the second coordinate to obtain a plurality of classes;
screening out from each class a set of points of the laser spot that are greater than a first threshold from a second matching occupancy grid, wherein the second matching occupancy grid for each laser spot is the occupancy grid closest to each laser spot determined based on the second coordinates;
confirming the matching degree of the point concentration laser points and the second matching occupied grids, and screening out laser points with poor matching degree; and/or the presence of a gas in the gas,
and when the total number of the laser points of the point set is smaller than a second threshold value and the ratio of the total number of the laser points of the point set to the total number of the laser points in the class to which the point set belongs is smaller than a third threshold value, screening out all the laser points in the point set smaller than the third threshold value.
3. The lidar positioning method of claim 2, wherein determining a degree of matching between the laser point in the set of confirmation points and the second matching occupancy grid, and screening out laser points with poor matching comprises:
confirming the number of laser points in the point set corresponding to the same second matching occupation grid, and screening out the laser points corresponding to the same second matching occupation grid when the number exceeds a fourth threshold value and the distance between two laser points which are farthest away in the laser points corresponding to the same second matching occupation grid exceeds a fifth threshold value; and/or the presence of a gas in the gas,
taking the laser points and a plurality of adjacent laser points of the laser points in the point set in the class to which the point set belongs as a feature calculation point set, and calculating the geometric features of the geometric figure to which the laser points in the point set belong based on the feature calculation point set; and calculating a reference geometric characteristic based on the second matching occupation grids of all the laser points in the characteristic calculation point set, and screening out the laser points with low similarity degree based on the similarity degree of the geometric characteristic and the reference geometric characteristic.
4. The lidar positioning method according to claim 3, wherein the calculating the geometric feature of the geometric figure to which the laser point in the point set belongs based on the feature calculation point set comprises: the normal vector and the curvature radius of the laser point are approximately solved by utilizing the eigenvalue and the eigenvector of the covariance of the characteristic calculation point set distribution;
calculating a reference geometric feature based on the second matching occupancy grid of all laser points in the feature calculation point set, including: approximately solving a reference normal vector and a reference curvature radius of the laser point based on the eigenvalue and eigenvector of the covariance of the distribution of the central points of the second matching occupied grids of all the laser points in the feature calculation point set;
the screening out laser points with low similarity degree based on the similarity degree of the geometric features and the reference geometric features comprises the following steps: when the included angle between the normal vector of the laser point and the reference direction quantity is larger than a sixth threshold value and smaller than a seventh threshold value, screening out the laser points which are larger than the sixth threshold value and smaller than the seventh threshold value;
and when the absolute value of the difference value between the curvature radius of the laser point and the reference curvature radius is larger than an eighth threshold value, screening out the laser points larger than the eighth threshold value.
5. The lidar positioning method of claim 2, wherein the screening out from each class a set of points of the laser spot greater than a first threshold from a first occupancy grid previously comprises: calculating a residual vector of the laser spot according to the second coordinate of the laser spot and the coordinate of the first occupied grid center point of the laser spot;
said screening out from each class a set of points of the laser spot larger than a first threshold from the first occupancy grid, comprising: and screening out a point set of the laser points with the residual vector modular length larger than a first threshold value from each class, wherein the laser point residual vector modular length is calculated by the residual vector of the laser points.
6. The lidar positioning method of claim 1, wherein matching the second set of laser points to a map to determine a final pose of the lidar previously comprises: calculating a residual vector of the laser spot according to the second coordinate of the laser spot and the coordinate of the first occupied grid center point of the laser spot;
the matching of the second set of laser points with the map to determine the final pose of the lidar previously included: and determining a weight coefficient of each laser point in the second laser point set, and determining the final pose of the laser radar according to the weight coefficient and the residual vector of each laser point in the second laser point set.
7. The lidar positioning method of claim 6, wherein the weight coefficient of each laser spot in the second set of laser spots is initially 1, and wherein determining the weight coefficient of each laser spot in the second set of laser spots comprises:
calculating the direction angle of the residual vector of the laser point according to the residual vector of the laser point in the second laser point set;
clustering laser points with similar residual vector directions based on the residual vector direction angles to obtain a plurality of subsets;
and matching the subsets in the subset pairs with the opposite residual vector directions with the raster map one by one to confirm the matching scores of the subsets, and updating the weight coefficients of all the laser points in the subsets based on the matching scores.
8. The lidar positioning method of claim 7, wherein the step of matching the subsets of the subset pairs with opposite residual vector directions to the grid map one by one comprises:
confirming the sum of the residual vector moduli of all laser points in each subset of the subset pair with the opposite residual vector directions, and calculating the sum of the residual vector moduli of all the laser points in the second laser point set;
confirming whether the ratio of the sum of the modulus lengths of the residual vectors of all the subsets in the subset pair with the opposite residual vector direction to the sum of the modulus lengths of the residual vectors of the second laser point set is larger than a ninth threshold value;
and when the ratio is greater than a ninth threshold value and the ratio of the number of the laser points in all the subsets in the subset pairs with opposite residual vector directions to the total number of the laser points in the second laser point set exceeds a tenth threshold value, matching the subsets in the subset pairs with opposite residual vector directions with the grid map one by one.
9. The lidar positioning method according to claim 1, wherein the screening out laser points with poor matching degree based on the first coordinate comprises:
determining an occupancy grid closest to each laser spot based on the first coordinates as a first matching occupancy grid for each laser spot;
laser spots that are greater than an eleventh threshold away from the first matching occupancy grid are screened out.
10. The lidar positioning method of claim 1, wherein estimating the predicted pose of the lidar comprises:
the predicted pose of the lidar is calculated based on other sensor information, or,
and estimating the linear velocity and the yaw angular velocity of the laser radar based on the final pose of the laser radar determined by the previous frames of laser point cloud data, and determining the predicted pose of the laser radar based on the linear velocity and the yaw angular velocity and the time interval between the current frame of laser point cloud data and the previous frame of laser point cloud data.
11. A lidar positioning apparatus comprising a memory and a processor; the memory has stored therein a computer program for execution by the processor to carry out the steps of the method according to any one of claims 1 to 10.
12. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method according to any one of claims 1 to 10.
CN202010561444.8A 2020-06-18 2020-06-18 Laser radar positioning method and related device thereof Pending CN111765882A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010561444.8A CN111765882A (en) 2020-06-18 2020-06-18 Laser radar positioning method and related device thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010561444.8A CN111765882A (en) 2020-06-18 2020-06-18 Laser radar positioning method and related device thereof

Publications (1)

Publication Number Publication Date
CN111765882A true CN111765882A (en) 2020-10-13

Family

ID=72721008

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010561444.8A Pending CN111765882A (en) 2020-06-18 2020-06-18 Laser radar positioning method and related device thereof

Country Status (1)

Country Link
CN (1) CN111765882A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112363158A (en) * 2020-10-23 2021-02-12 浙江华睿科技有限公司 Pose estimation method for robot, and computer storage medium
CN112539756A (en) * 2020-11-30 2021-03-23 深圳市银星智能科技股份有限公司 Long channel identification method and robot
CN112612034A (en) * 2020-12-24 2021-04-06 哈尔滨工业大学芜湖机器人产业技术研究院 Pose matching method based on laser frame and probability map scanning
CN112710299A (en) * 2020-12-04 2021-04-27 深圳市优必选科技股份有限公司 Repositioning method, repositioning device, terminal equipment and storage medium
CN112987027A (en) * 2021-01-20 2021-06-18 长沙海格北斗信息技术有限公司 Positioning method of AMCL algorithm based on Gaussian model and storage medium
CN113607173A (en) * 2021-09-14 2021-11-05 成都睿芯行科技有限公司 Robot laser positioning method based on FPGA
WO2022267285A1 (en) * 2021-06-25 2022-12-29 深圳市优必选科技股份有限公司 Robot pose determination method and apparatus, robot and storage medium
CN117359644A (en) * 2023-12-04 2024-01-09 四川福莫斯工业技术有限公司 Point laser emission point pose calibration method and system based on robot

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106650640A (en) * 2016-12-05 2017-05-10 浙江大学 Negative obstacle detection method based on local structure feature of laser radar point cloud
CN107590827A (en) * 2017-09-15 2018-01-16 重庆邮电大学 A kind of indoor mobile robot vision SLAM methods based on Kinect
US20180306922A1 (en) * 2017-04-20 2018-10-25 Baidu Online Network Technology (Beijing) Co., Ltd Method and apparatus for positioning vehicle
US20190205695A1 (en) * 2017-12-29 2019-07-04 Baidu Online Network Technology (Beijing) Co., Ltd Method and apparatus for determining matching relationship between point cloud data
CN110018496A (en) * 2018-01-10 2019-07-16 北京京东尚科信息技术有限公司 Obstacle recognition method and device, electronic equipment, storage medium
CN110221603A (en) * 2019-05-13 2019-09-10 浙江大学 A kind of long-distance barrier object detecting method based on the fusion of laser radar multiframe point cloud
WO2019178429A1 (en) * 2018-03-15 2019-09-19 Kaarta, Inc. Methods and systems for filtering point cloud map data for use with acquired scan frame data
CN110297224A (en) * 2019-08-01 2019-10-01 深圳前海达闼云端智能科技有限公司 Laser radar positioning method and device, robot and computing equipment
CN110530368A (en) * 2019-08-22 2019-12-03 浙江大华技术股份有限公司 A kind of robot localization method and apparatus
CN110645998A (en) * 2019-09-10 2020-01-03 上海交通大学 Dynamic object-free map segmentation establishing method based on laser point cloud
CN110969649A (en) * 2019-11-29 2020-04-07 上海有个机器人有限公司 Matching evaluation method, medium, terminal and device of laser point cloud and map
CN111077539A (en) * 2019-12-30 2020-04-28 华通科技有限公司 Bird detection system based on laser radar
US20200174131A1 (en) * 2018-11-29 2020-06-04 Waymo Llc Localization Initialization For Autonomous Vehicles
CN111275763A (en) * 2020-01-20 2020-06-12 深圳市普渡科技有限公司 Closed loop detection system, multi-sensor fusion SLAM system and robot

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106650640A (en) * 2016-12-05 2017-05-10 浙江大学 Negative obstacle detection method based on local structure feature of laser radar point cloud
US20180306922A1 (en) * 2017-04-20 2018-10-25 Baidu Online Network Technology (Beijing) Co., Ltd Method and apparatus for positioning vehicle
CN107590827A (en) * 2017-09-15 2018-01-16 重庆邮电大学 A kind of indoor mobile robot vision SLAM methods based on Kinect
US20190205695A1 (en) * 2017-12-29 2019-07-04 Baidu Online Network Technology (Beijing) Co., Ltd Method and apparatus for determining matching relationship between point cloud data
CN110018496A (en) * 2018-01-10 2019-07-16 北京京东尚科信息技术有限公司 Obstacle recognition method and device, electronic equipment, storage medium
WO2019178429A1 (en) * 2018-03-15 2019-09-19 Kaarta, Inc. Methods and systems for filtering point cloud map data for use with acquired scan frame data
US20200174131A1 (en) * 2018-11-29 2020-06-04 Waymo Llc Localization Initialization For Autonomous Vehicles
CN110221603A (en) * 2019-05-13 2019-09-10 浙江大学 A kind of long-distance barrier object detecting method based on the fusion of laser radar multiframe point cloud
CN110297224A (en) * 2019-08-01 2019-10-01 深圳前海达闼云端智能科技有限公司 Laser radar positioning method and device, robot and computing equipment
CN110530368A (en) * 2019-08-22 2019-12-03 浙江大华技术股份有限公司 A kind of robot localization method and apparatus
CN110645998A (en) * 2019-09-10 2020-01-03 上海交通大学 Dynamic object-free map segmentation establishing method based on laser point cloud
CN110969649A (en) * 2019-11-29 2020-04-07 上海有个机器人有限公司 Matching evaluation method, medium, terminal and device of laser point cloud and map
CN111077539A (en) * 2019-12-30 2020-04-28 华通科技有限公司 Bird detection system based on laser radar
CN111275763A (en) * 2020-01-20 2020-06-12 深圳市普渡科技有限公司 Closed loop detection system, multi-sensor fusion SLAM system and robot

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
BASSIRI, A.; ASGHARI OSKOEI, M.; BASIRI, A.: "Particle Filter and Finite Impulse Response Filter Fusion and Hector SLAM to Improve the Performance of Robot Positioning", JOURNAL OF ROBOTICS, vol. 2018, pages 1 - 9 *
王任栋;徐友春;齐尧;韩栋斌;李华: "一种鲁棒的城市复杂动态场景点云配准方法", 机器人, vol. 40, no. 03, pages 257 - 265 *
陈逍遥 任小玲 夏 邢 史政坤: "基于标记的多状态离群点去除算法", 国外电子测量技术, vol. 39, no. 1, pages 39 - 43 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112363158B (en) * 2020-10-23 2024-03-12 浙江华睿科技股份有限公司 Pose estimation method for robot, robot and computer storage medium
CN112363158A (en) * 2020-10-23 2021-02-12 浙江华睿科技有限公司 Pose estimation method for robot, and computer storage medium
CN112539756A (en) * 2020-11-30 2021-03-23 深圳市银星智能科技股份有限公司 Long channel identification method and robot
CN112710299A (en) * 2020-12-04 2021-04-27 深圳市优必选科技股份有限公司 Repositioning method, repositioning device, terminal equipment and storage medium
CN112710299B (en) * 2020-12-04 2024-05-17 深圳市优必选科技股份有限公司 Repositioning method, repositioning device, terminal equipment and storage medium
CN112612034A (en) * 2020-12-24 2021-04-06 哈尔滨工业大学芜湖机器人产业技术研究院 Pose matching method based on laser frame and probability map scanning
CN112612034B (en) * 2020-12-24 2023-10-13 长三角哈特机器人产业技术研究院 Pose matching method based on laser frame and probability map scanning
CN112987027A (en) * 2021-01-20 2021-06-18 长沙海格北斗信息技术有限公司 Positioning method of AMCL algorithm based on Gaussian model and storage medium
CN112987027B (en) * 2021-01-20 2024-03-15 长沙海格北斗信息技术有限公司 Positioning method of AMCL algorithm based on Gaussian model and storage medium
WO2022267285A1 (en) * 2021-06-25 2022-12-29 深圳市优必选科技股份有限公司 Robot pose determination method and apparatus, robot and storage medium
CN113607173B (en) * 2021-09-14 2023-10-20 成都睿芯行科技有限公司 Robot laser positioning method based on FPGA
CN113607173A (en) * 2021-09-14 2021-11-05 成都睿芯行科技有限公司 Robot laser positioning method based on FPGA
CN117359644B (en) * 2023-12-04 2024-02-27 四川福莫斯工业技术有限公司 Point laser emission point pose calibration method and system based on robot
CN117359644A (en) * 2023-12-04 2024-01-09 四川福莫斯工业技术有限公司 Point laser emission point pose calibration method and system based on robot

Similar Documents

Publication Publication Date Title
CN111765882A (en) Laser radar positioning method and related device thereof
US10140529B2 (en) Method, apparatus and device for detecting lane lines
CN110472553B (en) Target tracking method, computing device and medium for fusion of image and laser point cloud
EP3518001B1 (en) Method for increasing the reliability of determining the position of a vehicle on the basis of a plurality of detection points
CN109932713B (en) Positioning method, positioning device, computer equipment, readable storage medium and robot
CN111444767B (en) Pedestrian detection and tracking method based on laser radar
CN112880694B (en) Method for determining the position of a vehicle
CN110349188B (en) Multi-target tracking method, device and storage medium based on TSK fuzzy model
CN111340875A (en) Space moving target detection method based on three-dimensional laser radar
CN112651274B (en) Road obstacle detection device, road obstacle detection method, and recording medium
CN112947415A (en) Indoor path planning method based on meaning information of barrier
CN111929653B (en) Target detection and tracking method and system based on unmanned ship navigation radar
JP6459981B2 (en) Image processing apparatus, subject identification method, and program
CN114612665B (en) Pose estimation and dynamic vehicle detection method based on normal vector histogram features
CN110736456A (en) Two-dimensional laser real-time positioning method based on feature extraction in sparse environment
Jumanov et al. Optimization of micro-object identification based on detection and correction of distorted image points
CN113281718B (en) 3D multi-target tracking system and method based on laser radar scene flow estimation
CN114241448A (en) Method and device for obtaining heading angle of obstacle, electronic equipment and vehicle
CN117291936A (en) Point cloud segmentation method, device, equipment and medium
CN111735443B (en) Dense target track correlation method based on assignment matrix
CN110749325B (en) Flight path planning method and device
CN112556698A (en) Ship behavior change and abnormality recognition method, equipment and medium
CN110275895B (en) Filling equipment, device and method for missing traffic data
CN113011376B (en) Marine ship remote sensing classification method and device, computer equipment and storage medium
CN116863325A (en) Method for multiple target detection and related product

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20201225

Address after: C10, 1199 Bin'an Road, Binjiang District, Hangzhou City, Zhejiang Province, 310051

Applicant after: ZHEJIANG HUARAY TECHNOLOGY Co.,Ltd.

Address before: No.1187 Bin'an Road, Binjiang District, Hangzhou City, Zhejiang Province

Applicant before: ZHEJIANG DAHUA TECHNOLOGY Co.,Ltd.

CB02 Change of applicant information
CB02 Change of applicant information

Address after: 310051 8 / F, building a, 1181 Bin'an Road, Binjiang District, Hangzhou City, Zhejiang Province

Applicant after: Zhejiang Huarui Technology Co.,Ltd.

Address before: C10, 1199 Bin'an Road, Binjiang District, Hangzhou City, Zhejiang Province, 310051

Applicant before: ZHEJIANG HUARAY TECHNOLOGY Co.,Ltd.