首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 533 毫秒
1.
In a modern application of Stokes formula for geoid determination, regional terrestrial gravity is combined with long-wavelength gravity information supplied by an Earth gravity model. Usually, several corrections must be added to gravity to be consistent with Stokes formula. In contrast, here all such corrections are applied directly to the approximate geoid height determined from the surface gravity anomalies. In this way, a more efficient workload is obtained. As an example, in applications of the direct and first and second indirect topographic effects significant long-wavelength contributions must be considered, all of which are time consuming to compute. By adding all three effects to produce a combined geoid effect, these long-wavelength features largely cancel. The computational scheme, including two least squares modifications of Stokes formula, is outlined, and the specific advantages of this technique, compared to traditional gravity reduction prior to Stokes integration, are summarised in the conclusions and final remarks. AcknowledgementsThis paper was written whilst the author was a visiting scientist at Curtin University of Technology, Perth, Australia. The hospitality and fruitful discussions with Professor W. Featherstone and his colleagues are gratefully acknowledged.  相似文献   

2.
Gravimetric geoid determination by Stokes formula requires that the effects of topographic masses be removed prior to Stokes integration. This step includes the direct topographic and the downward continuation (DWC) effects on gravity anomaly, and the computations yield the co-geoid height. By adding the effect of restoration of the topography, the indirect effect on the geoid, the geoid height is obtained. Unfortunately, the computations of all these topographic effects are hampered by the uncertainty of the density distribution of the topography. Usually the computations are limited to a constant topographic density, but recently the effects of lateral density variations have been studied for their direct and indirect effects on the geoid. It is emphasised that the DWC effect might also be significantly affected by a lateral density variation. However, instead of computing separate effects of lateral density variation for direct, DWC and indirect effects, it is shown in two independent ways that the total geoid effect due to the lateral density anomaly can be represented as a simple correction proportional to the lateral density anomaly and the elevation squared of the computation point. This simple formula stems from the fact that the significant long-wavelength contributions to the various topographic effects cancel in their sum. Assuming that the lateral density anomaly is within 20% of the standard topographic density, the derived formula implies that the total effect on the geoid is significant at the centimetre level for topographic elevations above 0.66 km. For elevations of 1000, 2000 and 5000 m the effect is within ± 2.2, ± 8.8 and ± 56.8 cm, respectively. For the elevation of Mt. Everest the effect is within ± 1.78 m.  相似文献   

3.
This study emphasizes that the harmonic downward continuation of an external representation of the Earth’s gravity potential to sea level through the topographic masses implies a topographic bias. It is shown that the bias is only dependent on the topographic density along the geocentric radius at the computation point. The bias corresponds to the combined topographic geoid effect, i.e., the sum of the direct and indirect topographic effects. For a laterally variable topographic density function, the combined geoid effect is proportional to terms of powers two and three of the topographic height, while all higher order terms vanish. The result is useful in geoid determination by analytical continuation, e.g., from an Earth gravity model, Stokes’s formula or a combination thereof.  相似文献   

4.
Today the combination of Stokes formula and an Earth gravity model (EGM) for geoid determination is a standard procedure. However, the method of modifying Stokes formula varies from author to author, and numerous methods of modification exist. Most methods modify Stokes kernel, but the most widely applied method, the remove compute restore technique, removes the EGM from the gravity anomaly to attain a residual gravity anomaly under Stokes integral, and at least one known method modifies both Stokes kernel and the gravity anomaly. A general model for modifying Stokes formula is presented; it includes most of the well-known techniques of modification as special cases. By assuming that the error spectra of the gravity anomalies and the EGM are known, the optimum model of modification is derived based on the least-squares principle. This solution minimizes the expected mean square error (MSE) of all possible solutions of the general geoid model. A practical formula for estimating the MSE is also presented. The power of the optimum method is demonstrated in two special cases. AcknowledgementsThis paper was partly written whilst the author was a visiting scientist at The University of New South Wales, Sydney, Australia. He is indebted to Professor W. Kearsley and his colleagues, and their hospitality is acknowledged.  相似文献   

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

