CN107817003B - External parameter calibration method of distributed large-size space positioning system - Google Patents

External parameter calibration method of distributed large-size space positioning system Download PDF

Info

Publication number
CN107817003B
CN107817003B CN201610822085.0A CN201610822085A CN107817003B CN 107817003 B CN107817003 B CN 107817003B CN 201610822085 A CN201610822085 A CN 201610822085A CN 107817003 B CN107817003 B CN 107817003B
Authority
CN
China
Prior art keywords
calibration
matrix
point
transmitter
projection
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610822085.0A
Other languages
Chinese (zh)
Other versions
CN107817003A (en
Inventor
钟征宇
刘青
贾昆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Aerocomm Technology Co ltd
Original Assignee
Xi'an Aerocomm Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xi'an Aerocomm Technology Co ltd filed Critical Xi'an Aerocomm Technology Co ltd
Priority to CN201610822085.0A priority Critical patent/CN107817003B/en
Publication of CN107817003A publication Critical patent/CN107817003A/en
Application granted granted Critical
Publication of CN107817003B publication Critical patent/CN107817003B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C25/00Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass

Landscapes

  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

The invention relates to an external parameter calibration method of a distributed large-size space positioning system, which solves the difficult problems of complicated calibration process of a large-size space measurement system, sensitivity of system global optimization to an external parameter initial value, high cost and the like, and provides a method for realizing the initial calibration of the system by using a two-dimensional calibration rod by assuming that a virtual projection plane exists in front of a transmitter, and then realizing the global calibration of the system by using a Levenberg-Marquardt iterative optimization and a method for minimizing a target function. The method has the advantages of simple operation process and high calibration speed, and can effectively reduce the system cost while improving the system measurement precision and the working efficiency.

Description

External parameter calibration method of distributed large-size space positioning system
Technical Field
The invention relates to an external parameter calibration method for a distributed large-size space positioning system. The principle is as follows: firstly, a virtual projection plane is assumed to exist in front of a transmitter, a virtual transmitter perspective projection model is established, initial calibration of each calibration unit is achieved, and then overall calibration of the system is achieved through a method of Levenberg-Marquardt iterative optimization and objective function minimization.
Background
With the increasing requirements of assembly of large-scale product parts such as aviation, aerospace, ships, automobiles and the like and the requirements of large parts on precision positioning of butt joint pairing and real-time pose measurement and control, large-size space digital measurement systems and application technologies thereof are widely concerned by the industry and academia. At present, the mature large-size measurement technology applied at home and abroad mainly comprises a laser tracker measurement system, a large-view-field vision measurement system, a theodolite measurement system and the like.
The distributed large-size space positioning system is mainly composed of a plurality of measuring base stations, can monitor all parts of a large measured object at the same time, is high in real-time performance and large in measuring range, coordinates contradiction between measuring accuracy and large-size space by means of adding the base stations, and is wide in application prospect. Before measurement, the position relationship between the system world coordinate system and each measurement base station needs to be established, that is, the system external parameter calibration is performed, but the existing external calibration method usually uses other auxiliary measurement devices or needs to add additional constraint conditions, so that additional error sources are easily added, the operation process is complicated, the cost is high, and the field application is not facilitated.
Disclosure of Invention
Aiming at the defects of the prior art, the invention aims to provide an external parameter calibration method of a distributed large-size space positioning system, which realizes the quick and accurate calibration of the system, thereby realizing the high-precision and omnibearing large-size space measurement.
The technical scheme adopted by the invention for realizing the purpose is as follows:
1) two near-infrared lasers are installed inside a transmitter rotating platform, the transmitter rotating platform rotates anticlockwise at a constant speed around a rotating shaft under the driving of a servo motor, two fan-shaped laser planes with fixed inclination angles and an LED synchronous optical signal are sent to a measuring space in real time, and the LED synchronous optical signal is used as a zero position signal of a laser plane rotating period;
2) reasonably arranging transmitter networks according to the measurement space, setting fixed rotating speed of each transmitter respectively, sequentially increasing 100r/min for each transmitter by taking 2500r/min as a starting point, and numbering the transmitters in an anticlockwise direction;
3) placing a calibration rod, and establishing a coordinate system of each transmitter, wherein the specific method comprises the steps of taking the intersection point of two laser planes on a rotating shaft as the origin of the coordinate system, taking the intersection line of the laser plane 1 and a vertical Y-axis plane when the laser plane sweeps across a receiver on the left side of the calibration rod as the Z axis, taking the direction pointing to the receiver on the left side as the Z axis forward direction, and determining the X axis direction through a right-hand rule;
4) every two transmitters are used as a calibration unit, the calibration units are numbered, a transmitter coordinate system is used as a world coordinate system, calibration rod with fixed length is used for collecting calibration point data, and a calibration point set is established;
5) respectively calibrating the single calibration unit, and calculating the coordinates of the projection points of the calibration points on a virtual projection plane on the assumption that the virtual projection plane exists at a specific position in the forward direction perpendicular to the Z axis of the transmitter, wherein the model is as follows:
Figure GDA0003070101540000021
wherein P ═ (x y z): the coordinate of a projection point P of the calibration point on the virtual projection plane;
az: the horizontal angle of the calibration point under the transmitter coordinate system;
pa: the pitch angle of the calibration point under the transmitter coordinate system;
d: assuming that a virtual projection plane exists at z ═ d;
6) establishing a virtual transmitter perspective projection model, calculating the coordinates of the projection points of the calibration points on a virtual projection plane, and performing coordinate transformation of a projection point set through translation transformation and scaling transformation, wherein the model comprises the following steps:
p ═ MP and M ═ STS
Wherein, P: according to the projection point coordinates on the virtual projection plane obtained in the step 5;
p': projecting point coordinates on the virtual projection plane after coordinate transformation;
m: a coordinate transformation matrix of the point set;
s: the coordinate scaling transformation matrix of the set of points,
Figure GDA0003070101540000031
d0the average distance from each point in the point set to the image origin;
TS: the coordinate translation transformation matrix of the set of points,
Figure GDA0003070101540000032
(Tx Ty)Tis the centroid coordinate of the point set;
7) calculating an essential matrix, and decomposing the essential matrix to obtain possible solutions of a rotation matrix and a unit translation vector, wherein the method comprises the following specific steps:
I. calculating a coefficient matrix according to the projection point coordinates obtained in the step 6 after the coordinate transformation of the point set, taking the calibration unit one as an example, and using the model as follows:
Am=(a1 a2 … aN)T
wherein Am: coefficient matrix, aiIs the element of the ith in the coefficient matrix Am, and N is the number of the index points;
Figure GDA0003070101540000033
wherein d in the step 5 is 1;
Figure GDA0003070101540000034
the projection point coordinates of the ith calibration point in the first transmitter coordinate system;
Figure GDA0003070101540000035
the projection point coordinates of the ith calibration point in the second transmitter coordinate system;
decomposing the coefficient matrix Am by using singular values to obtain an essential matrix E;
decompose the essence matrix, the model is as follows:
E=USVT
R1=UZVTor R2=UZTVT
Figure GDA0003070101540000036
t1=V(0 0 1)TOr t2=-V(0 0 1)T
Determining possible solutions for the rotation matrix and unit translation vector, the model being as follows:
(R1|t1)
(R1|t2)
(R2|t1)
(R2|t2)
8) and (3) rejecting the pseudo solution in the step 7 by utilizing physical screening and Sampson distance minimum constraint to obtain a unique solution of a rotation matrix and a unit translation vector, wherein the specific steps are as follows:
I. physical screening: ensuring that the index point is in front of both transmitters, i.e. Z1> 0 and Z2>0;
Sampson distance minimum constraint: ensuring that, in the virtual perspective projection model, the projection point Pi2To a straight line L1Distance of (2) and projection point Pi1To the outer polar line L2The sum of the squared distances of (a) is minimal, the model is as follows:
Figure GDA0003070101540000041
9) and (3) determining a scale factor according to the unique solution of the rotation matrix and the unit translation vector obtained by screening in the step 8 and taking the length of the calibration rod as a constraint to obtain a translation matrix, wherein the model is as follows:
Figure GDA0003070101540000042
T=kt
wherein, k: a scale factor;
l: calibrating the theoretical length of the rod;
l': calculating the length of the calibration rod according to the rotation matrix and the unit translation vector in the step 8;
t: the unit translation vector obtained in the step 8;
t: translating the matrix;
10) and (3) sequentially calculating a rotation matrix and a translation matrix of each calibration unit in the steps 4-9, taking the rotation matrix and the translation matrix as initial values, and minimizing an objective function by using a Levenberg-Marquardt iterative optimization method to obtain a final rotation matrix and a final translation matrix, wherein the model is as follows:
Figure GDA0003070101540000051
wherein N is the number of the calibration points;
j is the number of times the calibration rod collects data;
l'jis the calculated value of the jth calibration rod length;
l is the theoretical length of the calibration rod.
The method avoids the use of other auxiliary measuring equipment, firstly assumes that a virtual projection plane exists in front of the transmitter, establishes a virtual transmitter perspective projection model, realizes initial calibration of the external parameters by a linear method, then adopts a Levenberg-Marquardt iterative optimization method, minimizes an objective function, realizes global calibration by a nonlinear method, and solves the problems of complicated calibration process, sensitivity of nonlinear optimization to the initial values of the external parameters and the like. The invention ensures that the solving parameters are globally optimal, the solving speed is high, and the optimization result is stable, thereby ensuring the calibration precision of the system parameters.
Drawings
FIG. 1 is a schematic system flow diagram;
FIG. 2 is a schematic diagram of a transmitter architecture;
FIG. 3 is a transmitter numbering diagram;
FIG. 4 is a schematic view of a calibration rod;
FIG. 5 is a schematic diagram of establishing a transmitter coordinate system;
fig. 6 is a schematic view of a virtual transmitter perspective projection model.
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings.
An external parameter calibration method of a distributed large-size space positioning system comprises the following steps:
1) as shown in fig. 2, two near-infrared lasers are installed inside a transmitter rotating platform, the transmitter rotating platform rotates counterclockwise around a rotating shaft at a constant speed under the driving of a servo motor, and sends two fan-shaped laser planes with fixed inclination angles and an LED synchronous optical signal to a measurement space in real time, and the LED synchronous optical signal is used as a zero position signal of a laser plane rotating period;
2) as shown in fig. 4, taking four transmitters as an example, a network of transmitters is reasonably arranged according to a measurement space, each transmitter is respectively set with a fixed rotation speed of 2500r/min, 2600r/min, 2700r/min and 2800r/min, and the transmitters are numbered in a counterclockwise direction, which are respectively a first transmitter, a second transmitter, a third transmitter and a fourth transmitter;
3) as shown in fig. 5, a calibration rod is placed in front of a transmitter network, and a coordinate system of each transmitter is established, wherein the specific method is that the intersection point of two laser planes on a rotating shaft is taken as the origin of the coordinate system, the direction of the two laser planes is taken as the forward direction of a Y axis downwards along the rotating shaft of the transmitter, the intersection line of the laser plane 1 which sweeps through a receiver on the left side of the calibration rod and a vertical Y axis plane is taken as a Z axis, the direction pointing to the receiver on the left side is taken as the forward direction of the Z axis, and the direction of the X axis is determined by;
4) numbering the calibration units, taking a transmitter coordinate system as a world coordinate system, taking every two transmitters as a calibration unit, and respectively: the method comprises the following steps that a first transmitter and a second transmitter are calibration units 1, a first transmitter and a third transmitter are calibration units 2, a first transmitter and a fourth transmitter are calibration units 3, calibration rod with fixed lengths are used for collecting calibration point data, and a calibration point set is established;
5) respectively calibrating the single calibration unit, and calculating the coordinates of the projection points of the calibration points on a virtual projection plane on the assumption that the virtual projection plane exists at a specific position in the forward direction perpendicular to the Z axis of the transmitter, wherein the model is as follows:
Figure GDA0003070101540000061
wherein P ═ (x y z): the coordinate of a projection point P of the calibration point on the virtual projection plane;
az: the horizontal angle of the calibration point under the transmitter coordinate system;
pa: the pitch angle of the calibration point under the transmitter coordinate system;
d: assuming that a virtual projection plane exists at z ═ d;
6) as shown in fig. 6, taking the calibration unit 1 as an example, a virtual transmitter perspective projection model is established, coordinates of a projection point of the calibration point on a virtual projection plane are calculated, and coordinate transformation of a projection point set is performed through translation transformation and scaling transformation, and the model is as follows:
p ═ MP and M ═ STS
Wherein, P: according to the projection point coordinates on the virtual projection plane obtained in the step 5;
p': projecting point coordinates on the virtual projection plane after coordinate transformation;
m: a coordinate transformation matrix of the point set;
s: the coordinate scaling transformation matrix of the set of points,
Figure GDA0003070101540000071
d0the average distance from each point in the point set to the image origin;
TS: the coordinate translation transformation matrix of the set of points,
Figure GDA0003070101540000072
(Tx Ty)Tis the centroid coordinate of the point set;
7) calculating an essential matrix, and decomposing the essential matrix to obtain possible solutions of a rotation matrix and a unit translation vector, wherein the method comprises the following specific steps:
I. calculating a coefficient matrix according to the projection point coordinates obtained in the step 6 after the coordinate transformation of the point set, taking the calibration unit one as an example, and using the model as follows:
Am=(a1 a2 … aN)T
wherein Am: coefficient matrix, aiIs the element of the ith in the coefficient matrix Am, and N is the number of the index points;
Figure GDA0003070101540000081
wherein d in the step 5 is 1;
Figure GDA0003070101540000082
the projection point coordinates of the ith calibration point in the first transmitter coordinate system;
Figure GDA0003070101540000083
the ith index point is in number twoProjection point coordinates under a transmitter coordinate system;
decomposing the coefficient matrix Am by using singular values to obtain an essential matrix E;
decompose the essence matrix, the model is as follows:
E=USVT
R1=UZVTor R2=UZTVT
Figure GDA0003070101540000084
t1=V(0 0 1)TOr t2=-V(0 0 1)T
Determining possible solutions for the rotation matrix and unit translation vector, the model being as follows:
(R1|t1)
(R1|t2)
(R2|t1)
(R2|t2)
8) and (3) rejecting the pseudo solution in the step 7 by utilizing physical screening and Sampson distance minimum constraint to obtain a unique solution of a rotation matrix and a unit translation vector, wherein the specific steps are as follows:
I. physical screening: ensuring that the index point is in front of both transmitters, i.e. Z1> 0 and Z2>0;
Sampson distance minimum constraint: ensuring that, in the virtual perspective projection model, the projection point Pi2To a straight line L1Distance of (2) and projection point Pi1To the outer polar line L2The sum of the squared distances of (a) is minimal, the model is as follows:
Figure GDA0003070101540000091
9) and (3) determining a scale factor according to the unique solution of the rotation matrix and the unit translation vector obtained by screening in the step 8 and taking the length of the calibration rod as a constraint to obtain a translation matrix, wherein the model is as follows:
Figure GDA0003070101540000092
T=kt
wherein, k: a scale factor;
l: calibrating the theoretical length of the rod;
l': calculating the length of the calibration rod according to the rotation matrix and the unit translation vector in the step 8;
t: the unit translation vector obtained in the step 8;
t: translating the matrix;
10) and (3) sequentially calculating a rotation matrix and a translation matrix of each calibration unit in the steps 4-9, taking the rotation matrix and the translation matrix as initial values, and minimizing an objective function by using a Levenberg-Marquardt iterative optimization method to obtain a final rotation matrix and a final translation matrix, wherein the model is as follows:
Figure GDA0003070101540000093
wherein N is the number of the calibration points;
j is the number of times the calibration rod collects data;
l'jis the calculated value of the jth calibration rod length;
l is the theoretical length of the calibration rod.

