On the sensitivity of the location, size and shape of the GPS ambiguity search space to certain changes in the stochastic model

In this contribution we analyse in a qualitative sense for the geometry-free model the dependency of the location, the size and the shape of the ambiguity search space on different factors of the stochastic model. For this purpose a rather general stochastic model is used. It includes time-correlation, cross-correlation, satellite elevation dependency and the use of an a priori weighted ionospheric model, having the ionosphere-fixed model and the ionosphere-float model as special cases. It is shown that the location is invariant for changes in the cofactor matrix of the phase observables. This also holds true for the cofactor matrix of the code observables in the ionosphere-float case. As for time-correlation and satellite elevation dependency, it is shown that they only affect the size of the search space, but not its shape and orientation. It is also shown that the least-squares ambiguities, their variance matrix and its determinant, for, respectively, the ionosphere-fixed model, the ionosphere-float model and the ionosphere-weighted model, are all related through the same scalar weighted mean, the weight of which is governed by the variance ratio of the ionospheric delays and the code observables. A closed-form expression is given for the area of the search space in which all contributing factors are easily recognized. From it one can infer by how much the area gets blown up when the ionospheric spatial decorrelation increases. This multiplication factor is largest when one switches from the ionosphere-fixed model to the ionosphere-float model, in which case it is approximately equal to the ratio of the standard deviation of phase with that of code. The area gives an indication of the number of grid points inside the search space. Received: 11 November 1996 / Accepted: 21 March 1997     

An integrated wavelet concept of physical geodesy

For the determination of the earth's gravity field many types of observations are nowadays available, including terrestrial gravimetry, airborne gravimetry, satellite-to-satellite tracking, satellite gradio-metry, etc. The mathematical connection between these observables on the one hand and gravity field and shape of the earth on the other is called the integrated concept of physical geodesy. In this paper harmonic wavelets are introduced by which the gravitational part of the gravity field can be approximated progressively better and better, reflecting an increasing flow of observations. An integrated concept of physical geodesy in terms of harmonic wavelets is presented. Essential tools for approximation are integration formulas relating an integral over an internal sphere to suitable linear combinations of observation functionals, i.e. linear functionals representing the geodetic observables. A scale discrete version of multiresolution is described for approximating the gravitational potential outside and on the earth's surface. Furthermore, an exact fully discrete wavelet approximation is developed for the case of band-limited wavelets. A method for combined global outer harmonic and local harmonic wavelet modelling is proposed corresponding to realistic earth's models. As examples, the role of wavelets is discussed for the classical Stokes problem, the oblique derivative problem, satellite-to-satellite tracking, satellite gravity gradiometry and combined satellite-to-satellite tracking and gradiometry. Received: 28 February 1997 / Accepted: 17 November 1997     

简便实用的GPS网平差模型

采用数值微商求偏导数的方法，从而使ＧＰＳ网平差的计算模型非常简洁。选用平面坐标和大地高作为平差参数，可以方便地简化为平面网或高程网，加入地面归心数据和常规观测数据也变得十分容易。对旋转角大的地方独立坐标系，提出了处理办法。经对实测数据的处理，证明了模型的正确性     

Generalized inverses of nonlinear mappings and the nonlinear geodetic datum problem

Motivated by the existing theory of the geometric characteristics of linear generalized inverses of linear mappings, an attempt is made to establish a corresponding mathematical theory for nonlinear generalized inverses of nonlinear mappings in finite- dimensional spaces. The theory relies on the concept of fiberings consisting of disjoint manifolds (fibers) in which the domain and range spaces of the mappings are partitioned. Fiberings replace the quotient spaces generated by some characteristic subspaces in the linear case. In addition to the simple generalized inverse, the minimum-distance and the x ₀-nearest generalized inverse are introduced and characterized, in analogy with the least-squares and the minimum-norm generalized inverses of the linear case. The theory is specialized to the geodetic mapping from network coordinates to observables and the nonlinear transformations (Baarda's S-transformations) between different solutions are defined with the help of transformation parameters obtained from the solution of nonlinear equations. In particular, the transformations from any solution to an x ₀-nearest solution (corresponding to Meissl's inner solution) are given for two- and three-dimensional networks for both the similarity and the rigid transformation case. Finally the nonlinear theory is specialized to the linear case with the help of the singular-value decomposition and algebraic expressions with specific geometric meaning are given for all possible types of generalized inverses. Received: 11 April 1996 / Accepted: 19 April 1997     

Evaluation of three ionospheric delay computation methods for ground-based GNSS receivers

