On the accurate numerical evaluation of geodetic convolution integrals

In the numerical evaluation of geodetic convolution integrals, whether by quadrature or discrete/fast Fourier transform (D/FFT) techniques, the integration kernel is sometimes computed at the centre of the discretised grid cells. For singular kernels—a common case in physical geodesy—this approximation produces significant errors near the computation point, where the kernel changes rapidly across the cell. Rigorously, mean kernels across each whole cell are required. We present one numerical and one analytical method capable of providing estimates of mean kernels for convolution integrals. The numerical method is based on Gauss-Legendre quadrature (GLQ) as efficient integration technique. The analytical approach is based on kernel weighting factors, computed in planar approximation close to the computation point, and used to convert non-planar kernels from point to mean representation. A numerical study exemplifies the benefits of using mean kernels in Stokes’s integral. The method is validated using closed-loop tests based on the EGM2008 global gravity model, revealing that using mean kernels instead of point kernels reduces numerical integration errors by a factor of ~5 (at a grid-resolution of 10 arc min). Analytical mean kernel solutions are then derived for 14 other commonly used geodetic convolution integrals: Hotine, Eötvös, Green-Molodensky, tidal displacement, ocean tide loading, deflection-geoid, Vening-Meinesz, inverse Vening-Meinesz, inverse Stokes, inverse Hotine, terrain correction, primary indirect effect, Molodensky’s G1 term and the Poisson integral. We recommend that mean kernels be used to accurately evaluate geodetic convolution integrals, and the two methods presented here are effective and easy to implement.     

CHECKING TRAVERSE COMPUTATIONS

Abstract

Traverse Computations must be Checked.—A traverse is a chain of points connected by angular and linear measurements. The check on observations is provided by the agreement, obtained in computations, between the terminals of the traverse (terminal bearings and terminal co-ordinates) taken as fixed. This check is not sufficient, however, to serve as a check on the computations. As a matter of principle, computations should be free of errors; there are no limits of tolerance in computational work except for discrepancies arising from inaccuracy of last figures. Secondly, errors in computation may occur that are not revealed by the traverse misclosures, not to speak of compensational errors, the field for which is very favourable in traverse work.     

Improved efficiency of maximum likelihood analysis of time series with temporally correlated errors

Most time series of geophysical phenomena have temporally correlated errors. From these measurements, various parameters are estimated. For instance, from geodetic measurements of positions, the rates and changes in rates are often estimated and are used to model tectonic processes. Along with the estimates of the size of the parameters, the error in these parameters needs to be assessed. If temporal correlations are not taken into account, or each observation is assumed to be independent, it is likely that any estimate of the error of these parameters will be too low and the estimated value of the parameter will be biased. Inclusion of better estimates of uncertainties is limited by several factors, including selection of the correct model for the background noise and the computational requirements to estimate the parameters of the selected noise model for cases where there are numerous observations. Here, I address the second problem of computational efficiency using maximum likelihood estimates (MLE). Most geophysical time series have background noise processes that can be represented as a combination of white and power-law noise, \(1/f^{\alpha }\) with frequency, f. With missing data, standard spectral techniques involving FFTs are not appropriate. Instead, time domain techniques involving construction and inversion of large data covariance matrices are employed. Bos et al. (J Geod, 2013. doi: 10.1007/s00190-012-0605-0) demonstrate one technique that substantially increases the efficiency of the MLE methods, yet is only an approximate solution for power-law indices >1.0 since they require the data covariance matrix to be Toeplitz. That restriction can be removed by simply forming a data filter that adds noise processes rather than combining them in quadrature. Consequently, the inversion of the data covariance matrix is simplified yet provides robust results for a wider range of power-law indices.     

正态-Gamma先验下粗差探测的Bayes方法

本文主要运用Bayes统计推断的基本原理,提出和建立将观测信息与正态——Gamma先验信息融合、基于综合信息判断粗差的Bayes方法。首先,根据Bayes统计推断的基本原理,建立了判断粗差的Bayes方法——后验概率法,然后针对测量平差实际,考虑未知参数的正态-Gamma先验信息,分别给出了非等权独立观测条件下基于均值漂移模型和方差膨胀模型的后验概率的具体计算公式,并给出了相应的粗差的Bayes估算方法。最后提出了基于后验概率进行粗差探测的实施过程和具体步骤。数值试验的结果表明,本文提出的粗差探测的Bayes方法对多个粗差的同时定位和定值是相当有效的。     

结合影响分析的抗隐差型Bayes粗差探测方法

