Next Article in Journal
Adaptive Subspace Signal Detection in Structured Interference Plus Compound Gaussian Sea Clutter
Next Article in Special Issue
Validating Precise Orbit Determination from Satellite-Borne GPS Data of Haiyang-2D
Previous Article in Journal
Combining Different Transformations of Ground Hyperspectral Data with Unmanned Aerial Vehicle (UAV) Images for Anthocyanin Estimation in Tree Peony Leaves
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Technical Note

Adaptive Kalman Filter for Real-Time Precise Orbit Determination of Low Earth Orbit Satellites Based on Pseudorange and Epoch-Differenced Carrier-Phase Measurements

1
Institute of Space Sciences, Shandong University, Weihai 264209, China
2
College of Geology Engineering and Geomatics, Chang’an University, Xi’an 710054, China
3
Beijing Urban Construction Exploration & Surveying Design and Research Institute, Beijing 100101, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(9), 2273; https://doi.org/10.3390/rs14092273
Submission received: 14 March 2022 / Revised: 3 May 2022 / Accepted: 6 May 2022 / Published: 8 May 2022
(This article belongs to the Special Issue Precision Orbit Determination of Satellites)

Abstract

:
Real-time precise orbit determination (POD) of low earth orbiters (LEOs) is crucial for orbit maintenance as well as autonomous operation for space missions. The Global Positioning System (GPS) has become the dominant technique for real-time precise orbit determination (POD) of LEOs. However, the observation conditions of near-earth space are more critical than those on the ground. Real-time POD accuracy can be seriously affected when the observation environment suffers from strong space events, i.e., a heavy solar storm. In this study, we proposed a reliable adaptive Kalman filter based on pseudorange and epoch-differenced carrier-phase measurements. This approach uses the epoch-differenced carrier phase to eliminate the ambiguities and thus reduces the significant number of unknown parameters. Real calculations demonstrate that four to five observed GPS satellites is sufficient to solve reliable position parameters. Furthermore, with accurate pseudorange and epoch-differenced carrier-phase-based reference orbits, orbital dynamic disturbance can be detected precisely and reliably with an adaptive Kalman filter. Analyses of Swarm-A POD show that sub-meter level real-time orbit solutions can be obtained when the observation conditions are good. For poor observation conditions such as the GRACE-A satellite on 8 September 2017, when fewer than five GPS satellites were observed for 14% of the observation time, 1–2 m orbital accuracy can still be achieved with the proposed approach.

1. Introduction

In recent years, low earth orbit (LEO) satellites have been widely used in scientific research as well as military, civil and other fields, such as gravimetry, altimetry, meteorology and other earth observation missions [1,2,3,4,5,6]. The operation of scientific missions of these LEO satellites requires precise post- or real-time orbits. Currently, the precise orbit determination (POD) of LEO satellites mainly depends on the GPS technique. GPS shows incomparable advantages in the POD of LEO satellites since it is an all-weather, high accuracy, low cost and continuous observation system. GPS together with other new global navigation satellite system, i.e., BeiDou navigation satellite system (BDS) and Galileo, will further improve the performance of LEO POD [7]. Usually, centimeter level precise orbits can be obtained using the dynamic or reduced-dynamic approach in post-processing mode, and relevant techniques have been extensively studied [8,9,10,11,12,13].
However, continued advancements in remote sensing technology along with a trend towards highly autonomous spacecraft provide a strong motivation for accurate real-time navigation of LEO satellites [14]. Autonomous navigation refers to the carrier being able to receive and process the spaceborne GPS data automatically. Position, velocity and other state information can therefore be generated and provided to the other instruments mounted on this satellite or spacecraft for autonomous operation of the whole system [15]. For example, many LEO satellites will be launched into space for navigation augmentation of Global Navigation Satellite System (GNSS) satellites in the future as planned. Such a huge constellation of LEO satellites can also serve as navigation satellites and transmit navigation signals to the ground. The broadcast ephemeris needs to be generated and broadcasted to the global users for positioning services. It will be a massive task for ground processing centers to process such a large quantity of spaceborne and ground tracking data. It is also not easy for the ground station to upload the ephemeris information to the LEO satellites, which is usually performed for the GNSS satellites since the LEOs move relatively fast (~7.9 km/s) and it is not easy to keep the data link uninterrupted. Therefore, one more practical and easy implementation method is for LEO satellites to generate real-time orbits automatically and independently.
Autonomous orbit determination is also the main approach for attitude and orbit control systems of LEO satellites. Compared to GNSS, LEO satellites have a much lower orbital altitude and suffer more significant atmospheric drag forces. This may cause a great burden to ground facilities if the orbit maintenance work of a large constellation of LEO satellites is performed on the ground. Therefore, autonomous orbit determination is necessary and urgent. Orbital anomalies can also be detected over time.
Finally, applications of real-time space weather monitoring and Earth observation would also benefit from accurate real-time LEO orbit information, including the onboard geocoding of high-resolution imagery, the open loop operation of altimeters and atmospheric sounding [16,17,18].
It is not easy to retrieve stable and high-accuracy real-time orbit information of LEO satellites and spacecrafts at all times. The performance of onboard real-time POD needs to consider the balance among computational efficiency, in-orbit processor resources and accuracy. The accuracy calculated from pseudorange observations is at the meter level. Sometimes it is necessary to introduce carrier-phase observations to improve the accuracy. Since gross errors usually appear more frequently in spaceborne observations than those from the ground, and the observation is often discontinuous, a large number of outliers may occur in kinematic orbit determination solutions. From this point of view, the dynamic model as well as the Kalman filter technique is generally applied for real-time POD in order to improve the reliability and stability of orbit solutions. Decimeter level orbital accuracy can be achieved with pseudorange and carrier-phase measurements when the observation conditions are good [14,19].
It is obvious that including the spaceborne carrier-phase measurements leads to the problem of resolving a large number of ambiguities since the continuous observation arc is short from the view of a LEO satellite to the GNSS. The limited capability of the onboard processor restricts the resolution of many ambiguities. Moreover, in the case of insufficient observation data or frequent loss of satellite tracking, the estimation of a large number of ambiguities together with the dynamic model parameters would also lead to the singularity of solutions, and sometimes the results are even inferior to those calculated with pseudorange-only observations. Indeed, we need to consider the balance between the capability of the onboard processor and the quality of observations to obtain optimal solutions.
Therefore, we propose to use epoch-differenced carrier-phase measurements to remove ambiguities, and pseudorange-measurement-based solutions are taken as the a priori orbits. When there is an insufficient number of observed satellites, observations from a small number of satellites are sufficient for epoch-differencing to derive the position increments, and the position at the current epoch is obtained by orbit stacking. This approach processed in real-time is different from that in post-processing mode when all observations are stacked to resolve the position parameters entirely [20]. The position increments are estimated epoch-by-epoch in real-time, and together with the known positions from the last epoch, the position at the current epoch is finally obtained in a combined adjustment. This approach shows advantages when the number of observations is insufficient to resolve all ambiguities reliably.
In addition, the LEO satellites are faced with more complex and harsh environments in the near-earth space than the ground, and the most significant effect is atmospheric drag forces. Solar storms can also degrade the orbit accuracy of LEO satellites significantly, which is mainly caused by the increased density of atmosphere [21]. The conventional dynamic model at a usual sampling interval is not sensitive to the dynamic disturbance, and the influence on real-time POD accuracy is at the decimeter to meter level. Such a disturbance is not easily detected using purely pseudorange observations, while the inclusion of epoch-differenced carrier-phase measurements can improve performance. In this study, we propose to apply the adaptive Kalman filter to detect the state disturbance by using an adaptive factor with the pseudorange and epoch-differenced carrier-phase measurements.
Following this introduction, the algorithms and models of adaptive Kalman filter based on pseudorange and epoch-differenced carrier-phase measurements are described in detail in Section 2. The data and processing strategies are introduced in Section 3. Finally, the real-time POD results of various schemes used to evaluate our approaches are analyzed in Section 4. The conclusions are drawn in Section 5.

