奇异值分解法在病态问题中的应用

用截断奇异值分解法和修正奇异值分解法对大地测量病态问题进行了处理,并与最小二乘的结果进行了比较,最后将两者方法进行结合同样对病态方程进行了处理,得到了结合奇异值分解的解,并与截断奇异值、修正奇异值分解的解以及真值进行比较,发现结合奇异值分解的解即修正奇异值截断法比截断奇异值和修正奇异值的解更加靠近于真值,修正奇异值截断法相比于截断奇异值和修正奇异值法对于病态方程抗干扰能力更强,更具有实际意义。     

利用多光谱卫星遥感和深度学习方法进行青藏高原积雪判识

青藏高原积雪对全球气候变化十分重要,针对已有积雪遥感判识方法中普遍采用的可见光与红外光谱数据易受复杂地形与高海拔影响,导致青藏高原地区积雪判识精度较低的问题,提出了一种基于多光谱遥感与地理信息数据特征级融合的积雪遥感判识方法:以风云三号卫星可见光与红外多光谱遥感资料与多要素地理信息作为数据源,由地面实测雪深数据与现有积雪产品交叉筛选出样本标签,构建并训练基于层叠去噪自编码器(SDAE)的特征融合与分类网络,从而有效辨识青藏高原遥感图像中的云、积雪以及无雪地表。经地面实测雪深数据验证,该方法分类精度显著高于使用相同数据源的FY-3A/MULSS积雪产品,略高于国际主流积雪产品MOD10A1与MYD10A1,并且年均云覆盖率最低。试验结果表明该方法可有效地减少云层对积雪判识的干扰,提升分类精度。     

利用压密方程计算分析西安市地面沉降特征

利用西安市地下水位监测资料,基于水文地质三维结构模型,在不考虑粘土层滞后压缩变形的理想情况下,采用压密方程计算获得西安市抽取地下水可引起的理论地面沉降量及沉降分布特征,结果表明:计算所得沉降显著区位于西安市西南部的高新区及东南部的曲江新区,沉降量呈由北向南递增的特征,与InSAR监测结果整体趋势具有较好的一致性;地面沉降分布特征受到活动地裂缝影响,沉降曲线呈近NE向偏转展布。研究结果可为沉降灾害预防研究提供宏观的参考信息。     

改进的距离约束最小二乘模糊度搜索算法

针对短基线最小二乘模糊度搜索算法搜索效率低的缺点,该文提出了一种改善模糊度搜索空间以提高模糊度搜索效率的方法。最小二乘搜索算法通过基线长度范围确定模糊度搜索空间,完成模糊度的搜索和确认。距离约束的短基线模糊度搜索空间是一个空心椭球,各模糊度之间具有相关性,影响模糊度搜索效率。该文通过最小二乘去相关调整方法降低模糊度之间的相关性,减少搜索次数。实验证明,该算法的搜索效率提高了25%以上,且随着基线长度的增加,搜索效率有更多的提升。     

Long-term soil moisture dynamics derived from GNSS interferometric reflectometry: a case study for Sutherland,South Africa

Soil moisture is a geophysical key observable for predicting floods and droughts, modeling weather and climate and optimizing agricultural management. Currently available in situ observations are limited to small sampling volumes and restricted number of sites, whereas measurements from satellites lack spatial resolution. Global navigation satellite system (GNSS) receivers can be used to estimate soil moisture time series at an intermediate scale of about 1000 m². In this study, GNSS signal-to-noise ratio (SNR) data at the station Sutherland, South Africa, are used to estimate soil moisture variations during 2008–2014. The results capture the wetting and drying cycles in response to rainfall. The GNSS Volumetric Water Content (VWC) is highly correlated (r ² = 0.8) with in situ observations by time-domain reflectometry sensors and is accurate to 0.05 m³/m³. The soil moisture estimates derived from the SNR of the L1 and L2P signals compared to the L2C show small differences with a RMSE of 0.03 m³/m³. A reduction in the SNR sampling rate from 1 to 30 s has very little impact on the accuracy of the soil moisture estimates (RMSE of the VWC difference 1–30 s is 0.01 m³/m³). The results show that the existing data of the global tracking network with continuous observations of the L1 and L2P signals with a 30-s sampling rate over the last two decades can provide valuable complementary soil moisture observations worldwide.     

