首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 593 毫秒
1.
Green's function for the boundary-value problem of Stokes's type with ellipsoidal corrections in the boundary condition for anomalous gravity is constructed in a closed form. The `spherical-ellipsoidal' Stokes function describing the effect of two ellipsoidal correcting terms occurring in the boundary condition for anomalous gravity is expressed in O(e 2 0)-approximation as a finite sum of elementary functions analytically representing the behaviour of the integration kernel at the singular point ψ=0. We show that the `spherical-ellipsoidal' Stokes function has only a logarithmic singularity in the vicinity of its singular point. The constructed Green function enables us to avoid applying an iterative approach to solve Stokes's boundary-value problem with ellipsoidal correction terms involved in the boundary condition for anomalous gravity. A new Green-function approach is more convenient from the numerical point of view since the solution of the boundary-value problem is determined in one step by computing a Stokes-type integral. The question of the convergence of an iterative scheme recommended so far to solve this boundary-value problem is thus irrelevant. Received: 5 June 1997 / Accepted: 20 February 1998  相似文献   

2.
The ellipsoidal Stokes problem is one of the basic boundary-value problems for the Laplace equation which arises in physical geodesy. Up to now, geodecists have treated this and related problems with high-order series expansions of spherical and spheroidal (ellipsoidal) harmonics. In view of increasing computational power and modern numerical techniques, boundary element methods have become more and more popular in the last decade. This article demonstrates and investigates the nullfield method for a class of Robin boundary-value problems. The ellipsoidal Stokes problem belongs to this class. An integral equation formulation is achieved, and existence and uniqueness conditions are attained in view of the Fredholm alternative. Explicit expressions for the eigenvalues and eigenfunctions for the boundary integral operator are provided. Received: 22 October 1996 / Accepted: 4 August 1997  相似文献   

3.
确定似大地水准面的Hotine-Helmert边值解算模型   总被引：1，自引：1，他引：0  
马健  魏子卿  任红飞 《测绘学报》2019,48(2):153-160
空间大地测量技术的发展使大地高的观测成为可能,从而为第二大地边值问题的研究带来了新的机遇,本文对基于Helmert第二压缩法的第二边值问题(简称为Hotine-Helmert边值问题)展开研究。首先介绍了地形直接、间接影响的定义与算法,然后推导了Hotine-Helmert边值问题的解算模型。Hotine-Helmert边值理论无须计算地形压缩对重力的次要间接影响,因而较Stokes-Helmert边值理论更简单。此外,文中引入了一种低阶修正的Hotine截断核函数,该核函数较传统的截断核函数能有效地改善似大地水准面的解算精度。为了验证本文构建的Hotine-Helmert边值解算模型的有效性和实用性,本文将EIGEN-6C4模型的前360阶作为参考模型,利用Hotine-Helmert边值解算模型构建了我国中部地区6°×4°范围、1.5′×1.5′分辨率的重力似大地水准面,其精度达到±4.8 cm。  相似文献   

4.
J. Li 《Journal of Geodesy》2005,79(1-3):64-70
Integral formulas are derived which can be used to convert the second-order radial gradient of the disturbing potential, as boundary values, into the disturbing potential, gravity anomaly and the deflection of the vertical. The derivations are based on the fundamental differential equation as the boundary condition in Stokes’s boundary-value problem and the modified Poisson integral formula in which the zero and first-degree spherical harmonics are excluded. The rigorous kernel functions, corresponding to the integral operators, are developed by the methods of integration.  相似文献   

5.
We propose a methodology for local gravity field modelling from gravity data using spherical radial basis functions. The methodology comprises two steps: in step 1, gravity data (gravity anomalies and/or gravity disturbances) are used to estimate the disturbing potential using least-squares techniques. The latter is represented as a linear combination of spherical radial basis functions (SRBFs). A data-adaptive strategy is used to select the optimal number, location, and depths of the SRBFs using generalized cross validation. Variance component estimation is used to determine the optimal regularization parameter and to properly weight the different data sets. In the second step, the gravimetric height anomalies are combined with observed differences between global positioning system (GPS) ellipsoidal heights and normal heights. The data combination is written as the solution of a Cauchy boundary-value problem for the Laplace equation. This allows removal of the non-uniqueness of the problem of local gravity field modelling from terrestrial gravity data. At the same time, existing systematic distortions in the gravimetric and geometric height anomalies are also absorbed into the combination. The approach is used to compute a height reference surface for the Netherlands. The solution is compared with NLGEO2004, the official Dutch height reference surface, which has been computed using the same data but a Stokes-based approach with kernel modification and a geometric six-parameter “corrector surface” to fit the gravimetric solution to the GPS-levelling points. A direct comparison of both height reference surfaces shows an RMS difference of 0.6 cm; the maximum difference is 2.1 cm. A test at independent GPS-levelling control points, confirms that our solution is in no way inferior to NLGEO2004.  相似文献   

