Zero-Velocity Update-Based GNSS/IMU Tightly Coupled Algorithm with the Constraint of the Earth’s Rotation Angular Velocity for Cableway Bracket Deformation Monitoring

Cableways have been widely used in industrial areas, cities, and scenic spots due to their advantages, such as being a convenient mode of transportation, time-saving, labor-saving, and low cost, as well as offering environmental protection. To ensure the safe operation of a cableway, based on the characteristic that the velocity of the cableway bracket is approximately zero in a static deformation monitoring environment, a deformation monitoring method called zero velocity update (ZUPT)-based GNSS/IMU tightly coupled algorithm with the constraint of the Earth’s rotation angular velocity was proposed. The proposed method can effectively solve the problem of a single GNSS being unable to output attitude, which is directly related to the status of wire ropes and cable cars. Meanwhile, ZUPT is used to restrain the Kalman filter’s divergence when IMU is stationary. However, the improvements of ZUPT on attitude are not obvious, so the constraint of the Earth’s rotation angular velocity was applied. The performance of the proposed method was evaluated through monitoring the cableway bracket of the Yimeng Mountain Tourism area in Shandong. Compared with the ZUPT-based GNSS/IMU tightly coupled algorithm (ZUPT-TC), the proposed method can further constrain the error accumulation of IMU while stationary and, therefore, it can provide reliable position and attitude information on cableway brackets.


