CN113607188A - Calibration system and method of multi-view-field star sensor based on theodolite cross-hair imaging - Google Patents
Calibration system and method of multi-view-field star sensor based on theodolite cross-hair imaging Download PDFInfo
- Publication number
- CN113607188A CN113607188A CN202110883235.XA CN202110883235A CN113607188A CN 113607188 A CN113607188 A CN 113607188A CN 202110883235 A CN202110883235 A CN 202110883235A CN 113607188 A CN113607188 A CN 113607188A
- Authority
- CN
- China
- Prior art keywords
- theodolite
- star sensor
- view
- field
- calibration
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
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
A theodolite cross-hair imaging-based multi-view-field star sensor calibration system and method comprises the following steps: modeling the process of imaging a cross hair of the theodolite on the star sensor image surface through a telescope objective lens and a star sensor lens; establishing a calibration track according to the field parameters, and acquiring cross images shot by the star sensor at different calibration positions and theodolite measurement angles; extracting a cross image into a cross point image coordinate, taking the cross point coordinate and a theodolite angle as input data, and constructing an optimization problem according to a cross imaging model to solve the rotation relation between a theodolite coordinate system and a view field coordinate system; establishing a rotation relation between theodolite coordinate systems through theodolite mutual aiming; and fusing the step-by-step results to obtain a multi-view-field star sensor structure parameter calibration result. The method keeps high calibration precision, uses a plurality of electronic theodolites to construct a calibration system, has free system configuration, and is applied to a multi-view-field star sensor system with large size, large weight and any view-field distribution structure.
Description
Technical Field
The invention belongs to the technical field of astronomical navigation, and provides a calibration system and method of a multi-view-field star sensor based on theodolite cross hair imaging.
Background
The measurement precision and the dynamic performance are two important performance indexes of the star sensor. The single-view field star sensor is limited by the size of a view field, the attitude measurement error of the single-view field star sensor around the direction of an optical axis is one order of magnitude higher than that of the other two axes, and the multi-view field star sensor can compensate the low-precision measurement axes of other view fields by using the measurement result of the high-precision measurement axes of one view field, so that the three-axis equal-precision attitude measurement is realized; because the field range of the multi-field star sensor is far larger than that of a single-field star sensor with the same specification, the number of observation stars of the multi-field star sensor is larger than that of the single-field star sensor, and the multi-field star sensor has higher integral attitude measurement precision; in addition, when the angular speed of the star sensor is increased, the detection sensitivity of stars and the like of each field of view is reduced differently, the number of observed stars is reduced, and the multi-field star sensor can realize the information fusion of star points of each field of view, and can still output effective postures as long as the total number of observed stars is more than three, so that the multi-field star sensor has higher dynamic performance than a single-field star sensor.
The key step of measuring the attitude of the multi-view-field star sensor system is to fuse star vectors in different view fields into a unified coordinate system. Due to the existence of processing and adjusting errors, the rotation relation between the coordinate systems of the fields of view of the multi-field-of-view star sensor, namely the structural parameters, deviates from the design values, and therefore the attitude measurement is inaccurate. Therefore, the accurate calibration of the structural parameters has important significance for ensuring the attitude measurement accuracy of the multi-view-field star sensor.
At present, three methods are commonly used for calibrating the structural parameters of the multi-field star sensor. The first method is an external field calibration method, namely, different fields of view of the multi-field-of-view star sensor are controlled to be exposed simultaneously, the rotation relation of each field of view coordinate system relative to a celestial coordinate system is respectively calculated, and the structural parameters of each field of view are obtained by using the celestial coordinate system as a middle coordinate system. The method is easily influenced by atmospheric refraction, so that the calibration result precision is poor, and the measurement precision requirement of the multi-view-field star sensor cannot be met. The second method is to calibrate the structural parameters by combining an electronic theodolite with a reference cube mirror, and is disclosed as a method for calibrating the optical axis pointing direction of a star sensor probe assembly (application No. CN2011104609570, published No. CN 102538825A). The method is characterized in that the reference cube mirror is fixedly connected to each field of view of a multi-field star sensor, the rotating relation between a field of view coordinate system and a prism coordinate system can be established when the square elements inside and outside the star sensor of each field of view are calibrated, the rotating relation between each prism coordinate system and a theodolite coordinate system can be established by using a plurality of theodolite sighting prisms, the rotating relation between each theodolite coordinate system can be established by the mutual sighting of the theodolite, and finally, the field of view structural parameters are obtained by fusion. The calibration precision of the method is limited by the processing precision of the reference cube mirror to a great extent, the processing precision of the reference cube mirror can only reach the angle classification level, and the calibration precision can not meet the use requirement of the multi-field star sensor when the structural parameter precision is higher. The third method uses a three-axis turntable to calibrate the structural parameters, and is disclosed in the patent of ' a calibration system and a calibration method of a multi-view field star sensor based on a three-axis turntable ' (application number CN2018113663006, publication number CN109459059A) '. However, the mounting size and the load capacity of a three-axis turntable are limited, the method is only suitable for calibrating a light and small multi-view-field star sensor system, the visual axis distribution of the multi-view-field star sensor also needs to have symmetry, and the method cannot be applied to the multi-view-field star sensor system with large volume, large weight or asymmetric visual axis distribution. Furthermore, applications such as patent application, application No. CN2010101884887, publication No.: CN 101858755A discloses a calibration method of a star sensor; application No.: CN2010105127866, publication No.: the three-probe star sensor pointing calibration method disclosed in CN 102032918A; application No.: CN2014102253866, publication no: CN104154928A discloses a method for calibrating installation errors of a built-in star sensor suitable for an inertial platform; application No.: CN2014103608057, publication no: CN105318891B discloses a calibration device for the installation error of a star sensor reference cube mirror; application No.: CN2018115226969, publication No.: CN 109506645A discloses a star sensor mounting matrix ground precision measurement method; application No.: CN2019107203476, publication No.: CN 110345970A discloses a method and a device for calibrating an optical navigation sensor; application No.: CN2019113583595, publication No.: the calibration method between the star sensor measurement coordinate system and the reference cubic mirror coordinate system disclosed in CN 111044077A, and the like, all use theodolite to calibrate the conversion relationship between the star sensor relevant parameters and the relevant coordinate system, but these methods all need to be used in combination with the reference cubic mirror, some methods also need to add additional reference sources, such as collimator tubes, high-precision rotating tables or calibration fields formed by calibration points, and the like, the calibration system is complicated to build, and the reference cubic mirror needs to be used to build the intermediate conversion coordinate system, the related coordinate systems are numerous, calibration errors are difficult to control in multiple times of coordinate system conversion, and the insufficient processing precision of the reference cubic mirror can also cause the precision reduction of the calibration result.
Disclosure of Invention
The invention aims to: in the existing multi-view field star sensor structure parameter calibration method, both the external field calibration method and the theodolite combined prism calibration method cannot provide a calibration result with high enough precision, and the calibration method based on the three-axis turntable can provide a high-precision calibration result but has strict requirements on the overall dimension, the overall weight and the visual axis structure of the multi-view field star sensor.
Aiming at the problems, the method provides a multi-view field star sensor structure parameter calibration method based on electronic theodolite cross hair imaging, the method uses a plurality of electronic theodolites to respectively correspond to each view field of the multi-view field star sensor, a luminous cross hair in each theodolite can be imaged on an image surface of the corresponding view field through an optical system, the rotation relation between a theodolite coordinate system and a view field coordinate system can be calibrated by establishing an imaging model and collecting images, the rotation relation between the theodolite coordinate systems is established by mutually aiming through the theodolites, and finally, the view field structure parameters are obtained through fusion. The calibration process is mainly based on the optical imaging principle, higher calibration precision is guaranteed, meanwhile, the used electronic theodolite is free in installation mode, and the size, the weight and the configuration of the multi-view-field star sensor are not limited.
The technical scheme of the invention is as follows:
the invention discloses a calibration system of a multi-view-field star sensor based on theodolite cross hair imaging, which comprises a marble platform (1), a multi-view-field star sensor (2), an electronic theodolite (3) and data acquisition equipment (4); the marble platform (1) is used for bearing all other equipment; the multi-view-field star sensor (2) is arranged on the marble platform (1), and the visual axis of each view field is parallel to the plane of the marble platform; the electronic theodolite (3) is placed in front of each view field of the multi-view-field star sensor, the optical axis of a telescope of the electronic theodolite is superposed with the view axes of the view fields, and the base of the electronic theodolite is fixed and then the support foot leveling theodolite is adjusted; the data acquisition equipment (4) is connected with the multi-view-field star sensor (2) and the theodolite (3) and is used for acquiring images of the multi-view-field star sensor (2) and measurement data of the theodolite (3).
The invention discloses a calibration method of a multi-view-field star sensor based on theodolite cross hair imaging, which adopts the calibration system and comprises the following steps:
modeling a process of imaging a cross hair of a theodolite on an image surface of a star sensor through a telescope objective and a star sensor lens;
establishing a calibration track, and collecting cross images shot by the star sensor at different calibration positions and theodolite measurement angles;
extracting a cross image into a cross point image coordinate by using an image processing method, preprocessing the acquired data, and constructing an optimization problem according to a cross imaging model to solve the rotation relation between a theodolite coordinate system and a view field coordinate system;
establishing a rotation relation between theodolite coordinate systems through theodolite mutual aiming;
and step five, fusing the stepping results to obtain a multi-view-field star sensor structure parameter calibration result.
Further, the step one further comprises the following steps: when the theodolite telescope is stable at a certain position, the luminous cross hairs clearly image on an image surface through the telescope objective lens and the star sensor lens; the telescope optical axis direction vector under the view field coordinate system is determined by the horizontal angle and the vertical angle of the theodolite and the installation angle of the theodolite coordinate system and the view field coordinate system; the vector of the principal ray direction of the cross-hair cross point under the field-of-view coordinate system is determined by the image coordinate of the cross point and the star sensor imaging model; the two vectors are parallel to each other and in opposite directions.
Further, the second step further includes the following steps: firstly, adjusting the brightness of cross hairs and the focusing position of a telescope, and establishing a calibration track according to the field condition of a star sensor; and sequentially rotating the theodolite to each calibration position in the calibration track, collecting the star sensor images for multiple times and measuring the horizontal angle and the vertical angle of the theodolite for multiple times.
Further, the third step further includes the following steps: averaging the horizontal angle and the vertical angle of the theodolite collected at each calibration position to obtain the accurate theodolite angle of the calibration position; extracting sub-pixel level coordinates of the intersection points of the collected star sensor images by using an image processing method, and averaging the coordinates of the intersection points belonging to the same calibration position to obtain accurate image coordinates of the intersection points of the calibration position; according to a cross-hair imaging model of the theodolite, using a telescope optical axis direction vector vSAnd cross point chief ray direction vector wSThe difference is a residual function, and a minimum optimization problem is constructed for the rotational Euler angle theta of the theodolite coordinate system and the single-view field coordinate system1,θ2,θ3Solving, and constructing a minimal value optimization problem as follows:
where m is the total number of calibration positions, #i＝(Hzi,Vi)TTheodolite angle, x, for the ith calibration positioni＝(xp,i,yp,i)TFor the coordinates of the intersection image at the ith calibration position, p ═ f, xc,yc,p1,p2,q1,q2)TIs the internal parameter of the star sensor, including the focal length f and the principal point (x)c,yc) Radial distortion factor (p)1,p2) And tangential distortion coefficient (q)1,q2) The residual function is expressed as
Wherein the content of the first and second substances,
and
the vectors v of the optical axis directions of the telescopes respectively corresponding to the ith calibration positionS,iAnd ith calibration position cross point chief ray direction vector wS,iScaling to z is a vector after 1 plane. Solving the obtained rotation Euler angle theta1,θ2,θ3Converted into a direction cosine matrix form and recorded as a view field number k
Further, the fourth step further includes the following steps: aiming cross hairs of the two theodolites at each other and recording the measurement angle to obtain the rotation relation of the coordinate systems of the two theodolites; mutually aiming all the theodolites in pairs to obtain the rotation relation between coordinate systems of all the theodolites, and recording the direction cosine matrix from the coordinate system of the theodolite No. k to the coordinate system of the theodolite No. l as
Further, the fifth step further includes the following steps: and fusing all the obtained rotation relations to obtain a multi-view-field star sensor structure parameter calibration result, wherein a direction cosine matrix from a k view field coordinate system to a l view field coordinate system is as follows:
the invention also applies the calibration method of the multi-view field star sensor based on the cross-hair imaging of the theodolite to the astronomical navigation device.
The technical scheme of the invention can realize the following beneficial technical effects:
the calibration method calibrates the rotation relation between the theodolite coordinate system and the view field coordinate system by using a cross-hair imaging mode, and the calibration precision of the optical imaging method is usually in the order of angular seconds and is far higher than that of a reference cube mirror method.
The calibration method is characterized in that the multi-view-field star sensor to be calibrated is arranged on a stable base, a calibration system is constructed by a plurality of electronic theodolites around the multi-view-field star sensor, a high-precision rotary table is not needed, the limitation of the high-precision rotary table on the size and the weight of a bearing object is removed, and the method can be applied to the multi-view-field star sensor system with large size and large weight.
The calibration method uses the electronic theodolite and a single star sensor to form the sub-calibration systems, the sub-calibration systems are independent from each other, the special properties such as symmetry and the like of the field distribution of the multi-field star sensor are not required, and the method can be applied to the multi-field star sensor system with any structure and any number of fields.
Drawings
FIG. 1 is a schematic structural diagram of a calibration system of a multi-field-of-view star sensor based on cross hair imaging of an electronic theodolite according to the present invention;
FIG. 2 is a theodolite telescope optical axis direction vector model;
FIG. 3 is an imaging model of a star sensor;
FIG. 4 is a cross hair image real shot local picture;
fig. 5 shows calibration trajectories of the theodolite-single-view sub-calibration system.
Detailed Description
As shown in fig. 1, a calibration system of a theodolite cross-hair imaging-based multi-view-field star sensor comprises a marble platform 1 for providing stable support, a multi-view-field star sensor 2 to be calibrated, a plurality of electronic theodolites 3 with auto-collimation function and data acquisition equipment 4. Marble platform 1 be used for bearing all the other all equipment, isolated external vibrations, marble platform 1 is arranged in to multi-view star sensor 2, each view field visual axis is parallel with the platform plane, theodolite 3 places before each view field of multi-view star sensor, its telescope optical axis and each view field visual axis coincidence, base position fixed back adjustment stabilizer blade flattening theodolite, data acquisition equipment 4 link to each other with multi-view star sensor 2 and theodolite 3 for gather multi-view star sensor 2's image and theodolite 3's measured data.
The invention uses the self luminous cross-hair of the theodolite as the reference source, combines the reference source and the measuring device into a whole, reduces the complexity of the construction of the calibration system, adopts a cross-hair imaging method in the calibration process, does not need to use a reference cube mirror to establish a middle conversion coordinate system, simplifies the coordinate system conversion process, and simultaneously reduces the large-scale leveling work of the equipment placed on the marble platform in the practical application because the marble platform has high flatness.
The invention discloses a calibration method of a multi-view-field star sensor based on theodolite cross hair imaging, which uses the system of the invention and comprises the following steps:
step one, modeling is carried out on the process that the cross hairs of the theodolite are imaged on the star sensor image surface through the telescope objective lens and the star sensor lens.
In the system, each theodolite and one of the fields of view of the multi-field star sensor form a sub-calibration system, and the sub-system is used for calibrating the rotation relationship between the theodolite and the corresponding field of view. When the telescope of the theodolite is stabilized at a certain position, the cross-hair lighted by the auto-collimation function is imaged on the star sensor image surface through the telescope objective lens and the star sensor lens, and the direction of the principal ray of the cross-hair cross point is coincided with the direction of the optical axis of the telescope. The horizontal angle, the vertical angle and the installation angle of the theodolite are combined, and the direction vector of the optical axis of the telescope can be calculated. And combining the image coordinates of the cross-hair cross point and the internal parameters of the star sensor to calculate the direction vector of the main ray of the cross point. The two vectors are parallel to each other and opposite in direction, and the two vectors are used as a central pivot to establish a cross hair imaging model.
Star sensor single view field coordinate system OS-XSYSZSAt the optical center of the lens corresponding to the field of view, XSAxis and YSThe axes being parallel to the grid direction of the image sensor, X, respectivelySThe axis pointing to the right of the image, YSThe axis pointing above the image, ZSThe axis points to the front of the lens along the direction of the optical axis of the lens. During the calibration process, all vectors are converted into the coordinate system for calculation.
As shown in FIG. 2, the vector v of the theodolite telescope in the direction of the optical axisTDirectly under the zero coordinate system of the theodolite. Zero coordinate system OT-XTYTZTFixedly connected to the base, Z thereofTOTXTThe plane being parallel to the local horizontal plane, YTThe axis extending in the local vertical direction downwards, ZTThe axis direction coincides with the telescope pointing direction when the horizontal angle is 0 degrees and the vertical angle is 90 degrees, XTThe axial direction is determined by the right hand rule. When the horizontal angle of the theodolite is Hz and the vertical angle is V, the zero coordinate system OT-XTYTZTDirection vector v of lower telescope optical axisTThe expression of (a) is:
suppose that three-axis Euler angles (YXZ rotation order) from the theodolite zero coordinate system to the star sensor single-view field measurement coordinate system are respectively theta1,θ2,θ3Then the direction cosine matrix corresponding to the conversion relation
Watch (A)The expression is as follows:
according to the direction cosine matrix, the star sensor view field coordinate system OS-XSYSZSDirection vector v of lower telescope optical axisSThe expression of (a) is:
since the length of the direction vector is fixed to 1 and has substantially only 2 degrees of freedom, v is set for convenience of useSScaling to the plane where z is 1, the expression is:
as shown in fig. 3, the relationship between the image coordinates of the cross-hair intersection and the direction vector of the principal ray of the intersection can be established according to the star sensor imaging model. The internal parameters of the star sensor used by the model comprise a focal length f and a principal point (x)c，yc) Coefficient of radial distortion p1，p2And tangential distortion coefficient q1，q2The calibration of the parameters is completed before the calibration of the structural parameters of the multi-view-field star sensor.
Cross point image coordinates (x)p，yp) The image coordinate system is established under an image coordinate system, the origin of the image coordinate system is positioned in the center of a pixel at the lower left corner of an image, the X axis and the Y axis are parallel to the grid direction of the image, the X axis points to the right of the image, and the Y axis points to the upper part of the image. Coordinate (x) of image plane coordinate system for correcting image distortion with origin located under image coordinate systemc，yc) Here, the X-axis and Y-axis directions are the same as the image coordinate system, and a length unit is used. Thus, the image plane coordinate (x) of the intersectionp1，yp1) Comprises the following steps:
xp1＝DX·(xp-xc)
yp1＝DY·(yp-yc)
where DX and DY are pixel lengths in X-axis and Y-axis directions, respectively.
From the image distortion model, the coordinates (x) on the actual image planep1，yp1) The corresponding distortion amounts (Δ x, Δ y) are:
wherein p is1，p2Is the radial distortion coefficient, q1，q2Is the tangential distortion coefficient. Then the undistorted image plane coordinate (x) of the intersection pointp2，yp2) Comprises the following steps:
xp2＝xp1-Δx
yp2＝yp2-Δy
constructing a direction vector w of a main ray of a cross point according to a small hole imaging model of a star sensorSComprises the following steps:
for convenient use, w isSAlso scaled to the z-1 plane, the expression:
and step two, establishing a calibration track, and collecting cross images shot by the star sensor at different calibration positions and theodolite measurement angles.
Lighting the autocollimation cross wire on the graticule of theodolite, focusing the telescope to infinity, and regulating the brightness of cross wire to make the starThe brightest part of the cross-hair image captured by the sensor is just not saturated in gray scale, as shown in fig. 4. Manually adjusting the horizontal angle and the vertical angle of the theodolite to make the image coordinate of the cross point approximately coincide with the coordinate of the principal point, and recording the horizontal angle and the vertical angle (Hz) at the moment0,V0)。
Horizontal and vertical angles (Hz) recorded as described above0,V0) Centered, calibration trajectories spaced 1 ° apart within a 14 ° field of view were constructed, as shown in fig. 5. And adjusting the angle of the theodolite to a first calibration position, keeping the angle of the theodolite unchanged, controlling the star sensor to shoot 100 images, and simultaneously recording the horizontal angle and the vertical angle of the theodolite for 100 times. And after the images and the angles are recorded, controlling the theodolite to move to the next calibration position, and repeating the recording process until the data recording at the last calibration position is finished.
And step three, extracting the cross image into a cross point image coordinate by using an image processing method, preprocessing the acquired data, and constructing an optimization problem according to a cross imaging model to solve the rotation relation between the theodolite coordinate system and the view field coordinate system.
And (3) eliminating coarse measurement results by using a 3 sigma criterion for 100 groups of horizontal angles and vertical angles collected from each calibration position, and respectively averaging the rest horizontal angles and the rest vertical angles to obtain actual horizontal angles and actual vertical angles of the corresponding calibration positions. For the ith calibration position, note as (Hz)i,Vi)。
And processing images shot by all the star sensors by using a quadric surface fitting method to obtain the sub-pixel image coordinates of the cross-hair cross points. The method firstly traverses the whole image, searches the pixel point with the maximum gray value, records the coordinate (x) of the pixel0,y0) As pixel level coordinates of the intersection. Further, the gradient and Hessian matrix at the intersection pixel level coordinates are calculated using a first order differential operator and a second order differential operator. The first order differential operator and the second order differential operator are constructed according to the following formulas:
convolving the image f (x, y) by using the operator to obtain the first-order partial derivative g of the imagex(x,y),gy(x, y) and second partial derivative Hxx(x,y),Hxy(x,y),Hyy(x, y), the expression is:
further, an image gradient g and a Hessian matrix H at the pixel level coordinate of the cross point can be obtained, and the expression is as follows:
fitting the gray scale change condition in the cross point neighborhood into a quadric according to the gradient g and the Hessian matrix H, and calculating the deviation s between the maximum point of the quadric and the current point, wherein the calculation formula is as follows:
the maximum point of the quadric surface is (x)0+s1,y0+s2) Is denoted as (x)p,yp) This coordinate corresponds to the sub-pixel level coordinate of the cross-hair like intersection.
And (3) eliminating coarse measurement results from the coordinates of the cross points extracted from the 100 images shot at each calibration position by using a 3 sigma criterion, and averaging the coordinates of all the remaining cross points to obtain the accurate cross point coordinates of the corresponding calibration position. For the ith calibration position, note as (x)p,i,yp,i)。
According to the model in the step one, the vectors of the directions of the optical axes of the telescope corresponding to the horizontal angle and the vertical angle
Cross point chief ray direction vector corresponding to cross point coordinate
Are substantially the same vector. To be provided with
And
the difference is a residual function, and a minimum optimization problem is constructed for the rotational Euler angle theta of the theodolite coordinate system and the single-view field coordinate system1,θ2,θ3Solving, and constructing a minimal value optimization problem as follows:
where m is the total number of calibration positions, #i＝(Hzi,Vi)TTheodolite angle, x, for the ith calibration positioni＝(xp,i,yp,i)TFor the coordinates of the intersection image at the ith calibration position, p ═ f, xc,yc,p1,p2,q1,q2)TIs the internal parameter of the star sensor, including the focal length f and the principal point (x)c,yc) Radial distortion factor (p)1,p2) And tangential distortion coefficient (q)1,q2) The residual function expression is:
wherein the content of the first and second substances,
and
the vectors v of the optical axis directions of the telescopes respectively corresponding to the ith calibration positionS,iAnd ith calibration position cross point chief ray direction vector wS,iScaling to z is a vector after 1 plane.
And solving the optimization problem by using a Levenberg-Marquardt nonlinear least square method to obtain a rotation Euler angle between a theodolite coordinate system and a single-view field coordinate system. For convenient use, the Euler angle theta is rotated1,θ2,θ3Conversion into a directional cosine matrix
The expression is as follows:
the calibration process is carried out on each theodolite-single-view sub-calibration system, and the obtained direction cosine matrix is recorded as the field number k
And step four, establishing a rotation relation between theodolite coordinate systems through theodolite mutual aiming.
And (3) lightening the auto-collimation cross-shaped wires on the theodolite reticle, and focusing the telescopes to an infinite distance to enable the telescopes of the two theodolites to be approximately in opposite positions. Observing from the telescope of one theodolite, manually adjusting the horizontal angle and the vertical angle of the theodolite to align the cross image formed by the other theodolite with the center of a cross reticle on a reticle of the telescope of the theodolite. The same alignment is performed for viewing from the other theodolite telescope and the operation is repeated a number of times.
Keeping the angles of the theodolites unchanged, and measuring the horizontal angle and the vertical angle of the two theodolites by data acquisition equipment for 100 times respectively.
Two groups of data collected during the cross-sighting of the theodolite are subjected to 3 sigma criterion to eliminate coarse measurement results, the rest data are averaged, and the results are respectively recorded as (Hz)A,VA) And (Hz)B,VB). The angle which needs to rotate around the Y axis is converted from the coordinate system of the theodolite A to the coordinate system of the theodolite B
Comprises the following steps:
for convenience of use, directional cosine matrices are used
Expressing the conversion relation, the expression is:
mutually aiming all the theodolites in pairs to obtain the rotation relation between coordinate systems of all the theodolites, and recording the direction cosine matrix from the coordinate system of the theodolite No. k to the coordinate system of the theodolite No. l as
And step five, fusing the stepping results to obtain a multi-view-field star sensor structure parameter calibration result.
And fusing all the obtained rotation relations to obtain a multi-view-field star sensor structure parameter calibration result, wherein a direction cosine matrix from a k view field coordinate system to a l view field coordinate system is as follows:
effects of the invention
In order to show the calibration effect of the method, the method is applied to calibrate the structural parameters of the airborne high-dynamic star sensor system with three view fields. The system takes a central No. 1 view field as a main view field, and the finally obtained structural parameter form is a direction cosine matrix converted from a No. 1 view field coordinate system to other view field coordinate systems.
The calibration process of the theodolite-single-view sub-calibration system is shown below by taking the view field No. 1 as an example. The internal parameters of the star sensor of the view field are as follows:
the star sensor images and theodolite angles at 145 calibration positions are collected by a sub-calibration system consisting of a No. 1 view field and a corresponding theodolite, and the data after pretreatment are shown in the following table:
obtaining a direction cosine matrix from a theodolite coordinate system to a view field coordinate system in the sub-calibration system through optimization solution
Comprises the following steps:
the calibration results of the No. 2 view field coordinate system and the corresponding theodolite coordinate system can be obtained by the same method
Comprises the following steps:
through 1 theodolite and 2 theodolite cross sights, it is:
calculating to obtain a direction cosine matrix from the coordinate system of the theodolite No. 1 to the coordinate system of the theodolite No. 2
Comprises the following steps:
furthermore, the direction cosine matrix from the No. 1 view field coordinate system to the No. 2 view field coordinate system
Comprises the following steps:
the same method can be used for solving the direction cosine matrix from the No. 1 view field coordinate system to the No. 3 view field coordinate system
Comprises the following steps:
and completing the multi-view-field structural parameter calibration of the airborne high-dynamic star sensor system.
Compared with the existing calibration method adopting the theodolite, the method does not need to use an additional reference source coordinate system and a reference cube mirror coordinate system, eliminates errors caused by multiple times of conversion of the coordinate system and errors caused by processing deviation of the reference cube mirror, and can ensure higher calibration precision. Compared with another common multi-view-field star sensor calibration method based on a three-axis turntable, the method adopts the theodolite which can be placed at will to construct the calibration system, removes the requirements of the high-precision three-axis turntable on the overall dimension, the whole weight and the view field distribution structure of the multi-view-field star sensor system when the calibration system is constructed, and can be applied to the multi-view-field star sensor system with large dimension, large weight and any view field distribution structure.
The foregoing shows and describes the general principles, essential features, and advantages of the invention. It will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, which are merely illustrative of the principles of the invention, but that various changes and modifications may be made without departing from the spirit and scope of the invention, which fall within the scope of the invention as claimed. The scope of the invention is defined by the appended claims and equivalents thereof.
Claims (8)
1. A calibration system of a multi-view field star sensor based on theodolite cross hair imaging comprises a marble platform (1), a multi-view field star sensor (2), an electronic theodolite (3) and data acquisition equipment (4); the method is characterized in that: the marble platform (1) is used for bearing all other equipment; the multi-view-field star sensor (2) is arranged on the marble platform (1), and the visual axis of each view field is parallel to the plane of the marble platform; the electronic theodolite (3) is placed in front of each view field of the multi-view-field star sensor, the optical axis of a telescope of the electronic theodolite is superposed with the view axes of the view fields, and the base of the electronic theodolite is fixed and then the support foot leveling theodolite is adjusted; the data acquisition equipment (4) is connected with the multi-view-field star sensor (2) and the theodolite (3) and is used for acquiring images of the multi-view-field star sensor (2) and measurement data of the theodolite (3).
2. A calibration method of a multi-view field star sensor based on theodolite cross hair imaging adopts the calibration system of claim 1, and is characterized in that: the method comprises the following steps:
modeling a process of imaging a cross hair of a theodolite on an image surface of a star sensor through a telescope objective and a star sensor lens;
establishing a calibration track, and collecting cross images shot by the star sensor at different calibration positions and theodolite measurement angles;
extracting a cross image into a cross point image coordinate by using an image processing method, preprocessing the acquired data, and constructing an optimization problem according to a cross imaging model to solve the rotation relation between a theodolite coordinate system and a view field coordinate system;
establishing a rotation relation between theodolite coordinate systems through theodolite mutual aiming;
and step five, fusing the stepping results to obtain a multi-view-field star sensor structure parameter calibration result.
3. The method for calibrating the multi-view-field star sensor based on the theodolite cross hair imaging as claimed in claim 2, wherein: the first step further comprises the following steps: when the theodolite telescope is stable at a certain position, the luminous cross hairs clearly image on an image surface through the telescope objective lens and the star sensor lens; the telescope optical axis direction vector under the view field coordinate system is determined by the horizontal angle and the vertical angle of the theodolite and the installation angle of the theodolite coordinate system and the view field coordinate system; the vector of the principal ray direction of the cross-hair cross point under the field-of-view coordinate system is determined by the image coordinate of the cross point and the star sensor imaging model; the two vectors are parallel to each other and in opposite directions.
4. The method for calibrating the multi-view-field star sensor based on the theodolite cross hair imaging as claimed in claim 2, wherein the method comprises the following steps: the second step further comprises the following steps: firstly, adjusting the brightness of cross hairs and the focusing position of a telescope, and establishing a calibration track according to the field condition of a star sensor; and sequentially rotating the theodolite to each calibration position in the calibration track, collecting the star sensor images for multiple times and measuring the horizontal angle and the vertical angle of the theodolite for multiple times.
5. The method for calibrating the multi-view-field star sensor based on the theodolite cross hair imaging as claimed in claim 2, wherein: the third step further comprises the following steps: averaging the horizontal angle and the vertical angle of the theodolite collected at each calibration position to obtain the accurate theodolite angle of the calibration position; extracting sub-pixel level coordinates of the intersection points of the collected star sensor images by using an image processing method, and averaging the coordinates of the intersection points belonging to the same calibration position to obtain accurate image coordinates of the intersection points of the calibration position; according to a cross-hair imaging model of the theodolite, using a telescope optical axis direction vector vSAnd crossDirection vector w of point chief raySThe difference is a residual function, and a minimum optimization problem is constructed for the rotational Euler angle theta of the theodolite coordinate system and the single-view field coordinate system1，θ2，θ3Solving, and constructing a minimal value optimization problem as follows:
where m is the total number of calibration positions, #i＝(Hzi，Vi)TTheodolite angle, x, for the ith calibration positioni＝(xp，i，yp，i)TFor the coordinates of the intersection image at the ith calibration position, p ═ f, xc，yc，p1，p2，q1，q2)TIs the internal parameter of the star sensor, including the focal length f and the principal point (x)c，yc) Radial distortion factor (p)1，p2) And tangential distortion coefficient (q)1，q2) The residual function is expressed as
Wherein the content of the first and second substances,
and
the vectors v of the optical axis directions of the telescopes respectively corresponding to the ith calibration positionS，iAnd ith calibration position cross point chief ray direction vector wS，iScaling to z is a vector after 1 plane. Solving the obtained rotation Euler angle theta1，θ2，θ3Converted into a direction cosine matrix form and recorded as a view field number k
6. The method for calibrating the multi-view-field star sensor based on the theodolite cross hair imaging as claimed in claim 2, wherein: the fourth step further comprises the following steps: aiming cross hairs of the two theodolites at each other and recording the measurement angle to obtain the rotation relation of the coordinate systems of the two theodolites; mutually aiming all the theodolites in pairs to obtain the rotation relation between coordinate systems of all the theodolites, and recording the direction cosine matrix from the coordinate system of the theodolite No. k to the coordinate system of the theodolite No. l as
7. The method for calibrating the multi-view-field star sensor based on the theodolite cross hair imaging as claimed in claim 2, wherein: the fifth step further comprises the following steps: and fusing all the obtained rotation relations to obtain a multi-view-field star sensor structure parameter calibration result, wherein a direction cosine matrix from a k view field coordinate system to a l view field coordinate system is as follows:
8. the method for calibrating the theodolite cross hair imaging-based multi-field star sensor as claimed in any one of claims 2 to 7 is applied to an astronomical navigation device.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110883235.XA CN113607188B (en) | 2021-08-02 | 2021-08-02 | Theodolite cross-hair imaging-based multi-view-field star sensor calibration system and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110883235.XA CN113607188B (en) | 2021-08-02 | 2021-08-02 | Theodolite cross-hair imaging-based multi-view-field star sensor calibration system and method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113607188A true CN113607188A (en) | 2021-11-05 |
| CN113607188B CN113607188B (en) | 2022-07-05 |
Family
ID=78306535
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110883235.XA Active CN113607188B (en) | 2021-08-02 | 2021-08-02 | Theodolite cross-hair imaging-based multi-view-field star sensor calibration system and method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113607188B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114577235A (en) * | 2022-01-28 | 2022-06-03 | 北京控制工程研究所 | Cross-scale calibration method and system for space extremely-high-precision pointing measuring instrument |
| CN117706959A (en) * | 2023-12-12 | 2024-03-15 | 北京控制工程研究所 | Non-cooperative intersection sensor view field adjusting method based on semi-physical test system |
Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102032918A (en) * | 2010-10-20 | 2011-04-27 | 郑州辰维科技股份有限公司 | Method for calibrating direction of three-probe start sensor |
| EP2397816A1 (en) * | 2010-06-18 | 2011-12-21 | Leica Geosystems AG | Method for verifying a surveying instrument's external orientation |
| CN103344256A (en) * | 2013-06-19 | 2013-10-09 | 哈尔滨工业大学 | Laboratory testing method for multi-field-of-view star sensor |
| CN104406607A (en) * | 2014-11-21 | 2015-03-11 | 北京航空航天大学 | Multi-visual field composite optical sensor calibration device and method |
| CN105318891A (en) * | 2014-07-25 | 2016-02-10 | 北京航天计量测试技术研究所 | Star sensor reference cube-prism installation error calibration apparatus |
| CN105737848A (en) * | 2014-12-10 | 2016-07-06 | 上海新跃仪表厂 | System-level star sensor star observation system and star observation method thereof |
| CN109459059A (en) * | 2018-11-21 | 2019-03-12 | 北京航天计量测试技术研究所 | A kind of star sensor outfield conversion benchmark measurement system and method |
| CN110006462A (en) * | 2019-05-23 | 2019-07-12 | 长春工业大学 | Star sensor on-orbit calibration method based on singular value decomposition |
| CN110345970A (en) * | 2019-08-06 | 2019-10-18 | 西安中科微星光电科技有限公司 | A kind of optical navigation sensor scaling method and its equipment |
| CN111024127A (en) * | 2019-12-27 | 2020-04-17 | 苏州大学 | Method and system for detecting inter-satellite angular position error of high-resolution dynamic satellite simulator |
| CN111044077A (en) * | 2019-12-25 | 2020-04-21 | 中国科学院国家空间科学中心 | Calibration method between star sensor measurement coordinate system and reference cube mirror coordinate system |
| CN112432705A (en) * | 2020-10-28 | 2021-03-02 | 同济大学 | Multispectral imaging system and method based on dynamic visual axis adjustment principle |
| CN112596257A (en) * | 2020-12-30 | 2021-04-02 | 中国科学院长春光学精密机械与物理研究所 | Optical axis calibration method of off-axis reflective optical lens |
| CN112802115A (en) * | 2020-12-26 | 2021-05-14 | 长光卫星技术有限公司 | Geometric calibration method and device for multi-focal-plane spliced large-view-field off-axis camera |
-
2021
- 2021-08-02 CN CN202110883235.XA patent/CN113607188B/en active Active
Patent Citations (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2397816A1 (en) * | 2010-06-18 | 2011-12-21 | Leica Geosystems AG | Method for verifying a surveying instrument's external orientation |
| CN102032918A (en) * | 2010-10-20 | 2011-04-27 | 郑州辰维科技股份有限公司 | Method for calibrating direction of three-probe start sensor |
| CN103344256A (en) * | 2013-06-19 | 2013-10-09 | 哈尔滨工业大学 | Laboratory testing method for multi-field-of-view star sensor |
| CN105318891A (en) * | 2014-07-25 | 2016-02-10 | 北京航天计量测试技术研究所 | Star sensor reference cube-prism installation error calibration apparatus |
| CN104406607A (en) * | 2014-11-21 | 2015-03-11 | 北京航空航天大学 | Multi-visual field composite optical sensor calibration device and method |
| CN105737848A (en) * | 2014-12-10 | 2016-07-06 | 上海新跃仪表厂 | System-level star sensor star observation system and star observation method thereof |
| CN109459059A (en) * | 2018-11-21 | 2019-03-12 | 北京航天计量测试技术研究所 | A kind of star sensor outfield conversion benchmark measurement system and method |
| CN110006462A (en) * | 2019-05-23 | 2019-07-12 | 长春工业大学 | Star sensor on-orbit calibration method based on singular value decomposition |
| CN110345970A (en) * | 2019-08-06 | 2019-10-18 | 西安中科微星光电科技有限公司 | A kind of optical navigation sensor scaling method and its equipment |
| CN111044077A (en) * | 2019-12-25 | 2020-04-21 | 中国科学院国家空间科学中心 | Calibration method between star sensor measurement coordinate system and reference cube mirror coordinate system |
| CN111024127A (en) * | 2019-12-27 | 2020-04-17 | 苏州大学 | Method and system for detecting inter-satellite angular position error of high-resolution dynamic satellite simulator |
| CN112432705A (en) * | 2020-10-28 | 2021-03-02 | 同济大学 | Multispectral imaging system and method based on dynamic visual axis adjustment principle |
| CN112802115A (en) * | 2020-12-26 | 2021-05-14 | 长光卫星技术有限公司 | Geometric calibration method and device for multi-focal-plane spliced large-view-field off-axis camera |
| CN112596257A (en) * | 2020-12-30 | 2021-04-02 | 中国科学院长春光学精密机械与物理研究所 | Optical axis calibration method of off-axis reflective optical lens |
Non-Patent Citations (10)
| Title |
|---|
| WEI, XG 等: ""Star sensor calibration based on integrated modelling with intrinsic and extrinsic parameters"", 《MEASUREMENT: JOURNAL OF THE INTERNATIONAL MEASUREMENT CONFEDERATION》 * |
| WEI, XG 等: ""Star sensor calibration based on integrated modelling with intrinsic and extrinsic parameters"", 《MEASUREMENT: JOURNAL OF THE INTERNATIONAL MEASUREMENT CONFEDERATION》, vol. 55, 30 December 2014 (2014-12-30), pages 117 - 124 * |
| ZHANG, CF 等: ""Optimized star sensors laboratory calibration"", 《APPLIED OPTICS》 * |
| ZHANG, CF 等: ""Optimized star sensors laboratory calibration"", 《APPLIED OPTICS》, vol. 57, no. 5, 30 December 2018 (2018-12-30), pages 1067 - 1074 * |
| 孙利 等: ""多视场星敏感器结构参数标定方法"", 《北京航空航天大学学报》 * |
| 孙利 等: ""多视场星敏感器结构参数标定方法"", 《北京航空航天大学学报》, vol. 41, no. 8, 30 December 2015 (2015-12-30), pages 1532 - 1536 * |
| 熊琨等: "基于三轴转台的多视场星敏感器标定方法", 《红外与激光工程》 * |
| 熊琨等: "基于三轴转台的多视场星敏感器标定方法", 《红外与激光工程》, no. 04, 13 December 2018 (2018-12-13), pages 242 - 247 * |
| 贺俊吉等: "结构光三维视觉检测中光条图像处理方法研究", 《北京航空航天大学学报》 * |
| 贺俊吉等: "结构光三维视觉检测中光条图像处理方法研究", 《北京航空航天大学学报》, no. 07, 30 August 2003 (2003-08-30), pages 31 - 35 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114577235A (en) * | 2022-01-28 | 2022-06-03 | 北京控制工程研究所 | Cross-scale calibration method and system for space extremely-high-precision pointing measuring instrument |
| CN114577235B (en) * | 2022-01-28 | 2024-03-15 | 北京控制工程研究所 | Cross-scale calibration method and system for spatial extremely high-precision pointing measuring instrument |
| CN117706959A (en) * | 2023-12-12 | 2024-03-15 | 北京控制工程研究所 | Non-cooperative intersection sensor view field adjusting method based on semi-physical test system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113607188B (en) | 2022-07-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Clarke et al. | The development of camera calibration methods and models | |
| CN104154928B (en) | Installation error calibrating method applicable to built-in star sensor of inertial platform | |
| US6731329B1 (en) | Method and an arrangement for determining the spatial coordinates of at least one object point | |
| CN107782293B (en) | Spacecraft equipment posture information measurement method based on six degree of freedom laser tracking target | |
| CN108759798B (en) | Method for realizing precision measurement of high-precision spacecraft | |
| CN113607188B (en) | Theodolite cross-hair imaging-based multi-view-field star sensor calibration system and method | |
| CN106468544B (en) | Satellite high-precision angle-measuring method based on photoelectric auto-collimator | |
| CN109682399B (en) | Precision verification method for position and pose measurement result of total station based on three-axis turntable | |
| CN110646016B (en) | Distributed POS calibration method and device based on theodolite and vision-assisted flexible base line | |
| CN111665023B (en) | Telescope distortion measuring device and method | |
| CN108154535B (en) | Camera calibration method based on collimator | |
| CN106643735A (en) | Indoor positioning method and device and mobile terminal | |
| CN114001756B (en) | Small-field-of-view star sensor outfield ground star finding method | |
| CN111915685A (en) | Zoom camera calibration method | |
| CN110068313B (en) | Digital zenith instrument orientation method based on projection transformation | |
| Yuan et al. | A precise calibration method for line scan cameras | |
| CN105424060B (en) | A kind of measurement method of aircraft star sensor and strapdown inertial measurement unit installation error | |
| CN108645392B (en) | Camera installation posture calibration method and device | |
| Shi et al. | Automatic astronomical survey method based on video measurement robot | |
| CN109489642B (en) | Dynamic measurement method for relative attitude of two cube mirrors under any spatial attitude | |
| Meyer et al. | Coordinates of lunar features | |
| CN115290008A (en) | Angle calibration algorithm of image measurement collimator | |
| Sekulic et al. | Daniel K. Inouye solar telescope optical alignment plan | |
| CN110866951B (en) | Method for correcting optical axis inclination of monocular camera | |
| CN109141385B (en) | Positioning method of total station instrument without leveling |
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 |