共查询到20条相似文献，搜索用时 15 毫秒
1.
The Scripps Orbit and Permanent Array Center (SOPAC) has completed development for the UNAVCO community of first-generation GPS Seamless Archive (GSAC) software. The GSAC is a virtual archive composed of an assembly of agencies and investigators exchanging information about their respective GPS-related data holdings in a well defined, cohesive manner. The superset of this published information is collected and ingested into centralized databases administered currently by two data brokers (Retailers), who make the data available to the public in a seamless manner. There are three user interfaces available: the interactive GSAC Wizard, a command-line Unix-style executable called gsac-client, and a front door HTTP service called the GSAC Retailer Service Interface. Each user interface provides access to the data collections of 6 different GPS archives (GSAC Wholesalers) in North America. Together these archives have published more than 2 million GPS data files pertaining to over 10,000 different geodetic monuments. These datasets are composed in large part of data collected by US scientists and their collaborators over the period 1986 to the present in Western North America and other tectonically active regions around the globe, as well as the holdings of two IGS global data centers. In this article, we describe how the three GSAC user interfaces provide the community a powerful set of tools for seamlessly mining information and collecting data files from a distributed network of GPS archives.The GPS Toolbox is a column dedicated to highlighting algorithms and source code utilized by GPS Engineers and scientists. If you have an interesting program or software package you would like to share with our readers, please pass it along; e-mail it to us at gps-toolbox@ngs.noaa.gov. To comment on any of the source code discussed here, or to download source code, visit our website at . This column is edited by Stephen Hilla, National Geodetic Survey, NOAA, Silver Spring, Maryland, and Mike Craymer, Geodetic Survey Division, Natural Resources Canada, Ottawa, Ontario, Canada. 相似文献
2.
A new theory for high-resolution regional geoid computation without applying Stokess formula is presented. Operationally, it uses various types of gravity functionals, namely data of type gravity potential (gravimetric leveling), vertical derivatives of the gravity potential (modulus of gravity intensity from gravimetric surveys), horizontal derivatives of the gravity potential (vertical deflections from astrogeodetic observations) or higher-order derivatives such as gravity gradients. Its algorithmic version can be described as follows: (1) Remove the effect of a very high degree/order potential reference field at the point of measurement (POM), in particular GPS positioned, either on the Earths surface or in its external space. (2) Remove the centrifugal potential and its higher-order derivatives at the POM. (3) Remove the gravitational field of topographic masses (terrain effect) in a zone of influence of radius r. A proper choice of such a radius of influence is 2r=4×104 km/n, where n is the highest degree of the harmonic expansion. (cf. Nyquist frequency). This third remove step aims at generating a harmonic gravitational field outside a reference ellipsoid, which is an equipotential surface of a reference potential field. (4) The residual gravitational functionals are downward continued to the reference ellipsoid by means of the inverse solution of the ellipsoidal Dirichlet boundary-value problem based upon the ellipsoidal Abel–Poisson kernel. As a discretized integral equation of the first kind, downward continuation is Phillips–Tikhonov regularized by an optimal choice of the regularization factor. (5) Restore the effect of a very high degree/order potential reference field at the corresponding point to the POM on the reference ellipsoid. (6) Restore the centrifugal potential and its higher-order derivatives at the ellipsoidal corresponding point to the POM. (7) Restore the gravitational field of topographic masses ( terrain effect) at the ellipsoidal corresponding point to the POM. (8) Convert the gravitational potential on the reference ellipsoid to geoidal undulations by means of the ellipsoidal Bruns formula. A large-scale application of the new concept of geoid computation is made for the Iran geoid. According to the numerical investigations based on the applied methodology, a new geoid solution for Iran with an accuracy of a few centimeters is achieved.Acknowledgments. The project of high-resolution geoid computation of Iran has been support by National Cartographic Center (NCC) of Iran. The University of Tehran, via grant number 621/3/602, supported the computation of a global geoid solution for Iran. Their support is gratefully acknowledged. A. Ardalan would like to thank Mr. Y. Hatam, and Mr. K. Ghazavi from NCC and Mr. M. Sharifi, Mr. A. Safari, and Mr. M. Motagh from the University of Tehran for their support in data gathering and computations. The authors would like to thank the comments and corrections made by the four reviewers and the editor of the paper, Professor Will Featherstone. Their comments helped us to correct the mistakes and improve the paper. 相似文献
3.
The value of ocean reflections of GPS signals to enhance satellite altimetry: data distribution and error analysis 总被引：4，自引：0，他引：4
The time and space distribution of general reflection altimetry from two satellites (senders and receivers) via the the oceans surface is examined with specific reference to GPS senders and two current receiving satellites. While a considerable enhancement of conventional altimeter coverage is possible in all configurations if the reflection signals can be used, repeating passes of these (with GPS senders) having reasonably small cycle times (days to tens of days) occur only if the receiving orbit is nearly polar. Results of an analysis of the fundamental geometry show that over a large range of reflection angles the error of recovered sea heights depends almost entirely on the errors in the delay signal and the radial error of the receiving satellite (using current estimates of GPS orbit accuracies). The most critical element is the precision of the delay measurement. Both it and the accuracy of the receiving orbit should be below the decimeter level for the technique to achieve its full potential. 相似文献
4.
This study makes an initial comparison of three GPS-like constellations. Starting with a simplified constellation of 25 GPS satellites as a reference, GPS(25), we determine what kinematic positioning improvements would result from a constellation comprising a Hi component of 16 GPS satellites (at roughly 16.8 earth radii) coupled with a Lo component of 49 GPS satellites (at roughly 2.1 earth radii). We also include a GPS constellation of 49 GPS satellites, GPS(49), which comprises orbits like the GPS(25) constellation. The GPS(49) and the Hi(16)/Lo(49) constellations have semi-major axes selected so that they have exactly the same average number of satellites above 7.5 degrees elevation (averaged over 24 hours). What motivated this study was a need to measure the benefits, to precision differential kinematic positioning methods (i.e., RTK), which result from the higher Doppler shifts (hence speedier integrated Doppler) generated by the Lo component. Quicker initial convergence was anticipated, of course. 相似文献
5.
On Helmert’s methods of condensation 总被引：2，自引：0，他引：2
B. Heck 《Journal of Geodesy》2003,77(3-4):155-170
Helmerts first and second method of condensation are reviewed and generalized in two respects: First, the point at which the effects of topographical and condensation masses are calculated may be situated on or outside the topographical surface; second, the depth of the condensation layer below the geoid is arbitrary. While the first extension permits the application of the generalized model to the evaluation of airborne and satellite data, the second one gives an additional degree of freedom which can be used to provide a smooth gravity field after reducing the observation data. The respective formulae are derived for the generalized condensation model in both planar and spherical approximation. A comparison of the planar and the spherical model shows some structural differences, which are primarily visible in the out-of-integral terms. Considering the respective formulae for the combined topographic–condensation reduction on the background of the density structure of the Earths lithosphere, the consequences for the residual gravity field are investigated; it is shown that the residual field after applying Helmerts second model of reduction is very rough, making this procedure unfavourable for downward continuation. Further considerations refer to the question of which sets of formulae should be used in geoid and quasigeoid determination. It is concluded that for high-precision applications the generalized spherical model, involving a depth of the condensation layer of between 20 and 30 km, should be superior to Helmerts second model of condensation, although it requires the direct calculation of the indirect effect, which is larger than in the case of Helmerts second method of condensation. 相似文献
7.
Summary The signals transmitted by Block II satellites of the Global Positioning System (GPS) can be degraded to limit the highest accuracy of the system (10 m or better point positioning) to authorized users. This mode of degraded operation is called Selective Availability (S/A). S/A involves the degradation in the quality of broadcast orbits and satellite clock dithering. We monitored the dithered satellite oscillator and investigated the effect of this clock dithering on high accuracy relative positioning. The effect was studied over short 3-meter and zero-baselines with two GPS receivers. The equivalent S/A effects for baselines ranging from 0 to >10,000 km can be examined with short test baselines if the receiver clocks are deliberately mis-synchronized by a known and varying amount. Our results show that the maximum effect of satellite clock dithering on GPS double difference phase residuals grows as a function of the clock synchronization error according to: S/A
effect
=0.04 cm/msec, and it increases as a function of baseline length like: S/A
effect
=0.014 cm/100 km. These are equations for maximum observed values of post-fit residuals due to S/A. The effect on GPS baselines is likely to be smaller than the 0.14 mm for a baseline separation of 100 km. We therefore conclude, for our limited data set, and for the level of S/A during our tests, that S/A clock dithering has negligible effect on all terrestrial GPS baselines if double difference processing techniques are employed and if the GPS receivers remain synchronized to better than 10 msec. S/A may constitute a problem, however, if accurate point processing is required, or if GPS receivers are not synchronized. We suggest and test two different methods to monitor satellite frequency offsets due to S/A. S/A modulates GPS carrier frequencies in the range of-2 Hz to +2 Hz over time periods of several minutes. The methods used in this paper to measure the satellite clock dither could be applied by the civilian GPS community to continuously monitor S/A clock dithering. The monitored frequencies may aid high accuracy point positioning applications in a postprocessing mode (Malys and Ortiz 1989), and differential GPS with poorly synchronized receivers (Feigl et al. 1991). 相似文献
8.
Baardas reliability measures for outliers, as well as sensitivity and separability measures for deformations, are functions of the lower bound of the non-centrality parameter (LBNP). This parameter, which is taken from Baardas well-known nomograms, is actually a non-centrality parameter of the cumulative distribution function (CDF) of the non-central 2-distribution yielding a complementary probability of the desired power of the test, i.e. probability of Type II error. It is investigated how the LBNP can be computed for desired probabilities (power of the test and significance level) and known degrees of freedom. Two recursive algorithms, namely bisection and the Newton algorithm, were applied to compute the LBNP after the definition of a stable and accurate algorithm for the computation of the corresponding CDF. Despite the fact that the recursive algorithms ensure some desired accuracy, it is presented numerically that the Newton algorithm has a faster convergence to the solution than the bisection algorithm. 相似文献
9.
The six-hourly values of the atmospheric angular momentum (AAM) functions computed by the U.S. National Meteorological Center (NMC) were used to estimate the effects of the atmospheric tides on the Earth's rotation. Variations of the equatorial components
1 and
2 of the AAM have periods close to gravitational tidesP
1 andK
1.The amplitudes of the detected variations in
1 and
2 functions have been found to be much larger than the theoretical ones, the reason of this amplification remains unexplained. According to theoretical formulations, these waves can be expressed only as retrograde motions. Because of frame effects, there is a correspondance between diurnal retrograde polar motion and precession-nutations and the atmospheric effect on polar motion cannot be detected from observations.The second part of this paper deals the effects of atmospheric tides in Earth rotation. High-frequency UT1 variations have been derived from VLBI and GPS techniques during the SEARCH'92 campaign (Study ofEarth-AtmosphereRapidCHanges) (Dickey et al. 1994). They have been compared to values derived by Ray et al. (1994) from global ocean tide model. The results obtained in the present paper show the existence of variations of thermal origin with an amplitude of about 1µs in Universal Time UT1. The agreement between observed and theoretical values is better when the determined thermal atmospheric tides are taken into account.Oceanic tidal signal explains a large part (60% of the signal variance) of the diurnal and sub-diurnal variations. Our results show that only a small part of the residuals (5%) accounts for the atmospheric tidal effects. The residual signal remains unexplained; it might be due to mismodelization of oceanic or atmospheric tides or effect of other geophysical phenomena. 相似文献
10.
Jussi Kääriäinen 《Journal of Geodesy》1992,66(3):281-283
Two water tube clinometers are installed in a deep mine in Lohja, Finland. One, lying in an E-W direction, is 177 m long and the other, lying in a N-S direction, is 62 m long. These instruments are used for tidal records and research. Because the station was on the zone of totality of the 1990 solar eclipse, the instruments offered a suitable tool for investigating whether there was any shielding effect on the direction of the vertical. During the eclipse, the resolution of the instruments was increased to 0.00001 and 0.00004, respectively. p ]Recordings and data evaluation are discussed, and the results obtained show that no gravitational shielding was found at the level of the above accuracy. 相似文献
11.
Ernst Heinrich Knickmeyer 《Journal of Geodesy》1994,68(2):71-76
The inspection of a pipeline from the inside becomes more and more widespread. An inspection tool, called Geometry Pig, carries a variety of sensors, among which a strapdown inertial system provides essential geometric information on the pipeline. In this paper the inspection pig is described as a surveying tool, which delivers information on centre line, shape and features of the pipe. After some background information it is explained, why the inertial system is the optimal instrument to deliver some structural parameters of the pipeline. The interaction between data evaluation and measurement procedure is discussed in view of the very few possibilities to vary the procedure. One of them is a forced rotation about the forward axis. Analogies are shown to traditional surveying procedures. Finally some results of actual surveys are presented. A positioning accuracy of 1:2500 is achieved with respect to tie points. 相似文献
12.
Time variations in the Earths gravity field at periods longer than 1 year, for degree-two spherical harmonics, C21, S21, and C20, are estimated from accurately measured Earth rotational variations. These are compared with predictions of atmospheric, oceanic, and hydrologic models, and with independent satellite laser ranging (SLR) results. There is remarkably good agreement between Earth rotation and model predictions of C21 and S21 over a 22-year period. After decadal signals are removed, Earth-rotation-derived interannual C20 variations are dominated by a strong oscillation of period about 5.6 years, probably due to uncertainties in wind and ocean current estimates. The model-predicted C20 agrees reasonably well with SLR observations during the 22-year period, with the exception of the recent anomaly since 1997/1998. 相似文献
13.
The Center for Orbit Determination in Europe (CODE) has been involved in the processing of combined GPS/GLONASS data during the International GLONASS Experiment (IGEX). The resulting precise orbits were analyzed using the program SORBDT. Introducing one satellites positions as pseudo-observations, the program is capable of fitting orbital arcs through these positions using an orbit improvement procedure based on the numerical integration of the satellites orbit and its partial derivative with respect to the orbit parameters. For this study, the program was enhanced to estimate selected parameters of the Earths gravity field. The orbital periods of the GPS satellites are —in contrast to those of the GLONASS satellites – 2:1 commensurable (P
Sid:P
GPS) with the rotation period of the Earth. Therefore, resonance effects of the satellite motion with terms of the geopotential occur and they influence the estimation of these parameters. A sensitivity study of the GPS and GLONASS orbits with respect to the geopotential coefficients reveals that the correlations between different geopotential coefficients and the correlations of geopotential coefficients with other orbit parameters, in particular with solar radiation pressure parameters, are the crucial issues in this context. The estimation of the resonant geopotential terms is, in the case of GPS, hindered by correlations with the simultaneously estimated radiation pressure parameters. In the GLONASS case, arc lengths of several days allow the decorrelation of the two parameter types. The formal errors of the estimates based on the GLONASS orbits are a factor of 5 to 10 smaller for all resonant terms.
AcknowledgmentsThe authors would like to thank all the organizations involved in the IGS and the IGEX campaign, in particular those operating an IGS or IGEX observation site and providing the indispensable data for precise orbit determination. 相似文献
14.
Wieslaw Kosek 《Journal of Geodesy》1993,67(1):1-9
In order to find short periodic oscillations in the Earth's rate of rotation, atmospheric angular momentum, solar activity the Maximum Entropy Spectral Analysis — MESA (Burg, 1967) has been applied. The MESA with moving autoregressive order has been introduced in order to detect more accurately periods of very weak short periodic variations. Oscillations with periods of about 75, 50, 27 and 18 days have been found in length of day — LOD, from which tidal oscillations were removed up to 35 days — LODR computed by the Center for Space Research — CSR from Lageos Laser Ranging data, in the axial component of atmospheric angular momentum —
3 determined by the U.S. National Meteorological Center — NMC and in the geomagnetic activity represented by the geomagnetic index —A
p (Lincoln, 1967). These oscillations computed by Ormsby band pass filter (Ormsby, 1961) are in a very good phase agreement in the case of oscillations with periods of 50 and 18 days in these 3 series. The MESA of the cross covariance estimations between LODR-
3,
3-A
p,A
p-LODR, LODR-FLUX,
3-FLUX, andA
p-FLUX has confirmed the existence of common oscillations with periods of 70, 50, 27 and 18 days. This indicates a possible relationship between solar activity and the short periodic exchange of angular momentum between the atmosphere and the solid Earth. 相似文献
15.
The treatment of the permanent tidal deformation of the Earth in GPS computation has been an almost unmentioned topic in the GPS literature. However, the ever increasing accuracy and the need to combine the GPS based coordinates with other methods requires a consistent way to handle the tides. Our survey shows that both the ITRF-xx coordinates and the GPS based coordinates are nowadays reduced to a non-tidal crust, conventionally defined using physically meaningless parameters. We propose to use instead the zero-crust concept which corresponds to concepts already accepted in the resolution of IAG in 1983 for gravimetric works. 相似文献
16.
J. Klokoník J. Kostelecký C.A. Wagner P. Schwintzer Ch. Förste R. Scharroo 《Journal of Geodesy》2005,78(7-8):405-417
Since the advent of CHAMP, the first in a series of low-altitude satellites being almost continuously and precisely tracked by GPS, a new generation of long-wavelength gravitational geopotential models can be derived. The accuracy evaluation of these models depends to a large extent on the comparison with external data of comparable quality. Here, two CHAMP-derived models, EIGEN-1S and EIGEN-2, are tested with independent long-term-averaged single satellite crossover (SSC) sea heights from three altimetric satellites (ERS-1, ERS-2 and Geosat). The analyses show that long-term averages of crossover residuals still are powerful data to test CHAMP gravity field models. The new models are tested in the spatial domain with the aid of ERS-1/-2 and Geosat SSCs, and in the spectral domain with latitude-lumped coefficient (LLC) corrections derived from the SSCs. The LLC corrections allow a representation of the satellite-orbit-specific error spectra per order of the models spherical harmonic coefficients. These observed LLC corrections are compared to the LLC projections from the models variance–covariance matrix. The excessively large LLC errors at order 2 found in the case of EIGEN-2 with the ERS data are discussed. The degree-dependent scaling factors for the variance-covariance matrices of EIGEN-1S and –2, applied to obtain more realistic error estimates of the solved-for coefficients, are compatible with the results found here. 相似文献
17.
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. 相似文献
18.
A new, high-resolution and high-precision geoid has been computed for the whole of Canada and part of the U.S., ranging from 35°N to about 90°N in latitude and 210°E to 320°E in longitude. The OSU91A geopotential model complete to degree and order 360 was combined with a 5 × 5 mean gravity anomaly grid and 1km × 1km topographical information to generate the geoid file. The remove-restore technique was adopted for the computation of terrain effects by Helmert's condensation reduction. The contribution of the local gravity data to the geoid was computed strictly by the 1D-FFT technique, which allows for the evaluation of the discrete spherical Stokes integral without any approximation, parallel by parallel. The indirect effects of up to second order were considered. The internal precision of the geoid, i.e. the contribution of the gravity data and the model coefficients noise, was also evaluated through error propagation by FFT. In a relative sense, these errors seem to agree quite well with the external errors and show clearly the weak areas of the geoid which are mostly due to insufficient gravity data coverage. Comparison of the gravimetric geoid with the GPS/levelling-derived geoidal heights of eight local GPS networks with a total of about 900 stations shows that the absolute agreement with respect to the GPS/levelling datum is generally better than 10 cm RMS and the relative agreement ranges, in most cases, from 4 to 1 ppm over short distances of about 20 to 100km, 1 to 0.5 ppm over distances of about 100 to 200 km, and 0.5 to 0.1 ppm for baselines of 200 to over 1000 km. Other existing geoids, such as UNB90, GEOID90 and GSD91, were also included in the comparison, showing that the new geoid achieves the best agreement with the GPS/levelling data.Presented at theIAG General Meeting, Beijing, P.R. China, Aug. 6–13, 1993 相似文献
19.
A method based on multi-antennae linked to a common GPS receiver is proposed. The goal of the technique is to improve height determination for baselines a few kilometres in length. The advantage of this technique resides in the elimination of relative clock parameters in the between-antenna single difference observations. Because single difference observations are free of clock errors more geometrical strength remains to determine the baseline components. This statement is valid as long as intercable biases can be carefully calibrated. For millimetre height determination, the intercable calibration must be done at the same level of accuracy. Under this assumption it is shown that in general the height standard deviation improves by a factor of about three compared to standard GPS data processing. With the proposed method, the effect of relative tropospheric zenith delay errors becomes a bit smaller (in absolute value), compared to standard data processing. To absorb this error, a relative tropospheric zenith delay parameter may be estimated. Even with this additional parameter in the solution the height standard deviation remains two times smaller than the results of standard processing techniques (without tropospheric zenith delay parameter), and at least five times smaller than in the results obtained from standard processing including one tropospheric zenith delay parameter. 相似文献
20.
Jürgen Kusche 《Journal of Geodesy》1996,70(4):183-187
This paper deals with position estimation and path control for Autonomous Guided Vehicles (AGV). To enable a vehicle or a mobile robot in following a continuous virtual path without human control, these techniques play an important role. The relationship between the vehicle's motion in 3-dimensional space and the shape of a curved surface is described. In particular, the introduction of a digital terrain model in dead reckoning is considered. Moreover, a possible nonlinear control is developed based on curvilinear path coordinates, and the proof for global stability is given. To achieve general validity, these topics are treated here independently of the cart's special mechanization (the configuration of steered wheels and driven wheels). Simulation studies are presented to illustrate the investigations. 相似文献