CN109242912A - Join scaling method, electronic equipment, storage medium outside acquisition device - Google Patents
Join scaling method, electronic equipment, storage medium outside acquisition device Download PDFInfo
- Publication number
- CN109242912A CN109242912A CN201810996832.1A CN201810996832A CN109242912A CN 109242912 A CN109242912 A CN 109242912A CN 201810996832 A CN201810996832 A CN 201810996832A CN 109242912 A CN109242912 A CN 109242912A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- coordinate
- acquisition device
- fitting
- movable platform
- 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 58
- 238000003860 storage Methods 0.000 title claims description 7
- 238000006243 chemical reaction Methods 0.000 claims abstract description 97
- 230000009466 transformation Effects 0.000 claims abstract description 30
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 20
- 239000003550 marker Substances 0.000 claims description 45
- 239000011159 matrix material Substances 0.000 claims description 34
- 238000013519 translation Methods 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 16
- 238000005457 optimization Methods 0.000 claims description 13
- 230000006870 function Effects 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000011065 in-situ storage Methods 0.000 claims description 3
- 238000009434 installation Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明提供采集装置外参标定方法,包括步骤:获取采集装置采集的标识信息;计算标识坐标系与采集装置坐标系的第一坐标转换关系;根据第一坐标转换关系得到标识坐标系下标识信息中的标记点在采集装置坐标系下对应的坐标点,采用拟合算法拟合坐标点,建立拟合坐标系,计算拟合坐标系与采集装置坐标系之间的第二坐标转换关系;根据预设轨迹对应的拟合算法拟合坐标点,生成拟合轨迹,计算拟合坐标系与可移动平台坐标系之间的第三坐标转换关系;通过第二坐标转换关系和第三坐标转换关系,计算采集装置坐标系与可移动平台坐标系之间的转换关系。本发明能够减小采集装置安装不平整带来的安装误差,通过优化标定结果使得标定结果更准确,标定误差小。
The invention provides a method for calibrating external parameters of a collection device, which comprises the steps of: acquiring identification information collected by the collection device; calculating a first coordinate conversion relationship between the identification coordinate system and the collection device coordinate system; obtaining identification information under the identification coordinate system according to the first coordinate conversion relationship The corresponding coordinate points of the marked points in the acquisition device coordinate system, the fitting algorithm is used to fit the coordinate points, the fitting coordinate system is established, and the second coordinate conversion relationship between the fitting coordinate system and the acquisition device coordinate system is calculated; The fitting algorithm corresponding to the preset trajectory fits the coordinate points, generates the fitting trajectory, and calculates the third coordinate transformation relationship between the fitted coordinate system and the movable platform coordinate system; through the second coordinate transformation relationship and the third coordinate transformation relationship , calculate the conversion relationship between the coordinate system of the acquisition device and the coordinate system of the movable platform. The invention can reduce the installation error caused by the uneven installation of the collecting device, and by optimizing the calibration result, the calibration result is more accurate and the calibration error is small.
Description
技术领域technical field
本发明涉及采集装置外参标定技术领域,尤其涉及采集装置外参标定方法、电子设备、存储介质。The invention relates to the technical field of external parameter calibration of a collection device, in particular to a method, electronic equipment and storage medium for the calibration of external parameters of a collection device.
背景技术Background technique
通常需要预先对采集装置进行标定以便获得准确的采集装置参数。采集装置参数包括内部参数和外部参数,采集装置内部参数是采集装置本身固有的参数,例如焦距、光心位置等,采集装置外参表示采集装置在世界坐标系(也称全局坐标系)中的位置和方向,从数学角度而言,采集装置外参表征了采集装置坐标系与指定的世界坐标系(也称为全局坐标系)之间的坐标变换关系。采集装置位姿估计是在采集装置与被摄物体的相对运动过程中估计采集装置的外参,而估计算法中使用到的采集装置内参数通常已在离线情况下通过采集装置标定方法确定,传统的采集装置-轮式里程计标定方式对轮式里程计有较高精度要求,而且误差较大。Usually, the acquisition device needs to be calibrated in advance in order to obtain accurate acquisition device parameters. The parameters of the acquisition device include internal parameters and external parameters. The internal parameters of the acquisition device are inherent parameters of the acquisition device, such as focal length, optical center position, etc. The external parameters of the acquisition device represent the acquisition device in the world coordinate system (also called the global coordinate system). Position and orientation, from a mathematical point of view, the external parameters of the acquisition device represent the coordinate transformation relationship between the coordinate system of the acquisition device and the specified world coordinate system (also called the global coordinate system). The pose estimation of the acquisition device is to estimate the external parameters of the acquisition device during the relative motion of the acquisition device and the object to be photographed, and the internal parameters of the acquisition device used in the estimation algorithm are usually determined offline by the acquisition device calibration method. The acquisition device-wheel odometer calibration method has high precision requirements for the wheel odometer, and the error is large.
发明内容SUMMARY OF THE INVENTION
为了克服现有技术的不足,本发明的目的在于提供采集装置外参标定方法,直接标定运动中心,实现客观自由度完整标定,优化标定结果,准确度高,误差小。In order to overcome the deficiencies of the prior art, the purpose of the present invention is to provide a method for calibrating the external parameters of the acquisition device, directly calibrating the motion center, realizing the complete calibration of the objective degrees of freedom, optimizing the calibration result, with high accuracy and small error.
本发明提供采集装置外参标定方法,包括以下步骤:The present invention provides a method for calibrating external parameters of a collection device, comprising the following steps:
获取标识信息,获取采集装置在若干采集时刻采集的信息,所述信息为可移动平台沿预设轨迹与标识相对移动过程中的标识信息;Acquiring identification information, acquiring information collected by the collection device at several collection moments, the information being the identification information during the relative movement of the movable platform along the preset track and the identification;
建立第一坐标转换关系,计算每个采集时刻标识坐标系与采集装置坐标系的转换关系,记为第一坐标转换关系;Establish a first coordinate conversion relationship, calculate the conversion relationship between the coordinate system of each collection moment and the coordinate system of the collection device, and record it as the first coordinate conversion relationship;
建立第二坐标转换关系,在所述标识坐标系上获取所述标识信息中的标记点,根据所述第一坐标转换关系得到所述标记点每个采集时刻在所述采集装置坐标系下对应的若干坐标点,采用拟合算法拟合所述坐标点,建立拟合坐标系,计算所述拟合坐标系与所述采集装置坐标系之间的转换关系,记为第二坐标转换关系;Establish a second coordinate conversion relationship, obtain the marking points in the marking information on the marking coordinate system, and obtain the corresponding marking points in the coordinate system of the collecting device at each acquisition moment of the marking points according to the first coordinate conversion relationship Several coordinate points of , use a fitting algorithm to fit the coordinate points, establish a fitting coordinate system, calculate the conversion relationship between the fitting coordinate system and the coordinate system of the acquisition device, and record it as the second coordinate conversion relationship;
建立第三坐标转换关系,根据所述预设轨迹对应的拟合算法拟合所述坐标点,生成拟合轨迹,计算所述拟合坐标系与可移动平台坐标系之间的转换关系,记为第三坐标转换关系;Establish a third coordinate conversion relationship, fit the coordinate points according to the fitting algorithm corresponding to the preset trajectory, generate a fitted trajectory, calculate the conversion relationship between the fitted coordinate system and the movable platform coordinate system, record is the third coordinate conversion relationship;
采集装置外参标定,通过第二坐标转换关系和第三坐标转换关系,计算所述采集装置坐标系与所述可移动平台坐标系之间的转换关系。In the external parameter calibration of the acquisition device, the conversion relationship between the coordinate system of the acquisition device and the coordinate system of the movable platform is calculated through the second coordinate conversion relationship and the third coordinate conversion relationship.
进一步地,所述建立第一坐标转换关系中,对所述标识信息进行处理,建立图像坐标系,得到所述标识在所述图像坐标系和所述标识坐标系上的坐标,并结合采集装置内参和畸变系数,通过位姿估计方法计算所述标识坐标系与所述采集装置坐标系之间的转换关系。Further, in the establishment of the first coordinate conversion relationship, the identification information is processed, an image coordinate system is established, the coordinates of the identification on the image coordinate system and the identification coordinate system are obtained, and the collection device is combined. The internal parameters and distortion coefficients are used to calculate the transformation relationship between the identification coordinate system and the acquisition device coordinate system through a pose estimation method.
进一步地,所述预设轨迹包括所述可移动平台以所述标识为参考系生成的旋转轨迹与直行轨迹。Further, the preset trajectory includes a rotation trajectory and a straight trajectory generated by the movable platform using the identifier as a reference frame.
进一步地,所述获取标识信息包括以下步骤:Further, the obtaining identification information includes the following steps:
获取旋转图像,获取在可移动平台按所述旋转轨迹原地旋转过程中采集装置采集的标识信息,获得若干旋转图像;Acquiring a rotating image, acquiring the identification information collected by the collecting device during the in-situ rotation of the movable platform according to the rotating trajectory, and obtaining several rotating images;
获取直行图像,获取在可移动平台按所述直行轨迹沿所述可移动平台坐标系的坐标轴直线行驶过程中采集装置采集的标识信息,获得若干直行图像。Acquiring straight images, acquiring identification information collected by the collecting device during the process of moving the movable platform in a straight line along the coordinate axis of the coordinate system of the movable platform according to the straight trajectory, and obtaining several straight images.
进一步地,所述建立第一坐标转换关系包括以下步骤:Further, the establishment of the first coordinate conversion relationship includes the following steps:
生成旋转位姿,对所述旋转图像进行处理,获得所述标记点的图像坐标,根据所述标记点的图像坐标及所述标记点对应所述标记板坐标系标识坐标系的坐标,并结合采集装置内参和畸变系数,计算所述标记板坐标系标识坐标系相对于所述采集装置坐标系下的位姿,获得旋转位姿;Generate a rotation pose, process the rotated image, obtain the image coordinates of the marker point, identify the coordinate system of the marker board coordinate system according to the image coordinates of the marker point and the marker point corresponding to the coordinate system, and combine Collecting the internal parameters and distortion coefficients of the device, calculating the pose of the marker plate coordinate system and the marker coordinate system relative to the coordinate system of the collecting device, and obtaining the rotational pose;
生成直行位姿,对所述直行图像进行处理,获得所述标记点的图像坐标,根据所述标记点的图像坐标及所述标记点对应所述标记板坐标系标识坐标系的坐标,并结合采集装置内参和畸变系数,计算所述标记板坐标系标识坐标系相对于采集装置坐标系下的位姿,获得直行位姿。Generate a straight pose, process the straight image, obtain the image coordinates of the marker, identify the coordinates of the coordinate system according to the image coordinates of the marker and the marker plate coordinate system corresponding to the marker, and combine The internal parameters and distortion coefficients of the acquisition device are acquired, and the pose of the marker coordinate system of the marker coordinate system relative to the coordinate system of the acquisition device is calculated to obtain the straight-ahead pose.
进一步地,所述建立第二坐标转换关系中,采用平面拟合算法拟合所述坐标点,生成拟合平面,在所述拟合平面内建立拟合坐标系。Further, in establishing the second coordinate conversion relationship, a plane fitting algorithm is used to fit the coordinate points, a fitting plane is generated, and a fitting coordinate system is established in the fitting plane.
进一步地,所述生成第二坐标转换关系包括以下步骤:Further, the generating of the second coordinate conversion relationship includes the following steps:
计算所述坐标点与所述拟合平面之间的距离,当所述距离大于阈值时,剔除所述坐标点,并重新拟合剩余坐标点;Calculate the distance between the coordinate point and the fitting plane, when the distance is greater than a threshold, remove the coordinate point, and re-fit the remaining coordinate points;
建立拟合坐标系,在所述拟合平面内取第一点和第二点,所述第一点和所述第二点在所述采集装置坐标系的坐标分别为第一坐标和第二坐标,根据所述第一坐标和所述第二坐标计算第一单位向量、第二单位向量、第三单位向量,所述第一单位向量为第一坐标点与第二坐标点差向量的单位向量,所述第三单位向量为所述拟合平面的法向量,所述第二单位向量为所述第一单位向量与所述第三单位向量的向量积,以所述第一点为坐标原点,分别以所述第一单位向量、第二单位向量、第三单位向量的方向为x轴、y轴、z轴的方向建立拟合坐标系;A fitting coordinate system is established, a first point and a second point are taken in the fitting plane, and the coordinates of the first point and the second point in the coordinate system of the acquisition device are the first coordinate and the second point respectively Coordinate, calculate the first unit vector, the second unit vector, and the third unit vector according to the first coordinate and the second coordinate, and the first unit vector is the unit vector of the difference vector between the first coordinate point and the second coordinate point , the third unit vector is the normal vector of the fitting plane, the second unit vector is the vector product of the first unit vector and the third unit vector, and the first point is the coordinate origin , and establish a fitting coordinate system with the directions of the first unit vector, the second unit vector, and the third unit vector as the directions of the x-axis, the y-axis, and the z-axis;
生成坐标转换关系,根据所述第一单位向量、第二单位向量、第三单位向量及所述第一坐标生成所述拟合坐标系与所述采集装置坐标系之间的转换关系。A coordinate conversion relationship is generated, and a conversion relationship between the fitting coordinate system and the acquisition device coordinate system is generated according to the first unit vector, the second unit vector, the third unit vector, and the first coordinates.
进一步地,所述建立第三坐标转换关系中,对所述拟合坐标系下沿所述旋转轨迹得到的坐标点进行椭圆拟合,记拟合坐标系中的椭圆中心为可移动平台的旋转中心,得到拟合平面坐标系与可移动平台坐标系之间的平移矩阵;对所述拟合坐标系下沿直行轨迹得到的坐标点进行直线拟合,得到拟合平面坐标系与可移动平台坐标系之间的旋转矩阵,得到所述拟合坐标系与可移动平台坐标系之间的转换关系。Further, in the establishment of the third coordinate conversion relationship, ellipse fitting is performed on the coordinate points obtained along the rotation trajectory under the fitting coordinate system, and the center of the ellipse in the fitting coordinate system is recorded as the rotation of the movable platform. center to obtain the translation matrix between the fitted plane coordinate system and the movable platform coordinate system; perform straight line fitting on the coordinate points obtained along the straight trajectory under the fitted coordinate system to obtain the fitted plane coordinate system and the movable platform The rotation matrix between the coordinate systems obtains the conversion relationship between the fitting coordinate system and the movable platform coordinate system.
进一步地,通过第二坐标转换关系、第三坐标转换关系计算可移动平台的旋转中心在所述采集装置坐标系下的坐标,得到采集装置坐标系与可移动平台坐标系间的平移矩阵;Further, calculating the coordinates of the rotation center of the movable platform under the coordinate system of the collection device through the second coordinate conversion relationship and the third coordinate conversion relationship, and obtaining a translation matrix between the coordinate system of the collection device and the movable platform coordinate system;
通过第三转换关系计算拟合平面坐标系与可移动平台坐标系之间的旋转矩阵,通过第二转换关系计算采集装置坐标系与可移动平台坐标系的旋转矩阵;Calculate the rotation matrix between the fitting plane coordinate system and the movable platform coordinate system through the third transformation relationship, and calculate the rotation matrix between the acquisition device coordinate system and the movable platform coordinate system through the second transformation relationship;
由所述平移矩阵和旋转矩阵计算所述采集装置坐标系与所述可移动平台坐标系之间的转换关系。The transformation relationship between the coordinate system of the acquisition device and the coordinate system of the movable platform is calculated from the translation matrix and the rotation matrix.
进一步地,还包括优化标定结果,通过所述采集装置坐标系与所述可移动平台坐标系之间的转换关系、所述第一坐标转换关系构建代价函数,将所述代价函数输入优化库进行计算,生成优化标定结果。Further, it also includes optimizing the calibration result, constructing a cost function through the conversion relationship between the coordinate system of the acquisition device and the coordinate system of the movable platform, and the first coordinate conversion relationship, and inputting the cost function into an optimization library for execution. Calculate to generate optimized calibration results.
进一步地,所述优化标定结果还包括建立里程计约束,输入优化库进行计算,生成优化标定结果。Further, the optimization calibration result also includes establishing odometer constraints, inputting the optimization library for calculation, and generating the optimization calibration result.
进一步地,还包括建立重投影约束、圆约束、直线约束中的一种或多种优化标定结果。Further, it also includes establishing one or more optimal calibration results among reprojection constraints, circle constraints, and line constraints.
一种电子设备,包括:处理器;An electronic device, comprising: a processor;
存储器;以及程序,其中所述程序被存储在所述存储器中,并且被配置成由处理器执行,所述程序包括用于执行上述采集装置外参标定方法。a memory; and a program, wherein the program is stored in the memory and is configured to be executed by a processor, the program comprising for executing the above-mentioned method for calibrating an external parameter of an acquisition device.
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行上述采集装置外参标定方法。A computer-readable storage medium on which a computer program is stored, the computer program being executed by a processor to perform the above-mentioned method for calibrating an external parameter of an acquisition device.
相比现有技术,本发明的有益效果在于:Compared with the prior art, the beneficial effects of the present invention are:
本发明提供采集装置外参标定方法,采集装置外参标定方法,包括以下步骤:获取标识信息,获取采集装置在若干采集时刻采集的信息,信息为可移动平台沿预设轨迹与标识相对移动过程中的标识信息;建立第一坐标转换关系,计算每个采集时刻标识坐标系与采集装置坐标系的转换关系,记为第一坐标转换关系;建立第二坐标转换关系,在标识坐标系上获取标识信息中的标记点,根据第一坐标转换关系得到标记点每个采集时刻在采集装置坐标系下对应的若干坐标点,采用拟合算法拟合坐标点,建立拟合坐标系,计算拟合坐标系与采集装置坐标系之间的转换关系,记为第二坐标转换关系;建立第三坐标转换关系,根据预设轨迹对应的拟合算法拟合坐标点,生成拟合轨迹,计算拟合坐标系与可移动平台坐标系之间的转换关系,记为第三坐标转换关系;采集装置外参标定,通过第二坐标转换关系和第三坐标转换关系,计算采集装置坐标系与可移动平台坐标系之间的转换关系。本发明涉及电子设备与可读存储介质,用于执行采集装置外参标定方法;本发明通过增加拟合平面进行采集装置外参标定,减小因采集装置安装不平整带来机械上的安装误差,通过优化标定结果使得标定结果更准确,不依赖于高精度的里程计,标定误差小。The invention provides a method for calibrating external parameters of a collection device, and the method for calibrating external parameters of a collection device, comprising the following steps: acquiring identification information, and acquiring information collected by a collection device at several collection moments, the information being the relative movement process of the movable platform and the identification along a preset track Establish the first coordinate conversion relationship, calculate the conversion relationship between the logo coordinate system and the acquisition device coordinate system at each collection moment, and record it as the first coordinate conversion relationship; establish the second coordinate conversion relationship, and obtain it from the logo coordinate system For the marked points in the identification information, according to the first coordinate conversion relationship, several coordinate points corresponding to the marked points in the coordinate system of the acquisition device are obtained at each acquisition moment, and a fitting algorithm is used to fit the coordinate points, establish a fitting coordinate system, and calculate the fitting. The conversion relationship between the coordinate system and the coordinate system of the acquisition device is recorded as the second coordinate conversion relationship; the third coordinate conversion relationship is established, the coordinate points are fitted according to the fitting algorithm corresponding to the preset trajectory, the fitting trajectory is generated, and the fitting is calculated. The transformation relationship between the coordinate system and the coordinate system of the movable platform is recorded as the third coordinate transformation relationship; for the external parameter calibration of the acquisition device, the coordinate system of the acquisition device and the movable platform are calculated through the second coordinate transformation relationship and the third coordinate transformation relationship. Transformation relationship between coordinate systems. The invention relates to an electronic device and a readable storage medium, which are used to perform a method for calibrating external parameters of a collection device; the invention performs external parameter calibration of the collection device by adding a fitting plane, so as to reduce mechanical installation errors caused by uneven installation of the collection device , by optimizing the calibration results to make the calibration results more accurate, independent of the high-precision odometer, and the calibration error is small.
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,并可依照说明书的内容予以实施,以下以本发明的较佳实施例并配合附图详细说明如后。本发明的具体实施方式由以下实施例及其附图详细给出。The above description is only an overview of the technical solution of the present invention. In order to understand the technical means of the present invention more clearly, and implement it according to the content of the description, the preferred embodiments of the present invention are described in detail below with the accompanying drawings. Specific embodiments of the present invention are given in detail by the following examples and the accompanying drawings.
附图说明Description of drawings
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:The accompanying drawings described herein are used to provide a further understanding of the present invention and constitute a part of the present application. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation of the present invention. In the attached image:
图1为本发明的采集装置外参标定方法流程图;Fig. 1 is the flow chart of the external parameter calibration method of the acquisition device of the present invention;
图2为本发明实施例中的aprilTag码示意图;Fig. 2 is the schematic diagram of aprilTag code in the embodiment of the present invention;
图3为本发明实施例中的标记板示意图;3 is a schematic diagram of a marker board in an embodiment of the present invention;
图4为本发明实施例中的采集装置及可移动平台运动示意图;4 is a schematic diagram of the movement of the acquisition device and the movable platform in the embodiment of the present invention;
图5为本发明实施例中的拟合平面坐标系、采集装置坐标系及可移动平台坐标系示意图一;5 is a schematic diagram 1 of a fitting plane coordinate system, a collection device coordinate system, and a movable platform coordinate system in an embodiment of the present invention;
图6为本发明实施例中的拟合平面坐标系、采集装置坐标系及可移动平台坐标系示意图二;6 is a second schematic diagram of a fitting plane coordinate system, a collection device coordinate system, and a movable platform coordinate system in an embodiment of the present invention;
图7为本发明实施例中的采集装置及可移动平台运动示意图。FIG. 7 is a schematic diagram of the movement of the acquisition device and the movable platform in the embodiment of the present invention.
图中:1、可移动平台;2、采集装置;3、标记板;4、采集装置坐标系;5、拟合平面坐标系;6、旋转轨迹;7、旋转点;8、执行轨迹;9、直行点;10、可移动平台坐标系。In the figure: 1. Movable platform; 2. Acquisition device; 3. Marker board; 4. Coordinate system of acquisition device; 5. Fitted plane coordinate system; 6. Rotation track; 7. Rotation point; 8. Execution track; 9 , Straight point; 10. Coordinate system of movable platform.
具体实施方式Detailed ways
下面,结合附图以及具体实施方式,对本发明做进一步描述,需要说明的是,在不相冲突的前提下,以下描述的各实施例之间或各技术特征之间可以任意组合形成新的实施例。The present invention will be further described below with reference to the accompanying drawings and specific embodiments. It should be noted that, on the premise of no conflict, the embodiments or technical features described below can be combined arbitrarily to form new embodiments. .
采集装置外参标定方法,如图1所示,包括以下步骤:The external parameter calibration method of the acquisition device, as shown in Figure 1, includes the following steps:
获取标识信息,标记点在同一平面且规则排列在标记板上,标记板置于固定工作平台上,采集装置安装于可移动平台的中心位置附近,朝向固定工作平台,获取采集装置在若干采集时刻采集的信息,信息为可移动平台沿预设轨迹与标识相对移动过程中的标识信息,标识信息存储若干标记点。所述的标识用以定位,可以是一个或多个具有具体形状的特征、图案、二维码,其可以被识别并从中提取坐标点。二维码可以是AprilTag或其他形式。所述的标识也可以是二维码和其他图案、色彩的组合。具体而言,所述的标识是明显的,且具备被识别特征的特性。当其采用二维码时,可以排列成特定形状的阵列。其二维码数量越多,定位越准确。所述的采集装置可以是相机,也可以是其他具有识别功能的传感器。本实施例中,以aprilTag码作为标记物,如图2所示,aprilTag码为黑色,背景通常为白色,aprilTag码的每个角点作为一个标记点,每个aprilTag码具有固定的ID号和朝向,如图3所示,若干aprilTag码在同一平面且规则排列在标记板上,如图4所示,本实施例中,在可移动平台1运动过程中,采集装置2至少能够观测到标记板3上的4个标记点,标记点之间的物理坐标关系已知,可以是任何区别背景,能够在图像中被识别且唯一确定的点,即每个标记点在一幅图像中具有确定的唯一坐标。Obtain the identification information, the marking points are on the same plane and regularly arranged on the marking plate, the marking plate is placed on the fixed working platform, the collecting device is installed near the center of the movable platform, facing the fixed working platform, and the collecting device is acquired at several collecting moments The collected information is the identification information during the relative movement of the movable platform along the preset track and the mark, and the identification information stores several marking points. The identification is used for positioning, and can be one or more features with specific shapes, patterns, and two-dimensional codes, which can be identified and coordinate points can be extracted therefrom. The QR code can be AprilTag or other forms. The identification can also be a combination of two-dimensional code and other patterns and colors. Specifically, the identification is obvious and has the characteristics of the identified feature. When it uses a two-dimensional code, it can be arranged in an array of a specific shape. The more QR codes, the more accurate the positioning. The collection device may be a camera, or may be other sensors with identification functions. In this embodiment, the aprilTag code is used as a marker. As shown in Figure 2, the aprilTag code is black, and the background is usually white. Each corner of the aprilTag code is used as a marker point, and each aprilTag code has a fixed ID number and Orientation, as shown in Figure 3, several aprilTag codes are on the same plane and regularly arranged on the marking plate, as shown in Figure 4, in this embodiment, during the movement of the movable platform 1, the collecting device 2 can at least observe the markings The 4 markers on the board 3, the physical coordinate relationship between the markers is known, can be any distinguishing background, can be identified in the image and uniquely determined points, that is, each marker has a certain number in an image. the unique coordinates of .
在一实施例中,优选的,预设轨迹包括可移动平台以标识为参考系生成的旋转轨迹与直行轨迹。优选的,获取标识信息包括以下步骤:In an embodiment, preferably, the preset trajectory includes a rotation trajectory and a straight trajectory generated by the movable platform with the identification as the reference frame. Preferably, acquiring the identification information includes the following steps:
获取旋转图像,获取在可移动平台按旋转轨迹原地旋转过程中采集装置采集的标识信息,获得若干旋转图像;本实施例中,使可移动平台原地旋转一周,在可移动平台旋转过程中,采集装置至少可观测到一个aprilTag码,采集装置获取图像。Obtain a rotating image, obtain identification information collected by the collection device during the in-situ rotation of the movable platform according to the rotation trajectory, and obtain several rotating images; , the acquisition device can observe at least one aprilTag code, and the acquisition device acquires the image.
获取直行图像,获取在可移动平台按直行轨迹沿可移动平台坐标系的坐标轴直线行驶过程中采集装置采集的标识信息,获得若干直行图像。本实施例中,使可移动平台沿可移动坐标系的x轴直线行驶一段距离,直行时,采集装置至少可观测到一个apriltag码,采集装置获取图像。A straight image is acquired, the identification information collected by the collecting device is acquired during the process of the movable platform traveling straight along the coordinate axis of the coordinate system of the movable platform according to the straight trajectory, and several straight images are obtained. In this embodiment, the movable platform is made to travel in a straight line along the x-axis of the movable coordinate system for a certain distance. When traveling straight, the acquisition device can observe at least one apriltag code, and the acquisition device acquires an image.
建立第一坐标转换关系,计算每个采集时刻标识坐标系与采集装置坐标系的转换关系,记为第一坐标转换关系。优选的,建立第一坐标转换关系中,对标识信息进行处理,建立图像坐标系,得到标识在图像坐标系和标识坐标系上的坐标,并结合采集装置内参和畸变系数,通过位姿估计方法计算标识坐标系与采集装置坐标系之间的转换关系。本实施例中,利用棋盘格标定板与MATLAB中的采集装置标定工具标定采集装置的内参及畸变参数,通过图像处理手段识别apriltag码上四个标记点的图像坐标根据标记点的图像坐标及标记点对应标记板坐标系的坐标,并结合采集装置内参和畸变系数,通过位姿估计方法如PNP位姿估计方法计算标记板坐标系与采集装置坐标系之间的转换关系。A first coordinate conversion relationship is established, and the conversion relationship between the coordinate system of the marker at each collection moment and the coordinate system of the collection device is calculated, which is recorded as the first coordinate conversion relationship. Preferably, in establishing the first coordinate conversion relationship, the identification information is processed, an image coordinate system is established, and the coordinates of the identification on the image coordinate system and the identification coordinate system are obtained, and combined with the internal parameters and distortion coefficients of the acquisition device, the pose estimation method is used. Calculate the conversion relationship between the identification coordinate system and the acquisition device coordinate system. In this embodiment, the internal parameters and distortion parameters of the acquisition device are calibrated by using the checkerboard calibration board and the acquisition device calibration tool in MATLAB, and the image coordinates of the four marked points on the apriltag code are identified by image processing methods according to the image coordinates and marks of the marked points. The points correspond to the coordinates of the coordinate system of the marker board, and combined with the internal parameters and distortion coefficients of the acquisition device, the transformation relationship between the coordinate system of the marker board and the coordinate system of the acquisition device is calculated by a pose estimation method such as the PNP pose estimation method.
在一实施例中,优选的,建立第一坐标转换关系包括以下步骤:In an embodiment, preferably, establishing the first coordinate conversion relationship includes the following steps:
生成旋转位姿,对旋转图像进行处理,获得标记点的图像坐标,根据标记点的图像坐标及标记点对应标记板坐标系标识坐标系的坐标,并结合采集装置内参和畸变系数,计算标记板坐标系标识坐标系相对于采集装置坐标系下的位姿,获得旋转位姿;旋转位姿包括位置和姿态,位置用平移向量表示,姿态用旋转向量表示,在旋转过程中,假设一共保存了M个旋转位姿。Generate the rotation pose, process the rotated image, and obtain the image coordinates of the marked point. According to the image coordinates of the marked point and the coordinate system of the marking plate corresponding to the marked point, the coordinates of the coordinate system are identified, and the internal parameters of the acquisition device and the distortion coefficient are used to calculate the marked plate. The coordinate system identifies the coordinate system relative to the pose in the coordinate system of the acquisition device, and obtains the rotation pose; the rotation pose includes the position and the attitude, the position is represented by a translation vector, and the attitude is represented by a rotation vector. During the rotation process, it is assumed that a total of M rotation poses.
生成直行位姿,对直行图像进行处理,获得标记点的图像坐标,根据标记点的图像坐标及标记点对应标记板坐标系标识坐标系的坐标,并结合采集装置内参和畸变系数,计算标记板坐标系标识坐标系相对于采集装置坐标系下的位姿,获得直行位姿。直行位姿包括位置和姿态,位置用平移向量表示,姿态用旋转向量表示,假设一共保存了N个直行位姿。Generate the straight pose, process the straight image, and obtain the image coordinates of the marked points. According to the image coordinates of the marked points and the coordinates of the marking plate coordinate system corresponding to the marked points, the coordinates of the coordinate system are identified, and the internal parameters of the acquisition device and the distortion coefficient are used to calculate the marked plate. The coordinate system identifies the coordinate system relative to the pose in the coordinate system of the acquisition device, and obtains the straight-line pose. The straight pose includes position and attitude. The position is represented by a translation vector, and the attitude is represented by a rotation vector. It is assumed that a total of N straight poses are saved.
在一实施例中,世界坐标系用符号Ow-XwYwZw表示,世界坐标系是定义在空间中的绝对坐标系,不会随着采集装置或者物体的移动而改变。可移动平台坐标系用符号Ob-XbYbZb表示,可移动平台的前进方向与x轴正方向平行。采集装置坐标系用符号Oc-XcYcZc表示。标记板坐标系用符号Om-XmYmZm表示。拟合平面坐标系用符号Op-XpYpZp表示。图像坐标系,用符号OI-XIYI表示。设某一旋转矩阵为aRb,它是一个3×3单位矩阵,表示相对于坐标系a下的坐标系b的旋转,通过罗德里格斯变化可得到该旋转矩阵对应的旋转向量aψb(aαb,aβb,aγb)T,T表示矩阵的转置,同样的,旋转向量亦可通过罗德里格斯公式得到旋转矩阵。设某一平移向量为atb,它是一个3×1的向量,表示相对于坐标系a下的坐标系b的平移关系。设某个坐标为api(axi,ayi,azi),表示坐标系a下第i个点的坐标。设第i(i=1,…,M+N)次获取图像时,采集装置总共观测到n个标记点,其中第k(k=1,…,n)个标记点的图像坐标表示为Iipk(uk,vk),即图像坐标系中第i获取的图像中的第k标记点坐标,标记板坐标系下坐标mipk(mixk,miyk,mizk),即标记板坐标系中第i获取的图像中的第k点坐标,得到采集装置坐标系与标定板坐标之间的转换为ctmi,cψmi,若在旋转过程得到的位姿也可表示为同样,在直行过程中得到的位姿也可表示为 In one embodiment, the world coordinate system is represented by the symbol O w -X w Y w Z w , and the world coordinate system is an absolute coordinate system defined in space, which will not change with the movement of the acquisition device or the object. The coordinate system of the movable platform is represented by the symbol O b -X b Y b Z b , and the forward direction of the movable platform is parallel to the positive direction of the x-axis. The coordinate system of the acquisition device is represented by the symbol O c -X c Y c Z c . The marker plate coordinate system is denoted by the notation O m -X m Y m Z m . The fitted plane coordinate system is denoted by the notation Op - X p Y p Z p . Image coordinate system, represented by the symbol O I -X I Y I. Let a rotation matrix be a R b , which is a 3×3 unit matrix, which represents the rotation relative to the coordinate system b under the coordinate system a, and the rotation vector a ψ corresponding to the rotation matrix can be obtained through the Rodrigues change b ( a α b , a β b , a γ b ) T , where T represents the transpose of the matrix. Similarly, the rotation vector can also be obtained by the Rodrigues formula to obtain the rotation matrix. Let a certain translation vector be a t b , which is a 3×1 vector, which represents the translation relationship with respect to the coordinate system b under the coordinate system a. Let a certain coordinate be a p i ( a x i , a y i , a z i ), which means the coordinates of the i-th point in the coordinate system a. It is assumed that when the image is acquired for the i-th (i=1,...,M+N) time, the acquisition device observes a total of n markers, wherein the image coordinates of the k-th (k=1,...,n) marker are denoted as Ii p k (u k ,v k ), that is, the coordinates of the kth marker point in the image obtained by the ith in the image coordinate system, and the coordinates in the marker plate coordinate system mi p k ( mi x k , mi y k , mi z k ) , that is, the coordinate of the kth point in the image obtained by the ith in the coordinate system of the marker plate, and the conversion between the coordinate system of the acquisition device and the coordinate of the calibration plate is obtained as ct mi , c ψ mi , if the pose obtained during the rotation process is also can be expressed as Similarly, the pose obtained in the process of going straight can also be expressed as
建立第二坐标转换关系,在标识坐标系上获取标识信息中的标记点,根据第一坐标转换关系得到标记点每个采集时刻在采集装置坐标系下对应的若干坐标点,采用拟合算法拟合坐标点,建立拟合坐标系,计算拟合坐标系与采集装置坐标系之间的转换关系,记为第二坐标转换关系;优选的,建立第二坐标转换关系中,采用平面拟合算法拟合坐标点,生成拟合平面,在拟合平面内建立拟合坐标系。优选的,生成第二坐标转换关系包括以下步骤:Establish a second coordinate conversion relationship, obtain the marked points in the identification information on the identification coordinate system, obtain a number of coordinate points corresponding to the marked point in the coordinate system of the acquisition device at each acquisition moment according to the first coordinate conversion relationship, and use a fitting algorithm to fit the points. Combine the coordinate points, establish a fitting coordinate system, calculate the conversion relationship between the fitting coordinate system and the acquisition device coordinate system, and record it as the second coordinate conversion relationship; preferably, in establishing the second coordinate conversion relationship, a plane fitting algorithm is used Fit the coordinate points, generate a fitting plane, and establish a fitting coordinate system in the fitting plane. Preferably, generating the second coordinate conversion relationship includes the following steps:
计算坐标点与拟合平面之间的距离,当距离大于阈值时,剔除坐标点,并重新拟合剩余坐标点;本实施例中阈值为0.0001m,当距离大于0.0001m时,判定坐标点为离散点,剔除离散点,并重新拟合剩余坐标点,直到不再存在离散点。Calculate the distance between the coordinate point and the fitting plane. When the distance is greater than the threshold, remove the coordinate point and re-fit the remaining coordinate points; in this embodiment, the threshold is 0.0001m, and when the distance is greater than 0.0001m, the coordinate point is determined to be Discrete points, remove the discrete points, and refit the remaining coordinate points until there are no more discrete points.
建立拟合平面坐标系,在拟合平面内取第一点和第二点,第一点以p1表示,第二点以p2表示,第一点和第二点在采集装置坐标系的坐标分别为第一坐标和第二坐标,第一坐标为第二坐标为根据第一坐标和第二坐标计算第一单位向量、第二单位向量、第三单位向量,第一单位向量为第一坐标点与第二坐标点差向量的单位向量,第三单位向量为拟合平面的法向量,第二单位向量为第一单位向量与第三单位向量的向量积,以第一点为坐标原点,分别以第一单位向量、第二单位向量、第三单位向量的方向为x轴、y轴、z轴的方向建立拟合坐标系;具体地,第一单位向量第三单位向量第二单位向量向量方向平行于p2p1方向,向量方向是平面法向量方向,p2p1位于拟合平面内,向量和向量相互垂直,根据正交原理得到向量,即向量方向分别垂直于向量和向量如图5所示,若可移动平台坐标系10的z轴与采集装置坐标系4的z轴同向,则以点p1为坐标原点,设定x轴正方向沿单位向量方向,y轴正方向沿单位向量方向,z轴正方向沿单位向量方向建立拟合平面坐标系5,得到拟合平面坐标系5与采集装置坐标系4的转换关系:Establish a fitting plane coordinate system, take the first point and the second point in the fitting plane, the first point is represented by p1, the second point is represented by p2, the coordinates of the first point and the second point in the coordinate system of the acquisition device are respectively is the first coordinate and the second coordinate, the first coordinate is The second coordinate is Calculate the first unit vector, the second unit vector, and the third unit vector according to the first coordinate and the second coordinate. The first unit vector is the unit vector of the difference vector between the first coordinate point and the second coordinate point, and the third unit vector is the fitting The normal vector of the plane, the second unit vector is the vector product of the first unit vector and the third unit vector, with the first point as the coordinate origin, and the directions of the first unit vector, the second unit vector, and the third unit vector are The directions of the x-axis, y-axis, and z-axis establish a fitting coordinate system; specifically, the first unit vector third unit vector second unit vector vector The direction is parallel to the p 2 p 1 direction, the vector direction is the plane normal vector direction, p 2 p 1 lies in the fitted plane, the vector and vector perpendicular to each other, according to the principle of orthogonality vector, i.e. vector The directions are perpendicular to the vector and vector As shown in FIG. 5 , if the z-axis of the movable platform coordinate system 10 and the z-axis of the acquisition device coordinate system 4 are in the same direction, the point p 1 is taken as the coordinate origin, and the positive direction of the x-axis is set along the unit vector direction, the positive y-axis is along the unit vector direction, the positive z-axis is along the unit vector The direction establishes the fitting plane coordinate system 5, and obtains the conversion relationship between the fitting plane coordinate system 5 and the acquisition device coordinate system 4:
如图6所示,若可移动平台坐标系10的z轴与采集装置坐标系4的z轴反向,则以点p1为坐标原点,设定x轴正方向沿单位向量方向,y轴负方向沿单位向量方向,z轴负方向沿单位向量方向建立拟合平面坐标系5,得到拟合平面坐标系5与采集装置坐标系4的转换关系:As shown in FIG. 6 , if the z-axis of the movable platform coordinate system 10 is opposite to the z-axis of the acquisition device coordinate system 4 , the point p1 is taken as the coordinate origin, and the positive direction of the x-axis is set along the unit vector direction, the negative y-axis is along the unit vector direction, the negative z-axis is along the unit vector The direction establishes the fitting plane coordinate system 5, and obtains the conversion relationship between the fitting plane coordinate system 5 and the acquisition device coordinate system 4:
其中,cRp,ctp分别表示采集装置坐标系相对拟合平面坐标系的旋转变化矩阵和平移矩阵。则拟合平面坐标系相对采集装置坐标系的旋转变化矩阵和平移矩阵为:Wherein, c R p and c t p respectively represent the rotation change matrix and translation matrix of the coordinate system of the acquisition device relative to the fitted plane coordinate system. Then the rotation change matrix and translation matrix of the fitted plane coordinate system relative to the acquisition device coordinate system are:
pRc=(cRp)T,ptc=-(pRc·ctp) p R c =( c R p ) T , p t c =-( p R c · c t p )
生成坐标转换关系,根据第一单位向量、第二单位向量、第三单位向量及第一坐标生成拟合平面坐标系与采集装置坐标系之间的转换关系。A coordinate conversion relationship is generated, and a conversion relationship between the fitted plane coordinate system and the acquisition device coordinate system is generated according to the first unit vector, the second unit vector, the third unit vector and the first coordinates.
建立第三坐标转换关系,根据预设轨迹对应的拟合算法拟合坐标点,生成拟合轨迹,计算拟合坐标系与可移动平台坐标系之间的转换关系,记为第三坐标转换关系;例如根据平面拟合算法拟合平面,根据椭圆拟合算法拟合椭圆,根据直线拟合算法拟合直线。如:通过旋转位姿和直行位姿分别将旋转图像和直行图像中标记点对应标记板坐标系的坐标转换为采集装置坐标系下的坐标,获得采集装置坐标系下的坐标点,由于可移动平台旋转得到的椭圆为2D椭圆,而采集的标记点为3D点,采用平面拟合算法拟合坐标点,坐标点坐标等效于平移向量,拟合平面方程为:Establish a third coordinate conversion relationship, fit the coordinate points according to the fitting algorithm corresponding to the preset trajectory, generate a fitting trajectory, calculate the conversion relationship between the fitted coordinate system and the movable platform coordinate system, and record it as the third coordinate conversion relationship ; For example, the plane is fitted according to the plane fitting algorithm, the ellipse is fitted according to the ellipse fitting algorithm, and the straight line is fitted according to the straight line fitting algorithm. For example, the coordinates of the marker plate coordinate system corresponding to the marker points in the rotating image and the straight image are converted into the coordinates under the coordinate system of the acquisition device through the rotation pose and the straight pose, respectively, and the coordinate points under the coordinate system of the acquisition device are obtained. The ellipse obtained by the platform rotation is a 2D ellipse, and the collected mark points are 3D points. The plane fitting algorithm is used to fit the coordinate points. The coordinates of the coordinate points are equivalent to the translation vector. The fitting plane equation is:
Ap·x+Bp·y+Cp·z+Dp=0A p x+B p y+C p z+D p =0
其中,(Ap,Bp,Cp)为平面法向量,则平面与x轴相交于与y轴相交于与z轴相交于Ap,Bp,Cp,Dp均可由拟合平面得到。还可以根据采集装置模型以及映射关系,将采集装置坐标系变换到其他根据需要确定的坐标系。Among them, (A p , B p , C p ) is the plane normal vector, then the plane intersects the x-axis at intersects the y-axis at intersects the z-axis at A p , B p , C p , D p can be obtained from the fitted plane. It is also possible to transform the coordinate system of the acquisition device into other coordinate systems determined as required according to the acquisition device model and the mapping relationship.
在一实施例中,优选的,建立第三坐标转换关系中,对拟合坐标系下沿旋转轨迹得到的坐标点进行椭圆拟合,记拟合坐标系中的椭圆中心为可移动平台的旋转中心,得到拟合平面坐标系与可移动平台坐标系之间的平移矩阵;对拟合坐标系下沿直行轨迹得到的坐标点进行直线拟合,得到拟合平面坐标系与可移动平台坐标系之间的旋转矩阵,得到拟合坐标系与可移动平台坐标系之间的转换关系。In one embodiment, preferably, in establishing the third coordinate conversion relationship, ellipse fitting is performed on the coordinate points obtained along the rotation trajectory under the fitting coordinate system, and the center of the ellipse in the fitting coordinate system is the rotation of the movable platform. center to obtain the translation matrix between the fitted plane coordinate system and the movable platform coordinate system; perform straight line fitting on the coordinate points obtained along the straight trajectory under the fitted coordinate system to obtain the fitted plane coordinate system and the movable platform coordinate system The rotation matrix between the fitting coordinate system and the movable platform coordinate system is obtained.
采集装置外参标定,通过第二坐标转换关系和第三坐标转换关系,计算采集装置坐标系与可移动平台坐标系之间的转换关系。优选的,通过第二坐标转换关系、第三坐标转换关系计算可移动平台的旋转中心在采集装置坐标系下的坐标,得到采集装置坐标系与可移动平台坐标系间的平移矩阵;通过第三转换关系计算拟合平面坐标系与可移动平台坐标系之间的旋转矩阵,通过第二转换关系计算采集装置坐标系与可移动平台坐标系的旋转矩阵;由平移矩阵和旋转矩阵计算采集装置坐标系与可移动平台坐标系之间的转换关系。如图7所示,根据采集装置坐标系4与拟合平面坐标系5之间的转换关系,设第i个点在采集装置坐标系4下的坐标为则其在拟合平面坐标系5下的坐标为:The external parameters of the acquisition device are calibrated, and the conversion relationship between the coordinate system of the acquisition device and the coordinate system of the movable platform is calculated through the second coordinate conversion relationship and the third coordinate conversion relationship. Preferably, the coordinates of the rotation center of the movable platform in the coordinate system of the acquisition device are calculated through the second coordinate conversion relationship and the third coordinate conversion relationship, so as to obtain a translation matrix between the coordinate system of the acquisition device and the coordinate system of the movable platform; The transformation relationship calculates the rotation matrix between the fitting plane coordinate system and the movable platform coordinate system, and calculates the rotation matrix between the acquisition device coordinate system and the movable platform coordinate system through the second transformation relationship; calculates the acquisition device coordinates from the translation matrix and the rotation matrix. The transformation relationship between the system and the coordinate system of the movable platform. As shown in FIG. 7 , according to the conversion relationship between the acquisition device coordinate system 4 and the fitting plane coordinate system 5 , the coordinates of the i-th point in the acquisition device coordinate system 4 are set as Then its coordinates in the fitting plane coordinate system 5 are:
ppi=pRc·cpi+ptc p p i = p R c · c p i + p t c
由于拟合前已去除离散点,旋转点7或直行点9在拟合平面坐标系5下的z坐标趋近于0,可认为在拟合平面坐标系5下的坐标点都位于x-y平面上,取拟合平面坐标系5下旋转点7的x、y轴坐标进行椭圆拟合,生成旋转轨迹6,椭圆拟合的方式采用椭圆拟合函数进行拟合,输出结果为椭圆中心在相应坐标系下的坐标,假设为(xc,yc),则该点在拟合平面坐标系下的坐标为(xc,yc,0),根据相对运动原理,认定椭圆中心为可移动平台的旋转中心,通过坐标变化得到旋转中心在采集装置坐标系下的坐标:Since the discrete points have been removed before fitting, the z-coordinate of the rotation point 7 or the straight point 9 in the fitting plane coordinate system 5 is close to 0, and it can be considered that the coordinate points in the fitting plane coordinate system 5 are all located on the xy plane , take the x and y-axis coordinates of the rotation point 7 under the fitting plane coordinate system 5 to perform ellipse fitting, and generate a rotation trajectory 6. The ellipse fitting method adopts the ellipse fitting function to fit, and the output result is that the center of the ellipse is at the corresponding coordinate The coordinates in the coordinate system, assuming that it is (x c , y c ), then the coordinates of the point in the fitting plane coordinate system are (x c , y c , 0). According to the principle of relative motion, the center of the ellipse is identified as the movable platform The rotation center of , and the coordinates of the rotation center in the coordinate system of the acquisition device are obtained through the coordinate change:
ctb=cRp·[xc,yc,0]T+ctp c t b = c R p ·[x c ,y c ,0] T + c t p
取拟合平面坐标5下直行点9的x、y轴坐标进行直线拟合,生成直线轨迹8,直线拟合的方式采用直线拟合函数进行拟合,输出结果为直线方程系数,直线方程表示为:Take the x and y-axis coordinates of the straight point 9 under the fitting plane coordinate 5 to perform straight line fitting, and generate a straight line trajectory 8. The straight line fitting method uses a straight line fitting function for fitting, and the output result is the straight line equation coefficient, which is represented by the straight line equation. for:
Bl·y=Al·x+Cl B l ·y=A l ·x+C l
根据直线方程系数Al,Bl得到直线与x轴的夹角,即前进方向与x轴之间的角度:According to the coefficients A l and B l of the straight line equation, the angle between the straight line and the x-axis is obtained, that is, the angle between the forward direction and the x-axis:
直线方向为可移动平台的前进方向,得到拟合平面坐标系与可移动平台坐标系之间的旋转矩阵:The straight line direction is the forward direction of the movable platform, and the rotation matrix between the fitted plane coordinate system and the movable platform coordinate system is obtained:
则采集装置坐标系与可移动平台坐标系之间的旋转关系:Then the rotation relationship between the coordinate system of the acquisition device and the coordinate system of the movable platform:
cRb=cRp·pRb c R b = c R p · p R b
即采集装置坐标系与可移动平台坐标系之间的转换关系:旋转矩阵为cRb,平移矩阵为ctb,为了方便下一过程优化,将旋转矩阵转化为旋转向量cψb(cαb,cβb,cγb)。That is, the conversion relationship between the coordinate system of the acquisition device and the coordinate system of the movable platform: the rotation matrix is c R b , and the translation matrix is c t b . In order to facilitate the next process optimization, the rotation matrix is converted into a rotation vector c ψ b ( c α b , c β b , c γ b ).
在一实施例中,优选的,还包括步骤优化标定结果,根据相对运动原理,通过采集装置坐标系与可移动平台坐标系之间的转换关系、第一坐标转换关系构建代价函数,将代价函数输入优化库进行计算,生成优化标定结果。优选的,优化标定结果还包括建立里程计约束,输入优化库进行计算,生成优化标定结果。优选的,还包括建立重投影约束、圆约束、直线约束中的一种或多种优化标定结果。采集装置坐标系和可移动平台坐标系之间的转换为:ctb(cxb,cyb),cψb(cαb,cβb,cγb),即采集装置光心与可移动平台的坐标之间的位置关系;旋转直行过程中,采集装置坐标系下标记板坐标系之间的转换为:为了区别旋转和直行过程,旋转过程中采集装置坐标系下标记板坐标系之间的转换为:直行过程中采集装置坐标系下标记板坐标系之间的转换为: 直线系数为:Al,Bl,圆的半径为r,优化变量总自由度为6M+6N+9。本实施例中,代价函数具体为:In one embodiment, it is preferable to further include steps to optimize the calibration result, according to the principle of relative motion, construct a cost function through the conversion relationship between the coordinate system of the acquisition device and the coordinate system of the movable platform, and the first coordinate conversion relationship, and convert the cost function into the cost function. Input the optimization library for calculation, and generate the optimization calibration result. Preferably, optimizing the calibration result further includes establishing an odometer constraint, inputting it into an optimization library for calculation, and generating an optimized calibration result. Preferably, it also includes establishing one or more optimal calibration results among reprojection constraints, circle constraints, and line constraints. The conversion between the coordinate system of the collecting device and the coordinate system of the movable platform is: c t b ( c x b , c y b ), c ψ b ( c α b , c β b , c γ b ), that is, the light of the collecting device The positional relationship between the center and the coordinates of the movable platform; in the process of rotating straight, the conversion between the coordinate systems of the marker plate under the coordinate system of the acquisition device is: In order to distinguish between the rotation and the straight travel process, the conversion between the coordinate systems of the marker plate under the coordinate system of the acquisition device during the rotation process is: The conversion between the coordinate systems of the marker plate under the coordinate system of the acquisition device during the straight travel is: The coefficients of the straight line are: A l , B l , the radius of the circle is r, and the total degree of freedom of the optimization variable is 6M+6N+9. In this embodiment, the cost function is specifically:
1、直线约束,可得到N个方程: 1. Straight line constraints, N equations can be obtained:
2、圆约束,可得到M个方程: 2. Circle constraints, M equations can be obtained:
3、重投影约束,可得到2(M+N)个方程:3. Reprojection constraints, 2(M+N) equations can be obtained:
设第i(i=1,…,M+N)次获取图像时,采集装置总共观测到n个标记点,其中第j(j=1,…,n)个标记点的图像坐标表示为Iipj(uj,vj),标记板坐标系下坐标mipk(mixk,miyk,mizk),得到采集装置坐标系与标定板坐标系之间的转换为ctmi,cψmi。重投影过程,首先计算采集装置坐标系下第j个标记点的坐标:Assuming that when the image is acquired for the ith (i=1,...,M+N) time, the acquisition device observes n markers in total, wherein the image coordinates of the jth (j=1,...,n) marker are denoted as Ii p j (u j , v j ), the coordinates mi p k ( mi x k , mi y k , mi z k ) under the coordinate system of the marker board, the conversion between the coordinate system of the acquisition device and the coordinate system of the calibration board is obtained as c t mi , c ψ mi . In the reprojection process, first calculate the coordinates of the jth marker point in the coordinate system of the acquisition device:
其次,通过畸变模型得到畸变后的图像坐标,理论上观测到的图像坐标应与重投影得到的图像坐标相等:Secondly, the distorted image coordinates are obtained through the distortion model. In theory, the observed image coordinates should be equal to the image coordinates obtained by reprojection:
dist=1+r·(k1+k2·r)dist=1+r·(k 1 +k 2 ·r)
4、里程计约束(轨迹约束),可得到6(M+N)个方程:4. Odometer constraint (trajectory constraint), 6(M+N) equations can be obtained:
对于里程计,可直接获取世界坐标系下各个时刻的位置和姿态,For the odometer, the position and attitude of each moment in the world coordinate system can be directly obtained,
即(wtb0,wRb0),…,(wtbi,wRbi),(wtbj,wRbj),得到第j个轨迹点相对于第i个轨迹点的旋转平移:That is ( wt b0 , w R b0 ), ..., ( wt bi , w R bi ), ( wt bj , w R bj ), the rotational translation of the j-th trajectory point relative to the i-th trajectory point is obtained:
平移:bitbj=(wRbi)T·(wtbj-wtbi)Translation: bit bj = ( w R bi ) T · ( w t bj - w t bi )
旋转:biRbj=(wRbi)T·wRbj Rotation: bi R bj = ( w R bi ) T · w R bj
对于采集装置,获取到采集装置坐标系和标记板坐标系之间的位姿(ctm0,cRm0),…,(ctmi,cRmi),(ctmj,cRmj),若将标定板坐标系固定,上述位姿等价于(c0tm,c0Rm),…,(citm,ciRm),(cjtm,cjRm),通过采集装置坐标系和可移动平台坐标系得到上述轨迹:For the acquisition device, the poses ( c t m0 , cR m0 ), ( c t mi , c R mi ), ( c t mj , c R mj ) between the coordinate system of the acquisition device and the marker plate coordinate system are obtained. , if the calibration plate coordinate system is fixed, the above pose is equivalent to ( c0 t m , c0 R m ), ..., ( ci t m , ci R m ), ( cj t m , cj R m ), through the acquisition device The coordinate system and the movable platform coordinate system get the above trajectory:
平移:bitbj=bRc·(ciRm·(bRc·cjRm)T·(-btc)+ciRm·(cjRm)T·(-cjtm)+citm)+btc Translation: bit t bj = b R c · ( ci R m · ( b R c · cj R m ) T · (- b t c )+ ci R m · ( cj R m ) T · (- cj t m ) + ci t m )+ b t c
旋转:biRbj=bRc·ciRm·(bRc·cjRm)T Rotation: bi R bj = b R c · ci R m · ( b R c · cj R m ) T
共得到9(M+N)个方程,而变量的自由度只有6N+6个,所以满足优化条件。将四个约束用优化库去求取最后结果,得到优化标定结果。A total of 9(M+N) equations are obtained, and the degrees of freedom of the variables are only 6N+6, so the optimization conditions are satisfied. The four constraints are used to obtain the final result with the optimization library, and the optimized calibration result is obtained.
一种电子设备,包括:处理器;An electronic device, comprising: a processor;
存储器;以及程序,其中程序被存储在存储器中,并且被配置成由处理器执行,程序包括用于执行上述采集装置外参标定方法。a memory; and a program, wherein the program is stored in the memory and is configured to be executed by the processor, the program includes a method for performing the above-mentioned external parameter calibration method of the acquisition device.
一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行上述采集装置外参标定方法。A computer-readable storage medium on which a computer program is stored, and the computer program is executed by a processor to perform the above-mentioned method for calibrating external parameters of an acquisition device.
本发明提供采集装置外参标定方法,采集装置外参标定方法,包括以下步骤:获取标识信息,获取采集装置在若干采集时刻采集的信息,信息为可移动平台沿预设轨迹与标识相对移动过程中的标识信息;建立第一坐标转换关系,计算每个采集时刻标识坐标系与采集装置坐标系的转换关系,记为第一坐标转换关系;建立第二坐标转换关系,在标识坐标系上获取标识信息中的标记点,根据第一坐标转换关系得到标记点每个采集时刻在采集装置坐标系下对应的若干坐标点,采用拟合算法拟合坐标点,建立拟合坐标系,计算拟合坐标系与采集装置坐标系之间的转换关系,记为第二坐标转换关系;建立第三坐标转换关系,根据预设轨迹对应的拟合算法拟合坐标点,生成拟合轨迹,计算拟合坐标系与可移动平台坐标系之间的转换关系,记为第三坐标转换关系;采集装置外参标定,通过第二坐标转换关系和第三坐标转换关系,计算采集装置坐标系与可移动平台坐标系之间的转换关系。本发明涉及电子设备与可读存储介质,用于执行采集装置外参标定方法;本发明通过增加拟合平面进行采集装置外参标定,减小因采集装置安装不平整带来机械上的安装误差,通过优化标定结果使得标定结果更准确,不依赖于高精度的里程计,标定误差小。The invention provides a method for calibrating external parameters of a collection device, and the method for calibrating external parameters of a collection device, comprising the following steps: acquiring identification information, and acquiring information collected by a collection device at several collection moments, the information being the relative movement process of the movable platform and the identification along a preset track Establish the first coordinate conversion relationship, calculate the conversion relationship between the logo coordinate system and the acquisition device coordinate system at each collection moment, and record it as the first coordinate conversion relationship; establish the second coordinate conversion relationship, and obtain it from the logo coordinate system For the marked points in the identification information, according to the first coordinate conversion relationship, several coordinate points corresponding to the marked points in the coordinate system of the acquisition device are obtained at each acquisition moment, and a fitting algorithm is used to fit the coordinate points, establish a fitting coordinate system, and calculate the fitting. The conversion relationship between the coordinate system and the coordinate system of the acquisition device is recorded as the second coordinate conversion relationship; the third coordinate conversion relationship is established, the coordinate points are fitted according to the fitting algorithm corresponding to the preset trajectory, the fitting trajectory is generated, and the fitting is calculated. The transformation relationship between the coordinate system and the coordinate system of the movable platform is recorded as the third coordinate transformation relationship; for the external parameter calibration of the acquisition device, the coordinate system of the acquisition device and the movable platform are calculated through the second coordinate transformation relationship and the third coordinate transformation relationship. Transformation relationship between coordinate systems. The invention relates to an electronic device and a readable storage medium, which are used to perform a method for calibrating external parameters of a collection device; the invention performs external parameter calibration of the collection device by adding a fitting plane, so as to reduce mechanical installation errors caused by uneven installation of the collection device , by optimizing the calibration results to make the calibration results more accurate, independent of the high-precision odometer, and the calibration error is small.
以上,仅为本发明的较佳实施例而已,并非对本发明作任何形式上的限制;凡本行业的普通技术人员均可按说明书附图所示和以上而顺畅地实施本发明;但是,凡熟悉本专业的技术人员在不脱离本发明技术方案范围内,利用以上所揭示的技术内容而做出的些许更动、修饰与演变的等同变化,均为本发明的等效实施例;同时,凡依据本发明的实质技术对以上实施例所作的任何等同变化的更动、修饰与演变等,均仍属于本发明的技术方案的保护范围之内。The above are only preferred embodiments of the present invention, and do not limit the present invention in any form; any person of ordinary skill in the industry can smoothly implement the present invention as shown in the accompanying drawings and above; however, any Those skilled in the art, without departing from the scope of the technical solution of the present invention, make use of the above-disclosed technical content to make some changes, modifications and equivalent changes of evolution are equivalent embodiments of the present invention; at the same time, Any alteration, modification and evolution of any equivalent changes made to the above embodiments according to the essential technology of the present invention still fall within the protection scope of the technical solution of the present invention.
Claims (14)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810996832.1A CN109242912A (en) | 2018-08-29 | 2018-08-29 | Join scaling method, electronic equipment, storage medium outside acquisition device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810996832.1A CN109242912A (en) | 2018-08-29 | 2018-08-29 | Join scaling method, electronic equipment, storage medium outside acquisition device |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109242912A true CN109242912A (en) | 2019-01-18 |
Family
ID=65068707
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810996832.1A Pending CN109242912A (en) | 2018-08-29 | 2018-08-29 | Join scaling method, electronic equipment, storage medium outside acquisition device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109242912A (en) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109829952A (en) * | 2019-02-11 | 2019-05-31 | 广东工业大学 | A kind of outer ginseng scaling method, device and the relevant device of positioning system |
CN110207722A (en) * | 2019-06-11 | 2019-09-06 | 远形时空科技(北京)有限公司 | A kind of automation calibration for cameras mileage system and method |
CN111753565A (en) * | 2019-03-27 | 2020-10-09 | 北京外号信息技术有限公司 | Method and electronic equipment for presenting information related to optical communication device |
CN112082491A (en) * | 2020-09-11 | 2020-12-15 | 苏州杰锐思智能科技股份有限公司 | Height detection method based on point cloud |
CN112184824A (en) * | 2019-07-05 | 2021-01-05 | 杭州海康机器人技术有限公司 | Camera external parameter calibration method and device |
CN112258590A (en) * | 2020-12-08 | 2021-01-22 | 杭州迦智科技有限公司 | Laser-based depth camera external parameter calibration method, device and storage medium thereof |
CN112567439A (en) * | 2020-11-09 | 2021-03-26 | 驭势(上海)汽车科技有限公司 | Method and device for determining traffic flow information, electronic equipment and storage medium |
WO2021056503A1 (en) * | 2019-09-29 | 2021-04-01 | 深圳市大疆创新科技有限公司 | Positioning method and apparatus for movable platform, movable platform, and storage medium |
CN112654886A (en) * | 2020-05-27 | 2021-04-13 | 华为技术有限公司 | External parameter calibration method, device, equipment and storage medium |
CN112884842A (en) * | 2020-12-25 | 2021-06-01 | 杭州迦智科技有限公司 | Synchronous calibration method for zero position of double rudder wheels and external parameters of acquisition device and movable platform |
CN113063353A (en) * | 2021-03-31 | 2021-07-02 | 深圳中科飞测科技股份有限公司 | Coordinate system establishing method, detection device, detection equipment and storage medium |
CN113252066A (en) * | 2020-02-13 | 2021-08-13 | 纳恩博(北京)科技有限公司 | Method and device for calibrating parameters of odometer equipment, storage medium and electronic device |
CN113538598A (en) * | 2021-07-21 | 2021-10-22 | 北京能创科技有限公司 | Active Stereo Vision System Calibration Method |
CN113814987A (en) * | 2021-11-24 | 2021-12-21 | 季华实验室 | Multi-camera robot hand-eye calibration method and device, electronic equipment and storage medium |
CN114538027A (en) * | 2020-11-26 | 2022-05-27 | 合肥欣奕华智能机器股份有限公司 | Full-automatic visual positioning transfer equipment and control method thereof |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103279949A (en) * | 2013-05-09 | 2013-09-04 | 浙江大学 | Operation method of self-positioning robot-based multi-camera parameter automatic calibration system |
CN103542981A (en) * | 2013-09-28 | 2014-01-29 | 大连理工大学 | Method for measuring rotary inertia through binocular vision |
CN105869150A (en) * | 2016-03-24 | 2016-08-17 | 杭州南江机器人股份有限公司 | Movable platform calibration device and calibration method based on visual recognition |
CN106705956A (en) * | 2017-02-28 | 2017-05-24 | 南京工程学院 | Rapid industrial robot tail end pose measuring device and measuring method thereof |
CN107993266A (en) * | 2017-12-05 | 2018-05-04 | 西北工业大学 | Light-field camera intrinsic parameter caliberating device and method |
-
2018
- 2018-08-29 CN CN201810996832.1A patent/CN109242912A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103279949A (en) * | 2013-05-09 | 2013-09-04 | 浙江大学 | Operation method of self-positioning robot-based multi-camera parameter automatic calibration system |
CN103542981A (en) * | 2013-09-28 | 2014-01-29 | 大连理工大学 | Method for measuring rotary inertia through binocular vision |
CN105869150A (en) * | 2016-03-24 | 2016-08-17 | 杭州南江机器人股份有限公司 | Movable platform calibration device and calibration method based on visual recognition |
CN106705956A (en) * | 2017-02-28 | 2017-05-24 | 南京工程学院 | Rapid industrial robot tail end pose measuring device and measuring method thereof |
CN107993266A (en) * | 2017-12-05 | 2018-05-04 | 西北工业大学 | Light-field camera intrinsic parameter caliberating device and method |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109829952A (en) * | 2019-02-11 | 2019-05-31 | 广东工业大学 | A kind of outer ginseng scaling method, device and the relevant device of positioning system |
CN111753565A (en) * | 2019-03-27 | 2020-10-09 | 北京外号信息技术有限公司 | Method and electronic equipment for presenting information related to optical communication device |
CN111753565B (en) * | 2019-03-27 | 2021-12-24 | 北京外号信息技术有限公司 | Method and electronic equipment for presenting information related to optical communication device |
CN110207722A (en) * | 2019-06-11 | 2019-09-06 | 远形时空科技(北京)有限公司 | A kind of automation calibration for cameras mileage system and method |
CN112184824A (en) * | 2019-07-05 | 2021-01-05 | 杭州海康机器人技术有限公司 | Camera external parameter calibration method and device |
WO2021056503A1 (en) * | 2019-09-29 | 2021-04-01 | 深圳市大疆创新科技有限公司 | Positioning method and apparatus for movable platform, movable platform, and storage medium |
CN113252066A (en) * | 2020-02-13 | 2021-08-13 | 纳恩博(北京)科技有限公司 | Method and device for calibrating parameters of odometer equipment, storage medium and electronic device |
CN113252066B (en) * | 2020-02-13 | 2024-04-09 | 纳恩博(北京)科技有限公司 | Calibration method and device for parameters of odometer equipment, storage medium and electronic device |
WO2021237520A1 (en) * | 2020-05-27 | 2021-12-02 | 华为技术有限公司 | Method and apparatus for calibrating extrinsics, and device and storage medium |
CN112654886A (en) * | 2020-05-27 | 2021-04-13 | 华为技术有限公司 | External parameter calibration method, device, equipment and storage medium |
CN112082491A (en) * | 2020-09-11 | 2020-12-15 | 苏州杰锐思智能科技股份有限公司 | Height detection method based on point cloud |
CN112567439A (en) * | 2020-11-09 | 2021-03-26 | 驭势(上海)汽车科技有限公司 | Method and device for determining traffic flow information, electronic equipment and storage medium |
CN114538027A (en) * | 2020-11-26 | 2022-05-27 | 合肥欣奕华智能机器股份有限公司 | Full-automatic visual positioning transfer equipment and control method thereof |
CN112258590A (en) * | 2020-12-08 | 2021-01-22 | 杭州迦智科技有限公司 | Laser-based depth camera external parameter calibration method, device and storage medium thereof |
CN112884842B (en) * | 2020-12-25 | 2024-03-22 | 杭州迦智科技有限公司 | External parameter synchronous calibration method for double rudder zero position and acquisition device and movable platform |
CN112884842A (en) * | 2020-12-25 | 2021-06-01 | 杭州迦智科技有限公司 | Synchronous calibration method for zero position of double rudder wheels and external parameters of acquisition device and movable platform |
CN113063353A (en) * | 2021-03-31 | 2021-07-02 | 深圳中科飞测科技股份有限公司 | Coordinate system establishing method, detection device, detection equipment and storage medium |
CN113063353B (en) * | 2021-03-31 | 2023-05-09 | 深圳中科飞测科技股份有限公司 | Coordinate system establishment method, detection device, detection equipment and storage medium |
CN113538598A (en) * | 2021-07-21 | 2021-10-22 | 北京能创科技有限公司 | Active Stereo Vision System Calibration Method |
CN113538598B (en) * | 2021-07-21 | 2022-03-25 | 北京能创科技有限公司 | Active stereo vision system calibration method |
CN113814987A (en) * | 2021-11-24 | 2021-12-21 | 季华实验室 | Multi-camera robot hand-eye calibration method and device, electronic equipment and storage medium |
CN113814987B (en) * | 2021-11-24 | 2022-06-03 | 季华实验室 | Multi-camera robot hand-eye calibration method, device, electronic device and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109242912A (en) | Join scaling method, electronic equipment, storage medium outside acquisition device | |
CN104596502B (en) | Object posture measuring method based on CAD model and monocular vision | |
CN112183171B (en) | Method and device for building beacon map based on visual beacon | |
CN109029257A (en) | Based on stereoscopic vision and the large-scale workpiece pose measurement system of structure light vision, method | |
CN111089569A (en) | A large-scale box measurement method based on monocular vision | |
CN102376089A (en) | Target correction method and system | |
CN105444696B (en) | A kind of binocular ranging method and its application based on perspective projection line measurement model | |
CN107941217A (en) | A kind of robot localization method, electronic equipment, storage medium, device | |
CN106990776A (en) | Robot goes home localization method and system | |
CN110136211A (en) | A workpiece positioning method and system based on active binocular vision technology | |
CN101901502B (en) | A global optimal registration method for multi-view point cloud data in optical 3D measurement | |
CN114310901B (en) | Coordinate system calibration method, device, system and medium for robot | |
CN113870366B (en) | Calibration method and calibration system of three-dimensional scanning system based on pose sensor | |
CN112958960A (en) | Robot hand-eye calibration device based on optical target | |
CN104167001B (en) | Large-visual-field camera calibration method based on orthogonal compensation | |
CN113393413B (en) | Water area measurement method and system based on monocular and binocular vision synergy | |
CN110992416A (en) | High-reflection-surface metal part pose measurement method based on binocular vision and CAD model | |
CN104123726A (en) | Blanking point based large forging measurement system calibration method | |
CN110363801A (en) | The corresponding point matching method of workpiece material object and workpiece three-dimensional CAD model | |
CN111145267B (en) | 360-degree panoramic view multi-camera calibration method based on IMU assistance | |
CN118493394A (en) | Method for actively sensing target distance by monocular vision robot | |
CN112700505A (en) | Binocular three-dimensional tracking-based hand-eye calibration method, equipment and storage medium | |
CN113155152A (en) | Camera and inertial sensor spatial relationship self-calibration method based on lie group filtering | |
CN104751452A (en) | Monocular camera calibration method based on any known movement | |
Van Toan et al. | A single 2-D LiDAR extrinsic calibration for autonomous mobile robots |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190118 |