GNSS observables for ionospheric estimation are commonly based on carrier-to-code leveling (CCL) and precise point positioning (PPP) methods. The CCL method is a geometry-free method which uses carrier phase to level pseudorange observation for decreasing multipath error and observation noise. However, the ionospheric observable based on the CCL has been proven to be affected by leveling errors. The leveling errors are caused by pseudorange multipath and intraday variation of receiver DCB. To obtain more accurate ionospheric observable, the PPP method takes advantage of precise satellite-to-ground range for retrieving slant total electron content and is less affected by the leveling errors. Previous studies have only proven that the ionospheric observables extracted by the two methods are affected by the leveling errors. The influence on ionospheric observable by the pseudorange inter-receiver satellite bias (IRSB) of the receiver has not been taken into consideration. Also, the magnitude of the differences between the ionospheric observables extracted by the two methods has also not been given. In this work, three methods, namely, the CCL, the conventional ionospheric-free PPP method which uses the ionospheric-free Hatch–Melbourne–Wubbena (HMW) function, and the University of Calgary (UOFC) PPP method, are selected to analyze and compare the differences of ionospheric observables and the global ionospheric maps, using a large number of measured data from international GNSS service global stations. Experimental results show that the accuracy of ionospheric observables obtained by the three methods is not only related to the leveling error, but also pseudorange IRSB. The IRSB of the receiver exerts a major effect on the ionospheric observables obtained by the CCL method and a minor effect on the ionospheric observables obtained by the HMW and UOFC methods. The accuracies in the latter case are similar and superior to those obtained by the CCL. The differences of the ionospheric observables obtained by the CCL and UOFC methods, or the CCL and HMW methods, are at decimeter level, whereas the difference of the ionospheric observables obtained by the UOFC and HMW methods is at centimeter level. The UOFC method presented the highest single-frequency pseudorange positioning accuracy using estimated global ionospheric products, followed by the HMW and the CCL methods which presented the lowest positioning accuracy.     

GPS三频非差观测数据周跳的自动探测与改正研究

在GPS三频非差观测数据的处理中,由于伪距噪声的影响,利用原始的伪距和载波相位观测数据估计的模糊度误差比较大,不能用于探测和改正周跳。对原始观测数据进行平滑或适当的组合处理,可降低观测噪声的影响。因此,本文选择合适的经过平滑或组合处理后的观测数据作为探测周跳的检验量,探测并改正单个频点上的周跳。在分析了一般周跳的特点并在研究双频周跳自动探测与改正方法基础上,提出了选取检验量的四条基本原则。最后,依此原则选取了三个观测值组合作为周跳检验量,利用该组检验量实现三频非差观测数据周跳的自动探测与改正。     

An Information fusion approach for PALSAR data to retrieve soil moisture

Estimation of vegetation covered soil moisture with satellite images is still a challenging task. Several models are available for soil moisture retrieval in which water cloud model (WCM) is most common. But, it requires an estimation of accurate vegetation parameterization. Thus, there is a need to develop such an approach for soil moisture retrieval which minimize these limitations. Therefore, this paper deals with the soil moisture retrieval using fully polarimetric SAR data by fusing the information from different bands. Various polarimetric indices and observables were critically analysed, and found that the index; SPAN (total scattered power) gives better information of vegetation cover as compared to other indices/observables. Based on this, WCM model has been modified using SPAN as parameter and soil moisture content were retrieved.     

ISO2002: an analytical stochastic model of multi-difference GPS carrier-phase data

The differencing technique is useful in global positioning system (GPS) positioning when two or more GPS receivers collect simultaneous observables from common satellites at each epoch, and all carrier-phase observables have the same normal distribution. An analytical probability distribution of the single-, double-, triple- and multi-difference GPS observables is obtained. This analytical model, called ISO2002, has a good matrix structure, in which I indicates the number of receivers, S indicates the number of observed satellites, and O indicates the number of epochs. The variance–covariance matrix can be expressed as the Kronecker product of several small matrices, so its inverse is equal to the Kronecker product of the inverses of these sub-matrices. Moreover, these small matrices are circulant or symmetric diagonal Toeplitz matrices, so their inverses have analytical solutions. The analytical model ISO2002 proposed to compute the inverse variance–covariance matrix is shown to be very effective.     

基于三频SNR数据探测北斗GEO卫星多径效应

多径误差是全球卫星导航系统（global navigation satellite system,GNSS）高精度实时定位的主要误差源。由于多径效应具有时变特性且与观测环境有关,难以有效探测与剔除,尤其是动态定位领域。鉴于多径效应对于同颗卫星不同频率观测值影响不同,本文采用三频信噪比（signal-noise ratio,SNR）频率间差分数据,设计算法探测北斗对地静止卫星（geostationary orbit satellite,GEO）多径误差,并依据开阔环境SNR频间差分数据序列进行统计分析,确定其探测阈值。实验结果表明,北斗GEO卫星易受多径误差影响,传统方法难以实时有效探测,本文提出的新方法可有效探测GEO多径误差,并基于设定的阈值对其数据进行降权或剔除,可应用于GNSS数据质量控制。     

