共查询到20条相似文献,搜索用时 15 毫秒
1.
The structure of normal matrices occurring in the problem of weighted least-squares spherical harmonic analysis of measurements
scattered on a sphere with random noises is investigated. Efficient algorithms for the formation of the normal matrices are
derived using fundamental relations inherent to the products of two surface spherical harmonic functions. The whole elements
of a normal matrix complete to spherical harmonic degree L are recursively obtained from its first row or first column extended to degree 2L with only O(L
4) computational operations. Applications of the algorithms to the formation of surface normal matrices from geoid undulations
and surface gravity anomalies are discussed in connection with the high-degree geopotential modeling.
Received: 22 March 1999 / Accepted: 23 December 1999 相似文献
2.
Downward continuation and geoid determination based on band-limited airborne gravity data 总被引:4,自引:3,他引:4
The downward continuation of the harmonic disturbing gravity potential, derived at flight level from discrete observations
of airborne gravity by the spherical Hotine integral, to the geoid is discussed. The initial-boundary-value approach, based
on both the direct and inverse solution to Dirichlet's problem of potential theory, is used. Evaluation of the discretized
Fredholm integral equation of the first kind and its inverse is numerically tested using synthetic airborne gravity data.
Characteristics of the synthetic gravity data correspond to typical airborne data used for geoid determination today and in
the foreseeable future: discrete gravity observations at a mean flight height of 2 to 6 km above mean sea level with minimum
spatial resolution of 2.5 arcmin and a noise level of 1.5 mGal. Numerical results for both approaches are presented and discussed.
The direct approach can successfully be used for the downward continuation of airborne potential without any numerical instabilities
associated with the inverse approach. In addition to these two-step approaches, a one-step procedure is also discussed. This
procedure is based on a direct relationship between gravity disturbances at flight level and the disturbing gravity potential
at sea level. This procedure provided the best results in terms of accuracy, stability and numerical efficiency. As a general
result, numerically stable downward continuation of airborne gravity data can be seen as another advantage of airborne gravimetry
in the field of geoid determination.
Received: 6 June 2001 / Accepted: 3 January 2002 相似文献
3.
The multiresolution character of collocation 总被引:3,自引:0,他引:3
C. Kotsakis 《Journal of Geodesy》2000,74(3-4):275-290
An interesting theoretical connection between the statistical (non-stochastic) collocation principle and the multiresolution/wavelet
framework of signal approximation is presented. The rapid developments in multiresolution analysis theory over the past few
years have provided very useful (theoretical and practical) tools for approximation and spectral studies of irregularly varying
signals, thus opening new possibilities for `non-stationary' gravity field modeling. It is demonstrated that the classic multiresolution
formalism according to Mallat's pioneering work lies at the very core of some of the general approximation principles traditionally
used in physical geodesy problems. In particular, it is shown that the use of a spatio-statistical (non-probabilistic) minimum
mean-square-error criterion for optimal linear estimation of deterministic signals, in conjunction with regularly gridded
data, always gives rise to a generalized multiresolution analysis in the Hilbert space L
2(R), under some mild constraints on the spatial covariance function and the power spectrum of the unknown field under consideration.
Using the theory and the actual approximation algorithms associated with statistical collocation, a new constructive framework
for building generalized multiresolution analyses in L
2(R) is presented, without the need for the usual dyadic restriction that exists in classic wavelet theory. The multiresolution and `non-stationary' aspects of the statistical collocation approximation
procedure are also discussed, and finally some conclusions and recommendations for future work are given.
Received: 26 January 1999 / Accepted: 16 August 1999 相似文献
4.
One of the most basic and important tools in optimal spectral gravity field modelling is the method of Wiener filtering.
Originally developed for applications in analogue signal analysis and communication engineering, Wiener filtering has become
a standard linear estimation technique of modern operational geodesy, either as an independent practical tool for data de-noising
in the frequency domain or as an integral component of a more general signal estimation methodology (input–output systems
theory). Its theoretical framework is based on the Wiener–Kolmogorov linear prediction theory for stationary random fields
in the presence of additive external noise, and thus it is closely related to the (more familiar to geodesists) method of
least-squares collocation with random observation errors. The main drawback of Wiener filtering that makes its use in many
geodetic applications problematic stems from the stationarity assumption for both the signal and the noise involved in the
approximation problem. A modified Wiener-type linear estimation filter is introduced that can be used with noisy data obtained
from an arbitrary deterministic field under the masking of non-stationary random observation errors. In addition, the sampling
resolution of the input data is explicitly taken into account within the estimation algorithm, resulting in a resolution-dependent
optimal noise filter. This provides a more insightful approach to spectral filtering techniques for noise reduction, since
the data resolution parameter has not been directly incorporated in previous formulations of frequency-domain estimation problems
for gravity field signals with discrete noisy data.
Received: 1 November 2000 / Accepted: 19 June 2001 相似文献
5.
H. Nahavandchi 《Journal of Geodesy》2002,76(6-7):345-352
It is suggested that a spherical harmonic representation of the geoidal heights using global Earth gravity models (EGM) might
be accurate enough for many applications, although we know that some short-wavelength signals are missing in a potential coefficient
model. A `direct' method of geoidal height determination from a global Earth gravity model coefficient alone and an `indirect'
approach of geoidal height determination through height anomaly computed from a global gravity model are investigated. In
both methods, suitable correction terms are applied. The results of computations in two test areas show that the direct and
indirect approaches of geoid height determination yield good agreement with the classical gravimetric geoidal heights which
are determined from Stokes' formula. Surprisingly, the results of the indirect method of geoidal height determination yield
better agreement with the global positioning system (GPS)-levelling derived geoid heights, which are used to demonstrate such
improvements, than the results of gravimetric geoid heights at to the same GPS stations. It has been demonstrated that the
application of correction terms in both methods improves the agreement of geoidal heights at GPS-levelling stations. It is
also found that the correction terms in the direct method of geoidal height determination are mostly similar to the correction
terms used for the indirect determination of geoidal heights from height anomalies.
Received: 26 July 2001 / Accepted: 21 February 2002 相似文献
6.
Simplified techniques for high-degree spherical harmonic synthesis are extended to include gravitational potential second
derivatives with respect to latitude.
Received: 23 July 2001 / Accepted: 12 April 2002
Acknowledgement. The authors would like to thank Christian Tscherning for recommending Laplace's equation as an accuracy test. Our use of
Legendre's differential equation, as the most direct means for extending our simplified synthesis methods to second-order
derivatives, was a direct result of this suggestion.
Correspondence to: S. A. Holmes 相似文献
7.
The recovery of a full set of gravity field parameters from satellite gravity gradiometry (SGG) is a huge numerical and computational
task. In practice, parallel computing has to be applied to estimate the more than 90 000 harmonic coefficients parameterizing
the Earth's gravity field up to a maximum spherical harmonic degree of 300. Three independent solution strategies (preconditioned
conjugate gradient method, semi-analytic approach, and distributed non-approximative adjustment), which are based on different
concepts, are assessed and compared both theoretically and on the basis of a realistic-as-possible numerical simulation regarding
the accuracy of the results, as well as the computational effort. Special concern is given to the correct treatment of the
coloured noise characteristics of the gradiometer. The numerical simulations show that the three methods deliver nearly identical
results—even in the case of large data gaps in the observation time series. The newly proposed distributed non-approximative
adjustment approach, which is the only one of the three methods that solves the inverse problem in a strict sense, also turns
out to be a feasible method for practical applications.
Received: 17 December 2001 / Accepted: 17 July 2002
Acknowledgments. We would like to thank Prof. W.-D. Schuh, Institute of Theoretical Geodesy, University of Bonn, for providing us with the
serial version of the PCGMA algorithm, which forms the basis for the parallel PCGMA package developed at our institute. This
study was partially performed in the course of the GOCE project `From E?tv?s to mGal+', funded by the European Space Agency
(ESA) under contract No. 14287/00/NL/DC.
Correspondence to: R. Pail 相似文献
8.
A synthetic Earth for use in geodesy 总被引:1,自引:0,他引:1
R. Haagmans 《Journal of Geodesy》2000,74(7-8):503-511
A synthetic Earth and its gravity field that can be represented at different resolutions for testing and comparing existing
and new methods used for global gravity-field determination are created. Both the boundary and boundary values of the gravity
potential can be generated. The approach chosen also allows observables to be generated at aircraft flight height or at satellite
altitude. The generation of the synthetic Earth shape (SES) and gravity-field quantities is based upon spherical harmonic
expansions of the isostatically compensated equivalent rock topography and the EGM96 global geopotential model. Spherical
harmonic models are developed for both the synthetic Earth topography (SET) and the synthetic Earth potential (SEP) up to
degree and order 2160 corresponding to a 5′×5′ resolution. Various sets of SET, SES and SEP with boundary geometry and boundary
values at different resolutions can be generated using low-pass filters applied to the expansions. The representation is achieved
in point sets based upon refined triangulation of a octahedral geometry projected onto the chosen reference ellipsoid. The
filter cut-offs relate to the sampling pattern in order to avoid aliasing effects. Examples of the SET and its gravity field
are shown for a resolution with a Nyquist sampling rate of 8.27 degrees.
Received: 6 August 1999 / Accepted: 26 April 2000 相似文献
9.
The Cartesian moments of the mass density of a gravitating body and the spherical harmonic coefficients of its gravitational
field are related in a peculiar way. In particular, the products of inertia can be expressed by the spherical harmonic coefficients
of the gravitational potential as was derived by MacCullagh for a rigid body. Here the MacCullagh formulae are extended to
a deformable body which is restricted to radial symmetry in order to apply the Love–Shida hypothesis. The mass conservation
law allows a representation of the incremental mass density by the respective excitation function. A representation of an
arbitrary Cartesian monome is always possible by sums of solid spherical harmonics multiplied by powers of the radius. Introducing
these representations into the definition of the Cartesian moments, an extension of the MacCullagh formulae is obtained. In
particular, for excitation functions with a vanishing harmonic coefficient of degree zero, the (diagonal) incremental moments
of inertia also can be represented by the excitation coefficients. Four types of excitation functions are considered, namely:
(1) tidal excitation; (2) loading potential; (3) centrifugal potential; and (4) transverse surface stress. One application
of the results could be model computation of the length-of-day variations and polar motion, which depend on the moments of
inertia.
Received: 27 July 1999 / Accepted: 24 May 2000 相似文献
10.
Fast spherical collocation: theory and examples 总被引:2,自引:4,他引:2
It has long been known that a spherical harmonic analysis of gridded (and noisy) data on a sphere (with uniform error for
a fixed latitude) gives rise to simple systems of equations. This idea has been generalized for the method of least-squares
collocation, when using an isotropic covariance function or reproducing kernel. The data only need to be at the same altitude
and of the same kind for each latitude. This permits, for example, the combination of gravity data at the surface of the Earth
and data at satellite altitude, when the orbit is circular. Suppose that data are associated with the points of a grid with
N values in latitude and M values in longitude. The latitudes do not need to be spaced uniformly. Also suppose that it is required to determine the
spherical harmonic coefficients to a maximal degree and order K. Then the method will require that we solve K systems of equations each having a symmetric positive definite matrix of only N × N. Results of simulation studies using the method are described.
Received: 18 October 2001 / Accepted: 4 October 2002
Correspondence to: F. Sansò 相似文献
11.
H. Nahavandchi 《Journal of Geodesy》2000,74(6):488-496
The direct topographical correction is composed of both local effects and long-wavelength contributions. This implies that
the classical integral formula for determining the direct effect may have some numerical problems in representing these different
signals. On the other hand, a representation by a set of harmonic coefficients of the topography to, say, degree and order
360 will omit significant short-wavelength signals. A new formula is derived by combining the classical formula and a set
of spherical harmonics. Finally, the results of this solution are compared with the Moritz topographical correction in a test
area.
Received: 27 July 1998 / Accepted: 29 March 2000 相似文献
12.
A unified approach to the Clenshaw summation and the recursive computation of very high degree and order normalised associated Legendre functions 总被引:11,自引:7,他引:11
Spherical harmonic expansions form partial sums of fully normalised associated Legendre functions (ALFs). However, when evaluated
increasingly close to the poles, the ultra-high degree and order (e.g. 2700) ALFs range over thousands of orders of magnitude.
This causes existing recursion techniques for computing values of individual ALFs and their derivatives to fail. A common
solution in geodesy is to evaluate these expansions using Clenshaw's method, which does not compute individual ALFs or their
derivatives. Straightforward numerical principles govern the stability of this technique. Elementary algebra is employed to
illustrate how these principles are implemented in Clenshaw's method. It is also demonstrated how existing recursion algorithms
for computing ALFs and their first derivatives are easily modified to incorporate these same numerical principles. These modified
recursions yield scaled ALFs and first derivatives, which can then be combined using Horner's scheme to compute partial sums,
complete to degree and order 2700, for all latitudes (except at the poles for first derivatives). This exceeds any previously
published result. Numerical tests suggest that this new approach is at least as precise and efficient as Clenshaw's method.
However, the principal strength of the new techniques lies in their simplicity of formulation and implementation, since this
quality should simplify the task of extending the approach to other uses, such as spherical harmonic analysis.
Received: 30 June 2000 / Accepted: 12 June 2001 相似文献
13.
In October 1998 the IGEX field campaign, the first coordinated international effort to monitor GLONASS satellites on global
basis, was started. Currently about 40 institutions worldwide support this effort either by providing GLONASS tracking data
or in operating related data and analysis centers. The increasing quality and consistency of the calculated GLONASS orbits
(about 25 cm early in 2000), even after the end of the official IGEX field campaign, are shown. Particular attention is drawn
to the combination of precise ephemerides in order to generate a robust, reliable and complete IGEX orbits product. Some problems
in modeling the effect of solar radiation pressure on GLONASS satellites are demonstrated. Finally, the expected benefits
and prospects of the upcoming International GLOnass Service-Pilot Project (IGLOS-PP) of the International GPS Service (IGS)
are discussed in more detail.
Received: 17 August 2000 / Accepted: 12 April 2001 相似文献
14.
Local geoid determination combining gravity disturbances and GPS/levelling: a case study in the Lake Nasser area, Aswan, Egypt 总被引:1,自引:0,他引:1
C. C. Tscherning Awar Radwan A. A. Tealeb S. M. Mahmoud M. Abd El-Monum Ramdan Hassan I. El-Syaed K. Saker 《Journal of Geodesy》2001,75(7-8):343-348
The use of GPS for height control in an area with existing levelling data requires the determination of a local geoid and
the bias between the local levelling datum and the one implicitly defined when computing the local geoid. If only scarse gravity
data are available, the heights of new data may be collected rapidly by determining the ellipsoidal height by GPS and not
using orthometric heights. Hence the geoid determination has to be based on gravity disturbances contingently combined with
gravity anomalies. Furthermore, existing GPS/levelling data may also be used in the geoid determination if a suitable general
gravity field modelling method (such as least-squares collocation, LSC) is applied. A comparison has been made in the Aswan
Dam area between geoids determined using fast Fourier transform (FFT) with gravity disturbances exclusively and LSC using
only the gravity disturbances and the disturbances combined with GPS/levelling data. The EGM96 spherical harmonic model was
in all cases used in a remove–restore mode. A total of 198 gravity disturbances spaced approximately 3 km apart were used,
as well as 35 GPS/levelling points in the vicinity and on the Aswan Dam. No data on the Nasser Lake were available. This gave
difficulties when using FFT, which requires the use of gridded data. When using exclusively the gravity disturbances, the
agreement between the GPS/levelling data were 0.71 ± 0.17 m for FFT and 0.63 ± 0.15 for LSC. When combining gravity disturbances
and GPS/levelling, the LSC error estimate was ±0.10 m. In the latter case two bias parameters had to be introduced to account
for a possible levelling datum difference between the levelling on the dam and that on the adjacent roads.
Received: 14 August 2000 / Accepted: 28 February 2001 相似文献
15.
Some modifications of Stokes' formula that account for truncation and potential coefficient errors 总被引:1,自引:1,他引:0
Stokes' formula from 1849 is still the basis for the gravimetric determination of the geoid. The modification of the formula,
originating with Molodensky, aims at reducing the truncation error outside a spherical cap of integration. This goal is still
prevalent among various modifications. In contrast to these approaches, some least-squares types of modification that aim
at reducing the truncation error, as well as the error stemming from the potential coefficients, are demonstrated. The least-squares
estimators are provided in the two cases that (1) Stokes' kernel is a priori modified (e.g. according to Molodensky's approach)
and (2) Stokes' kernel is optimally modified to minimize the global mean square error. Meissl-type modifications are also
studied. In addition, the use of a higher than second-degree reference field versus the original (Pizzetti-type) reference
field is discussed, and it is concluded that the former choice of reference field implies increased computer labour to achieve
the same result as with the original reference field.
Received: 14 December 1998 / Accepted: 4 October 1999 相似文献
16.
Simplified formulae for the BIQUE estimation of variance components in disjunctive observation groups 总被引:7,自引:1,他引:6
General rigorous and simplified formulae are reported for the best invariant quadratic unbiased estimates of the variance–covariance
components, which can be applied to all least-squares adjustments with the general linear stochastic model. Simplified procedures
are given for two cases frequently recurring in geodetic applications: uncorrelated groups of correlated or uncorrelated observations,
with more than one variance component in each group.
Received: 19 November 1998 / Accepted: 21 March 2000 相似文献
17.
L. E. Sjöberg 《Journal of Geodesy》2001,75(5-6):283-290
The topographic and atmospheric effects of gravimetric geoid determination by the modified Stokes formula, which combines
terrestrial gravity and a global geopotential model, are presented. Special emphasis is given to the zero- and first-degree
effects. The normal potential is defined in the traditional way, such that the disturbing potential in the exterior of the
masses contains no zero- and first-degree harmonics. In contrast, it is shown that, as a result of the topographic masses,
the gravimetric geoid includes such harmonics of the order of several centimetres. In addition, the atmosphere contributes
with a zero-degree harmonic of magnitude within 1 cm.
Received: 5 November 1999 / Accepted: 22 January 2001 相似文献
18.
An efficient algorithm is proposed for gravity field recovery from Gravity Field and Steady-State Ocean Circulation Explorer
(GOCE) satellite gravity gradient observations. The mathematical model is formulated in the time domain, which allows the
inclusion of realistic observational noise models. The algorithm combines the iterative solution of the normal equations,
using a Richardson-type iteration scheme, with the fast computation of the right-hand side of the normal equations in each
iteration step by a suitable approximation of the design matrix. The convergence of the iteration is investigated, error estimates
are provided, and the unbiasedness of the method is proved. It is also shown that the method does not converge to the solution
of the normal equations. The performance of the approach for white noise and coloured noise is demonstrated along a simulated
GOCE orbit up to spherical harmonic degree and order 180. The results also indicate that the approximation error may be neglected.
Received: 30 November 1999 / Accepted: 31 May 2000 相似文献
19.
A method is presented for performing global spherical harmonic computation by two-dimensional Fourier transformations. The method goes back to old literature (Schuster 1902) and tackles the problem of non-orthogonality of Legendre-functions, when discretized on an equi-angular grid. Both analysis and synthesis relations are presented, which link the spherical harmonic spectrum to a two-dimensional Fourier spectrum. As an alternative, certain functions of co-latitude are introduced, which are orthogonal to discretized Legendre functions. Several independent Fourier approaches for spherical harmonic computation fit into our general scheme. 相似文献
20.
E. W. Grafarend 《Journal of Geodesy》2001,75(7-8):363-390
In a comparison of the solution of the spherical horizontal and vertical boundary value problems of physical geodesy it is
aimed to construct downward continuation operators for vertical deflections (surface gradient of the incremental gravitational
potential) and for gravity disturbances (vertical derivative of the incremental gravitational potential) from points on the
Earth's topographic surface or of the three-dimensional (3-D) Euclidean space nearby down to the international reference sphere
(IRS). First the horizontal and vertical components of the gravity vector, namely spherical vertical deflections and spherical
gravity disturbances, are set up. Second, the horizontal and vertical boundary value problem in spherical gravity and geometry
space is considered. The incremental gravity vector is represented in terms of vector spherical harmonics. The solution of
horizontal spherical boundary problem in terms of the horizontal vector-valued Green function converts vertical deflections
given on the IRS to the incremental gravitational potential external in the 3-D Euclidean space. The horizontal Green functions
specialized to evaluation and source points on the IRS coincide with the Stokes kernel for vertical deflections. Third, the
vertical spherical boundary value problem is solved in terms of the vertical scalar-valued Green function. Fourth, the operators
for upward continuation of vertical deflections given on the IRS to vertical deflections in its external 3-D Euclidean space
are constructed. Fifth, the operators for upward continuation of incremental gravity given on the IRS to incremental gravity
to the external 3-D Euclidean space are generated. Finally, Meissl-type diagrams for upward continuation and regularized downward
continuation of horizontal and vertical gravity data, namely vertical deflection and incremental gravity, are produced.
Received: 10 May 2000 / Accepted: 26 February 2001 相似文献