首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
On the adjustment of combined GPS/levelling/geoid networks   总被引:12,自引:7,他引:5  
A detailed treatment of adjustment problems in combined global positioning system (GPS)/levelling/geoid networks is given. The two main types of `unknowns' in this kind of multi-data 1D networks are usually the gravimetric geoid accuracy and a 2D spatial field that describes all the datum/systematic distortions among the available height data sets. An accurate knowledge of the latter becomes especially important when we consider employing GPS techniques for levelling purposes with respect to a local vertical datum. Two modelling alternatives for the correction field are presented, namely a pure deterministic parametric model, and a hybrid deterministic and stochastic model. The concept of variance component estimation is also proposed as an important statistical tool for assessing the actual gravimetric geoid noise level and/or testing a priori determined geoid error models. Finally, conclusions are drawn and recommendations for further study are suggested. Received: 9 September 1998 / Accepted: 8 June 1999  相似文献   

2.
Summary A datum change between two geodetic systems with points in common may be derived in three stages; slight adjustments of coordinates to make the networks of common points geometrically similar in the two systems; a scale factor to make them geometrically congruent; finally, an orthogonal transformation to swing them into coincidence. The geometrical concept is developed of a “datum screw”, not arbitrarily chosen as is the “origin” or “datum point” of a geodetic survey, but intrinsic to the geometry. The conditions under which it degenerates to a simple “datum shift” are discussed. Differential and other formulae for changes of spheroid and of datum are given, together with a set of tables of coefficients.  相似文献   

3.
 The weighted Procrustes algorithm is presented as a very effective tool for solving the three-dimensional datum transformation problem. In particular, the weighted Procrustes algorithm does not require any initial datum parameters for linearization or any iteration procedure. As a closed-form algorithm it only requires the values of Cartesian coordinates in both systems of reference. Where there is some prior information about the variance–covariance matrix of the two sets of Cartesian coordinates, also called pseudo-observations, the weighted Procrustes algorithm is able to incorporate such a quality property of the input data by means of a proper choice of weight matrix. Such a choice is based on a properly designed criterion matrix which is discussed in detail. Thanks to the weighted Procrustes algorithm, the problem of incorporating the stochasticity measures of both systems of coordinates involved in the seven parameter datum transformation problem [conformal group ℂ7(3)] which is free of linearization and any iterative procedure can be considered to be solved. Illustrative examples are given. Received: 7 January 2002 / Accepted: 9 September 2002 Correspondence to: E. W. Grafarend  相似文献   

4.
Robust estimation of geodetic datum transformation   总被引:18,自引:1,他引:17  
Y. Yang 《Journal of Geodesy》1999,73(5):268-274
The robust estimation of geodetic datum transformation is discussed. The basic principle of robust estimation is introduced. The error influence functions of the robust estimators, together with those of least-squares estimators, are given. Particular attention is given to the robust initial estimates of the transformation parameters, which should have a high breakdown point in order to provide reliable residuals for the following estimation. The median method is applied to solve for robust initial estimates of transformation parameters since it has the highest breakdown point. A smooth weight function is then used to improve the efficiency of the parameter estimates in successive iterative computations. A numerical example is given on a datum transformation between a global positioning system network and the corresponding geodetic network in China. The results show that when the coordinates are contaminated by outliers, the proposed method can still give reasonable results. Received: 25 September 1997 / Accepted: 1 March 1999  相似文献   

5.
A review of recent progress and current activities towards an improved formulation and solution of geodetic boundary value problems is given. Improvements stimulated and required by the dramatic changes of the real world of geodetic measurements are focused upon. Altimetry–gravimetry problems taking into account various scenarios of non-homogeneous data coverage are discussed in detail. Other problems are related to free geodetic datum parameters, most of all the vertical datum, overdetermination or additional constraints imposed by satellite geodetic observations or models. Some brief remarks are made on pseudo-boundary value problems for geoid determination and on purely gravitational boundary-value problems. Received: 17 March 1999 / Accepted: 19 April 1999  相似文献   

6.
 Modern dual-frequency global positioning system (GPS) receivers are capable of providing direct measurements of both L1 C/A (C1) and P code (P1) without the use of the Y-codes under Anti-Spoofing. A discrepancy or bias between the C1 and P1 measurements from these receivers has however been of concern to operators and users of GPS reference networks. For the purpose of modeling and estimation, the nature and characteristics of the discrepancy must be investigated. The research results presented indicate that the discrepancy between the C1 and P1 measurements contains two different types of components: one is of constant type while another is time variant. A method has been developed for their modeling and estimation. The residual C1–P1 time series after a satellite-dependent bias removal agree at a few-centimeter level, indicating the effectiveness of the proposed model. This allows the C1–P1 discrepancy, both constant and non-constant components, to be removed from GPS reference network solutions. Numerical results are provided to support the analysis. Received: 12 October 1999 / Accepted: 11 May 2000  相似文献   

7.
Adaptive Kalman Filtering for INS/GPS   总被引:69,自引:0,他引:69  
After reviewing the two main approaches of adaptive Kalman filtering, namely, innovation-based adaptive estimation (IAE) and multiple-model-based adaptive estimation (MMAE), the detailed development of an innovation-based adaptive Kalman filter for an integrated inertial navigation system/global positioning system (INS/GPS) is given. The developed adaptive Kalman filter is based on the maximum likelihood criterion for the proper choice of the filter weight and hence the filter gain factors. Results from two kinematic field tests in which the INS/GPS was compared to highly precise reference data are presented. Results show that the adaptive Kalman filter outperforms the conventional Kalman filter by tuning either the system noise variance–covariance (V–C) matrix `Q' or the update measurement noise V–C matrix `R' or both of them. Received: 14 September 1998 / Accepted: 21 December 1998  相似文献   

8.
It is shown that also in a rank deficient Gauss-Markov model higher weights of the observations automatically improve the precision of the estimated parameters as long as they are computed in thesame datum. However, the amount of improvement in terms of the trace of the dispersion matrix isminimum for the so-called “free datum” which corresponds to the pseudo-inverse normal equations matrix. This behaviour together with its consequences is discussed by an example with special emphasis on geodetic networks for deformation analysis.  相似文献   

9.
 This paper suggests that potential coefficient models of the Earth's gravitational potential be used to calculate height anomalies which are then reduced to geoid undulations where such quantities are needed for orthometric height determination and vertical datum definition through a potential coefficient realization of the geoid. The process of the conversion of the height anomaly into a geoid undulation is represented by a height anomaly gradient term and the usual N–ζ term that is dependent on elevation and the Bouguer anomaly. Using a degree 360 expansion of 30′ elevations and the OSU91A potential coefficient model, a degree 360 representation of the correction terms was computed. The magnitude of N–ζ reached –3.4 m in the Himalaya Mountains with smaller, but still significant, magnitudes in other mountainous regions. Received: 6 May 1996; Accepted: 30 October 1996  相似文献   

10.
A Quaternion-based Geodetic Datum Transformation Algorithm   总被引:1,自引:1,他引:1  
This paper briefly introduces quaternions to represent rotation parameters and then derives the formulae to compute quaternion, translation and scale parameters in the Bursa–Wolf geodetic datum transformation model from two sets of co-located 3D coordinates. The main advantage of this representation is that linearization and iteration are not needed for the computation of the datum transformation parameters. We further extend the formulae to compute quaternion-based datum transformation parameters under constraints such as the distance between two fixed stations, and develop the corresponding iteration algorithm. Finally, two numerical case studies are presented to demonstrate the applications of the derived formulae.  相似文献   