Claims (1)

1. An external parameter calibration method of a distributed large-size space positioning system is characterized by comprising the following steps:
1) two near-infrared lasers are installed inside a transmitter rotating platform, the transmitter rotating platform rotates anticlockwise at a constant speed around a rotating shaft under the driving of a servo motor, two fan-shaped laser planes with fixed inclination angles and an LED synchronous optical signal are sent to a measuring space in real time, and the LED synchronous optical signal is used as a zero position signal of a laser plane rotating period;
2) reasonably arranging transmitter networks according to the measurement space, setting fixed rotating speed of each transmitter respectively, sequentially increasing 100r/min for each transmitter by taking 2500r/min as a starting point, and numbering the transmitters in an anticlockwise direction;
3) placing a calibration rod, and establishing a coordinate system of each transmitter, wherein the specific method comprises the steps of taking the intersection point of two laser planes on a rotating shaft as the origin of the coordinate system, taking the intersection line of the laser plane 1 and a vertical Y-axis plane when the laser plane sweeps across a receiver on the left side of the calibration rod as the Z axis, taking the direction pointing to the receiver on the left side as the Z axis forward direction, and determining the X axis direction through a right-hand rule;
4) every two transmitters are used as a calibration unit, the calibration units are numbered, a transmitter coordinate system is used as a world coordinate system, calibration rod with fixed length is used for collecting calibration point data, and a calibration point set is established;
5) respectively calibrating the single calibration unit, and calculating the coordinates of the projection points of the calibration points on a virtual projection plane on the assumption that the virtual projection plane exists at a specific position in the forward direction perpendicular to the Z axis of the transmitter, wherein the model is as follows:
Figure FDA0003070101530000011
wherein P ═ (x y z): the coordinate of a projection point P of the calibration point on the virtual projection plane;
az: the horizontal angle of the calibration point under the transmitter coordinate system;
pa: the pitch angle of the calibration point under the transmitter coordinate system;
d: assuming that a virtual projection plane exists at z ═ d;
6) establishing a virtual transmitter perspective projection model, calculating the coordinates of the projection points of the calibration points on a virtual projection plane, and performing coordinate transformation of a projection point set through translation transformation and scaling transformation, wherein the model comprises the following steps:
p ═ MP and M ═ STS
Wherein, P: according to the projection point coordinates on the virtual projection plane obtained in the step 5;
p': projecting point coordinates on the virtual projection plane after coordinate transformation;
m: a coordinate transformation matrix of the point set;
s: the coordinate scaling transformation matrix of the set of points,
Figure FDA0003070101530000021
d0the average distance from each point in the point set to the image origin;
TS: the coordinate translation transformation matrix of the set of points,
Figure FDA0003070101530000022
(Tx Ty)Tis the centroid coordinate of the point set;
7) calculating an essential matrix, and decomposing the essential matrix to obtain possible solutions of a rotation matrix and a unit translation vector, wherein the method comprises the following specific steps:
I. calculating a coefficient matrix according to the projection point coordinates obtained in the step 6 after the coordinate transformation of the point set, taking the calibration unit one as an example, and using the model as follows:
Am=(a1 a2 … aN)T
wherein Am: coefficient matrix, aiIs the element of the ith in the coefficient matrix Am, and N is the number of the index points;
Figure FDA0003070101530000023
wherein d in the step 5 is 1;
Figure FDA0003070101530000024
the projection point coordinates of the ith calibration point in the first transmitter coordinate system;
Figure FDA0003070101530000025
the ith index point is sent in the second numberProjection point coordinates under a projector coordinate system;
decomposing the coefficient matrix Am by using singular values to obtain an essential matrix E;
decompose the essence matrix, the model is as follows:
E=USVT
R1=UZVTor R2=UZTVT
Figure FDA0003070101530000026
t1=V(001)TOr t2=-V(001)T
Determining possible solutions for the rotation matrix and unit translation vector, the model being as follows:
(R1|t1)
(R1|t2)
(R2|t1)
(R2|t2)
8) and (3) rejecting the pseudo solution in the step 7 by utilizing physical screening and Sampson distance minimum constraint to obtain a unique solution of a rotation matrix and a unit translation vector, wherein the specific steps are as follows:
I. physical screening: ensuring that the index point is in front of both transmitters, i.e. Z1> 0 and Z2>0;
Sampson distance minimum constraint: ensuring that, in the virtual perspective projection model, the projection point Pi2To a straight line L1Distance of (2) and projection point Pi1To the outer polar line L2The sum of the squared distances of (a) is minimal, the model is as follows:
Figure FDA0003070101530000031
9) and (3) determining a scale factor according to the unique solution of the rotation matrix and the unit translation vector obtained by screening in the step 8 and taking the length of the calibration rod as a constraint to obtain a translation matrix, wherein the model is as follows:
Figure FDA0003070101530000032
T=kt
wherein, k: a scale factor;
l: calibrating the theoretical length of the rod;
l': calculating the length of the calibration rod according to the rotation matrix and the unit translation vector in the step 8;
t: the unit translation vector obtained in the step 8;
t: translating the matrix;
10) and (3) sequentially calculating a rotation matrix and a translation matrix of each calibration unit in the steps 4-9, taking the rotation matrix and the translation matrix as initial values, and minimizing an objective function by using a Levenberg-Marquardt iterative optimization method to obtain a final rotation matrix and a final translation matrix, wherein the model is as follows:
Figure FDA0003070101530000041
wherein N is the number of the calibration points;
j is the number of times the calibration rod collects data;
l'jis the calculated value of the jth calibration rod length;
l is the theoretical length of the calibration rod.
CN201610822085.0A 2016-09-14 2016-09-14 External parameter calibration method of distributed large-size space positioning system Active CN107817003B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610822085.0A CN107817003B (en) 2016-09-14 2016-09-14 External parameter calibration method of distributed large-size space positioning system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610822085.0A CN107817003B (en) 2016-09-14 2016-09-14 External parameter calibration method of distributed large-size space positioning system

