引力场梯度张量的非奇异公式推导

传统的引力场梯度计算公式在两极附近存在奇异性,需要换用其他的非奇异计算公式。从奇异性产生的原因入手,并结合勒让德函数的有关性质,推导了一组新的计算公式。实际计算验证了该公式的正确性和有效性。     

Determination of the potential of homogeneous polyhedral bodies using line integrals

For the determination of the potential of irregular inhomogeneous bodies they can be decomposed into (polyhedral) parts of homogeneous density. Efficient formulas for the computation of the gravitational potential (and its first and second derivatives) of homogeneous polyhedral bodies are presented. They are obtained using a transformation of the volume integral into line integrals. The most important property of the solution is that all ten quantities under consideration (potential, 3 components of the gravitation vector, 6 components of the tensor of the second derivatives) can be represented by using only two different line integrals. Furthermore, all coordinate transformations needed in the evaluation are chosen in such a way that they do not appear in the final result. The consequence, favorable for efficient programming, is that the same transcendental expressions along each edge of the polyhedron are needed for all ten quantities; even the same linear combinations of them for individual surfaces are appearing in different formulas. The expressions obtained are probably the simplest possible, which is also reflected in the fact that for the special case of a right rectangular prism they may easily be specialized to the usual well-known formulas. Received 28 Juni 1994; Accepted 13 September 1996     

Comparative study on two methods for calculating the gravitational potential of a prism

The determination of the gravitational potential of a prism plays an important role in physical geodesy and geophysics. However, there are few literatures that provide accurate approaches for determining the gravitational potential of a prism. Discrete element method can be used to determine the gravitational potential of a prism, and can approximate the true gravitational potential values with sufficient accuracy (the smaller each element is, the more accurate the result is). Although Nagy’s approach provided a closed expression, one does not know whether it is valid, due to the fact that this approach has not been confirmed in literatures. In this paper, a study on the comparison of Nagy’s approach with discrete element method is presented. The results show that Nagy’s formulas for determining the gravitational potential of a prism are valid in the domain both inside and outside the prism.     

Comparative study on two methods for calculating the gravitational potential of a prism

The determination of the gravitational potential of a prism plays an important role in physical geodesy and geophysics. However, there are few literatures that provide accurate approaches for determining the gravitational potential of a prism. Discrete element method can be used to determine the gravitational potential of a prism, and can approximate the true gravitational potential values with sufficient accuracy (the smaller each element is, the more accurate the result is). Although Nagy’s approach provi...     

Optimized formulas for the gravitational field of a tesseroid

Various tasks in geodesy, geophysics, and related geosciences require precise information on the impact of mass distributions on gravity field-related quantities, such as the gravitational potential and its partial derivatives. Using forward modeling based on Newton’s integral, mass distributions are generally decomposed into regular elementary bodies. In classical approaches, prisms or point mass approximations are mostly utilized. Considering the effect of the sphericity of the Earth, alternative mass modeling methods based on tesseroid bodies (spherical prisms) should be taken into account, particularly in regional and global applications. Expressions for the gravitational field of a point mass are relatively simple when formulated in Cartesian coordinates. In the case of integrating over a tesseroid volume bounded by geocentric spherical coordinates, it will be shown that it is also beneficial to represent the integral kernel in terms of Cartesian coordinates. This considerably simplifies the determination of the tesseroid’s potential derivatives in comparison with previously published methodologies that make use of integral kernels expressed in spherical coordinates. Based on this idea, optimized formulas for the gravitational potential of a homogeneous tesseroid and its derivatives up to second-order are elaborated in this paper. These new formulas do not suffer from the polar singularity of the spherical coordinate system and can, therefore, be evaluated for any position on the globe. Since integrals over tesseroid volumes cannot be solved analytically, the numerical evaluation is achieved by means of expanding the integral kernel in a Taylor series with fourth-order error in the spatial coordinates of the integration point. As the structure of the Cartesian integral kernel is substantially simplified, Taylor coefficients can be represented in a compact and computationally attractive form. Thus, the use of the optimized tesseroid formulas particularly benefits from a significant decrease in computation time by about 45 % compared to previously used algorithms. In order to show the computational efficiency and to validate the mathematical derivations, the new tesseroid formulas are applied to two realistic numerical experiments and are compared to previously published tesseroid methods and the conventional prism approach.     

A comparison of the tesseroid,prism and point-mass approaches for mass reductions in gravity field modelling

The calculation of topographic (and iso- static) reductions is one of the most time-consuming operations in gravity field modelling. For this calculation, the topographic surface of the Earth is often divided with respect to geographical or map-grid lines, and the topographic heights are averaged over the respective grid elements. The bodies bounded by surfaces of constant (ellipsoidal) heights and geographical grid lines are denoted as tesseroids. Usually these ellipsoidal (or spherical) tesseroids are replaced by “equivalent” vertical rectangular prisms of the same mass. This approximation is motivated by the fact that the volume integrals for the calculation of the potential and its derivatives can be exactly solved for rectangular prisms, but not for the tesseroids. In this paper, an approximate solution of the spherical tesseroid integrals is provided based on series expansions including third-order terms. By choosing the geometrical centre of the tesseroid as the Taylor expansion point, the number of non-vanishing series terms can be greatly reduced. The zero-order term is equivalent to the point-mass formula. Test computations show the high numerical efficiency of the tesseroid method versus the prism approach, both regarding computation time and accuracy. Since the approximation errors due to the truncation of the Taylor series decrease very quickly with increasing distance of the tesseroid from the computation point, only the elements in the direct vicinity of the computation point have to be separately evaluated, e.g. by the prism formulas. The results are also compared with the point-mass formula. Further potential refinements of the tesseroid approach, such as considering ellipsoidal tesseroids, are indicated.     

近区地形直接与间接影响的棱柱模型算法

基于Helmert第二压缩法进行边值解算时需要计算地形压缩对重力的直接影响和对（似）大地水准面的间接影响。计算近区直接、间接影响的传统积分算法仍是二重积分形式。该算法以网格中心点处的积分核作为网格积分核的平均值的计算模式在一定程度上引入了近似误差。另外,直接、间接影响的传统积分算法在中央区存在奇异性,需单独计算中央网格地形影响,因而增加了计算的复杂性。为此,本文推导了近区地形直接、间接影响的棱柱模型公式,一方面提高了地形影响的计算精度;另一方面中央区不存在奇异性,从而简化了计算过程。为避免棱柱模型存在的平面近似误差,可使用顾及地球曲率的棱柱模型算法计算地形影响。最后通过试验得出结论,在（似）大地水准面精度要求较高的应用中,应尽量使用顾及地球曲率的棱柱模型算法计算地形影响。     

Gravity recovery using COSMIC GPS data: application of orbital perturbation theory

 COSMIC is a joint Taiwan–US mission to study the atmosphere using the Global Positioning System (GPS) occultation technique. Improved formulas are developed for the radial, along-track, and cross-track perturbations, which are more accurate than the commonly used order-zero formulas. The formulas are used to simulate gravity recovery using the geodetic GPS data of COSMIC in the operational phase. Results show that the EGM96 model can be improved up to degree 26 using 1 year of COSMIC data. TOPEX/POSEIDON altimeter data are used to derive a temporal gravity variation. COSMIC cannot reproduce this gravity variation perfectly because of data noise and orbital configuration, but the recovered field clearly shows the gravity signature due to mass movement in an El Ni?o. Received: 3 March 2000 / Accepted: 10 November 2000     

基于逆Vening-Meinesz公式的测高重力中央区效应精密计算

为提高利用逆Vening-Meinesz公式反演测高重力中央区效应的精度,视中央区为矩形域,将垂线偏差分量表示成双二次多项式插值形式,引入非奇异变换,推导出了重力异常的计算公式。以低纬度区域2′×2′的垂线偏差实际数据为背景场进行了计算,结果表明,当中央区包含4个网格时,传统公式与推导出的重力异常计算公式误差的最大值大于1 mGal。推导出的公式可为高精度测高重力中央区效应的计算提供理论依据。     

An Evaluation of Fractal Methods for Characterizing Image Complexity

Previously, we developed an integrated software package called ICAMS (Image Characterization and Modeling System) to provide specialized spatial analytical functions for interpreting remote sensing data. This paper evaluates three fractal dimension measurement methods that have been implemented in ICAMS: isarithm, variogram, and a modified version of triangular prism. To provide insights into how the fractal methods compare with conventional spatial techniques in measuring landscape complexity, the performance of two spatial autocorrelation methods, Moran's I and Geary's C, is also evaluated. Results from analyzing 25 simulated surfaces having known fractal dimensions show that both the isarithm and triangular prism methods can accurately measure a range of fractal surfaces. The triangular prism method is most accurate at estimating the fractal dimension of surfaces having higher spatial complexity, but it is sensitive to contrast stretching. The variogram method is a comparatively poor estimator for all surfaces, particularly those with high fractal dimensions. As with the fractal techniques, spatial autocorrelation techniques have been found to be useful for measuring complex images, but not images with low dimensionality. Fractal measurement methods, as well as spatial autocorrelation techniques, can be applied directly to unclassified images and could serve as a tool for change detection and data mining.     

Integral formulas for geopotential coefficient determination from gravity anomalies versus gravity disturbances

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.     

利用余弦变换的重力梯度Parker正演方法及其应用

重力测量数据存在地形数据产生的高频分量的影响,高精度地形数据正演重力梯度也能较好地反映重力局部高频特征。为获得高精度重力梯度数据,实现基准梯度数据库精确快速构建,研究了利用数字高程模型正演重力梯度的频率域快速计算方法,推导出基于余弦变换的Parker正演重力梯度理论公式。数值实验结果表明,余弦变换频率域正演方法平均绝对误差可达到0.5E左右精度要求,与傅里叶变换正演方法相比误差可减小3dB左右,与棱柱法等空间域正演方法相比,该方法计算规模小,速度优势明显。     

On the ellipsoidal corrections to gravity anomalies computed using the inverse Stokes integral

The formulas of the ellipsoidal corrections to the gravity anomalies computed using the inverse Stokes integral are derived. The corrections are given in the integral formulas and expanded in the spherical harmonics series. If a coefficient model such as the OSU91A is given, the corrections can be easily computed. Received: 19 August 1996 / Accepted: 28 September 1998     

The Eighth International Cartographic Conference and Fifth General Assembly of the International Cartographic Association: A Report

The utility of a stepped statistical surface for either choropleth style or absolute quantity themes depends on how accurately readers can estimate prism heights, and whether prism volumes affect those estimates. Testing of a stepped surface through comparison of state pairs without benefit of legend revealed that most readers respond to prism heights, not volumes. The association of values with heights is consistent for a variety of data areas and for three different themes. Even for an absolute quantity, which logically can be represented by volume, the surface, scaled by height, conveyed magnitudes with as much accuracy as scaled circles used for the same data. Altogether, the results show a stepped surface scaled by height to be a versatile device, and suggest that a surface representing values by volumes might be misinterpreted.     

Fourier geoid determination with irregular data

The fast Fourier transform (FFT) and, recently, the fast Hartley transform (FHT) have been extensively used by geodesists for efficient geoid determination. For this kind of efficiency, data must be given on a regular grid and, consequently, a pre-processing step of interpolation is required when only point measurements are available. This paper presents a way of computing a grid of geoid undulations N without explicitly gridding the data. The method is applicable to all FFT or FHT techniques of geoid or terrain effects determination, and it works with planar as well as spherical formulas. This method can be used not only for, e.g., computing a grid of undulations from irregular gravity anomalies g but it also lends itself to other applications, such as the gridding of gravity anomalies and, since the contribution of each data point is computed individually, the update of N- or g-grids as soon as new point measurements become available. In the case that there are grid cells which contain no measurements, the results of gravity interpolation or geoid estimation can be drastically improved by incorporating into the procedure a frequency-domain interpolating function. In addition to numerical results obtained using a few simple interpolating functions, the paper presents briefly the mathematical formulas for recovering missing grid values and for transforming values from one grid to another which might be rotated and/or scaled with respect to the first one. The geodetic problems where these techniques may find applications are pointed out throughout the paper.Presented at theIAG General Meeting, Beijing, P.R. China, Aug. 6–13, 1993     

海洋重力似大地水准面与区域测高似大地水准面的拟合问题

研究了将陆地重力似大地水准面与GPS水；住似大地水准面拟合的处理方法推广到海洋的问题．首先从理论上证明了当存在海面地形．则海洋大地水准面与似大地水准面不重合．导出了在海洋上大地水；住面差距与高程异常之间差值的公式．由此给出了求定平均海面相对于区域高程基准的正常高以及测高似大地水准面的计算公式。由于测高平均海面与GPS大地高有相近的精度．提出了将海洋重力似大地水准面与区域测高似大地水准面拟合的处理方法．并利用当前最新的海面地形模型和测高平均海面模型做了数值估计。     

A comparison of methods for the inversion of airborne gravity data

Four integral-based methods for the inversion of gravity disturbances, derived from airborne gravity measurements, into the disturbing potential on the Bjerhammar sphere and the Earths surface are investigated and compared with least-squares (LS) collocation. The performance of the methods is numerically investigated using noise-free and noisy observations, which have been generated using a synthetic gravity field model. It is found that advanced interpolation of gravity disturbances at the nodes of higher-order numerical integration formulas significantly improves the performance of the integral-based methods. This is preferable to the commonly used one-point composed Newton–Cotes integration formulas, which intrinsically imply a piecewise constant interpolation over a patch centered at the observation point. It is shown that the investigated methods behave similarly for noise-free observations, but differently for noisy observations. The best results in terms of root-mean-square (RMS) height-anomaly errors are obtained when the gravity disturbances are first downward continued (inverse Poisson integral) and then transformed into potential values (Hotine integral). The latter has a strong smoothing effect, which damps high-frequency errors inherent in the downward-continued gravity disturbances. An integral method based on the single-layer representation of the disturbing potential shows a similar performance. This representation has the advantage that it can be used directly on surfaces with non-spherical geometry, whereas classical integral-based methods require an additional step if gravity field functionals have to be computed on non-spherical geometries. It is shown that defining the single-layer density on the Bjerhammar sphere gives results with the same quality as obtained when using the Earths topography as support for the single-layer density. A comparison of the four integral-based methods with LS collocation shows that the latter method performs slightly better in terms of RMS height-anomaly errors.     

Analytical computation of gravity effects for polyhedral bodies

On the basis of recent analytical results we derive new formulas for computing the gravity effects of polyhedral bodies which are expressed solely as function of the coordinates of the vertices of the relevant faces. We thus prove that such formulas exhibit no singularity whenever the position of the observation point is not aligned with an edge of a face. In the opposite case, the contribution of the edge to the potential to its first-order derivative and to the diagonal entries of the second-order derivative is deemed to be zero on the basis of some claims which still require a rigorous mathematical proof. In contrast with a common statement in the literature, it is proved that only the off-diagonal entries of the second-order derivative of the potential do exhibit a noneliminable singularity when the observation point is aligned with an edge of a face. The analytical provisions on the range of validity of the derived formulas have been fully confirmed by the Matlab $^{\textregistered }$ program which has been coded and thoroughly tested by computing the gravity effects induced by real asteroids at arbitrarily placed observation points.     

Effects of topographic–isostatic masses on gravitational functionals at the Earth’s surface and at airborne and satellite altitudes

Topographic–isostatic masses represent an important source of gravity field information, especially in the high-frequency band, even if the detailed mass-density distribution inside the topographic masses is unknown. If this information is used within a remove-restore procedure, then the instability problems in downward continuation of gravity observations from aircraft or satellite altitudes can be reduced. In this article, integral formulae are derived for determination of gravitational effects of topographic–isostatic masses on the first- and second-order derivatives of the gravitational potential for three topographic–isostatic models. The application of these formulas is useful for airborne gravimetry/gradiometry and satellite gravity gradiometry. The formulas are presented in spherical approximation by separating the 3D integration in an analytical integration in the radial direction and 2D integration over the mean sphere. Therefore, spherical volume elements can be considered as being approximated by mass-lines located at the centre of the discretization compartments (the mass of the tesseroid is condensed mathematically along its vertical axis). The errors of this approximation are investigated for the second-order derivatives of the topographic–isostatic gravitational potential in the vicinity of the Earth’s surface. The formulas are then applied to various scenarios of airborne gravimetry/gradiometry and satellite gradiometry. The components of the gravitational vector at aircraft altitudes of 4 and 10 km have been determined, as well as the gravitational tensor components at a satellite altitude of 250 km envisaged for the forthcoming GOCE (gravity field and steady-state ocean-circulation explorer) mission. The numerical computations are based on digital elevation models with a 5-arc-minute resolution for satellite gravity gradiometry and 1-arc-minute resolution for airborne gravity/gradiometry.