GPS相位观测值的模糊度函数法与最小二乘法的等价性及模糊度函数法的改进

本文主要讨论了ＧＰＳ相位观测工函数法处理与最小二乘法处理（如双闪差分相位观测值的最小二乘处理方法）的等价性，并提出了一种改进的模糊度函数模型，该模型具有更好的求解特性，在此基础上，给出了模糊度函数法及改进的模糊度函数法的精度评定方法。最后通过对实测数据的处理，证明了理论和方法的有效性。     

北斗卫星伪距码偏差特性及其影响分析

在北斗导航卫星伪距码偏差特性分析的基础上,建立了倾斜地球同步轨道卫星（IGSO）和中轨卫星（MEO）的伪距码偏差多项式改正模型;并利用星间单差宽巷小数周一致性,分析建立北斗地球同步轨道卫星（GEO）卫星伪距码偏差改正模型。采用武汉大学北斗试验网、中国陆态网络和MGEX网不同位置、不同类型接收机观测数据进行分析验证,结果表明,北斗卫星伪距码偏差特性与观测值频率、卫星类型相关,所有GEO和IGSO卫星变化规律相同,所有MEO卫星变化规律相同,与接收机类型、测站位置和观测时间无关,偏差值大小随卫星高度角变化,其变化规律稳定,可以采用建立的两类改正模型（GEO/IGSO和MEO）进行修正。通过偏差修正后的伪距无电离层组合的残差、双频SPP以及单频PPP三个方面验证了伪距码偏差改正模型的正确性。     

A synthetic Earth for use in geodesy

 A synthetic Earth and its gravity field that can be represented at different resolutions for testing and comparing existing and new methods used for global gravity-field determination are created. Both the boundary and boundary values of the gravity potential can be generated. The approach chosen also allows observables to be generated at aircraft flight height or at satellite altitude. The generation of the synthetic Earth shape (SES) and gravity-field quantities is based upon spherical harmonic expansions of the isostatically compensated equivalent rock topography and the EGM96 global geopotential model. Spherical harmonic models are developed for both the synthetic Earth topography (SET) and the synthetic Earth potential (SEP) up to degree and order 2160 corresponding to a 5′×5′ resolution. Various sets of SET, SES and SEP with boundary geometry and boundary values at different resolutions can be generated using low-pass filters applied to the expansions. The representation is achieved in point sets based upon refined triangulation of a octahedral geometry projected onto the chosen reference ellipsoid. The filter cut-offs relate to the sampling pattern in order to avoid aliasing effects. Examples of the SET and its gravity field are shown for a resolution with a Nyquist sampling rate of 8.27 degrees. Received: 6 August 1999 / Accepted: 26 April 2000     

小波分析在单频GPS快速精密定位中的应用研究

单频GPS快速定位方程是严重病态的,应用最小二乘原理得到的模糊度浮点解大大偏离其准确值,应用LAMBDA方法难以正确地固定模糊度。本文将GPS载波相位双差观测量在不同的小波空间和尺度空间进行分解和重构,去除高频测量噪声,可减小测量噪声对GPS快速定位中病态方程解的影响,提高模糊度浮点解的精度,缩小模糊度搜索空间。实验表明,对于GPS短基线,仅利用1min左右的单频载波观测数据,经过基于haar、db4、coif4和sym4小波的5尺度小波变换后,可获得较准确的模糊度浮点解,应用LAMBDA法可正确地固定模糊度,达到厘米级定位精度。     

Precision, volume and eigenspectra for GPS ambiguity estimation based on the time-averaged satellite geometry

In this contribution we consider the time-averaged GPS single-baseline model and study in a qualitative sense its relation with the geometry-free model and the geometry-based model. The least-squares estimators of the model are derived and their properties discussed. Special attention is given to the ambiguity search space, since it plays such a crucial role in the problem of integer ambiguity estimation and validation. Easy-to-evaluate, closed-form expressions are presented for the volumes of the ambiguity search spaces that belong to the geometry-free model, the single-epoch geometry-based model and the time-averaged model. By means of an eigenvalue analysis, the geometry of the ambiguity search spaces is revealed and its impact on the search for the integer least-squares ambiguities discussed. Received: 3 April 1996; Accepted: 6 January 1997     

RINEX_HO: second- and third-order ionospheric corrections for RINEX observation files