Introduction
Cableways are transportation tools that have functions such as freight and passenger transportation.Compared with other transportation equipment, a cableway has the advantages of being a convenient mode of transportation, time-saving, labor-saving, and low cost, as well as offering environmental protection, and it is therefore widely used in tourism [1,2].Nowadays, the number of in-use passenger cableways in China has reached 1148, and more than 4000 patents for passenger cableways have been disclosed.However, due to the relatively late development and incomplete technology of cableway safety monitoring in China, accidents related to cableways are constantly occurring in the Yunnan Xishuangbanna Cableway, Hunan Tianmenshan Cableway, Shanxi Huashan Cableway, Nanjing Zijinshan Cableway, and other areas.In this case, safety monitoring technology is a must to ensure the safe operation of the cableway [3].Generally, the support device of the cableway is composed of a steel wire rope and a cableway bracket.It is acknowledged that the occurrence of cableway accidents is often due to damage to the steel wire rope, and the cableway bracket is used to support the normal operation of the steel wire rope, so potential faults and abnormal states of the cableway can be detected in a timely manner by monitoring the three-dimensional deformation and attitude changes in the cableway bracket [4].Therefore, it is imperative to monitor the deformation and attitude changes in the cableway bracket.
Periodic contact measurement is a common method used in monitoring the deformation and attitude changes in cableway brackets, but it has poor timeliness and low efficiency.In addition, periodic contact measurement cannot reflect the operational status of the cableway in real time.In order to improve the accuracy and real-time performance of monitoring, dual-or multi-antenna GNSSs can be used to obtain the deformation and attitude changes in the cableway bracket [5,6].Bock proposed that GPS can be used to monitor surface deformation [7].Subsequently, GPS technology was applied to deformation monitoring of natural gas pipelines and achieved sub-centimeter-level monitoring accuracy [8].In [9], GPS was found to be able to monitor building vibration frequencies of 1.2 Hz, and it was proven that the use of GNSS for deformation monitoring is actually feasible.What is more, GNSS was also successfully applied in bridge structural deformation and deflection monitoring [10], dam deformation monitoring [11], and volcanic turbulence and post-earthquake deformation monitoring [12].Unfortunately, high-precision attitude is difficult to obtain directly through GNSS due to the structural defects in the design of the cableway bracket.
The inertial measurement unit (IMU) is based on Newton's mechanical laws of relative inertial space, using inertial components such as gyroscopes and accelerometers to sense the acceleration and angular velocity of the carrier during motion.Then, the position, velocity, and attitude navigation information of the carrier can be obtained by integral calculation.The IMU has the advantages of strong independence, anti-interference ability, and high sampling frequency, and is widely used in integrated navigation [13][14][15].However, extended time sampling of INS in a stationary state will lead to drifting error [16].Therefore, the outstanding short-term relative measurement accuracy of INS in integrated navigation applications was studied and analyzed [17][18][19].The GNSS/IMU integrated method can fully utilize the accuracy of GNSS positioning, as well as the high sampling frequency and attitude acquisition ability of IMU [20], so it can provide accurate monitoring data for cableway operation.In the environment of dynamic GNSS/IMU deformation monitoring, the feasibility of integrated navigation for track deformation monitoring has been verified, which means that integrated navigation can perform high-precision deformation monitoring [21][22][23].Furthermore, road flatness measurement and track irregularity measurement can also be carried out through integrated navigation [24].Different from dynamic measurement, GNSS/IMU static measurement has attracted great attention because of its unique properties in deformation monitoring.In [25], it was proven that integrated positioning can monitor high-frequency small deformation (such as seismic waves) and meet the accuracy requirements, which also provides theoretical preparation for the GNSS/IMU integrated monitoring of structural deformation.The evaluation and testing of the monitoring accuracy of the GNSS/IMU tightly coupled algorithm in complex environments using a vibration table in [26] also demonstrated that the GNSS/IMU combination can improve accuracy and robustness compared to GNSS alone.Aiming to address the problem of using fixed Earth gravity values and ignoring Earth's rotation in traditional pre-integration algorithms, an IMU pre-integration algorithm considering Earth's rotation and gravity changes was proposed, and it was shown to improve positioning accuracy by 32.41% and yaw accuracy by 4.23% [27].Moreover, a GNSS/IMU robust adaptive filtering algorithm based on innovation was proposed, and it significantly improved positioning and attitude accuracy compared to extended Kalman filtering [28].
Currently, cableway deformation monitoring is mainly achieved through the establishment of an automated and information-based integrated platform for remote monitoring [29,30].In [31], a monitoring system for passenger cableways was preliminarily designed that can monitor and diagnose cableway information.Afterward, a set of software for parallel cableway control and monitoring was developed and successfully applied [32].With the further development of information technology, automatic, informationized, and networked monitoring has become the mainstream.Wireless networks and cloud services were used to improve the informatization level and monitoring efficiency of cableway deformation monitoring [33].Meanwhile, cableway deformation monitoring based on multi-sensor fusion has generated wide interest, and there have been many studies on the monitoring of wire ropes in cableways.It was found that regular condition monitoring of steel wire ropes based on non-destructive investigation and evaluation techniques can improve safety and optimal rope life [34].Lagerev et al. analyzed the wire rope of a cableway and observed the deformation of the wire rope through computer simulation of the overall time-varying reliability index [35].Sukhorukov et al. used magnetic flux leakage testing (MFL) in automatic wire rope monitoring systems, which can perform automatic monitoring and save time and effort [36].In addition, by combining GPR and Differential GPS (DGPS), high-frequency and high-precision data for cableway deformation monitoring can be obtained to establish a corresponding cableway deformation monitoring system [37].In terms of maintenance, fuzzy inference systems can be used to simulate the failure modes of large-scale cableway transportation systems, thereby achieving maintenance automation [38].In addition, there are relatively few studies conducted on deformation monitoring of other parts of the cableway.In [39], the safety detection device for the cableway bracket was thoroughly studied, and the fault monitoring process was successfully analyzed.Zhang et al. used acoustic emission technology to monitor the condition and diagnose faults of the bearings in the cableway [40].
Overall, current cableway deformation monitoring mainly focuses on the development of relevant automated monitoring platforms or research on the deformation monitoring of related parts such as cable cars and wire ropes.However, the cableway bracket used to support the normal operation of the steel wire rope has not been given enough attention, as potential faults and abnormal states of the cableway will have a significant impact on the three-dimensional deformation and attitude changes in the cableway bracket.Therefore, the goal of this paper is to propose a GNSS/IMU tightly coupled algorithm for cableway bracket deformation monitoring and accident avoidance.To overcome the problem of Kalman filter divergence and drifting error caused by static IMU measurement, ZUPT and the constraint of the Earth's rotation angular velocity are applied to the measurement model and state model of the Kalman filter for the GNSS/IMU tightly coupled algorithm.
The remainder of this paper is organized as follows: Section 2 lists the formulas of the dynamic and observation equations of the GNSS/IMU tightly coupled algorithm and extended Kalman filter (EKF).In Section 3, the ZUPT and the constraint of the Earth's rotation angular velocity are derived and described.The set of test results and data analysis are given in Section 4. Finally, the paper is concluded in Section 5.

Extended Kalman Filter for GNSS/IMU Tightly Coupled Algorithm
The EKF is widely considered as the best linear filter when the system state and measurement noise are both zero-mean Gaussian white noise.In this section, the dynamic and observation equations of the GNSS/IMU tightly coupled algorithm and the EKF are provided.

