首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 546 毫秒
1.
The GEOID96 high-resolution geoid height model for the United States   总被引:4,自引:0,他引:4  
The 2 arc-minute × 2 arc-minute geoid model (GEOID96) for the United States supports the conversion between North American Datum 1983 (NAD 83) ellipsoid heights and North American Vertical Datum 1988 (NAVD 88) Helmert heights. GEOID96 includes information from global positioning system (GPS) height measurements at optically leveled benchmarks. A separate geocentric gravimetric geoid, G96SSS, was first calculated, then datum transformations and least-squares collocation were used to convert from G96SSS to GEOID96. Fits of 2951 GPS/level (ITRF94/NAVD 88) benchmarks to G96SSS show a 15.1-cm root mean square (RMS) around a tilted plane (0.06 ppm, 178 azimuth), with a mean value of −31.4 cm (15.6-cm RMS without plane). This mean represents a bias in NAVD 88 from global mean sea level, remaining nearly constant when computed from subsets of benchmarks. Fits of 2951 GPS/level (NAD 83/NAVD 88) benchmarks to GEOID96 show a 5.5-cm RMS (no tilts, zero average), due primarily to GPS error. The correlated error was 2.5 cm, decorrelating at 40 km, and is due to gravity, geoid and GPS errors. Differences between GEOID96 and GEOID93 range from −122 to +374 cm due primarily to the non-geocentricity of NAD 83. Received: 28 July 1997 / Accepted: 2 September 1998  相似文献   

2.
 Four different implementations of Stokes' formula are employed for the estimation of geoid heights over Sweden: the Vincent and Marsh (1974) model with the high-degree reference gravity field but no kernel modifications; modified Wong and Gore (1969) and Molodenskii et al. (1962) models, which use a high-degree reference gravity field and modification of Stokes' kernel; and a least-squares (LS) spectral weighting proposed by Sj?berg (1991). Classical topographic correction formulae are improved to consider long-wavelength contributions. The effect of a Bouguer shell is also included in the formulae, which is neglected in classical formulae due to planar approximation. The gravimetric geoid is compared with global positioning system (GPS)-levelling-derived geoid heights at 23 Swedish Permanent GPS Network SWEPOS stations distributed over Sweden. The LS method is in best agreement, with a 10.1-cm mean and ±5.5-cm standard deviation in the differences between gravimetric and GPS geoid heights. The gravimetric geoid was also fitted to the GPS-levelling-derived geoid using a four-parameter transformation model. The results after fitting also show the best consistency for the LS method, with the standard deviation of differences reduced to ±1.1 cm. For comparison, the NKG96 geoid yields a 17-cm mean and ±8-cm standard deviation of agreement with the same SWEPOS stations. After four-parameter fitting to the GPS stations, the standard deviation reduces to ±6.1 cm for the NKG96 geoid. It is concluded that the new corrections in this study improve the accuracy of the geoid. The final geoid heights range from 17.22 to 43.62 m with a mean value of 29.01 m. The standard errors of the computed geoid heights, through a simple error propagation of standard errors of mean anomalies, are also computed. They range from ±7.02 to ±13.05 cm. The global root-mean-square error of the LS model is the other estimation of the accuracy of the final geoid, and is computed to be ±28.6 cm. Received: 15 September 1999 / Accepted: 6 November 2000  相似文献   

3.
 Two numerical techniques are used in recent regional high-frequency geoid computations in Canada: discrete numerical integration and fast Fourier transform. These two techniques have been tested for their numerical accuracy using a synthetic gravity field. The synthetic field was generated by artificially extending the EGM96 spherical harmonic coefficients to degree 2160, which is commensurate with the regular 5 geographical grid used in Canada. This field was used to generate self-consistent sets of synthetic gravity anomalies and synthetic geoid heights with different degree variance spectra, which were used as control on the numerical geoid computation techniques. Both the discrete integration and the fast Fourier transform were applied within a 6 spherical cap centered at each computation point. The effect of the gravity data outside the spherical cap was computed using the spheroidal Molodenskij approach. Comparisons of these geoid solutions with the synthetic geoid heights over western Canada indicate that the high-frequency geoid can be computed with an accuracy of approximately 1 cm using the modified Stokes technique, with discrete numerical integration giving a slightly, though not significantly, better result than fast Fourier transform. Received: 2 November 1999 / Accepted: 11 July 2000  相似文献   

