CN116295525A - Multi-IMU external parameter calibration method - Google Patents
Multi-IMU external parameter calibration method Download PDFInfo
- Publication number
- CN116295525A CN116295525A CN202310252703.2A CN202310252703A CN116295525A CN 116295525 A CN116295525 A CN 116295525A CN 202310252703 A CN202310252703 A CN 202310252703A CN 116295525 A CN116295525 A CN 116295525A
- Authority
- CN
- China
- Prior art keywords
- imu
- calibrated
- representing
- imus
- angular velocity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000013519 translation Methods 0.000 claims abstract description 12
- 238000005295 random walk Methods 0.000 claims abstract description 9
- 230000001133 acceleration Effects 0.000 claims description 16
- 239000011159 matrix material Substances 0.000 claims description 13
- 230000005484 gravity Effects 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 230000005653 Brownian motion process Effects 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 description 2
- 230000004927 fusion Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000000272 proprioceptive effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
-
- 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
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
Abstract
The invention relates to a multi-IMU external parameter calibration method, which comprises the following steps: constructing an IMU output model; calibrating random walk and noise of each IMU on the IMU array; recording IMU array data; taking any IMU from the IMU array as a reference IMU, taking a coordinate system of the reference IMU as a reference coordinate system, establishing a nonlinear least square problem about rotation quaternion between the reference IMU and other IMUs to be calibrated, and solving; generating virtual angular velocity observation values of the rest IMUs to be calibrated according to the calibrated relative rotation external parameters; and establishing a nonlinear least square problem between the reference IMU and the rest IMUs to be calibrated with respect to translation according to the virtual angular velocity observation values of the rest IMUs to be calibrated, and solving. The invention is independent of environment and external sensors, and can rapidly complete calibration.
Description
Technical Field
The invention relates to the technical field of automatic driving, in particular to a multi-IMU external parameter calibration method.
Background
In the field of simultaneous localization and mapping (SLAM), multi-inertial sensor fusion is a hot topic. The inertial measurement unit IMU can output an observed value at a high frequency as a proprioceptive sensor, but has a problem of posture drift in long-time positioning. Thus, the IMU may cooperate with external sensors (e.g., lidar and cameras) to provide global observations. Because of the advantages of compact size, low cost, and the like of microelectromechanical system (MEMS) IMUs, SLAM systems can add more IMUs for master-slave backup or to improve positioning accuracy.
Most SLAM studies ignore the extrinsic accuracy of multiple IMUs, but it plays an important role in fusion algorithms, and most vision-multiple IMU systems assume that the extrinsic between each IMU and sensor module is perfectly calibrated. However, according to simulation experiments, if the sensor external parameters with enough accuracy cannot be ensured, only one IMU is used, and the positioning accuracy is better.
To date, existing external calibration methods for multiple IMUs require accurate self-trajectories, relying on expensive turrets or external sensors (e.g., cameras). While these algorithms may perform well in certain environments, they are limited to the external environment and additional sensor devices.
Disclosure of Invention
The invention aims to solve the technical problem of providing a multi-IMU external parameter calibration method which is independent of environment and external sensors and can quickly complete calibration.
The technical scheme adopted for solving the technical problems is as follows: the multi-IMU external parameter calibration method comprises the following steps:
an IMU output model is built, wherein the output model comprises random walk and noise;
calibrating random walk and noise of each IMU on the IMU array;
recording IMU array data;
taking any IMU from the IMU array as a reference IMU, taking a coordinate system of the reference IMU as a reference coordinate system, establishing a nonlinear least square problem about rotation quaternion between the reference IMU and other IMUs to be calibrated, and solving and completing rotation external parameter calibration of the other IMUs to be calibrated relative to the reference IMU;
generating virtual angular velocity observation values of the rest IMUs to be calibrated according to the calibrated relative rotation external parameters;
and establishing a nonlinear least square problem about translation between the reference IMU and the other IMUs to be calibrated according to the virtual angular velocity observation values of the other IMUs to be calibrated, and solving and completing the translation external parameter calibration of the other IMUs to be calibrated relative to the reference IMU.
The IMU output model isWherein I represents a sensor coordinate system, W represents a world coordinate system, < >>And->Respectively representing a triaxial angular velocity value output by a gyroscope and a triaxial linear acceleration value output by an accelerometer in the IMU at the time t; omega (t) and a (t) represent the observed value of the gyroscope and the observed value of the accelerometer in the IMU at time t, respectively; b g (t) and b a (t) representing the deviation of the observed value of the gyroscope and the deviation of the observed value of the accelerometer in the IMU at time t, respectively, and modeling as a random walk; η (eta) g (t) and eta a (t) modeling as wiener processes, representing noise of the observations of the gyroscopes and noise of the observations of the accelerometers in the IMU at time t, respectively; r represents a rotation matrix in a world coordinate system; g represents a gravity vector.
When the IMU array data is recorded, the IMU array data is excited to each axis in an 8-shaped mode, and 200 frames of data are recorded, wherein the data comprise angular velocity and linear acceleration.
The nonlinear least squares problem between the reference IMU and the remaining IMUs with respect to rotation is:wherein, A represents a reference IMU, B represents an IMU to be calibrated, B q A representing a rotation quaternion between the reference IMU and the remaining IMUs to be calibrated, T representing a time set of samples, < ->Is the residual of the angular velocity observation,χ represents the system state quantity of the IMU to be estimated, +.>A first covariance matrix is represented by,wherein sigma g Representing the variance of the gravity estimation +.>Representing the variance of the bias estimate, Δt represents the time sampling interval, I 3 Representing the identity matrix of 3*3.
The virtual angular velocity observation value of the IMU to be calibrated passesCalculated, wherein->Representing the virtual angular velocity observation, freq represents the sampling frequency.
The nonlinear least square problem about translation between the reference IMU and the remaining IMUs to be calibrated is:wherein,, B p A representing translation parameters between the reference IMU and the remaining IMUs to be calibrated, < >>Is the residual of the angular velocity and linear acceleration observations,/-> A p B Representing the position of the reference IMU->A second covariance matrix is represented by a second covariance matrix,
σ a representing the variance of the acceleration noise,representing the variance of the acceleration bias.
Advantageous effects
Due to the adoption of the technical scheme, compared with the prior art, the invention has the following advantages and positive effects: the invention does not need IMU turntable equipment and the like, and does not depend on an external sensor, and compared with the prior method, the invention can achieve better calibration accuracy within 2 seconds by only collecting 200 frames of data. Compared with the prior art, the method has the advantages of higher speed and higher applicability under the same data, and can be used for successfully calibrating various IMU array experiments.
Drawings
FIG. 1 is a flow chart of a multi-IMU extrinsic parameter calibration method according to an embodiment of the present invention.
Detailed Description
The invention will be further illustrated with reference to specific examples. It is to be understood that these examples are illustrative of the present invention and are not intended to limit the scope of the present invention. Further, it is understood that various changes and modifications may be made by those skilled in the art after reading the teachings of the present invention, and such equivalents are intended to fall within the scope of the claims appended hereto.
The embodiment of the invention relates to a multi-IMU external parameter calibration method, as shown in fig. 1, comprising the following steps:
step 1, an IMU output model is built, the IMU generally comprises a gyroscope and an accelerometer, wherein the gyroscope outputs a three-axis angular velocity value, the accelerometer outputs a three-axis linear acceleration value, and the output model is as follows:
wherein I represents a sensor coordinate system, W represents a world coordinate system,and->Respectively representing a triaxial angular velocity value output by a gyroscope and a triaxial linear acceleration value output by an accelerometer in the IMU at the time t; omega (t) and a (t) respectively represent the observed value of a gyroscope and the observed value of an accelerometer in the IMU at the moment t, and R represents a rotation matrix under a world coordinate system; g represents a gravity vector.
b g (t) and b a (t) represents the deviation of the observed value of the gyroscope and the deviation of the observed value of the accelerometer in the IMU at time t, respectively, and is modeled as a random walk:
wherein,,and->Noise estimates representing angular velocity and acceleration bias, respectively,/->And->Respectively representing angular velocity and acceleration bias variance, I 3 Representing the identity matrix of 3*3.
η g (t) and eta a (t) noise representing the gyroscope observations and the accelerometer observations in the IMU at time t, respectively, modeled as a wiener process:
wherein sigma g Sum sigma a The variances of the angular velocities are respectively represented.
And 2, calibrating random walk and noise parameters of each IMU on the IMU array through Kalibr.
And 3, recording IMU array data, wherein the recorded data comprises angular velocity and linear acceleration, the IMU array data can be excited to each axis by 8 words, the recorded data volume is about 200 frames of data, for example, the IMU with the frequency of 100Hz is required to record 2 seconds of data.
And 4, taking any IMU from the IMU array as a reference IMU, taking a coordinate system of the reference IMU as a reference coordinate system, establishing a nonlinear least square problem about rotation quaternion between the reference IMU and the other IMU to be calibrated, and completing the rotation external parameter calibration of the other IMU to be calibrated relative to the reference IMU by solving the nonlinear least square problem. For example: the A-th IMU is taken as a reference IMU, the B-th IMU is taken as an IMU to be calibrated, and rotation quaternion between the A-th IMU and the B-th IMU is defined B q A Nonlinear least squares problem of (2):
and 5, generating virtual angular velocity observation values of the rest IMUs to be calibrated according to the calibrated relative rotation external parameters, wherein the virtual angular velocity observation values are generated by taking the B-th IMU as an example as follows:
wherein,,representing the virtual angular velocity observation, freq represents the sampling frequency.
And 6, establishing a nonlinear least square problem about translation between the reference IMU and the other IMUs to be calibrated according to the virtual angular velocity observation values of the other IMUs to be calibrated, and solving and completing the translation external parameter calibration of the other IMUs to be calibrated relative to the reference IMU. For example: the A-th IMU is taken as a reference IMU, the B-th IMU is taken as an IMU to be calibrated, and translation is defined between the A-th IMU and the B-th IMU B p A Nonlinear least squares problem of (2):
wherein,, A p B representing the location of the reference IMU.
in this embodiment, the nonlinear least squares problem can be solved using Google Ceres or OpenSLAM g2o tools.
It is not difficult to find that the invention does not need IMU turntable equipment and the like, and meanwhile, does not depend on an external sensor, compared with the existing method, the invention can achieve better calibration accuracy within 2 seconds by only collecting 200 frames of data. Compared with the prior art, the method has the advantages of higher speed and higher applicability under the same data, and can be used for successfully calibrating various IMU array experiments.
Claims (6)
1. The multi-IMU external parameter calibration method is characterized by comprising the following steps of:
an IMU output model is built, wherein the output model comprises random walk and noise;
calibrating random walk and noise of each IMU on the IMU array;
recording IMU array data;
taking any IMU from the IMU array as a reference IMU, taking a coordinate system of the reference IMU as a reference coordinate system, establishing a nonlinear least square problem about rotation quaternion between the reference IMU and other IMUs to be calibrated, and solving and completing rotation external parameter calibration of the other IMUs to be calibrated relative to the reference IMU;
generating virtual angular velocity observation values of the rest IMUs to be calibrated according to the calibrated relative rotation external parameters;
and establishing a nonlinear least square problem about translation between the reference IMU and the other IMUs to be calibrated according to the virtual angular velocity observation values of the other IMUs to be calibrated, and solving and completing the translation external parameter calibration of the other IMUs to be calibrated relative to the reference IMU.
2. The multi-IMU extrinsic reference calibration method according to claim 1, wherein said IMU output model isWherein I represents a sensor coordinate system, W represents a world coordinate system, < >>Andrespectively representing a triaxial angular velocity value output by a gyroscope and a triaxial linear acceleration value output by an accelerometer in the IMU at the time t; omega (t) and a (t) represent the observed value of the gyroscope and the observed value of the accelerometer in the IMU at time t, respectively; b g (t) and b a (t) representing the deviation of the observed value of the gyroscope and the deviation of the observed value of the accelerometer in the IMU at time t, respectively, and modeling as a random walk; η (eta) g (t) and eta a (t) modeling as wiener processes, representing noise of the observations of the gyroscopes and noise of the observations of the accelerometers in the IMU at time t, respectively; r represents a rotation matrix in a world coordinate system; g represents a gravity vector.
3. The method of claim 1, wherein the recording of IMU array data is performed by exciting each axis in an 8-word manner, and the total recording of 200 frames of data, including angular velocity and linear acceleration.
4. The multi-IMU extrinsic calibration method according to claim 2, wherein a nonlinear least squares problem between the reference IMU and the remaining IMUs with respect to rotation is:wherein, A represents a reference IMU, B represents an IMU to be calibrated, B q A representing the rotation quaternion between the reference IMU and the remaining IMUs to be calibrated, T representing the time set of samples,is the residual of the angular velocity observation, +.>χ represents the system state quantity of the IMU to be estimated, +.>Representing a first covariance matrix,>wherein sigma g Representing the variance of the gravity estimation +.>Representing the variance of the bias estimate, Δt represents the time sampling interval, I 3 Representing the identity matrix of 3*3.
6. The multi-IMU extrinsic calibration method according to claim 5, wherein a nonlinear least squares problem between the reference IMU and the remaining IMUs to be calibrated with respect to translation is:wherein,, B p A representing translation parameters between the reference IMU and the remaining IMUs to be calibrated, < >>Is the residual of the angular velocity and linear acceleration observations,/-> A p B Representing the position of the reference IMU->Representing a second covariance matrix,>wherein sigma a Representing acceleration noise variance->Representing the variance of the acceleration bias.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310252703.2A CN116295525A (en) | 2023-03-16 | 2023-03-16 | Multi-IMU external parameter calibration method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310252703.2A CN116295525A (en) | 2023-03-16 | 2023-03-16 | Multi-IMU external parameter calibration method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116295525A true CN116295525A (en) | 2023-06-23 |
Family
ID=86799106
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310252703.2A Pending CN116295525A (en) | 2023-03-16 | 2023-03-16 | Multi-IMU external parameter calibration method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116295525A (en) |
-
2023
- 2023-03-16 CN CN202310252703.2A patent/CN116295525A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Accelerometer for mobile robot positioning | |
Akeila et al. | Reducing low-cost INS error accumulation in distance estimation using self-resetting | |
EP1760431A1 (en) | Vehicle comprising an inertial navigation system with a plurality of Kalman filters | |
CN110956665B (en) | Bidirectional calculation method, system and device for turning track of vehicle | |
CN108007477B (en) | Inertial pedestrian positioning system error suppression method based on forward and reverse filtering | |
US20180066943A1 (en) | Fused Sensor Ensemble for Navigation and Calibration Process Therefor | |
Chen et al. | Improving inertial sensor by reducing errors using deep learning methodology | |
Ariffin et al. | Low cost MEMS gyroscope and accelerometer implementation without Kalman Filter for angle estimation | |
CN107990901B (en) | User direction positioning method based on sensor | |
Amirsadri et al. | Practical considerations in precise calibration of a low-cost MEMS IMU for road-mapping applications | |
US20180180420A1 (en) | Method and System for Improving Inertial Measurement Unit Sensor Signals | |
WO2018182528A1 (en) | Trajectory estimation system and method | |
Yuan et al. | Indoor pedestrian navigation using miniaturized low-cost MEMS inertial measurement units | |
CN113218389B (en) | Vehicle positioning method, device, storage medium and computer program product | |
CN116295525A (en) | Multi-IMU external parameter calibration method | |
CN101949702B (en) | Quick self-testing method for GNSS PVT quality by using MEMS accelerometer | |
CN113959464B (en) | Gyroscope-assisted accelerometer field calibration method and system | |
Kowalczuk et al. | Modelling an accelerometer for robot position estimation | |
CN114137592A (en) | Switching method and system for fusion positioning of multi-source sensor | |
Vaispacher et al. | Application of linearized Kalman filter in integration of navigation systems | |
KR20230048324A (en) | Extended dead reckoning accuracy | |
CN114184190A (en) | Inertial/odometer integrated navigation system and method | |
Akeila et al. | A self-resetting method for reducing error accumulation in INS-based tracking | |
CN113075665A (en) | Underwater positioning method, underwater vehicle navigation device and computer readable storage medium | |
CN110082805A (en) | A kind of 3 D locating device and method |
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 |