Initial Alignment for SINS Based on Pseudo-Earth Frame in Polar Regions

An accurate initial alignment must be required for inertial navigation system (INS). The performance of initial alignment directly affects the following navigation accuracy. However, the rapid convergence of meridians and the small horizontal component of rotation of Earth make the traditional alignment methods ineffective in polar regions. In this paper, from the perspective of global inertial navigation, a novel alignment algorithm based on pseudo-Earth frame and backward process is proposed to implement the initial alignment in polar regions. Considering that an accurate coarse alignment of azimuth is difficult to obtain in polar regions, the dynamic error modeling with large azimuth misalignment angle is designed. At the end of alignment phase, the strapdown attitude matrix relative to local geographic frame is obtained without influence of position errors and cumbersome computation. As a result, it would be more convenient to access the following polar navigation system. Then, it is also expected to unify the polar alignment algorithm as much as possible, thereby further unifying the form of external reference information. Finally, semi-physical static simulation and in-motion tests with large azimuth misalignment angle assisted by unscented Kalman filter (UKF) validate the effectiveness of the proposed method.


Introduction
With the development of worldwide transportation, the importance of polar navigation is increasing. The polar navigation technology can ensure the safety and reliability of vehicles when they sail in polar regions. However, due to the special geographical and electromagnetic characteristics [1], common navigation methods such as satellite navigation, radio navigation and geomagnetic navigation do not always work efficiently in polar regions [2,3]. As a result, inertial navigation system (INS), which is a highly autonomous and self-contained navigation system, becomes an optimal choice for polar navigation [4,5].
However, in polar regions, the rapid convergence of meridians makes the traditional north-oriented mechanization ineffective because the characteristic that the y-axis always points toward the north would result in both mechanical and mathematical difficulties in a platform instrument or strapdown attitude matrix [6,7]. As a result, the wander navigation system [8], grid navigation system [9], and transverse navigation system [10] have been designed and developed to solve the problem. Nonetheless, any one of these navigation systems does not have the ability of global navigation. Consequently, a combination of them is usually applied to achieve the global navigation [11]. To benefit from the understandability and familiarity, the local geographic frame is being now widely used as the navigation frame in mid-low latitude regions. As a result, the combination of the north-oriented navigation system and a wander navigation system, or a grid navigation system, or a transverse degree of observability of azimuth [25]. Hence, the external velocity reference is preferred selection for the observation information of filter. For velocity reference information, there are usually two forms: velocity information projected in body frame (V b ) and velocity information projected in the navigation frame (V n ). However, polar alignment assisted by V b would also be difficult to achieve due to the smaller horizontal component of rotation of Earth. As a result, V n should be selected as the observation information of filter to achieve the polar alignment [19,21]. This means that polar alignment with the grid frame or transverse frame needs different forms of velocity reference information due to the fact that the different attitude reference frames are employed in the grid alignment or transverse alignment. As a result, the grid alignment or transverse alignment would be disadvantageous to achieve a unified polar alignment algorithm and to unify the form of external reference information further.
According to above analysis and considering the features of global navigation of SINS, a polar alignment algorithm, which can directly obtain the attitude matrix relative to local geographic frame when the initial alignment is finished, is expected to be proposed. The transformation of attitude matrix is not influenced by position errors. Moreover, a unified polar alignment algorithm is also expected to be explored as much as possible. Then the form of external observation information of filter can also be unified for polar alignment in future. Furthermore, it would also be advantage to implement the acquirement of the observation information in a form from the external navigation device. On the other hand, since the definition of coordinate system just like that in the traditional Earth frame, the transverse frame is more immediate and convenient [19]. Therefore, an improved SINS mechanization based on pseudo-Earth frame [26], which is a generalized transverse Earth frame, is proposed in this paper. Then, a novel alignment algorithm based on pseudo-Earth frame and backward process is proposed to achieve the initial alignment in polar regions.
The purpose of this research is expected to provide a new idea to solve the problem of polar alignment and explores a unified polar alignment algorithm as much as possible. The rest of this paper is organized as follows: the proposed polar alignment scheme based on pseudo-Earth frame and backward process is presented in Section 2. In Section 3, the dynamic error modeling with large azimuth misalignment angle is designed. The nonlinear model of polar alignment is given in Section 4. The semi-physical static simulation and in-motion tests with the proposed method are presented in Section 5, and Section 6 is the conclusion.

