共查询到20条相似文献，搜索用时 906 毫秒
1.
在评估重力场模型计算空间扰动引力精度时,对模型截断误差常采用阶方差方法。文中将6种经典的重力异常阶方差模型与现有超高阶重力场模型的阶方差进行比较,TSD模型与重力场模型的差值最小。根据重力异常阶方差模型TSD,文中分析不同高度、不同阶次利用重力场模型计算空中扰动引力时截断误差的影响。实验结果表明:36阶模型截断误差最大径向和水平方向分别为26.455 1mGal、25.946 3mGal;360阶模型截断误差最大径向和水平方向分别为9.969 0mGal、9.960 9 mGal;2160阶模型截断误差最大径向和水平方向分别为2.538 5 mGal、2.538 1mGal;2160阶模型计算空中扰动引力时,即使在低空附近,截断误差在2.5mGal以内,计算高度超过5km,截断误差可以忽略;超过400km的高度,都可以用36阶模型计算,截断误差在1mGal以内。 相似文献
2.
Y. Hagiwara 《Journal of Geodesy》1972,46(4):453-466
A general formula giving Molodenskii coefficientsQ
n of the truncation errors for the geoidal height is introduced in this paper. A relation betweenQ
n andq
n, Cook’s truncation function, is also obtained. Cook (1951) has treated the truncation errors for the deflection of the vertical
in the Vening Meinesz integration. Molodenskii et al. (1962) have also derived the truncation error formulas for the deflection
of the vertical. It is proved in this paper that these two formulas are equivalent. 相似文献
3.
Recurrence relations have been derived for truncation error coefficients of the extended Stokes' function and its partial
derivatives required in the computation of the disturbing gravity vector at any elevation above the earth's surface. The corresponding
formulae, the example of values of the truncation error coefficients for H=30.1 km and ψ0=30∘ and the estimations of truncation error are given in this article.
Received: 26 January 1996 / Accepted: 11 June 1997 相似文献
4.
Improved convergence rates for the truncation error in gravimetric geoid determination 总被引：2，自引：2，他引：0
When Stokes's integral is used over a spherical cap to compute a gravimetric estimate of the geoid, a truncation error results
due to the neglect of gravity data over the remainder of the Earth. Associated with the truncation error is an error kernel
defined over these two complementary regions. An important observation is that the rate of decay of the coefficients of the
series expansion for the truncation error in terms of Legendre polynomials is determined by the smoothness properties of the
error kernel. Previously published deterministic modifications of Stokes's integration kernel involve either a discontinuity
in the error kernel or its first derivative at the spherical cap radius. These kernels are generalised and extended by constructing
error kernels whose derivatives at the spherical cap radius are continuous up to an arbitrary order. This construction is
achieved by smoothly continuing the error kernel function into the spherical cap using a suitable degree polynomial. Accordingly,
an improved rate of convergence of the spectral series representation of the truncation error is obtained.
Received: 21 April 1998 / Accepted: 4 October 1999 相似文献
5.
Local geoid determination from airborne vector gravimetry 总被引：3，自引：2，他引：1
Methods are illustrated to compute the local geoid using the vertical and horizontal components of the gravity disturbance vector derived from an airborne GPS/inertial navigation system. The data were collected by the University of Calgary in a test area of the Canadian Rocky Mountains and consist of multiple parallel tracks and two crossing tracks of accelerometer and gyro measurements, as well as precise GPS positions. Both the boundary-value problem approach (Hotines integral) and the profiling approach (line integral) were applied to compute the disturbing potential at flight altitude. Cross-over adjustments with minimal control were investigated and utilized to remove error biases and trends in the estimated gravity disturbance components. Final estimation of the geoid from the vertical gravity disturbance included downward continuation of the disturbing potential with correction for intervening terrain masses. A comparison of geoid estimates to the Canadian Geoid 2000 (CGG2000) yielded an average standard deviation per track of 14 cm if they were derived from the vertical gravity disturbance (minimally controlled with a cross-over adjustment), and 10 cm if derived from the horizontal components (minimally controlled in part with a simulated cross-over adjustment). Downward continuation improved the estimates slightly by decreasing the average standard deviation by about 0.5 cm. The application of a wave correlation filter to both types of geoid estimates yielded significant improvement by decreasing the average standard deviation per track to 7.6 cm. 相似文献
6.
Brian Farelly 《Journal of Geodesy》1991,65(2):92-101
Six sources of error in the use of Fourier methods for the conversion of geoid heights to gravity anomalies are considered. The errors due to spherical approximation are unimportant. The errors due to approximations in Stokes' integral may be eliminated by use of the gravity coating rather than the gravity anomaly. The chord-to-arc error and the truncation error may be reduced by using a reference field. Tapering of the edges of the measurement window reduces the truncation error. The long-wavelength components of the high degree spherical harmonics cause small offsets in the resulting gravity anomalies. The errors due to the plane approximation can be reduced by appropriate choice of map projection and area of integration. 相似文献
7.
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. 相似文献
8.
Computation of spherical harmonic coefficients and their error estimates using least-squares collocation 总被引：4，自引：0，他引：4
C. C. Tscherning 《Journal of Geodesy》2001,75(1):12-18
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 相似文献
9.
分别采用基于梯度、基于泊松积分和基于快速傅里叶变换（FFT）的地面重力向上延拓方案,并提出交叉检验方法估计地面重力数据误差及其空中误差传播,对毛乌素测区GT-2A航空重力测量系统采集的空中测线数据进行外符合精度评价。对比结果表明：地面重力格网插值误差和代表性误差对空中点的影响达到0.66~0.92 mGal（1 Gal=1×10-2 m/s2）,航空重力数据误差估计必须扣除这一影响;基于泊松积分和基于FFT的地面重力向上延拓方法能够客观评价航空重力观测值的外符合精度,二者表现相当;扣除地面重力误差影响后,在包含残余边界效应的情况下,毛乌素测区GT-2A航空重力空中测线重力扰动的外符合精度优于1.42 mGal。 相似文献
10.
Meissel-Stokes核函数应用于区域大地水准面分析 总被引：1，自引：0，他引：1
为提高区域大地水准面计算精度,基于EGM2008地球重力场位系数模型分析Meissel-Stokes核函数、截断误差系数以及截断误差。选取实验区,采用移去-恢复法评价Meissel-Stokes核函数计算大地水准面的精度。结果表明:Meissel-Stokes核函数及其截断误差系数收敛速度快;截断误差小且稳定。在积分半径不易扩展的情况下,应用Meissel-Stokes核函数计算区域大地水准面,比标准Stokes计算大地水准面精度略高。 相似文献
11.
Richard H. Rapp 《Journal of Geodesy》1980,54(2):149-163
A gravimetric geoid computed using different techniques has been compared to a geoid derived from Geos-3 altimeter data in
two 30°×30° areas: one in the Tonga Trench area and one in the Indian Ocean. The specific techniques used were the usual Stokes
integration (using 1°×1° mean anomalies) with the Molodenskii truncation procedure; a modified Stokes integration with a modified
truncation method; and computations using three sets of potential coefficients including one complete to degree 180. In the
Tonga Trench area the standard deviation of the difference between the modified Stokes’ procedure and the altimeter geoid
was ±1.1 m while in the Indian Ocean area the difference was ±0.6 m. Similar results were found from the 180×180 potential
coefficient field. However, the differences in using the usual Stokes integration procedure were about a factor of two greater
as was predicted from an error analysis.
We conclude that there is good agreement at the ±1 m level between the two types of geoids. In addition, systematic differences
are at the half-meter level. The modified Stokes procedure clearly is superior to the usual Stokes method although the 180×180
solution is of comparable accuracy with the computational effort six times less than the integration procedures. 相似文献
12.
G. L. Ramillien 《Journal of Geodesy》2017,91(8):881-895
A radial integration of spherical mass elements (i.e. tesseroids) is presented for evaluating the six components of the second-order gravity gradient (i.e. second derivatives of the Newtonian mass integral for the gravitational potential) created by an uneven spherical topography consisting of juxtaposed vertical prisms. The method uses Legendre polynomial series and takes elastic compensation of the topography by the Earth’s surface into account. The speed of computation of the polynomial series increases logically with the observing altitude from the source of anomaly. Such a forward modelling can be easily applied for reduction of observed gravity gradient anomalies by the effects of any spherical interface of density. An iterative least-squares inversion of measured gravity gradient coefficients is also proposed to estimate a regional set of juxtaposed topographic heights. Several tests of recovery have been made by considering simulated gradients created by idealistic conical and irregular Great Meteor seamount topographies, and for varying satellite altitudes and testing different levels of uncertainty. In the case of gravity gradients measured at a GOCE-type altitude of \(\sim \)300 km, the search converges down to a stable but smooth topography after 10–15 iterations, while the final root-mean-square error is \(\sim \)100 m that represents only 2 % of the seamount amplitude. This recovery error decreases with the altitude of the gravity gradient observations by revealing more topographic details in the region of survey. 相似文献
13.
M. K. Paul 《Journal of Geodesy》1973,47(4):413-425
Neglecting distant zones in the computation of geoidal height using Stokes' formula gives rise to some truncation error. This
truncation error is expressible as a weighted summation of the zonal harmonic components of the gravity anomaly. Making use
of the well-known properties of Legendre polynomials, a compact method of computing these theoretical coefficients has been
developed in this paper. 相似文献
14.
Lars E. Sjöberg 《Journal of Geodesy》1989,63(2):213-221
Integral formulas are derived for the determination of geopotential coefficients from gravity anomalies and gravity disturbances
over the surface of the Earth. First order topographic corrections to spherical formulas are presented. In addition new integral
formulas are derived for the determination of the external gravity field from surface gravity.
Taking advantage of modern satellite positioning techniques, it is suggested that, in general, the external gravity field
as well as individual coefficients are better determined from gravity disturbances than from gravity anomalies. 相似文献
15.
Gravity gradient modeling using gravity and DEM 总被引：2，自引：0，他引：2
A model of the gravity gradient tensor at aircraft altitude is developed from the combination of ground gravity anomaly data
and a digital elevation model. The gravity data are processed according to various operational solutions to the boundary-value
problem (numerical integration of Stokes’ integral, radial-basis splines, and least-squares collocation). The terrain elevation
data are used to reduce free-air anomalies to the geoid and to compute a corresponding indirect effect on the gradients at
altitude. We compare the various modeled gradients to airborne gradiometric data and find differences of the order of 10–20 E
(SD) for all gradient tensor elements. Our analysis of these differences leads to a conclusion that their source may be primarily
measurement error in these particular gradient data. We have thus demonstrated the procedures and the utility of combining
ground gravity and elevation data to validate airborne gradiometer systems. 相似文献
16.
基于物理大地测量边值问题的解,利用一阶边界算子定义,推导重力异常Δg、单层密度μ、大地水准面高N,垂线偏差ε、扰动重力δg等扰动场元的解。利用球谐函数的正交特性,通过对核函数的算子运算,可以得到上述扰动场元的有关逆变换公式。相对经典物理大地测量公式应用的边界面条件,笔者将含有因子r的对应扰动场元反演关系的公式称为广义积分公式。针对常用的重力异常Δg、大地水准面高N,垂线偏差ε、扰动重力δg计算,重点分析它们之间的变换关系,给出利用某个选定扰动场元计算其他扰动场元的广义积分公式。同时,通过对积分边界面的讨论,分析经典公式与广义积分公式的差异和联系。最后,给出所有外部扰动场元与核函数映射的关系表。 相似文献
17.
A formula for computing the gravity disturbance from the second radial derivative of the disturbing potential 总被引：6，自引：0，他引：6
J. Li 《Journal of Geodesy》2002,76(4):226-231
A formula for computing the gravity disturbance and gravity anomaly from the second radial derivative of the disturbing potential
is derived in detail using the basic differential equation with spherical approximation in physical geodesy and the modified
Poisson integral formula. The derived integral in the space domain, expressed by a spherical geometric quantity, is then converted
to a convolution form in the local planar rectangular coordinate system tangent to the geoid at the computing point, and the
corresponding spectral formulae of 1-D FFT and 2-D FFT are presented for numerical computation.
Received: 27 December 2000 / Accepted: 3 September 2001 相似文献
18.
Christopher Jekeli 《Journal of Geodesy》1980,54(2):137-147
Errors are considered in the outer zone contribution to oceanic undulation differences as obtained from a set of potential
coefficients complete to degree 180. It is assumed that the gravity data of the inner zone (a spherical cap), consisting of
either gravity anomalies or gravity disturbances, has negligible error. This implies that error estimates of the total undulation
difference are analyzed. If the potential coefficients are derived from a global field of 1°×1° mean anomalies accurate to
εΔg=10 mgal, then for a cap radius of 10°, the undulation difference error (for separations between 100 km and 2000 km) ranges
from 13 cm to 55 cm in the gravity anomaly case and from 6 cm to 36 cm in the gravity disturbance case. If εΔg is reduced to 1 mgal, these errors in both cases are less than 10 cm. In the absence of a spherical cap, both cases yield
identical error estimates: about 68 cm if εΔg=1 mgal (for most separations) and ranging from 93 cm to 160 cm if εΔg=10 mgal. Introducing a perfect 30-degree reference field, the latter errors are reduced to about 110 cm for most separations. 相似文献
19.
20.