6.
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.  相似文献   

7.
Solving the geodetic boundary-value problem (GBVP) for the precise determination of the geoid requires proper use of the fundamental equation of physical geodesy as the boundary condition given on the geoid. The Stokes formula and kernel are the result of spherical approximation of this fundamental equation, which is a violation of the proper relation between the observed quantity (gravity anomaly) and the sought function (geoid). The violation is interpreted here as the improper formulation of the boundary condition, which implies the spherical Stokes kernel to be in error compared with the proper kernel of integral transformation. To remedy this error, two correction kernels to the Stokes kernel were derived: the first in both closed and spectral forms and the second only in spectral form. Contributions from the first correction kernel to the geoid across the globe were [−0.867 m, +1.002 m] in the low-frequency domain implied by the GRIM4-S4 purely satellite-derived geopotential model. It is a few centimeters, on average, in the high-frequency domain with some exceptions of a few meters in places of high topographical relief and sizable geological features in accordance with the EGM96 combined geopotential model. The contributions from the second correction kernel to the geoid are [−0.259 m, +0.217 m] and [−0.024 m, +0.023 m] in the low- and high-frequency domains, respectively.  相似文献   

8.
 This paper generalizes the Stokes formula from the spherical boundary surface to the ellipsoidal boundary surface. The resulting solution (ellipsoidal geoidal height), consisting of two parts, i.e. the spherical geoidal height N 0 evaluated from Stokes's formula and the ellipsoidal correction N 1, makes the relative geoidal height error decrease from O(e 2) to O(e 4), which can be neglected for most practical purposes. The ellipsoidal correction N 1 is expressed as a sum of an integral about the spherical geoidal height N 0 and a simple analytical function of N 0 and the first three geopotential coefficients. The kernel function in the integral has the same degree of singularity at the origin as the original Stokes function. A brief comparison among this and other solutions shows that this solution is more effective than the solutions of Molodensky et al. and Moritz and, when the evaluation of the ellipsoidal correction N 1 is done in an area where the spherical geoidal height N 0 has already been evaluated, it is also more effective than the solution of Martinec and Grafarend. Received: 27 January 1999 / Accepted: 4 October 1999  相似文献   

9.
    
When the values of gravity anomalies are given at the geoid, Ag can be calculated at altitude by application of Poisson’s integral theorem. The process requires integration of Δg multiplied by the Poisson kernel function over the entire globe. It is common practice to add to the kernel function terms that will ensure removal of any zeroth and first order components of Δg that may be present. The effects of trancating the integration at the boundary of a spherical cap of earth central half angle ψo have been analyzed using an adaptation of Molodenskii’s procedure. The extension process without removal terms retains the correct effects of inaccuracies in the constant term of the gravity reference model used in the definition of Δg. Furthermore, the effects of ignoring remote zones or unmapped areas in the integration process are very much smaller for the extension without removal terms than for the commonly used formula with removal terms. For these reasons the Poisson vertical extension process without removal terms is to be preferred over the extension with the zeroth order term removal. Truncation of this process at the point recommended for the Stokes integration, namely, the first zero crossing of the Stokes kernel function, leaves negligible truncation errors.  相似文献   

10.
On the accurate numerical evaluation of geodetic convolution integrals   总被引：3，自引：2，他引：1  
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.  相似文献   

11.
Geodetic surface integrals play an important role in the numerical solution of geodetic boundary-value problems. In many cases they can be evaluated using fast methods in the frequency domain (FFT). However, this is not possible in general, because the domain of integration may be non-trivial (as is the surface of the Earth), the kernel function may not be of convolution type, or the data distribution may be heterogeneous. Therefore, fast evaluation strategies are also required in the space domain. They are more difficult to design because only one property is left where a more or less fast evaluation strategy can be built upon: the potential type of the kernel function. Consequently, the idea is not to replace well-established frequency domain techniques, but to supplement them. Our approach to this problem goes in two directions: (1) we use advanced cubature methods where the integration nodes automatically densify in the vicinity of the evaluation points; (2) we use powerful computer hardware, namely MIMD computers with distributed memory. This enables us to evaluate geodetic surface integrals of any practical complexity in reasonable time and accuracy. This is shown in a numerical example. Received: 7 May 1996 / Accepted:17 March 1997  相似文献   

12.
A new theory for high-resolution regional geoid computation without applying Stokess formula is presented. Operationally, it uses various types of gravity functionals, namely data of type gravity potential (gravimetric leveling), vertical derivatives of the gravity potential (modulus of gravity intensity from gravimetric surveys), horizontal derivatives of the gravity potential (vertical deflections from astrogeodetic observations) or higher-order derivatives such as gravity gradients. Its algorithmic version can be described as follows: (1) Remove the effect of a very high degree/order potential reference field at the point of measurement (POM), in particular GPS positioned, either on the Earths surface or in its external space. (2) Remove the centrifugal potential and its higher-order derivatives at the POM. (3) Remove the gravitational field of topographic masses (terrain effect) in a zone of influence of radius r. A proper choice of such a radius of influence is 2r=4×104 km/n, where n is the highest degree of the harmonic expansion. (cf. Nyquist frequency). This third remove step aims at generating a harmonic gravitational field outside a reference ellipsoid, which is an equipotential surface of a reference potential field. (4) The residual gravitational functionals are downward continued to the reference ellipsoid by means of the inverse solution of the ellipsoidal Dirichlet boundary-value problem based upon the ellipsoidal Abel–Poisson kernel. As a discretized integral equation of the first kind, downward continuation is Phillips–Tikhonov regularized by an optimal choice of the regularization factor. (5) Restore the effect of a very high degree/order potential reference field at the corresponding point to the POM on the reference ellipsoid. (6) Restore the centrifugal potential and its higher-order derivatives at the ellipsoidal corresponding point to the POM. (7) Restore the gravitational field of topographic masses ( terrain effect) at the ellipsoidal corresponding point to the POM. (8) Convert the gravitational potential on the reference ellipsoid to geoidal undulations by means of the ellipsoidal Bruns formula. A large-scale application of the new concept of geoid computation is made for the Iran geoid. According to the numerical investigations based on the applied methodology, a new geoid solution for Iran with an accuracy of a few centimeters is achieved.Acknowledgments. The project of high-resolution geoid computation of Iran has been support by National Cartographic Center (NCC) of Iran. The University of Tehran, via grant number 621/3/602, supported the computation of a global geoid solution for Iran. Their support is gratefully acknowledged. A. Ardalan would like to thank Mr. Y. Hatam, and Mr. K. Ghazavi from NCC and Mr. M. Sharifi, Mr. A. Safari, and Mr. M. Motagh from the University of Tehran for their support in data gathering and computations. The authors would like to thank the comments and corrections made by the four reviewers and the editor of the paper, Professor Will Featherstone. Their comments helped us to correct the mistakes and improve the paper.  相似文献   

13.
Efficient numerical computation of integrals defined on closed surfaces in ℝ3 with non-integrable point singularities that arise in physical geodesy is discussed. The method is based on the use of polar coordinates and the definition of integrals with non-integrable point singularities as Hadamard finite part integrals. First the behavior of singular integrals under smooth parameter transformations is studied, and then it is shown how they can be reduced to absolutely integrable functions over domains in ℝ2. The correction terms that usually arise if the substitution rule is formally applied, in contrast to absolutely integrable functions, are calculated. It is shown how to compute the regularized integrals efficiently, and, numerical efforts for various orders of singularity are compared. Finally, efficient numerical integration methods are discussed for integrals of functions that are defined as singular integrals, a task that typically arises in Galerkin boundary element methods. Received: 15 April 1997 / Accepted: 7 May 1998  相似文献   

14.
The upward-downward continuation of a harmonic function like the gravitational potential is conventionally based on the direct-inverse Abel-Poisson integral with respect to a sphere of reference. Here we aim at an error estimation of the “planar approximation” of the Abel-Poisson kernel, which is often used due to its convolution form. Such a convolution form is a prerequisite to applying fast Fourier transformation techniques. By means of an oblique azimuthal map projection / projection onto the local tangent plane at an evaluation point of the reference sphere of type “equiareal” we arrive at a rigorous transformation of the Abel-Poisson kernel/Abel-Poisson integral in a convolution form. As soon as we expand the “equiareal” Abel-Poisson kernel/Abel-Poisson integral we gain the “planar approximation”. The differences between the exact Abel-Poisson kernel of type “equiareal” and the “planar approximation” are plotted and tabulated. Six configurations are studied in detail in order to document the error budget, which varies from 0.1% for points at a spherical height H=10km above the terrestrial reference sphere up to 98% for points at a spherical height H = 6.3×106km. Received: 18 March 1997 / Accepted: 19 January 1998  相似文献   

15.
A new local existence and uniqueness theorem is obtained for the scalar geodetic boundary-value problem in spherical coordinates. The regularities H α and H 1+α are assumed for the boundary data g (gravity) and v (gravitational potential) respectively. Received: 27 July 1998 / Accepted: 19 April 1999  相似文献   

16.
The target of the spheroidal Gauss–Listing geoid determination is presented as a solution of the spheroidal fixed–free two-boundary value problem based on a spheroidal Bruns' transformation (“spheroidal Bruns' formula”). The nonlinear spheroidal Bruns' transform (nonlinear spheroidal Bruns' formula), the spheroidal fixed part and the spheroidal free part of the two-boundary value problem are derived. Four different spheroidal gravity models are treated, in particular to determine whether they pass the test to fit to the postulate of a level ellipsoidal gravity field, namely of Somigliana–Pizzetti type. Received: 4 May 1999 / Accepted: 21 May 1999  相似文献   

17.
借助以地心参考椭球面为边界面的第二大地边值问题的理论,基于Helmert空间的Neumann边值条件,给定Helmert扰动位的椭球解表达式,并详细推导第二类勒让德函数及其导数的递推关系、Helmert扰动位函数的椭球积分解以及类椭球Hotine积分核函数的实用计算公式,便于后续椭球域第二大地边值问题的实际研究。  相似文献   

18.
The differential equations which generate a general conformal mapping of a two-dimensional Riemann manifold found by Korn and Lichtenstein are reviewed. The Korn–Lichtenstein equations subject to the integrability conditions of type vectorial Laplace–Beltrami equations are solved for the geometry of an ellipsoid of revolution (International Reference Ellipsoid), specifically in the function space of bivariate polynomials in terms of surface normal ellipsoidal longitude and ellipsoidal latitude. The related coefficient constraints are collected in two corollaries. We present the constraints to the general solution of the Korn–Lichtenstein equations which directly generates Gau?–Krüger conformal coordinates as well as the Universal Transverse Mercator Projection (UTM) avoiding any intermediate isometric coordinate representation. Namely, the equidistant mapping of a meridian of reference generates the constraints in question. Finally, the detailed computation of the solution is given in terms of bivariate polynomials up to degree five with coefficients listed in closed form. Received: 3 June 1997 / Accepted: 17 November 1997  相似文献   

19.
P. Holota 《Journal of Geodesy》1997,71(10):640-651
In this paper the linear gravimetric boundary-value problem is discussed in the sense of the so-called weak solution. For this purpose a Sobolev weight space was constructed for an unbounded domain representing the exterior of the Earth and quantitative estimates were deduced for the trace theorem and equivalent norms. In the generalized formulation of the problem a special decomposition of the Laplace operator was used to express the oblique derivative in the boundary condition which has to be met by the solution. The relation to the classical formulation was also shown. The main result concerns the coerciveness (ellipticity) of a bilinear form associated with the problem under consideration. The Lax-Milgram theorem was used to decide about the existence, uniqueness and stability of the weak solution of the problem. Finally, a clear geometrical interpretation was found for a constant in the coerciveness inequality, and the convergence of approximation solutions constructed by means of the Galerkin method was proved. Received: 21 June 1996 / Accepted: 14 April 1997  相似文献   

20.
The solution of the linear Molodensky problem by analytical continuation to point level is numerically the most convenient of all the theoretically equivalent solutions. It is obtained by successively applying the same integral operator and it does not depend explicitly on the terrain inclination. However, its dependence on the computation point restricts somehow the computational efficiency. The use of the Fourier transform for the evaluation of the integral operator in planar approximation lessens significantly the burden of computations. Using this spectral approach, the problem has been reformulated and solved in the frequency domain. Moreover, it is shown that the solution can be easily split into two steps: (a) “downward” continuation to sea level, which is independent of the computation point, and (b) “upward” continuation from sea to point level, using the values computed at sea level. Such a treatment not only simplifies the formulas and increases the numerical efficiency but also clarifies the physical interpretation and the theoretical equivalence of the continuation solution with respect to the other solution types. Numerical tests have been performed to investigate which terms in the Molodensky series are of significance for geoid and deflection computations. The practical difficulty of differences in the grid spacings of gravity and height data has been overcome by frequency domain interpolation. Presented at theXIX IUGG General Assembly, Vancouver, B.C., August 9–22, 1987.  相似文献   

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

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