共查询到20条相似文献,搜索用时 125 毫秒
1.
Iterative vector methods for computing geodetic latitude and height from rectangular coordinates 总被引:4,自引:4,他引:4
J. Pollard 《Journal of Geodesy》2002,76(1):36-40
Two iterative vector methods for computing geodetic coordinates (φ, h) from rectangular coordinates (x, y, z) are presented. The methods are conceptually simple, work without modification at any latitude and are easy to program. Geodetic
latitude and height can be calculated to acceptable precision in one iteration over the height range from −106 to +109 m.
Received: 13 December 2000 / Accepted: 13 July 2001 相似文献
2.
G. Even-Tzur 《Journal of Geodesy》2002,76(8):455-461
The performance of geodetic monitoring networks is heavily influenced by the configuration of the measured GPS vectors. As
an effective design of the GPS measurements will decrease GPS campaign costs and increase the accuracy and reliability of
the entire network, the identification of the preferred GPS vectors for measurement has been highlighted as a core problem
in the process of deformation monitoring. An algorithm based on a sensitivity analysis of the network, as dependent upon a
postulated velocity field, is suggested for the selection of the optimal GPS vectors. Relevant mathematical and statistical
concepts are presented as the basis for an improved method of vector configuration design. A sensitivity analysis of the geodetic
geodynamic network in the north of Israel is presented, where the method is examined against two deformation models, the Simple
Transform Fault and the Locked Fault. The proposed method is suggested as a means for the improvement of the design of monitoring
networks, a common practice worldwide.
Received: 30 July 2001 / Accepted: 3 June 2002
Acknowledgments. It is my pleasant duty to thank the Survey of Israel and Dr. E. Ostrovsky for providing the variance–covariance matrix of
the G1 network in northern Israel. I would like to thank the reviewers of this paper for their constructive and helpful remarks. 相似文献
3.
R. Lehmann 《Journal of Geodesy》2000,74(3-4):327-334
The definition and connection of vertical datums in geodetic height networks is a fundamental problem in geodesy. Today,
the standard approach to solve it is based on the joint processing of terrestrial and satellite geodetic data. It is generalized
to cases where the coverage with terrestrial data may change from region to region, typically across coastlines. The principal
difficulty is that such problems, so-called altimetry–gravimetry boundary-value problems (AGPs), do not admit analytical solutions
such as Stokes' integral. A numerical solution strategy for the free-datum problem is presented. Analysis of AGPs in spherical
and constant radius approximation shows that two of them are mathematically well-posed problems, while the classical AGP-I
may be ill posed in special situations.
Received: 2 December 1998 / Accepted: 30 November 1999 相似文献
4.
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 相似文献
5.
H. Vermeille 《Journal of Geodesy》2002,76(8):451-454
The transformation from geocentric coordinates to geodetic coordinates is usually carried out by iteration. A closed-form
algebraic method is proposed, valid at any point on the globe and in space, including the poles, regardless of the value of
the ellipsoid's eccentricity.
Received: 14 August 2000 / Accepted: 26 June 2002 相似文献
6.
The Somigliana–Pizzetti gravity field (the International gravity formula), namely the gravity field of the level ellipsoid
(the International Reference Ellipsoid), is derived to the sub-nanoGal accuracy level in order to fulfil the demands of modern
gravimetry (absolute gravimeters, super conducting gravimeters, atomic gravimeters). Equations (53), (54) and (59) summarise
Somigliana–Pizzetti gravity Γ(φ,u) as a function of Jacobi spheroidal latitude φ and height u to the order ?(10−10 Gal), and Γ(B,H) as a function of Gauss (surface normal) ellipsoidal latitude B and height H to the order ?(10−10 Gal) as determined by GPS (`global problem solver'). Within the test area of the state of Baden-Württemberg, Somigliana–Pizzetti
gravity disturbances of an average of 25.452 mGal were produced. Computer programs for an operational application of the new
international gravity formula with (L,B,H) or (λ,φ,u) coordinate inputs to a sub-nanoGal level of accuracy are available on the Internet.
Received: 23 June 2000 / Accepted: 2 January 2001 相似文献
7.
The latest gravimetric geoid model for Japan, JGEOID2000, was successfully combined with the nationwide net of GPS at benchmarks,
yielding a new hybrid geoid model for Japan, GSIGEO2000. The least-squares collocation (LSC) method was applied as an interpolation
for fitting JGEOID2000 to the GPS/leveling geoid undulations. The GPS/leveling geoid undulation data were reanalyzed in advance,
in terms of three-dimensional positions from GPS and orthometric heights from leveling. The new hybrid geoid model is, therefore,
compatible with the new Japanese geodetic reference frame. GSIGEO2000 was evaluated internally and independently and the precision
was estimated at 4 cm throughout nearly the whole region.
Received: 15 October 2001 / Accepted: 27 March 2002
Acknowledgments. Messrs. Toshio Kunimi and Tadashi Saito at the Third Geodetic Division of the Geographical Survey Institute (GSI) mainly
carried out the computations of most of the updated leveled heights. With regard to the reanalysis of GPS data, the discussions
with Messrs. Yuki Hatanaka and Shoichi Matsumura of GSI were of great help in building the analysis strategy. Messrs. Kazuyuki
Tanaka and Hiromi Shigematsu collaborated in the preparatory stages of GPS data computation. The authors' thanks are extended
to these colleagues. Some plots were made by GMT software (Wessel and Smith 1991).
Correspondence to: Y. Kuroishi 相似文献
8.
A solution to the downward continuation effect on the geoid determined by Stokes' formula 总被引:2,自引:1,他引:2
L.E. Sjöberg 《Journal of Geodesy》2003,77(1-2):94-100
The analytical continuation of the surface gravity anomaly to sea level is a necessary correction in the application of Stokes'
formula for geoid estimation. This process is frequently performed by the inversion of Poisson's integral formula for a sphere.
Unfortunately, this integral equation corresponds to an improperly posed problem, and the solution is both numerically unstable,
unless it is well smoothed, and tedious to compute. A solution that avoids the intermediate step of downward continuation
of the gravity anomaly is presented. Instead the effect on the geoid as provided by Stokes' formula is studied directly. The
practical solution is partly presented in terms of a truncated Taylor series and partly as a truncated series of spherical
harmonics. Some simple numerical estimates show that the solution mostly meets the requests of a 1-cm geoid model, but the
truncation error of the far zone must be studied more precisely for high altitudes of the computation point. In addition,
it should be emphasized that the derived solution is more computer efficient than the detour by Poisson's integral.
Received: 6 February 2002 / Accepted: 18 November 2002
Acknowledgements. Jonas ?gren carried out the numerical calculations and gave some critical and constructive remarks on a draft version of
the paper. This support is cordially acknowledged. Also, the thorough work performed by one unknown reviewer is very much
appreciated. 相似文献
9.
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 相似文献
10.
Dual algebraic formulation of differential GPS 总被引:1,自引:0,他引:1
A new approach to differential GPS is presented. The corresponding theoretical framework calls on elementary concepts of
algebraic graph theory. The notion of double difference, which is related to that of closure in the sense of Kirchhoff, is
revisited in this context. The Moore–Penrose pseudo-inverse of the closure operator plays a key role in the corresponding
dual formulation. This approach, which is very attractive from a conceptual point of view, sheds a new light on the Teunissen
formulation.
Received: 21 May 2002 / Accepted: 18 November 2002
Correspondence to: A. Lannes 相似文献
11.
The standard analytical approach which is applied for constructing geopotential models OSU86 and earlier ones, is based on
reducing the boundary value equation to a sphere enveloping the Earth and then solving it directly with respect to the potential
coefficients
n,m
. In an alternative procedure, developed by Jekeli and used for constructing the models OSU91 and EGM96, at first an ellipsoidal
harmonic series is developed for the geopotential and then its coefficients
n,m
e
are transformed to the unknown
n,m
. The second solution is more exact, but much more complicated. The standard procedure is modified and a new simple integral
formula is derived for evaluating the potential coefficients. The efficiency of the standard and new procedures is studied
numerically. In these solutions the same input data are used as for constructing high-degree parts of the EGM96 models. From
two sets of
n,m
(n≤360,|m|≤n), derived by the standard and new approaches, different spectral characteristics of the gravity anomaly and the geoid undulation
are estimated and then compared with similar characteristics evaluated by Jekeli's approach (`etalon' solution). The new solution
appears to be very close to Jekeli's, as opposed to the standard solution. The discrepancies between all the characteristics
of the new and `etalon' solutions are smaller than the corresponding discrepancies between two versions of the final geopotential
model EGM96, one of them (HDM190) constructed by the block-diagonal least squares (LS) adjustment and the other one (V068)
by using Jekeli's approach. On the basis of the derived analytical solution a new simple mathematical model is developed to
apply the LS technique for evaluating geopotential coefficients.
Received: 12 December 2000 / Accepted: 21 June 2001 相似文献
12.
The New Hebrides experiment consisted of setting up a pair of DORIS beacons in remote tropical islands in the southwestern
Pacific, between 1993 and 1997. Because of orbitography requirements on TOPEX/Poséidon, the beacons were only transmitting
to SPOT satellites. Root-mean-square (RMS) scatters at the centimeter level on the latitude and vertical components were achieved,
but 2-cm RMS scatters affected the longitude component. Nevertheless, results of relative velocity (123 mm/year N250°) are
very consistent with those obtained using the global positioning system (GPS) (126 mm/yr N246°). The co-seismic step (12 mm
N60°) related to the Walpole event (M
W = 7.7) is consistent with that derived from GPS (10 mm N30°) or from the centroid moment tensor (CMT) of the quake (12 mm
N000°).
Received: 19 November 1999 / Accepted: 17 May 2000 相似文献
13.
A global plate motion model is established based on the ITRF97 velocity fields and geological model NUVEL1. Sub-plate models
are estimated by using the velocity fields derived from 45 global positioning system (GPS) sites under the ITRF97 reference
frame in China. Comparisons between space geodesy and geological models are given. It is found that the Euler vector of the
AFRC–EURA pair has an obvious discrepancy between space geodetic and geological models. The motion patterns of tectonic blocks
in China predicted by GPS are consistent with those of geological data on the whole.
Received: 9 November 2000 / Accepted: 17 September 2001 相似文献
14.
The identification of mean semi-major axes (suitably defined) for satellite orbits to satisfy a variety of requirements for
geodesy, geophysics and oceanography, in terms of repeat orbits (with orbital resonances), is investigated. Various options
for the definition of semi-major axis, from the viewpoint of satellite dynamics, are described. Simple simulations of the
expected resonant changes in inclination are presented, and tools for the analysis of orbit resonances to extract certain
lumped harmonic coefficients of the geopotential (e.g. from the very precise CHAMP orbit) are resurrected. Finally, a preliminary
example of the 46th-order resonance analysis possible for CHAMP, based on the mean orbital elements produced by GFZ (GeoForschungs
Zentrum) for ephemeris prediction, is presented.
Received: 10 July 2001 / Accepted: 17 July 2002
Correspondence to: J. Klokočník at Ondřejov Observatory
Acknowledgements. We thank Prof. Dr. Ch. Reigber, Dr. P. Schwintzer, Dr. T. Gruber and Dr. R. K?nig from GFZ Potsdam for various consultations
and discussions, and for the CHAMP two-line mean elements. This investigation was performed under the aegis of CEDR (Center
for Earth's Dynamics Research, Prague-Ondřejov); it has been supported by project LN00A005 (provided by the Ministry of Education
of the Czech Republic) and by grant A 3004 of the Grant Agency of the Academy of Sciences of the Czech Republic. 相似文献
15.
Analysis of the first year of Earth rotation parameters with a sub-daily resolution gained at the CODE processing center of the IGS 总被引:4,自引:4,他引:0
The solutions of the CODE Analysis Center submitted to the IGS, the International Global Position System (GPS) Service for
Geodynamics, are based on three days of observation of about 80–100 stations of the IGS network. The Earth rotation parameters
(ERPs) are assumed to vary linearly over the three days with respect to an a priori model. Continuity at the day boundaries
as well as the continuity of the first derivatives are enforced by constraints. Since early April 1995 CODE has calculated
a new ERP series with an increased time resolution of 2 hours. Again continuity is enforced at the 2-hours-interval boundaries.
The analysis method is described, particularly how to deal with retrograde diurnal terms in the ERP series which may not be estimated with satellite geodetic methods. The results obtained from the first year of data covered by the time series
(time interval from 4 April 1995 to 30 June 1996) are also discussed. The series is relatively homogeneous in the sense of
the used orbit model and the a priori model for the ERPs. The largest source of excitation at daily and sub-daily periods
is likely to be the effect of the ocean tides. There is good agreement between the present results and Topex/Poseidon ocean
tide models, as well as with models based on Very Long Baseline Interferometry (VLBI) and Satellite Laser Ranging (SLR) data.
Non-oceanic periodic variations are also observed in the series. Their origin is most probably a consequence of the GPS solution
strategy; other possible sources are the atmospheric tides.
Received: 13 July 1999 / Accepted: 21 March 2000 相似文献
16.
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ò 相似文献
17.
Z. Martinec 《Journal of Geodesy》2003,77(1-2):41-49
Three independent gradiometric boundary-value problems (BVPs) with three types of gradiometric data, {Γ
rr
}, {Γ
r
θ,Γ
r
λ} and {Γθθ−Γλλ,Γθλ}, prescribed on a sphere are solved to determine the gravitational potential on and outside the sphere. The existence and
uniqueness conditions on the solutions are formulated showing that the zero- and the first-degree spherical harmonics are
to be removed from {Γ
r
θ,Γ
r
λ} and {Γθθ−Γλλ,Γθλ}, respectively. The solutions to the gradiometric BVPs are presented in terms of Green's functions, which are expressed in
both spectral and closed spatial forms. The logarithmic singularity of the Green's function at the point ψ=0 is investigated
for the component Γ
rr
. The other two Green's functions are finite at this point. Comparisons to the paper by van Gelderen and Rummel [Journal of
Geodesy (2001) 75: 1–11] show that the presented solution refines the former solution.
Received: 3 October 2001 / Accepted: 4 October 2002 相似文献
18.
Global mean sea surface heights (SSHs) and gravity anomalies on a 2′×2′ grid were determined from Seasat, Geosat (Exact Repeat Mission and Geodetic Mission), ERS-1 (1.5-year mean of 35-day, and
GM), TOPEX/POSEIDON (T/P) (5.6-year mean) and ERS-2 (2-year mean) altimeter data over the region 0∘–360∘ longitude and –80∘–80∘ latitude. To reduce ocean variabilities and data noises, SSHs from non-repeat missions were filtered by Gaussian filters
of various wavelengths. A Levitus oceanic dynamic topography was subtracted from the altimeter-derived SSHs, and the resulting
heights were used to compute along-track deflection of the vertical (DOV). Geoidal heights and gravity anomalies were then
computed from DOV using the deflection-geoid and inverse Vening Meinesz formulae. The Levitus oceanic dynamic topography was
added back to the geoidal heights to obtain a preliminary sea surface grid. The difference between the T/P mean sea surface
and the preliminary sea surface was computed on a grid by a minimum curvature method and then was added to the preliminary
grid. The comparison of the NCTU01 mean sea surface height (MSSH) with the T/P and the ERS-1 MSSH result in overall root-mean-square
(RMS) differences of 5.0 and 3.1 cm in SSH, respectively, and 7.1 and 3.2 μrad in SSH gradient, respectively. The RMS differences
between the predicted and shipborne gravity anomalies range from 3.0 to 13.4 mGal in 12 areas of the world's oceans.
Received: 26 September 2001 / Accepted: 3 April 2002
Correspondence to: C. Hwang
Acknowledgements. This research is partly supported by the National Science Council of ROC, under grants NSC89-2611-M-009-003-OP2 and NSC89-2211-E-009-095.
This is a contribution to the IAG Special Study Group 3.186. The Geosat and ERS1/2 data are from NOAA and CERSAT/France, respectively.
The T/P data were provided by AVISO. The CLS and GSFC00 MSS models were kindly provided by NASA/GSFC and CLS, respectively.
Drs. Levitus, Monterey, and Boyer are thanked for providing the SST model. Dr. T. Gruber and two anonymous reviewers provided
very detailed reviews that improved the quality of this paper. 相似文献
19.
A methodology for precise determination of the fundamental geodetic parameter w
0, the potential value of the Gauss–Listing geoid, as well as its time derivative 0, is presented. The method is based on: (1) ellipsoidal harmonic expansion of the external gravitational field of the Earth
to degree/order 360/360 (130 321 coefficients; http://www.uni-stuttgard.de/gi/research/ index.html projects) with respect
to the International Reference Ellipsoid WGD2000, at the GPS positioned stations; and (2) ellipsoidal free-air gravity reduction
of degree/order 360/360, based on orthometric heights of the GPS-positioned stations. The method has been numerically tested
for the data of three GPS campaigns of the Baltic Sea Level project (epochs 1990.8,1993.4 and 1997.4). New w
0 and 0 values (w
0=62 636 855.75 ± 0.21 m2/s2, 0=−0.0099±0.00079 m2/s2 per year, w
0/&γmacr;=6 379 781.502 m,0/&γmacr;=1.0 mm/year, and &γmacr;= −9.81802523 m2/s2) for the test region (Baltic Sea) were obtained. As by-products of the main study, the following were also determined: (1)
the high-resolution sea surface topography map for the Baltic Sea; (2) the most accurate regional geoid amongst four different
regional Gauss–Listing geoids currently proposed for the Baltic Sea; and (3) the difference between the national height datums
of countries around the Baltic Sea.
Received: 14 August 2000 / Accepted: 19 June 2001 相似文献
20.
Ionospheric variation may be considered as a stationary time series under quiet conditions. However, the disturbance of a
stationary random process from stationarity results in the bias of corresponding samples from the stationary observations,
and in the change of statistical model parameters of the process. From a general mathematical aspect, a new method is presented
for monitoring ionospheric variations, based on the characteristic of time-series observation of GPS, and an investigation
of the statistical properties of the estimated auto-covariance of the random ionospheric delay when changing the number of
samples in the time series is carried out. A preliminary scheme for monitoring ionospheric delays is proposed.
Received: 18 August 2000 / Accepted: 12 April 2001 相似文献