2. Methods

2.1. Preprocessing of Pseudorange and Carrier-Phase Measurements

Generally, there are more cycle slips on spaceborne observations than those on the ground since the continuous observation time from a GNSS to a LEO satellite is much shorter (15–25 min). It is also a challenge for a fast-moving LEO satellite to capture the GNSS signals continuously, as gross errors may occur more frequently than ground receivers. Therefore, accurate data preprocessing is important. In real-time POD, the on-the-fly data editing algorithm is generally used to preprocess the pseudorange and carrier-phase measurements. In the filtering process, for each GNSS satellite, the predicted LEO orbits from orbital dynamic integration and the GNSS orbits from the broadcast ephemeris are used to calculate the receiver clock errors, and the pseudorange residuals of each satellite are checked for whether they exceed a given threshold value. If so, the pseudorange observations of this satellite are eliminated. The detailed algorithms can be found in [14]. The carrier-phase measurements are dealt with in a similar way. When forming the epoch-differenced carrier-phase observation equations with the predicted LEO orbits and GNSS orbital positions, most error sources can be eliminated. The remaining unknowns are only the difference of receiver clock offsets between consecutive epochs. At a sampling interval of 10 s, the difference can be at 1.2 dm for the BlackJack Receiver [22] and 1.8 dm for the Global Navigation Satellite System Occultation Sounder (GNOS) spaceborne receiver [13]. Therefore, it can be concluded that for a commonly used spaceborne GPS receiver, the effects of receiver clock offsets are generally at the decimeter level. If the residual of the epoch-differenced carrier-phase measurements of a certain satellite at a certain frequency is at the meter level, there must be cycle slips larger than 1 cycle.
The determination of the position of cycle slip starts from the first epoch. From the sequential process, the cycle slip occurring on the current epoch can be detected. Large cycle slips should be repaired to under 1 cycle. Then, the effect of cycle slips below 1 cycle on positioning is at the centimeter level and can be neglected in real-time POD. In addition, errors from cycle slips do not pass on to the next epoch and therefore have no effect on the following solutions.

2.2. Adaptive Kalman Filter Based on Pseudorange and Epoch-Differenced Carrier-Phase Measurements

2.2.1. Observation Model

The ionosphere-free (IF) observation equation is expressed as follows:
d P C r j = x r x j ρ r j d x + y r y j ρ r j d y + z r z j ρ r j d z + c δ d t r
d L C r j = x r x j ρ r j d x + y r y j ρ r j d y + z r z j ρ r j d z + B r j + c δ d t r
where d P C r j and d L C r j are IF observed minus calculated (o-c) pseudorange and carrier-phase measurements, respectively, for satellite j and receiver r . x r , y r , z r is the a priori receiver position at three directions, d x , d y , d z is the corresponding position corrections. ρ is the geometry distance. x j , y j , z j is satellite fixed position at three directions. c is the speed of light. δ d t r is the receiver clock offset correction in the estimation. B is the ambiguity parameter. The unknowns are 3 positional correction parameters, receiver clock offsets and ambiguities. The tropospheric delay is not considered since the LEO orbit is usually above the troposphere. In order to eliminate the ambiguity parameter, the following normal equation at epoch k 1 is derived for epoch-differencing:
δ L C r , k 1 j = L C r , k 1 j ρ r , k 1 j c d t r , k 1 + c d t r , k 1 j
where ρ r , k 1 j is calculated based on positions at epoch k 1 , and the ambiguity parameter B is eliminated in the differencing process.
Next is the determination of the weight for pseudorange and carrier-phase measurements. For the two types of observations, we have the following covariance matrix for observational noise:
D w = D P D P Δ Φ D Δ Φ P D Δ Φ
where D P and D Δ Φ are the variance for pseudorange and epoch-differenced carrier-phase measurements, respectively. For simplicity, the correlation between the pseudorange and the differenced carrier phase is neglected, and therefore D P Δ Φ and D Δ Φ P are both set to 0. The initial standard deviation for the pseudorange and epoch-differenced carrier phase is set to 1 and 0.01, respectively. Therefore, the weight of the carrier-phase measurements is about 104 times the pseudorange measurements. The pseudorange-based solution provides the absolute reference orbit, and the carrier-phase-based solution dominates the final orbit estimates.

2.2.2. Dynamic Model