6.
Recent papers in the geodetic literature promote the reduction of gravity for geoid determination according to the Helmert condensation technique where the entire reduction is made in place before downward continuation. The alternative approach, primarily developed by Moritz, uses two evaluation points, one at the Earths surface, the other on the (co-)geoid, for the direct topographic effect. Both approaches are theoretically legitimate and the derivations in each case make use of the planar approximation and a Lipschitz condition on height. Each method is re-formulated from first principles, yielding equations for the direct effect that contain only the spherical approximation. It is shown that neither method relies on a linear relationship between gravity anomalies and height (as claimed by some). Numerical tests, however, show that the practical implementations of these two approaches yield significant differences. Computational tests were performed in three areas of the USA, using 1×1 grids of gravity data and 30×30 grids of height data to compute the gravimetric geoid undulation, and GPS/leveled heights to compute the geometric geoid undulation. Using the latter as a control, analyses of the gravimetric undulations indicate that while in areas with smooth terrain no substantial differences occur between the gravity reduction methods, the Moritz–Pellinen (MP) approach is clearly superior to the Vanicek–Martinec (VM) approach in areas of rugged terrain. In theory, downward continuation is a significant aspect of either approach. Numerically, however, based on the test data, neither approach benefited by including this effect in the areas having smooth terrain. On the other hand, in the rugged, mountainous area, the gravimetric geoid based on the VM approach was improved slightly, but with the MP approach it suffered significantly. The latter is attributed to an inability to model the downward continuation of the Bouguer anomaly accurately in rugged terrain. Applying the higher-order, more accurate gravity reduction formulas, instead of their corresponding planar and linear approximations, yielded no improvement in the accuracy of the gravimetric geoid undulation based on the available data.  相似文献   

7.
The method of analytical downward continuation has been used for solving Molodensky’s problem. This method can also be used to reduce the surface free air anomaly to the ellipsoid for the determination of the coefficients of the spherical harmonic expansion of the geopotential. In the reduction of airborne or satellite gradiometry data, if the sea level is chosen as reference surface, we will encounter the problem of the analytical downward continuation of the disturbing potential into the earth, too. The goal of this paper is to find out the topographic effect of solving Stoke’sboundary value problem (determination of the geoid) by using the method of analytical downward continuation. It is shown that the disturbing potential obtained by using the analytical downward continuation is different from the true disturbing potential on the sea level mostly by a −2πGρh 2/p. This correction is important and it is very easy to compute and add to the final results. A terrain effect (effect of the topography from the Bouguer plate) is found to be much smaller than the correction of the Bouguer plate and can be neglected in most cases. It is also shown that the geoid determined by using the Helmert’s second condensation (including the indirect effect) and using the analytical downward continuation procedure (including the topographic effect) are identical. They are different procedures and may be used in different environments, e.g., the analytical downward continuation procedure is also more convenient for processing the aerial gravity gradient data. A numerical test was completed in a rough mountain area, 35°<ϕ<38°, 240°<λ<243°. A digital height model in 30″×30″ point value was used. The test indicated that the terrain effect in the test area has theRMS value ±0.2−0.3 cm for geoid. The topographic effect on the deflections of the vertical is around1 arc second.  相似文献   

8.
Assuming that the gravity anomaly and disturbing potential are given on a reference ellipsoid, the result of Sjöberg (1988, Bull Geod 62:93–101) is applied to derive the potential coefficients on the bounding sphere of the ellipsoid to order e 2 (i.e. the square of the eccentricity of the ellipsoid). By adding the potential coefficients and continuing the potential downward to the reference ellipsoid, the spherical Stokes formula and its ellipsoidal correction are obtained. The correction is presented in terms of an integral over the unit sphere with the spherical approximation of geoidal height as the argument and only three well-known kernel functions, namely those of Stokes, Vening-Meinesz and the inverse Stokes, lending the correction to practical computations. Finally, the ellipsoidal correction is presented also in terms of spherical harmonic functions. The frequently applied and sometimes questioned approximation of the constant m, a convenient abbreviation in normal gravity field representations, by e 2/2, as introduced by Moritz, is also discussed. It is concluded that this approximation does not significantly affect the ellipsoidal corrections to potential coefficients and Stokes formula. However, whether this standard approach to correct the gravity anomaly agrees with the pure ellipsoidal solution to Stokes formula is still an open question.  相似文献   

9.
Geoid determination using one-step integration   总被引：1，自引：1，他引：0  
P. Novák 《Journal of Geodesy》2003,77(3-4):193-206
A residual (high-frequency) gravimetric geoid is usually computed from geographically limited ground, sea and/or airborne gravimetric data. The mathematical model for its determination from ground gravity is based on the transformation of observed discrete values of gravity into gravity potential related to either the international ellipsoid or the geoid. The two reference surfaces are used depending on height information that accompanies ground gravity data: traditionally orthometric heights determined by geodetic levelling were used while GPS positioning nowadays allows for estimation of geodetic (ellipsoidal) heights. This transformation is usually performed in two steps: (1) observed values of gravity are downward continued to the ellipsoid or the geoid, and (2) gravity at the ellipsoid or the geoid is transformed into the corresponding potential. Each of these two steps represents the solution of one geodetic boundary-value problem of potential theory, namely the first and second or third problem. Thus two different geodetic boundary-value problems must be formulated and solved, which requires numerical evaluation of two surface integrals. In this contribution, a mathematical model in the form of a single Fredholm integral equation of the first kind is presented and numerically investigated. This model combines the solution of the first and second/third boundary-value problems and transforms ground gravity disturbances or anomalies into the harmonically downward continued disturbing potential at the ellipsoid or the geoid directly. Numerical tests show that the new approach offers an efficient and stable solution for the determination of the residual geoid from ground gravity data.  相似文献   

10.
 The downward continuation of the harmonic disturbing gravity potential, derived at flight level from discrete observations of airborne gravity by the spherical Hotine integral, to the geoid is discussed. The initial-boundary-value approach, based on both the direct and inverse solution to Dirichlet's problem of potential theory, is used. Evaluation of the discretized Fredholm integral equation of the first kind and its inverse is numerically tested using synthetic airborne gravity data. Characteristics of the synthetic gravity data correspond to typical airborne data used for geoid determination today and in the foreseeable future: discrete gravity observations at a mean flight height of 2 to 6 km above mean sea level with minimum spatial resolution of 2.5 arcmin and a noise level of 1.5 mGal. Numerical results for both approaches are presented and discussed. The direct approach can successfully be used for the downward continuation of airborne potential without any numerical instabilities associated with the inverse approach. In addition to these two-step approaches, a one-step procedure is also discussed. This procedure is based on a direct relationship between gravity disturbances at flight level and the disturbing gravity potential at sea level. This procedure provided the best results in terms of accuracy, stability and numerical efficiency. As a general result, numerically stable downward continuation of airborne gravity data can be seen as another advantage of airborne gravimetry in the field of geoid determination. Received: 6 June 2001 / Accepted: 3 January 2002  相似文献   

11.
Least-squares collocation and Stokes integral formula, as implemented using the Fast Fourier Technique, handle the harmonic downward continuation problem quite differently. FFT furthermore requires gridded data, amplifying the difference of methods.We have in this paper studied numerically the effects of downward continuation and gridding in a mountainous area in central Norway. Topographically smoothed data were used in order to reduce these effects. Despite the smoothing, it was found that the vertical gravity gradient had values up to -11 mgal/km. The corresponding differences between geoid heights and the height anomalies at altitude reached 12 cm.The differences between geoid heights obtained using collocation or FFT with gravity data at terrain level or sea level showed differences between the values of up to 10 cm r.m.s. A part of this difference was a consequence of different data areas used in the FFT and collocation solution, though.Major discrepancies between the solutions were found in areas where the topographic smoothing could not be applied (deep fjords with no depth information in the used DTM) or where there seemed to be gross errors in the data.We conclude that proper handling of harmonic continuation is important, even when we as here have used a 1 km resolution DTM for the calculation of topographic effects. The effect of data gridding, required for the FFT method, seems not to be as serious as the need to limit the data distribution area, required when least squares collocation is used with randomly distributed data.  相似文献   

12.
Ellipsoidal geoid computation   总被引：1，自引：1，他引：0  
Modern geoid computation uses a global gravity model, such as EGM96, as a third component in a remove–restore process. The classical approach uses only two: the reference ellipsoid and a geometrical model representing the topography. The rationale for all three components is reviewed, drawing attention to the much smaller precision now needed when transforming residual gravity anomalies. It is shown that all ellipsoidal effects needed for geoid computation with millimetric accuracy are automatically included provided that the free air anomaly and geoid are calculated correctly from the global model. Both must be consistent with an ellipsoidal Earth and with the treatment of observed gravity data. Further ellipsoidal corrections are then negligible. Precise formulae are developed for the geoid height and the free air anomaly using a global gravity model, given as spherical harmonic coefficients. Although only linear in the anomalous potential, these formulae are otherwise exact for an ellipsoidal reference Earth—they involve closed analytical functions of the eccentricity (and the Earths spin rate), rather than a truncated power series in e2. They are evaluated using EGM96 and give ellipsoidal corrections to the conventional free air anomaly ranging from –0.84 to +1.14 mGal, both extremes occurring in Tibet. The geoid error corresponding to these differences is dominated by longer wavelengths, so extrema occur elsewhere, rising to +766 mm south of India and falling to –594 mm over New Guinea. At short wavelengths, the difference between ellipsoidal corrections based only on EGM96 and those derived from detailed local gravity data for the North Sea geoid GEONZ97 has a standard deviation of only 3.3 mm. However, the long-wavelength components missed by the local computation reach 300 mm and have a significant slope. In Australia, for example, such a slope would amount to a 600-mm rise from Perth to Cairns.  相似文献   

13.
We investigate the stability of a discrete downward continuation problem for geoid determination when the surface gravity observations are harmonically continued from the Earth's surface to the geoid. The discrete form of Poisson's integral is used to set up the system of linear algebraic equations describing the problem. The posedness of the downward continuation problem is then expressed by means of the conditionality of the matrix of a system of linear equations. The eigenvalue analysis of this matrix for a particularly rugged region of the Canadian Rocky Mountains shows that the discrete downward continuation problem is stable once the topographical heights are discretized with a grid step of size 5 arcmin or larger. We derive two simplified criteria for analysing the conditionality of the discrete downward continuation problem. A comparison with the proper eigenvalue analysis shows that these criteria provide a fairly reliable view into the conditionality of the problem.The compensation of topographical masses is a possible way how to stabilize the problem as the spectral contents of the gravity anomalies of compensated topographical masses may significantly differ from those of the original free-air gravity anomalies. Using surface gravity data from the Canadian Rocky Mountains, we investigate the efficiency of highly idealized compensation models, namely the Airy-Heiskanen model, the Pratt-Hayford model, and Helmert's 2nd condensation technique, to dampen high-frequency oscillations of the free-air gravity anomalies. We show that the Airy-Heiskanen model reduces high-frequencies of the data in the most efficient way, whereas Helmert's 2nd condensation technique in the least efficient way. We have found areas where a high-frequency part of the surface gravity data has been completely removed by adopting the Airy-Heiskanen model which is in contrast to the nearly negligible dampening effect of Helmert's 2nd condensation technique. Hence, for computation of the geoid over the Canadian Rocky Mountains, we recommend the use of the Airy-Heiskanen compensation model to reduce the gravitational effect of topographical masses.In addition, we propose to solve the discrete downward continuation problem by means of a simple Jacobi's iterative scheme which finds the solution without determining and storing the matrix of a system of equations. By computing the spectral norm of the matrix of a system of equations for the topographical 5 × 5 heights from a region of the Canadian Rocky Mountains, we rigorously show that Jacobi's iterations converge to the solution; that the problem was well posed then ensures that the solution is not contaminated by large roundoff errors. On the other hand, we demonstrate that for a rugged mountainous region of the Rocky Mountains the discrete downward continuation problem becomes ill-conditioned once the grid step size of both the surface observations and the solution is smaller than 1 arcmin. In this case, Jacobi's iterations converge very slowly which prevents their use for searching the solution due to accumulating roundoff errors.  相似文献   

14.
Local geoid determination from airborne vector gravimetry   总被引：3，自引：2，他引：1  
Methods are illustrated to compute the local geoid using the vertical and horizontal components of the gravity disturbance vector derived from an airborne GPS/inertial navigation system. The data were collected by the University of Calgary in a test area of the Canadian Rocky Mountains and consist of multiple parallel tracks and two crossing tracks of accelerometer and gyro measurements, as well as precise GPS positions. Both the boundary-value problem approach (Hotines integral) and the profiling approach (line integral) were applied to compute the disturbing potential at flight altitude. Cross-over adjustments with minimal control were investigated and utilized to remove error biases and trends in the estimated gravity disturbance components. Final estimation of the geoid from the vertical gravity disturbance included downward continuation of the disturbing potential with correction for intervening terrain masses. A comparison of geoid estimates to the Canadian Geoid 2000 (CGG2000) yielded an average standard deviation per track of 14 cm if they were derived from the vertical gravity disturbance (minimally controlled with a cross-over adjustment), and 10 cm if derived from the horizontal components (minimally controlled in part with a simulated cross-over adjustment). Downward continuation improved the estimates slightly by decreasing the average standard deviation by about 0.5 cm. The application of a wave correlation filter to both types of geoid estimates yielded significant improvement by decreasing the average standard deviation per track to 7.6 cm.  相似文献   