The Mechanization of SINS with Pseudo-Earth Frame
The pseudo-Earth frame (p-frame) can be obtained by revolving Earth-Centered Earth-Fixed frame (e-frame). As shown in Figure 1, the origin of the p-frame (o p x p y p z p ) is still at the Earth's mass center; N p denotes the pseudo-North Pole; the o p x p axis overlaps with previous o e z e axis; the o p y p axis is pointing from the center of Earth to the projection point of initial position of vehicle S on the equatorial plane; and the o p z p axis is perpendicular to the o p y p axis and conforms to the right hand rule, while laying in the equatorial plane. The detailed definition of the p-frame can be found in [26]. Supposing the initial longitude, latitude and altitude of vehicle are λ 0 , L 0 and h 0 , respectively. Then, the rotating matrix between the p-frame and e-frame can be written as:  For the definitions of pseudo-Earth coordinate system, they are similar to the ones in traditional Earth coordinate system [26]. ) is also similar to the East-North-Up geographic frame (namely, g -frame), which is shown in Figure  1. The initial pseudo position of vehicle can also be given by: According to above definitions, the pseudo-Earth frame is a generalized transverse Earth frame built with the initial position of vehicle in fact [17]. Therefore, the pseudo navigation mechanization of SINS is also similar to the one with transversal Earth coordinate system. In pseudo navigation mechanization, the p g -frame is selected as navigation frame ( n -frame). Only the angular velocity of p -frame relative to inertial frame ( i -frame) differs from that in the local-level coordinate system. As a result, the pseudo navigation equations can be written by [17,26]: where (  ) denotes the skew symmetric matrix of (  ); c ab  denotes the rotation angular rate of frame b relative to frame a represented in the c frame; sin sin cos ;  denotes the angular rate of Earth rotation;  denotes the outputs of gyro; b f denotes the specific force; and n g denotes the gravity vector in the n-frame, For the definitions of pseudo-Earth coordinate system, they are similar to the ones in traditional Earth coordinate system [26]. The position of vehicle is determined by pseudo longitude λ p , pseudo latitude L p , and pseudo altitude h p . The pseudo geographic frame o g p x g p y g p z g p (namely, E p N p U p ) is also similar to the East-North-Up geographic frame (namely, g-frame), which is shown in Figure 1. The initial pseudo position of vehicle can also be given by: According to above definitions, the pseudo-Earth frame is a generalized transverse Earth frame built with the initial position of vehicle in fact [17]. Therefore, the pseudo navigation mechanization of SINS is also similar to the one with transversal Earth coordinate system. In pseudo navigation mechanization, the g p -frame is selected as navigation frame (n-frame). Only the angular velocity of p-frame relative to inertial frame (i-frame) differs from that in the local-level coordinate system. As a result, the pseudo navigation equations can be written by [17,26]: . .
where ( * ×) denotes the skew symmetric matrix of ( * ); ω c ab denotes the rotation angular rate of frame b relative to frame a represented in the c frame; ω n ip = −Ω sin λ p − Ω sin L p cos λ p Ω cos L p cos λ p T ; Ω denotes the angular rate of Earth rotation; ω n pn ib denotes the outputs of gyro; f b denotes the specific force; and g n denotes the gravity vector in the n-frame, g n ≈ 0 0 −g T . Since the pseudo-Earth frame is a generalized transverse Earth frame and the pseudo latitude at the initial position (namely, L p0 ) is always zero from Equation (2), the pseudo-Earth frame can solve the problem of attitude reference frame of polar alignment. In addition, the polar alignment algorithm based on pseudo-Earth frame is also easy to implement because it is similar to the polar alignment algorithm with the transverse Earth frame. However, if the strapdown attitude matrix relative to g-frame is obtained by direct transformation with both the pseudo attitude matrix and the relationship of them, the attitude matrix would be still influenced by the position errors. The cumbersome derivation and computational complexity must also be confronted. Nevertheless, by computation, it is easily find that the initial transformation matrix between g p -frame and g-frame at initial position is a constant matrix and never change with the initial position of vehicle. As a result, it would be very advantageous to achieve the transformation of attitude matrix between g-frame and g p -frame without the influence of position errors and the cumbersome computation. The initial transformation matrix C g p g (0) can also be given by: Substituting Equations (1) and (2) into Equation (6), the initial transformation matrix C g p g (0) can be rewritten by: From Equation (7), we can obtain that the g p -frame is identical to the South-East-Up geographic frame at initial position. As a result, if the pseudo attitude matrix of vehicle at initial position is obtained when the initial alignment is completed, it would be equal to the attitude matrix relative to local geographic frame is obtained without influence of transformation errors and cumbersome computation. Moreover, the coarse alignment method needs also not to be redesigned, and the traditional coarse alignment algorithm can be applied directly. Therefore, the pseudo navigation mechanization would be simpler to solve the problem of attitude reference frame for polar alignment. Nevertheless, the pseudo attitude matrix of vehicle at initial position is expected to be obtained when the initial alignment based on pseudo navigation mechanization is completed.
Remark. If the initial latitude of vehicle relative to e-frame is 90 • , the initial longitude of vehicle can choose the arbitrary value with [−180 • , 180 • ]. Then, the pseudo-Earth frame is built based on the selected longitude. The initial azimuth angle can still be obtained by the initial pseudo azimuth angle from Equation (7). Nevertheless, we should notice that the azimuth angle relative to g-frame in pole is the angle between the vehicle direction and the selected longitude.