The Kalman filter is used to solve the state parameters in real-time considering the dynamic models. Based on pseudorange and epoch-differenced carrier-phase measurements, the ambiguity parameter is removed. The unknown parameters can be expressed as:
x = r , r ˙ , C R , C D , a , c d t r
which comprises 3-dimensional positions and velocity vectors r and r ˙ in an Earth-fixed reference frame; one solar radiation pressure coefficient C R ; one atmospheric drag coefficient C D ; three empirical acceleration vectors a = a R , a T , a N in radial, along and cross directions and the receiver clock offset parameter c d t r .
First the Kalman filter is initialized. We set X ^ 0 = X 0 r e f , P ^ 0 = P 0 r e f , where X ^ 0 is the initial state estimate, and here it is obtained from the kinematic solution X 0 r e f ; P ^ 0 is the initial state covariance, and the setting value P 0 r e f can be seen in Table 3 in Section 3.2.
The predicted state including position and velocity can be calculated based on dynamic integration using a 4th-order Runge–Kutta scheme with Richardson extrapolation, and the prediction of empirical accelerations can be derived as the following:
a ¯ k = e t k t k 1 / τ a ^ k 1
C ¯ R , C ¯ D , c δ t ¯ k = C ^ R , C ^ D , c δ t ^ k 1
where τ is the correlation time of empirical accelerations. a ^ k 1 is the empirical acceleration at epoch time t k 1 , and a ¯ k is the predicted value at epoch time t k . The predicted values of coefficient C R and C D and the predicted value of the receiver clock offset c δ t r do not change in the filter propagation process.
Firstly, we need to calculate the state transition matrix, Φ k , k 1 = X k r e f / X k 1 r e f , where X k r e f means the reference orbit, which is calculated by integration of the satellite dynamics:
Φ k , k 1 = ϕ x ϕ C R , C D ϕ a 0 0 I 1 0 0 0 0 M 0 0 0 0 I 2 12 + n × 12 + n
where ϕ x = r / r r / r ˙ r ˙ / r r ˙ / r ˙ 6 × 6 , ϕ C R , C D = r / C R r / C D r ˙ / C R r ˙ / C D 6 × 2 , ϕ a = r / a R , T , N r ˙ / a R , T , N 6 × 3 , M is a 3 × 3 diagonal matrix, and the ith element m i = e t k t k 1 / τ , I 1 is 2 × 2 unit matrix, I 2 is the n + 1 × n + 1 unit matrix.
Then, the predicted state X ¯ k and state transition matrix Φ k at epoch k can be calculated based on the state estimate X ^ k 1 and covariance matrix P ^ k 1 at epoch k 1 , and we have
X ¯ k = X k r e f + Φ k , k 1 X ^ k 1 X k 1 r e f
where X k r e f is the reference orbit at epoch k calculated based on observation information. Thus, the covariance matrix P ¯ k for the predicted state X ¯ k is expressed as:
P ¯ k = Φ k , k 1 P ^ k 1 Φ k , k 1 T + Q k
with P ^ k - 1 as the covariance for the estimates at time t k 1 , Q k as the covariance increases due to the accumulated effect of process noise:
Q k = 0 0 0 0 0 Q a 0 0 0 0 Q δ t 0 0 0 0 0 12 + n × 12 + n
where Q a is covariance for the three empirical accelerations and is a 3 × 3 diagonal matrix, the ith element is expressed as q i = σ i 2 1 m i 2 and σ i 2 i = 1 , 3 is the corresponding steady-state variance. For the clock offset, we have Q δ t = σ δ t 2 τ δ t t k t k 1 , where σ δ t 2 and τ δ t are the steady-state variance and the correlation time of receiver clock offset, respectively.
The measurement is therefore updated, and the gain matrix K ¯ k and state estimates X ^ k together with the covariance matrix P ^ k can be calculated:
K ¯ k = P ¯ k H k T H k P ¯ k H k T + R k 1
X ^ k = X ¯ k + K ¯ k L k
P ^ k = I K ¯ k H k P ¯ k I K ¯ k H k T + K k R k K k T
where R k is the covariance matrix for observation noise. Equation (14) can guarantee the non-negative property of the covariance matrix. H is the design matrix, and L k is the observed minus calculated pseudorange and epoch-differenced carrier-phase measurements. If n satellites are observed at epoch k , we have the following expression for H and L k :
H = e 1 x k e 1 y k e 1 z k 1 e 1 x k e 1 y k e 1 z k 1 e 2 x k e 2 y k e 2 z k 1 e 2 x k e 2 y k e 2 z k 1 e n x k e n y k e n z k 1 e n x k e n y k e n z k 1 , L k = d P C r , k 1 d L C r , k 1 δ L C r , k 1 1 d P C r , k 2 d L C r , k 2 δ L C r , k 1 2 d P C r , k n d L C r , k n δ L C r , k 1 n
Usually, the extended Kalman filter (EKF) is used for kinematic positioning of a fast-moving carrier. This is to solve the problem of easy divergence of the conventional linear Kalman filter for LEO POD [14,23]. However, the EKF also leads to the problem of re-initialization of the state transition matrix Φ k , k 1 , which consumes the spaceborne computing resources. This may cause a problem when resources are limited. This is also not convenient for programming. Generally, the reference orbit can be obtained by integration of the satellite dynamics and can be denoted as X k r e f . Here, for sequential dynamic model disturbance detection, we need to construct a robust reference orbit for the subsequent Kalman filter, which is calculated using only the observation information in a kinematic navigation solution, and we have:
X ˜ k = X k r e f + x k
where x k = H k T P k H k 1 H k T P k L k , P k is the weight for observations, P k = R k 1 . In contrast to the previous study where the kinematic solution is only from the pseudorange observations, here the inclusion of epoch-differenced carrier-phase measurements can obtain an improved reference orbit; therefore, the filter is not easy to diverge.
It can be seen that the LEO position can be solved with only 4 sequential observed GPS satellites with epoch-differenced observations as the unknowns are three positional parameters and one receiver clock offset. This is also demonstrated in [24,25,26]. In this approach, the pseudorange measurements can provide the absolute reference for orbit solutions. The position increments calculated with epoch-differenced carrier-phase measurements can be added to the filter estimated orbits, and therefore high accuracy orbits are able to be passed on to the next epoch. This technique manifests its advantage when there is an insufficient number of observed satellites. The position increments dominate the process of high accuracy orbit inheritance since the orbit solutions during the period of poor observation conditions can still be estimated with position increments at these epochs and filter estimated high accuracy orbit at a historical epoch.
In addition, the observation conditions in near-earth space are more challenging than those on the ground. LEO satellites suffer more atmosphere drag forces from solar storms impact. Under such circumstances, the dynamic disturbance can cause significant impact on LEO POD. Therefore, we propose to refine the predicted orbit based on an adaptive Kalman filter. This is achieved by adjusting the predicted state covariance:
P ¯ k = P ¯ k / α k = Φ k P ^ k 1 Φ k T + Q k / α k
where P ¯ k is the equivalent predicted covariance. The predicted state covariance P ¯ k need to be adjusted when there is a discrepancy between the predicted state X ¯ k and the estimated state X ^ k . α k is the adaptive factor. In order to successfully implement this approach, one needs to consider the observation conditions. When the number of observed and valid satellites is larger than 4, the following equation is used to construct the dynamic model error discriminant statistics:
α k = 1 Δ X ¯ k c c / Δ X ¯ k Δ X ¯ k > c
where Δ X ¯ k = X ˜ k X ¯ k t r ( Σ X ¯ k ) , X ˜ k X ¯ k = Δ x 2 + Δ y 2 + Δ z 2 1 2 , X ˜ k is the observation information-based kinematic orbit and is shown in Equation (16). Δ x , Δ y , Δ z means the state discrepancies in x , y , and z directions. t r refers to solving the trace.
The determination of constant c is related to the precision of the dynamic model and the observation quality. In this study, c is set to 2.5 based on numerous empirical experiments. When the number of satellites is lower than 4, we cannot obtain reliable solutions from observation-only information; therefore, Equations (17) and (18) are not applicable, and α k is set to 1. In this case, the inclusion of epoch-differenced carrier-phase measurements plays a minor role when the observation conditions are not good.
The advantage of applying the pseudorange and epoch-differenced carrier-phase measurements to adjust the state discrepancies is that a more accurate reference orbit can be obtained to distinguish the state discrepancy. The effect of this method can be reviewed intuitively through the following figure (Figure 1):
The blue line represents the reference orbit, while the red line represents the true orbit in the case when there is orbit dynamic disturbance. The blue dots indicate the reference orbit based on pseudorange-only measurements, while the red dots indicate the pseudorange and epoch-differenced carrier-phase-based orbits. It can be inferred that the orbit calculated with the pseudorange measurements lacks sufficient accuracy to detect the dynamic orbit anomaly reliably. With the inclusion of the carrier-phase measurements, the orbits tend to be more accurate and are distributed around the “true” orbits; therefore, the orbit dynamics anomaly can be reliably detected. Usually the epoch-differenced observations can induce error accumulation problems. In this case, the reference orbit is updated by using the final orbit solutions at each certain period of time, and the accumulation errors are cleared. When there is an insufficient number of observations, the reference orbits are updated at every epoch with the help of orbit dynamics.
However, the disadvantage of this approach is that we need to store the orbital solution information of current epoch. This may become critical when there are limit storage resources onboard this satellite.