The Dynamic and Observation Equations of GNSS/INS Tightly Coupled Algorithm
The GNSS/INS tight combination is the process of inputting GNSS observations and predicted INS observations together into a Kalman filter, and obtaining the position and navigation results of the carrier through filtering fusion.Compared with loose combination, tight combination has better positioning accuracy and robustness, especially when the number of satellites is small or loss of lock occurs, it can significantly improve the availability of GNSS under harsh observation conditions.In addition, owing to the fact that tight combinations only need one filter for optimal estimation, they have higher computational efficiency in practical work.In this paper, the state model of the EKF based on the IMU error model for the GNSS/IMU tightly coupled system adopts an 18-dimensional state: where X is the state error vector of GNSS/INS integrated navigation; δψ e , δv e , and δr e represent the three-dimensional attitude error, velocity error, and position error, respectively; b b a and b b g are acceleration deviation and gyro drift in the carrier coordinate frame, respectively; dt r is the receiver clock error; d .t r is the receiver clock error drift; and ∇∆N is the double difference ambiguity.
The dynamic model of the GNSS/INS tightly coupled navigation system can be expressed as: where X k is the system state vector; Φ k/k−1 is the state transition matrix; ω ie is the Earth's rotation angular velocity; L and H are the geographic latitude and altitude, respectively; (V e , V n , V u ) are the velocities of the carrier in the navigation coordinate system; R M and R N are the curvature radii of the prime and meridian circles at the location of the carrier; ( f e , f n , f u ) are the specific force outputs of the accelerometer; C n b is the coordinate conversion matrix from the navigation coordinate system to the carrier coordinate system; F 44 and F 55 are the zero-bias error state transition matrices of accelerometers and gyroscopes based on first-order Gaussian Markov processes; and t α is the clock drift correlation time coefficients based on first-order Gaussian Markov processes.w k−1 is the process noise vector, and its mean and covariance are as follows: where E[•] is the mean and Q k is the process noise covariance matrix.
Generally, the measurement model established based on the difference between the pseudorange, carrier phase, and Doppler measurement values of the GNSS and INS are as follows: where Z k represents the measurement vector; H k is the measurement matrix; D and E represent the design matrices of position error vectors and velocity error vectors, respectively; C represents the error conversion matrix between the spatial Cartesian coordinate system and the geodetic coordinate system; v k is the measurement noise vector; v k and w k−1 are uncorrelated white noise sequences, and their mean and covariances are as follows: where E[•] is the mean, and R k is the measurement noise covariance matrix.In addition, owing to the short baseline differential positioning model used in this paper, errors such as ionospheric delay and tropospheric delay can be reduced or eliminated.

Extended Kalman Filter
Generally, a Kalman filter includes two parts, a state update and a measurement update, and a linear mathematical model is required.Unfortunately, in practical applications, most systems such as the measurement model in GNSS/INS tightly integrated navigation are nonlinear.In this case, the nonlinear equation needs to be linearized, and the EKF is a widely used method that involves Taylor expansion of nonlinear functions.
State update of EKF: where X k/k−1 and P k/k−1 are the predicted state vector and its covariance; X k−1 and P k−1 represent the state update vector and its covariance.Measurement update steps of EKF: (a) The gain matrix K k : (b) The innovation vector (c) The state update vector X k : (d) The covariance of state update vector P k : where K k describes the weights between new measurements and predicted information in the system dynamics model.

ZUPT-Based GNSS/IMU Tightly Coupled Algorithm with the Constraint of the Earth's Rotation Angular Velocity
In this section, in order to obtain reliable position and attitude results using the GNSS/IMU tightly coupled algorithm, the zero-velocity update was added to the measurement model, and the constraint of the Earth's angular velocity was added to the state model.