4.
 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  相似文献   

5.
In order to study the Baltic Sea Level change and to unify national height systems a two week GPS campaign was performed in the region in Autumn 1990. Parties from Denmark, Finland, Germany, Poland and Sweden carried out GPS measurements at 26 tide gauges along the Baltic sea and 8 VLBI and SLR fiducial stations with baseline lengths ranging from 230 km to 1600 km. The observations were processed in the network mode with the Bernese version 3.3 software using orbit improvement techniques. To get rid of the scale error introduced by the ionospheric refraction from single-frequency data, the local models of the ionosphere were estimated using L4 observations. The tropospheric zenith corrections were also considered. The preliminary results show average root mean square (RMS) errors of about ±3 cm in the horizontal position and ±7 cm in the vertical position relative to the Potsdam SLR station in ITRF89 system. After transformation of the GPS results to geoid heights using the levelled heights, an absolute comparison with gravimetric geoid heights using the least squares modification of Stokes' formula (LSMS), the modified Molodensky and the NKG Scandinavian geoid 1989 (NGK-89) models gives a standard deviation of the difference of ±7cm to ±9cm for the NKG-89 model and of ±9cm to ±30cm for the LSMS and the modified Molodensky model. The Swedish height system is found to be about 8-37cm higher than those of the other Baltic countries for NKG-89 model.  相似文献   

6.
Bathymetry data from Sognefjord, Norway, have been included in a terrain model, and their influence on the geoid has been calculated. The test area, located in the western part of Norway, was chosen due to its deep fjords and high mountains. Inclusion of bathymetry data in the terrain model altered the computed gravimetric geoid by as much as a few decimeters. The effect was detectable to a distance of more than 100 km. All calculated geoids, both with and without bathymetry data in the terrain model, fit the geoidal heights determined by available Global Positioning System (GPS) and levelling heights at the sub-decimetre level. Contrary to expectations, the accuracy in geoid prediction was reduced when using bathymetric data. The geoid changes were largest over the fjord where no GPS points were located. Different methods on the same area [isostatic and Residual Terrain Model (RTM)-terrain reductions] showed differences of approximately 1 m. Rigorous distinction between quasigeoid and geoid was found to be essential in this kind of area. Received: 12 May 1997 / Accepted 7 May 1998  相似文献   

7.
Fast and accurate relative positioning for baselines less than 20 km in length is possible using dual-frequency Global Positioning System (GPS) receivers. By measuring orthometric heights of a few GPS stations by differential levelling techniques, the geoid undulation can be modelled, which enables GPS to be used for orthometric height determination in a much faster and more economical way than terrestrial methods. The geoid undulation anomaly can be very useful for studying tectonic structure. GPS, levelling and gravity measurements were carried out along a 200-km-long highly undulating profile, at an average elevation of 4000 m, in the Ladak region of NW Himalaya, India. The geoid undulation and gravity anomaly were measured at 28 common GPS-levelling and 67 GPS-gravity stations. A regional geoid low of nearly −4 m coincident with a steep negative gravity gradient is compatible with very recent findings from other geophysical studies of a low-velocity layer 20–30 km thick to the north of the India–Tibet plate boundary, within the Tibetan plate. Topographic, gravity and geoid data possibly indicate that the actual plate boundary is situated further north of what is geologically known as the Indus Tsangpo Suture Zone, the traditionally supposed location of the plate boundary. Comparison of the measured geoid with that computed from OSU91 and EGM96 gravity models indicates that GPS alone can be used for orthometric height determination over the Higher Himalaya with 1–2 m accuracy. Received: 10 April 1997 / Accepted: 9 October 1998  相似文献   