Backward Process for Pseudo Navigation System
With electronic technology development, the capabilities of the navigation computer become large. The backward process of data has been proposed for initial alignment [27][28][29]. By utilizing the saved data of inertial measurement unit (IMU) with the opposite time sequence, the backward process can accelerate the process of alignment. Nevertheless, here the requirement that the pseudo attitude matrix of vehicle at initial position is obtained at the end of alignment phase is expected to be realized by utilizing the saved IMU data. If the pseudo starpdown attitude matrix and the pseudo position in backward process are identical to the forward ones at the same time, the pseudo position with backward process can return to initial position. Then, the initial attitude matrix relative to g-frame can be obtained from Equation (7). On the other hand, the backward process can still be applied to decrease the alignment time.
From Equations (3)-(5), the discrete form of forward pseudo navigation equations can be given as follows: where k is the discrete time index for forward computation, and k = 0, 1, 2, · · · m; A k denotes the computed value of A at time t k , while A k−1 denotes the computed value of A at time t k−1 ; T denotes the sampling period; I denotes the identity matrix. Then, through further transposition, the backward pseudo navigation equations can also be obtained: Moreover, assuming the discrete time index for backward computation is l(l = 0, 1, 2, · · · m), the backward pseudo navigation equations can be written as follows: where the superscript "←" denotes this term is for the backward process.
Defining ξ denotes the ξ th sample with the normal time sequence (ξ = 1, 2, · · · m) and η denotes the η th sample with the opposite time sequence (η = 1, 2, · · · m). According to the forward computed process and backward computed process above, the corresponding relations of time index can be shown in Figure 2.  According to the analysis above, the backward strapdown attitude matrix should equal to the forward one at the same time. As a result, the corresponding terms of Equations (11) and (14) should also be equal. Comparing Equation (11) and Equation (14), we can obtain: Considering the corresponding relations of time index from Figure 2, we further have: Obviously, k (19) and (20) can be rewritten by: On the other hand, it is also expected that the backward pseudo position could be identical to the forward ones at the same time. Thus, comparing the corresponding terms of Equations (13) and (16), we can similarly obtain: Then, we can further have: From Equation (23), we can conclude that the sign of the pseudo velocity in backward process should be opposite to the corresponded one in forward process. Therefore, after taking opposite, Equation (12) should equal to Equation (15), correspondingly. Then we can get: From Equation (24), we have: According to the analysis above, the backward strapdown attitude matrix should equal to the forward one at the same time. As a result, the corresponding terms of Equations (11) and (14) should also be equal. Comparing Equation (11) and Equation (14), we can obtain: Considering the corresponding relations of time index from Figure 2, we further have: Obviously, ξ = k, η = l. Combining Equations (19) and (20) can be rewritten by: On the other hand, it is also expected that the backward pseudo position could be identical to the forward ones at the same time. Thus, comparing the corresponding terms of Equations (13) and (16), we can similarly obtain: Then, we can further have: From Equation (23), we can conclude that the sign of the pseudo velocity in backward process should be opposite to the corresponded one in forward process. Therefore, after taking opposite, Equation (12) should equal to Equation (15), correspondingly. Then we can get: From Equation (24), we have: From Equations (25) and (26), the signs of the outputs of accelerometer and the gravity in backward process should be the same as the forward ones. Moreover, substituting Equation (23) into Equation (27), we can obtain: From Equation (28), we can conclude that the sign of the angular rate of rotation of Earth in backward process needs to be taken negative. In addition, combining with Equations (21), (28), and (29), the outputs of gyro in backward process can be obtained by: From Equation (30), the sign of outputs of gyro should be opposite to the forward one in the backward process.
According to above derivation, the expectation that the pseudo starpdown attitude matrix and the pseudo position in backward process should be the same as the forward ones at the same time can be realized. While the pseudo velocity of backward process is opposite to the forward one at the same time. Moreover, the signs of the outputs of gyro and the angular rate of rotation of Earth in backward process should be opposite to the corresponded ones with the forward process. The others are all the same as the forward ones. As a result, the expectation that the pseudo attitude matrix of vehicle at initial position is obtained at the end of alignment phase can be realized by backward process with pseudo navigation mechanization. Further, the initial transformation matrix C g p g (0) can be utilized to obtain the strapdown attitude matrix relative to g-frame with a simple transformation logic and without influence of position errors.