ZUPT
When the carrier is in a stationary state, the velocity of the carrier is zero.Inertial navigation errors can be corrected through IMU zero-velocity detection technology [41].ZUPT utilizes the calculated velocity of the IMU in the carrier as the observation of the system velocity error, corrects other error quantities, and improves the integrated navigation results in a stationary state without the need for external sensors.Therefore, it is an effective, inexpensive, and easy-to-implement technology.Therefore, when the carrier is judged to be in a zero-velocity state, the zero-velocity constraint equation will be added to the observation Equation ( 5) to form the zero-velocity observation equation: The updated state estimation value is used to provide feedback to and correct the navigation parameter error inside the IMU, completing the static zero-velocity correction.The precision of the zero-velocity constraint equation is extremely high.By giving a small prior variance, the Kalman filtering parameter of the zero-velocity additional equation is used to solve for strong constraints, achieving the goal of correcting the accumulation of static state errors.
In the static deformation monitoring environment of the cableway, the constraint method of the zero-velocity update (ZUPT) is particularly useful for the static GNSS/IMU tightly coupled deformation monitoring data.It was found that by constraining the threeaxis velocity of inertial navigation, the divergence of the Kalman filter can be effectively suppressed and the corresponding pose data can be successfully calculated.However, this technology cannot suppress the drift of the combined positioning and attitude error.

Constraint of the Earth's Rotation Angular Velocity
Although ZUPT has a certain inhibitory effect on the Kalman filter divergence when the carrier is stationary, it is just the constraint on the measurement model.There may still be position divergence and yaw drift during the calculation process, because the rotation of the navigation system caused by the curvature of the Earth's surface in the state of the velocity update in the inertial navigation mechanical arrangement: ) where ω n ie is the angular velocity; the subscript ie represents the direction of rotation (the rotation of the Earth-fixed coordinate system relative to the inertial coordinate system); the superscript n represents the corresponding coordinate system (which, here, is the navigation coordinate system); C n e is the coordinate conversion matrix from the navigation coordinate system to the Earth-fixed coordinate system; and ω ie is the Earth's rotation angular velocity.
To address the problem mentioned above, based on the carrier being stationary, a further constraint on the Earth's angular velocity is applied to the state model under the ZUPT-based tightly coupled solution, that is, to constrain the yaw direction component of the Earth's rotation angular velocity by setting: As a result, the state transition matrix (3) with the constraint of the Earth's rotation angular velocity (CERAV) is: According to the descriptions in Sections 3.1 and 3.2, the flowchart for practical implementation of the ZUPT-based GNSS/IMU tightly coupled algorithm with CERAV is shown in Figure 1.

Experiments and Results Analysis
In this section, the proposed method, called the ZUPT-based GNSS/IMU tightly coupled algorithm with CERAV (ZUPT-TC-CERAV), is compared with the ZUPT-based GNSS/IMU tightly coupled algorithm (ZUPT-TC).Due to the divergence of the Kalman filter, the GNSS/IMU tightly coupled algorithm without constraint cannot provide position and attitude solutions, so it is not considered.In addition, the GNSS alone (which can only provide the position not affected by IMU divergence, but cannot provide attitude) is used to demonstrate the effectiveness of the proposed method.
Experimental data were collected on the crossbar of the cableway in the Yimeng Mountain Tourism area on 23 September 2022, over a duration of approximately 15 min.The environment of the cableway bracket is shown in Figure 2. The tourist sightseeing cable car station is located in the southern area of the scenic area.The pick-up station is located at the top of the main peak tower in the southern area, and the drop-off station is located at Waterfall Square.The one-way passenger capacity of the sightseeing cable car is 1668 people/h, and each carriage can accommodate 160 people.The inclined length of the route is 1835.5 m.

