Automatic Estimation of Dynamic Lever Arms for a Position and Orientation System

An inertially stabilized platform (ISP) is generally equipped with a position and orientation system (POS) to isolate attitude disturbances and to focus surveying sensors on interesting targets. However, rotation of the ISP will result in a time-varying lever arm between the measuring center of the inertial measurement unit (IMU) and the phase center of the Global Positioning System (GPS) antenna, making it difficult to measure and provide compensation. To avoid the complexity of manual measurement and improve surveying efficiency, we propose an automatic estimation method for the dynamic lever arm. With the aid of the ISP encoder data, we decompose the variable lever arm into two constant lever arms to be estimated on line. With a complete 21-dimensional state Kalman filter, we accurately and simultaneously accomplish navigation and dynamic lever arm calibration. Our observability analysis provides a valuable insight into the conditions under which the lever arms can be estimated, and we use the error distribution method to reveal which error sources are the most influential. The simulation results demonstrate that the dynamic lever arm can be estimated to within [0.0104; 0.0110; 0.0178] m, an accuracy that is equivalent to the positioning accuracy of Carrier-phase Differential GPS (CDGPS).


Introduction
Airborne earth observation takes the aircraft as the platform, and it uses remote sensing load such as synthesize aperture radar (SAR) and charge coupled device (CCD) array camera to acquire a wide range, highly accurate and multilayered space-time information of global surface and deep earth [1]. A position and orientation system (POS) is a dedicated strapdown inertial navigation system (SINS)/Global Positioning System (GPS) integrated system for airborne remote sensing [2]. As a direct georeferencing (DG) system, the POS provides the ability to directly relate the data collected by a remote sensing system to earth without using traditional ground-based measurements [3]. To improve mapping accuracy and efficiency, the POS is required to provide the position information to centimeter-level accuracy, and orientation data to sub-arcminute accuracy in either real-time navigation or post-mission, such as the current state-of-the-art airborne system POS AV610 [3,4]. This is primarily accomplished by the inertial measurement unit (IMU) integrated with a Carrier-phase Differential GPS (CDGPS).
The key algorithm for POS is still SINS/GPS integrated navigation [5,6]. For the loosely-coupled integration method discussed here, the positions and velocities of the SINS and GPS solutions are compared, and the resulting differences are used as the measurement values for a Kalman filter [7,8]. Due to the extremely high positioning accuracy of a CDGPS, typically 1-2 cm plus 1 ppm of baseline separation in real-time kinematic (RTK) or post-processing [9,10], many overlooked errors must be treated seriously. Among these, the time and space asynchronies are the dominant error sources. The time asynchrony means the time delay of the GPS position and the velocity information with respect to the SINS data in real-time systems, which will badly influence the accuracy of the SINS/GPS navigation solution. Based on one pulse per second (1PPS) of the GPS, which can deliver very high time accuracies such as 20 ns, the time delay can be significantly attenuated by the professional time synchronization technology [11]. The space asynchrony means the lever arm between the GPS antenna and the IMU. Generally, the GPS antenna is mounted on the roof of the plane for better GPS satellite visibility [7,12], while the IMU is installed inside the plane cabin or under the aircraft belly as close to the surveying sensors as possible to weaken the influence of flexural deformation. Therefore, the phase center of the GPS antenna and the measuring center of the IMU are separated by a lever arm. Since the lever arm in length is usually several meters, it severely deteriorates measurement values under angular excitations if no compensation is provided.
Many efforts have been devoted to improving SINS/GPS integration accuracy by attenuating the lever arm uncertainty. One method is manual measurement by a total station [13,14], which can obtain the length of the lever arm with sufficient accuracy, but suffers from complicated operation. Furthermore, it is a difficult task in itself to determine the phase center of a GPS antenna and the sensing center of an IMU. Repeated measurements are inevitable if the GPS antenna or IMU is replaced or remounted. The other method is automatic identification of the lever arm by the SINS/GPS integration navigation [1,[15][16][17][18][19][20], in which the lever arm is extended as an error state of the Kalman filter to be estimated, while the vehicle maneuvers. The latter method is simple and effective, but practical only if the lever arm is constant.
To maintain imaging quality, aerial surveying sensors are required to move in a straight line with a given direction, but a plane may deviate from the ideal motion trajectory, due to gusts, air turbulence, or flight control errors [7,21]. Consequently, an inertially stabilized platform (ISP) with three gimbals is usually used to isolate the angular motion interference and to focus the surveying sensor onto the desired targets, such as the general ISP PAV80 and GSM3000 [7,22,23]. The use of ISP makes the lever arm problem more complicated. Generally, the surveying sensor and the IMU are rigidly mounted close together on the ISP, but the measuring center of the IMU cannot be ensured to coincide with the rotating center of the ISP. When the ISP rotates, the lever arm between the IMU center and the GPS antenna center changes; for this reason, we call it a dynamic lever arm. Since the lever arm is not constant in this situation, traditional lever arm estimation methods do not work. In [7,12], the dynamic lever arm is decomposed into two constant lever arms to be compensated in flying missions, but they are still manually measured in advance. In order to avoid complex manual measurements and to improve surveying efficiency, we propose an automatic approach to estimating two separate lever arms that provide equivalent compensation.
The rest of this paper is organized as follows: Section 2 describes the system model and the problems to be solved. Section 3 presents the proposed integrated navigation approach and analyzes the corresponding observability in detail. Section 4 verifies our method through simulation tests, and conclusions are drawn in Section 5.

System Description
As shown in Figure 1, the IMU and surveying sensors are usually fixed on the inner gimbal of the ISP. Thus, the mounting angle and lever arm between the two are constant, reducing the accuracy requirements for the ISP. The GPS antenna is mounted on the roof of the plane for better satellite visibility. The related coordinate systems are defined as follows: n-frame, the navigation frame, chosen as the local level east-north-up (ENU) coordinate; b-frame, the IMU body frame, implicitly predefined by the calibrated sensitive axes of the inertial sensors, with the origin located at the sensitive center of IMU (point O I ), and axes pointing right, forward, and upward; b c -frame, the carrier frame, fixed to the flight vehicle and originating at the rotating center of the ISP (point O O ). The axis directions of b c -frame do not rotate with the gimbals, and coincide with the  ISP only if the encoder data of three gimbals are all zeros; b s -frame, the surveying sensor frame, rigidly fixed to the surveying sensor with the origin located at the surveying center (point O S ).
Based on the established frames, the related lever arms are defined as follows:     L is uncertain and time-varying, it is referred to as a dynamic lever arm.
Since the attitude matrix n b C can be obtained from the SINS solution, the transformation matrix c n b C from c b -frame to the n-frame can be obtained as follows: in which the direction cosine matrix C is formed with the encoder data of the ISP [7]: where: Three lever arms-L b OI , L b c OG , and L b IS -are constant, while lever arm L b c IG changes with the ISP rotation. Because the lever arm L b c IG is uncertain and time-varying, it is referred to as a dynamic lever arm.
Since the attitude matrix C n b can be obtained from the SINS solution, the transformation matrix C n b c from b c -frame to the n-frame can be obtained as follows: in which the direction cosine matrix C b b c is formed with the encoder data of the ISP [7]: where: the angles ψ T , θ T , and γ T are encoder data obtained from the inner, middle, and outer gimbals of the ISP, respectively.
Referring to Equation (2), the angular rate of the b-frame with respect to the b c -frame, denoted as ω b b c b , can be calculated as follows [7]: where . θ T , . γ T , and . ψ T are the angular rates of roll, pitch, and heading gimbals of the ISP, respectively, which are obtained by the time differentiation of the corresponding encoder data.
According to the lever arm effect [24], if the position and velocity parameters at point O I are obtained from navigation solutions, and the lever arm L b IS is known, the corresponding information at point O S can be calculated as follows: where p S and p I denote the position information consisting of latitude, longitude, altitude; v n S and v n I are the velocity data at point O S and point O I in the n-frame; L and h are the local latitude and height; R N and R M denote the transverse and meridian radius of curvature. The relative angular rate where v n O and v n G denote the velocities at point O O and point O G in the n-frame, and the relative angular rate ω b c eb c can be obtained by: In a normal integrated navigation mechanism, the velocity and/or position differences between the SINS and GPS solutions are generally taken as the measurement values. To improve the measurement accuracy, the dynamic lever arm L b c IG must be compensated precisely even if it is difficult. Based on Equations (7) and (8), the dynamic lever arm can be compensated as follows: Combining Equations (9) and (11) yields: This means that the dynamic lever arm L b c IG can be decomposed into two constant lever arms L b c OG and L b OI , because of the fact that the ISP rotating center O O is fixed to both the b-frame and the b c -frame. From the geometric location relationship shown in Figure 1, we can deduce the same results as Equation (12). Therefore, if the two constant lever arms L b c OG and L b OI are precisely measured in advance, the dynamic lever arm L b c IG can be compensated as Equation (12). However, as mentioned before, manual measurement is complex, tedious, and inefficient. To avoid these drawbacks and enhance the efficiency of surveying, we propose a method to automatically estimate the lever arm and to accomplish the integrated navigation simultaneously.

SINS/GPS Integrated Navigation Model
As shown in Equation (11), if the SINS/GPS integrated navigation model can be designed to estimate the lever arms L b c OG and L b OI , the complex dynamic lever arm L b c IG can be compensated automatically, and the tedious work of manual measurement can be eliminated. Combining the classical error model of SINS and the constant lever arms L b OI and L b c OG , we define a 21-dimensional error state vector of SINS/GPS integration as: .
where F I is a 15 × 15 transition matrix based on typical SINS error model as detailed in [25]. The symbols ε b w and ∇ b w are the noises of gyros and accelerometers, and 0 i×j denotes an i × j zero matrix. For loosely coupled SINS/GPS integrated navigation, there are two measurement models to be selected in practical POS, including position measurements and position/velocity measurements. The associated measurement equations are: where w gp and w gv represent the position and velocity noises of the GPS. Since Equation (11) can be rewritten as: In this paper, we depict the position measurement equation as an example, while the position/velocity measurement equations have the same results.

Observability Analysis for Lever Arms
Usually, an observability analysis for an integrated navigation system is conducted to determine whether a specific error state is observable, and under what conditions it can be estimated [26]. According to global observability theory [27], a state is said to be observable if it can be determined uniquely from the measurements for a finite time and under a possible condition. Therefore, we investigate the observability of lever arms directly, depending on the measurement equation in this section.
The position measurement in Equation (15) can be rewritten as an error model: Clearly, this is not a deterministic system, but a stochastic case. However, the related theoretical analysis is still useful for revealing the observable conditions and guiding the design of verification trajectories. Ignoring the driving noise and measurement noise temporarily, Equation (18) itself determines the observability of the associated errors δp I , L b OI , and L b c OG . It can be seen that the acceleration or deceleration excitations have little effect on lever arms L b OI and L b c OG , and their observabilities are determined by direction cosine matrices C n b c and C b c b . If only the ISP rotates, C b c b changes while C n b c remains constant; thus, the lever arm L b OI will be distinguished. In this situation, the SINS position error cannot be accurately estimated for the unobserved lever arm L b c OG . Considering the construction of attitude matrix C b c b , the rotation of one ISP gimbal can reveal two components of the lever arm L b OI perpendicular to the rotating axis. Consequently, rotations of at least two gimbals are needed to estimate all three components.
If only the attitude of the vehicle changes, C n b c varies while C b c b stays the same, and the combined lever arm OI can be estimated. Thus, position measurement accuracy can be improved in spite of the unseparated lever arm L b OI . It can be seen that the change in C n b c is crucial to determining the final positioning accuracies. Similarly, improvement of all three dimensional position accuracies depends on attitude changing in at least two directions of the vehicle.
In summary, the estimate of constant lever arms L b OI and L b c OG requires the change of direction cosine matrices C b c b and C n b c . The greater the change of the matrixes, the stronger the observability of the corresponding lever arms. Certainly, the final estimation accuracies of associated lever arms are limited by the measurement error noise w gp .
In aerial mapping missions, the aircraft usually maneuvers intentionally before entering the sensing area to enhance the observability of the POS [17,20], especially to improve the heading accuracy. During the mapping process, long straight portions of the flight will cause heading error growth; thus, periodical aircraft maneuvers are required to separate the navigation errors from inertial sensor errors [18]. Common flight maneuvers include circling, coordinated turn, or "8" shaped flight, in which at least two directions of attitude change occur, such as in the heading and roll axes. Thus, in those cases without ISP, the lever arm L b c IG is usually treated as a constant to be estimated in these planned maneuver processes, and satisfactory results can be obtained [1,17]. As long as any two gimbals of ISP are controlled to rotate during these auxiliary maneuvering processes, the two constant lever arms L b OI and L b c OG are solvable. Therefore, the observable conditions for the dynamic lever arm estimation are easy to satisfy in practical applications.
It should be noted that the auxiliary matrix C b c b is constructed from the encoder data of the ISP, as shown in Equation (2). Since there are some unavoidable errors in the encoder angle measurements, their influences should be evaluated. Denoting the angle measurement errors as µ b c , the introduced position errors are deduced as follows: The term (L b c OG × µ b c ) determines the magnitude of the position errors that are introduced by encoder errors. Assuming that the lever arm L b c OG is 2 m and the encoder angle errors of the ISP are 3 arcminutes, the related position errors will be at the 2 mm level. The commonly used ISP PAV80 [22] and GSM3000 [23] can provide such angular accuracy. In this situation, the influence of encoder angle errors is limited and much smaller than the positioning accuracy of typical CDGPS. However, if the lever arm L b c OG is too large, the interrelated influence will not be negligible. The residual encoder error of the ISP can be modelled as a zero mean white noises process, which has the same effect on automatic estimation method and manual measurement method of lever arm. Thus, the GPS antenna and the ISP should be installed as close as possible, otherwise higher performance ISP is needed.
The time delay of the GPS information relative to SINS is another major error source. Since the GPS receiver and the SINS unit are two separate subsystems, the frequency difference between the 1PPS and SINS sampling clock could cause the time delay to be a variable (usually a ramp or triangular wave). Thus, the method of taking the time delay as a constant error state to be estimated will not work. The effective way is to discipline the SINS sampling clock with 1PPS [28], or to record the delayed time with a counter for software compensation [11]. Even so, it is difficult for the time delay to be perfectly compensated. In addition, although the same sampling clock is used in the IMU, the time-asynchrony between each inertial sensor is inevitable, due to the differences of physical characteristics [29], which makes the time synchronization more complicated. Denoting the time-asynchrony error as δt, its influence on position measurement is: Assuming that δt is 0.1 ms and the normal flight velocity is 50 m/s, the introduced position error is about 5 mm. Clearly, larger time-asynchrony error is intolerable at this speed. After time compensation as mentioned earlier, the residual time-asynchrony error can be modelled as a zero mean white noises process for analysis.

Error Distribution Analysis
The purpose of error distribution analysis for stochastic systems is to evaluate the final estimation accuracy of the concerned states, and to distinguish which error sources primarily contribute to the concerned estimation uncertainties. Monte-Carlo simulations are appropriate to determine error distributions, but suffer from a large computational burden. Based on the law of large numbers [30], Monte-Carlo simulations necessitate a large number of single-factor simulation tests and very long times to establish the error distributions if many error sources are involved. Considering linear error propagation, covariance simulation programs [31,32] are commonly used to provide numerical time histories depicting the accuracy of a given system configuration in terms of the covariance of its associated linearized error state vector, avoiding large numbers of simulations. According to the established SINS/GPS Kalman filter, we designed a covariance analysis method to get the error distribution budgets, especially focusing on which error sources contribute to the estimation uncertainty of lever arms and the final position errors.
Three system models are involved in the covariance analysis method: the real world model, the full-order design filter, and the reduced-order filter model [33]. The real world model, also known as true model, describes the behavior of the actual system, to the best of one's knowledge and ability. The full-order filter is dedicated to yielding the estimate for all state errors of the real world. The reduced-order filter model is designed for filtering calculation in actual system, in which some weakly observable states have been excluded from the full-order filter to ease computational burden.
In this paper, the real world model is defined in Appendix A, the corresponding full-order filter is depicted in Appendix B, and the reduced-order filter model is designed in Section 3.1. The state vector of the reduced-order filter is composed of the first 21 elements of the full-order filter state. In addition, the full-order filter model contains three types of measurement noise, the GPS positioning noise, the ISP encoder noise and the time-asynchrony noise, and only the first type is selected as the measurement noise of the reduced-order filter model.
For the reduced-order filter model, after the state equation and measurement equation are established and discretized, we update the standard linear Kalman model in two calculation loops [34]: the state filtering loop and covariance calculation loop. For ease of depiction, the discretized state filtering loop is listed here:x while the covariance calculation loop is updated as follows: where the hat label "ˆ" denotes the calculated value in the reduced-order filter model, distinguishing it from that of the full-order filter model; the subscript k − 1 and k denote the last and current cycle, while k/k − 1 expresses the propagation value;Φ k,k−1 ,Γ k−1 andĤ k are discretized values from the state equation and measurement equation;Q k−1 andR k are the previously set covariance values of process noises and measurement noises;P k is the current covariance matrix recursively calculated from the initial valueP 0 ; andK k is the calculated filtering gain matrix. It can be seen that the covariance calculation loop is independent of the state filtering loop. Once the filter is designed and the trajectory is selected, the associated matricesΦ k,k−1 ,P 0 ,Q k−1 ,R k andĤ k are uniquely determined, so is the gain matrixK k . The filtering program only needs to be run once, so that the gain matrixK k can be obtained at every measurement update. It should be recorded and substituted into the full-order filter model for covariance analysis.
In the full-order filter model, the state equation and measurement equation are extended as follows: x where the superscript "*" denotes the matrix involved in the full-order filter model. It can be seen that there are three types of errors that affect the estimation results: the initial state errors x * 0 , the process noises W * k−1 , and the measurement noises V * k . Based on the linear error model, the estimated error of a concerned state at any moment can be expressed as a linear cumulative sum: where S * is the propagation matrix of initial errors; n, l, and m are the nonzero element numbers of the initial state, process noise, and measurement noise, respectively; x * i,W j or x * i,V j denotes the estimated error of ith state caused by the jth process noise W j or the jth measurement noise V j . For the sake of simplicity, the current cycle k is omitted in Equation (28).
The filtering gain matrixK k is the crucial link between the reduced-order filter model and the full-order filter model, which can be extended into the full-order filter model as follows: Since the gain matrix K * k of full-order filter model is known, the covariance calculation loop is updated as follows: Firstly, based on linear propagation characteristics of associated errors, and setting Q * = 0 and R * = 0 in Equations (30) and (31), the mean square deviation propagation matrix of initial errors can be calculated as follows: with the initial value S * 0 = I. Secondly, setting P * 0 = 0 and R * = 0 in Equations (30) and (31), the propagation matrix of the jth process noise W j is obtained by: with the initial value A j 0 = 0, where e jj denotes the 39 × 39 matrix with e(j, j) = 1 and all else zeros. If there are l-dimensional nonzero process noises, Equations (33) and (34) should be executed l times with different e jj to obtain all of the propagation matrices.
Thirdly, setting P * 0 = 0 and Q * = 0 in Equations (30) and (31), the propagation matrix of the jth measurement noise V j is obtained by: with the initial value B j 0 = 0 and m calculations required for m-dimensional measurement noises. Therefore, the estimation accuracy of x * i in Equation (28) can be evaluated by the corresponding mean square deviation σ xi : where σ 2 x * j0 , Q jj , and R jj are the deviations of the jth initial state error, jth process noise, and jth measurement noise, respectively, while A j ii and B j ii are the ith diagonal elements of covariance matrices A j k and B j k , respectively. Based on the linear error model, the mean square deviation σ xi,xj of linear combination of x * i and x * j can be expressed as σ 2 xi + σ 2 xj , and so on. Therefore, according to the above covariance analyzing method, the error distribution budgets of the full order filter or real-world model can be constructed from the calculation of the reduced order filter. The complexity and computational burden of simulation are greatly reduced.

Simulation Conditions
Considering the extremely high accuracy requirements for a lever arm in POS, it is difficult to find another more accurate benchmark to evaluate the estimation results. As an expedient, we carried out a series of simulation tests instead of real tests to verify the proposed method. Considering the effect of the time-asynchrony error, we assume that the aircraft flies at 50 m/s uniform speed, and the angular motion excitations are added into the simulation trajectory. The initial values needed for the trajectory are given as: To verify the observability analysis results in Section 3.2, two simulation trajectories with different rotation sequences are designed as shown in Table 1. The total time of one simulation trajectory is 150 s, and the related rotation angular rates in two trajectories are all set to 10 • /s. As shown in Table 1, the only difference between the two trajectories is that steps 2 and 3 are exchanged, which means that the change of matrices C b c b and C n b c occur in a different time sequence. The rotation angles of ISP gimbals and vehicular attitudes are shown in Figure 2.
Referring to the system POS AV610 produced by Applanix [4,35], we simulated an aviation-grade SINS and a CDGPS receiver with sensor specifications, including initial errors, process noises, and measurement noises listed in Table 2. Without loss of generality, the IMU data is sampled at 100 Hz, the ISP encoder data outputs at 20 Hz, and the GPS positioning information at 1 Hz.  Referring to the system POS AV610 produced by Applanix [4,35], we simulated an aviation-grade SINS and a CDGPS receiver with sensor specifications, including initial errors, process noises, and measurement noises listed in Table 2. Without loss of generality, the IMU data is sampled at 100 Hz, the ISP encoder data outputs at 20 Hz, and the GPS positioning information at 1 Hz.

Simulation Results of Lever Arm Estimation
We use the reduced-order filter model with parameters set as in Appendix C to estimate the lever arms and position errors in two different trajectories. The simulation results of Trajectory 1 are illustrated in Figure 3, while the simulation results of Trajectory 2 are given in Figure 4.

Simulation Results of Lever Arm Estimation
We use the reduced-order filter model with parameters set as in Appendix C to estimate the lever arms and position errors in two different trajectories. The simulation results of Trajectory 1 are illustrated in Figure 3, while the simulation results of Trajectory 2 are given in Figure 4. Figure 3a shows the estimation errors of lever arm L b OI , Figure 3b,c illustrate the estimation errors of lever arm L b c OG and dynamic lever arm L b c IG , while Figure 3d gives the position errors of SINS/GPS integrated navigation. As shown in Figure 3a, the heading gimbal of ISP rotates by 90 • at 30 s, and the x and y components of lever arm L b OI are estimated correctly. When the roll gimbal of the ISP rotates by 30 • at 60 s, the remaining component L bz OI is identified quickly. It can be seen that the rotations of ISP gimbals can distinguish the lever arm L b OI , but have no effect on lever arm L b c OG . Therefore, there are still significant position errors in integrated navigation before 90 s as shown in Figure 3d. Once the vehicular attitude changes at 90 s and 120 s, the combined lever arm OI is estimated as in Figure 3c and the position errors converge immediately as shown in Figure 3d. Since the lever arm L b OI is identified, the constant lever arm L b c OG is determined with the convergence of lever arm L b c IG .    To evaluate lever arm estimation accuracy and positioning accuracy, we carried out 100 Monte-Carlo simulations, 50 times for each trajectory. The final lever arm estimation results at the end of every test are shown in Figure 5, and Figure 6 shows the final positioning results. In these Figures, the terms "Trj1" is the abbreviation for Trajectory 1, "Trj2" is the abbreviation for Trajectory 2, As shown in Figure 4c, the combined lever arm L b IG is identified when the vehicular heading turns at 30 s and the roll rotates at 60 s. Clearly, the change of attitude matrix C n b c is helpful to estimate the equivalent lever arm L b IG . Even though the constant lever arms L b OI and L b c OG cannot be separated, it does not affect the positioning accuracy of SINS/GPS, as shown in Figure 4d. As the ISP gimbals begin to rotate at 90 s and 120 s, the lever arms L b OI and L b c OG are further separated, as shown in Figure 4a,b. To evaluate lever arm estimation accuracy and positioning accuracy, we carried out 100 Monte-Carlo simulations, 50 times for each trajectory. The final lever arm estimation results at the end of every test are shown in Figure 5, and Figure 6 shows the final positioning results. In these Figures, the terms "Trj1" is the abbreviation for Trajectory 1, "Trj2" is the abbreviation for Trajectory 2, and "RMS" are the envelopes of root mean square error statistics. To evaluate lever arm estimation accuracy and positioning accuracy, we carried out 100 Monte-Carlo simulations, 50 times for each trajectory. The final lever arm estimation results at the end of every test are shown in Figure 5, and Figure 6 shows the final positioning results. In these Figures, the terms "Trj1" is the abbreviation for Trajectory 1, "Trj2" is the abbreviation for Trajectory 2, and "RMS" are the envelopes of root mean square error statistics.  Table 2.
In manual measurements, the lever arms can be measured with an accuracy of 1 mm level by a laser total station instrument [17], if the phase center of the GPS antenna, the sensing center of the IMU, and the rotating center of the ISP are all marked out without error. It is true that the automatic estimation method presented here cannot accomplish such accuracy, which is restricted by the accuracy of CDGPS. However, in many practical situations, it is difficult to accurately locate these three centers, especially for the IMU. As the sensitive centers of the three-axis accelerometers do not coincide, the IMU actually have a virtual sensing center determined by the calibration algorithm of the inner arms [36], and the calibration accuracy is greatly influenced by the accuracy of inertial  Table 2.
In manual measurements, the lever arms can be measured with an accuracy of 1 mm level by a laser total station instrument [17], if the phase center of the GPS antenna, the sensing center of the IMU, and the rotating center of the ISP are all marked out without error. It is true that the automatic estimation method presented here cannot accomplish such accuracy, which is restricted by the accuracy of CDGPS. However, in many practical situations, it is difficult to accurately locate these three centers, especially for the IMU. As the sensitive centers of the three-axis accelerometers do not coincide, the IMU actually have a virtual sensing center determined by the calibration algorithm of the inner arms [36], and the calibration accuracy is greatly influenced by the accuracy of inertial sensors. If the sensing centers cannot be marked out physically, the accuracy of manual measuring cannot be guaranteed.
The purpose of the lever arm compensation is to obtain a higher positioning accuracy of SINS/GPS integration. Assuming that the lever arm has been manually measured with an accuracy of 1 mm, under the same simulation conditions, we obtain the positioning results of the SINS/GPS integration as shown in Figure 7.
As shown in Figures 5 and 6, the final estimation accuracy for the equivalent dynamic lever arm b IG L is [0.0104; 0.0110; 0.0178] m (RMS), and the positioning accuracy is [0.0123; 0.0128; 0.0181] m (RMS), which is nearly the same as the positioning accuracy of CDGPS listed in Table 2.
In manual measurements, the lever arms can be measured with an accuracy of 1 mm level by a laser total station instrument [17], if the phase center of the GPS antenna, the sensing center of the IMU, and the rotating center of the ISP are all marked out without error. It is true that the automatic estimation method presented here cannot accomplish such accuracy, which is restricted by the accuracy of CDGPS. However, in many practical situations, it is difficult to accurately locate these three centers, especially for the IMU. As the sensitive centers of the three-axis accelerometers do not coincide, the IMU actually have a virtual sensing center determined by the calibration algorithm of the inner arms [36], and the calibration accuracy is greatly influenced by the accuracy of inertial sensors. If the sensing centers cannot be marked out physically, the accuracy of manual measuring cannot be guaranteed.
The purpose of the lever arm compensation is to obtain a higher positioning accuracy of SINS/GPS integration. Assuming that the lever arm has been manually measured with an accuracy of 1 mm, under the same simulation conditions, we obtain the positioning results of the SINS/GPS integration as shown in Figure 7.  The positioning accuracy of Figure 7 is [0.0071; 0.0068; 0.0078] m (RMS), which is slightly higher than that of the automatic estimation method. Even so, compared with the state-of-art product, POS AV610 [4], which can accomplish typical positioning accuracy of 0.02 m (horizontal)/0.05 m (vertical) by post-processing, the accuracy of the proposed automatic estimation method is still acceptable. Thus, the automatic estimation method is an effective alternative approach for certain situations, in which the measuring conditions are not available or the dynamic lever arm cannot be measured accurately and efficiently.

Error Distribution Simulation Results
The error distribution simulation tests are devoted to determining which error sources primarily affect the estimation accuracies of specific parameters. For sake of simplicity, the associated error sources in full-order model are classified into eight categories, as shown in Table 3, with the same RMS values as Table 2. The noteworthy parameters are chosen as the norm of lever arm L b OI , with RMS error denoted as σ L b OI ; the norm of lever arm L b c OG with RMS error denoted as σ L b c OG ; and the norm of positioning error with RMS error denoted as σ|δp|. Under the Trajectory 1, the effects of each category on the estimation accuracies are illustrated in Figure 8.    As shown in Figure 8, the initial errors and process noises have little effect on the positioning and lever arm estimation accuracy, while the final positioning accuracy are mainly determined by As shown in Figure 8, the initial errors and process noises have little effect on the positioning and lever arm estimation accuracy, while the final positioning accuracy are mainly determined by the CDGPS error, consistent with the previous analysis. In addition, the time-asynchrony error and ISP encoder error have important influences on the final accuracy, and thus should be treated seriously. We note that the results of Figure 8 are based on input error values shown in Table 2. According to the linear error model, if the input errors expand n times, their effects will also increase n times.

Conclusions
The ISP equipped in a POS results in a variable lever arm, making it difficult to compensate. Instead of inefficient manual measurement strategies, our proposed method uses the encoder data from ISP gimbals to decompose the dynamic lever arm into two constant lever arms to be estimated automatically. Based on the established 21-dimensional state Kalman filter, our error analysis and simulation tests provide some meaningful conclusions. First, rotation of the ISP and an attitude change of the vehicle are necessary to determine the dynamic lever arm. Second, encoder errors of the ISP have a limited effect on estimation accuracy, so that general turntables can be accepted if the lever arm L b c OG is small enough. Third, precise time compensation is necessitated in both the lever arm estimation method and the manual measuring method. Last, the estimation accuracy of the dynamic lever arm mainly depends on the positioning accuracy of the CDGPS. Therefore, the proposed method can be used as an effective substitute for manual measurement with acceptable accuracy.
Based on the real world model, the state vector of the full-order filter is selected as: where the superscript "*" denotes the matrix involved in full-order filter model. The accompanying state equation of the full-order filter is designed as: where: .