8.
The well-known International Association of Geodesy (IAG) approach to the atmospheric geoid correction in connection with Stokes' integral formula leads to a very significant bias, of the order of 3.2 m, if Stokes' integral is truncated to a limited region around the computation point. The derived truncation error can be used to correct old results. For future applications a new strategy is recommended, where the total atmospheric geoid correction is estimated as the sum of the direct and indirect effects. This strategy implies computational gains as it avoids the correction of direct effect for each gravity observation, and it does not suffer from the truncation bias mentioned above. It can also easily be used to add the atmospheric correction to old geoid estimates, where this correction was omitted. In contrast to the terrain correction, it is shown that the atmospheric geoid correction is mainly of order H of terrain elevation, while the term of order H 2 is within a few millimetres. Received: 20 May 1998 / Accepted: 19 April 1999  相似文献   

9.
 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.  相似文献   

10.
 The AUSGeoid98 gravimetric geoid model of Australia has been computed using data from the EGM96 global geopotential model, the 1996 release of the Australian gravity database, a nationwide digital elevation model, and satellite altimeter-derived marine gravity anomalies. The geoid heights are on a 2 by 2 arc-minute grid with respect to the GRS80 ellipsoid, and residual geoid heights were computed using the 1-D fast Fourier transform technique. This has been adapted to include a deterministically modified kernel over a spherical cap of limited spatial extent in the generalised Stokes scheme. Comparisons of AUSGeoid98 with GPS and Australian Height Datum (AHD) heights across the continent give an RMS agreement of ±0.364 m, although this apparently large value is attributed partly to distortions in the AHD. Received: 10 March 2000 / Accepted: 21 February 2001  相似文献   

11.
The short wavelength geoid undulations, caused by topography, amount to several decimeters in mountainous areas. Up to now these effects are computed by means of digital terrain models in a grid of 100–500m. However, for many countries these data are not yet available or their collection is too expensive. This problem can be overcome by considering the special behaviour of the gravity potential along mountain slopes. It is shown that 90 per cent of the topographic effects are represented by a simple summation formula, based on the average height differences and distances between valleys and ridges along the geoid profiles, δN=[30.H.D.+16.(H−H′).D] in mm/km, (error<10%), whereH, H′, D are estimated in a map to the nearest 0.2km. The formula is valid for asymmetric sides of valleys (H, H′) and can easily be corrected for special shapes. It can be used for topographic refinement of low resolution geoids and for astrogeodetic projects. The “slope method” was tested in two alpine areas (heights up to 3800m, astrogeodetic deflection points every 170km 2) and resulted in a geoid accuracy of ±3cm. In first order triangulation networks (astro points every 1000km 2) or for gravimetric deflections the accuracy is about 10cm per 30km. Since a map scale of 1∶500.000 is sufficient, the method is suitable for developing countries, too.  相似文献   

12.
A new gravimetric geoid model, USGG2009 (see Abbreviations), has been developed for the United States and its territories including the Conterminous US (CONUS), Alaska, Hawaii, Guam, the Commonwealth of the Northern Mariana Islands, American Samoa, Puerto Rico and the US Virgin Islands. USGG2009 is based on a 1′ × 1′ gravity grid derived from the NGS surface gravity data and the DNSC08 altimetry-derived anomalies, the SRTM-DTED1 3′′ DEM for its topographic reductions, and the global geopotential model EGM08 as a reference model. USGG2009 geoid heights are compared with control values determined at 18,398 Bench Marks over CONUS, where both the ellipsoidal height above NAD 83 and the Helmert orthometric height above NAVD 88 are known. Correcting for the ellipsoidal datum difference, this permits a comparison of the geoid heights to independent data. The standard deviation of the differences is 6.3 cm in contrast to 8.4 cm for its immediate predecessor— USGG2003. To minimize the effect of long-wavelength errors that are known to exist in NAVD88, these comparisons were made on a state-by-state basis. The standard deviations of the differences range from 3–5 cm in eastern states to about 6–9 cm in the more mountainous western states. If the GPS/Bench Marks-derived geoid heights are corrected by removing a GRACE-derived estimate of the long-wavelength NAVD88 errors before the comparison, the standard deviation of their differences from USGG2009 drops to 4.3 cm nationally and 2–4 cm in eastern states and 4–8 in states with a maximum error of 26.4 cm in California and minimum of −32.1 cm in Washington. USGG2009 is also compared with geoid heights derived from 40 tide-gauges and a physical dynamic ocean topography model in the Gulf of Mexico; the mean of the differences is 3.3 cm and their standard deviation is 5.0 cm. When USGG2009-derived deflections of the vertical are compared with 3,415 observed surface astro-geodetic deflections, the standard deviation of the differences in the N–S and E–W components are 0.87′′ and 0.94′′, respectively.  相似文献   

13.
The second Baltic Sea Level (BSL) GPS campaign was run for one week in June 1993. Data from 35 tide gauge sites and five fiducial stations were analysed, for three fiducial stations (Onsala, Mets?hovi and Wettzell) fixed at the ITRF93 system. On a time-scale of 5 days, precision was several parts in 109 for the horizontal and vertical components. Accuracies were about 1 cm in comparison with the International GPS Geodynamical Service (IGS) coordinates in three directions. To connect the Swedish and the Finnish height systems, our numerical application utilises three approaches: a rigorous approach, a bias fit and a three-parameter fit. The results between the Swedish RH70 and the Finnish N 60 systems are estimated to −19.3 ± 6.5, −17 ± 6 and −15 ± 6 cm, respectively, by the three approaches. The results of the three indirect methods are in an agreement with those of a direct approach from levelling and gravity measurements. Received: 3 April 1996 / Accepted: 4 August 1997  相似文献   

14.
In precise geoid determination by Stokes formula, direct and primary and secondary indirect terrain effects are applied for removing and restoring the terrain masses. We use Helmert's second condensation method to derive the sum of these effects, together called the total terrain effect for geoid. We develop the total terrain effect to third power of elevation H in the original Stokes formula, Earth gravity model and modified Stokes formula. It is shown that the original Stokes formula, Earth gravity model and modified Stokes formula all theoretically experience different total terrain effects. Numerical results indicate that the total terrain effect is very significant for moderate topographies and mountainous regions. Absolute global mean values of 5–10 cm can be reached for harmonic expansions of the terrain to degree and order 360. In another experiment, we conclude that the most important part of the total terrain effect is the contribution from the second power of H, while the contribution from the third power term is within 9 cm. Received: 2 September 1996 / Accepted: 4 August 1997  相似文献   

15.
 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  相似文献   

16.
A synthetic [simulated] Earth gravity model (SEGM) of the geoid, gravity and topography has been constructed over Australia specifically for validating regional gravimetric geoid determination theories, techniques and computer software. This regional high-resolution (1-arc-min by 1-arc-min) Australian SEGM (AusSEGM) is a combined source and effect model. The long-wavelength effect part (up to and including spherical harmonic degree and order 360) is taken from an assumed errorless EGM96 global geopotential model. Using forward modelling via numerical Newtonian integration, the short-wavelength source part is computed from a high-resolution (3-arc-sec by 3-arc-sec) synthetic digital elevation model (SDEM), which is a fractal surface based on the GLOBE v1 DEM. All topographic masses are modelled with a constant mass-density of 2,670 kg/m3. Based on these input data, gravity values on the synthetic topography (on a grid and at arbitrarily distributed discrete points) and consistent geoidal heights at regular 1-arc-min geographical grid nodes have been computed. The precision of the synthetic gravity and geoid data (after a first iteration) is estimated to be better than 30 μ Gal and 3 mm, respectively, which reduces to 1 μ Gal and 1 mm after a second iteration. The second iteration accounts for the changes in the geoid due to the superposed synthetic topographic mass distribution. The first iteration of AusSEGM is compared with Australian gravity and GPS-levelling data to verify that it gives a realistic representation of the Earth’s gravity field. As a by-product of this comparison, AusSEGM gives further evidence of the north–south-trending error in the Australian Height Datum. The freely available AusSEGM-derived gravity and SDEM data, included as Electronic Supplementary Material (ESM) with this paper, can be used to compute a geoid model that, if correct, will agree to in 3 mm with the AusSEGM geoidal heights, thus offering independent verification of theories and numerical techniques used for regional geoid modelling.Electronic Supplementary Material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s00190-005-0002-z  相似文献   