Publications (2)

Publication Number Publication Date
CN107817003A CN107817003A (en) 2018-03-20
CN107817003B true CN107817003B (en) 2021-07-06

Family

ID=61600603

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610822085.0A Active CN107817003B (en) 2016-09-14 2016-09-14 External parameter calibration method of distributed large-size space positioning system

Country Status (1)

Country Link
CN (1) CN107817003B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109884658B (en) * 2019-03-04 2020-11-27 北京工业大学 Laser tracker station position method based on laser tracker multi-station position measuring system
CN110006459B (en) * 2019-04-30 2020-11-10 西安交通大学 Self-correcting system and method for R-LATs measurement network in workshop dynamic operation environment
CN110654571B (en) * 2019-11-01 2023-10-20 西安航通测控技术有限责任公司 Nondestructive testing robot system and method for surface defects of aircraft skin
CN111622114B (en) * 2020-05-08 2021-08-10 中铁大桥科学研究院有限公司 Bridge segment prefabrication construction line shape error adjusting method
CN113158846B (en) * 2021-04-06 2024-02-20 中航航空电子有限公司 Anti-interference sight tracking method for head-mounted waveguide display device

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20080091891A (en) * 2007-04-10 2008-10-15 삼성중공업 주식회사 The automatic calibration method in robot based multi-laser vision system
CN102252663A (en) * 2011-04-27 2011-11-23 西安交通大学 Field calibration method for local area space positioning system
CN105320271A (en) * 2014-07-10 2016-02-10 精工爱普生株式会社 HMD calibration with direct geometric modeling
CN105353999A (en) * 2014-05-27 2016-02-24 空中客车集团有限公司 Method for projecting virtual data and device enabling said projection

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1118685C (en) * 2001-10-15 2003-08-20 天津大学 In-situ calibration device for large-size 3D space measurement and its measuring method
EP1724609A1 (en) * 2005-05-18 2006-11-22 Leica Geosystems AG Method of determining postion of receiver unit
CN101307489A (en) * 2008-02-02 2008-11-19 中国科学院物理研究所 Halogen-containing high temperature superconductor crystal and method for making same
CN100575865C (en) * 2008-07-18 2009-12-30 西安交通大学 Space-location method based on double-rotating laser plane transmitter network
CN201757972U (en) * 2010-04-15 2011-03-09 沈阳瑞克上成广告传媒有限公司 Holographic suspension virtual image system
JP2015506121A (en) * 2011-11-16 2015-02-26 ストーウェル・リンネル・エルエルシー System and method for 3D projection mapping with robot controlled objects
CZ304495B6 (en) * 2012-12-12 2014-05-28 ÄŚVUT v Praze, Fakulta strojnĂ­ Device for optical measurement and/or optical calibration of a body position within a space
EP2787323A1 (en) * 2013-04-05 2014-10-08 Leica Geosystems AG Measuring device with function for calibrating a display image position of an electronic reticule
US9189711B2 (en) * 2013-10-28 2015-11-17 Eastman Kodak Company Method for aligning imaging systems
TW201528048A (en) * 2014-01-03 2015-07-16 Egismos Technology Corp Image-based virtual interactive device and method thereof

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20080091891A (en) * 2007-04-10 2008-10-15 삼성중공업 주식회사 The automatic calibration method in robot based multi-laser vision system
CN102252663A (en) * 2011-04-27 2011-11-23 西安交通大学 Field calibration method for local area space positioning system
CN105353999A (en) * 2014-05-27 2016-02-24 空中客车集团有限公司 Method for projecting virtual data and device enabling said projection
CN105320271A (en) * 2014-07-10 2016-02-10 精工爱普生株式会社 HMD calibration with direct geometric modeling

Also Published As

Publication number Publication date
CN107817003A (en) 2018-03-20

Similar Documents

Publication Publication Date Title
CN107817003B (en) External parameter calibration method of distributed large-size space positioning system
CN106226780B (en) More rotor-wing indoor positioning systems and implementation method based on scanning laser radar
CN103591891B (en) Fine control field precision source-tracing method for indoor space measuring and locating system
CN110470304B (en) High-precision target positioning and speed measuring method based on unmanned aerial vehicle photoelectric platform
CN112325796A (en) Large-scale workpiece profile measuring method based on auxiliary positioning multi-view point cloud splicing
CN112525164A (en) Method for detecting deformation of super high-rise building based on unmanned aerial vehicle oblique photography technology
CN104515478A (en) Automatic three-dimensional measuring method and automatic three-dimensional measuring system for high-precision blade of aviation engine
CN111811483B (en) Multi-camera networking positioning calibration method for airplane digital assembly
CN109238247A (en) A kind of six degree of freedom measurement method towards the complicated scene of large space
CN108917723A (en) A kind of pose on-line measurement system and method for the docking of cylindrical bay section
CN104880215B (en) The method of multisensor platform camera laser integration calibration
CN104422425A (en) Irregular-outline object space attitude dynamic measuring method
CN107727118B (en) Method for calibrating GNC subsystem equipment attitude measurement system in large aircraft
CN111208497A (en) Airborne laser radar system adjustment processing method
CN207456381U (en) Improve the device of laser tracker measurement accuracy
CN110871824B (en) Method and system for monitoring surrounding environment of track
CN106228564A (en) The outer parameter two step associating online calibration method of many mesh camera and system
CN105953755A (en) Networking large-size space positioning system on-site calibration method
CN106908078B (en) On-site calibration method of networked large-size space positioning system
CN115792800A (en) Grid search-based double-station three-dimensional cross positioning method
CN113177987B (en) Outfield global calibration method and system of visual tracking measurement system
Wu et al. A novel precise guiding method for visual guiding theodolite measurement in volume space
CN115183981A (en) Wind tunnel coordinate shafting calibration method and calibration device based on coding marks
CN110148218A (en) A kind of method of high-volume airborne lidar point cloud data global optimization
CN106908018A (en) A kind of intrinsic parameter scaling method of double-rotating laser plane transmitter

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant