首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The AUSGeoid09 model of the Australian Height Datum   总被引：8，自引：6，他引：2  
AUSGeoid09 is the new Australia-wide gravimetric quasigeoid model that has been a posteriori fitted to the Australian Height Datum (AHD) so as to provide a product that is practically useful for the more direct determination of AHD heights from Global Navigation Satellite Systems (GNSS). This approach is necessary because the AHD is predominantly a third-order vertical datum that contains a ~1 m north-south tilt and ~0.5 m regional distortions with respect to the quasigeoid, meaning that GNSS-gravimetric-quasigeoid and AHD heights are inconsistent. Because the AHD remains the official vertical datum in Australia, it is necessary to provide GNSS users with effective means of recovering AHD heights. The gravimetric component of the quasigeoid model was computed using a hybrid of the remove-compute-restore technique with a degree-40 deterministically modified kernel over a one-degree spherical cap, which is superior to the remove-compute-restore technique alone in Australia (with or without a cap). This is because the modified kernel and cap combine to filter long-wavelength errors from the terrestrial gravity anomalies. The zero-tide EGM2008 global gravitational model to degree 2,190 was used as the reference field. Other input data are ~1.4 million land gravity anomalies from Geoscience Australia, 1′ × 1′ DNSC2008GRA altimeter-derived gravity anomalies offshore, the 9′′ × 9′′ GEODATA-DEM9S Australian digital elevation model, and a readjustment of Australian National Levelling Network (ANLN) constrained to the CARS2006 mean dynamic ocean topography model. To determine the numerical integration parameters for the modified kernel, the gravimetric component of AUSGeoid09 was compared with 911 GNSS-observed ellipsoidal heights at benchmarks. The standard deviation of fit to the GNSS-AHD heights is ±222 mm, which dropped to ±134 mm for the readjusted GNSS-ANLN heights showing that careful consideration now needs to be given to the quality of the levelling data used to assess gravimetric quasigeoid models. The publicly released version of AUSGeoid09 also includes a geometric component that models the difference between the gravimetric quasigeoid and the zero surface of the AHD at 6,794 benchmarks. This a posteriori fitting used least-squares collocation (LSC) in cross-validation mode to determine a correlation length of 75 km for the analytical covariance function, whereas the noise was taken from the estimated standard deviation of the GNSS ellipsoidal heights. After this LSC surface fitting, the standard deviation of fit reduced to ±30 mm, one-third of which is attributable to the uncertainty in the GNSS ellipsoidal heights.  相似文献   

2.
基于最小二乘配置法向下延拓航空重力的过程中,由于协方差矩阵严重病态,影响延拓结果的稳定性和精度。针对这一问题,提出了航空重力向下延拓的最小二乘配置Tikhonov正则化法。基于全球协方差函数模型建立航空重力数据与地面重力数据的协方差关系,引入基于广义交叉验证法,选择正则化参数的Tikhonov正则化法改善协方差矩阵的病态性,抑制观测噪声对延拓结果的放大影响。基于EGM2008重力场模型,设计了山区、丘陵和海域3种不同地形区域的航空重力数据向下延拓的仿真实验,实验结果验证了该方法的有效性。  相似文献   

3.
On the basis of gravity field model （EIGEN_CG01C）, together with multi-altimeter data, the improved deflection of the vertical gridded in 2＇×2＇ in China marginal sea and gridded in 5＇×5＇ in the global sea was determined by using the weighted method of along-track least squares, and the accuracy is better than 1.2^# in China marginal sea. As for the quality of the deflection of the vertical, it meets the challenge for the gravity field of high resolution and accuracy, it shows that, compared with the shipboard gravimetry in the sea, the accuracy of the gravity anomalies computed with the marine deflection of the vertical by inverse Vening-Meinesz formula is 7.75 m.s ^-2.  相似文献   

4.
 Equations expressing the covariances between spherical harmonic coefficients and linear functionals applied on the anomalous gravity potential, T, are derived. The functionals are the evaluation functionals, and those associated with first- and second-order derivatives of T. These equations form the basis for the prediction of spherical harmonic coefficients using least-squares collocation (LSC). The equations were implemented in the GRAVSOFT program GEOCOL. Initially, tests using EGM96 were performed using global and regional sets of geoid heights, gravity anomalies and second-order vertical gravity gradients at ground level and at altitude. The global tests confirm that coefficients may be estimated consistently using LSC while the error estimates are much too large for the lower-order coefficients. The validity of an error estimate calculated using LSC with an isotropic covariance function is based on a hypothesis that the coefficients of a specific degree all belong to the same normal distribution. However, the coefficients of lower degree do not fulfil this, and this seems to be the reason for the too-pessimistic error estimates. In order to test this the coefficients of EGM96 were perturbed, so that the pertubations for a specific degree all belonged to a normal distribution with the variance equal to the mean error variance of the coefficients. The pertubations were used to generate residual geoid heights, gravity anomalies and second-order vertical gravity gradients. These data were then used to calculate estimates of the perturbed coefficients as well as error estimates of the quantities, which now have a very good agreement with the errors computed from the simulated observed minus calculated coefficients. Tests with regionally distributed data showed that long-wavelength information is lost, but also that it seems to be recovered for specific coefficients depending on where the data are located. Received: 3 February 2000 / Accepted: 23 October 2000  相似文献   

5.
 Two numerical techniques are used in recent regional high-frequency geoid computations in Canada: discrete numerical integration and fast Fourier transform. These two techniques have been tested for their numerical accuracy using a synthetic gravity field. The synthetic field was generated by artificially extending the EGM96 spherical harmonic coefficients to degree 2160, which is commensurate with the regular 5′ geographical grid used in Canada. This field was used to generate self-consistent sets of synthetic gravity anomalies and synthetic geoid heights with different degree variance spectra, which were used as control on the numerical geoid computation techniques. Both the discrete integration and the fast Fourier transform were applied within a 6∘ spherical cap centered at each computation point. The effect of the gravity data outside the spherical cap was computed using the spheroidal Molodenskij approach. Comparisons of these geoid solutions with the synthetic geoid heights over western Canada indicate that the high-frequency geoid can be computed with an accuracy of approximately 1 cm using the modified Stokes technique, with discrete numerical integration giving a slightly, though not significantly, better result than fast Fourier transform. Received: 2 November 1999 / Accepted: 11 July 2000  相似文献   

6.
Mean gravity anomalies, deflections of the vertical, and a geopotential model complete to degree and order180 are combined in order to determine geoidal heights in the area bounded by [34°≦ϕ≤42°, 18°≦λ≦28°]. Moreover, employing point gravity anomalies simultaneously with the above data, an attempt is made to predict deflections of the vertical in the same area. The method used in the computations is least squares collocation. Using empirical covariance functions for the data, the suitable errors for the different sources of observations, and the optimum cap radius around each point of evaluation, an accuracy better than±0.60m for geoidal heights and±1″.5 for deflections of the vertical is obtained taking into account existing systematic effects. This accuracy refers to the comparison between observed and predicted values.  相似文献   

7.
LSC法（最小二乘配置法）因能融合不同种类重力观测数据确定大地水准面的特性而受到广泛关注,但由于协方差矩阵存在病态性,微小的观测误差将被协方差矩阵的小奇异值放大,导致计算的配置结果不稳定且精度偏低。本文提出Tikhonov_LSC法,即在LSC法中引入Tikhonov正则化算法,基于GCV法选择协方差矩阵的正则化参数,利用正则化参数修正协方差矩阵的小奇异值,以抑制其对观测误差的放大影响。基于Tikhonov_LSC法计算大地水准面,能有效提高其稳定性和精度。通过以EGM2008重力场模型分别计算山区、丘陵和海域重力异常作为基础数据确定相应区域大地水准面的实验,验证了该方法的有效性。  相似文献   

8.
The least squares collocation algorithm for estimating gravity anomalies from geodetic data is shown to be an application of the well known regression equations which provide the mean and covariance of a random vector (gravity anomalies) given a realization of a correlated random vector (geodetic data). It is also shown that the collocation solution for gravity anomalies is equivalent to the conventional least-squares-Stokes' function solution when the conventional solution utilizes properly weighted zero a priori estimates. The mathematical and physical assumptions underlying the least squares collocation estimator are described.  相似文献   

9.
The accuracy of the gravity field approximation depends on the amount of the available data and their distribution as well as on the variation of the gravity field. The variation of the gravity field in the Greek mainland, which is the test area in this study, is very high (the variance of point free air gravity anomalies is 3191.5mgal 2). Among well known reductions used to smooth the gravity field, the complete isostatic reduction causes the best possible smoothing, however remain strong local anomalies which disturb the homogeneity of the gravity field in this area. The prediction of free air gravity anomalies using least squares collocation and regional covariance function is obtained within a ±4 ... ±19mgal accuracy depending on the local peculiarities of the free air gravity field. By taking into account the topography and its isostatic compensation with the usual remove-restore technique, the accuracy of the prediction mentioned obove was increased by about a factor of 4 and the prediction results become quite insensitive to the covariance function used (local or regional). But when predicting geoidal heights, in spite of using the smoothed field, the prediction results remain still depend on the covariance function used in such a way that differences up to about 50cm/100km result between relative geoidal heights computed with regional or local covariance functions.  相似文献   

10.
A synthetic Earth for use in geodesy   总被引：1，自引：0，他引：1  
 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  相似文献   

11.
This study demonstrates that in mountainous areas the use of residual terrain model (RTM) data significantly improves the accuracy of vertical deflections obtained from high-degree spherical harmonic synthesis. The new Earth gravitational model EGM2008 is used to compute vertical deflections up to a spherical harmonic degree of 2,160. RTM data can be constructed as difference between high-resolution Shuttle Radar Topography Mission (SRTM) elevation data and the terrain model DTM2006.0 (a spherical harmonic terrain model that complements EGM2008) providing the long-wavelength reference surface. Because these RTM elevations imply most of the gravity field signal beyond spherical harmonic degree of 2,160, they can be used to augment EGM2008 vertical deflection predictions in the very high spherical harmonic degrees. In two mountainous test areas—the German and the Swiss Alps—the combined use of EGM2008 and RTM data was successfully tested at 223 stations with high-precision astrogeodetic vertical deflections from recent zenith camera observations (accuracy of about 0.1 arc seconds) available. The comparison of EGM2008 vertical deflections with the ground-truth astrogeodetic observations shows root mean square (RMS) values (from differences) of 3.5 arc seconds for ξ and 3.2 arc seconds for η, respectively. Using a combination of EGM2008 and RTM data for the prediction of vertical deflections considerably reduces the RMS values to the level of 0.8 arc seconds for both vertical deflection components, which is a significant improvement of about 75%. Density anomalies of the real topography with respect to the residual model topography are one factor limiting the accuracy of the approach. The proposed technique for vertical deflection predictions is based on three publicly available data sets: (1) EGM2008, (2) DTM2006.0 and (3) SRTM elevation data. This allows replication of the approach for improving the accuracy of EGM2008 vertical deflection predictions in regions with a rough topography or for improved validation of EGM2008 and future high-degree spherical harmonic models by means of independent ground truth data.  相似文献   

12.
An inverse Poisson integral technique has been used to determine a gravity field on the geoid which, when continued by analytic free space methods to the topographic surface, agrees with the observed field. The computation is performed in three stages, each stage refining the previous solution using data at progressively increasing resolution (1o×1o, 5′×5′, 5/8′×5/8′) from a decreasing area of integration. Reduction corrections are computed at 5/8′×5/8′ granularity by differencing the geoidal and surface values, smoothed by low-pass filtering and sub-sampled at 5′ intervals. This paper discusses 1o×1o averages of the reduction corrections thus obtained for 172 1o×1o squares in western North America. The 1o×1o mean reduction corrections are predominantly positive, varying from −3 to +15mgal, with values in excess of 5mgal for 26 squares. Their mean andrms values are +2.4 and 3.6mgal respectively and they correlate well with the mean terrain corrections as predicted byPellinen in 1962. The mean andrms contributions from the three stages of computation are: 1o×1o stage +0.15 and 0.7mgal; 5′×5′ stage +1.0 and 1.6mgal; and 5/8′×5/8′ stage +1.3 and 1.8mgal. These results reflect a tendency for the contributions to become larger and more systematically positive as the wavelengths involved become shorter. The results are discussed in terms of two mechanisms; the first is a tendency for the absolute values of both positive and negative anomalies to become larger when continued downwards and, the second, a non-linear rectification, due to the correlation between gravity anomaly and topographic height, which results in the values continued to a level surface being systematically more positive than those on the topography.  相似文献   

13.
A global geopotential model, like EGM2008, is not capable of representing the high-frequency components of Earth’s gravity field. This is known as the omission error. In mountainous terrain, omission errors in EGM2008, even when expanded to degree 2,190, may reach amplitudes of 10 cm and more for height anomalies. The present paper proposes the utilisation of high-resolution residual terrain model (RTM) data for computing estimates of the omission error in rugged terrain. RTM elevations may be constructed as the difference between the SRTM (Shuttle Radar Topography Mission) elevation model and the DTM2006.0 spherical harmonic topographic expansion. Numerical tests, carried out in the German Alps with a precise gravimetric quasigeoid model (GCG05) and GPS/levelling data as references, demonstrate that RTM-based omission error estimates improve EGM2008 height anomaly differences by 10 cm in many cases. The comparisons of EGM2008-only height anomalies and the GCG05 model showed 3.7 cm standard deviation after a bias-fit. Applying RTM omission error estimates to EGM2008 reduces the standard deviation to 1.9 cm which equates to a significant improvement rate of 47%. Using GPS/levelling data strongly corroborates these findings with an improvement rate of 49%. The proposed RTM approach may be of practical value to improve quasigeoid determination in mountainous areas without sufficient regional gravity data coverage, e.g., in parts of Asia, South America or Africa. As a further application, RTM omission error estimates will allow refined validation of global gravity field models like EGM2008 from GPS/levelling data.  相似文献   

14.
航空重力向下延拓是病态问题,而广义补偿最小二乘法可以很好地克服病态性。研究了基于广义补偿最小二乘法的逆Poisson积分的航空重力向下延拓模型,并设计实验方案,将EGM2008地球位模型计算的重力异常作为仿真实验的数据,分别用最小二乘、Tikhonov正则化、广义补偿最小二乘3种方法求解,对其精度及仿真效果进行比较。结果表明,广义补偿最小二乘方法精度高,对仿真效果有显著提高。  相似文献   

15.
如果最小二乘拟合推估法被应用在重力异常、高程异常等的内插中 ,当观测值中含有粗差时 ,由此拟合的协方差函数就不能精确表征其统计性质。本文先从协方差函数的拟合过程入手 ,通过分析传统的协方差函数拟合法的无抗差性 ,提出了协方差函数的抗差拟合法  相似文献   

16.
协方差函数的抗差拟合   总被引：3，自引：2，他引：3  
刘念 《测绘科学》2001,26(3):25-28
如果最小二乘拟合推估法被应用在重力异常、高程异常等的内插中 ,当观测值中含有粗差时 ,由此拟合的协方差函数就不能精确表征其统计性质。本文先从协方差函数的拟合过程入手 ,通过分析传统的协方差函数拟合法的无抗差性 ,提出了协方差函数的抗差拟合法  相似文献   

17.
Summary Basically two different evaluation methods are available to compute geoid heights from residual gravity anomalies in the inner zone: numerical integration and least squares collocation.If collocation is not applied to a global gravity data set, as is usually the case in practice, its result will not be equal to the numerical integration result. However, the cross covariance function between geoid heights and gravity anomalies can be adapted such that the geoid contribution is computed only from a small gravity area up to a certain distance o from the computation point. Using this modification, identical results are obtained as from numerical integration.Applying this modification makes the results less dependent on the covariance function used. The difference between numerical integration and collocation is mainly caused by the implicitly extrapolated residual gravity anomaly values, outside the original data area. This extrapolated signal depends very much on the covariance function used, while the interpolated values within the original data area depend much less on it.As a sort of by-product, this modified collocation formula also leads to a new combination technique of numerical integration and collocation, in which the optimizing practical properties of both methods are fully exploited.Numerical examples are added as illustration.  相似文献   

18.
High-resolution EIGEN6C4 and EGM2008 Bouguer gravity data of 2190 degree spherical harmonic over the Singhbhum-Orissa Craton, India, have been generated from the International Centre for Global Earth Models. The Bouguer gravity anomaly difference maps of (i) in situ and EIGEN6C4, (ii) in situ and EGM2008 and iii) EIGEN6C4 and EGM2008 of the study area are compared. It reveals that EIGEN6C4 has lesser systematic error than EGM2008. However, from different profile plots of Bouguer gravity, east–west horizontal derivative and north–south horizontal derivative anomalies of the in situ, EIGEN6C4 and EGM2008, it is observed that most of the signatures of lithounits and geological structural elements are delineated very well by EGM2008 and match 94–98% with those of EIGEN6C4. Further, the Bouguer gravity, east–west horizontal derivative and north–south horizontal derivative anomalies of EGM2008 data over the study area have been used effectively for identifying various lithounits and geological structural elements.  相似文献   

19.
A detailed gravimetric geoid in the North Atlantic Ocean, named DGGNA-77, has been computed, based on a satellite and gravimetry derived earth potential model (consisting in spherical harmonic coefficients up to degree and order 30) and mean free air surface gravity anomalies (35180 1°×1° mean values and 245000 4′×4′ mean values). The long wavelength undulations were computed from the spherical harmonics of the reference potential model and the details were obtained by integrating the residual gravity anomalies through the Stokes formula: from 0 to 5° with the 4′×4′ data, and from 5° to 20° with the 1°×1° data. For computer time reasons the final grid was computed with half a degree spacing only. This grid extends from the Gulf of Mexico to the European and African coasts. Comparisons have been made with Geos 3 altimetry derived geoid heights and with the 5′×5′ gravimetric geoid derived byMarsh andChang [8] in the northwestern part of the Atlantic Ocean, which show a good agreement in most places apart from some tilts which porbably come from the satellite orbit recovery.  相似文献   

20.
In support of requirements for the U.S. Air Force Cambridge Research Laboratories, gravity anomalies have been upward continued to several elevations in different areas of the United States. One area was 340 to 400 N in latitude and 960 to 1030 W in longitude, generally called the Oklahoma area. The computations proceeded from 26, 032 point anomalies to the prediction of mean anomalies in 14, 704, 2.5′×2.5′ blocks and 9,284, 5′×5′ blocks. These anomalies were upward continued along 28 profiles at 5′ intervals for every 30′ in latitude and longitude. These anomalies at elevations were meaned in various patterns to form mean 30′×30″, 10×10, 50×50 blocks. Comparisons were then made to the corresponding ground values. The results of these comparisons lead to practical recommendations on the arrangement of flight profiles in airborne gravimetry.  相似文献   

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

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