3. Materials

3.1. Data

The spaceborne GPS data from Swarm-A and Gravity Recovery and Climate Experiment (GRACE)-A satellites collected on 8 September 2017 were used for the experimental analysis. The constellation configuration parameters are listed in Table 1:
A strong solar geomagnetic storm occurred on 8 September 2017, and the post-processing POD of LEO satellites were seriously affected according to [21,28] due to the large dynamic model error caused by increased atmospheric density in the thermosphere. In this study, the data on that day were chosen for real-time analysis. In order to overcome the effect of increased atmospheric density on precise POD, the conventional method is to introduce more frequent dynamic model parameters. Here, we analyzed the approach proposed in this study to assess the state disturbance of real-time POD.

3.2. Processing Strategy

The parameter settings for real-time POD are as follows. Due to the limited power of the onboard processor unit, the order of the gravity model is set to 70 × 70. There is only 2–3 cm improvement in the POD with a higher order of gravity model which would increase the amount of computation significantly. The settings of the other dynamic model parameters are listed in Table 2. As to the observation model, undifferenced dual-frequency IF LC and PC combined observations are used, and the a priori constraints for pseudorange and carrier-phase measurements are 1 m and 0.02 cycles, respectively. The cut-off elevation is 1° for spaceborne observations. The sampling interval is 10 s. The ambiguities are estimated as float values when normal carrier-phase measurements are used.
The settings of initial variance, steady-state variance and correlation time are shown in Table 3. These values were set according to numerous experimental tests for optimal solutions.

4. Results and Analysis

