SINGLE STAR-BASED ORIENTATION METHOD USING DUAL-AXIS LEVEL SENSOR

Abstract

Disclosed is a single star-based orientation method using a dual-axis level sensor, which includes a calibration process and an actual calculation process.

Claims

1. A single star-based orientation method using a dual-axis level sensor, comprising: (S1) denoting any two orthogonal side surfaces of a hexahedron on a star sensor as a first side surface and a second side surface, respectively; denoting a component of a vector ν.sub.right on an X.sub.n-axis of a quasi-horizontal reference system as t.sub.4; and denoting a component of the vector ν.sub.right on a Z.sub.n-axis of the quasi-horizontal reference system as t.sub.6; wherein the vector ν.sub.right is obtained by rotating a normal vector of the second side surface around the Z.sub.n-axis of the quasi-horizontal reference system by a preset azimuth angle θ.sub.dir.sup.right; (S2) calculating the t.sub.4 and t.sub.6 through the following equations: $\{\begin{array}{c}{t}_4=-\left(\sin \alpha \sin \left(\Delta {\theta }_x^{\mathrm{pitch}}-{\theta }_x^{\mathrm{pitch}}\right)+\cos \alpha \sin \left(\Delta {\theta }_y^{\mathrm{pitch}}-{\theta }_y^{\mathrm{pitch}}\right)\right)\\ t_6=\left(\cos \alpha \sin \left(\Delta {\theta }_x^{\mathrm{pitch}}-{\theta }_x^{\mathrm{pitch}}\right)-\sin \alpha \sin \left(\Delta {\theta }_y^{\mathrm{pitch}}-{\theta }_y^{\mathrm{pitch}}\right)\right)\end{array};$ wherein α is an angle between an X.sub.s axis of a dual-axis level sensor reference system and an X-axis of a hexahedron reference system on the star sensor; θ.sub.x.sup.pitch is a measurement reading of the dual-axis level sensor on the X.sub.s-axis; θ.sub.y.sup.pitch is a measurement reading of the dual-axis level sensor on a Y.sub.s-axis of the dual-axis level sensor reference system; Δθ.sub.x.sup.pitch is a bias of the dual-axis level sensor on the X.sub.s-axis; and Δθ.sub.y.sup.pitch is a bias of the dual-axis level sensor on the Y.sub.s-axis; (S3) denoting an observation vector of a single astronomical object in the hexahedron reference system as ν.sub.PRI; and denoting a reference vector of the single astronomical object in an inertial reference system as ν.sub.GND, wherein the reference vector of the single astronomical object is obtained from a star catalog, and the ν.sub.GNDand ν.sub.PRI satisfy the following equation:
R.sub.z(θ.sub.z).Math.R.sub.x(θ.sub.x).Math.R.sub.y(θ.sub.y).Math.ν.sub.GND=ν.sub.PRI; wherein R.sub.y(θ.sub.y) denotes that the reference vector ν.sub.GND rotates around a Y-axis of the hexahedron reference system by an azimuth angle θ.sub.y; R.sub.x(θ.sub.x) denotes that the reference vector ν.sub.GND rotates around the X-axis of the hexahedron reference system by a pitch angle θ.sub.x; and R.sub.z(θ.sub.z) denotes that the reference vector ν.sub.GND rotates around a Z-axis of the hexahedron reference system by a roll angle θ.sub.z; (S4) calculating the pitch angle θ.sub.x and the roll angle θ.sub.z according to the t.sub.4 and t.sub.6 through the following equations: $\{\begin{array}{c}{\theta }_x=-\arcsin t_6\\ {\theta }_z=-\arcsin \left(t_4/\cos {\theta }_x\right)\end{array};$ (S5) letting ν.sub.GND0=R.sub.x(−θ.sub.z).Math.R.sub.z(−θ.sub.z).Math.ν.sub.PRI and obtaining ν.sub.GND=R.sub.y(−θ.sub.y).Math.ν.sub.GND0; and calculating the azimuth angle θ.sub.y through the following equation to complete a single star-based orientation: ${\theta }_y=\arctan 2\left(\frac{v1.Math.v6-v3.Math.v4}{v{4}^2+v{6}^2},\frac{v1.Math.v4+v3.Math.v6}{v{4}^2+v{6}^2}\right);$ wherein [ν1 ν2 ν3].sup.T=ν.sub.GND; ν.sub.1-ν3 are components of the ν.sub.GND on three axes of the inertial reference system, respectively; [ν4 ν5 ν6].sup.T=ν.sub.GND0; and ν4-ν6 are components of the ν.sub.GND0 on the X.sub.n-axis, the Y.sub.n-axis and the Z.sub.n-axis, respectively.

2. The single star-based orientation method of claim 1, wherein in the step (S1), the preset azimuth angle θ.sub.dir.sup.right is set through steps of: defining the X.sub.n-axis of the quasi-horizontal reference system as a projection of the X-axis of the hexahedron reference system on a horizontal plane; defining a Y.sub.n-axis of the quasi-horizontal reference system such that the Y.sub.n-axis is orthogonal to the X.sub.n-axis on the horizontal plane and a Z.sub.n-axis of the quasi-horizontal reference system satisfies a right-hand rule, and has a general upward pointing direction; measuring, by a first theodolite and a second theodolite, a normal pitch angle θ.sub.pitch.sup.front of the first side surface and a normal pitch angle θ.sub.pitch.sup.right of the second side surface, respectively; wherein in the quasi-horizontal reference system, a normal vector ν.sub.front of the first side surface is expressed as: $v_{\mathrm{front}}=R_y\left({\theta }_{\mathrm{pitch}}^{front}\right).Math.\left[\begin{array}{c}0\\ 0\\ 1\end{array}\right]=\left[\begin{array}{c}vf1\\ 0\\ v\mathrm{f3}\end{array}\right];$ wherein R.sub.y(θ.sub.pitch.sup.front) denotes that a vector [0 0 1] rotates around a Y.sub.n-axis of the quasi-horizontal reference system by the normal pitch angle θ.sub.pitch.sup.front; and νƒ1, 0 and νƒ3 are components of the normal vector ν.sub.front on the X.sub.n-axis, Y.sub.n-axis and Z.sub.n-axis, respectively; a normal vector ν.sub.right of the second side surface without rotating by the preset azimuth angle θ.sub.dir.sup.right is expressed as: ${\stackrel{\_}{v}}_{\mathrm{right}}=R_y\left({\theta }_{\mathrm{pitch}}^{\mathrm{right}}\right).Math.\left[\begin{array}{c}0\\ 0\\ 1\end{array}\right]=\left[\begin{array}{c}vr1\\ 0\\ vr3\end{array}\right];$ wherein R.sub.y(θ.sub.pitch.sup.right) denotes that the vector [0 0 1] rotates around the Y.sub.n-axis of the quasi-horizontal reference system by the normal pitch angle θ.sub.pitch.sup.right; and νr1, 0 and νr3 are components of the normal vector ν.sub.right on the X.sub.n-axis, Y.sub.n-axis and Z.sub.n-axis, respectively; supposing that the vector ν.sub.right is obtained by rotating the vector ν.sub.right by the preset azimuth angle θ.sub.dir.sup.right and expressing the vector ν.sub.right as follows: $v_{\mathrm{right}}=R_z\left({\theta }_{\mathrm{dir}}^{\mathrm{right}}\right).Math.\left[\begin{array}{c}vr1\\ 0\\ vr3\end{array}\right]=\left[\begin{array}{c}\mathrm{vr}1.Math.\cos {\theta }_{\mathrm{dir}}^{\mathrm{right}}\\ -\mathrm{vr}1.Math.\sin {\theta }_{\mathrm{dir}}^{\mathrm{right}}\\ vr3\end{array}\right];$ wherein R.sub.z(θ.sub.dir.sup.right) denotes that a vector [νr1 0 νr3] rotates around the Z.sub.n-axis by the preset azimuth angle θ.sub.dir.sup.right; and according to an orthogonal relationship between the normal vector ν.sub.front and the normal vector ν.sub.right, obtaining the following equation: ${\theta }_{\mathrm{dir}}^{\mathrm{right}}=a\cos \left(\frac{-\mathrm{vf}3.Math.\mathrm{vr}3}{\mathrm{vf}1.Math.\mathrm{vr}1}\right).$

3. The single star-based orientation method of claim 2, wherein an axis of a measuring lens barrel of the first theodolite coincides with a normal of the first side surface of the hexahedron, and an axis of a measuring lens barrel of the second theodolite coincides with a normal of the second side surface of the hexahedron.

4. The single star-based orientation method for orientation of claim 1, wherein in the step (S2), the angle α, the bias Δθ.sub.x.sup.pitch and the bias Δθ.sub.y.sup.pitch are obtained through steps of: defining the X.sub.n-axis of the quasi-horizontal reference system as a projection of the X-axis of the hexahedron reference system on a horizontal plane; defining a Y.sub.n-axis of the quasi-horizontal reference system such that the Y.sub.n-axis is orthogonal to the X.sub.n-axis on the horizontal plane and the Z.sub.n-axis of the quasi-horizontal reference system satisfies a vertically-upward right-hand rule; allowing the dual-axis level sensor reference system to be coplanar with the hexahedron reference system; and expressing a direction vector of the X.sub.s-axis in the hexahedron reference system as [cos α 0−sin α] and a direction vector of the Y.sub.s-axis in the hexahedron reference system as [−sin α 0−cos α]; denoting a normal vector of the first side surface as ν.sub.front; and obtaining a normal vector ν.sub.up, of a third side surface orthogonal to both the first side surface and the second side surface by cross product, expressed as:
ν.sub.up=ν.sub.right×ν.sub.front; obtaining a transformation matrix R.sub.GND2PRI from the quasi-horizontal reference system to the hexahedron reference system, expressed as:
R.sub.GND2PRI=[ν.sub.frontν.sub.rightν.sub.up]; obtaining a transpose matrix of the transformation matrix R.sub.GND2PRI, expressed as: $R_{PRI2GND}={R^T}_{\mathrm{GND}2\mathrm{PRI}}=\left[\begin{array}{ccc}{t}_1& t_2& t_3\\ t_4& t_5& t_6\\ t_7& t_8& t_9\end{array}\right];$ wherein T represents a transposition operation; and t.sub.1-t.sub.9 are elements of the transformation matrix R.sub.GND2PRI; acquiring a direction vector ν.sub.x, of the X.sub.s-axis in the quasi-horizontal reference system, expressed as: $v_{X_s}=\left[\begin{array}{c}{t}_1\cos \alpha -t_3\sin \alpha \\ t_4\cos \alpha -t_6\sin \alpha \\ t_7\cos \alpha -t_9\sin \alpha \end{array}\right];$ acquiring a direction vector ν.sub.y, of the Y.sub.s-axis in the quasi-horizontal reference system, expressed as: $v_{Y_s}=\left[\begin{array}{c}-t_1\sin \alpha -t_3\cos \alpha \\ -t_4\sin \alpha -t_6\cos \alpha \\ -t_7\sin \alpha -t_9\cos \alpha \end{array}\right];$ building relations between t.sub.4 and t.sub.6 according to the direction vector ν.sub.x, and the direction vector ν.sub.y, expressed as: $\{\begin{array}{c}\sin \left(\Delta {\theta }_x^{\mathrm{pitch}}-{\theta }_x^{\mathrm{pitch}}\right)=t_4\cos \alpha -t_6\sin \alpha \\ \sin \left(\Delta {\theta }_y^{\mathrm{pitch}}-{\theta }_y^{\mathrm{pitch}}\right)=-t_4\sin \alpha -t_6\cos \alpha \end{array};$ according to the relations between t.sub.4 and t.sub.6, obtaining a relation between the measurement reading θ.sub.x.sup.pitch and the bias Δθ.sub.x.sup.pitch and a relation between the measurement reading θ.sub.y.sup.pitch and the bias Δθ.sub.y.sup.pitch, expressed as: $\{\begin{array}{c}{\theta }_x^{\mathrm{pitch}}=\Delta {\theta }_x^{\mathrm{pitch}}-\arcsin \left(t_4\cos \alpha -t_6\sin \alpha \right)\\ {\theta }_y^{\mathrm{pitch}}=\Delta {\theta }_y^{\mathrm{pitch}}-\arcsin \left(-t_4\sin \alpha -t_6\cos \alpha \right)\end{array};$ changing a two-axis level degree of the dual-axis level sensor to obtain (θ.sub.x.sup.pitch, θ.sub.y.sup.pitch) and (t.sub.4, t.sub.6) under different two-axis level degrees; and plugging the (θ.sub.x.sup.pitch, θ.sub.y.sup.pitch) and (t.sub.4, t.sub.6) under different two-axis level degrees into the relation between the measurement reading θ.sub.x.sup.pitch and the bias Δθ.sub.x.sup.pitch and the relation between the measurement reading θ.sub.y.sup.pitch and the bias Δθ.sub.y.sup.pitch to obtain the angle α, the bias Δθ.sub.x.sup.pitch and the bias Δθ.sub.y.sup.pitch.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0041] FIG. 1 is a flow chart of a single star-based orientation method using a dual-axis level sensor;

[0042] FIG. 2 schematically depicts a relationship between a dual-axis level sensor reference system, a hexahedron reference system and a quasi-horizontal reference system; and

[0043] FIG. 3 schematically depicts a calibration of the dual-axis level sensor reference system and the quasi-horizontal reference system.

[0044] In the drawings: 1, servo star sensor; 2, hexahedron reference; 3, first theodolite; 4, second theodolite; and 5, data recording computer.

DETAILED DESCRIPTION OF EMBODIMENTS

[0045] The disclosure will be clearly and completely described below with reference to the accompanying drawings and embodiments.

[0046] As shown in the flow chart in FIG. 1, a single star-based orientation method using a dual-axis level sensor includes a first calibrating process, a second calibrating process and an actual calculating process. A key of the single star-based orientation is to project an observed single astronomical object onto a “virtual” horizontal plane (quasi-horizontal plane) via a dual-axis level sensor to obtain an azimuth angle of the projection, which is then compared with that in an inertial reference system. An azimuth angle difference is a current measured azimuth angle.

[0047] Specifically, the actual calculating process of the single star-based orientation method is provided, which is described as follows.

[0048] (S1) Any two orthogonal side surfaces of a hexahedron on a star sensor are denoted as a first side surface and a second side surface, respectively. A component of a vector ν.sub.right on an X.sub.n-axis of a quasi-horizontal reference system is denoted as t.sub.4. A component of the vector ν.sub.right on a Z.sub.n-axis of the quasi-horizontal reference system is denoted as t.sub.6. The vector ν.sup.right is obtained by rotating a normal vector of the second side surface ν.sub.right around the Z.sub.n-axis of the quasi-horizontal reference system by a preset azimuth angle θ.sub.dir.sup.right.

[0049] (S2) The t.sub.4 and t.sub.6 are calculated through the following equations:

[00013]$\{\begin{array}{c}{t}_4=-\left(\sin \alpha \sin \left(\Delta {\theta }_x^{\mathrm{pitch}}-{\theta }_x^{\mathrm{pitch}}\right)+\cos \alpha \sin \left(\Delta {\theta }_y^{\mathrm{pitch}}-{\theta }_y^{\mathrm{pitch}}\right)\right)\\ t_6=\left(\cos \mathrm{\alpha sin}\left(\Delta {\theta }_x^{\mathrm{pitch}}-{\theta }_x^{\mathrm{pitch}}\right)-\sin \alpha \sin \left(\Delta {\theta }_y^{\mathrm{pitch}}-{\theta }_y^{\mathrm{pitch}}\right)\right)\end{array}$

[0050] wherein α is an angle between an X.sub.s-axis of a dual-axis level sensor reference system and an X-axis of a hexahedron reference system on the star sensor. θ.sub.x.sup.pitch is a measurement reading of the dual-axis level sensor on the X.sub.s-axis; θ.sub.y.sup.pitch is a measurement reading of the dual-axis level sensor on a Y.sub.s-axis of the dual-axis level sensor reference system. Δθ.sub.x.sup.pitch is a bias of the dual-axis level sensor on the X.sub.s-axis. Δθ.sub.y.sup.pitch is a bias of the dual-axis level sensor on the Y.sub.s-axis.

[0051] (S3) An observed vector of a single astronomical object in the hexahedron reference system is denoted as ν.sub.PRI. A reference vector of the single astronomical object in an inertial reference system is denoted as ν.sub.GND, where the reference vector of the single astronomical object is obtained from a star catalog, and the ν.sub.GND and ν.sub.PRI satisfy the following equation:

R.sub.z(θ.sub.z).Math.R.sub.x(θ.sub.x).Math.R.sub.y(θ.sub.y).Math.ν.sub.GND=ν.sub.PRI

[0052] where R.sub.y(θ.sub.y) denotes that the reference vector ν.sub.GND rotates around a Y-axis of the hexahedron reference system by an azimuth angle θ.sub.y. R.sub.y(θ.sub.x) denotes that the reference vector ν.sub.GND rotates around the X-axis of the hexahedron reference system by a pitch angle θ.sub.x. R.sub.z(θ.sub.z) denotes that the reference vector ν.sub.GND rotates around a Z-axis of the hexahedron reference system by a roll angle θ.sub.z.

[0053] (S4) The pitch angle θ.sub.x and the roll angle θ.sub.z are calculated according to the t.sub.4 and t.sub.6 through the following equations:

[00014]$\{\begin{array}{c}{\theta }_x=-\arcsin t_6\\ {\theta }_z=-\arcsin \left(t_4/\cos {\theta }_x\right)\end{array}.$

[0054] It should be noted that a transformation matrix R.sub.GND2PRI from the quasi-horizontal reference system to the hexahedron reference system is described according to a rotation sequence of rotating around the Y-axis by the θ.sub.y, rotating around the X-axis by the θ.sub.x and rotating around the Z-axis by the θ.sub.z (2-1-3 rotation) to obtain relations as:

[00015]$\{\begin{array}{c}{t}_4=\cos {\theta }_x\sin {\theta }_z\\ t_6=\sin {\theta }_x\end{array}.$

[0055] The equations for computing the pitch angle θ.sub.x and the roll angle θ.sub.z are obtained through an arcsin function in accordance with the above relations.

[0056] (S5) ν.sub.GND0=R.sub.x(−θ.sub.x).Math.R.sub.z (−θ.sub.z).Math.ν.sub.PRI is made. According to a rotation sequence of rotating around the Z-axis by the θ.sub.z, rotating around the X-axis by the θ.sub.x and rotating around the Y-axis by the θ.sub.z (3-1-2 rotation), ν.sub.GND R.sub.y (−θ.sub.y).Math.ν.sub.GND0 is obtained. The azimuth angle θ.sub.y is calculated through the following equation to complete a single star-based orientation:

[00016]${\theta }_y=\arctan 2\left(\frac{v1.Math.v6-v3.Math.v4}{v{4}^2+v{6}^2},\frac{v1.Math.v4-v3.Math.v6}{v{4}^2+v{6}^2}\right)$

[0057] where [ν1 ν2 ν3].sup.T=ν.sub.GND.Math.ν.sub.1-ν3 are components of the ν.sub.GND on the X-axis, the Y-axis and the Z-axis, respectively. [ν4 ν5 ν6].sup.T=ν.sub.GND0. ν4-ν6 are components of the ν.sub.GND0 on the X.sub.n-axis, the Y.sub.n-axis and the Z.sub.n-axis, respectively.

[0058] A first calibrating process shown in FIG. 1 is a calibration of a transformational relationship between the quasi-horizontal reference system and the hexahedron reference system.

[0059] The quasi-horizontal reference system is constructed according to an angle information of a normal vector of the first side surface ν.sub.front and the normal vector of the second side surface (considering that the first side surface is a front side surface and the second side surface is a right side surface). Referring to FIG. 2, the X.sub.n-axis of the quasi-horizontal reference system is defined as a projection of the X-axis of the hexahedron reference system on a horizontal plane. A Y.sub.n-axis of the quasi-horizontal reference system is defined such that the Y.sub.n-axis of the quasi-horizontal reference system is constrained by orthogonality and lies in the horizontal plane, namely, the Y.sub.n-axis is orthogonal to the X.sub.n-axis on the horizontal plane and a Z.sub.n-axis of the quasi-horizontal reference system satisfies a right-hand rule, and has a general upward pointing direction.

[0060] A normal pitch angle θ.sub.pitch.sup.front of the first side surface and a normal pitch angle θ.sub.pitch.sup.front of the second side surface are measured by a first theodolite and a second theodolite, respectively. An axis of a measuring lens barrel of the first theodolite coincides with a normal of the first side surface of the hexahedron, and an axis of a measuring lens barrel of the second theodolite coincides with a normal of the second side surface of the hexahedron. In the quasi-horizontal reference system, the normal vector ν.sub.front is expressed as

[00017]$v_{\mathrm{front}}=R_y\left({\theta }_{\mathrm{pitch}}^{\mathrm{front}}\right).Math.\left[\begin{array}{c}0\\ 0\\ 1\end{array}\right]=\left[\begin{array}{c}vf1\\ 0\\ vf3\end{array}\right]$

[0061] where R.sub.y(θ.sub.pitch.sup.front) denotes that a vector [0 0 1] rotates around a Y.sub.n-axis by the normal pitch angle θ.sub.pitch.sup.front. νƒ1, 0 and νƒ3 are components of the second normal vector ν.sub.front on the X.sub.n-axis, the Y.sub.n-axis and the Z.sub.n-axis, respectively. The normal vector ν.sub.right without rotating by the preset azimuth angle θ.sub.dir.sup.right is expressed as:

[00018]${\stackrel{\_}{v}}_{\mathrm{right}}=R_y\left({\theta }_{\mathrm{pitch}}^{\mathrm{right}}\right).Math.\left[\begin{array}{c}0\\ 0\\ 1\end{array}\right]=\left[\begin{array}{c}vr1\\ 0\\ vr3\end{array}\right],$

[0062] where R.sub.y(θ.sub.pitch.sup.right) denotes that the vector [00 1] rotates around the Y.sub.n-axis of the quasi-horizontal reference system by the normal pitch angle θ.sub.pitch.sup.right. νr1, 0 and νr3 are components of the normal vector ν.sub.right on the X.sub.n-axis, Y.sub.n-axis and Z.sub.n-axis, respectively.

[0063] The vector ν.sub.right is obtained by rotating the vector ν.sub.right by the preset azimuth and the vector ν.sub.right is expressed as follows:

[00019]$v_{\mathrm{right}}=R_z\left({\theta }_{\mathrm{dir}}^{\mathrm{right}}\right).Math.\left[\begin{array}{c}vr1\\ 0\\ vr3\end{array}\right]=\left[\begin{array}{c}vr1.Math.\cos {\theta }_{\mathrm{dir}}^{\mathrm{right}}\\ -vr1.Math.\sin {\theta }_{\mathrm{dir}}^{\mathrm{right}}\\ vr3\end{array}\right]$

[0064] where R.sub.z(θ.sub.dir.sup.right) denotes that a vector [νr1 0 νr3] rotates around the Z.sub.n-axis by the preset azimuth angle θ.sub.dir.sup.right.

[0065] According to an orthogonal relationship between the normal vector ν.sub.front and the normal vector ν.sub.right (ν.sub.front ⊥ ν.sub.right), the following equation is obtained:

[00020]${\theta }_{\mathrm{dir}}^{\mathrm{right}}=\mathrm{acos}\left(\frac{-vf3.Math.vr3}{vf1.Math.vr1}\right).$

[0066] A third normal vector ν.sub.up of a third surface is obtained according to a cross product, expressed as

ν.sub.up=ν.sub.right×ν.sub.front.

[0067] Therefore, the transformation matrix R.sub.GND2PRI is expressed as:

R.sub.GND2PRI=[ν.sub.frontν.sub.rightν.sub.up].

[0068] A second calibrating process shown in FIG. 1 is a calibration of a transformational relationship between the dual-axis level sensor reference system and the quasi-horizontal reference system.

[0069] Referring to FIG. 2, an angle between the X.sub.s-axis and the X-axis is assumed as α. Since the dual-axis level sensor reference system is coplanar with the hexahedron reference system, an angle between the Y.sub.s-axis and the X-axis is α+π/2. A direction vector of the X.sub.s-axis in the hexahedron reference system is expressed as [cos α 0−sin α] and a direction vector of the Y.sub.s-axis in the hexahedron reference system is expressed as [−sin α 0−cos α].

[0070] A transpose matrix of the transformation matrix R.sub.GND2PRI is obtained, expressed as follows:

[00021]$R_{\mathrm{PRI}2\mathrm{GND}}=R_{\mathrm{GND}2\mathrm{PRI}}^T=\left[\begin{array}{ccc}{t}_1& t_2& t_3\\ t_4& t_5& t_6\\ t_7& t_8& t_9\end{array}\right]$

[0071] where T represents a transposition operation. t.sub.1-t.sub.9 are elements of the transformation matrix R.sub.GND2PRI.

[0072] A direction vector ν.sub.x, of the X.sub.s-axis in the quasi-horizontal reference system is acquired, expressed as follows:

[00022]$v_{X_s}=\left[\begin{array}{c}{t}_1\cos \alpha -t_3\sin \alpha \\ t_4\cos \alpha -t_6\sin \alpha \\ t_7\cos \alpha -t_9\sin \alpha \end{array}\right].$

[0073] A direction vector ν.sub.y, of the Y.sub.s-axis in the quasi-horizontal reference system is acquired, expressed as follows:

[00023]$v_{Y_s}=\left[\begin{array}{c}-t_1\sin \alpha -t_3\cos \alpha \\ -t_4\sin \alpha -t_6\cos \alpha \\ -t_7\sin \alpha -t_9\cos \alpha \end{array}\right].$

[0074] The measurement reading of the direction vector ν.sub.x, on the dual-axis level sensor is θ.sub.x.sup.pitch. The measurement reading of the direction vector ν.sub.y, on the dual-axis level sensor is θ.sub.y.sup.pitch. The bias of the dual-axis level sensor on the Xs-axis is Δθ.sub.x.sup.pitch. The bias of the dual-axis level sensor on the Ys-axis is Δθ.sub.y.sup.pitch

[0075] Relations between t.sub.4 and t.sub.6 according to the direction vector ν.sub.x, and the direction vector ν.sub.y, is built, expressed as follows:

[00024]$\{\begin{array}{c}\sin \left(\Delta {\theta }_x^{\mathrm{pitch}}-{\theta }_x^{\mathrm{pitch}}\right)=t_4\cos \alpha -t_6\sin \alpha \\ \sin \left(\Delta {\theta }_y^{\mathrm{pitch}}-{\theta }_y^{\mathrm{pitch}}\right)=-t_4\sin \alpha -t_6\cos \alpha \end{array}.$

[0076] According to the relations between t.sub.4 and t.sub.6, a relation between the measurement reading θ.sub.x.sup.pitch and the bias Δθ.sub.x.sup.pitch and a relation between the measurement reading θ.sub.y.sup.pitch and the bias Δθ.sub.y.sup.pitch are obtained, expressed as:

[00025]$\{\begin{array}{c}{\theta }_x^{\mathrm{pitch}}=\Delta {\theta }_x^{\mathrm{pitch}}-\arcsin \left(t_4\cos \alpha -t_6\sin \alpha \right)\\ {\theta }_y^{\mathrm{pitch}}=\Delta {\theta }_y^{\mathrm{pitch}}-\arcsin \left(-t_4\sin \alpha -t_6\cos \alpha \right)\end{array}.$

[0077] The star sensor is placed under a calibration device, as shown in FIG. 3. A two-axis level degree of the dual-axis level sensor is changed to obtain (θ.sub.x.sup.pitch, θ.sub.y.sup.pitch) and (t.sub.4, t.sub.6) under different two-axis level degrees and then multiple mapping relations of (θ.sub.x.sup.pitch, θ.sub.x.sup.pitch).fwdarw.(t.sub.4, t.sub.6) are obtained. The (θ.sub.x.sup.pitch, θ.sub.y.sup.pitch) and (t.sub.4, t.sub.6) under different two-axis level degrees into the relation between the measurement reading θ.sub.x.sup.pitch and the bias Δθ.sub.x.sup.pitch and the relation between the measurement reading θ.sub.y.sup.pitch and the bias Δθ.sub.y.sup.pitch are plugged to obtain the angle α, the bias Δθ.sub.x.sup.pitch and the bias Δθ.sub.y.sup.pitch.

[0078] In the orientation method provided herein, the hexahedron reference system is measured by means of the autocollimation of the theodolite, and the quasi-horizontal reference system is measured by the dual-axis level sensor. Then structural parameters of the quasi-horizontal reference system and the hexahedron reference system can be calibrated according to the pitch angle measured by the theodolite and the reading of the dual-axis level sensor. Accordingly, with the help of the theodolite, the collection of all calibrated data can be completed, simplifying the calibration process.

[0079] In summary, the single star-based orientation method provided herein can be described as follows.

[0080] The first calibrating process is performed as follows. Normal vectors of two orthogonal side surfaces on the hexahedron of the star sensor are measured to obtain the transformational relationship between the quasi-horizontal reference system and the hexahedron reference system. The X.sub.n-axis of the quasi-horizontal reference system is defined as a projection of the X-axis of the hexahedron reference system in a horizontal plane. The Y.sub.n-axis is orthogonal to the X.sub.n-axis on the horizontal plane and the Z.sub.n-axis of the quasi-horizontal reference system satisfies a vertically-upward right-hand rule.

[0081] The second calibrating process is performed as follows. By means of the mapping relation between the measurement reading of the dual-axis level sensor and normal measurements of the theodolites, the transformational relationship between the dual-axis level sensor reference system and the quasi-horizontal reference system is obtained. After the second calibrating process, the angle between an X.sub.s-axis of a dual-axis level sensor reference system and an X-axis of a hexahedron reference system, the bias of the dual-axis level sensor on the X.sub.s-axis, and the bias of the dual-axis level sensor on the Y.sub.s-axis are obtained.

[0082] The actual calculating process is performed as follows. By means of an attitude information of the X.sub.s-axis and the Y.sub.s-axis obtained in the second calibrating process, the vectors expressed in the quasi-horizontal reference system is obtained through the observation vector of the single astronomical object. An azimuth angle difference between the quasi-horizontal reference system and the inertial reference system, that is orientation information, is obtained by combining with the reference vector of the single astronomical object in the inertial reference system from the star catalog.

[0083] This application eliminates an uncertainty of solution when solving for an attitude with a single vector by means of the information from the dual-axis level sensor. Furthermore, a quasi-horizontal reference system is constructed through the 2-1-3 rotation. Vectors of the single astronomical object are projected onto the quasi-horizontal reference system. By comparing the azimuth angle difference between a projected azimuth angle and an azimuth angle in the inertial reference system, the single star-based orientation of is completed.

[0084] Mentioned above are merely some embodiments of this disclosure, which are not intended to limit the disclosure. It should be understood that any changes, replacements and modifications made by those killed in the art without departing from the spirit and scope of this disclosure should fall within the scope of the present disclosure defined by the appended claims.