When GNSS receivers capable of collecting dual-frequency data are available, it is possible to eliminate the first-order ionospheric effect in the data processing through the ionosphere-free linear combination. However, the second- and third-order ionospheric effects still remain. The first-, second- and third-order ionospheric effects are directly proportional to the total electron content (TEC), although the second- and third-order effects are influenced, respectively, by the geomagnetic field and the maximum electron density. In recent years, the international scientific community has given more attention to these kinds of effects and some works have shown that for high precision GNSS positioning these effects have to be taken into consideration. We present a software tool called RINEX_HO that was developed to correct GPS observables for second- and third-order ionosphere effects. RINEX_HO requires as input a RINEX observation file, then computes the second- and third-order ionospheric effects, and applies the corrections to the original GPS observables, creating a corrected RINEX file. The mathematical models implemented to compute these effects are presented, as well as the transformations involving the earth’s magnetic field. The use of TEC from global ionospheric maps and TEC calculated from raw pseudorange measurements or pseudoranges smoothed by phase is also investigated.     

星载GNSS-R海浪有效波高反演模型构建EI北大核心CSCD

海浪有效波高是描述海况的重要参数,在海浪和海洋动力学的预测中起着重要的作用。然而,卫星雷达高度计和浮标等传统方法难以实现高时空分辨率的SWH估计。星载GNSS-R为估计SWH提供了一种思路。本文提出了一种基于归一化积分延迟波形反演星载GNSS-R海浪SWH的方法。首先,对星载GNSS-R延迟多普勒图进行去噪处理和数据过滤以便对DDM数据进行严格质量控制。然后,从DDM中提取NIDW,并基于NIDW计算4个GNSS-R观测值(即归一化积分时延波形的前沿斜率和前沿波形值之和,归一化积分时延波形的后沿斜率和后沿波形值之和)。随后,基于这4个观测值建立了星载GNSS-R海浪SWH反演经验模型。最后,分别将ERA5和AVISO SWH数据产品作为验证数据,并将反演模型的SWH估计结果与ERA5和AVISO SWH数据产品进行比较和分析。试验结果表明,当采用ERA5 SWH数据作为验证数据时,4个观测值反演SWH的均方根误差和相关系数分别优于0.66 m和0.65;当采用AVISO SWH数据作为验证数据时,4个观测值反演SWH的均方根误差和相关系数分别优于0.68 m和0.70。进一步表明了本文建模方法在星载GNSS-R SWH估计方面具有可行性和可靠性。     

Polarimetric formulation of the visibility function equation including cross-polar antenna patterns

The European Space Agency's Soil Moisture and Ocean Salinity (SMOS) mission will be the first one using two-dimensional aperture synthesis radiometry for Earth observation. This study presents the formulation that relates instrument observables and brightness temperature maps including cross-polar antenna voltage patterns, which may be also different from element to element. Finally, the radiometric accuracy degradation if cross-polar patterns are neglected in the image reconstruction is studied.     

Rapid GPS ambiguity resolution for short and long baselines

 A method of quick initial carrier cycle ambiguity resolution is described. The method applies to high-quality dual-band global positioning system observations. Code measurements on both frequencies must be available. The rapidity of the method is achieved through smoothing pseudoranges by phase observables and forming linear combinations between the phase observables. Two cases are investigated. Case 1: ionospheric bias is neglected (short distances); and case 2: the bias is taken into account (longer distances, more than, say, 10 km). The method was tested on six baselines, from 1 to 31 km long. In most cases, single-epoch ambiguity resolution was achieved. Received: 6 October 1999 / Accepted: 4 March 2002     

Examining the C1-P1 Pseudorange Bias

Use of GPS tracking data from different dual-frequency receiver types (cross-correlating vs. codeless) has revealed satellite-dependent biases in pseudorange observables P1 (Y-code) and C1 (C/A, Clear Acquisition code). These biases can have a direct effect on clock estimates, carrier phase bias fixing, and other parameters estimated in GPS data processing. A set of satellite-specific compensatory pseudorange offsets is calculated, and each is applied to a wee of daily global network analyses in which satlellite, receiver, atmospheric, and Earth rotation parameters are estimated. Results from these analyses are then compared to those from corresponding baseline cases in which no biases were applied. There is also some evidence that suggests that the pseudorange biases differ even among codeless receiver models. Hence, a second set of offsets is computed on a different basis, and compared with the baseline model in a similar manner. A preliminary examination of C1-P1 variations over time is presented. Finally, recommendations are made for the use of the calculated offsets, and consideration is given to a future dissemination of updates to these values as necessary. ? 2001 John Wiley & Sons, Inc.