Models and algorithms mentioned above were used to evaluate the performance of real-time POD under poor observation conditions, i.e., solar storm weather. In order to validate the feasibility of algorithms, real-time POD of Swarm-A satellite on a normal day was performed in advance.
The spaceborne observation data collected on 24 April 2020 were used for experimental analysis. The corresponding GPS broadcast ephemeris was obtained from the IGS analysis center. In contrast to post-processing, where all observations are stacked and the ambiguities and dynamic model parameters are solved in a batch-processing mode, the observation data were processed epoch-by-epoch in real time. The ambiguity was resolved as a float constant sequentially. The single epoch ambiguity resolution accuracy is around the meter level using broadcast ephemeris.
Five kinds of schemes were designed to evaluate the performance of the approach proposed in this study:
Scheme 1: real-time POD based on pseudorange and carrier-phase measurements with Kalman filter;
Scheme 2: real-time POD based on pseudorange-only measurements with Kalman filter;
Scheme 3: real-time POD based on pseudorange-only measurements with adaptive Kalman filter;
Scheme 4: real-time POD based on pseudorange and epoch-differenced carrier-phase measurements with Kalman filter;
Scheme 5: real-time POD based on pseudorange and epoch-differenced carrier-phase measurements with adaptive Kalman filter.
Figure 2, Figure 3, Figure 4 and Figure 5 show real-time POD results of Swarm-A on a normal day. The corresponding statistical solutions are listed in Table 4. It can be seen that overall accuracy is around the sub-meter level in the radial, along and cross directions. Since the observation conditions of Swarm-A are good, a sufficient number of observed satellites provide good conditions for reliable ambiguity resolution with pseudorange and carrier-phase measurements. The accuracy is 1.1, 1.2, and 1.0 m for the radial, along and cross component, respectively. Comparing results of S2 and S3 as well as S4 and S5, it can be seen that the adaptive filter can improve the POD accuracy based on either pseudorange or epoch-differenced carrier-phase measurements, and the improvement with the epoch-differenced carrier phase is more significant than that with pseudorange measurements. The average improvement in three components is 30.2% when comparing solutions from S5 and S4, while it is only 15.5% comparing S3 and S2. This is due to more accurate observational solutions from epoch-differenced carrier-phase measurements; therefore, the adaptive filter can be more successfully applied.
Real-time POD on the day of a solar storm was performed. The results are shown in Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 and Table 5. The results in Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 show the differences between our solutions and the reference orbits. Statistical results in Table 5 are corresponding root mean square (RMS) values. The reference orbit for Swarm-A and GRACE-A are taken from the GeoforschungsZentrum (GFZ) and European Space Agency (ESA) official orbit. The reference orbits are all calculated based on the reduced-dynamic approach, and the accuracy can be better than 5 cm [9,11]. Please note that unlike Figure 5 and Figure 8, the time scales in the horizontal axes in Figure 6, Figure 7, Figure 9 and Figure 10 are plotted compactly since two schemes of results are displayed in these figures. The orbital solutions should use the attitude data to change the GPS antenna phase center to the mass center. For the GRACE satellite, the antenna phase center attitude-four-element file is used. For the Swarm satellite, the Common Data Format (CDF) attitude data file first should be converted in advance (the CDF converter can be found at https://spdf.gsfc.nasa.gov/pub/software/cdf, accessed on 20 May 2021).
From Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10 and Table 5, we can see that in good observation conditions such as Swarm-A, POD results from S1 are better than S2, S3 and S4. However, for GRACE-A, the POD results from S1 are almost the same as S2, S3 and S4. There are two reasons for this. The first is that the Swarm and GRACE IF carrier-phase measurements have a noise level of about 9 mm and are much lower than pseudorange observations. This is why S1 results of Swarm-A are better than S2, S3 and S4. The noise level is calculated based on geometry-free combined observations; details can be found in [28]. The noise levels of spaceborne pseudorange and carrier-phase measurements for other LEO satellites have also been evaluated, i.e., the overall precision of the L1 and L2 measurements are 3.5 mm and 0.8 m for the P1 and P2 measurements of the Luojia-1A satellite [32]. It can be seen that our solutions for the evaluation of spaceborne measurement noise are generally at the same level as other studies.
The second, for GRACE-A, is that the observation conditions are very poor due to the increased atmosphere density at a much lower orbit altitude (~360 km) during the period of a strong solar geomagnetic storm. The conventional approach is not able to estimate the atmospheric drag parameters reliably. Furthermore, the performance of the spaceborne receiver is also seriously affected, as fewer than five satellites are observed during the two periods, and an insufficient number of observations is not beneficial to the resolution of excessive carrier-phase ambiguities or the dynamic parameters. Therefore, POD results of S1 are not good, although the measurement noise of S1 is much lower.
For Swarm-A solutions from S2, the pseudorange noise is at the meter level, and the final orbit solutions are around 1.6, 1.9 and 1.7 m. Using the adaptive Kalman filter in S3 can improve the accuracy by only about 13.6%. This is because pseudorange observations are not sensitive to the state disturbance. The improvement mainly lies at the period of the strong solar storm, which is between 2–6 and 12–18 hours. Comparing S2 and S4, the inclusion of epoch-differenced carrier-phase measurements can improve accuracy by about 22.4% for the radial, along and cross components. When the adaptive filter is applied, the large state disturbance is effectively suppressed, and the accuracy is improved by about 36.5% when comparing S4 and S5.
For GRACE-A, the observation conditions are not good, and on 8 September 2017 there were fewer than five observed satellites for 14% of the observation time (Nsat. in Figure 8). It is difficult to apply the adaptive filter to detect the state disturbance reliably. The orbit accuracy is around 2.2 m for the three components after applying the adaptive filter using the pseudorange measurements. When the epoch-differenced carrier-phase measurements are included, high accuracy orbit solutions can also be passed to the next epoch using the epoch-differenced orbit increments when there are only about four tracked satellites. At the same time, the orbit dynamic model error can also be precisely detected, and the orbit accuracy is improved by about 31.1% when the adaptive filter is applied.
Solutions of adaptive factors for Swarm-A and GRACE-A are also displayed in Figure 11 to demonstrate the above analyses. POD solutions from S3 are not sensitive to the state disturbance due to the limited accuracy of pseudorange measurement-based solutions. The inclusion of epoch-differenced carrier-phase measurements can detect the state distance more reliably. This can be viewed by comparing the left and right parts of subplots in Figure 11. More dense dots indicate that more adaptive factors are involved in the adjustment. Furthermore, it can also be seen that the adaptive filter with observations of Swarm-A performed much better than that of GRACE-A due to its good observation conditions, and the more applied adaptive factors indicate that more state distances are reliably detected.
In sum, the pseudorange and epoch-differenced carrier-phase measurements can achieve one-meter-level real-time orbit accuracy, and the accuracy can be further improved to the sub-meter level when the adaptive filter is applied under good observation conditions. However, the adaptive Kalman filter plays a minor role of suppressing the dynamic model error for pseudorange-only solutions when there is an insufficient number of observed satellites. This situation can be improved with the inclusion of epoch-differenced carrier-phase measurements. Application of the epoch-differenced carrier-phase-based adaptive Kalman filter to various LEO real-time POD missions would be worthwhile.

5. Discussion

This study discusses the performance of spaceborne pseudorange and epoch-differenced carrier-phase measurement-based real-time POD of LEO satellites. This approach shows an apparent advantage when the observation conditions are poor. This is achieved by reducing the number of ambiguities when making the epoch-differenced carrier-phase measurements. However, epoch-differencing requires satellites to be observed continuously. When the number of continuous observed satellites is not sufficient, i.e., fewer than 4, the Kalman filter is re-initialized. At this moment, pseudorange observations play a major role, and SPP is applied in the filter re-initialization process. It is therefore difficult to detect the state disturbance, and the adaptive filter is not applicable. It is still a challenge to deal with the problem of insufficient observations with existing approaches, and this needs to be explored in the future.
From Figure 6, Figure 7, Figure 9 and Figure 10, we can see that the orbit errors in the radial, along, cross directions all have significant periodic errors, and the period is the same as the orbital period. This is because in real-time processing, the periodic phenomenon appears when the float ambiguity is resolved epoch-by-epoch, unlike the reduced-dynamic-based post-processing in which the periodic errors are no longer significant after the dynamic smoothing.

6. Conclusions

In this study, we assessed the LEO real-time POD performance with pseudorange and epoch-differenced carrier-phase measurements. The epoch-differenced carrier-phase-based real-time POD models and algorithms are proposed to deal with LEO real-time POD under poor observation conditions. Five schemes including the pseudorange and carrier-phase-based, pseudorange-only, pseudorange and epoch-differenced carrier-phase-based real-time POD and the application of the adaptive filter were designed to evaluate the performance of the proposed approach.
Firstly, the real-time POD of Swarm-A on a normal day was calculated to evaluate the basic performance of the proposed approach. About 1.5 m accuracy real-time orbital solutions in three components were obtained with pseudorange-only observations. With the inclusion of epoch-differenced carrier-phase measurements, there is a slight improvement by about 18% in the average accuracy is achieved. More importantly, the adaptive filter can both improve the accuracy with either pseudorange-based observations or the inclusion of epoch-differenced carrier-phase measurements, and the improvement is more significant with the latter.
Then, the performances of the proposed approach on Swarm-A and GRACE-A real-time POD on a day with strong solar storm were evaluated. Real calculations show that the pseudorange and epoch-differenced carrier-phase measurement-based orbital solutions can achieve almost the same accuracy as those calculated with IF pseudorange and carrier-phase measurements but with much less computation, and this approach is adaptive when there are few (only four) observed satellites. When the observation conditions are good, 1–2 m real-time orbit accuracy in radial, along and cross directions can be obtained with the proposed approach. The adaptive Kalman filter can further improve the accuracy to the sub-meter level. The pseudorange and epoch-differenced carrier-phase measurement-based solution is more sensitive to the state disturbance than the pseudorange-only observations, and the adaptive Kalman filter can be successfully applied based on such observation combinations. Furthermore, the adaptive filter can contribute more significantly to the improvement of real-time POD of LEO satellites when there is a sufficient number of observed satellites. Results demonstrate that the adaptive filter can still improve the real-time orbital accuracy by 25.6%, 31.1% and 28.0%, respectively, in the radial, along and cross directions even in poor observation conditions, i.e., fewer than five GPS satellites are observed in 14% of the full observation time.

Author Contributions

Conceptualization, M.L. and T.X.; methodology, M.L.; software, Y.S. and K.W.; formal analysis, M.L. and T.X; investigation, X.F. and D.W.; resources, Y.S. and K.W.; writing—original draft preparation, M.L.; writing—review and editing, M.L. and Y.S.; funding acquisition, T.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (2020YFB0505800 and 2020YFB0505804), Postdoctoral Science Foundation of China (2021M691902) and National Natural Science Foundation of China (Grant No. 41874032).

Acknowledgments

We are very grateful to the International GNSS Service (IGS) for providing the broadcast ephemeris (ftp://igs.ign.fr/pub/igs/data/, accessed on 1 May, 2021). The Swarm and GRACE spaceborne GPS data and official orbits provided by the European Space Agency (ESA) as part of the Copernicus program and GFZ, which can be accessed at ftp://swarm-diss.eo.esa.int/Level1b/, accessed on 2 May 2021 and https://isdc.gfz-potsdam.de/grace-isdc/, accessed on 2 May 2021, are greatly acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Reigber, C.; Luhr, H.; Schwintzer, P.; Wickert, J. Earth Observation with CHAMP: Results from Three Years in Orbit; Springer: Berlin, Germany, 2005. [Google Scholar]
  2. Massmann, F.-H.; Beerer, J.; Tapley, B.; Reigber, C. GRACE Mission Status and Future Plans, (Geophysical Research Abstracts, Vol. 9, 07022, 2007). In Proceedings of the General Assembly European Geosciences Union (EGU), Vienna, Austria, 15–20 April 2007. [Google Scholar]
  3. Kornfeld, R.P.; Arnold, B.W.; Gross, M.A.; Dahya, N.T.; Klipstein, W.M.; Gath, P.; Bettadpur, S. GRACE-FO: The Gravity Recovery and Climate Experiment Follow-On Mission. J. Spacecr. Rockets 2019, 56, 931–951. [Google Scholar] [CrossRef]
  4. Garcia, E.S.; Sandwell, D.T.; Smith, W. Retracking CryoSat-2, Envisat and Jason-1 radar altimetry waveforms for improved gravity field recovery. Geophys. J. Int. 2014, 196, 2. [Google Scholar] [CrossRef] [Green Version]
  5. Hwang, C.; Tseng, T.P.; Lin, T. Precise orbit determination for the FORMOSAT-3/COSMIC satellite mission using GPS. J. Geod. 2009, 83, 477–489. [Google Scholar] [CrossRef]
  6. Olsen, N.; Eigil, F.C.; Rune, F. The Swarm Satellite Constellation Application and Research Facility (SCARF) and Swarm data products. Earth Planets Space 2013, 65, 1189–1200. [Google Scholar] [CrossRef]
  7. Yuan, Y.; Mi, X.; Zhang, B. Initial assessment of single-and dual-frequency BDS-3 RTK positioning. Satell. Navig. 2020, 1, 7. [Google Scholar] [CrossRef]
  8. Bock, H.; Jäggi, A.; Meyer, U.; Visser, P.; Ijssel, J.; Helleputte, T.V.; Markus, H.; Hugentobler, U. GPS-derived orbits for the GOCE satellite. J. Geod. 2011, 85, 807–818. [Google Scholar] [CrossRef] [Green Version]
  9. Kang, Z.; Tapley, B.; Bettadpur, S.; Ries, J.; Nagel, P.; Pastor, R. Precise orbit determination for the GRACE mission using only GPS data. J. Geod. 2006, 80, 322–331. [Google Scholar] [CrossRef]
  10. Montenbruck, O.; Hackel, S.; Jaggi, A. Precise orbit determination of the Sentinel-3A altimetry satellite using ambiguity-fixed GPS carrier phase observations. J. Geod. 2018, 92, 711–726. [Google Scholar] [CrossRef] [Green Version]
  11. Montenbruck, O.; Hackel, S.; van den Ijssel, J.; Arnold, D. Reduced dynamic and kinematic precise orbit determination for the Swarm mission from 4 years of GPS tracking. GPS Solut. 2018, 22, 79. [Google Scholar] [CrossRef]
  12. Yoon, Y.T.; Eineder, M.; Yague-Martinez, N.; Montenbruck, O. TerraSAR-X precise trajectory estimation and quality assessment. IEEE Trans. Geosci. Remote Sens. 2009, 47, 1859–1868. [Google Scholar] [CrossRef]
  13. Li, M.; Li, W.; Shi, C.; Jiang, K.; Guo, X.; Dai, X.; Meng, X.; Yang, Z.; Yang, G. Precise orbit determination of the Fengyun-3C satellite using onboard GPS and BDS observations. J. Geod. 2017, 91, 1313–1327. [Google Scholar] [CrossRef] [Green Version]
  14. Montenbruck, O.; Ramos-Bosch, P. Precision real-time navigation of LEO satellites using global positioning system measurements. GPS Solut. 2008, 12, 187–198. [Google Scholar] [CrossRef]
  15. Yang, Y.X.; Mao, Y.; Sun, B.J. Basic performance and future developments of BeiDou global navigation satellite system. Satell. Navig. 2020, 1, 1. [Google Scholar] [CrossRef] [Green Version]
  16. Haas, L.; Pittelkau, M. Real-time high accuracy GPS onboard orbit determination for use on remote sensing satellites. In Proceedings of the ION GPS 99, Institute of Navigation, Nashville, TN, USA, 14–17 September 1999; pp. 829–836. [Google Scholar]
  17. Loiselet, M.; Stricker, N.; Menard, Y.; Luntama, J.P. GRAS—MetOp’s GPS-based atmospheric sounder. ESA Bull. 2000, 102, 38–44. [Google Scholar]
  18. Williams, J.; Lightsey, E.G.; Yoon, S.P.; Schutz, R.E. Testing of the ICESat BlackJack GPS receiver engineering model. In Proceedings of the ION-GPS-2002, Portland, OR, USA, 24–27 September 2002. [Google Scholar]
  19. Feng, Y. An alternative orbit integration algorithm for GPS based precise LEO autonomous navigation. GPS Solut. 2001, 5, 1–11. [Google Scholar] [CrossRef]
  20. Chen, J.P.; Wang, J.X. Kinematic precise orbit determination of low earth orbiter based on epoch-difference strategy. J. Geod. Geodyn. 2007, 27, 57–61. [Google Scholar]
  21. Zhang, K.; Li, X.; Xiong, C.; Meng, X.; Li, X.; Yuan, Y.; Zhang, X. The Influence of Geomagnetic Storm of 7 September 2017 on the Swarm Precise Orbit Determination. J. Geophys. Res. Space Phys. 2019, 124, 6971–6984. [Google Scholar] [CrossRef]
  22. Montenbruck, O.; Kroes, R. In-flight performance analysis of the CHAMP BlackJack GPS Receiver. GPS Solut. 2003, 7, 74–86. [Google Scholar] [CrossRef]
  23. Tang, J.; Wang, H.; Chen, Q.; Chen, Z.; Zheng, J.; Cheng, H.; Liu, L. A time-efficient implementation of Extended Kalman Filter for sequential orbit determination and a case study for onboard application. Adv. Space Res. 2018, 62, 343–358. [Google Scholar] [CrossRef]
  24. Salazar, D.; Hernandez-Pajares, M.; Juan-Zornoza, J.M.; Sanz-Subirana, J.; Aragon-Angel, A. EVA: GPS-based extended velocity and acceleration determination. J. Geod. 2011, 85, 329–340. [Google Scholar] [CrossRef] [Green Version]
  25. Li, M.; Xu, T.; Lu, B.; He, K. Multi-GNSS precise orbit positioning for airborne gravimetry over Antarctica. GPS Solut. 2019, 23, 53. [Google Scholar] [CrossRef]
  26. Li, M.; Xu, T.; Flechtner, F.; Foerste, C.; Lu, B.; He, K. Improving the Performance of Multi-GNSS (Global Navigation Satellite System) Ambiguity Fixing for Airborne Kinematic Positioning over Antarctica. Remote Sens. 2019, 11, 992. [Google Scholar] [CrossRef] [Green Version]
  27. Wagner, C.; Mcadoo, D.; Klokocnik, J. Degradation of Geopotential Recovery from Short Repeat-Cycle Orbits: Application to GRACE Monthly Fields. J. Geod. 2006, 80, 94–103. [Google Scholar] [CrossRef]
  28. Li, M.; Xu, T.; Ge, H.; Guan, M.; Yang, H.; Fang, Z.; Gao, F. LEO-Constellation-Augmented BDS Precise Orbit determination Considering Spaceborne Observational Errors. Remote Sens. 2021, 13, 3189. [Google Scholar] [CrossRef]
  29. Förste, C.; Bruinsma, S.; Shako, R.; Marty, J.C.; Flechtner, F.; Abrykosov, O.; Dahle, C.; Lemoine, J.-M.; Neumayer, K.-H.; Biancale, R. EIGEN-6—A New Combined Global Gravity Field Model Including GOCE Data from the Collaboration of GFZ Potsdam and GRGS Toulouse (Geophysical Research Abstracts Vol.13, EGU2011-3242-2, 2011). In Proceedings of the General Assembly European Geosciences Union, Vienna, Austria, 3–8 April 2011. [Google Scholar]
  30. Petit, G.; Luzum, B. IERS Conventions 2010; No.36 in IERS Technical Note; Verlag des Bundesamts für Kartographie und Geodäsie: Frankfurt am Main, Germany, 2010. [Google Scholar]
  31. Lyard, F.; Lefevre, F.; Letellier, T.; Francis, O. Modelling the global ocean tides: Modern insights from FES2004. Ocean Dyn. 2006, 56, 394–415. [Google Scholar] [CrossRef]
  32. Wang, L.; Chen, R.; Li, D.; Zhang, G.; Shen, X.; Yu, B.; Wu, C.; Xie, S.; Zhang, P.; Li, M.; et al. Initial assessment of the LEO based navigation signal augmentation system from Luojia-1A satellite. Sensors 2018, 18, 3919. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Diagram of comparison of orbit discrepancies detected with different types of observations (left: pseudorange; right: pseudorange and epoch-differenced carrier-phase).
Figure 1. Diagram of comparison of orbit discrepancies detected with different types of observations (left: pseudorange; right: pseudorange and epoch-differenced carrier-phase).
Remotesensing 14 02273 g001
Figure 2. Swarm-A real-time POD results on a normal day of Scheme 1.
Figure 2. Swarm-A real-time POD results on a normal day of Scheme 1.
Remotesensing 14 02273 g002
Figure 3. Swarm-A real-time POD results on a normal day of Schemes 2 (left) and 3 (right).
Figure 3. Swarm-A real-time POD results on a normal day of Schemes 2 (left) and 3 (right).
Remotesensing 14 02273 g003
Figure 4. Swarm-A real-time POD results on a normal day of Schemes 4 (left) and 5 (right).
Figure 4. Swarm-A real-time POD results on a normal day of Schemes 4 (left) and 5 (right).
Remotesensing 14 02273 g004
Figure 5. Swarm-A real-time POD results of Scheme 1.
Figure 5. Swarm-A real-time POD results of Scheme 1.
Remotesensing 14 02273 g005
Figure 6. Swarm-A real-time POD results of Schemes 2 (left) and 3 (right).
Figure 6. Swarm-A real-time POD results of Schemes 2 (left) and 3 (right).
Remotesensing 14 02273 g006
Figure 7. Swarm-A real-time POD results of Schemes 4 (left) and 5 (right).
Figure 7. Swarm-A real-time POD results of Schemes 4 (left) and 5 (right).
Remotesensing 14 02273 g007
Figure 8. GRACE-A real-time POD results of Scheme 1.
Figure 8. GRACE-A real-time POD results of Scheme 1.
Remotesensing 14 02273 g008
Figure 9. GRACE-A real-time POD results of Schemes 2 (left) and 3 (right).
Figure 9. GRACE-A real-time POD results of Schemes 2 (left) and 3 (right).
Remotesensing 14 02273 g009
Figure 10. GRACE-A real-time POD results of Schemes 4 (left) and 5 (right).
Figure 10. GRACE-A real-time POD results of Schemes 4 (left) and 5 (right).
Remotesensing 14 02273 g010
Figure 11. The adaptive factors for Swarm-A and GRACE-A real-time POD solutions from Schemes 3 and 5 on the day of a solar storm. The two upper subplots are for Swarm-A while the two lower subplots are for GRACE-A. The two left subplots are from Scheme 3 (pseudorange-only) and the right are from Scheme 5 (pseudorange and epoch-differenced carrier-phase-based).
Figure 11. The adaptive factors for Swarm-A and GRACE-A real-time POD solutions from Schemes 3 and 5 on the day of a solar storm. The two upper subplots are for Swarm-A while the two lower subplots are for GRACE-A. The two left subplots are from Scheme 3 (pseudorange-only) and the right are from Scheme 5 (pseudorange and epoch-differenced carrier-phase-based).
Remotesensing 14 02273 g011
Table 1. Swarm and GRACE Constellation configurations.
Table 1. Swarm and GRACE Constellation configurations.
LEO ConstellationSwarmGRACE
SatelliteSwarm-A/B/CGRACE-A/B
AltitudeA/C: ~460 km, B: ~510 km~500 km
InclinationA/C: 87.35°, B: 87.75°89.5°
Orbit typeCircular near-polar orbitsCircular near-polar orbits
Repeat cycle7–10 monthsA sparse repeat track of 61 revolutions every 4 days [27]
GoalGeomagnetic observationDetection of the Earth gravity variations
Spaceborne observationsGPSGPS
Sampling interval10 s10 s
Table 2. Parameter settings of dynamic models.
Table 2. Parameter settings of dynamic models.
Dynamic ModelSetting
Earth gravity modelEIGEN-6C (70 × 70) [29]
N-bodyJPL DE405
Solid tide and pole tideIERS 2010 [30]
Ocean tideFES 2004 [31]
RelativelyIERS 2010
Solar radiation pressureMacro Model [11] for both Swarm-A and GRACE-A satellites
Atmospheric dragStatic Harris–Priester density model, fixed superficial area, estimating the drag parameter C D every 4 h.
Empirical accelerationsFirst order Gauss–Markov model, piecewise periodical terms in the along, cross and radial components
Table 3. Settings of the process-noise-related parameters.
Table 3. Settings of the process-noise-related parameters.
ParameterInitial VarianceSteady State VarianceCorrelation Time
Position (m)1.0--
Velocity (m/s)1.0--
Receiver clock offset (m)500.050.030.0
Empirical force acceleration in radial (nm/s2)100.0200.02000.0
Empirical force acceleration in track (nm/s2)400.0800.02000.0
Empirical force acceleration in normal (nm/s2)200.0400.02000.0
Table 4. Statistical results of Swarm-A orbit solutions on a normal day from different schemes (unit: m).
Table 4. Statistical results of Swarm-A orbit solutions on a normal day from different schemes (unit: m).
SchemeS1S2S3S4S5
Swarm-ARadial1.091.381.111.080.70
Along1.231.581.391.170.88
Cross0.971.451.181.220.78
Table 5. Statistical results of Swarm-A and GRACE-A orbit solutions from different schemes (unit: m).
Table 5. Statistical results of Swarm-A and GRACE-A orbit solutions from different schemes (unit: m).
SchemeS1S2S3S4S5
Swarm-ARadial1.211.611.301.220.91
Along1.451.921.821.571.06
Cross1.041.711.391.400.86
GRACE-ARadial2.252.862.111.941.35
Along2.293.142.412.461.75
Cross2.362.802.212.011.47
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, M.; Xu, T.; Shi, Y.; Wei, K.; Fei, X.; Wang, D. Adaptive Kalman Filter for Real-Time Precise Orbit Determination of Low Earth Orbit Satellites Based on Pseudorange and Epoch-Differenced Carrier-Phase Measurements. Remote Sens. 2022, 14, 2273. https://doi.org/10.3390/rs14092273

AMA Style

Li M, Xu T, Shi Y, Wei K, Fei X, Wang D. Adaptive Kalman Filter for Real-Time Precise Orbit Determination of Low Earth Orbit Satellites Based on Pseudorange and Epoch-Differenced Carrier-Phase Measurements. Remote Sensing. 2022; 14(9):2273. https://doi.org/10.3390/rs14092273

Chicago/Turabian Style

Li, Min, Tianhe Xu, Yali Shi, Kai Wei, Xianming Fei, and Dixing Wang. 2022. "Adaptive Kalman Filter for Real-Time Precise Orbit Determination of Low Earth Orbit Satellites Based on Pseudorange and Epoch-Differenced Carrier-Phase Measurements" Remote Sensing 14, no. 9: 2273. https://doi.org/10.3390/rs14092273

APA Style

Li, M., Xu, T., Shi, Y., Wei, K., Fei, X., & Wang, D. (2022). Adaptive Kalman Filter for Real-Time Precise Orbit Determination of Low Earth Orbit Satellites Based on Pseudorange and Epoch-Differenced Carrier-Phase Measurements. Remote Sensing, 14(9), 2273. https://doi.org/10.3390/rs14092273

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop