CN113569393B - Robust absolute orientation algorithm - Google Patents
Robust absolute orientation algorithm Download PDFInfo
- Publication number
- CN113569393B CN113569393B CN202110789936.7A CN202110789936A CN113569393B CN 113569393 B CN113569393 B CN 113569393B CN 202110789936 A CN202110789936 A CN 202110789936A CN 113569393 B CN113569393 B CN 113569393B
- Authority
- CN
- China
- Prior art keywords
- translation vector
- weight
- rotation matrix
- algorithm
- calculating
- 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.)
- Active
Links
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 45
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 238000013519 translation Methods 0.000 claims abstract description 44
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000005457 optimization Methods 0.000 claims description 13
- 230000001351 cycling effect Effects 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 3
- 230000002829 reductive effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C1/00—Measuring angles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
-
- 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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/16—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using electromagnetic waves other than radio waves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Electromagnetism (AREA)
- Computer Hardware Design (AREA)
- Multimedia (AREA)
- Complex Calculations (AREA)
Abstract
The invention relates to a robust absolute orientation algorithm, which comprises the following steps of S1, constructing an absolute orientation problem solving model based on a weighting coefficient; step S2, setting the weight of each coordinate point asStep S3, deducing an analytic solution algorithm of an absolute directional problem solving model, and calculating a rotation matrix R and a translation vector t corresponding to the weight in the step S2; step S4, calculating the superposition error e i =||Rr' A,i ‑r' B,i After updating the weight by using the S estimation method, entering a step S3, and calculating a rotation matrix and a translation vector corresponding to the current weight; step S5, continuously looping steps S4 and S3, and outputting a rotation matrix R and a translation vector t when the iteration number of the S estimation method exceeds the set maximum iteration number. According to the invention, the coarse difference point can be automatically identified through the construction of the absolute directional problem solving model and the application of the updating weight of the S estimation method, and then through the analytic algorithm of the absolute directional problem solving model, the robustness of the algorithm is improved, and the solving precision of the algorithm is ensured.
Description
Technical Field
The invention relates to the field of photogrammetry, in particular to a robust absolute orientation algorithm.
Background
From the three-dimensional coordinates of a given n spatial points in two coordinate systems A, B, the problem of the rotation matrix R and translation vector t of coordinate system a relative to coordinate system B is solved, referred to as the absolute orientation problem. Absolute orientation is an important and fundamental problem in computer vision, and has wide application in the fields of point cloud matching, photogrammetry and the like.
The existing solving algorithm of the absolute orientation problem mainly comprises two types of an analysis algorithm and an iterative optimization algorithm. The analysis algorithm mainly comprises an absolute orientation algorithm based on an orthogonal matrix, an absolute orientation algorithm based on a quaternion and the like, and the main ideas of the algorithms are to solve and obtain an optimal solution of the absolute orientation problem through a matrix optimization technology. The iterative optimization algorithm comprises nonlinear optimization methods such as a beam adjustment algorithm, and the main idea of the algorithm is to use numerical optimization solving algorithms such as a least square method, and the like, and the residual error is reduced by continuous iteration until the residual error is smaller than a threshold value, so that an approximate optimal solution of the absolute orientation problem is obtained.
In solving the absolute orientation problem using actual data, measurement of three-dimensional space points often has observation noise and gross errors. For example, in point cloud registration, when three-dimensional point cloud data obtained by scanning in different directions are registered, rough differences exist in corresponding coordinates due to the reasons of limited measurement precision, different point cloud densities, partial point cloud data missing, characteristic point matching errors and the like. The traditional two resolving algorithms are solved under ideal conditions, the influence of coarse difference is not considered, and non-steady statistic such as arithmetic mean value and the like are used in the solving process, so that the robustness of the two resolving methods is not strong. The absolute orientation algorithm based on the MF estimator builds an objective function capable of inhibiting the influence of errors by modeling the errors and the coarse factors, so that the influence of the coarse is reduced to a certain extent, but the algorithm does not give an analytical solution, and the solving efficiency is limited.
Disclosure of Invention
Aiming at the problems, the invention provides a robust absolute orientation algorithm which can automatically identify rough difference points and improve algorithm robustness.
The technical scheme adopted for solving the technical problems is as follows: a robust absolute orientation algorithm comprising the steps of:
s1, constructing an absolute orientation problem solving model based on a weighting coefficient,
let n spatial points have three-dimensional coordinates r in coordinate systems A and B A,i And r B,i (i=1, 2, …, n) solving the model by a weighted method using an absolute directional problem in weighted form
Is used as an estimate of the rotation matrix and translation vector, where ω i More than or equal to 0 represents the weight corresponding to the ith coordinate point andSO (3) represents a third-order orthogonal matrix set with a determinant of 1, and I represents a 2-norm;
s2, setting the weight of each coordinate pointI.e. < ->
S3, deducing an analytic solution algorithm of an absolute directional problem solving model to obtain an optimal solution of a formula (2), calculating a rotation matrix R and a translation vector t corresponding to the weight value in S2,
the weighted centroid expression in the A, B coordinate system is recorded asAnd->(2) In the model of the formula, the solution of the rotation matrix and the translation vector is coupled, and new coordinates of an ith space point are introduced for decoupling>And->Ideal case->Wherein R is true And t true True values of rotation matrix and translation vector respectively, and then
r' B,i =R true r' A,i , (3)
And therefore, the optimization model is solved first
Obtaining a rotation matrix R corresponding to the current weight * Through again
Obtaining translation vector t in displacement parameter * Thereby realizing the aim of optimizing parameter decoupling,
due to
First two itemsIndependent of the rotation matrix, therefore, the optimization model of equation (4) is equivalent to
Wherein the method comprises the steps oftr () is the trace function of the matrix, and for solving the formula (7), singular value decomposition is performed on the matrix M to obtain +.>Wherein U and V are 3×3 orthogonal matrices, σ j 0 (j=1, 2, 3) is the singular value of matrix M, when
R * =V T U, (8)
(7) The maximum value is obtained and the maximum value is obtained,
to obtain R * Corresponding translation vector, deriving t from (5)The derivative is set to 0, and the translation direction corresponding to the weight value is obtainedThe amount is
S4, calculating the superposition error e i =||Rr' A,i -r' B,i After updating the weight by using the S estimation method, entering a step S3, calculating a rotation matrix corresponding to the current weight according to the formula (8), calculating a superposition error, and calculating a translation vector corresponding to the current weight according to the formula (9);
s5, continuously cycling the steps S4 and S3, and outputting a rotation matrix R and a translation vector t when the iteration number of the S estimation method exceeds the set maximum iteration number.
Preferably, the S estimation method in step S4 is specifically,
normalizing the error when the cycle number is 1Wherein->For the purpose of the scale estimation,represented |e i The median, 0.6754, is a constant in the S estimation, using
When the number of loops is greater than 1, the new weight is usedNormalization of the observation error>Using
As a new weight.
Compared with the prior art, the invention has the following beneficial effects:
1. according to the invention, the coarse difference points can be automatically identified through the construction of the absolute directional problem solving model and the application of the updating weight of the S estimation method, and then through the analytic algorithm of the absolute directional problem solving model, the influence caused by the coarse difference points can be effectively restrained, and under the condition of more coarse difference points, the solving precision can be ensured, and the robustness of the absolute directional algorithm is improved;
2. the invention provides an iterative weight calculation method based on an S estimation method, which can automatically identify rough difference points and automatically calculate proper weights of all coordinate points to realize the inhibition of the rough difference points.
Drawings
FIG. 1 is a graph showing the average solution error of a rotation matrix and a translation vector as a function of the percentage increase in the coarse difference;
wherein, the graph (a) corresponds to the rotation matrix error and the graph (b) corresponds to the translation vector error.
Detailed Description
The present invention will be described in detail with reference to fig. 1, wherein the exemplary embodiments and descriptions of the present invention are provided for explaining the present invention, but not limiting the present invention, and the absolute orientation algorithm of the present invention is the Proposed method algorithm.
A robust absolute orientation algorithm comprising the steps of:
s1, constructing an absolute orientation problem solving model based on a weighting coefficient,
let n spatial points have three-dimensional coordinates r in coordinate systems A and B A,i And r B,i (i=1, 2, …, n) solving the model by a weighted method using an absolute directional problem in weighted form
Is used as an estimate of the rotation matrix and translation vector, where ω i More than or equal to 0 represents the weight corresponding to the ith coordinate point andSO (3) represents a third-order orthogonal matrix set with determinant 1, and I represents 2 norms, if the weight is selected to satisfy the i coordinate point, omega is greater in the degree of the rough difference i The principle that the difference point is close to 0 is that the influence of the rough difference point on the optimization objective function of the formula (2) can be effectively reduced, so that the correct rotation matrix and translation vector can be obtained by solving the optimization model of the formula (2);
s2, setting the weight of each coordinate pointI.e. < ->
S3, deducing an analytic solution algorithm of an absolute directional problem solving model to obtain an optimal solution of a formula (2), calculating a rotation matrix R and a translation vector t corresponding to the weight value in S2,
the weighted centroid expression in the A, B coordinate system is recorded asAnd->(2) In the model of the formula, the solution of the rotation matrix and the translation vector is coupled, and new coordinates of an ith space point are introduced for decoupling>And->Ideal case->Wherein R is true And t true True values of rotation matrix and translation vector respectively, and then
r' B,i =R true r' A,i , (3)
And therefore, the optimization model is solved first
Obtaining a rotation matrix R corresponding to the current weight * Through again
Obtaining translation vector t in displacement parameter * Thereby realizing the aim of optimizing parameter decoupling,
due to
First two itemsIndependent of the rotation matrix, therefore, the optimization model of equation (4) is equivalent to
Wherein the method comprises the steps oftr () is the trace function of the matrix, and for solving the formula (7), singular value decomposition is performed on the matrix M to obtain +.>Wherein U and V are 3×3 orthogonal matrices, σ j 0 (j=1, 2, 3) is the singular value of matrix M, when
R * =V T U, (8)
(7) The maximum value is obtained and the maximum value is obtained,
to obtain R * Corresponding translation vector, deriving t from (5)Let the derivative be 0 to obtain a translation vector corresponding to the weight value as
S4, calculating the superposition error e i =||Rr' A,i -r' B,i After updating the weight by using the S estimation method, the method enters a step S3, calculates a rotation matrix corresponding to the current weight according to the formula (8), calculates the superposition error, calculates a translation vector corresponding to the current weight according to the formula (9),
normalizing the error when the cycle number is 1Wherein->For the purpose of the scale estimation,represented |e i The median, 0.6754, is a constant in the S estimation, using
When the number of loops is greater than 1, the new weight is usedNormalization of the observation error>Using
As a new weight;
s5, continuously cycling the steps S4 and S3, and outputting a rotation matrix R and a translation vector t when the iteration number of the S estimation method exceeds the set maximum iteration number.
In order to verify the effectiveness of the robust absolute orientation problem solving algorithm provided by the invention, the robustness is compared with a classical AbsPose algorithm based on an orthogonal matrix.
In the following simulation experiments, the spatial point coordinates and translation vector coordinates were in mm. The rotation matrix in the simulation test is a randomly generated SO (3) matrix, and three components of the translation vector t are in [0,1000 ]]The random number selected in the group. The coordinates of 20 ideal space points in the A coordinate system are uniformly distributed in [ -30,30]×[-30,30]×[-30,30]And (3) rotating and translating the ideal space point coordinate in the A coordinate system to obtain the ideal space point coordinate in the B coordinate system. To test robustness, the ideal spatial point coordinates in the A, B coordinate system are added to a mean value of 0, gaussian noise with variance of 0.2, then m (m=1,..10.) points were randomly selected in the B coordinate system, and 0,50 were added, respectively]Is used to simulate the coarse difference. The error of the rotation matrix is defined asWherein r is ij,true And r ij Elements of the ith row and jth column representing the rotation matrix true and resolved values, respectively, the translation vector error is defined as +.>Wherein t is true And t represent translation vector true and resolved values, respectively. For each m, independent experiments were performed 1000 times.
Fig. 1 shows the variation of the average solution error of the rotation matrix and the translation vector with the increase of the percentage of the coarse difference point, wherein the graph (a) corresponds to the rotation matrix error, the graph (b) corresponds to the translation vector error, the dots are the solution result of the absPose algorithm, and the square dots are the solution result of the algorithm provided by the invention. As can be seen from the graph, as the number of coarse difference points increases, both the rotation matrix error and the translation vector error of the abswise algorithm show a growing trend, and the growing speed is far greater than that of the robust absolute orientation problem solving algorithm provided by the invention. Under the condition that the percentages of the rough difference points are the same, the algorithm provided by the invention can greatly improve the measurement accuracy of the rotation matrix and the translation vector. When the rough difference point percentage is 50%, the average translation vector relative error obtained by the AbsPose algorithm is as high as 74.0%, and the translation vector relative error of the algorithm is only 2.3%, so that the solution accuracy of one order of magnitude is improved compared with the AbsPose. Therefore, the algorithm provided by the invention can effectively inhibit the influence caused by the rough difference points, can still ensure the solving precision under the condition of more rough difference points, and has good robustness.
The foregoing has described in detail the technical solutions provided by the embodiments of the present invention, and specific examples have been applied to illustrate the principles and implementations of the embodiments of the present invention, where the above description of the embodiments is only suitable for helping to understand the principles of the embodiments of the present invention; meanwhile, as for those skilled in the art, according to the embodiments of the present invention, there are variations in the specific embodiments and the application scope, and the present description should not be construed as limiting the present invention.
Claims (2)
1. A robust absolute orientation algorithm, characterized by: the method comprises the following steps:
s1, constructing an absolute orientation problem solving model based on a weighting coefficient,
let n spatial points have three-dimensional coordinates r in coordinate systems A and B A,i And r B,i (i=1, 2, …, n) solving the model by a weighted method using an absolute directional problem in weighted form
Is used as an estimate of the rotation matrix and translation vector, where ω i More than or equal to 0 represents the weight corresponding to the ith coordinate point andrepresenting a set of third-order orthogonal matrices of determinant 1, representation of 2 norms;
s2, setting the weight of each coordinate pointI.e. < ->
S3, deducing an analytic solution algorithm of an absolute directional problem solving model to obtain an optimal solution of a formula (2), calculating a rotation matrix R and a translation vector t corresponding to the weight value in S2,
the weighted centroid expression in the A, B coordinate system is recorded asAnd->(2) In the model of the formula, the solution of the rotation matrix and the translation vector is coupled, and new coordinates of an ith space point are introduced for decoupling>Andideal case->Wherein R is true And t true True values of rotation matrix and translation vector respectively, and then
r' B,i =R true r' A,i , (3)
And therefore, the optimization model is solved first
Obtaining a rotation matrix R corresponding to the current weight * Through again
Obtaining translation vector t in displacement parameter * Thereby realizing the aim of optimizing parameter decoupling,
due to
First two itemsIndependent of the rotation matrix, therefore, the optimization model of equation (4) is equivalent to
Wherein the method comprises the steps oftr () is the trace function of the matrix, and for solving the formula (7), singular value decomposition is performed on the matrix M to obtain +.>Wherein U and V are 3×3 orthogonal matrices, σ j 0 (j=1, 2, 3) is the singular value of matrix M, when
R * =V T U, (8)
(7) The maximum value is obtained and the maximum value is obtained,
to obtain R * Corresponding translation vector, deriving t from (5)Let the derivative be 0 to obtain a translation vector corresponding to the weight value as
S4, calculating the superposition error e i =||Rr' A,i -r' B,i After updating the weight by using the S estimation method, entering a step S3, calculating a rotation matrix corresponding to the current weight according to the formula (8), calculating a superposition error, and calculating a translation vector corresponding to the current weight according to the formula (9);
s5, continuously cycling the steps S4 and S3, and outputting a rotation matrix R and a translation vector t when the iteration number of the S estimation method exceeds the set maximum iteration number.
2. A robust absolute orientation algorithm according to claim 1, characterized in that: the S estimation method in step S4 specifically includes,
normalizing the error when the cycle number is 1Wherein->For the scale estimator +.>Represented |e i Median, 0.6754 is a constant in the S estimation, using
When the number of loops is greater than 1, the new weight is usedNormalization of the observation error>Using
As a new weight.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110789936.7A CN113569393B (en) | 2021-07-13 | 2021-07-13 | Robust absolute orientation algorithm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110789936.7A CN113569393B (en) | 2021-07-13 | 2021-07-13 | Robust absolute orientation algorithm |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113569393A CN113569393A (en) | 2021-10-29 |
CN113569393B true CN113569393B (en) | 2023-11-21 |
Family
ID=78164609
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110789936.7A Active CN113569393B (en) | 2021-07-13 | 2021-07-13 | Robust absolute orientation algorithm |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113569393B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115031927B (en) * | 2022-05-31 | 2023-07-07 | 西安电子科技大学 | High-precision positioning method for centroid of elliptical Gaussian distribution light spot |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106339603A (en) * | 2016-09-08 | 2017-01-18 | 武汉大学 | Relative orientation method based on axial angle vector |
CN107450045A (en) * | 2017-07-13 | 2017-12-08 | 中国人民解放军空军空降兵学院 | DOA estimation method based on the secondary weighted algorithms of FOCUSS |
FR3087387A1 (en) * | 2018-10-19 | 2020-04-24 | Compagnie Generale Des Etablissements Michelin | METHOD FOR SIMULATING THE TIME EVOLUTION OF A REAL-TIME PHYSICAL SYSTEM |
CN111784770A (en) * | 2020-06-28 | 2020-10-16 | 河北工业大学 | Three-dimensional attitude estimation method in disordered grabbing based on SHOT and ICP algorithm |
CN112818529A (en) * | 2021-01-22 | 2021-05-18 | 中国人民解放军国防科技大学 | Constrained space approximate orthogonal space filling test design method |
CN113051739A (en) * | 2021-03-17 | 2021-06-29 | 电子科技大学长三角研究院(衢州) | Robustness self-adaptive processing method based on sparse constraint |
-
2021
- 2021-07-13 CN CN202110789936.7A patent/CN113569393B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106339603A (en) * | 2016-09-08 | 2017-01-18 | 武汉大学 | Relative orientation method based on axial angle vector |
CN107450045A (en) * | 2017-07-13 | 2017-12-08 | 中国人民解放军空军空降兵学院 | DOA estimation method based on the secondary weighted algorithms of FOCUSS |
FR3087387A1 (en) * | 2018-10-19 | 2020-04-24 | Compagnie Generale Des Etablissements Michelin | METHOD FOR SIMULATING THE TIME EVOLUTION OF A REAL-TIME PHYSICAL SYSTEM |
CN111784770A (en) * | 2020-06-28 | 2020-10-16 | 河北工业大学 | Three-dimensional attitude estimation method in disordered grabbing based on SHOT and ICP algorithm |
CN112818529A (en) * | 2021-01-22 | 2021-05-18 | 中国人民解放军国防科技大学 | Constrained space approximate orthogonal space filling test design method |
CN113051739A (en) * | 2021-03-17 | 2021-06-29 | 电子科技大学长三角研究院(衢州) | Robustness self-adaptive processing method based on sparse constraint |
Non-Patent Citations (1)
Title |
---|
基于协方差矩阵估计的稳健Capon波束形成算法;陈明建;罗景青;龙国庆;;火力与指挥控制(第10期);98-101+106 * |
Also Published As
Publication number | Publication date |
---|---|
CN113569393A (en) | 2021-10-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111145227B (en) | Iterative integral registration method for space multi-view point cloud of underground tunnel | |
CN111795686A (en) | Method for positioning and mapping mobile robot | |
CN111798491A (en) | Maneuvering target tracking method based on Elman neural network | |
CN108801131B (en) | Filtering method and system for Beidou high-frequency deformation monitoring data | |
CN113569393B (en) | Robust absolute orientation algorithm | |
CN107870001A (en) | A kind of magnetometer bearing calibration based on ellipsoid fitting | |
CN112257722B (en) | Point cloud fitting method based on robust nonlinear Gaussian-Hermer model | |
CN110032709B (en) | Positioning and estimation method for abnormal point in geographic coordinate conversion | |
CN110900610A (en) | Industrial robot calibration method based on LM algorithm and particle filter algorithm optimization | |
CN111354040B (en) | Optical satellite image block adjustment method based on Partial EIV model | |
CN110375772A (en) | The ring laser stochastic error modeling of adaptive Kalman filter and compensation method | |
CN111553954B (en) | Online luminosity calibration method based on direct method monocular SLAM | |
CN108595373B (en) | Uncontrolled DEM registration method | |
CN112747773B (en) | Method for improving precision of gyroscope based on Allan variance and random polynomial | |
CN116091625A (en) | Binocular vision-based reference mark pose estimation method | |
CN113513978B (en) | High-precision measurement method and system for relative change of end face pose in high-temperature and low-temperature environment | |
CN114648544A (en) | Sub-pixel ellipse extraction method | |
CN112902844B (en) | Equipment circular port modeling method based on laser tracking | |
CN114092545A (en) | Self-adaptive grid searching method suitable for spherical target fitting | |
CN109270560B (en) | Multi-dimensional gross error positioning and value fixing method for area elevation abnormal data | |
CN108614804B (en) | Regularization Kalman filtering method based on signal-to-noise ratio test | |
CN112068092A (en) | Robust weighted observation fusion square root UKF filtering method for high-precision trajectory real-time orbit determination | |
CN115900511B (en) | Nonlinear separable least square-based magnetic dipole target positioning method | |
CN114140541B (en) | Parameter calibration method of multi-line structured light weld tracking sensor | |
CN115808160A (en) | Coordinate conversion method based on Rodrigue formula |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |