[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

Search Results (298)

Search Parameters:
Keywords = precise orbit determination

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
26 pages, 13225 KiB  
Article
A New and Tested Ionospheric TEC Prediction Method Based on SegED-ConvLSTM
by Yuanhang Liu, Yingkui Gong, Hao Zhang, Ziyue Hu, Guang Yang and Hong Yuan
Remote Sens. 2025, 17(5), 885; https://doi.org/10.3390/rs17050885 - 2 Mar 2025
Viewed by 337
Abstract
Total electron content (TEC) serves as a key parameter characterizing ionospheric conditions. Accurate prediction of TEC plays a crucial role in improving the precision of Global Navigation Satellite Systems (GNSS). However, existing research have predominantly emphasized spatial variations in the ionosphere, neglecting the [...] Read more.
Total electron content (TEC) serves as a key parameter characterizing ionospheric conditions. Accurate prediction of TEC plays a crucial role in improving the precision of Global Navigation Satellite Systems (GNSS). However, existing research have predominantly emphasized spatial variations in the ionosphere, neglecting the periodic changes of the ionosphere with the diurnal cycle. In this paper, we propose a TEC prediction model, which simultaneously considers both spatial and temporal characteristics to extract spatiotemporal features of ionospheric distribution. Additionally, we integrate several space weather element datasets into the prediction model framework, allowing the generation of multiple space weather feature values that represent the influence of space weather on the ionosphere at different latitudes and longitudes. Moreover, we apply Gaussian process regression (GPR) interpolation to geomagnetic data to characterize impact on the ionosphere, thereby enhancing the prediction accuracy. We compared our model with traditional image-based models such as convolutional neural networks (CNNs), convolutional long short-term memory networks (ConvLSTMs), a self-attention mechanism-integrated ConvLSTM (SAM-ConvLSTM) model, and one-day predicted ionospheric products (C1PG) provided by the Center for Orbit Determination in Europe (CODE). We also examined the effect of using different numbers of space weather feature values in these models. Our model outperforms the comparison models in terms of prediction error metrics, including mean absolute error (MAE), root mean square error (RMSE), correlation coefficient (CC), and the structural similarity index (SSIM). Furthermore, we analyzed the influence of different batch sizes on model training accuracy to find the best performance of each model. In addition, we investigated the model performance during geomagnetic quiet periods, where our model provided the most accurate predictions and demonstrates higher prediction accuracy in the equatorial anomaly region. We also analyzed the prediction performance of all models during space weather events. The results indicate that the proposed model is the least affected during geomagnetic storms and demonstrates superior prediction performance compared to other models. This study presents a more stable and high-performance spatiotemporal prediction model for TEC. Full article
Show Figures

Figure 1

Figure 1
<p>The time series of different space weather event elements and total TEC unit per hour during 2012–2024.</p>
Full article ">Figure 2
<p>The sliding window method used in this study divides continuous time segments into 60-day periods, with each period allocated as follows: 40 days for the training set and 20 days for the test set.</p>
Full article ">Figure 3
<p>SegED-ConvLSTM architecture.</p>
Full article ">Figure 4
<p>Global TEC distribution at the same time points for two consecutive days.</p>
Full article ">Figure 5
<p>ConvLSTM Structure.</p>
Full article ">Figure 6
<p>ED-ConvLSTM structure.</p>
Full article ">Figure 7
<p>Model Prediction Under Space Weather Anomaly Events.</p>
Full article ">Figure 8
<p>2014-05-12 Comparison between GroundTruth, TS-ConvLSTM-6, SegED-ConvLSTM-1, SegED-ConvLSTM-6.</p>
Full article ">Figure 9
<p>Absolute error distribution of TEC predictions for 2014-05-12 across different models.</p>
Full article ">Figure 10
<p>Comparison of ionospheric TEC predictions: ground truth vs. COPG and SegED-ConvLSTM-6 models.</p>
Full article ">Figure 11
<p>Comparison of SegED-ConvLSTM-6 and COPG at 12 May 2014.</p>
Full article ">
29 pages, 10206 KiB  
Article
Finite-Time Control for Satellite Formation Reconfiguration and Maintenance in LEO: A Nonlinear Lyapunov-Based SDDRE Approach
by Majid Bakhtiari, Amirhossein Panahyazdan and Ehsan Abbasali
Aerospace 2025, 12(3), 201; https://doi.org/10.3390/aerospace12030201 - 28 Feb 2025
Viewed by 333
Abstract
This paper introduces a nonlinear Lyapunov-based Finite-Time State-Dependent Differential Riccati Equation (FT-SDDRE) control scheme, considering actuator saturation constraints and ensuring that the control system operates within safe operational limits designed for satellite reconfiguration and formation-keeping in low Earth orbit (LEO) missions. This control [...] Read more.
This paper introduces a nonlinear Lyapunov-based Finite-Time State-Dependent Differential Riccati Equation (FT-SDDRE) control scheme, considering actuator saturation constraints and ensuring that the control system operates within safe operational limits designed for satellite reconfiguration and formation-keeping in low Earth orbit (LEO) missions. This control approach addresses the challenges of reaching the relative position and velocity vectors within a defined timeframe amid various orbital perturbations. The proposed approach guarantees precise formation control by utilizing a high-fidelity relative motion model that incorporates all zonal harmonics and atmospheric drag, which are the primary environmental disturbances in LEO. Additionally, the article presents an optimization methodology to determine the most efficient State-Dependent Coefficient (SDC) form regarding fuel consumption. This optimization process minimizes energy usage through a hybrid genetic algorithm and simulated annealing (HGASA), resulting in improved performance. In addition, this paper includes a sensitivity analysis to identify the optimized SDC parameterization for different satellite reconfiguration maneuvers. These maneuvers encompass radial, along-track, and cross-track adjustments, each with varying baseline distances. The analysis provides insights into how different parameterizations affect reconfiguration performance, ensuring precise and efficient control for each type of maneuver. The finite-time controller proposed here is benchmarked against other forms of SDRE controllers, showing reduced error margins. To further assess the control system’s effectiveness, an input saturation constraint is integrated, ensuring that the control system operates within safe operational limits, ultimately leading to the successful execution of the mission. Full article
(This article belongs to the Section Astronautics & Space Science)
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of ECI and LVLH frames used in relative motion analysis.</p>
Full article ">Figure 2
<p>Schematic diagram of the deputy satellite relative to the target satellite.</p>
Full article ">Figure 3
<p>The 3D sketch of deputy satellite trajectory in LVLH frame.</p>
Full article ">Figure 4
<p>The introduced model’s position accuracy compared to the ERM Model.</p>
Full article ">Figure 5
<p>The introduced model’s velocity accuracy compared to the ERM Model.</p>
Full article ">Figure 6
<p>The variation in the optimal values of <math display="inline"><semantics> <mrow> <mi>β</mi> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>γ</mi> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <mi>δ</mi> </mrow> </semantics></math> in the scenarios with radial motion only (<b>a</b>), along-track motion only (<b>b</b>), and cross-track motion only (<b>c</b>).</p>
Full article ">Figure 7
<p>Block diagram of the satellite formation flying control and optimization process.</p>
Full article ">Figure 8
<p>Uncontrolled motion of deputy satellites with respect to the target satellite.</p>
Full article ">Figure 9
<p>The process of cost reduction in optimization through the HGASA.</p>
Full article ">Figure 10
<p>The absolute sum of the magnitudes of the control forces for four deputy satellites using three types of SDC and the optimized form.</p>
Full article ">Figure 11
<p>The 3D sketch of deputy satellites’ formation reconfiguration and maintenance trajectory utilizing the Lyapunov-based FT-SDDRE method (In this figure, The filled circles indicate the deputy satellites’ initial positions, while the hollow circles represent their positions at the end of the mission).</p>
Full article ">Figure 12
<p>The 2D sketch of deputy satellites’ formation reconfiguration and maintenance trajectory utilizing the Lyapunov-based FT-SDDRE method in different perspectives (In this figure, The filled circles indicate the deputy satellites’ initial positions, while the hollow circles represent their positions at the end of the mission).</p>
Full article ">Figure 13
<p>The position tracking of deputy satellites utilizing the Lyapunov-based FT-SDDRE method.</p>
Full article ">Figure 14
<p>The velocity tracking of deputy satellites utilizing the Lyapunov-based SDDRE method.</p>
Full article ">Figure 15
<p>The relative distance among deputy satellites in the projected circular orbit formation after utilizing the Lyapunov-based FT-SDDRE method.</p>
Full article ">Figure 16
<p>The relative distance among deputy satellites in <math display="inline"><semantics> <mrow> <mi>y</mi> <mo>−</mo> <mi>z</mi> </mrow> </semantics></math> plane in the projected circular orbit formation after utilizing the Lyapunov-based FT-SDDRE method.</p>
Full article ">Figure 17
<p>The absolute position tracking error of the Lyapunov-based FT-SDDRE method compared to the classical SDRE and finite-time STM approach by deputy satellites.</p>
Full article ">Figure 18
<p>The absolute velocity tracking error of the Lyapunov-based FT-SDDRE method compared to the classical SDRE and finite-time STM approach by deputy satellites.</p>
Full article ">Figure 19
<p>The control force generated by Lyapunov-based FT-SDDRE controller for deputy satellites.</p>
Full article ">
17 pages, 3270 KiB  
Review
Progress in Atmospheric Density Inversion Based on LEO Satellites and Preliminary Experiments for SWARM-A
by Xiaoyu Bian, Cunying Xiao, Shuli Song and Mengjun Wu
Remote Sens. 2025, 17(5), 793; https://doi.org/10.3390/rs17050793 - 24 Feb 2025
Viewed by 136
Abstract
The vigorous development of Low Earth Orbit (LEO) satellite constellation programs imposes higher requirements for the accuracy of satellite orbit determination. Significant variations in atmospheric density within the operational region of LEO satellites are primary factors influencing their orbital decay and operational lifespan. [...] Read more.
The vigorous development of Low Earth Orbit (LEO) satellite constellation programs imposes higher requirements for the accuracy of satellite orbit determination. Significant variations in atmospheric density within the operational region of LEO satellites are primary factors influencing their orbital decay and operational lifespan. This article first summarizes the research advancements in atmospheric density inversion utilizing LEO satellites, comparing and analyzing the principles of various algorithms, factors affecting accuracy, as well as the advantages and disadvantages associated with different acquisition methods. Subsequently, we introduce recent progress in enhancing atmospheric density inversion algorithms and data analysis applications based on LEO satellites. The SWARM-A satellite, equipped with a high-precision GPS receiver and accelerometer, was employed to invert atmospheric density using both semi-long axis attenuation and accelerometer methodologies. The inversion results were compared against empirical models to validate their reliability; specifically, the correlation coefficient between the semi-long axis attenuation method and nrlmsise00 reached 0.9158, while that between the accelerometer method and nrlmsise00 attained 0.9204. Notably, the inversion accuracy achieved by the accelerometer slightly surpasses that of the semi-long axis attenuation method. These findings provide valuable support for predicting large air tightness based on LEO satellite orbit data inversions and for adjusting operational orbits to ensure successful execution of satellite missions. Full article
(This article belongs to the Section Atmospheric Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>Search for articles using keywords “LEO Satellite” and “Low Earth Orbit Satellite” in the Web of Science database from 2015 to 2024 and keyword co-occurrence/clustering maps created using CiteSpace6.3.1 software in order to analyze research centered on LEO satellites and find the most important factors in the research of LEO satellites. In this figure, all words that co-occurred with LEO satellite keywords more than 50 times are displayed and automatically clustered according to the research field.</p>
Full article ">Figure 2
<p>Schematic diagram of inversion of orbit data and accelerometer data. The orbit data adopts the specific steps of semi-long-axis attenuation method and energy conservation method, and the accelerometer data are calibrated and inverted.</p>
Full article ">Figure 3
<p>The product of SWARM-A satellite damping factor and effective area is usually expressed by B, the anti-ballistic coefficient. The anti-ballistic coefficient of SWARM-A fluctuates between 2.56 and 2.58.</p>
Full article ">Figure 4
<p>On 28 April 2024, the atmospheric density inversion results of the SWARM-A satellite were obtained, with red indicating the semi-long axis attenuation method and blue indicating the accelerometer inversion. The atmospheric density varied from 1.6 to 3.2, and there was a minimum value at 15 UT (h) in the morning and noon.</p>
Full article ">Figure 5
<p>(<b>a</b>) On 28 April 2024, the atmospheric density inversion results were compared with the NRLMISISE00 model. Orange represents the Nrlmisise00 simulated density, green represents the semi-long axis attenuation method for atmospheric density inversion, and blue represents the accelerometer inversion. (<b>b</b>) We take the hourly average of atmospheric density on 28 April 2024. Blue represents the accelerometer average, green is the hourly average corresponding to the semi-long axis attenuation method, and orange is the hourly average density of the Nrlmisise00 model.</p>
Full article ">Figure 6
<p>We linearly fit the semi-long axis attenuation and velocity meter inversion results in <a href="#remotesensing-17-00793-f005" class="html-fig">Figure 5</a>b with the empirical model, where (<b>a</b>) represents the fitting results of the accelerometer and nrlmisise00, and (<b>b</b>) represents the fitting effect of the semi-long axis attenuation and nrlmisise00.</p>
Full article ">
21 pages, 11665 KiB  
Article
Influences of Discontinuous Attitudes on GNSS/LEO Integrated Precise Orbit Determination Based on Sparse or Regional Networks
by Yuanxin Wang, Baoqi Sun, Kan Wang, Xuhai Yang, Zhe Zhang, Minjian Zhang and Meifang Wu
Remote Sens. 2025, 17(4), 712; https://doi.org/10.3390/rs17040712 - 19 Feb 2025
Viewed by 178
Abstract
A uniformly distributed global ground network is essential for the accurate determination of GNSS orbit and clock parameters. However, achieving an ideal ground network is often difficult. When limited to a sparse or regional network of ground stations, the integration of LEO satellites [...] Read more.
A uniformly distributed global ground network is essential for the accurate determination of GNSS orbit and clock parameters. However, achieving an ideal ground network is often difficult. When limited to a sparse or regional network of ground stations, the integration of LEO satellites can substantially enhance the accuracy of GNSS Precise Orbit Determination (POD). In practical processing, discontinuities with complicated gaps can occur in LEO attitude quaternions, particularly when working with a restricted observation network. This hampers the accuracy of determining GNSS/LEO integrated orbits. To address this, an investigation was conducted using data from seven LEO satellites, including those from Sentinel-3, GRACE-FO, and Swarm, to evaluate integrated POD performance under sparse or regional station conditions. Particular focus was placed on addressing attitude discontinuities. Four scenarios were analyzed, encompassing both continuous data availability and one-, two-, and three-hour interruptions after one hour of continuous data availability. The results showed that the proposed quaternion rotation matrix interpolation method is reliable for the integrated POD of GNSSs and LEOs with strict attitude control. Full article
Show Figures

Figure 1

Figure 1
<p>The rotation angles between the transformed body-fixed frame and the RTA frame for GRACE-C (<b>left</b>) and GRACE-D (<b>right</b>) on DOY 130 of 2021.</p>
Full article ">Figure 2
<p>The rotation angle between the transformed body-fixed frame and the RTA frame for Sentinel-3A (<b>left</b>) and Sentinel-3B (<b>right</b>).</p>
Full article ">Figure 3
<p>Sparse network with 22 stations selected for the integrated POD.</p>
Full article ">Figure 4
<p>Regional network with 7 stations selected for the integrated POD.</p>
Full article ">Figure 5
<p>ODOPs of regional stations and LEOs for GPS satellite per epoch.</p>
Full article ">Figure 6
<p>ODOPs of sparse stations and LEOs for GPS satellite per epoch.</p>
Full article ">Figure 7
<p>The accumulated distribution of GPS satellite orbital errors, including whether or not ERP and geocenter were estimated based on the sparse network. The average orbital accuracies for the radial, along-track, cross-track, and 1D directions for each test are given in the legend.</p>
Full article ">Figure 8
<p>The average orbital accuracy in the radial, along-track, cross-track, and 1D directions for each LEO satellite in the integrated POD based on the sparse network, with and without ERP estimation.</p>
Full article ">Figure 9
<p>The accumulated distribution of GPS satellite orbital errors based on sparse network under different LEO attitude quaternion situations. The average orbital accuracies for the radial, along-track, cross-track, and 1D directions for each test are given in the legend in millimeters.</p>
Full article ">Figure 10
<p>The 1D deviations in the orbital errors of GRACE satellites in integrated POD based on sparse network for different LEO attitude quaternion situations.</p>
Full article ">Figure 11
<p>The 1D deviations in the orbital errors of Swarm satellites in integrated POD based on sparse network for different LEO attitude quaternion situations.</p>
Full article ">Figure 12
<p>The 1D deviations in the orbital errors of Sentinel-3 satellites in integrated POD based on sparse network for different LEO attitude quaternion situations.</p>
Full article ">Figure 13
<p>The 1D mean time series in the orbital errors of GRACE-C satellites in integrated POD based on sparse network for different LEO attitude quaternion situations.</p>
Full article ">Figure 14
<p>The 1D mean time series in the orbital errors of Swarm-A satellites in integrated POD based on sparse network for different LEO attitude quaternion situations.</p>
Full article ">Figure 15
<p>The 1D mean time series in the orbital errors of Sentinel-3B satellites in integrated POD based on sparse network for different LEO attitude quaternion situations.</p>
Full article ">Figure 16
<p>The accumulated distribution of the orbital errors of the GPS satellite for different LEO attitude quaternion situations. The average orbit accuracy for the radial, along-track, cross-track, and 1D directions for each test is given in the legend in millimeters.</p>
Full article ">Figure 17
<p>The 1D deviations in the orbital errors of GRACE satellites in integrated POD based on regional network for different LEO attitude quaternion situations.</p>
Full article ">Figure 18
<p>The 1D deviations in the orbital errors of Swarm satellites in integrated POD based on regional network for different LEO attitude quaternion situations.</p>
Full article ">Figure 19
<p>The 1D deviations in the orbital errors of Sentinel-3 satellites in integrated POD based on regional network for different LEO attitude quaternion situations.</p>
Full article ">
17 pages, 3795 KiB  
Review
Comprehensive Analysis of HY-2B/2C/2D Satellite-Borne GPS Data Quality and Reduced-Dynamic Precise Orbit Determination
by Xin Jin, Guangzhe Wang, Jinyun Guo, Hailong Peng, Yongjun Jia and Xiaotao Chang
Aerospace 2025, 12(2), 102; https://doi.org/10.3390/aerospace12020102 - 30 Jan 2025
Viewed by 492
Abstract
The deployment of the HY-2B/2C/2D satellite constellation marks a significant advancement in China’s marine dynamic environmental satellite program, forming a robust three-satellite network. All satellites are equipped with the “HY2_Receiver”, an indigenous technological achievement. Precise orbit determination using this receiver is critical for [...] Read more.
The deployment of the HY-2B/2C/2D satellite constellation marks a significant advancement in China’s marine dynamic environmental satellite program, forming a robust three-satellite network. All satellites are equipped with the “HY2_Receiver”, an indigenous technological achievement. Precise orbit determination using this receiver is critical for monitoring dynamic oceanic parameters such as sea surface wind fields and heights. This study presents a detailed analysis and comparison of the GPS data quality from the HY-2B/2C/2D satellites, emphasizing the impact of phase center variation (PCV) model corrections on orbit accuracy, with a particular focus on high-precision reduced-dynamic orbit determination. The experimental results demonstrate that the GPS data from the satellites exhibit consistent satellite visibility and minimal multipath errors, confirming the reliability and stability of the receivers. Incorporating PCV model corrections significantly enhances orbit accuracy, achieving improvements of approximately 0.3 cm. Compared to DORIS-derived orbits from the Centre National d’Études Spatiales (CNES), the GPS-derived reduced-dynamic orbits consistently reach radial accuracies of 1.5 cm and three-dimensional accuracies of 3 cm. Furthermore, validation using Satellite Laser Ranging (SLR) data confirms orbit accuracies better than 3.5 cm, with 3D root mean square (RMS) accuracies exceeding 3 cm in the radial (R), along-track (T), and cross-track (N) directions. Notably, the orbit determination accuracy remains consistent across all satellites within the HY-2B/2C/2D constellation. This comprehensive analysis highlights the consistent and reliable performance of the indigenous “HY2_Receiver” in supporting high-precision orbit determination for the HY-2B/2C/2D constellation, demonstrating its capability to meet the rigorous demands of marine dynamic environmental monitoring. Full article
Show Figures

Figure 1

Figure 1
<p>Satellite and its payloads.</p>
Full article ">Figure 2
<p>Attitude variation features of HY-2B/2C/2D satellites.</p>
Full article ">Figure 3
<p>Change in F10.7 values.</p>
Full article ">Figure 4
<p>Number of single epoch observation satellites.</p>
Full article ">Figure 5
<p>Multipath error ((<b>a</b>) HY-2B satellite MP1, (<b>b</b>) HY-2B satellite MP2, (<b>c</b>) HY-2C satellite MP1, (<b>d</b>) HY-2C satellite MP2, (<b>e</b>) HY-2D satellite MP1, (<b>f</b>) HY-2D satellite MP2).</p>
Full article ">Figure 6
<p>PCV model ((<b>a</b>) HY-2B Residual method PCV model, (<b>b</b>) HY-2B Direct method PCV model, (<b>c</b>) HY-2C Residual method PCV model, (<b>d</b>) HY-2C direct method PCV model, (<b>e</b>) HY-2D Residual method PCV model, (<b>f</b>) HY-2D Direct method PCV model).</p>
Full article ">Figure 7
<p>Intercomparison between RD and CNES orbits in the R, T, and N directions.</p>
Full article ">Figure 8
<p>NP data quantity of SLR station.</p>
Full article ">
18 pages, 5909 KiB  
Communication
High-Speed Target Location Based on Photoelectric Imaging and Laser Ranging with Fast Steering Mirror Deflection
by Kui Shi, Hongtao Yang, Jianwei Peng, Yingjun Ma and Hongwei Zhang
Photonics 2025, 12(2), 108; https://doi.org/10.3390/photonics12020108 - 24 Jan 2025
Viewed by 653
Abstract
There is an increasing number of spacecrafts in orbit, and the collision impact of high-speed moving targets, such as space debris, can cause fatal damage to these spacecrafts. It has become increasingly important to rapidly and accurately locate high-speed moving targets in space. [...] Read more.
There is an increasing number of spacecrafts in orbit, and the collision impact of high-speed moving targets, such as space debris, can cause fatal damage to these spacecrafts. It has become increasingly important to rapidly and accurately locate high-speed moving targets in space. In this study, we designed a visible-light telephoto camera for observing high-speed moving targets and a laser rangefinder for measuring the precise distance of these targets, and we proposed a method of using fast steering mirror deflection to quickly direct the emitted laser towards such targets and measure the distance. Based on the principle of photographic imaging and the precise distance of targets, a collinear equation and a spatial target location model based on the internal and external orientation elements of the camera and the target distance were established, and the principle of target location and the method for calculating target point coordinates were determined. We analyzed the composition of target point location error and derived an equation for calculating such errors. Based on the actual values of various error components and the error synthesis theory, the accuracy of target location was calculated to be 26.5 m when the target distance is 30 km (the relative velocity is 8 km/s and the velocity component perpendicular to the camera’s optical axis is less than 3.75 km/s). This study provides a theoretical basis and a method for solving the practical needs of quickly locating high-speed moving targets in space and proposes specific measures to improve target location accuracy. Full article
Show Figures

Figure 1

Figure 1
<p>Principle of photographic imaging of targets.</p>
Full article ">Figure 2
<p>High-resolution optical system.</p>
Full article ">Figure 3
<p>MTF curve of the optical system.</p>
Full article ">Figure 4
<p>Schematic diagram of the light shading angle in the optical system.</p>
Full article ">Figure 5
<p>Structural composition of the optical objective lens.</p>
Full article ">Figure 6
<p>Image of the camera’s optical objective lens.</p>
Full article ">Figure 7
<p>Structural composition diagram of the target location system.</p>
Full article ">Figure 8
<p>Optical system of the lens for laser collimation.</p>
Full article ">Figure 9
<p>Optical system of the laser-receiving lens.</p>
Full article ">Figure 10
<p>Image of the optical objective lens used in the laser-receiving system.</p>
Full article ">Figure 11
<p>Composition of the high-speed deflection mirror system.</p>
Full article ">Figure 12
<p>Two-dimensional FSM.</p>
Full article ">Figure 13
<p>Schematic diagram of the cooperative target observation and location experiment.</p>
Full article ">Figure 14
<p>An image of the cooperative target, satellite <span class="html-italic">A</span>, captured in orbit.</p>
Full article ">
19 pages, 8948 KiB  
Article
Differential Code Bias Estimation and Accuracy Analysis Based on CSES Onboard GPS and BDS Observations
by Jiawen Pang, Fuying Zhu and Shang Wu
Remote Sens. 2025, 17(3), 374; https://doi.org/10.3390/rs17030374 - 23 Jan 2025
Viewed by 447
Abstract
An accurate estimation of Differential Code Bias (DCB) is essential for high-precision applications of the Global Navigation Satellite System (GNSS) and for the precise determination of GNSS-derived total electron content (TEC). This study leverages BeiDou Navigation Satellite System (BDS) and Global Positioning System [...] Read more.
An accurate estimation of Differential Code Bias (DCB) is essential for high-precision applications of the Global Navigation Satellite System (GNSS) and for the precise determination of GNSS-derived total electron content (TEC). This study leverages BeiDou Navigation Satellite System (BDS) and Global Positioning System (GPS) dual-frequency observations of the China Seismo-electromagnetic Satellite (CSES) from day of the year (DOY) 201 to DOY 232 in 2018, we evaluate the quality of CSES onboard GNSS observations, improve the data preprocessing method, and use the least-squares to estimate DCBs for both GNSS satellites and CSES receivers. A comprehensive analysis of the estimation accuracy is presented, revealing that DCBs for BDS satellites, derived from joint BDS and GPS observations, exhibit superior consistency compared to those from single BDS observations. Notably, the stability of DCBs for the CSES BDS receiver as well as for BDS GEO, IGSO, and MEO satellites has been significantly enhanced by 70%, 14%, 22%, and 23%, respectively. Conversely, the consistency of GPS satellite DCBs estimated from joint observations shows a decline when compared to the DCB products from the Center for Orbit Determination in Europe (CODE) and the Chinese Academy of Sciences (CAS). When fewer than nine satellites are tracked daily and nighttime observations are under 25%, estimation errors increase. The optimal DCB estimation is achieved with a cutoff elevation angle set at 10°, with monthly mean DCB values for CSES GPS and BDS receivers determined to be −2.193 ns and −1.099 ns, respectively, accompanied by root mean square errors (RMSEs) of 0.10 ns and 0.31 ns. The highest accuracy of DCBs estimated by the single-GPS scheme is corroborated by examining the occurrence of negative vertical total electron content (VTEC) percentages. Full article
Show Figures

Figure 1

Figure 1
<p>The statistical result of the number of GPS, BDS, and BDS + GPS satellites before (<b>upper part</b>) and after (<b>lower part</b>) preprocessing.</p>
Full article ">Figure 2
<p>MD and STD of GPS C1WC2W DCB estimates relative to CODE on DOYs 201–232/2018.</p>
Full article ">Figure 3
<p>The monthly MDs between estimated GPS DCB and DCB products from the CODE and CAS.</p>
Full article ">Figure 4
<p>The comparative analysis of the average percentage of nighttime observations (<b>upper part</b>) and the average number of tracked satellites (<b>lower part</b>) with the monthly MDs of GPS DCBs.</p>
Full article ">Figure 5
<p>The monthly STD of GPS DCBs estimated by two solution schemes.</p>
Full article ">Figure 6
<p>The monthly MDs between estimated BDS DCBs and DCB products from the DLR and CAS.</p>
Full article ">Figure 7
<p>The monthly STD of BDS DCBs estimated by two solution schemes.</p>
Full article ">Figure 8
<p>(<b>a</b>) the CSES GOR <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>D</mi> <mi>C</mi> <mi>B</mi> </mrow> <mrow> <mi>r</mi> <mi>G</mi> <mi>P</mi> <mi>S</mi> </mrow> </msub> </mrow> </semantics></math> time series and (<b>b</b>) the CSES GOR <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>D</mi> <mi>C</mi> <mi>B</mi> </mrow> <mrow> <mi>r</mi> <mi>B</mi> <mi>D</mi> <mi>S</mi> </mrow> </msub> </mrow> </semantics></math> time series. The straight lines represent the monthly mean <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>D</mi> <mi>C</mi> <mi>B</mi> </mrow> <mrow> <mi>r</mi> <mi>G</mi> <mi>P</mi> <mi>S</mi> </mrow> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>D</mi> <mi>C</mi> <mi>B</mi> </mrow> <mrow> <mi>r</mi> <mi>B</mi> <mi>D</mi> <mi>S</mi> </mrow> </msub> </mrow> </semantics></math> for different estimation methods.</p>
Full article ">Figure 9
<p>RMSE time series of CSES receiver DCBs at different cutoff elevations. (<b>a</b>) circles represent GPS receivers, and (<b>b</b>) squares represent BDS receivers.</p>
Full article ">Figure 10
<p>RMSE of DCBs for CSES receivers with different cutoff elevations. (<b>a</b>) represent GPS receivers, and (<b>b</b>) represent BDS receivers.</p>
Full article ">Figure 11
<p>VTEC time series for DOY 201/2018 estimated based on different scenarios.</p>
Full article ">
22 pages, 1097 KiB  
Article
Efficient AOA Estimation and NLOS Signal Utilization for LEO Constellation-Based Positioning Using Satellite Ephemeris Information
by Junqi Guo and Yang Wang
Appl. Sci. 2025, 15(3), 1080; https://doi.org/10.3390/app15031080 - 22 Jan 2025
Viewed by 634
Abstract
As large-scale low Earth orbit (LEO) constellations continue to expand, the potential of their signal strength for positioning applications should be fully leveraged. For high-precision angle of arrival (AOA) estimation, current spectrum search algorithms are computationally expensive. To address this, we propose a [...] Read more.
As large-scale low Earth orbit (LEO) constellations continue to expand, the potential of their signal strength for positioning applications should be fully leveraged. For high-precision angle of arrival (AOA) estimation, current spectrum search algorithms are computationally expensive. To address this, we propose a method that downscales the 2D joint spectrum search algorithm by incorporating satellite ephemeris a priori information. The proposed algorithm efficiently and accurately determines the azimuth and elevation angles of NLOS (non-line-of-sight) signals. Furthermore, an NLOS virtual satellite construction method is introduced for integrating NLOS satellite data into the positioning system using previously estimated azimuth and elevation angles. Simulation experiments, conducted with a uniform planar array antenna in environments containing both LOS (line-of-sight) and NLOS signals, demonstrate the effectiveness of the proposed solution. The results show that the azimuth determination algorithm reduces computational complexity without sacrificing accuracy, while the NLOS virtual satellite construction method significantly enhances positioning accuracy in NLOS environments. The geometric dilution of precision (GDOP) improved significantly, decreasing from values exceeding 10 to an average of less than 1.42. Full article
Show Figures

Figure 1

Figure 1
<p>Visual representation for GDOP [<a href="#B29-applsci-15-01080" class="html-bibr">29</a>].</p>
Full article ">Figure 2
<p>Diagram of signal incidence on an antenna array.</p>
Full article ">Figure 3
<p>Spectrum plot of the MUSIC algorithm with different angle search steps.</p>
Full article ">Figure 4
<p>View to demonstrate the angular relationship between the direct and reflected signal. Red and blue line represent reflect and direct paths, respectively.</p>
Full article ">Figure 5
<p>Spectrum plot of the reduced-dimension MUSIC algorithm with different angle search step.</p>
Full article ">Figure 6
<p>NLOS virtual satellite inversion diagram: Red lines represent NLOS paths, blue lines represent unavailable LOS paths, and dashed lines represent direct paths of NLOS inverted virtual satellites.</p>
Full article ">Figure 7
<p>GDOP scatter plot at different times under the condition of no virtual satellites below 60° NLOS.</p>
Full article ">Figure 8
<p>Number of observable satellites under different virtual satellite settings below 60° NLOS condition: no virtual satellites, half of the satellites are set as virtual, one-quarter of the satellites are set as virtual.</p>
Full article ">Figure 9
<p>GDOP under different virtual satellite settings below 60° NLOS condition: half of the satellites are set as virtual, one-quarter of the satellites are set as virtual.</p>
Full article ">Figure 10
<p>GDOP with under different virtual satellite settings below 60° (10/20) NLOS condition: half of the satellites are set as virtual, one-quarter of the satellites are set as virtual.</p>
Full article ">
18 pages, 6373 KiB  
Article
Comparisons and Analyses of Thermospheric Mass Densities Derived from Global Navigation Satellite System–Precise Orbit Determination and an Ionization Gauge–Orbital Neutral Atmospheric Detector Onboard a Spherical Satellite at 520 km Altitude
by Yujiao Jin, Xianguo Zhang, Maosheng He, Yongping Li, Xiangguang Meng, Jiangzhao Ai, Bowen Wang, Xinyue Wang and Yueqiang Sun
Remote Sens. 2025, 17(1), 98; https://doi.org/10.3390/rs17010098 - 30 Dec 2024
Viewed by 589
Abstract
Thermospheric mass densities are investigated to explore their responses to solar irradiance and geomagnetic activity during the period from 31 October to 7 November 2021. Utilizing data from the Global Navigation Satellite System (GNSS) payload and an ionization gauge mounted on the Orbital [...] Read more.
Thermospheric mass densities are investigated to explore their responses to solar irradiance and geomagnetic activity during the period from 31 October to 7 November 2021. Utilizing data from the Global Navigation Satellite System (GNSS) payload and an ionization gauge mounted on the Orbital Neutral Atmospheric Detector (OAD) payload onboard the QQ-Satellite, thermospheric mass densities are derived through two independent means: precise orbit determination (POD) and pressure measurements. For the first time, observations of these two techniques are compared and analyzed in this study to demonstrate similarities and differences. Both techniques exhibit similar spatial–temporal variations, with clear dependences on local solar time (LT). However, the hemispheric asymmetry is almost absent in simulations from the NRLMSISE-00 and DTM94 models compared with observations. At high latitudes, density enhancements of observations and simulations are shown, characterized by periodic bulge structures. In contrast, only the OAD-derived densities exhibit wave-like disturbances that propagate from two poles to lower latitudes during geomagnetic storm periods, suggesting a connection to traveling atmospheric disturbances (TADs). Over the long term, thermospheric mass densities derived from the two means of POD and the OAD show good agreements, yet prominent discrepancies emerge during specific periods and under different space-weather conditions. We propose possible interpretations as well as suggestions for utilizing these two means. Significantly, neutral winds should be considered in both methods, particularly at high latitudes and under storm conditions. Full article
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) The QQ-Satellite, and (<b>b</b>) the projections of the QQ-Satellite orbits in a local solar time–latitude frame.</p>
Full article ">Figure 2
<p>Thermospheric mass densities at (<b>a</b>–<b>d</b>) dawn and (<b>e</b>–<b>h</b>) dusk derived from (<b>a</b>,<b>e</b>) NRLMSISE-00, (<b>b</b>,<b>f</b>) DTM94, (<b>c</b>,<b>g</b>) POD, and (<b>d</b>,<b>h</b>) OAD. The (<b>a</b>,<b>e</b>) black solid line is the daily F10.7 index, and the (<b>a</b>,<b>e</b>) black dashed line is ten times of the daily Kp index.</p>
Full article ">Figure 3
<p>The evolutions of (<b>a</b>) thermospheric mass densities derived from two means and (<b>b</b>) their ratios, along with the evolutions of (<b>c</b>) F10.7 and (<b>d</b>) Kp indices.</p>
Full article ">Figure 4
<p>Normal distributions of ratios displayed in <a href="#remotesensing-17-00098-f003" class="html-fig">Figure 3</a>b for (<b>a</b>) all data, and for three periods defined in <a href="#sec3dot2-remotesensing-17-00098" class="html-sec">Section 3.2</a>: (<b>b</b>) increasing period, (<b>c</b>) decreasing period, and (<b>d</b>) transition period.</p>
Full article ">Figure 5
<p>(<b>a</b>,<b>c</b>,<b>e</b>) Thermospheric mass densities derived from POD and OAD, and (<b>b</b>,<b>d</b>,<b>f</b>) their ratios. The O<sub>2</sub>, O<sub>3</sub>, O<sub>50</sub>, O<sub>51</sub>, O<sub>68</sub>, and O<sub>69</sub> indicate the accumulated numbers of the QQ-Satellite’s orbits since 31 October 2021. And the shadows represent the increasing periods.</p>
Full article ">Figure 6
<p>Similar plot to <a href="#remotesensing-17-00098-f005" class="html-fig">Figure 5</a>, except for different satellite orbits. (<b>a</b>,<b>c</b>,<b>e</b>) Thermospheric mass densities derived from POD and OAD, and (<b>b</b>,<b>d</b>,<b>f</b>) their ratios. And the shadows represent the decreasing periods.</p>
Full article ">Figure 7
<p>Similar plot to <a href="#remotesensing-17-00098-f005" class="html-fig">Figure 5</a>, except for different satellite orbits. (<b>a</b>,<b>c</b>,<b>e</b>,<b>g</b>) Thermospheric mass densities derived from POD and OAD, and (<b>b</b>,<b>d</b>,<b>f</b>,<b>h</b>) their ratios. And the shadows represent the transition periods.</p>
Full article ">Figure 8
<p>The evolutions of (<b>a</b>) peak densities in the orbital bulge structures from simulations and observations and their (<b>b</b>) enhancements, along with the evolution of (<b>c</b>) the AE index.</p>
Full article ">Figure 9
<p>Spatial distributions of peak density in bulge structures at (<b>a</b>–<b>d</b>) dawn (0–8 LT) and at (<b>e</b>–<b>h</b>) dusk (16–24 LT), derived from (<b>a</b>,<b>e</b>) NRLMSISE-00, (<b>b</b>,<b>f</b>) DTM94, (<b>c</b>,<b>g</b>) POD, and (<b>d</b>,<b>h</b>) OAD.</p>
Full article ">Figure 10
<p>Statistical distributions in latitudes for bulge peaks derived from (<b>a</b>) NRLMSISE-00, (<b>b</b>) DTM94, (<b>c</b>) POD, and (<b>d</b>) OAD.</p>
Full article ">Figure 11
<p>Ratios of OAD densities during storm periods on 3 November to 5 November relative to a quiet period on 2 November at (<b>a</b>) dawn, 4–8 LT and at (<b>b</b>) dusk, 16–20 LT.</p>
Full article ">
19 pages, 6733 KiB  
Article
Real-Time Orbit Determination of Micro–Nano Satellite Using Robust Adaptive Filtering
by Jing Chen, Xiaojun Jin, Cong Hou, Likai Zhu, Zhaobin Xu and Zhonghe Jin
Sensors 2024, 24(24), 7988; https://doi.org/10.3390/s24247988 - 14 Dec 2024
Viewed by 663
Abstract
Low-performing GPS receivers, often used in challenging scenarios such as attitude maneuver and attitude rotation, are frequently encountered for micro–nano satellites. To address these challenges, this paper proposes a modified robust adaptive hierarchical filtering algorithm (named IARKF). This algorithm leverages robust adaptive filtering [...] Read more.
Low-performing GPS receivers, often used in challenging scenarios such as attitude maneuver and attitude rotation, are frequently encountered for micro–nano satellites. To address these challenges, this paper proposes a modified robust adaptive hierarchical filtering algorithm (named IARKF). This algorithm leverages robust adaptive filtering to dynamically adjust the distribution of innovation vectors and employs a fading memory weighted method to estimate measurement noise in real time, thereby enhancing the filter’s adaptability to dynamic environments. A segmented adaptive filtering strategy is introduced, allowing for flexible parameter adjustment in different dynamic scenarios. A micro–nano satellite equipped with a miniaturized dual-frequency GPS receiver is employed to demonstrate precise orbit determination capabilities. On-orbit GPS data from the satellite, collected in two specific scenarios—slow rotation and Earth-pointing stabilization—are analyzed to evaluate the proposed algorithm’s ability to cope with weak GPS signals and satellite attitude instability as well as to assess the achievable orbit determination accuracy. The results show that, compared to traditional Extended Kalman Filters (EKF) and other improved filtering algorithms, the IARKF performs better in reducing post-fit residuals and improving orbit prediction accuracy, demonstrating its superior robustness. The three-axes orbit determination internal consistency precision can reach the millimeter level. This work explores a feasible approach for achieving high-performance orbit determination in micro–nano satellites. Full article
(This article belongs to the Section Navigation and Positioning)
Show Figures

Figure 1

Figure 1
<p>Flowchart of IARKF strategy.</p>
Full article ">Figure 2
<p>Comparison of posterior residuals, number of positioning satellites, and GDOP.</p>
Full article ">Figure 3
<p>Comparison of posterior residuals in the sky vision.</p>
Full article ">Figure 4
<p>Comparison of posterior residual scatter.</p>
Full article ">Figure 5
<p>Filtering level variation curve and distribution map.</p>
Full article ">Figure 6
<p>Comparison of errors for overlapping arc.</p>
Full article ">Figure 7
<p>Comparison of posterior residuals, number of positioning satellites, and GDOP.</p>
Full article ">Figure 8
<p>Comparison of posterior residuals in the sky vision.</p>
Full article ">Figure 9
<p>Comparison of posterior residual scatter.</p>
Full article ">Figure 10
<p>Filtering level variation curve and distribution map.</p>
Full article ">Figure 11
<p>Comparison of errors for overlapping arc.</p>
Full article ">
18 pages, 6204 KiB  
Article
An Integrity Monitoring Method for Navigation Satellites Based on Multi-Source Observation Links
by Jie Xin, Dongxia Wang and Kai Li
Remote Sens. 2024, 16(23), 4574; https://doi.org/10.3390/rs16234574 - 6 Dec 2024
Viewed by 621
Abstract
The BeiDou-3 navigation satellite system (BDS-3) has officially provided positioning, navigation, and timing (PNT) services to global users since 31 July 2020. With the application of inter-satellite link technology, global integrity monitoring becomes possible. Nevertheless, the content of integrity monitoring is still limited [...] Read more.
The BeiDou-3 navigation satellite system (BDS-3) has officially provided positioning, navigation, and timing (PNT) services to global users since 31 July 2020. With the application of inter-satellite link technology, global integrity monitoring becomes possible. Nevertheless, the content of integrity monitoring is still limited by the communication capacity of inter-satellite links and the layout of ground monitoring stations. Low earth orbit (LEO) satellites have advantages in information-carrying rate and kinematic velocity and can be used as satellite-based monitoring stations for navigation satellites. Large numbers of LEO satellites can provide more monitoring data than ground monitoring stations and make it easier to obtain full-arc observation data. A new challenge of redundant data also arises. This study constructs multi-source observation links with satellite-to-ground, inter-satellite, and satellite-based observation data, proposes an integrity monitoring method with optimization of observation links, and verifies the performance of integrity monitoring with different observation links. The experimental results show four findings. (1) Based on the integrity status of BDS-3, the proposed system-level integrity mode can realize full-arc anomaly diagnosis in information and signals according to the observation conditions of the target satellite. Apart from basic navigation messages and satellite-based augmentation messages, autonomous messages and inter-satellite ranging data can be used to evaluate the state of the target satellite. (2) For a giant LEO constellation, only a small number of LEO satellites need to be selected to construct a minimum satellite-based observation unit that can realize multiple returns of navigation messages and reduce the redundancy of observation data. With the support of 12 and 30 LEO satellites, the minimum number of satellite-based observation links is 1 and 4, respectively, verifying that a small amount of LEO satellites could be used to construct a minimum satellite-based observation unit. (3) A small number of LEO satellites can effectively improve the observation geometry of the target satellite. An orbit determination observation unit, which consists of chosen satellite-to-ground and/or satellite-based observation links based on observation geometry, is proposed to carry out fast calculations of satellite orbit. If the orbit determination observation unit contains 6 satellite-to-ground monitoring links and 6/12/60 LEO satellites, the value of satellite position dilution of precision (SPDOP) is 38.37, 24.60, and 15.71, respectively, with a 92.95%, 95.49%, and 97.12% improvement than the results using 6 satellite-to-ground monitoring links only. (4) LEO satellites could not only expand the resolution of integrity parameters in real time but also augment the service accuracy of the navigation satellite system. As the number of LEO satellites increases, the area where UDRE parameters can be solved in real time is constantly expanding to a global area. The service accuracy is 0.93 m, 0.88 m, and 0.65 m, respectively, with augmentation of 6, 12, and 60 LEO satellites, which is an 8.9%, 13.7%, and 36.3% improvement compared with the results of regional service. LEO satellites have practical application values by improving the integrity monitoring of navigation satellites. Full article
Show Figures

Figure 1

Figure 1
<p>Design of integrity monitoring system for the BDS-3 satellite.</p>
Full article ">Figure 2
<p>Design of integrity monitoring system with support of multi-source observation links.</p>
Full article ">Figure 3
<p>The chosen valid satellite-based observation links (<math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>≤</mo> <msub> <mi>γ</mi> <mn>0</mn> </msub> </mrow> </semantics></math>).</p>
Full article ">Figure 4
<p>The chosen valid satellite-based observation links (<math display="inline"><semantics> <mrow> <mi>γ</mi> <mo>&gt;</mo> <msub> <mi>γ</mi> <mn>0</mn> </msub> </mrow> </semantics></math>).</p>
Full article ">Figure 5
<p>The chosen valid satellite-to-ground observation links.</p>
Full article ">Figure 6
<p>Processing of integrity monitoring system.</p>
Full article ">Figure 7
<p>Divided grid points and chosen ground monitoring stations.</p>
Full article ">Figure 8
<p>Multiple numbers in scenario 3.</p>
Full article ">Figure 9
<p>Coverage rate of four-multiple numbers in scenario 3.</p>
Full article ">Figure 10
<p>Coverage rate of more than four-multiple numbers in scenario 3. (The grids in the blue area can be monitored by no less than four LEO satellites; the grids in the yellow area can be monitored by less than four LEO satellites).</p>
Full article ">Figure 11
<p>Number of available links and SPDOP values in scenario 1.</p>
Full article ">Figure 12
<p>Number of available links and SPDOP values in scenario 2.</p>
Full article ">Figure 13
<p>Number of available links and SPDOP values in scenario 3.</p>
Full article ">Figure 14
<p>Number of available links and SPDOP values in scenario 4.</p>
Full article ">Figure 15
<p>Number of available links and PDOP values of BJFS in scenario 4.</p>
Full article ">
19 pages, 5444 KiB  
Article
Two-Dimensional Directions Determination for GNSS Spoofing Source Based on MEMS-Based Dual-GNSS/INS Integration
by Chengzhong Zhang, Dingjie Wang and Jie Wu
Remote Sens. 2024, 16(23), 4568; https://doi.org/10.3390/rs16234568 - 5 Dec 2024
Viewed by 789
Abstract
Satellite navigation spoofing is a major challenge in the field of satellite/inertial integrated navigation security. To effectively enhance the anti-spoofing capability of a low-cost GNSS/MEMS-SINS integrated navigation system, this paper proposes a method integrating a dual-antenna global navigation satellite system (GNSS) and a [...] Read more.
Satellite navigation spoofing is a major challenge in the field of satellite/inertial integrated navigation security. To effectively enhance the anti-spoofing capability of a low-cost GNSS/MEMS-SINS integrated navigation system, this paper proposes a method integrating a dual-antenna global navigation satellite system (GNSS) and a micro-inertial measurement unit (MIMU) to determine the two-dimensional directions of spoofing signal sources. The proposed method evaluates whether the single-difference carrier-phase measurements conform to the corresponding directions given in ephemeris files and employs micro-inertial navigation technology to determine the two-dimensional directions of the signal source. Based on a set of short-baseline dual-station measurements, the accuracy of the proposed method in determining the two-dimensional azimuths of satellites in synchronous orbits is verified, and the deviation from the real value is evaluated. The experimental results show that the proposed method can effectively identify the spoofed satellite signals while providing high-precision direction information at three different distances: 100 m, 10 km, and 36,000 km. The two-dimensional angle errors do not exceed 0.2 rad, 0.05 rad, and 0.01 rad, respectively. Full article
Show Figures

Figure 1

Figure 1
<p>Typical application scenario for the proposed method. The dotted line represents the trajectory of the unmanned vehicle during satellite signal deception.</p>
Full article ">Figure 2
<p>Integrated navigation model.</p>
Full article ">Figure 3
<p>Single-differenced carrier phase between two antennas on the baseline vector-source plane.</p>
Full article ">Figure 4
<p>Signal source position in the body coordinate system.</p>
Full article ">Figure 5
<p>Domains and standardised points of inter-station single-difference integer ambiguities on B<sub>1</sub> and B<sub>3.</sub> The square area represents the range of values for carrier phase ambiguities N<sub>1</sub> and N<sub>3</sub>. Due to the integer constraints for N<sub>1</sub> and N<sub>3</sub>, the feasible candidates for (N<sub>1</sub>, N<sub>3</sub>) can be illustrated as the dots in the square area.</p>
Full article ">Figure 6
<p>Experimental device.</p>
Full article ">Figure 7
<p>BDS real signal two-dimensional signal source determining errors. (<b>a</b>) Two-dimensional direction-finding error of BeiDou-2 in this experiment; (<b>b</b>) two-dimensional direction-finding error of BeiDou-3 in this experiment; (<b>c</b>) yaw angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line of BeiDou-2 in this experiment; (<b>d</b>) yaw angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line of BeiDou-3 in this experiment; (<b>e</b>) pitch angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line of BeiDou-2 in this experiment; and (<b>f</b>) pitch angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line of BeiDou-3 in this experiment.</p>
Full article ">Figure 7 Cont.
<p>BDS real signal two-dimensional signal source determining errors. (<b>a</b>) Two-dimensional direction-finding error of BeiDou-2 in this experiment; (<b>b</b>) two-dimensional direction-finding error of BeiDou-3 in this experiment; (<b>c</b>) yaw angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line of BeiDou-2 in this experiment; (<b>d</b>) yaw angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line of BeiDou-3 in this experiment; (<b>e</b>) pitch angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line of BeiDou-2 in this experiment; and (<b>f</b>) pitch angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line of BeiDou-3 in this experiment.</p>
Full article ">Figure 8
<p>Horizontal reference trajectories of the carrier.</p>
Full article ">Figure 9
<p>Near-range spoofing signal two-dimensional signal source determining errors. (<b>a</b>) Two-dimensional direction-finding error for the short-distance spoofing signal source; (<b>b</b>) yaw angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line for the short-distance spoofing signal source; and (<b>c</b>) pitch angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line for the short-distance spoofing signal source.</p>
Full article ">Figure 10
<p>MIMU attitude error curves.</p>
Full article ">Figure 11
<p>Simulation of spoofing signal two-dimensional signal source determining errors. (<b>a</b>) Simulation of two-dimensional direction-finding errors for the spoofing signal source; (<b>b</b>) simulation of the yaw angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line for the spoofing signal source; and (<b>c</b>) simulation of the pitch angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line for the spoofing signal source.</p>
Full article ">Figure 11 Cont.
<p>Simulation of spoofing signal two-dimensional signal source determining errors. (<b>a</b>) Simulation of two-dimensional direction-finding errors for the spoofing signal source; (<b>b</b>) simulation of the yaw angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line for the spoofing signal source; and (<b>c</b>) simulation of the pitch angle error and <math display="inline"><semantics> <mrow> <mn>3</mn> <mi>σ</mi> </mrow> </semantics></math> line for the spoofing signal source.</p>
Full article ">
34 pages, 14204 KiB  
Article
A Novel Algorithm for Precise Orbit Determination Using a Single Satellite Laser Ranging System Within a Single Arc for Space Surveillance and Tracking
by Dong-Gu Kim, Sang-Young Park and Eunji Lee
Aerospace 2024, 11(12), 989; https://doi.org/10.3390/aerospace11120989 - 29 Nov 2024
Viewed by 886
Abstract
A satellite laser ranging (SLR) system uses lasers to measure the range from ground stations to space objects with millimeter-level precision. Recent advances in SLR systems have increased their use in space surveillance and tracking (SST). The problem we are addressing, the precise [...] Read more.
A satellite laser ranging (SLR) system uses lasers to measure the range from ground stations to space objects with millimeter-level precision. Recent advances in SLR systems have increased their use in space surveillance and tracking (SST). The problem we are addressing, the precise orbit determination (POD) using one-dimensional range observations within a single arc, is challenging owing to infinite solutions because of limited observability. Therefore, general orbit determination algorithms struggle to achieve reasonable accuracy. The proposed algorithm redefines the cost value for orbit determination by leveraging residual tendencies in the POD process. The tendencies of residuals are quantified as R-squared values using Fourier series fitting to determine velocity vector information. The algorithm corrects velocity vector errors through the grid search method and least squares (LS) with a priori information. This approach corrects all six dimensions of the state vectors, comprising position and velocity vectors, utilizing only one dimension of the range observations. Simulations of three satellites using real data validate the algorithm. In all cases, the errors of the two-line element data (three-dimensional position error of 1 km and velocity error of 1 m/s, approximately) used as the initial values were reduced by tens of meters and the cm/s level, respectively. The algorithm outperformed the general POD algorithm using only the LS method, which does not effectively reduce errors. This study offers a more efficient and accurate orbit determination method, which improves the safety, cost efficiency, and effectiveness of space operations. Full article
(This article belongs to the Section Astronautics & Space Science)
Show Figures

Figure 1

Figure 1
<p>An example of a geometrically infinite number of orbits with range values observed at a single station within a single arc. (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi>y</mi> </mrow> <mrow> <mi>o</mi> <mo>,</mo> <msub> <mrow> <mi>t</mi> </mrow> <mrow> <mi>i</mi> </mrow> </msub> </mrow> </msub> </mrow> </semantics></math> is the actual observation value and <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>y</mi> </mrow> <mrow> <mi>c</mi> <msub> <mrow> <mo>,</mo> <mi>t</mi> </mrow> <mrow> <mi>i</mi> </mrow> </msub> </mrow> </msub> </mrow> </semantics></math> is the calculated observation value at observation time <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>t</mi> </mrow> <mrow> <mi>i</mi> </mrow> </msub> </mrow> </semantics></math>. <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="bold-italic">X</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msub> </mrow> </semantics></math> is the given initial value with error, and <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="bold-italic">X</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </mrow> </semantics></math> is the true state vector at <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>t</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> <mo>.</mo> </mrow> </semantics></math> <math display="inline"><semantics> <mrow> <msub> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">X</mi> </mrow> <mo>^</mo> </mover> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </mrow> </semantics></math> is the estimated state vector by POD using <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>y</mi> </mrow> <mrow> <mi>o</mi> <mo>,</mo> <msub> <mrow> <mi>t</mi> </mrow> <mrow> <mi>i</mi> </mrow> </msub> </mrow> </msub> </mrow> </semantics></math>).</p>
Full article ">Figure 2
<p>Simplified flowchart of RSR method.</p>
Full article ">Figure 3
<p>Example of unwanted overshoot effect: approximation of square wave by 30-term Fourier series [<a href="#B34-aerospace-11-00989" class="html-bibr">34</a>].</p>
Full article ">Figure 4
<p>Correlation of LS and position error (solid line) with dynamics (dotted line) and measurement error (dashed line).</p>
Full article ">Figure 5
<p>Examples of the bow-tie effect depending on the inadequacies of the LS dominated by the influence of the velocity vector error in the initial state.</p>
Full article ">Figure 6
<p>Meaning of SST and SSE for the R-squared value.</p>
Full article ">Figure 7
<p>Comparison of random search, grid search, and proposed grid search methods.</p>
Full article ">Figure 8
<p>Algorithm process of proposed grid search method for velocity vector searching.</p>
Full article ">Figure 9
<p>Flowchart of RSR POD algorithm.</p>
Full article ">Figure 10
<p>Range observations, which are the laser photon’s two-way travel time divided in half and multiplied by the speed of light, of Cryosat-2, Sentinel-3A, and Saral satellites utilized in the simulations.</p>
Full article ">Figure 11
<p>The comparison of errors between the SP-3 and TLE data in the GCRF for approximately 1 month around the time of observation in each case; (<b>a</b>) 3D position and velocity errors in all cases, (<b>b</b>) velocity vector errors from the Cryosat-2 satellite’s TLE near December 2018, (<b>c</b>) velocity vector errors from the Sentinel-3A satellite’s TLE near June 2018, (<b>d</b>) velocity vector errors from the Saral satellite’s TLE near August 2018.</p>
Full article ">Figure 12
<p>Results of the RSR POD for Case 1 (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mrow> <mi>v</mi> </mrow> </msub> </mrow> </semantics></math> is the velocity vector error computed from SP-3, assumed to be the true error value); (<b>a</b>) optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>x</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> among the gridded values <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>x</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> </mrow> </semantics></math> on the x-axis, (<b>b</b>) optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>y</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> among the gridded values <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>y</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> </mrow> </semantics></math> on the y-axis, (<b>c</b>) optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>z</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> among the gridded values <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>z</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> </mrow> </semantics></math> on the z-axis, (<b>d</b>) all cost values <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi>J</mi> </mrow> <mrow> <mi>R</mi> <mi>S</mi> <mi>R</mi> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> <mo>(</mo> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> </mrow> </semantics></math>) and minimum cost <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi>J</mi> </mrow> <mrow> <mi>R</mi> <mi>S</mi> <mi>R</mi> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> <mo>(</mo> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math>) in the k-th loop.</p>
Full article ">Figure 13
<p>Calculated residuals <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>19</mn> </mrow> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> with the first given initial state from TLE (left y-axis). <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>19</mn> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> using the optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> of the first, second, and third loops (right y-axis), the fitted curve <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="bold-italic">f</mi> </mrow> <mrow> <mn mathvariant="bold">6</mn> </mrow> </msub> <mo>(</mo> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>19</mn> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> <mo>)</mo> </mrow> </semantics></math> using the sixth Fourier series function for residuals <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>19</mn> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> in the RSR POD algorithm process for Case 1, and the difference <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>θ</mi> </mrow> <mrow> <mi>i</mi> </mrow> </msub> </mrow> </semantics></math> between the dependent variable of the fitted curve at observation time <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>t</mi> </mrow> <mrow> <mi>i</mi> </mrow> </msub> </mrow> </semantics></math> and the residual <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mi>i</mi> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> as data used in the fitting.</p>
Full article ">Figure 14
<p>POD results using only LS compared to that using RSR in the same condition for Case 1; (<b>a</b>) 3D position and velocity errors, (<b>b</b>) residuals (only LS: left y-axis, RSR: right y-axis).</p>
Full article ">Figure 15
<p>Results of the RSR POD for Case 2 (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mrow> <mi>v</mi> </mrow> </msub> </mrow> </semantics></math> is the velocity vector error computed from SP-3, assumed to be the true error value); (<b>a</b>) optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>x</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> among the gridded values <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>x</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> </mrow> </semantics></math> on the x-axis, (<b>b</b>) optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>y</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> among the gridded values <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>y</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> </mrow> </semantics></math> on the y-axis, (<b>c</b>) optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>z</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> among the gridded values <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>z</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> </mrow> </semantics></math> on the z-axis, (<b>d</b>) all cost values <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi>J</mi> </mrow> <mrow> <mi>R</mi> <mi>S</mi> <mi>R</mi> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> <mo>(</mo> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> </mrow> </semantics></math>) and minimum cost <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi>J</mi> </mrow> <mrow> <mi>R</mi> <mi>S</mi> <mi>R</mi> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> <mo>(</mo> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math>) in the first, third, and fifth loops.</p>
Full article ">Figure 16
<p>Calculated residuals <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>18</mn> </mrow> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> with the first given initial state from TLE (left y-axis). <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>18</mn> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> using the optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> of the first, third, and fifth loops (right y-axis), the fitted curve <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="bold-italic">f</mi> </mrow> <mrow> <mn mathvariant="bold">6</mn> </mrow> </msub> <mo>(</mo> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>18</mn> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> <mo>)</mo> </mrow> </semantics></math> using the sixth Fourier series function for residuals <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>18</mn> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> in the RSR POD algorithm process for Case 2, and the difference <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>θ</mi> </mrow> <mrow> <mi>i</mi> </mrow> </msub> </mrow> </semantics></math> between the dependent variable of the fitted curve at observation time <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>t</mi> </mrow> <mrow> <mi>i</mi> </mrow> </msub> </mrow> </semantics></math> and the residual <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mi>i</mi> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> as data used in the fitting.</p>
Full article ">Figure 17
<p>POD results of only LS compared to RSR in the same condition for Case 2; (<b>a</b>) 3D position and velocity errors, (<b>b</b>) residuals (only LS: left y-axis, RSR: right y-axis).</p>
Full article ">Figure 18
<p>Results of the RSR POD for Case 3 (<math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mrow> <mi>v</mi> </mrow> </msub> </mrow> </semantics></math> is the velocity vector error computed from SP-3, which is assumed to be the true error value); (<b>a</b>) optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>x</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> among the gridded values <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>x</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> </mrow> </semantics></math> on the x-axis, (<b>b</b>) optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>y</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> among the gridded values <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>y</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> </mrow> </semantics></math> on the y-axis, (<b>c</b>) optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>z</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> among the gridded values <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>v</mi> </mrow> <mrow> <mi>z</mi> </mrow> </msub> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> </mrow> </semantics></math> on the z-axis, (<b>d</b>) all cost values <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi>J</mi> </mrow> <mrow> <mi>R</mi> <mi>S</mi> <mi>R</mi> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> <mo>(</mo> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> </mrow> </semantics></math>) and minimum cost <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi>J</mi> </mrow> <mrow> <mi>R</mi> <mi>S</mi> <mi>R</mi> </mrow> <mrow> <mi>k</mi> </mrow> </msubsup> <mo>(</mo> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math>) in the first, second, third, fourth, and fifth loops.</p>
Full article ">Figure 19
<p>Calculated residuals <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>23</mn> </mrow> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> with the first given initial state from TLE (left y-axis). <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>23</mn> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> using the optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> of the first, third, and fifth loops (right y-axis), the fitted curve <math display="inline"><semantics> <mrow> <msub> <mrow> <mi mathvariant="bold-italic">f</mi> </mrow> <mrow> <mn mathvariant="bold">6</mn> </mrow> </msub> <mo>(</mo> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>23</mn> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> <mo>)</mo> </mrow> </semantics></math> using the sixth Fourier series function for residuals <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>23</mn> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> in the RSR POD algorithm process for Case 3, and the difference <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>θ</mi> </mrow> <mrow> <mi>i</mi> </mrow> </msub> </mrow> </semantics></math> between the dependent variable of the fitted curve at observation time <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>t</mi> </mrow> <mrow> <mi>i</mi> </mrow> </msub> </mrow> </semantics></math> and the residual <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi mathvariant="bold-italic">ϵ</mi> </mrow> <mrow> <mi>i</mi> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> as data used in the fitting.</p>
Full article ">Figure 20
<p>POD results of only LS compared to RSR in the same condition for Case 3; (<b>a</b>) 3D position and velocity errors, (<b>b</b>) residuals (only LS: left y-axis, RSR: right y-axis).</p>
Full article ">Figure 21
<p>Results compared to SP-3 by propagating the state vector <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">X</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>t</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <mn>0</mn> </mrow> </msubsup> </mrow> </semantics></math> given initially through the TLE, the estimated state vector <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">X</mi> </mrow> <mo>^</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>t</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <mn>0</mn> </mrow> </msubsup> </mrow> </semantics></math> through the only LS method, and the estimated state vector <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi mathvariant="bold-italic">X</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <msub> <mrow> <mi>t</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msub> </mrow> <mrow> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> calculated through the RSR POD algorithm for 24 h; (<b>a</b>) Case 1: Cryosat-2 satellite, (<b>b</b>) Case 2: Sentinel-3A satellite, (<b>c</b>) Case 3: Saral satellite.</p>
Full article ">Figure 22
<p>RMS of residuals per gridded vector of the first, second, and third loops for Case 1.</p>
Full article ">Figure 23
<p>Calculated residuals <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi>ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>11</mn> </mrow> <mrow> <msubsup> <mrow> <mi>e</mi> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> with the first given initial state from TLE (left y-axis). <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi>ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>11</mn> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> using the optimal solution <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>k</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </semantics></math> of the first, third, and fifth loops (right y-axis), the fitted curve <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>f</mi> </mrow> <mrow> <mn>4</mn> </mrow> </msub> <mo>(</mo> <msubsup> <mrow> <mi>ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>11</mn> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> <mo>)</mo> </mrow> </semantics></math> using the fourth Fourier series function for residuals <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi>ϵ</mi> </mrow> <mrow> <mn>1</mn> <mo>:</mo> <mn>11</mn> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> in the RSR POD algorithm process for 1/2 arc length of Case 3, and the difference <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>θ</mi> </mrow> <mrow> <mi>i</mi> </mrow> </msub> </mrow> </semantics></math> between the dependent variable of the fitted curve at observation time <math display="inline"><semantics> <mrow> <msub> <mrow> <mi>t</mi> </mrow> <mrow> <mi>i</mi> </mrow> </msub> </mrow> </semantics></math> and the residual <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi>ϵ</mi> </mrow> <mrow> <mi>i</mi> </mrow> <mrow> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mi>n</mi> <mo>*</mo> </mrow> </msubsup> </mrow> </msubsup> </mrow> </semantics></math> as data used in the fitting.</p>
Full article ">Figure 24
<p>All cost values, R-squared, <math display="inline"><semantics> <mrow> <msubsup> <mrow> <mi>J</mi> </mrow> <mrow> <mi>R</mi> <mi>S</mi> <mi>R</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msubsup> <mo>(</mo> <msubsup> <mrow> <mover accent="true"> <mrow> <mi>e</mi> </mrow> <mo>¯</mo> </mover> </mrow> <mrow> <mi>v</mi> </mrow> <mrow> <mn>1</mn> </mrow> </msubsup> </mrow> </semantics></math>), in the first loop of the RSR POD for 1/4 arc length of Case 3.</p>
Full article ">Figure 25
<p>Bow-tie effect following velocity error for MEO satellite, LAGEOS-1; (<b>a</b>) not-corrected velocity error case (error level: meter), (<b>b</b>) corrected velocity error case (error level: centimeter).</p>
Full article ">Figure 26
<p>Results of POD by EKF and UKF for Cases 1, 2, and 3; (<b>a</b>) 3D position and velocity errors for Case 1, (<b>b</b>) 3D position and velocity errors for Case 2, (<b>c</b>) 3D position and velocity errors for Case 3, and (<b>d</b>) residuals for three cases.</p>
Full article ">
18 pages, 6209 KiB  
Article
Impact of Latency and Continuity of GNSS Products on Filter-Based Real-Time LEO Satellite Clock Determination
by Meifang Wu, Kan Wang, Jinqian Wang, Wei Xie, Jiawei Liu, Beixi Chen, Yulong Ge, Ahmed El-Mowafy and Xuhai Yang
Remote Sens. 2024, 16(22), 4315; https://doi.org/10.3390/rs16224315 - 19 Nov 2024
Viewed by 758
Abstract
High-precision Low Earth Orbit (LEO) satellite clocks are essential for LEO-augmented Positioning, Navigation, and Timing (PNT) services. Nowadays, high-precision LEO satellite clocks can be determined in real-time using a Kalman filter either onboard or on the ground, as long as the GNSS observations [...] Read more.
High-precision Low Earth Orbit (LEO) satellite clocks are essential for LEO-augmented Positioning, Navigation, and Timing (PNT) services. Nowadays, high-precision LEO satellite clocks can be determined in real-time using a Kalman filter either onboard or on the ground, as long as the GNSS observations collected onboard LEO satellites can be transmitted to the ground in real-time. While various real-time and high-precision GNSS products are available nowadays in the latter case, their continuity and latencies in engineering reality are not as perfect as expected and will lead to unignorable impacts on the precision of the real-time LEO satellite clocks. In this study, based on real observations of Sentinel-3B, the impacts of different latencies and continuity of the real-time GNSS products on LEO real-time clocks are determined and discussed for two scenarios, namely the “epoch estimation” and “arc estimation” scenarios. The former case refers to the traditional filter-based processing epoch-by-epoch, and the latter case connects LEO satellite clocks from different rounds of filter-based processing under a certain arc length. The two scenarios lead to the “end-loss” and “mid-gap” situations. Latencies of the real-time GNSS products are discussed for the cases of orbit-only latency, clock-only latency, and combined forms, and different handling methods for the missing GNSS satellite clocks are discussed and compared. Results show that the real-time LEO satellite clock precision is very sensitive to the precision of real-time GNSS satellite clocks, and prediction of the latter becomes essential in case of their latencies. For the “end-loss” situation, with a latency of 30 to 120 s for the GNSS real-time clocks, the LEO satellite clock precision is reduced from about 0.2 to 0.28–0.57 ns. Waiting for the GNSS products in case of their short latencies and predicting the LEO satellite clocks instead could be a better option. For “arc-estimation”, when the gap of GNSS real-time products increases from 5 to 60 min, the real-time LEO clock precision decreases from 0.26 to 0.32 ns. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The continuity of the received real-time GPS satellite clocks from four analysis centers.</p>
Full article ">Figure 2
<p>Epoch-processing (<b>left</b>) and arc-processing (<b>right</b>) of real-time LEO satellite clock determination based on the kinematic model using the Kalman filter.</p>
Full article ">Figure 3
<p>Processing timelines of the two methods for real-time LEO satellite clocks determination based on the kinematic model and using the Kalman filter.</p>
Full article ">Figure 4
<p>(<b>Top</b>): Sentinel-3B satellite clocks determined using the Kalman filter-based kinematic model with the CNES real-time products without latency (blue line) and using the BLS-based reduced-dynamic model with the CODE final products (red line); (<b>Bottom</b>): differences between the two sets of clocks (real-time and final) in the top panel after unifying the time reference.</p>
Full article ">Figure 5
<p>Sentinel-3B clocks determined using the Kalman filter-based kinematic model and CNES real-time products with different latencies for Scenario A. The clocks using CNES real-time products without loss at the end are used as the reference.</p>
Full article ">Figure 6
<p>Sentinel-3B clocks determined using the Kalman filter-based kinematic model and CNES real-time products with different latencies for Scenario A. The BLS clocks using the CODE final products are used as the reference.</p>
Full article ">Figure 7
<p>Sentinel-3B clocks in the five scenarios (see <a href="#remotesensing-16-04315-t003" class="html-table">Table 3</a>) with a latency of 120 s compared with the clocks estimated using the CNES real-time products without latencies.</p>
Full article ">Figure 8
<p>Sentinel-3B clocks in the five scenarios (with a latency of 120 s, see <a href="#remotesensing-16-04315-t003" class="html-table">Table 3</a>) and the clocks without latencies in the CNES products, compared with the BLS clocks using the CODE final products.</p>
Full article ">Figure 9
<p>The MDEVs of the satellite clocks (<b>left</b>) and the clock estimation errors and prediction errors over 120 s (<b>right</b>) for G16 and G30.</p>
Full article ">Figure 10
<p>Sentinel-3B satellite clocks determined using the Kalman filter-based kinematic model with the CNES real-time products without latency (black line); using the BLS-based reduced-dynamic model with the CODE final products (red line); the connected near-real-time clocks estimated in a Kalman filter using the CNES real-time products without gaps (blue line).</p>
Full article ">Figure 11
<p>Differences between the C and K clocks (<b>top</b>) and C and R clocks (<b>bottom</b>) for Sentinel-3B. The explanations of C, K, and R clocks are given in <a href="#remotesensing-16-04315-f010" class="html-fig">Figure 10</a>.</p>
Full article ">Figure 12
<p>Near-real-time Sentinel-3B clocks based on a Kalman filter-based kinematic model, using CNES products with different “mid-gap” durations. (<b>Top</b>): Sentinel-3B clocks for the entire day; (<b>Bottom</b>): Sentinel-3B clocks from the first epoch affected to the last epoch affected.</p>
Full article ">Figure 13
<p>Near-real-time Sentinel-3B clock errors using the Kalman filter-based kinematic model with the CNES real-time products suffering from different “mid-gap” durations, compared with those using continuous CNES real-time products. (<b>Top</b>): Sentinel-3B clocks for the entire day; (<b>Bottom</b>): Sentinel-3B clocks from the first epoch affected to the last epoch affected.</p>
Full article ">Figure 14
<p>Near-real-time Sentinel-3B clock errors determined using the Kalman filter-based kinematic model with CNES products of different “mid-gap” durations, compared with the BLS clocks determined using the CODE final products. (<b>Top</b>): Sentinel-3B clocks for the entire day; (<b>Bottom</b>): Sentinel-3B clocks from the first epoch affected to the last epoch affected.</p>
Full article ">
17 pages, 5464 KiB  
Article
Geographically-Informed Modeling and Analysis of Platform Attitude Jitter in GF-7 Sub-Meter Stereo Mapping Satellite
by Haoran Xia, Xinming Tang, Fan Mo, Junfeng Xie and Xiang Li
ISPRS Int. J. Geo-Inf. 2024, 13(11), 413; https://doi.org/10.3390/ijgi13110413 - 15 Nov 2024
Viewed by 862
Abstract
The GF-7 satellite, China’s inaugural sub-meter-level stereoscopic mapping satellite, has been deployed for a wide range of applications, including natural resource investigation, environmental monitoring, fundamental surveying, and the development of global geospatial information resources. The satellite’s stable platform and reliable imaging systems are [...] Read more.
The GF-7 satellite, China’s inaugural sub-meter-level stereoscopic mapping satellite, has been deployed for a wide range of applications, including natural resource investigation, environmental monitoring, fundamental surveying, and the development of global geospatial information resources. The satellite’s stable platform and reliable imaging systems are crucial for achieving high-quality imaging and precise attitude measurements. However, the satellite’s operation is affected by both internal and external factors, which induce vibrations in the satellite platform, thereby affecting image quality and mapping accuracy. To address this challenge, this paper proposes a novel method for constructing a satellite platform vibration model based on geographic location information. The model is developed by integrating composite data from star sensors and gyroscopes (gyro) with subsatellite point location data. The experimental methodology involves the composite processing of gyro data and star sensor optical axis angles, integration of the processed data through time-matching and normalization, and denoising of the integrated data, followed by trigonometric fitting to capture the periodic characteristics of platform vibrations. The positions of the satellite substellar points are determined from the satellite orbit data. A rigorous geometric imaging model is then used to construct a vibration model with geographic location correlation in combination with the satellite subsatellite point positions. The experimental results demonstrate the following: (1) Over the same temporal range, there is a significant convergence in the waveform similarities between the gyro data and the star sensor optical axis angles, indicating a strong correlation in the jitter information; (2) The platform vibration exhibits a robust correlation with the satellite’s geographic location along its orbit. Specifically, the model reveals that the GF-7 satellite experiences the maximum vibration amplitude between 5° S and 20° S latitude during its ascending phase, and the minimum vibration amplitude between 5° N and 20° N latitude during the descending phase. The model established in this study offers theoretical support for optimizing satellite attitude and mitigating platform vibrations. Full article
Show Figures

Figure 1

Figure 1
<p>Comparative analysis of real and ideal vibration information from the star sensor and the gyro.</p>
Full article ">Figure 2
<p>Flow chart for vibration modeling method of the GF-7 satellite.</p>
Full article ">Figure 3
<p>Time consistency analysis of gyro data and star sensor data.</p>
Full article ">Figure 4
<p>Time offset analysis in composite data from the gyro and star sensor. Time Offset is the temporal difference between the two datasets. T1 is a specific time point.</p>
Full article ">Figure 5
<p>Schematic of the moving average filter.</p>
Full article ">Figure 6
<p>Satellite orbital data for calculating geographic locations of subsatellite points.</p>
Full article ">Figure 7
<p>Filtered gyro data across multiple satellite tracks: (<b>a</b>) Track 016396; (<b>b</b>) Track 016416; (<b>c</b>) Track 016431; (<b>d</b>) Track 016446.</p>
Full article ">Figure 8
<p>Optical axis clamping angle measurements from star sensor.</p>
Full article ">Figure 9
<p>Composite analysis of gyro and star sensor pinch angle data.</p>
Full article ">Figure 10
<p>Denoising results using moving average filter.</p>
Full article ">Figure 11
<p>Results of trigonometric function fitting for vibration data analysis. The blue color indicates composite data. The red curve represents the fitted result.</p>
Full article ">Figure 12
<p>Geographic distribution of satellite orbital paths.</p>
Full article ">Figure 13
<p>Variations in satellite flutter amplitude relative to geographic location.</p>
Full article ">
Back to TopTop