共查询到20条相似文献，搜索用时 31 毫秒
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. 相似文献
2.
利用小波理论讨论了奇异积分的计算问题,用实例说明了一维情况下计算的快速性和准确性。与Fourier变换法相比,列举了小波理论的优点,肯定了在局部重力场研究中小波理论用于实际的重要性。 相似文献
3.
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 相似文献
4.
With the advent of geodetic satellite missions mapping almost globally the Earth’s gravitational field, new methods and theoretical approaches have been developed and investigated to fully exploit the potential of their new observables. Besides estimating values of numerical coefficients in harmonic series of the gravitational potential, new applications emerged such as data validation and combination. In this contribution, new integral transformations are presented which transform principal components of the terrestrial deflection of the vertical onto disturbing satellite-to-satellite tracking and gradiometric data at altitude. Using spherical approximation, necessary integral kernel functions are derived in both spectral and closed forms. The behaviour of isotropic kernel functions is studied and the new integral transformations are tested in a closed-loop simulation using synthetic terrestrial and satellite data synthesized from a global gravitational model. New integral transformations can be used for data validation and combination purposes. 相似文献
5.
Geoid determination using one-step integration 总被引：1，自引：1，他引：0
P. Novák 《Journal of Geodesy》2003,77(3-4):193-206
A residual (high-frequency) gravimetric geoid is usually computed from geographically limited ground, sea and/or airborne gravimetric data. The mathematical model for its determination from ground gravity is based on the transformation of observed discrete values of gravity into gravity potential related to either the international ellipsoid or the geoid. The two reference surfaces are used depending on height information that accompanies ground gravity data: traditionally orthometric heights determined by geodetic levelling were used while GPS positioning nowadays allows for estimation of geodetic (ellipsoidal) heights. This transformation is usually performed in two steps: (1) observed values of gravity are downward continued to the ellipsoid or the geoid, and (2) gravity at the ellipsoid or the geoid is transformed into the corresponding potential. Each of these two steps represents the solution of one geodetic boundary-value problem of potential theory, namely the first and second or third problem. Thus two different geodetic boundary-value problems must be formulated and solved, which requires numerical evaluation of two surface integrals. In this contribution, a mathematical model in the form of a single Fredholm integral equation of the first kind is presented and numerically investigated. This model combines the solution of the first and second/third boundary-value problems and transforms ground gravity disturbances or anomalies into the harmonically downward continued disturbing potential at the ellipsoid or the geoid directly. Numerical tests show that the new approach offers an efficient and stable solution for the determination of the residual geoid from ground gravity data. 相似文献
6.
R. Lehmann 《Journal of Geodesy》1997,71(9):533-540
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 相似文献
7.
Numerical calculation of weakly singular surface integrals 总被引：1，自引：0，他引：1
Roland Klees 《Journal of Geodesy》1996,70(11):781-797
We consider semi-analytical and purely numerical integration methods for weakly singular integrals with point singularities on curved smooth surfaces. The methods can be applied to many practical computations in Geodesy, e.g. terrain corrections, Stokes' and Hotines' integral, surface potentials, and the solution of geodetic boundary value problems using integral equations. Current numerical integration techniques are reviewed. The most important semi-analytical and purely numerical techniques are described. Test calcualtions are done and the techniques are compared as regards accuracy and computational efficiency. Semi-analytical methods, which are based on some regularizing parameter transformations, are superior to purely numerical techniques. The best choice are modified polar coordinates defined in the parameter domain with the singularity as pole. Triangular coordinates show similar performance if carefully tuned. Extrapolation techniques and adaptive subdivision techniques behave poorly as regards accuracy and numerical efficiency. Standard integration techniques, which ignore the singularity, completely fail. 相似文献
8.
Optimized formulas for the gravitational field of a tesseroid 总被引：7，自引：3，他引：4
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. 相似文献
9.
The three-dimensional (3-D) resection problem is usually solved by first obtaining the distances connecting the unknown point P{X,Y,Z} to the known points Pi{Xi,Yi,Zi}i=1,2,3 through the solution of the three nonlinear Grunert equations and then using the obtained distances to determine the position {X,Y,Z} and the 3-D orientation parameters {,, }. Starting from the work of the German J. A. Grunert (1841), the Grunert equations have been solved in several substitutional steps and the desire as evidenced by several publications has been to reduce these number of steps. Similarly, the 3-D ranging step for position determination which follows the distance determination step involves the solution of three nonlinear ranging (`Bogenschnitt') equations solved in several substitution steps. It is illustrated how the algebraic technique of Groebner basis solves explicitly the nonlinear Grunert distance equations and the nonlinear 3-D ranging (`Bogenschnitt') equations in a single step once the equations have been converted into algebraic (polynomial) form. In particular, the algebraic tool of the Groebner basis provides symbolic solutions to the problem of 3-D resection. The various forward and backward substitution steps inherent in the classical closed-form solutions of the problem are avoided. Similar to the Gauss elimination technique in linear systems of equations, the Groebner basis eliminates several variables in a multivariate system of nonlinear equations in such a manner that the end product normally consists of a univariate polynomial whose roots can be determined by existing programs e.g. by using the roots command in Matlab.Acknowledgments.The first author wishes to acknowledge the support of JSPS (Japan Society of Promotion of Science) for the financial support that enabled the completion of the write-up of the paper at Kyoto University, Japan. The author is further grateful for the warm welcome and the good working atmosphere provided by his hosts Professors S. Takemoto and Y. Fukuda of the Department of Geophysics, Graduate School of Science, Kyoto University, Japan. 相似文献
10.
J. Feltens 《Journal of Geodesy》2009,83(2):129-137
The vector-based algorithm to transform Cartesian (X, Y, Z ) into geodetic coordinates (, λ, h) presented by Feltens (J Geod, 2007, doi:) has been extended for triaxial ellipsoids. The extended algorithm is again based on simple formulae and has successfully
been tested for the Earth and other celestial bodies and for a wide range of positive and negative ellipsoidal heights. 相似文献
11.
Solutions of the linearized geodetic boundary value problem for an ellipsoidal boundary to order e
3
The geodetic boundary value problem (GBVP) was originally formulated for the topographic surface of the Earth. It degenerates to an ellipsoidal problem, for example when topographic and downward continuation reductions have been applied. Although these ellipsoidal GBVPs possess a simpler structure than the original ones, they cannot be solved analytically, since the boundary condition still contains disturbing terms due to anisotropy, ellipticity and centrifugal components in the reference potential. Solutions of the so-called scalar-free version of the GBVP, upon which most recent practical calculations of geoidal and quasigeoidal heights are based, are considered. Starting at the linearized boundary condition and presupposing a normal field of Somigliana–Pizzetti type, the boundary condition described in spherical coordinates is expanded into a series with respect to the flattening f of the Earth. This series is truncated after the linear terms in f, and first-order solutions of the corresponding GBVP are developed in closed form on the basis of spherical integral formulae, modified by suitable reduction terms. Three alternative representations of the solution are discussed, implying corrections by adding a first-order non-spherical term to the solution, by reducing the boundary data, or by modifying the integration kernel. A numerically efficient procedure for the evaluation of ellipsoidal effects, in the case of the linearized scalar-free version of the GBVP, involving first-order ellipsoidal terms in the boundary condition, is derived, utilizing geopotential models such as EGM96. 相似文献
12.
A function having some properties of a wavelet and being harmonic around a given point in R
3 is defined, and three models showing the local relationships between the disturbing density, the disturbing potential and
the disturbing gravity are established by using the function as the kernel function of the integrals in the models. The local
relationship has two meanings. One is that we can evaluate with a high accuracy the integrals in the models by using mainly
high-accuracy and high-resolution data in a local area. The other is that we can obtain a stable solution with high resolution
when we invert the integrals in the models because of the rapid decrease of the kernel function of the integrals. As a result,
with these models we evaluate one quantity with high resolution, in a band limited by the maximum degree of a set of geopotential
coefficients or by the resolution (spacing) of the local data, from another quantity (or quantities) in a local area, and
the resulting solution is stable.
Received: 6 April 1998 / Accepted: 16 June 1999 相似文献
13.
为解决世界各国高程基准差异的问题,提出联合卫星重力场模型、地面重力数据、GNSS大地高、局部高程基准的正高或正常高,按大地边值问题法确定局部高程基准重力位差的方法。首先推导了利用传统地面"有偏"重力异常确定高程基准重力位差的方法;接着利用改化Stokes核函数削弱"有偏"重力异常的影响,并联合卫星重力场模型和地面"有偏"重力数据,得到独立于任何局部高程基准的重力水准面,以此来确定局部高程基准重力位差;最后利用GNSS+水准数据和重力大地水准面确定了美国高程基准与全球高程基准W0的重力位差为-4.82±0.05 m2s-2。 相似文献
14.
Peiliang Xu 《Journal of Geodesy》2003,77(5-6):247-256
The previous work of Xu on discrete nonlinear filtering is extended to continuous systems. The new results are summarized as follows: (1) a second-order unbiased prediction of the true state governed by a vector stochastic differential equation is worked out; (2) a set of coupled differential equations for a new truncated second-order nonlinear filter and its variance–covariance matrix are derived from the frequentist point of view. The new filter is proved to be unbiased to the second-order approximation; and, most importantly, (3) comparison of the new filtering and accuracy results with the literature on nonlinear filtering has indicated that more than 40 years of nonlinear filtering of continuous systems may have foundational problems.Acknowledgments.This work is supported by a Grant-in-Aid for Scientific Research (C13640422). The author thanks Prof J.A.R. Blais, Prof A. Dermanis and Prof B. Schaffrin for their constructive comments. 相似文献
15.
An operational algorithm for computation of terrain correction (or local gravity field modeling) based on application of closed-form solution of the Newton integral in terms of Cartesian coordinates in multi-cylindrical equal-area map projection of the reference ellipsoid is presented. Multi-cylindrical equal-area map projection of the reference ellipsoid has been derived and is described in detail for the first time. Ellipsoidal mass elements with various sizes on the surface of the reference ellipsoid are selected and the gravitational potential and vector of gravitational intensity (i.e. gravitational acceleration) of the mass elements are computed via numerical solution of the Newton integral in terms of geodetic coordinates {,,h}. Four base- edge points of the ellipsoidal mass elements are transformed into a multi-cylindrical equal-area map projection surface to build Cartesian mass elements by associating the height of the corresponding ellipsoidal mass elements to the transformed area elements. Using the closed-form solution of the Newton integral in terms of Cartesian coordinates, the gravitational potential and vector of gravitational intensity of the transformed Cartesian mass elements are computed and compared with those of the numerical solution of the Newton integral for the ellipsoidal mass elements in terms of geodetic coordinates. Numerical tests indicate that the difference between the two computations, i.e. numerical solution of the Newton integral for ellipsoidal mass elements in terms of geodetic coordinates and closed-form solution of the Newton integral in terms of Cartesian coordinates, in a multi-cylindrical equal-area map projection, is less than 1.6×10–8 m2/s2 for a mass element with a cross section area of 10×10 m and a height of 10,000 m. For a mass element with a cross section area of 1×1 km and a height of 10,000 m the difference is less than 1.5×10–4m2/s2. Since 1.5× 10–4 m2/s2 is equivalent to 1.5×10–5m in the vertical direction, it can be concluded that a method for terrain correction (or local gravity field modeling) based on closed-form solution of the Newton integral in terms of Cartesian coordinates of a multi-cylindrical equal-area map projection of the reference ellipsoid has been developed which has the accuracy of terrain correction (or local gravity field modeling) based on the Newton integral in terms of ellipsoidal coordinates.Acknowledgments. This research has been financially supported by the University of Tehran based on grant number 621/4/859. This support is gratefully acknowledged. The authors are also grateful for the comments and corrections made to the initial version of the paper by Dr. S. Petrovic from GFZ Potsdam and the other two anonymous reviewers. Their comments helped to improve the structure of the paper significantly. 相似文献
16.
Closed form solutions for transforming 3D Cartesian to geodetic coordinates reduce the problem to finding the real solutions
of the fourth degree latitude equation or variations of it. By using symbolic tools (Sturm–Habicht coefficients and subresultants
mainly) we study the methods (and polynomials) proposed by Vermeille and Borkowski to solve this problem. For Vermeille approach,
the region where it cannot be applied is completely characterized. For those points it is shown how to transform 3D Cartesian
to geodetic coordinates and a new method for solving Vermeille equation for those cases not yet covered is introduced. Concerning
Borkowski’s approach, the symbolic analysis produces a complete characterization of the singular cases (i.e. where multiple
roots appear). 相似文献
17.
Toshio Fukushima 《Journal of Geodesy》2018,92(2):123-130
In order to accelerate the spherical harmonic synthesis and/or analysis of arbitrary function on the unit sphere, we developed a pair of procedures to transform between a truncated spherical harmonic expansion and the corresponding two-dimensional Fourier series. First, we obtained an analytic expression of the sine/cosine series coefficient of the \(4 \pi \) fully normalized associated Legendre function in terms of the rectangle values of the Wigner d function. Then, we elaborated the existing method to transform the coefficients of the surface spherical harmonic expansion to those of the double Fourier series so as to be capable with arbitrary high degree and order. Next, we created a new method to transform inversely a given double Fourier series to the corresponding surface spherical harmonic expansion. The key of the new method is a couple of new recurrence formulas to compute the inverse transformation coefficients: a decreasing-order, fixed-degree, and fixed-wavenumber three-term formula for general terms, and an increasing-degree-and-order and fixed-wavenumber two-term formula for diagonal terms. Meanwhile, the two seed values are analytically prepared. Both of the forward and inverse transformation procedures are confirmed to be sufficiently accurate and applicable to an extremely high degree/order/wavenumber as \(2^{30}\,{\approx }\,10^9\). The developed procedures will be useful not only in the synthesis and analysis of the spherical harmonic expansion of arbitrary high degree and order, but also in the evaluation of the derivatives and integrals of the spherical harmonic expansion. 相似文献
18.
提出了一种用于Stokes积分和Hotine积分直接离散求和的快速算法。该算法将积分核表达为计算点纬度、流动点纬度和两点间经度差的函数,充分利用核函数的对称性,相同纬度的所有计算点只需计算一组核函数,计算次数远少于普通离散求和。基于EGM2008地球重力位模型的模拟实验表明,快速算法的计算效率远高于普通算法,有效解决了离散求和计算速度太慢的数值问题,且保留了球面积分的特性,可取代一维FFT用于计算Stokes积分和Hotine积分。 相似文献
19.
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. 相似文献
20.
利用非奇变换，将地形校正诸奇异积分转化为一组非奇异积分。理论分析和数值计算都表明，奇异积分非奇异后，可有效地提高地形校正中央区积分的精确度。 相似文献