11.
On aggregation in spatial econometric modelling   总被引:1,自引:1,他引:0  
The spatial aggregation problem – also termed the modifiable areal unit problem – has attracted regular attention in spatial statistics and econometrics. In this study econometric aggregation analysis is used to investigate the formal composition of meso-areal parameters given micro-areal underlying relations with spatial dependence. Impact on stochastic terms (possible meso-areal spatial autocorrelation) is also studied. Finally consequences for meso-areal estimation are derived, the general finding having been that spatial aggregation leads to meso-region specific parameter values, with the estimation problems this implies.  相似文献   

12.
Summary From a two-dimensional network adjustment it is well understood that the one orientation unknown of a theodolite frame is estimable, once the orientation datum parameter, e.g., one azimuth, is fixed. In three-dimensional networks the problem of estimability of three orientation unknowns inherent in a theodolite frame is more complex. Here we prove that not only the classical horizontal orientation unknown is estimable (up to the datum degrees of freedom), but also astronomical longitude and astronomical latitude which can be considered as two additional orientation unknowns of the theodolite frame moving with respect to an earth-fixed equatorial frame of reference. Thus the theodolite instrument can be considered—at least theoretically—a gradiometer measuring the variation of the directional parameters of the gravity vector from one point to another. Or up to the datum degrees of freedom astronomical longitude and astronomical latitude can be determined from only theodolite observations between exclusively terrestrial points. M?nicke (1982), has shown that despite the refraction problem the method works sufficiently well in practice.  相似文献   

13.
The crossover adjustment plays a central role in the processing of satellite altimeter measurements. The usual procedure is to form sea surface height differences at crossover points, solve for the radial orbit error (with due attention to the singular nature of the estimation problem) and then to construct altimetric sea-level maps using the mean sea surface heights at the crossover points. Our approach is very different, to make direct use of measurements at crossover points without differencing and to estimate simultaneously orbit parameters, mean sea surface height and sea surface height variability in a single, unified adjustment. The technique is suited for repeat data over an area small enough that adjoining passes may be considered to be parallel and to permit the solution of a set of linear equations of dimension equal to the number of crossover points. The size of the numerical problem is almost independent of the number of repeat cycles of the altimeter mission. Explicit recognition is given to the rank defect of the least-squares estimation problem; we show that, for an orbit model with r parameters, the rank defect of the local crossover problem is exactly r 2. The defect may be overcome by choosing an appropriate set of constraints – either giving a best fit of mean sea surface heights to a reference surface, or minimising orbit parameters, or a minimum norm solution in which both mean sea surface heights and orbit parameters are minimised. There is no need to choose a reference pass, all passes are treated equally and data gaps are easily accommodated. Numerical results are presented for the south-western Indian Ocean, based on the first 2 years of altimeter data from the Geosat Exact Repeat Mission. Received: 31 May 1996 / Accepted: 19 April 1997  相似文献   

14.
 The analysis of lunar laser ranging (LLR) data enables the determination of many parameters of the Earth–Moon system, such as lunar gravity coefficients, reflector and station coordinates which contribute to the realisation of the International Terrestrial Reference Frame 2000 (ITRF 2000), Earth orientation parameters [EOPs, which contribute to the global EOP solutions at the International Earth Rotation Service (IERS)] or quantities which parameterise relativistic effects in the solar system. The big advantage of LLR is the long time span of lunar observations (1970–2000). The accuracy of the normal points nowadays is about 1 cm.  The capability of LLR to determine tidal parameters is investigated. In principle, it could be assumed that LLR would contribute greatly to the investigation of tidal effects, because the Moon is the most important tide-generating body. In this respect some special topics such as treatment of the permanent tide and the effect of atmospheric loading are addressed and results for the tidal parameters h 2 and l 2 as well as values for the eight main tides are given. Received: 14 August 2000 / Accepted: 15 October 2001  相似文献   

