共查询到20条相似文献，搜索用时 31 毫秒
1.
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. 相似文献
2.
Summary The problem to detect configurational defects in geodetic networks is solved by a graph-theoretical algorithm, here applied
to triangular geodetic networks and being presented as a computer program in the Appendix. Based on an analysis of the incidence
matrix the algorithm detects, for instance, missing vertical directions which cause two types of deficiencies. In case of
only vertical direction measurements from one point to another, but no counter vertical direction measurements backwards,
the rank deficiency of the first type is identified. Furtheron if there are “bare” points with no vertical direction measurements
at all, the rank deficiency of the second type is found. The algorithm has proved a rank deficiency of 4+13=17 in theSW Finland triangular network which before has been found as the surprizing rank defect ofTAGNET 3d-operational adjustment. 相似文献
3.
The use of sampling-based Monte Carlo methods for the computation and propagation of large covariance matrices in geodetic applications is investigated. In particular, the so-called Gibbs sampler, and its use in deriving covariance matrices by Monte Carlo integration, and in linear and nonlinear error propagation studies, is discussed. Modifications of this technique are given which improve in efficiency in situations where estimated parameters are highly correlated and normal matrices appear as ill-conditioned. This is a situation frequently encountered in satellite gravity field modelling. A synthetic experiment, where covariance matrices for spherical harmonic coefficients are estimated and propagated to geoid height covariance matrices, is described. In this case, the generated samples correspond to random realizations of errors of a gravity field model.
AcknowledgementsThe authors are indebted to Pieter Visser and Pavel Ditmar for providing simulation output that was used in the GOCE error generation experiments. Furthermore, the NASA/NIMA/OSU team is acknowledged for providing public ftp access to the EGM96 error covariance matrix. The two anonymous reviewers are thanked for their valuable comments. 相似文献
4.
Summary
Riemann polar/normal coordinates are the constituents to generate the oblique azimuthal projection of geodesic type, here applied to the reference ellipsoid of revolution (biaxial ellipsoid).Firstly we constitute a minimal atlas of the biaxial ellipsoid built on {ellipsoidal longitude, ellipsoidal latitude} and {metalongitude, metalatitude}. TheDarboux equations of a 1-dimensional submanifold (curve) in a 2-dimensional manifold (biaxial ellipsoid) are reviewed, in particular to represent geodetic curvature, geodetic torsion and normal curvature in terms of elements of the first and second fundamental form as well as theChristoffel symbols. The notion of ageodesic anda geodesic circle is given and illustrated by two examples. The system of twosecond order ordinary differential equations of ageodesic (Lagrange portrait) is presented in contrast to the system of twothird order ordinary differential equations of ageodesic circle (Proofs are collected inAppendix A andB). A precise definition of theRiemann mapping/mapping of geodesics into the local tangent space/tangent plane has been found.Secondly we computeRiemann polar/normal coordinates for the biaxial ellipsoid, both in theLagrange portrait (Legendre series) and in theHamilton portrait (Lie series).Thirdly we have succeeded in a detailed deformation analysis/Tissot distortion analysis of theRiemann mapping. The eigenvalues — the eigenvectors of the Cauchy-Green deformation tensor by means of ageneral eigenvalue-eigenvector problem have been computed inTable 3.1 andTable 3.2 (1, 2 = 1) illustrated inFigures 3.1, 3.2 and3.3. Table 3.3 contains the representation ofmaximum angular distortion of theRiemann mapping. Fourthly an elaborate global distortion analysis with respect toconformal Gau-Krüger, parallel Soldner andgeodesic Riemann coordinates based upon theAiry total deformation (energy) measure is presented in a corollary and numerically tested inTable 4.1. In a local strip [-l
E,l
E] = [-2°, +2°], [b
S,b
N] = [-2°, +2°]Riemann normal coordinates generate the smallest distortion, next are theparallel Soldner coordinates; the largest distortion by far is met by theconformal Gau-Krüger coordinates. Thus it can be concluded that for mapping of local areas of the biaxial ellipsoid surface the oblique azimuthal projection of geodesic type/Riemann polar/normal coordinates has to be favored with respect to others. 相似文献
5.
The resolution of a nonlinear parametric adjustment model is addressed through an isomorphic geometrical setup with tensor
structure and notation, represented by a u-dimensional “model surface” embedded in a flat n-dimensional “observational space”.
Then observations correspond to the observational-space coordinates of the pointQ, theu initial parameters correspond to the model-surface coordinates of the “initial” pointP, and theu adjusted parameters correspond to the model-surface coordinates of the “least-squares” point
. The least-squares criterion results in a minimum-distance property implying that the vector
Q must be orthogonal to the model surface. The geometrical setup leads to the solution of modified normal equations, characterized
by a positive-definite matrix. The latter contains second-order and, optionally, thirdorder partial derivatives of the observables
with respect to the parameters. This approach significantly shortens the convergence process as compared to the standard (linearized)
method. 相似文献
6.
J. Feltens 《Journal of Geodesy》2008,82(8):493-504
Vector-based algorithms for the computation of azimuth, elevation and the ellipsoidal normal unit vector from 3D Cartesian
coordinates are presented. As a by-product, the formulae for the ellipsoidal normal vector can also be used to iteratively
transform rectangular Cartesian coordinates (X, Y, Z) into geodetic coordinates (φ, λ, h) for a height range from −5600 km to 108 km. Comparisons with existing methods indicate that the new transformation can compete with them. 相似文献
7.
C-D. Zhang H.T. Hsu X.P. Wu S.S. Li Q.B. Wang H.Z. Chai L. Du 《Journal of Geodesy》2005,79(8):413-420
The algorithm to transform from 3D Cartesian to geodetic coordinates is obtained by solving the equation of the Lagrange parameter.
Numerical experiments show that geodetic height can be recovered to 0.5 mm precision over the range from −6×106 to 1010 m.
Electronic Supplementary Material: Supplementary material is available in the online version of this article at 相似文献
8.
The weighted Procrustes algorithm is presented as a very effective tool for solving the three-dimensional datum transformation
problem. In particular, the weighted Procrustes algorithm does not require any initial datum parameters for linearization
or any iteration procedure. As a closed-form algorithm it only requires the values of Cartesian coordinates in both systems
of reference. Where there is some prior information about the variance–covariance matrix of the two sets of Cartesian coordinates,
also called pseudo-observations, the weighted Procrustes algorithm is able to incorporate such a quality property of the input
data by means of a proper choice of weight matrix. Such a choice is based on a properly designed criterion matrix which is
discussed in detail. Thanks to the weighted Procrustes algorithm, the problem of incorporating the stochasticity measures
of both systems of coordinates involved in the seven parameter datum transformation problem [conformal group ℂ7(3)] which is free of linearization and any iterative procedure can be considered to be solved. Illustrative examples are
given.
Received: 7 January 2002 / Accepted: 9 September 2002
Correspondence to: E. W. Grafarend 相似文献
9.
《测量评论》2013,45(66):166-174
AbstractThe computation of geographical coordinates in a geodetic triangulation is usually carried out using Puissant's method, in which the assumption is made the sphere radius ν (the radius of curvature of the spheroid perpendicular to the meridian) not only touches the spheroid along the whole small circle of latitude ?,but also, since ρ (the radius of curvature in meridian) is very nearly equal to ν it makes such close contact with the spheroid that the lengths of sides and angles of a geodetic triangle may be considered identical on both sphere and spheroid. 相似文献
10.
Wavelet evaluation of the Stokes and Vening Meinesz integrals 总被引：1，自引：0，他引：1
The wavelet transform is a powerful tool in evaluating some singular geodetic integrals. Due to its localization properties in both the time (space) and frequency (scale) domains, and because the kernels of some geodetic integrals have singular points and decay smoothly and quickly away from the singularities, many wavelet transform coefficients of the kernels become zeros or negligible, and only a small number of wavelet transform coefficients are significant. It is thus possible to significantly compress the kernels of these integrals on a wavelet basis by neglecting the zero coefficients and the small coefficients below a certain threshold. Therefore, wavelets provide a convenient way of efficiently evaluating these integrals in terms of fast computation and savings of computer memory. A modified algorithm for the wavelet evaluation of Stokes' integral is presented. The same modified algorithm is applied to the evaluation of the Vening Meinesz integral, whose kernel has a stronger singularity than does Stokes' kernel. Numerical examples illustrate the efficiency and accuracy of the wavelet methods.Acknowledgments.The author express their sincere thanks to Dr. Salamonwicz for providing his PhD thesis. E-mail correspondence between the authors and Dr. Barthelmes and Dr. Benciolini contributed to the work. R. Benciolini and the other two anonymous reviewers are thanked for their constructive comments. Support for this research was provided by research grants to Dr. Sideris from the Natural Sciences and Engineering Reserch Council of Canada (NSERC) and the Geomatics for Informed Decisions (GEOIDE) Network of Centres of Excellence. The MATLAB Wavelet Toolbox package was used as the platform to develop the software in this project. 相似文献
11.
Béla Paláncz Joseph L. Awange Piroska Zaletnyik Robert H. Lewis 《Journal of Geodesy》2010,84(1):79-95
A fundamental task in geodesy is solving systems of equations. Many geodetic problems are represented as systems of multivariate
polynomials. A common problem in solving such systems is improper initial starting values for iterative methods, leading to
convergence to solutions with no physical meaning, or to convergence that requires global methods. Though symbolic methods
such as Groebner bases or resultants have been shown to be very efficient, i.e., providing solutions for determined systems
such as 3-point problem of 3D affine transformation, the symbolic algebra can be very time consuming, even with special Computer
Algebra Systems (CAS). This study proposes the Linear Homotopy method that can be implemented easily in high-level computer languages like C++ and Fortran that are faster than CAS by at
least two orders of magnitude. Using Mathematica, the power of Homotopy is demonstrated in solving three nonlinear geodetic problems: resection, GPS positioning, and affine
transformation. The method enlarging the domain of convergence is found to be efficient, less sensitive to rounding of numbers,
and has lower complexity compared to other local methods like Newton–Raphson. 相似文献
12.
Marcin Ligas 《Journal of Geodesy》2012,86(4):249-256
A new method of transforming Cartesian to geodetic (or planetographic) coordinates on a triaxial ellipsoid is presented. The
method is based on simple reasoning coming from essentials of vector calculus. The reasoning results in solving a nonlinear
system of equations for coordinates of the point being the projection of a point located outside or inside a triaxial ellipsoid
along the normal to the ellipsoid. The presented method has been compared to a vector method of Feltens (J Geod 83:129–137,
2009) who claims that no other methods are available in the literature. Generally, our method turns out to be more accurate, faster
and applicable to celestial bodies characterized by different geometric parameters. The presented method also fits to the
classical problem of converting Cartesian to geodetic coordinates on the ellipsoid of revolution. 相似文献
13.
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. 相似文献
14.
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. 相似文献
15.
16.
The gravitational attraction of any polygonally shaped vertical prism with inclined top and bottom faces 总被引：2，自引：0，他引：2
D. A. Smith 《Journal of Geodesy》2000,74(5):414-420
A closed formula for computing the gravitational attraction of a general vertical prism with N+2 faces (N faces are vertical planes, the other two are the inclined top and bottom planes) in Cartesian coordinates is presented. In
addition, the special case of a triangular prism is discussed. Algebraic differences and overlooked singularity conditions
of a previously published formula of this computation (which was only for the triangular special case) were identified and
are also presented.
Received: 22 March 1999 / Accepted: 2 February 2000 相似文献
17.
Younghee Kwak Mathis Bloßfeld Ralf Schmid Detlef Angermann Michael Gerstl Manuela Seitz 《Journal of Geodesy》2018,92(9):1047-1061
The Celestial Reference System (CRS) is currently realized only by Very Long Baseline Interferometry (VLBI) because it is the space geodetic technique that enables observations in that frame. In contrast, the Terrestrial Reference System (TRS) is realized by means of the combination of four space geodetic techniques: Global Navigation Satellite System (GNSS), VLBI, Satellite Laser Ranging (SLR), and Doppler Orbitography and Radiopositioning Integrated by Satellite. The Earth orientation parameters (EOP) are the link between the two types of systems, CRS and TRS. The EOP series of the International Earth Rotation and Reference Systems Service were combined of specifically selected series from various analysis centers. Other EOP series were generated by a simultaneous estimation together with the TRF while the CRF was fixed. Those computation approaches entail inherent inconsistencies between TRF, EOP, and CRF, also because the input data sets are different. A combined normal equation (NEQ) system, which consists of all the parameters, i.e., TRF, EOP, and CRF, would overcome such an inconsistency. In this paper, we simultaneously estimate TRF, EOP, and CRF from an inter-technique combined NEQ using the latest GNSS, VLBI, and SLR data (2005–2015). The results show that the selection of local ties is most critical to the TRF. The combination of pole coordinates is beneficial for the CRF, whereas the combination of \(\varDelta \hbox {UT1}\) results in clear rotations of the estimated CRF. However, the standard deviations of the EOP and the CRF improve by the inter-technique combination which indicates the benefits of a common estimation of all parameters. It became evident that the common determination of TRF, EOP, and CRF systematically influences future ICRF computations at the level of several \(\upmu \)as. Moreover, the CRF is influenced by up to \(50~\upmu \)as if the station coordinates and EOP are dominated by the satellite techniques. 相似文献
18.
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. 相似文献
19.
Computing geodetic coordinates from geocentric coordinates 总被引：1，自引：1，他引：0
H. Vermeille 《Journal of Geodesy》2004,78(1-2):94-95
A closed-form algebraic method to transform geocentric coordinates to geodetic coordinates has previously been proposed. The validity domain of latitude and height formulae in the vicinity of the Earths core is specified. A new expression of longitude is proposed, excluding indetermination and sensitivity to round-off error around the ±180 degrees longitude discontinuity. 相似文献
20.
Contribution of Starlette,Stella, and AJISAI to the SLR-derived global reference frame 总被引：3，自引：3，他引：0
Krzysztof Sośnica Adrian Jäggi Daniela Thaller Gerhard Beutler Rolf Dach 《Journal of Geodesy》2014,88(8):789-804
The contribution of Starlette, Stella, and AJISAI is currently neglected when defining the International Terrestrial Reference Frame, despite a long time series of precise SLR observations and a huge amount of available data. The inferior accuracy of the orbits of low orbiting geodetic satellites is the main reason for this neglect. The Analysis Centers of the International Laser Ranging Service (ILRS ACs) do, however, consider including low orbiting geodetic satellites for deriving the standard ILRS products based on LAGEOS and Etalon satellites, instead of the sparsely observed, and thus, virtually negligible Etalons. We process ten years of SLR observations to Starlette, Stella, AJISAI, and LAGEOS and we assess the impact of these Low Earth Orbiting (LEO) SLR satellites on the SLR-derived parameters. We study different orbit parameterizations, in particular different arc lengths and the impact of pseudo-stochastic pulses and dynamical orbit parameters on the quality of the solutions. We found that the repeatability of the East and North components of station coordinates, the quality of polar coordinates, and the scale estimates of the reference are improved when combining LAGEOS with low orbiting SLR satellites. In the multi-SLR solutions, the scale and the \(Z\) component of geocenter coordinates are less affected by deficiencies in solar radiation pressure modeling than in the LAGEOS-1/2 solutions, due to substantially reduced correlations between the \(Z\) geocenter coordinate and empirical orbit parameters. Eventually, we found that the standard values of Center-of-mass corrections (CoM) for geodetic LEO satellites are not valid for the currently operating SLR systems. The variations of station-dependent differential range biases reach 52 and 25 mm for AJISAI and Starlette/Stella, respectively, which is why estimating station-dependent range biases or using station-dependent CoM, instead of one value for all SLR stations, is strongly recommended. This clearly indicates that the ILRS effort to produce CoM corrections for each satellite, which are site-specific and depend on the system characteristics at the time of tracking, is very important and needs to be implemented in the SLR data analysis. 相似文献