17.
The objective of this study is to evaluate two approaches, which use different representations of the Earth’s gravity field for downward continuation (DC), for determining Helmert gravity anomalies on the geoid. The accuracy of these anomalies is validated by 1) analyzing conformity of the two approaches; and 2) converting them to geoid heights and comparing the resulting values to GPS-leveling data. The first approach (A) consists of evaluating Helmert anomalies at the topography and downward-continuing them to the geoid. The second approach (B) downward-continues refined Bouguer anomalies to the geoid and transforms them to Helmert anomalies by adding the condensed topographical effect. Approach A is sensitive to the DC because of the roughness of the Helmert gravity field. The DC effect on the geoid can reach up to 2 m in Western Canada when the Stokes kernel is used to convert gravity anomalies to geoid heights. Furthermore, Poisson’s equation for DC provides better numerical results than Moritz’s equation when the resulting geoid models are validated against the GPS-leveling. On the contrary, approach B is significantly less sensitive to the DC because of the smoothness of the refined Bouguer gravity field. In this case, the DC (Poisson’s and Moritz’s) contributes only at the decimeter level to the geoid model in Western Canada. The maximum difference between the geoid models from approaches A and B is about 5 cm in the region of interest. The differences may result from errors in the DC such as numerical instability. The standard deviations of the hHN for both approaches are about 8 cm at the 664 GPS-leveling validation stations in Western Canada.  相似文献   

18.
 The analysis of lunar laser ranging (LLR) data enables the determination of many parameters of the Earth–Moon system, such as lunar gravity coefficients, reflector and station coordinates which contribute to the realisation of the International Terrestrial Reference Frame 2000 (ITRF 2000), Earth orientation parameters [EOPs, which contribute to the global EOP solutions at the International Earth Rotation Service (IERS)] or quantities which parameterise relativistic effects in the solar system. The big advantage of LLR is the long time span of lunar observations (1970–2000). The accuracy of the normal points nowadays is about 1 cm.  The capability of LLR to determine tidal parameters is investigated. In principle, it could be assumed that LLR would contribute greatly to the investigation of tidal effects, because the Moon is the most important tide-generating body. In this respect some special topics such as treatment of the permanent tide and the effect of atmospheric loading are addressed and results for the tidal parameters h 2 and l 2 as well as values for the eight main tides are given. Received: 14 August 2000 / Accepted: 15 October 2001  相似文献   

19.
Improvements in height datum transfer expected from the GOCE mission   总被引:1,自引:1,他引:1  
 One of the aims of the Earth Explorer Gravity Field and Steady-State Ocean Circulation (GOCE) mission is to provide global and regional models of the Earth's gravity field and of the geoid with high spatial resolution and accuracy. Using the GOCE error model, simulation studies were performed in order to estimate the accuracy of datum transfer in different areas of the Earth. The results showed that with the GOCE error model, the standard deviation of the height anomaly differences is about one order of magnitude better than the corresponding value with the EGM96 error model. As an example, the accuracy of the vertical datum transfer from the tide gauge of Amsterdam to New York was estimated equal to 57 cm when the EGM96 error model was used, while in the case of GOCE error model this accuracy was increased to 6 cm. The geoid undulation difference between the two places is about 76.5 m. Scaling the GOCE errors to the local gravity variance, the estimated accuracy varied between 3 and 7 cm, depending on the scaling model. Received: 1 March 2000 / Accepted: 21 February 2001  相似文献   

20.
About half a million marine gravity measurements over a 30×30 area centered on Japan have been processed and adjusted to produce a new free-air gravity map from a 5′×5′ grid. This map seems to have a better resolution than those previously published as measured by its correlation with bathymetry. The grid was used together with a high-degree and -order spherical harmonics geopotential model to compute a detailed geoid with two methods: Stokes integral and collocation. Comparisons with other available geoidal surfaces derived either from gravity or from satellite altimetry were made especially to test the ability of this new geoid at showing the sea surface topography as mapped by the Topex/Poseidon satellite. Over 2 months (6 cycles) the dynamic topography at ascending passes in the region (2347N and 123147E) was mapped to study the variability of the Kuroshio current. Received: 15 July 1994 / Accepted: 17 February 1997  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号