The Proposed Polar Alignment Scheme with Pseudo-Earth frame
From the analysis above, the alignment algorithm based on pseudo-Earth frame and backward process can be applied to solve the problem of polar alignment. The attitude matrix relative to local geographic frame is obtained when the initial alignment is completed. Moreover, it is also easy to implement without the influence of position errors and the cumbersome computation since the initial transformation matrix is a fixed known constant matrix and never influenced by initial position of vehicle from Equation (7).
Supposing the IMU data are saved in first forward fine alignment and the sampling index of saved data sequence is 1 to m. Moreover, the computed pseudo attitude matrix C g p b (m) 1 is also obtained at the end of first forward process. Subsequently, the backward process of the filter estimation is carried out from m to 1, and the pseudo attitude matrix C g p b (0) 2 is obtained at the end of first backward process. Then, the forward process and backward process based on pseudo-Earth frame is repeatedly conducted until the initial alignment is completed. Comparing with the traditional application of the backward process, one and only difference is that the proposed alignment algorithm should end with the backward process the filter estimation from m to 1. Meanwhile, the computed pseudo attitude matrix of time 0 (namely, C g p b (0) n ) is obtained. As a result, the pseudo attitude matrix at initial position is obtained by the backward process. The attitude matrix relative to g-frame is also obtained from Equation (7). Then, the navigation parameters relative to g-frame can be transformed to the corresponding reference frame of the following polar navigation system by the existing transformation program without any additional workloads and complexities. In addition, since the reciprocating process spends some time, the IMU data also need to be saved. The sampling index of saved data in the delay duration is represented from m to s. According to the transformational initial navigation parameters, the navigation calculation is conducted from the first sample to the sth sample by the selected polar navigation system, such as the wander navigation system, grid navigation system and transverse navigation system and so on. Then, the SINS accesses the navigation phase. The proposed polar alignment scheme based on pseudo-Earth frame and backward process is shown in Figure 3. sample to the sth sample by the selected polar navigation system, such as the wander navigation system, grid navigation system and transverse navigation system and so on. Then, the SINS accesses the navigation phase. The proposed polar alignment scheme based on pseudo-Earth frame and backward process is shown in Figure 3. In addition, the implementation of the proposed alignment algorithm is also illustrated in Figure  4. As shown in Figure 4, the traditional coarse alignment algorithm can be still applied to the proposed polar alignment process. Only an initialization with simple transformation is conducted to access the next fine alignment associated with pseudo-Earth frame. Moreover, since the initial pseudo latitude is always zero, the alignment model with the pseudo-Earth frame would be always equivalent to the traditional alignment model of equatorial regions. As a result, the initial alignment with the pseudo-Earth frame can solve the attitude reference problem of polar alignment.
Furthermore, since the initial transformation matrix (0) p g g C is a fixed known constant matrix, the backward process of SINS is proposed here, thereby obtaining the pseudo attitude matrix of vehicle at initial position at the end of initial alignment. As a result, the initial transformation matrix can be cleverly applied to obtain the strapdown attitude matrix relative to g -frame with a simple transformation and without influence of position errors. Thus, there is less cumbersome derivation and low computational complexity in the process of transformation, thereby decreasing the complexity of program and the burden of computer. Moreover, the backward process can still accelerate the filtering estimation and decrease the alignment time. On the other hand, the strapdown attitude matrix relative to g -frame is obtained at the end of alignment phase. It would also be more convenient to access any one of the following polar navigation systems. Because the existing transformation program of navigation parameters in global navigation system can be utilized directly, it would make it possible to achieve the unification of polar alignment algorithm. Therefore, the proposed polar alignment algorithm has superior performance. In addition, the implementation of the proposed alignment algorithm is also illustrated in Figure 4. As shown in Figure 4, the traditional coarse alignment algorithm can be still applied to the proposed polar alignment process. Only an initialization with simple transformation is conducted to access the next fine alignment associated with pseudo-Earth frame. Moreover, since the initial pseudo latitude is always zero, the alignment model with the pseudo-Earth frame would be always equivalent to the traditional alignment model of equatorial regions. As a result, the initial alignment with the pseudo-Earth frame can solve the attitude reference problem of polar alignment. Furthermore, since the initial transformation matrix C g p g (0) is a fixed known constant matrix, the backward process of SINS is proposed here, thereby obtaining the pseudo attitude matrix of vehicle at initial position at the end of initial alignment. As a result, the initial transformation matrix can be cleverly applied to obtain the strapdown attitude matrix relative to g-frame with a simple transformation and without influence of position errors. Thus, there is less cumbersome derivation and low computational complexity in the process of transformation, thereby decreasing the complexity of program and the burden of computer. Moreover, the backward process can still accelerate the filtering estimation and decrease the alignment time. On the other hand, the strapdown attitude matrix relative to g-frame is obtained at the end of alignment phase. It would also be more convenient to access any one of the following polar navigation systems. Because the existing transformation program of navigation parameters in global navigation system can be utilized directly, it would make it possible to achieve the unification of polar alignment algorithm. Therefore, the proposed polar alignment algorithm has superior performance. accelerate the filtering estimation and decrease the alignment time. On the other hand, the strapdown attitude matrix relative to g -frame is obtained at the end of alignment phase. It would also be more convenient to access any one of the following polar navigation systems. Because the existing transformation program of navigation parameters in global navigation system can be utilized directly, it would make it possible to achieve the unification of polar alignment algorithm. Therefore, the proposed polar alignment algorithm has superior performance.

