Ionospheric and receiver DCB-constrained multi-GNSS single-frequency PPP integrated with MEMS inertial measurements

Single-frequency precise point positioning (SF-PPP) is a potential precise positioning technique due to the advantages of the high accuracy in positioning after convergence and the low cost in operation. However, there are still challenges limiting its applications at present, such as the long convergence time, the low reliability, and the poor satellite availability and continuity in kinematic applications. In recent years, the achievements in the dual-frequency PPP have confirmed that its performance can be significantly enhanced by employing the slant ionospheric delay and receiver differential code bias (DCB) constraint model, and the multi-constellation Global Navigation Satellite Systems (GNSS) data. Accordingly, we introduce the slant ionospheric delay and receiver DCB constraint model, and the multi-GNSS data in SF-PPP modular together. In order to further overcome the drawbacks of SF-PPP in terms of reliability, continuity, and accuracy in the signal easily blocking environments, the inertial measurements are also adopted in this paper. Finally, we form a new approach to tightly integrate the multi-GNSS single-frequency observations and inertial measurements together to ameliorate the performance of the ionospheric delay and receiver DCB-constrained SF-PPP. In such model, the inter-system bias between each two GNSS systems, the inter-frequency bias between each two GLONASS frequencies, the hardware errors of the inertial sensors, the slant ionospheric delays of each user-satellite pair, and the receiver DCB are estimated together with other parameters in a unique Kalman filter. To demonstrate its performance, the multi-GNSS and low-cost inertial data from a land-borne experiment are analyzed. The results indicate that visible positioning improvements in terms of accuracy, continuity, and reliability can be achieved in both open-sky and complex conditions while using the proposed model in this study compared to the conventional GPS SF-PPP.     

A two-step ionospheric modeling algorithm considering the impact of GLONASS pseudo-range inter-channel biases

The Global Navigation Satellite System presents a plausible and cost-effective way of computing the total electron content (TEC). But TEC estimated value could be seriously affected by the differential code biases (DCB) of frequency-dependent satellites and receivers. Unlike GPS and other satellite systems, GLONASS adopts a frequency-division multiplexing access mode to distinguish different satellites. This strategy leads to different wavelengths and inter-frequency biases (IFBs) for both pseudo-range and carrier phase observations, whose impacts are rarely considered in ionospheric modeling. We obtained observations from four groups of co-stations to analyze the characteristics of the GLONASS receiver P1P2 pseudo-range IFB with a double-difference method. The results showed that the GLONASS P1P2 pseudo-range IFB remained stable for a period of time and could catch up to several meters, which cannot be absorbed by the receiver DCB during ionospheric modeling. Given the characteristics of the GLONASS P1P2 pseudo-range IFB, we proposed a two-step ionosphere modeling method with the priori IFB information. The experimental analysis showed that the new algorithm can effectively eliminate the adverse effects on ionospheric model and hardware delay parameters estimation in different space environments. During high solar activity period, compared to the traditional GPS + GLONASS modeling algorithm, the absolute average deviation of TEC decreased from 2.17 to 2.07 TECu (TEC unit); simultaneously, the average RMS of GPS satellite DCB decreased from 0.225 to 0.219 ns, and the average deviation of GLONASS satellite DCB decreased from 0.253 to 0.113 ns with a great improvement in over 55%.     

Tropospheric refractivity and zenith path delays from least-squares collocation of meteorological and GNSS data

Precise positioning requires an accurate a priori troposphere model to enhance the solution quality. Several empirical models are available, but they may not properly characterize the state of troposphere, especially in severe weather conditions. Another possible solution is to use regional troposphere models based on real-time or near-real time measurements. In this study, we present the total refractivity and zenith total delay (ZTD) models based on a numerical weather prediction (NWP) model, Global Navigation Satellite System (GNSS) data and ground-based meteorological observations. We reconstruct the total refractivity profiles over the western part of Switzerland and the total refractivity profiles as well as ZTDs over Poland using the least-squares collocation software COMEDIE (Collocation of Meteorological Data for Interpretation and Estimation of Tropospheric Pathdelays) developed at ETH Zürich. In these two case studies, profiles of the total refractivity and ZTDs are calculated from different data sets. For Switzerland, the data set with the best agreement with the reference radiosonde (RS) measurements is the combination of ground-based meteorological observations and GNSS ZTDs. Introducing the horizontal gradients does not improve the vertical interpolation, and results in slightly larger biases and standard deviations. For Poland, the data set based on meteorological parameters from the NWP Weather Research and Forecasting (WRF) model and from a combination of the NWP model and GNSS ZTDs shows the best agreement with the reference RS data. In terms of ZTD, the combined NWP-GNSS observations and GNSS-only data set exhibit the best accuracy with an average bias (from all stations) of 3.7 mm and average standard deviations of 17.0 mm w.r.t. the reference GNSS stations.     

Group delay variations of GPS transmitting and receiving antennas

GPS code pseudorange measurements exhibit group delay variations at the transmitting and the receiving antenna. We calibrated C1 and P2 delay variations with respect to dual-frequency carrier phase observations and obtained nadir-dependent corrections for 32 satellites of the GPS constellation in early 2015 as well as elevation-dependent corrections for 13 receiving antenna models. The combined delay variations reach up to 1.0 m (3.3 ns) in the ionosphere-free linear combination for specific pairs of satellite and receiving antennas. Applying these corrections to the code measurements improves code/carrier single-frequency precise point positioning, ambiguity fixing based on the Melbourne–Wübbena linear combination, and determination of ionospheric total electron content. It also affects fractional cycle biases and differential code biases.     

Error analysis of the 3D similarity coordinate transformation

The 3D similarity coordinate transformation with the Gauss–Helmert error model is investigated. The first-order error analysis of an analytical least-squares solution to this problem is developed in detail. While additive errors are assumed in the translation and scale estimates, a 3 × 1 multiplicative error vector is defined to effectively parameterize the rotation matrix estimation error. The propagation of the errors in the coordinate measurements to the errors in the estimated transformation parameters is derived step-by-step, and the formulae for calculating the variance–covariance matrix of the estimated parameters are presented.     

A Novel Approach of Despeckling SAR Images Using Nonlocal Means Filtering

Synthetic aperture radar (SAR) is a day and night, all weather satellite imaging technology. Inherent property of SAR image is speckle noise which produces granular patterns in the image. Speckle noise occurs due to the interference of backscattered echo from earth’s rough surface. There are various speckle reduction techniques in spatial domain and transform domain. Non local means filtering (NLMF) is the technique used for denoising which uses Gaussian weights. In NLMF algorithm, the filtering is performed by taking the weighted mean of all the pixels in a selected search area. The weight given to the pixel is based on the similarity measure calculated as the weighted Euclidean distance over the two windows. Non local means filtering smoothes out homogeneous areas but edges are not preserved. So a discontinuity adaptive weight is used in order to preserve heterogeneous areas like edges. This technique is called as discontinuity adaptive non local means filtering and is well-adapted and robust in the case of Additive White Gaussian Noise (AWGN) model. But speckle is a multiplicative random noise and hence Euclidean distance is not a good choice. This paper presents evaluation results of using different distance measures for improving the accuracy of the Non local means filtering technique. The results are verified using real and synthetic images and from the results it can be concluded that the usage of Manhattan distance improves the accuracy of NLMF technique. Non local approach is used as a preprocessing or post processing technique for many denoising algorithms. So improving NLMF technique would help improving many of the existing denoising techniques.     

Proximal Remote Sensing of Herbicide and Drought Stress in Field Grown Colocasia and Sweet Potato Plants by Sunlight-Induced Chlorophyll Fluorescence Imaging