通过分析隐差现象产生的原因,将基于观测误差的后验概率法与基于Kullback-Leiber距离的影响分析相结合,并利用逐个搜索法的思想,给出了一种粗差探测的抗隐差型Bayes方法。实验结果表明,采用该方法进行粗差探测切实可行,不但有效地发现了被掩盖的粗差,而且运算简单,效率较高。     

Computing components of the gravity field induced by distant topographic masses and condensed masses over the entire Earth using the 1-D FFT approach

 A new method for computing gravitational potential and attraction induced by distant, global masses on a global scale has been developed. The method uses series expansions and the well known one-dimensional fast Fourier transform (1-D FFT) method. It has been proven to be significantly faster than quadrature while being equally accurate. Various quantities were studied to cover the two primary applications of the Stokes–Helmert scheme of modeling effects. These two applications (or paths), given the names R/r/D and R/D/r, are briefly discussed, although the primary objective of the paper is to provide computational information to either path, rather than choosing one path as preferable to the other. It is further shown that the impact of masses outside a 4-degree cap can impact the absolute computation of the geoid at more than 1 cm, and should therefore be included in all local geoid computations seeking that accuracy. Received: 13 December 2000 / Accepted: 3 September 2001     

Development of the second-order derivatives of the Earth’s potential in the local north-oriented reference frame in orthogonal series of modified spherical harmonics

The second-order derivatives of the Earth’s potential in the local north-oriented reference frame are expanded in series of modified spherical harmonics. Linear relations are derived between the spectral coefficients of these series and the spectrum of the geopotential. On the basis of these relations, recurrence procedures are developed for evaluating the geopotential coefficients from the spectrum of each derivative and, inversely, for simulating the latter from a known geopotential model. Very simple structure of the derived expressions for the derivatives is convenient for estimating the geopotential coefficients by the least-squares procedure, at a certain step of processing satellite gradiometry data. Due to the orthogonality of the new series, the quadrature formula approach can be also applied, which allows avoidance of aliasing errors caused by the series truncation. The spectral coefficients of the derivatives are evaluated on the basis of the derived relations from the geopotential models EGM96 and EIGEN-CG01C at a mean orbital sphere of the GOCE satellite. Various characteristics of the spectra are studied corresponding to the EGM96 model. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Ultra-tight GPS/INS/PL integration: a system concept and performance analysis

The architecture of the ultra-tight GPS/INS/PL integration is the key to its successful performance; the main feature of this architecture is the Doppler feedback to the GPS receiver tracking loops. This Doppler derived from INS, when integrated with the carrier tracking loops, removes the Doppler due to vehicle dynamics from the GPS/PL signal thereby achieving a significant reduction in the carrier tracking loop bandwidth. The bandwidth reduction provides several advantages such as: improvement in anti-jamming performance, and increase in post correlated signal strength which in turn increases the dynamic range and accuracy of measurements. Therefore, any degradation in the derived Doppler estimates will directly affect the tracking loop bandwidth and hence its performance. The quadrature signals from the receiver correlator, I (in-phase) and Q (quadrature), form the measurements, whereas the inertial sensor errors, position, velocity and attitude errors form the states of the complementary Kalman filter. To specify a reliable measurement model of the filter for this type of integrated system, a good understanding of GPS/PL signal characteristics is essential. It is shown in this paper that phase and frequency errors are the variables that relate the measurements and the states in the Kalman filter. The main focus of this paper is to establish the fundamental mathematical relationships that form the measurement model, and to show explicitly how the system error states are related to the GPS/PL signals. The derived mathematical relationships encapsulated in a Kalman filter, are tested by simulation and shown to be valid.

---

THE HEIGHT CLOSING ERROR IN MULTIPLEX BRIDGES CAUSED BY NON-FLAT DIAPOSITIVES

Abstract

A recent investigation into the flatness of Multiplex diapositive slides has shown that flatness errors occur ranging from 0 to 0·03–0·04 mm. referred to the flat projector stage. In a first attempt to ascertain the effects of errors of this kind on Multiplex bridges, the flatness values of two sets of nine diapositives each were measured using a simple interferometric method and the vertical (wants of correspondence) and horizontal parallaxes introduced by these flatness errors were subjected to computational bridging. The resulting height errors at the end of the two strips proved to be of noticeable size, as large as +0·6 and ?0·9 mm. respectively. Indicative as these figures may already be, it seems useful to abstract the investigation from the vagaries of the individual case and to put it on a more general footing.     

粗差探测的Bayes方法

讨论粗差探测的Bayes方法。首先根据Bayes统计推断的基本原理,建立判断粗差的Bayes方法———后验概率法,然后针对测量平差实际,分别给出非等权独立观测条件下基于均值漂移模型和方差膨胀模型的后验概率的计算公式,最后结合一边角网算例,验证本文方法的效果。     

Navigation message demodulation for GPS receiver on-board spinning rockets