Using a regional numerical weather prediction model for GNSS positioning over Brazil

The global navigation satellite system (GNSS) can provide centimeter positioning accuracy at low costs. However, in order to obtain the desired high accuracy, it is necessary to use high-quality atmospheric models. We focus on the troposphere, which is an important topic of research in Brazil where the tropospheric characteristics are unique, both spatially and temporally. There are dry regions, which lie mainly in the central part of the country. However, the most interesting area for the investigation of tropospheric models is the wet region which is located in the Amazon forest. This region substantially affects the variability of humidity over other regions of Brazil. It provides a large quantity of water vapor through the humidity convergence zone, especially for the southeast region. The interconnection and large fluxes of water vapor can generate serious deficiencies in tropospheric modeling. The CPTEC/INPE (Center for Weather Forecasting and Climate Studies/Brazilian Institute for Space Research) has been providing since July 2012 a numerical weather prediction (NWP) model for South America, known as Eta. It has yield excellent results in weather prediction but has not been used in GNSS positioning. This NWP model was evaluated in precise point positioning (PPP) and network-based positioning. Concerning PPP, the best positioning results were obtained for the station SAGA, located in Amazon region. Using the NWP model, the 3D RMS are less than 10 cm for all 24 h of data, whereas the values reach approximately 60 cm for the Hopfield model. For network-based positioning, the best results were obtained mainly when the tropospheric characteristics are critical, in which case an improvement of up to 7.2 % was obtained in 3D RMS using NWP models.     

Generating GPS satellite fractional cycle bias for ambiguity-fixed precise point positioning

With the development of precise point positioning (PPP), the School of Geodesy and Geomatics (SGG) at Wuhan University is now routinely producing GPS satellite fractional cycle bias (FCB) products with open access for worldwide PPP users to conduct ambiguity-fixed PPP solution. We provide a brief theoretical background of PPP and present the strategies and models to compute the FCB products. The practical realization of the two-step (wide-lane and narrow-lane) FCB estimation scheme is described in detail. With GPS measurements taken in various situations, i.e., static, dynamic, and on low earth orbit (LEO) satellites, the quality of FCB estimation and the effectiveness of PPP ambiguity resolution (AR) are evaluated. The comparison with CNES FCBs indicated that our FCBs had a good consistency with the CNES ones. For wide-lane FCB, almost all the differences of the two products were within ±0.05 cycles. For narrow-lane FCB, 87.8 % of the differences were located between ±0.05 cycles, and 97.4 % of them were located between ±0.075 cycles. The experimental results showed that, compared with conventional ambiguity-float PPP, the averaged position RMS of static PPP can be improved from (3.6, 1.4, 3.6) to (2.0, 1.0, 2.7) centimeters for ambiguity-fixed PPP. The average accuracy improvement in the east, north, and up components reached 44.4, 28.6, and 25.0 %, respectively. A kinematic, ambiguity-fixed PPP test with observation of 80 min achieved a position accuracy of better than 5 cm at the one-sigma level in all three coordinate components. Compared with the results of ambiguity-float, kinematic PPP, the positioning biases of ambiguity-fixed PPP were improved by about 78.2, 20.8, and 65.1 % in east, north, and up. The RMS of LEO PPP test was improved by about 23.0, 37.0, and 43.0 % for GRACE-A and GRACE-B in radial, tangential, and normal directions when AR was applied to the same data set. These results demonstrated that the SGG FCB products can be produced with high quality for users anywhere around the world to carry out ambiguity-fixed PPP solutions.     

Modified compressive sensing approach for GNSS signal reception in the presence of interference