Pseudo Navigation Error Equations with Large Azimuth Misalignment Angle
Considering that an accurate coarse alignment of azimuth is difficult to obtain in polar regions, the dynamic error modeling with large azimuth misalignment angle is designed in this section. Defining the pseudo misalignment angle between computed navigational frame (n -frame) and n-frame is Here, the pseudo-pitch and pseudo-roll misalignment angles (namely, φ xp and φ yp ) are small misalignment angles. While the pseudo-azimuth angle should be considered as a large misalignment angle. Then, the rotating matrix between n -frame and n-frame can be written as [30]: cφ yp cφ zp − sφ yp sφ xp sφ zp cφ yp sφ zp + sφ yp sφ xp cφ zp −sφ yp cφ xp −cφ xp sφ zp cφ xp cφ zp sφ xp sφ yp cφ zp + cφ yp sφ xp sφ zp sφ yp sφ zp − cφ yp sφ xp cφ zp cφ yp cφ xp    (31) where sφ jp and cφ jp represent the sine function and cosine function, respectively, and j = x, y, z.
Defining the angular velocity of n -frame relative to n-frame is ω n nn , and the relationship between ω n nn and φ p can be given by [30]: From Equation (32), we have: where Only considering that the pseudo-azimuth angle is a large misalignment angle here, Equations (31) and (34) can be rewritten by: Due to the existence of misalignment angles and gyro drifts, the differential equation of actual computed pseudo strapdown attitude matrix can be written as: where where ε b denotes the gyro drifts, δω n in denotes the calculated error of ω n in . Since C n n = C n n (ω n nn ×) into Equation (40), we can obtain: Multiplying the left by C n n and the right by C b n , we further have: Then, Equation (42) can be rewritten by: Substituting Equation (38) and ω b nb = ω b ib − C b n ω n in into Equation (43), we can obtain: From Equation (33), the pseudo attitude error equation with large azimuth misalignment angle can be obtained by: where δω n in = δω n ip + δω n pn (46) T ω 11 = −Ω cos λ p δλ p ω 12 = −Ω cos L p cos λ p δL p + Ω sin L p sin λ p δλ p ω 13 − Ω sin L p cos λ p δL p − Ω cos L p sin λ p δλ p (47) Since the measurement errors and calculation errors exist, the differential equation of actual computed pseudo velocity can be written as: , ω n pn = ω n pn + δω n pn , g n = g n + δg n , and ∇ b denotes the zero-bias of accelerometer.
Then the pseudo velocity error equation can be obtained as: Ignoring δg n and second order small quantities, substituting Equations (4) and (49) into Equation (50), can obtain: Defining L p = L p + δL p , λ p = λ p + δλ p , h p = h p + δh p , we can obtain the pseudo position error equations from Equation (5). They are can be expressed as:

The Design of Nonlinear Filter Model for Polar Alignment
Aiming at the nonlinear error models derived above, a suitable nonlinear filtering should be selected to replace the standard Kalman filter (KF). Because the extended KF (EKF) has to bear a number of drawbacks such as cumbersome derivation, evaluation of Jacobian matrices, and larger linearized errors and so on [31,32], as a promising substitute for EKF [33], unscented KF (UKF) is employed to accomplish the polar alignment here.
Assuming the measurement errors of accelerometer and gyro are mainly composed of the constant error and zero mean Gaussian white noise. Then we have Consequently, the filtering state with pseudo frame can be set as: where δP p denotes the pseudo position errors, and δP p = δL p δλ p δh p T . According to the nonlinear pseudo error equations in Section 3, the state equation of the filtering model of initial alignment can be obtained as: where 0 3×1 denotes the zero matrix of indicated dimension. Further, the general nonlinear system dynamic model can be written by: .
where the specific expressions of f (X, t) can refer to Equation (54), and W(t) is the white system process noise with zero mean and covariance Q(t).
In addition, the differences of the pseudo velocity projected in g p -frame between the SINS and external navigation device should be selected as the observations of filter here since the velocity information projected in body frame (V b ) would be invalid caused by the smaller horizontal component of rotation of Earth. Then, the measurement equation can also be given by: where H = 0 3×6 I 3×3 0 3×6 is the observation matrix, and V(t) is the white noise with zero mean and covariance R(t).
On the other hand, the block diagram of polar alignment assisted by UKF is also shown in Figure 5. For application of UKF to initial alignment with large azimuth misalignment angle, the specific derivation and implement can be found in [30,31,33]. On the other hand, the block diagram of polar alignment assisted by UKF is also shown in Figure  5. For application of UKF to initial alignment with large azimuth misalignment angle, the specific derivation and implement can be found in [30,31,33].

Simulation and Test
In order to confirm the performance of the proposed polar alignment scheme, both semi-physical static simulation and in-motion tests, which are assisted by UKF under the large azimuth misalignment angle, are conducted in this section.

Semi-Physical Static Simulation
Due to the restriction of geography, the semi-physical static simulation is conducted. As attributed to the rotation of Earth, the real static angular velocity measured by gyro at different geographical positions can be exactly known. The angular velocity measured by gyro can be expressed as:

Simulation and Test
In order to confirm the performance of the proposed polar alignment scheme, both semi-physical static simulation and in-motion tests, which are assisted by UKF under the large azimuth misalignment angle, are conducted in this section.

Semi-Physical Static Simulation
Due to the restriction of geography, the semi-physical static simulation is conducted. As attributed to the rotation of Earth, the real static angular velocity measured by gyro at different geographical positions can be exactly known. The angular velocity measured by gyro can be expressed as: where ω b ib denotes the measured angular velocity; C b g = (C g b ) T and C g b denotes the attitude matrix of IMU relative to g-frame; ω g ie = 0 Ω cos L Ω sin L T , L is the latitude in e-frame. Therefore, the static outputs of gyro in polar regions can be simulated by the measured angular velocity in non-polar regions from Equation (57). The static experimental data are collected with fiber optical gyroscope (FOG) SINS, which fixed in a tri-axis high-precision turntable. The attitude reference of IMU is provided by turntable. The gyro drift and the accelerometer bias of test IMU are less 0.01 • /h, 100 µg respectively. The local longitude and latitude is 126.6778 • and 45.7778 • . The pitch, roll and azimuth of IMU relative to the g-frame are 0 • , 0 • and 30 • respectively. The static outputs of IMU for 1200 s are collected, and the update frequency of IMU is 100 Hz. Moreover, according to Equation (57), the original gyro's outputs and the modified gyro's outputs with 85 • latitude can be shown in Figure 6. Then, the static alignment with the pseudo-Earth frame is carried out, which is assisted by UKF under the condition of large azimuth misalignment angle. In this simulation, the initial pseudo misalignment angles xp  , yp  and zp  are chosen as 1°, 1° and 40°, respectively. Since the position of the vehicle never changes in the static alignment, the strapdown attitude matrix relative g -frame can be obtained directly from Equation (7) without the influence of position errors and the cumbersome computation. As a result, the application of backward process to polar alignment has not been presented here. Through the pseudo attitude matrix and the Equation (7), the static alignment results relative g -frame is obtained directly and is shown in Figure 7. Then, the static alignment with the pseudo-Earth frame is carried out, which is assisted by UKF under the condition of large azimuth misalignment angle. In this simulation, the initial pseudo misalignment angles φ xp , φ yp and φ zp are chosen as 1 • , 1 • and 40 • , respectively. Since the position of the vehicle never changes in the static alignment, the strapdown attitude matrix relative g-frame can be obtained directly from Equation (7) without the influence of position errors and the cumbersome computation. As a result, the application of backward process to polar alignment has not been presented here. Through the pseudo attitude matrix and the Equation (7), the static alignment results relative g-frame is obtained directly and is shown in Figure 7.
As shown in Figure 7, the estimated errors of level misalignment angles and azimuth misalignment angle are convergent with time. The level alignment errors can converge rapidly. However, the convergent time of azimuth error is about 1000 s. Because the horizontal component of rotation of Earth decreases to zero gradually with the increase of latitude, the degree of observability of azimuth would become smaller in polar regions. More time is required for azimuth alignment in polar regions. Moreover, it would also be difficult to achieve self-alignment near the pole. Even so, a better alignment result is still obtained with the proposed pseudo navigation mechanization at 85 • latitude. The alignment error of azimuth is −0.3607 • in 1200 s. Therefore, the proposed pseudo navigation mechanization can solve the problem of traditional attitude reference frame in polar regions caused by SINS mechanization.
Then, the static alignment with the pseudo-Earth frame is carried out, which is assisted by UKF under the condition of large azimuth misalignment angle. In this simulation, the initial pseudo misalignment angles xp  , yp  and zp  are chosen as 1°, 1° and 40°, respectively. Since the position of the vehicle never changes in the static alignment, the strapdown attitude matrix relative g -frame can be obtained directly from Equation (7) without the influence of position errors and the cumbersome computation. As a result, the application of backward process to polar alignment has not been presented here. Through the pseudo attitude matrix and the Equation (7), the static alignment results relative g -frame is obtained directly and is shown in Figure 7. As shown in Figure 7, the estimated errors of level misalignment angles and azimuth misalignment angle are convergent with time. The level alignment errors can converge rapidly. However, the convergent time of azimuth error is about 1000 s. Because the horizontal component of rotation of Earth decreases to zero gradually with the increase of latitude, the degree of observability of azimuth would become smaller in polar regions. More time is required for azimuth alignment in polar regions. Moreover, it would also be difficult to achieve self-alignment near the pole. Even so, a better alignment result is still obtained with the proposed pseudo navigation mechanization at 85° latitude. The alignment error of azimuth is −0.3607° in 1200 s. Therefore, the proposed pseudo navigation mechanization can solve the problem of traditional attitude reference frame in polar regions caused by SINS mechanization.

In-Motion Alignment Tests
Due to the restriction of geography, the actual test data of in-motion alignment are difficult to obtain. Consequently, only a vehicle test of low latitude, which is assisted by UKF under the condition of large azimuth misalignment angle, is firstly conducted to verify the feasibility of the proposed alignment scheme. Then, the polar alignment tests with the simulation data are carried out to validate the superior performance of the proposed method. In addition, in the process of polar alignment, the external velocity observation information projected in the navigation frame (namely, g p -frame) should be employed to fulfill in-motion alignment tests since the horizontal component of rotation of Earth is very small in and near pole regions.
Firstly, the vehicle test of low latitude is executed by equipped with a receiver of differential GPS and two different grades fiber optical gyroscope (FOG) SINSs as shown in Figure 8. The gyro constant drift and the accelerometer constant bias of test IMU are less 0.01 • /h, 100 µg respectively. The local longitude and latitude are 126.6778 • and 45.7778 • , respectively. In total, 250 s vehicle data were collected. The initial pseudo misalignment angles φ xp , φ yp and φ zp are chosen as 1 • , 1 • and 20 • , respectively. In addition, since this paper mainly focuses on the design of polar alignment scheme, the application of backward process to accelerate the alignment process more has not been presented. Only a backward process is conducted to implement polar alignment and accelerate the alignment process here. For the application of backward process to decrease the alignment time further, more details can be found in [27][28][29]. Here, the filtering estimation of initial alignment firstly works with test data by the forward process. Subsequently, the filtering estimation is carried out with the opposite time sequence. The alignment errors and the calculated pseudo position are shown in Figures 9  and 10, respectively.
As shown in Figure 9, the estimated errors of three-axis misalignment angles are convergent with time. The pseudo position can also return to the initial position with the backward process when the initial alignment is completed from Figure 10. Although there are still some differences between the initial value and the final value of pseudo position, and the drift errors of pseudo position are unavoidable because the external velocity information is selected as the observation information of filter in the process of initial alignment, the influence of the drift errors of pseudo position on the initial transformation matrix is very small and can be neglected since the initial pseudo latitude is always zero. Therefore, the pseudo strapdown attitude matrix of vehicle at initial position is obtained at the end of alignment phase. Then, the initial attitude matrix relative to g-frame can be obtained directly with Equation (7). There is less cumbersome derivation and low computational complexity in the process of transformation since the initial transformation matrix is a fixed known constant matrix. Consequently, we can conclude that the proposed polar alignment scheme would be feasible for in-motion alignment and has less complexity. constant drift and the accelerometer constant bias of test IMU are less 0.01°/h, 100 μg respectively. The local longitude and latitude are 126.6778° and 45.7778°, respectively. In total, 250 s vehicle data were collected. The initial pseudo misalignment angles xp  , yp  and zp  are chosen as 1°, 1° and 20°, respectively. In addition, since this paper mainly focuses on the design of polar alignment scheme, the application of backward process to accelerate the alignment process more has not been presented. Only a backward process is conducted to implement polar alignment and accelerate the alignment process here. For the application of backward process to decrease the alignment time further, more details can be found in [27][28][29]. Here, the filtering estimation of initial alignment firstly works with test data by the forward process. Subsequently, the filtering estimation is carried out with the opposite time sequence. The alignment errors and the calculated pseudo position are shown in Figures 9 and 10, respectively.   constant drift and the accelerometer constant bias of test IMU are less 0.01°/h, 100 μg respectively. The local longitude and latitude are 126.6778° and 45.7778°, respectively. In total, 250 s vehicle data were collected. The initial pseudo misalignment angles xp  , yp  and zp  are chosen as 1°, 1° and 20°, respectively. In addition, since this paper mainly focuses on the design of polar alignment scheme, the application of backward process to accelerate the alignment process more has not been presented. Only a backward process is conducted to implement polar alignment and accelerate the alignment process here. For the application of backward process to decrease the alignment time further, more details can be found in [27][28][29]. Here, the filtering estimation of initial alignment firstly works with test data by the forward process. Subsequently, the filtering estimation is carried out with the opposite time sequence. The alignment errors and the calculated pseudo position are shown in Figures 9 and 10, respectively.   As shown in Figure 9, the estimated errors of three-axis misalignment angles are convergent with time. The pseudo position can also return to the initial position with the backward process when the initial alignment is completed from Figure 10. Although there are still some differences between the initial value and the final value of pseudo position, and the drift errors of pseudo position are Next, the superior performance of the proposed method for polar alignment is verified by simulation data due to restriction of geography. In polar in-motion alignment simulations, the outputs of gyro and accelerometer are generated by a strapdown INS simulator. The constant and random drifts of gyro are 0.01 • /h and 0.001 • /h, respectively, and the constant and random biases of the accelerometer are 100 µg and 10 µg, respectively. The initial values of attitude relative to g-frame are all set to be zero, and a constant velocity of 10m/s in the along-vehicle direction is applied all the scenarios. According the definitions of pseudo-Earth frame, the vehicle would sail along the pseudo equator with the constant velocity of 10 m/s. Two groups of in-motion data with the different initial latitudes 87 • and 90 • are obtained. The initial longitudes are all selected as 126 • . The length of test data with the initial latitude 87 • is 600 s, while the length of test data with the initial latitude 90 • is 1200 s. The update frequency of IMU data is 100 Hz. In the process of initial alignment, the initial pseudo position errors L p ,λ p , and h p are all set to be 1 m. The initial pseudo velocity errors are all 0.1 m/s. The initial pseudo misalignment angles φ xp , φ yp , and φ zp are chosen as 1 • , 1 • , and 70 • , respectively. The alignment results of 87 • latitude and 90 • latitude are shown in Figures 11 and 12, respectively. The calculated pseudo positions with the pseudo navigation mechanization are also shown in Figures 13  and 14. In addition, a comparison with the traditional method based on transverse frame at 87 • latitude is also carried out. The final alignment result relative to g-frame based transverse frame is obtained by direct transformation with the calculated transverse position from transverse SINS calculation. The mean and variance of 30 alignment errors in 1200 s are shown in Table 1.         As shown in Figures 11 and 12, at the end of the first forward process, the azimuth alignment errors of 87° and 90° are −13.16° in 600 s and −18.5° in 1200 s, respectively. The alignment errors of azimuth are still larger. While the alignment errors of azimuth are −2.885° and −4.736°, respectively, at the end of the backward process. As a result, the proposed alignment scheme with the backward process extends the length of dada equivalently and can utilize the data adequately to decrease the alignment time. Moreover, the alignment performance of the proposed algorithm is also very favorable. The alignment accuracy can meet the requirement of polar alignment. Therefore, the proposed alignment algorithm can solve the problem of polar alignment caused by SINS mechanization.
On the other hand, in Figures 13 and 14, it is clear that the pseudo position returns to the initial position, when the initial alignment is completed. Although there are still some differences, the influence of the drift errors of pseudo position on the initial transformation matrix   As shown in Figures 11 and 12, at the end of the first forward process, the azimuth alignment errors of 87 • and 90 • are −13.16 • in 600 s and −18.5 • in 1200 s, respectively. The alignment errors of azimuth are still larger. While the alignment errors of azimuth are −2.885 • and −4.736 • , respectively, at the end of the backward process. As a result, the proposed alignment scheme with the backward process extends the length of dada equivalently and can utilize the data adequately to decrease the alignment time. Moreover, the alignment performance of the proposed algorithm is also very favorable. The alignment accuracy can meet the requirement of polar alignment. Therefore, the proposed alignment algorithm can solve the problem of polar alignment caused by SINS mechanization.
On the other hand, in Figures 13 and 14, it is clear that the pseudo position returns to the initial position, when the initial alignment is completed. Although there are still some differences, the influence of the drift errors of pseudo position on the initial transformation matrix C g p g (0) is very small and can be neglected since the initial pseudo latitude is always zero. Therefore, the initial transformation matrix C g p g (0) can be utilized to obtain the strapdown attitude matrix with a simple transformation logic and without influence of position errors. Moreover, with the simple transformation logic, there have less the cumbersome derivation and low computational complexity in the process of transformation, thereby decreasing the complexity of program and the burden of computer. Here, we should notice that the azimuth angle relative to g-frame is the angle between the vehicle direction and the 126 • longitude when the initial latitude is 90 • . In addition, in Table 1, we can see that the mean and variance of alignment errors based on the proposed algorithm are also smaller as compared with the traditional method. This is because the proposed algorithm eliminates the influence of position errors on transformational accuracy. Therefore, the proposed algorithm has superior performance. On the other hand, since the strapdown attitude matrix relative to g-frame is obtained at the end of alignment phase, it would also be more convenient to access any one of the following polar navigation systems without any additional workloads and complexities. Because the existing transformation program of navigation parameters in global navigation system can be utilized directly, the proposed polar alignment algorithm would be advantageous to unify the forms of polar alignment and the external observation information of filter.
Furthermore, comparing Figures 11 and 12, more time is required for polar alignment with the initial latitude 90 • . This is the other key problem, which must be confronted in polar alignment. The degree of observability of azimuth is smaller in polar regions since the horizontal component of rotation of Earth decreases to zero gradually with the increase of latitude. As a result, the future efforts will focus on the improvement of the degree of observability of azimuth, thereby improving the polar alignment performance more.

Conclusions
In this paper, from the perspective of global inertial navigation, a generalized transverse Earth frame built based on the initial position of vehicle (namely, pseudo-Earth frame) is proposed as the attitude reference frame for polar alignment. Then, a polar alignment algorithm based on the pseudo-Earth frame and backward process is proposed to implement polar alignment. Through the proposed algorithm, the attitude matrix relative to local geographic frame is obtained without influence of position errors when the initial alignment is completed, and there is less cumbersome derivation and low computational complexity. As a result, it is very convenient to access any one of the following polar navigation systems without any additional workloads and complexities. Moreover, it would also be advantageous to unify the forms of polar alignment and the external observation information of filter, thereby further unifying the implemented form of the reference information acquired from the external navigation device. In addition, the backward process can still accelerate the filtering estimation and decrease the alignment time. Finally, the semi-physical static simulation and in-motion tests with large azimuth misalignment angle assisted by UKF are carried out to verify the proposed polar alignment algorithm, and the results exhibit the superior performance of the proposed alignment algorithm in polar regions.