Chlorophyll fluorescence is an indicator of plant photosynthetic activity and has been used to monitor the health status of vegetation. Several studies have exploited the application of red/far-red chlorophyll fluorescence ratio in detecting the impact of various types of stresses in plants. Recently, sunlight-induced chlorophyll fluorescence imaging has been used to detect and discriminate different stages of mosaic virus infection in potted cassava plants with a multi-spectral imaging system (MSIS). In this study, the MSIS is used to investigate the impact of drought and herbicide stress in field grown crop plants. Towards this control and treatment groups of colocasia and sweet potato plants were grown in laterite soil beds and the reflectance images of these crop plants were recorded up to 14-days of treatment at the Fraunhofer lines of O₂ B at 687 nm and O₂ A at 759.5 nm and the off-lines at 684 and 757.5 nm. The recorded images were analyzed using the Fraunhofer Line Discrimination technique to extract the sunlight-induced chlorophyll fluorescence component from the reflectance images of the plant leaves. As compared to the control group, the chlorophyll fluorescence image ratio (F ₆₈₇/F ₇₆₀) in the treatment groups of both the plant varieties shows an increasing trend with increase in the extent of stress. Further, the F ₆₈₇/F ₇₆₀ ratio was found to correlate with the net photosynthetic rate (Pn) and stomatal conductance (g_s) of leaves. The correlation coefficient (R ²) for the relationship of F ₆₈₇/F ₇₆₀ ratio with Pn were found to be 0.78, 0.79 and 0.78, respectively for the control, herbicide treated and drought treated colocasia plants, while these were 0.77, 0.86 and 0.88, respectively for sweet potato plants. The results presented show the potential of proximal remote sensing and the application F ₆₈₇/F ₇₆₀ fluorescence image ratio for effective monitoring of stress-induced changes in field grown plants.     

基于交互多模型的水下目标跟踪方法

According to the requirements of real-time performance and reliability in underwater maneuvering target tracking as well as clarifying motion features of the underwater target, an interacting multiple model algorithm based on fuzzy logic inference （FIMM） is proposed. Maneuvering patterns of the target are represented by model sets, including the constant velocity model （CA）, the Singer mode~, and the nearly constant speed horizontal-turn model （HT） in FIMM technology. The simulation results show that compared to conventional IMM, the reliability and real-time performance of underwater target tracking can be improved by FIMM algorithm.     

Bringing GPS into harsh environments for fully automated deformation monitoring

Engineering projects that require deformation monitoring frequently utilize geodetic sensors to measure displacements of target points located in the deformation zone. In situations where control stations and targets are separated by a kilometer or more, GPS can offer higher precision position updates at more frequent intervals than can normally be achieved using total station technology. For large-scale deformation projects requiring the highest precision, it is therefore advisable to use a combination of the two sensors. In response to the need for high precision, continuous GPS position updates in harsh deformation monitoring environments, a software has been developed that employs triple-differenced carrier-phase measurements in a delayed-state Kalman filter. Two data sets were analyzed to test the capabilities of the software. In the first test, a GPS antenna was displaced using a translation stage to mimic slow deformation. In the second test, data collected at a large open pit mine were processed. It was shown that the delayed-state Kalman filter developed could detect millimeter-level displacements of a GPS antenna. The actual precision attained depends upon the amount of process noise infused at each epoch to accommodate the antenna displacements. Higher process noise values result in quicker detection times, but at the same time increase the noise in the solutions. A slow, 25 mm displacement was detected within 30 min of the full displacement with sigma values in E, N and U of ±10 mm or better. The same displacement could also be detected in less than 5 h with sigma values in E, N and U of ±5 mm or better. The software works best for detecting long period deformations (e.g., 20 mm per day or less) for which sigma values of 1–2 mm are attained in all three solution components. It was also shown that the triple-differenced carrier-phase observation can be used to significantly reduce the effects of residual tropospheric delay that would normally plague double-differenced observations in harsh GPS environments.

Geodetic intersection on the ellipsoid

Through each of two known points on the ellipsoid a geodesic is passing in a known azimuth. We solve the problem of intersection of the two geodesics. The solution for the latitude is obtained as a closed formula for the sphere plus a small correction, of the order of the eccentricity of the ellipsoid, which is determined by numerical integration. The solution is iterative. Once the latitude is obtained, the longitude is determined without iteration.