We describe a demodulation scheme for the navigation message of GPS receivers on spin-stabilized rockets. Doppler frequencies due to fast and complex dynamics, in particular high-rate spin, cause errors in carrier frequency tracking. The effects of such errors on navigation message demodulation are described through theoretical analysis and numerical simulation. A demodulation scheme that includes a frequency estimator is proposed to account for frequency tracking errors. It is demonstrated that demodulation performance is degraded 5 dB due to frequency uncertainty. Simulation results showed that a demodulator which includes maximum likelihood (ML) frequency estimator achieves near-optimal symbol error rate under these conditions. Demodulation with ML estimator achieves a bit error rate below 10^?5 for a C/N ₀?=?35 dB–Hz, for spin rates below 2.7?Hz, and a rocket radius smaller than 1 m. For the cases in which computational capabilities of the on-board GPS receiver is insufficient to implement the demodulator with ML estimator, frequency estimation methods with low complexity were also tested through numerical simulation. The proposed Kay and Quinn-Fernandes combination achieves a bit error rate below 10^?5 for a C/N ₀?=?37 dB–Hz while requiring 1/10 of processing time.     

Gross Error Compensation for Gravity Field Analysis Based on Kinematic Orbit Data

This paper aims at a comparative study of several measures to compensate for gross errors in kinematic orbit data. It starts with a simulation study on the influence of a single outlier in the orbit data on the gravity field solution. It is shown that even a single outlier can degrade the resulting gravity field solution considerably. To compensate for outliers, two different strategies are investigated: wavelet filters, which detect and eliminate gross errors, and robust estimators, which due to an iterative downweighting gradually ignore those observations that lead to large residuals. Both methods are applied in the scope of the analysis of a 2-year kinematic CHAMP (challenging minisatellite payload) orbit data set. In various real data studies, robust estimators outperform wavelet filters in terms of resolution of the derived gravity field solution. This superior performance is at the cost of computational load, as robust estimators are implemented iteratively and require the solution of large sets of linear equations several times.     

基于方差膨胀模型的粗差探测Bayes方法

基于方差膨胀模型提出并建立了用观测信息的同时利用先验信息判断粗差的Bayes方法。首先,根据Bayes统计推断的基本原理,建立了判断粗差的Bayes方法——后验概率法;然后针对测量平差实际,考虑未知参数的两种先验信息分别给出了非等权独立观测条件下基于该模型的后验概率的计算公式;最后对模拟算例进行了计算和分析。试验结果表明,用给出的探测粗差的Bayes方法是切实可行的。     

A multi-technique approach for characterizing the SVN49 signal anomaly, part 2: chip shape analysis

Due to a satellite internal reflection at the L5 test payload, the SVN49 (PRN1) GPS satellite exhibits a static multipath on the L1 and L2 signals, which results in elevation-dependent tracking errors for terrestrial receivers. Using a 30-m high-gain antenna, code and carrier phase measurements as well as raw in-phase and quadrature radio frequency samples have been collected during a series of zenith passes in mid-April 2010 to characterize the SVN49 multipath and its impact on common users. Following an analysis of the receiver tracking data and the IQ constellation provided in Part 1 of this study, the present Part 2 provides an in-depth investigation into chip shapes for the L1 and L2 signals. A single reflection model is found to be compatible with the observed chip shape distortions and key parameters for an elevation dependent multipath model are derived. A good agreement is found between multipath parameters derived independently from raw IQ-samples and measurements of a so-called Vision Correlator. The chip shapes and their observed variation with elevation can be used to predict the multipath response of different correlator types within a tracking receiver. The multipath model itself is suitable for implementation in a signal simulator and thus enables laboratory testing of actual receiver hardware.     

Spheroidal equal angular DEMs: The specificity of morphometric treatment

Digital elevation models (DEMs) are commonly constructed using two main types of regular grids: plane square grids and spheroidal equal angular grids. Methods and algorithms intended for plane square‐gridded DEMs should not be directly applied to spheroidal equal angular DEMs. This is because these grids have fundamentally different geometry. However, some researchers continue to apply square‐grid algorithms to spheroidal equal angular DEMs. It seems appropriate to consider once again the specifity of morphometric treatment of spheroidal equal angular DEMs. This article, first, demonstrates possibilities of direct calculation of local, nonlocal, and combined morphometric variables from spheroidal equal angular DEMs exemplified by slope gradient, catchment area, and topographic index. Second, the article shows computational errors when algorithms for plane square‐gridded DEMs are unreasonably applied to spheroidal equal angular DEMs. The study is exemplified by two DEMs. A medium‐resolution DEM of a relatively small, high‐mountainous area (Mount Elbrus) was extracted from the SRTM1 DEM. A low‐resolution DEM of a vast region with the diverse topography (the central and western regions of Kenya) was extracted from the SRTM30_PLUS DEM. The results show that application of square‐grid methods to spheroidal equal angular DEMs leads to substantial computational errors in models of morphometric variables.     