Experiments and Results Analysis
In this section, the proposed method, called the ZUPT-based GNSS/IMU tightly coupled algorithm with CERAV (ZUPT-TC-CERAV), is compared with the ZUPT-based GNSS/IMU tightly coupled algorithm (ZUPT-TC).Due to the divergence of the Kalman filter, the GNSS/IMU tightly coupled algorithm without constraint cannot provide position and attitude solutions, so it is not considered.In addition, the GNSS alone (which can only provide the position not affected by IMU divergence, but cannot provide attitude) is used to demonstrate the effectiveness of the proposed method.
Experimental data were collected on the crossbar of the cableway in the Yimeng Mountain Tourism area on 23 September 2022, over a duration of approximately 15 min.The environment of the cableway bracket is shown in Figure 2. The tourist sightseeing cable car station is located in the southern area of the scenic area.The pick-up station is located at the top of the main peak tower in the southern area, and the drop-off station is located at Waterfall Square.The one-way passenger capacity of the sightseeing cable car is 1668 people/h, and each carriage can accommodate 160 people.The inclined length of the route is 1835.5 m.The environment of the cableway bracket is shown in Figure 2. The tourist sightseeing cable car station is located in the southern area of the scenic area.The pick-up station is located at the top of the main peak tower in the southern area, and the drop-off station is located at Waterfall Square.The one-way passenger capacity of the sightseeing cable car is 1668 people/h, and each carriage can accommodate 160 people.The inclined length of the route is 1835.5 m.The experimental platform consists of two Hi-Target receivers and a SPAN-IGM-A1 IMU.The update rate of the Hi-Target receiver is 1 Hz, and that of the SPAN-IGM-A1 IMU is 200 Hz.The technical parameters for the SPAN-IGM-A1 IMU are listed in Table 1, and the experimental platform is shown in Figure 3.The experimental platform consists of two Hi-Target receivers and a SPAN-IGM-A1 IMU.The update rate of the Hi-Target receiver is 1 Hz, and that of the SPAN-IGM-A1 IMU is 200 Hz.The technical parameters for the SPAN-IGM-A1 IMU are listed in Table 1, and the experimental platform is shown in Figure 3.  Figure 4 is a schematic diagram of the position relationship between the navigation coordinate and the cableway coordinate.In this paper, the baseline has been transformed from the navigation coordinate (E-N-U) to the cableway coordinate (X-Y-Z) in order to facilitate the exploration of position and attitude.The number of satellites and dilution of precision (DOP) are shown in Figures 5 and 6.It is obvious that the satellite observation environment is good.from the navigation coordinate (E-N-U) to the cableway coordinate (X-Y-Z) in order to facilitate the exploration of position and attitude.The number of satellites and dilution of precision (DOP) are shown in Figures 5 and 6.It is obvious that the satellite observation environment is good.Figures 7 and 8, respectively, contain the position of GNSS, ZUPT-TC, and ZUPT-TC-CERAV in navigation coordinates and cableway coordinates, and it can be seen that the position solution can be obtained by applying ZUPT, but the problem of filtering divergence still exists.Especially after turning the coordinates, the position of ZUPT-TC diverged in all three directions.Meanwhile, the position of ZUPT-TC-CERAV is relatively stable owing to the constraint of the state model, which also leads to the difference in the initial position of ZUPT-TC and ZUPT-TC-CERAV.In addition, the phenomenon that the  Figures 7 and 8, respectively, contain the position of GNSS, ZUPT-TC, and ZUPT-TC-CERAV in navigation coordinates and cableway coordinates, and it can be seen that the position solution can be obtained by applying ZUPT, but the problem of filtering divergence still exists.Especially after turning the coordinates, the position of ZUPT-TC diverged in all three directions.Meanwhile, the position of ZUPT-TC-CERAV is relatively stable owing to the constraint of the state model, which also leads to the difference in the Figures 7 and 8, respectively, contain the position of GNSS, ZUPT-TC, and ZUPT-TC-CERAV in navigation coordinates and cableway coordinates, and it can be seen that the position solution can be obtained by applying ZUPT, but the problem of filtering divergence still exists.Especially after turning the coordinates, the position of ZUPT-TC diverged in all three directions.Meanwhile, the position of ZUPT-TC-CERAV is relatively stable owing to the constraint of the state model, which also leads to the difference in the initial position of ZUPT-TC and ZUPT-TC-CERAV.In addition, the phenomenon that the position of ZUPT-TC-CERAV basically matches the position of the GNSS, which is not affected by IMU divergence, verifies the effectiveness of ZUPT-TC-CERAV again.The attitude of ZUPT-TC and ZUPT-TC-CERAV are provided in Figure 9. Obviously, compared to ZUPT-TC-CERAV, the trends of roll and pitch obtained by ZUPT-TC are relatively stable, while that of yaw shows divergence.At the same time, the trends of roll, pitch, and yaw obtained by ZUPT-TC-CERAV are all stable.The reason for this could be that the drifting error of yaw in the stationary state is significantly suppressed by adding   The attitude of ZUPT-TC and ZUPT-TC-CERAV are provided in Figure 9. Obviously, compared to ZUPT-TC-CERAV, the trends of roll and pitch obtained by ZUPT-TC are relatively stable, while that of yaw shows divergence.At the same time, the trends of roll, pitch, and yaw obtained by ZUPT-TC-CERAV are all stable.The reason for this could be that the drifting error of yaw in the stationary state is significantly suppressed by adding The attitude of ZUPT-TC and ZUPT-TC-CERAV are provided in Figure 9. Obviously, compared to ZUPT-TC-CERAV, the trends of roll and pitch obtained by ZUPT-TC are relatively stable, while that of yaw shows divergence.At the same time, the trends of roll, pitch, and yaw obtained by ZUPT-TC-CERAV are all stable.The reason for this could be that the drifting error of yaw in the stationary state is significantly suppressed by adding the constraint of Earth's rotation angular velocity.To further illustrate the difference between ZUPT-TC and ZUPT-TC-CERAV, the position error and attitude error are shown in Figures 10 and 11.It is indicated that the position error of ZUPT-TC shows a divergent trend in the X, Y, and Z directions, while the errors of the ZUPT-TC-CERAV solution are all relatively stable, and the difference in the position error between GNSS and ZUPT-TC-CERAV is small.Similarly, the same divergent trend caused by the stationary IMU is shown in the yaw error of ZUPT-TC, and ZUPT-TC-CERAV can effectively solve this problem.To further illustrate the difference between ZUPT-TC and ZUPT-TC-CERAV, the position error and attitude error are shown in Figures 10 and 11.It is indicated that the position error of ZUPT-TC shows a divergent trend in the X, Y, and Z directions, while the errors of the ZUPT-TC-CERAV solution are all relatively stable, and the difference in the position error between GNSS and ZUPT-TC-CERAV is small.Similarly, the same divergent trend caused by the stationary IMU is shown in the yaw error of ZUPT-TC, and ZUPT-TC-CERAV can effectively solve this problem.From the description in Table 2 and the analysis of the experimental results, we can observe that ZUPT-TC-CERAV is able to overcome the problem of Kalman filter divergence and drifting error caused by the static IMU measurement, in contrast to ZUPT-TC, and its reliable performance of positioning and attitude determination is demonstrated by the minor difference from single GNSS.

