共查询到19条相似文献，搜索用时 109 毫秒
1.
应用文献 [1 ]推导出的球谐系数与椭球谐系数的转换关系 ,给出了椭球界面下Neumann边值问题的积分解 相似文献
2.
应用文献[1]推导出的球谐系数与椭球谐系数的转换关系,给出了椭球界面下Neumann边值问题的积分解. 相似文献
3.
本文结合椭球域的特点导出了三类椭球域边值问题的准格林函数解。这些解的主项分别是椭球域泊松积分，纽曼积分和司托克斯积分；次项是Ｏ（ｅ＾２）量级的非谐和级数改正项，可由位系数模型按相应的公式算得。 相似文献
4.
本文采用物理方法讨论不同椭球变换的问题，论证了重力位低阶球谐系数与大地坐标系转换参数的关系，说明几何法转换模型（即大地坐标微分公式）不过是取至二阶次的谐函微分公式，为了提高转换精度，应该采用较高阶次的物理法微分公式（即球面正交多项式）作为不同系统的转换模型。 相似文献
5.
扰动重力梯度的非奇异表示 总被引：5，自引：0，他引：5
在局部指北坐标系中用地心球坐标来表示扰动重力梯度张量，当计算点趋近于两极时，由于Legendre函数的一阶和二阶导数以及分母上所含余纬的正弦函数，将导致扰动重力梯度张量的计算出现无穷大。因此，本文引入了Legendre函数的一阶和二阶导数以及 无奇异性的计算公式，并且进一步推导了 无奇异性的计算公式。在将Legendre函数的一阶和二阶导数以及 、 无奇异性的计算公式代入到扰动重力梯度张量各分量的求解中时，又充分考虑了m等于0，1，2以及其它量时的复杂情况，建立了扰动重力梯度张量各分量无奇异性的详细计算模型。通过模拟实验表明，本文所建立的详细计算模型不仅能够完全满足当前卫星重力梯度张量计算的精度要求，而且模型稳定、可靠、易于编程实现。 相似文献
6.
针对10 800阶地形球谐系数模型构建的计算精度、计算稳定性及超大规模运算量等问题,本文首先通过比较分析矩形离散积分方法、Gauss-Legendre积分方法和Driscoll/Healy积分方法,验证了Driscoll/Healy积分方法具有更高的计算精度。然后,给出了改进的Belikov公式,将完全正常化缔和勒让德函数在全球范围内以优于10-12精度高效、稳定递推至10 800阶;提出了联合FFT和基于OpenMP的多核并行方法求解超高阶球谐系数的优化策略,显著提高了球谐系数模型构建的计算效率。最后,利用Earth2014＿TBI与全球海底地形模型STO＿IEU2020格网数据建立了全球10 800阶地形球谐系数模型sph.10 800＿IEU,总体精度与Earth2014＿TBI2014.shc相当,在试验海域的相对精度略优于模型Earth2014＿TBI2014.shc。 相似文献
7.
8.
借助以地心参考椭球面为边界面的第二大地边值问题的理论,基于Helmert空间的Neumann边值条件,给定Helmert扰动位的椭球解表达式,并详细推导第二类勒让德函数及其导数的递推关系、Helmert扰动位函数的椭球积分解以及类椭球Hotine积分核函数的实用计算公式,便于后续椭球域第二大地边值问题的实际研究。 相似文献
9.
深入研究了利用卫星重力梯度数据确定地球重力场模型的球谐分析方法，导出了由重力梯度张量的球函数展开系数确定扰动位球谐系数的实用解算模型。模拟试算结果验证了本文算法的有效性和实用性。 相似文献
10.
针对不同模型方法在低空扰动引力计算中的适用问题,该文选取我国某山地区域,分析比较球谐位系数模型法、点质量模型法和Stokes积分法在低空不同高度的扰动引力计算精度及效率;并且分析误差来源和个别改进办法。结果表明:点质量模型计算低空扰动引力精度较高,且速度最快;去奇异点的Stokes积分法可以解决低空积分时数值溢出的问题,但精度较低;球谐位系数模型法原理简单,但计算速度最慢。 相似文献
11.
This paper is devoted to the spherical and spheroidal harmonic expansion of the gravitational potential of the topographic masses in the most rigorous way. Such an expansion can be used to compute gravimetric topographic effects for geodetic and geophysical applications. It can also be used to augment a global gravity model to a much higher resolution of the gravitational potential of the topography. A formulation for a spherical harmonic expansion is developed without the spherical approximation. Then, formulas for the spheroidal harmonic expansion are derived. For the latter, Legendre’s functions of the first and second kinds with imaginary variable are expanded in Laurent series. They are then scaled into two real power series of the second eccentricity of the reference ellipsoid. Using these series, formulas for computing the spheroidal harmonic coefficients are reduced to surface harmonic analysis. Two numerical examples are presented. The first is a spherical harmonic expansion to degree and order 2700 by taking advantage of existing software. It demonstrates that rigorous spherical harmonic expansion is possible, but the computed potential on the geoid shows noticeable error pattern at Polar Regions due to the downward continuation from the bounding sphere to the geoid. The second numerical example is the spheroidal expansion to degree and order 180 for the exterior space. The power series of the second eccentricity of the reference ellipsoid is truncated at the eighth order leading to omission errors of 25 nm (RMS) for land areas, with extreme values around 0.5 mm to geoid height. The results show that the ellipsoidal correction is 1.65 m (RMS) over land areas, with maximum value of 13.19 m in the Andes. It shows also that the correction resembles the topography closely, implying that the ellipsoidal correction is rich in all frequencies of the gravity field and not only long wavelength as it is commonly assumed. 相似文献
12.
Gravity data observed on or reduced to the ellipsoid are preferably represented using ellipsoidal harmonics instead of spherical harmonics. Ellipsoidal harmonics, however, are difficult to use in practice because the computation of the associated Legendre functions of the second kind that occur in the ellipsoidal harmonic expansions is not straightforward. Jekeli’s renormalization simplifies the computation of the associated Legendre functions. We extended the direct computation of these functions—as well as that of their ratio—up to the second derivatives and minimized the number of required recurrences by a suitable hypergeometric transformation. Compared with the original Jekeli’s renormalization the associated Legendre differential equation is fulfilled up to much higher degrees and orders for our optimized recurrences. The derived functions were tested by comparing functionals of the gravitational potential computed with both ellipsoidal and spherical harmonic syntheses. As an input, the high resolution global gravity field model EGM2008 was used. The relative agreement we found between the results of ellipsoidal and spherical syntheses is 10?14, 10?12 and 10?8 for the potential and its first and second derivatives, respectively. Using the original renormalization, this agreement is 10?12, 10?8 and 10?5, respectively. In addition, our optimized recurrences require less computation time as the number of required terms for the hypergeometric functions is less. 相似文献
13.
Computations of Fourier coefficients and related integrals of the associated Legendre functions with a new method along with their application to spherical harmonics analysis and synthesis are presented. The method incorporates a stable three-step recursion equation that can be processed separately for each colatitudinal Fourier wavenumber. Recursion equations for the zonal and sectorial modes are derived in explicit single-term formulas to provide accurate initial condition. Stable computations of the Fourier coefficients as well as the integrals needed for the projection of Legendre functions are demonstrated for the ultra-high degree of 10,800 corresponding to the resolution of one arcmin. Fourier coefficients, computed in double precision, are found to be accurate to 15 significant digits, indicating that the normalized error is close to the machine round-off error. The orthonormality, evaluated with Fourier coefficients and related integrals, is shown to be accurate to O(10?15) for degrees and orders up to 10,800. The Legendre function of degree 10,800 and order 5,000, synthesized from Fourier coefficients, is accurate to the machine round-off error. Further extension of the method to even higher degrees seems to be realizable without significant deterioration of accuracy. The Fourier series is applied to the projection of Legendre functions to the high-resolution global relief data of the National Geophysical Data Center of the National Oceanic and Atmospheric Administration, and the spherical harmonic degree variance (power spectrum) of global relief data is discussed. 相似文献
14.
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 相似文献
15.
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. 相似文献
16.
Least-squares collocation may be used for the estimation of spherical harmonic coefficients and their error and error correlations
from GOCE data. Due to the extremely large number of data, this requires the use of the so-called method of Fast Spherical
Collocation (FSC) which requires that data is gridded equidistantly on each parallel and have the same uncorrelated noise
on the parallel. A consequence of this is that error-covariances will be zero except between coefficients of the same signed
order (i.e., the same order and the same coefficient type C–C or S–S). If the data distribution and the characteristics of the data noise are symmetric with respect to the equator, then, within
a given order and coefficient type, the error-covariances amongst coefficients whose degrees are of different parity also
vanish. The deviation from this “ideal” pattern has been studied using data-sets of second order radial derivatives of the
anomalous potential. A total number of points below 17,000 were used having an equi-angular or an equal area distribution
or being associated with points on a realistic GOCE orbit but close to the nodes of a grid. Also the data were considered
having a correlated or an uncorrelated noise and three different signal covariance functions. Grids including data or not
including data in the polar areas were used. Using the functionals associated with the data, error estimates of coefficients
and error-correlations between coefficients were calculated up to a maximal degree and order equal to 90. As expected, for
the data-distributions with no data in the polar areas the error-estimates were found to be larger than when the polar areas
contained data. In all cases it was found that only the error-correlations between coefficients of the same order were significantly
different from zero (up to 88%). Error-correlations were significantly larger when data had been regarded as having non-zero
error-correlations. Also the error-correlations were largest when the covariance function with the largest signal covariance
distance was used. The main finding of this study was that the correlated noise has more pronounced impact on gridded data
than on data distributed on a realistic GOCE orbit. This is useful information for methods using gridded data, such as FSC. 相似文献
17.
Spherical harmonic modelling to ultra-high degree of Bouguer and isostatic anomalies 总被引：9，自引：4，他引：5
The availability of high-resolution global digital elevation data sets has raised a growing interest in the feasibility of obtaining their spherical harmonic representation at matching resolution, and from there in the modelling of induced gravity perturbations. We have therefore estimated spherical Bouguer and Airy isostatic anomalies whose spherical harmonic models are derived from the Earth’s topography harmonic expansion. These spherical anomalies differ from the classical planar ones and may be used in the context of new applications. We succeeded in meeting a number of challenges to build spherical harmonic models with no theoretical limitation on the resolution. A specific algorithm was developed to enable the computation of associated Legendre functions to any degree and order. It was successfully tested up to degree 32,400. All analyses and syntheses were performed, in 64 bits arithmetic and with semi-empirical control of the significant terms to prevent from calculus underflows and overflows, according to IEEE limitations, also in preserving the speed of a specific regular grid processing scheme. Finally, the continuation from the reference ellipsoid’s surface to the Earth’s surface was performed by high-order Taylor expansion with all grids of required partial derivatives being computed in parallel. The main application was the production of a 1′ × 1′ equiangular global Bouguer anomaly grid which was computed by spherical harmonic analysis of the Earth’s topography–bathymetry ETOPO1 data set up to degree and order 10,800, taking into account the precise boundaries and densities of major lakes and inner seas, with their own altitude, polar caps with bedrock information, and land areas below sea level. The harmonic coefficients for each entity were derived by analyzing the corresponding ETOPO1 part, and free surface data when required, at one arc minute resolution. The following approximations were made: the land, ocean and ice cap gravity spherical harmonic coefficients were computed up to the third degree of the altitude, and the harmonics of the other, smaller parts up to the second degree. Their sum constitutes what we call ETOPG1, the Earth’s TOPography derived Gravity model at 1′ resolution (half-wavelength). The EGM2008 gravity field model and ETOPG1 were then used to rigorously compute 1′ × 1′ point values of surface gravity anomalies and disturbances, respectively, worldwide, at the real Earth’s surface, i.e. at the lower limit of the atmosphere. The disturbance grid is the most interesting product of this study and can be used in various contexts. The surface gravity anomaly grid is an accurate product associated with EGM2008 and ETOPO1, but its gravity information contents are those of EGM2008. Our method was validated by comparison with a direct numerical integration approach applied to a test area in Morocco–South of Spain (Kuhn, private communication 2011) and the agreement was satisfactory. Finally isostatic corrections according to the Airy model, but in spherical geometry, with harmonic coefficients derived from the sets of the ETOPO1 different parts, were computed with a uniform depth of compensation of 30?km. The new world Bouguer and isostatic gravity maps and grids here produced will be made available through the Commission for the Geological Map of the World. Since gravity values are those of the EGM2008 model, geophysical interpretation from these products should not be done for spatial scales below 5 arc minutes (half-wavelength). 相似文献
18.
Fourier transform summation of Legendre series and D-functions 总被引：4，自引：1，他引：3
T. Risbo 《Journal of Geodesy》1996,70(7):383-396
The relation between D- and d-functions, spherical harmonic functions and Legendre functions is reviewed. Dmatrices and irreducible representations of the rotation group O(3) and SU(2) group are briefly reviewed. Two new recursive methods for calculations of D-matrices are presented. Legendre functions are evaluated as part of this scheme. Vector spherical harmonics in the form af generalized spherical harmonics are also included as well as derivatives of the spherical harmonics. The special dmatrices evaluated for argument equal to/2 offer a simple method of calculating the Fourier coefficients of Legendre functions, derivatives of Legendre functions and vector spherical harmonics. Summation of a Legendre series or a full synthesis on the unit sphere of a field can then be performed by transforming the spherical harmonic coefficients to Fourier coefficients and making the summation by an inverse FFT (Fast Fourier Transform). The procedure is general and can also be applied to evaluate derivatives of a field and components of vector and tensor fields. 相似文献