15.
The application of Stokes’s formula to determine the geoid height requires that topographic and atmospheric masses be mathematically removed prior to Stokes integration. This corresponds to the applications of the direct topographic and atmospheric effects. For a proper geoid determination, the external masses must then be restored, yielding the indirect effects. Assuming an ellipsoidal layering of the atmosphere with 15% increase in its density towards the poles, the direct atmospheric effect on the geoid height is estimated to be −5.51 m plus a second-degree zonal harmonic term with an amplitude of 1.1 cm. The indirect effect is +5.50 m and the total geoid correction thus varies between −1.2 cm at the equator to 1.9 cm at the poles. Finally, the correction needed to the atmospheric effect if Stokes’s formula is used in a spherical approximation, rather than an ellipsoidal approximation, of the Earth varies between 0.3 cm and 4.0 cm at the equator and pole, respectively.  相似文献   

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

17.
Prior to Stokes integration, the gravitational effect of atmospheric masses must be removed from the gravity anomaly g. One theory for the atmospheric gravity effect on the geoid is the well-known International Association of Geodesy approach in connection with Stokes integral formula. Another strategy is the use of a spherical harmonic representation of the topography, i.e. the use of a global topography computed from a set of spherical harmonics. The latter strategy is improved to account for local information. A new formula is derived by combining the local contribution of the atmospheric effect computed from a detailed digital terrain model and the global contribution computed from a spherical harmonic model of the topography. The new formula is tested over Iran and the results are compared with corresponding results from the old formula which only uses the global information. The results show significant differences. The differences between the two formulas reach 17 cm in a test area in Iran.  相似文献   

18.
Geoid and quasigeoid modelling from gravity anomalies by the method of least squares modification of Stokes’s formula with additive corrections is adapted for the usage with gravity disturbances and Hotine’s formula. The biased, unbiased and optimum versions of least squares modification are considered. Equations are presented for the four additive corrections that account for the combined (direct plus indirect) effect of downward continuation (DWC), topographic, atmospheric and ellipsoidal corrections in geoid or quasigeoid modelling. The geoid or quasigeoid modelling scheme by the least squares modified Hotine formula is numerically verified, analysed and compared to the Stokes counterpart in a heterogeneous study area. The resulting geoid models and the additive corrections computed both for use with Stokes’s or Hotine’s formula differ most in high topography areas. Over the study area (reaching almost 2 km in altitude), the approximate geoid models (before the additive corrections) differ by 7 mm on average with a 3 mm standard deviation (SD) and a maximum of 1.3 cm. The additive corrections, out of which only the DWC correction has a numerically significant difference, improve the agreement between respective geoid or quasigeoid models to an average difference of 5 mm with a 1 mm SD and a maximum of 8 mm.  相似文献   

19.
The solutions of four ellipsoidal approximations for the gravimetric geoid are reviewed: those of Molodenskii et al., Moritz, Martinec and Grafarend, and Fei and Sideris. The numerical results from synthetic tests indicate that Martinec and Grafarends solution is the most accurate, while the other three solutions contain an approximation error which is characterized by the first-degree surface spherical harmonic. Furthermore, the first 20 degrees of the geopotential harmonic series contribute approximately 90% of the ellipsoidal correction. The determination of a geoid model from the generalized Stokes scheme can accurately account for the ellipsoidal effect to overcome the first-degree surface spherical harmonic error regardless of the solution used.  相似文献   

20.
利用空中平均重力异常确定区域大地水准面   总被引：3，自引：0，他引：3  
提出了直接利用空中平均重力异常计算区域大地水准面的方法。模拟计算的结果表明, 该方法与传统的利用地面平均空间重力异常确定的大地水准面精度相当, 但其显著优点是勿需空中重力异常的向下解析延拓, 从而可以避免延拓误差对大地水准面精化的影响。  相似文献   

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

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