15.
Terrestrial free-air gravity anomalies form a most essential data source in the framework of gravity field determination. Gravity anomalies depend on the datums of the gravity, vertical, and horizontal networks as well as on the definition of a normal gravity field; thus gravity anomaly data are affected in a systematic way by inconsistencies of the local datums with respect to a global datum, by the use of a simplified free-air reduction procedure and of different kinds of height system. These systematic errors in free-air gravity anomaly data cause systematic effects in gravity field related quantities like e.g. absolute and relative geoidal heights or height anomalies calculated from gravity anomaly data. In detail it is shown that the effects of horizontal datum inconsistencies have been underestimated in the past. The corresponding systematic errors in gravity anomalies are maximum in mid-latitudes and can be as large as the errors induced by gravity and vertical datum and height system inconsistencies. As an example the situation in Australia is evaluated in more detail: The deviations between the national Australian horizontal datum and a global datum produce a systematic error in the free-air gravity anomalies of about −0.10 mgal which value is nearly constant over the continent  相似文献   

16.
Three-dimensional ray tracing through a numerical weather model has been applied to a global precise point positioning (PPP) campaign for modeling both the elevation angle- and azimuth-dependence of the tropospheric delay. Rather than applying the ray-traced slant delays directly, the delay has been parameterized in terms of slant factors, which are applied in a similar manner to traditional mapping functions, but which can account for the azimuthal asymmetry of the delay. Five strategies are considered: (1) Vienna Mapping Functions 1 (VMF1) and estimation of a residual zenith delay parameter; (2) VMF1, estimation of a residual zenith delay and estimation of two tropospheric gradient parameters; (3) three-dimensional ray-traced slant factors and estimation of a residual zenith delay; (4) using only ray-traced slant factors and no estimation of any tropospheric parameters and; (5) using both ray-traced slant factors and estimating a residual zenith delay and two tropospheric gradient parameters. The use of the ray-traced slant factors (solution 3) showed a 3.8% improvement in the repeatability of the up component when compared to the assumption of a symmetric atmosphere (solution 1), while the estimation of two tropospheric gradient parameters gave the best results showing an 7.6% improvement over solution 1 in the up component. Solution 4 performed well in the horizontal domain, allowing for sub-centimeter repeatability but the up component was degraded due to deficiencies in the modeling of the zenith delay, particularly for stations located at equatorial latitudes. The magnitude of the differences in the mean coordinates between solution 2 and solution 3, and the strong correlation with the differences between the north component and the ray-traced gradients (coefficient of correlation of 0.83), as well as the impact of observation geometry on the gradient solution indicate that the use of the ray-traced slant factors could have an implication on the realization of reference frames. The estimated tropospheric products from the PPP solutions were compared to those derived from ray tracing. For the zenith delay, a root mean square (RMS) of 5.4 mm was found, while for the gradient terms, a correlation coefficient of 0.46 for the N–S and 0.42 for the E–W was found for the north–south and east–west components, suggesting that there are still important differences in the gradient parameters which could be due to either errors in the NWM or to non-tropospheric error sources leaking into the PPP-estimated gradients.  相似文献   

17.
 Since the beginning of the International Global Navigation Satellite System (GLONASS) Experiment, IGEX, in October 1998, the Center for Orbit Determination in Europe (CODE) has acted as an analysis center providing precise GLONASS orbits on a regular basis. In CODE's IGEX routine analysis the Global Positioning System (GPS) orbits and Earth rotation parameters are introduced as known quantities into the GLONASS processing. A new approach is studied, where data from the IGEX network are combined with GPS observations from the International GPS Service (IGS) network and all parameters (GPS and GLONASS orbits, Earth rotation parameters, and site coordinates) are estimated in one processing step. The influence of different solar radiation pressure parameterizations on the GLONASS orbits is studied using different parameter subsets of the extended CODE orbit model. Parameterization with three constant terms in the three orthogonal directions, D, Y, and X (D = direction satellite–Sun, Y = direction of the satellite's solar panel axis), and two periodic terms in the X-direction, proves to be adequate for GLONASS satellites. As a result of the processing it is found that the solar radiation pressure effect for the GLONASS satellites is significantly different in the Y-direction from that for the GPS satellites, and an extensive analysis is carried out to investigate the effect in detail. SLR observations from the ILRS network are used as an independent check on the quality of the GLONASS orbital solutions. Both processing aspects, combining the two networks and changing the orbit parameterization, significantly improve the quality of the determined GLONASS orbits compared to the orbits stemming from CODE's IGEX routine processing. Received: 10 May 2000 / Accepted: 9 October 2000  相似文献   

18.
Least-squares collocation may be used for the estimation of spherical harmonic coefficients and their error and error correlations from GOCE data. Due to the extremely large number of data, this requires the use of the so-called method of Fast Spherical Collocation (FSC) which requires that data is gridded equidistantly on each parallel and have the same uncorrelated noise on the parallel. A consequence of this is that error-covariances will be zero except between coefficients of the same signed order (i.e., the same order and the same coefficient type CC or SS). If the data distribution and the characteristics of the data noise are symmetric with respect to the equator, then, within a given order and coefficient type, the error-covariances amongst coefficients whose degrees are of different parity also vanish. The deviation from this “ideal” pattern has been studied using data-sets of second order radial derivatives of the anomalous potential. A total number of points below 17,000 were used having an equi-angular or an equal area distribution or being associated with points on a realistic GOCE orbit but close to the nodes of a grid. Also the data were considered having a correlated or an uncorrelated noise and three different signal covariance functions. Grids including data or not including data in the polar areas were used. Using the functionals associated with the data, error estimates of coefficients and error-correlations between coefficients were calculated up to a maximal degree and order equal to 90. As expected, for the data-distributions with no data in the polar areas the error-estimates were found to be larger than when the polar areas contained data. In all cases it was found that only the error-correlations between coefficients of the same order were significantly different from zero (up to 88%). Error-correlations were significantly larger when data had been regarded as having non-zero error-correlations. Also the error-correlations were largest when the covariance function with the largest signal covariance distance was used. The main finding of this study was that the correlated noise has more pronounced impact on gridded data than on data distributed on a realistic GOCE orbit. This is useful information for methods using gridded data, such as FSC.  相似文献   

19.
Transformation of coordinates between two horizontal geodetic datums   总被引:2,自引:1,他引:2  
The following topics are discussed in this paper: the geocentric coordinate system and its different realizations used in geodetic practice; the definition of a horizontal geodetic datum (reference ellipsoid) and its positioning and orientation with respect to the geocentric coordinate system; positions on a horizontal datum and errors inherent in the process of positioning; and distortions of geodetic networks referred to a horizontal datum. The problem of determining transformation parameters between a horizontal datum and the geocentric coordinate system from known positions is then analysed. It is often found necessary to transform positions from one horizontal datum to another. These transformations are normally accomplished through the geocentric coordinate system and they include the transformation parameters of the two datums as well as the representation of the respective network distortions. Problems encountered in putting these transformations together are pointed out.  相似文献   

20.
A new estimator for VLBI baseline length repeatability   总被引:1,自引:1,他引:0  
O. Titov 《Journal of Geodesy》2009,83(11):1041-1049
The goal of this paper is to introduce a more effective technique to approximate for the “repeatability–baseline length” relationship that is used to evaluate the quality of geodetic VLBI results. Traditionally, this relationship is approximated by a quadratic function of baseline length over all baselines. The new model incorporates the mean number of observed group delays of the reference radio sources (i.e. estimated as global parameters) used in the estimation of each baseline. It is shown that the new method provides a better approximation of the “repeatability–baseline length” relationship than the traditional model. Further development of the new approach comes down to modeling the repeatability as a function of two parameters: baseline length and baseline slewing rate. Within the framework of this new approach the station vertical and horizontal uncertainties can be treated as a function of baseline length. While the previous relationship indicated that the station vertical uncertainties are generally 4–5 times larger than the horizontal uncertainties, the vertical uncertainties as determined by the new method are only larger by a factor of 1.44 over all baseline lengths.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号