利用GOCE模拟观测反演重力场的Torus法

在介绍Torus方法反演地球重力场模型的基本原理和方法的基础上,基于圆环面上均匀分布的卫星引力梯度模拟观测值解算了200阶次的地球重力场模型,在无误差情况下,Torus方法解算模型的阶误差RMS小于10-16,验证了该方法的严密性。利用61dGOCE卫星轨道上无误差的模拟引力梯度观测值解算了200阶次的地球重力场模型,分析了格网化误差、极空白对解算精度的影响,迭代3次后,在不考虑低次系数情况下,模型的大地水准面阶误差和累积误差均较小,最大值仅为0.022mm和0.099mm。在沿轨卫星引力梯度模拟数据中加入5mE/Hz1/2的白噪声,基于Torus方法和空域最小二乘法解算了200阶次的地球重力场模型,Torus方法的精度略低于空域最小二乘法的精度,在不考虑低次项的情况下,两种方法解算模型的大地水准面阶误差最大值分别为1.58cm和1.45cm,累积误差最大值分别为6.37cm和5.55cm。但由于采用了二维快速傅里叶技术和块对角最小二乘法,极大地提高了计算效率。本文数值结果说明Torus方法是一种独立有效的方法,可用于GOCE任务海量卫星引力梯度观测值反演重力场的快速解算。     

基于空间双曲交会的GNSS伪距单点定位算法

针对GNSS卫星导航中的伪距单点定位,提出一种不需要测站坐标近似值的非迭代算法。该算法将GNSS伪距导航定位方程转化为空间双曲定位方程,给出具体的解算步骤,研究了空间双曲定位方程的解(有两解),利用GNSS伪距导航定位的特点可消除多值性,从而实现无初值GNSS伪距单点定位。该算法与Bancroft算法相比,通过星间单差,与测站有关的公共误差项被消去,提高了定位精度;与传统的迭代算法相比,提高了计算效率,而且不需要测站坐标初值。最后通过IGS监测站实测数据对3种算法进行比较,验证了算法的有效性。     

Modeling errors in upward continuation for INS gravity compensation

Accurate upward continuation of gravity anomalies supports future precision, free-inertial navigation systems, since the latter cannot by themselves sense the gravitational field and thus require appropriate gravity compensation. This compensation is in the form of horizontal gravity components. An analysis of the model errors in upward continuation using derivatives of the standard Pizzetti integral solution (spherical approximation) shows that discretization of the data and truncation of the integral are the major sources of error in the predicted horizontal components of the gravity disturbance. The irregular shape of the data boundary, even the relatively rough topography of a simulated mountainous region, has only secondary effect, except when the data resolution is very high (small discretization error). Other errors due to spherical approximation are even less important. The analysis excluded all measurement errors in the gravity anomaly data in order to quantify just the model errors. Based on a consistent gravity field/topographic surface simulation, upward continuation errors in the derivatives of the Pizzetti integral to mean altitudes of about 3,000 and 1,500 m above the mean surface ranged from less than 1 mGal (standard deviation) to less than 2 mGal (standard deviation), respectively, in the case of 2 arcmin data resolution. Least-squares collocation performs better than this, but may require significantly greater computational resources.     

星地无线电双向时间比对法误差源分析

本文简要介绍了星地无线电双向时间比对法的基本原理,给出其地心惯性系中计算模型,分析了各种误差源对星地间相对钟差的影响,计算了在几种计算精度下卫星和地面站位置和速度误差所需要达到的精度.结果表明:影响星地无线电双向法精度的主要误差源为设备时延误差和电离层延迟误差,而由于卫星和地面站位置不准确引起的误差相对较小;对于电离层...     

A Bayesian approach to the detection of gross errors based on posterior probability

Existing methods for gross error detection, based on the mean shift model or the variance inflation model, have hardly considered or taken advantage of the potential prior information on the unknown parameters. This paper puts forward a Bayesian approach for gross error detection when prior information on the unknown parameters is available. Firstly, based on the basic principle of Bayesian statistical inference, the Bayesian method—posterior probability method—for the detection of gross errors is established. Secondly, considering either non-informative priors or normal-gamma priors on the unknown parameters, the computational formula of the posterior probability is given for both the mean shift model and the variance inflation model, respectively, under the condition of unequal weight and independent observations. Finally, as an example, a triangulation network is computed and analyzed, which shows that the method given here is feasible.