Discussion
The cableway bracket used to support the normal operation of the steel wire rope plays an important role in the detection of potential faults and abnormal states, as the cableway will have a significant impact on the three-dimensional deformation and attitude changes in the cableway bracket when an accident occurs.Although a single GNSS algorithm can provide a certain degree of accuracy in position information for deformation monitoring, it cannot provide attitude changes.Therefore, the GNSS/IMU tightly coupled algorithm is a must to obtain high-precision position and attitude data for cableway bracket deformation monitoring.However, when the IMU is stationary for a long time, it can cause significant drift errors.As a result, the GNSS/IMU tightly coupled algorithm without constraint cannot provide position and attitude solutions.
To overcome this problem, the ZUPT and the constraint of the Earth's rotation angular velocity are applied to the measurement model and the state model of the Kalman filter.By adding the zero-velocity constraint equation to the observation equation, ZUPT-TC can obtain position and attitude solutions, but there will still be drift in the position and yaw.Furthermore, the constraint of the Earth's rotation angular velocity is added to the state model, and then the position and attitude that match the actual situation are obtained.However, the absolute accuracy of the proposed method is still unavailable (because the reliable solution cannot be obtained), so we can only obtain its relative accuracy by comparing its solution with that of single GNSS.

3. 3 .
ZUPT-Based GNSS/IMU Tightly Coupled Algorithm with CERAV According to the descriptions in Sections 3.1 and 3.2, the flowchart for practical implementation of the ZUPT-based GNSS/IMU tightly coupled algorithm with CERAV is shown in Figure 1.Sensors 2023, 23, 9862 10 of 17

Figure 2 .
Figure 2. The environment of cableway bracket.Figure 2. The environment of cableway bracket.

Figure 2 .
Figure 2. The environment of cableway bracket.Figure 2. The environment of cableway bracket.

Figure 3 .
Figure 3. Experimental platform used in the test.

Figure 3 .
Figure 3. Experimental platform used in the test.

Figure 4
Figure4is a schematic diagram of the position relationship between the navigation coordinate and the cableway coordinate.In this paper, the baseline has been transformed from the navigation coordinate (E-N-U) to the cableway coordinate (X-Y-Z) in order to facilitate the exploration of position and attitude.The number of satellites and dilution of precision (DOP) are shown in Figures5 and 6.It is obvious that the satellite observation environment is good.

Figure 4 .
Figure 4. Position relationship between navigation coordinate and cableway coordinate.Figure 4. Position relationship between navigation coordinate and cableway coordinate.

Table 2 .
The RMSs of position and attitude errors of ZUPT-TC and ZUPT-TC-CERAV.