Pre-processing traditional navigation signals in global navigation satellite system (GNSS) receivers includes the conversion of an analog-to-digital sample and acquisition following the basic principle of Nyquist sampling theory. This condition inevitably increases the system computation time and cost of a modern wideband receiver. In recent years, the compressive sensing (CS) approach has been proven to effectively reduce the number of measurement samples required for digital signal acquisition systems. This method gives new potential to this modern design. In this study, a modified compressive sensing algorithm for the acquisition of a GNSS signal that is contaminated by an interfering signal is presented. The proposed method attempts to combine CS demodulation and the subspace projecting method to enhance GNSS signal acquisition performance with interference present. First, the received signal is sub-sampled and aliased from a compressive sampling process. This operation maintains the restricted isometry property (RIP) condition of the second sampling process using a Toeplitz-structured sensing matrix, which replaces a conventional random sensing matrix. The matrix ensures that distances between desired signals on the set of sparse space are not influenced by the sampling process. Next, the interference is eliminated through the orthogonal feature between the interference signal and the desired signal using the subspace projecting method. This also preserves the RIP of the projecting matrix to ensure that the original structure of the linear projection of the signal is preserved. After this, an iterative least-square method is utilized to recover the correlator output from the reception samples taken by the CS demodulator. In addition, the signal detection performance in the presence of co-channel interference using a CS demodulator is analyzed and evaluated. Finally, the relationships between signal detection probability, compressive factor and signal bandwidth are also illustrated. Several numerical results are presented to verify the theory.     

Weak GPS signal tracking using FFT discriminator in open loop receiver

An open loop tracking architecture, which tracks GPS signals under weak and challenging conditions, is analyzed. The in-phase and quadrature-phase integration pair is regarded as a single tone complex signal. An FFT-based method is used as a frequency discriminator to estimate the Doppler frequency residual of the single tone signal. Another FFT-based method applies complex squaring to eliminate the effect of the navigation data bits polarities. The performance of the FFT-based discriminators is assessed in three criteria. Those criteria are the signal strength and dynamic range that can be tracked and the accuracy of the estimated Doppler frequency. In addition, the performance of the discriminators is analyzed to provide the theoretical and simulated peak detection probability. The results indicate that the FFT discriminator can track signals about 5 dB weaker than the signals that can be tracked by the complex squared FFT discriminator. In a quasi-static environment, the Doppler frequency residual can be assumed to be around zero, which can enable the FFT-based discriminators to track signals with approximately 2 dB less power. Moreover, the performance of the FFT-based discriminators is compared with the performance of two other frequency discriminators, namely the fast–slow and power-based. The comparison results indicate that these two frequency discriminators give higher frequency estimation accuracy, but they have a narrower dynamic range.     

Precise orbit determination of BeiDou constellation: method comparison

Chinese BeiDou navigation satellite system is in official service as a regional constellation with five geostationary earth orbit (GEO) satellites, five inclined geosynchronous satellite orbit (IGSO) satellites and four medium earth orbit (MEO) satellites. There are mainly two methods for precise orbit determination of the BeiDou constellation found in the current literatures. One is the independent single-system method, where only BeiDou observations are used without help from other GNSS systems. The other is the two-step GPS-assisted method where in the first step, GPS data are used to resolve some common parameters, such as station coordinates, receiver clocks and zenith tropospheric delay parameters, which are then introduced as known quantities in BeiDou processing in the second step. We conduct a thorough performance comparison between the two methods. Observations from the BeiDou experimental tracking stations and the IGS Multi-GNSS Experiment network from January 1 to March 31, 2013, are processed with the Positioning and Navigation Data Analyst (PANDA) software. The results show that for BeiDou IGSO and MEO satellites, the two-step GPS-assisted method outperforms the independent single-system method in both internal orbit overlap precision and external satellite laser ranging validation. For BeiDou GEO satellites, the two methods show close performances. Zenith tropospheric delays estimated from the first method are very close to those estimated from GPS precise point positioning in the second method, with differences of several millimeters. Satellite clock estimates from the two methods show similar performances when assessing the stability of the BeiDou on board clocks.