首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
The rigorous determination of orthometric heights   总被引：1，自引：2，他引：1  
The main problem of the rigorous definition of the orthometric height is the evaluation of the mean value of the Earth’s gravity acceleration along the plumbline within the topography. To find the exact relation between rigorous orthometric and Molodensky’s normal heights, the mean gravity is decomposed into: the mean normal gravity, the mean values of gravity generated by topographical and atmospheric masses, and the mean gravity disturbance generated by the masses contained within geoid. The mean normal gravity is evaluated according to Somigliana–Pizzetti’s theory of the normal gravity field generated by the ellipsoid of revolution. Using the Bruns formula, the mean values of gravity along the plumbline generated by topographical and atmospheric masses can be computed as the integral mean between the Earth’s surface and geoid. Since the disturbing gravity potential generated by masses inside the geoid is harmonic above the geoid, the mean value of the gravity disturbance generated by the geoid is defined by applying the Poisson integral equation to the integral mean. Numerical results for a test area in the Canadian Rocky Mountains show that the difference between the rigorously defined orthometric height and the Molodensky normal height reaches ∼0.5 m.  相似文献   

2.
蒋涛  党亚民  郭春喜  陈斌  章传银 《测绘学报》2022,51(8):1757-1767
2020珠峰高程测量,首次确定并发布了基于国际高程参考系统(IHRS)的珠峰正高。在珠峰地区实现国际高程参考系统,采用的方案是建立珠峰区域高精度重力大地水准面。利用地球重力场谱组合理论和基于数据驱动的谱权确定方法,测试优选参考重力场模型及其截断阶数和球冠积分半径等关键参数,联合航空和地面重力等数据建立了珠峰区域重力似大地水准面模型,61点高精度GNSS水准高程异常检核表明,模型精度达3.8 cm,加入航空重力数据后模型精度提升幅度达51.3%。提出顾及高差改正的峰顶高程异常内插方法,采用顾及地形质量影响的高程异常——大地水准面差距转换改正严密公式,使用峰顶实测地面重力数据,基于国际高程参考系统定义的重力位值W0和GRS80参考椭球,最终确定了国际高程参考系统中的高精度珠峰峰顶大地水准面差距。  相似文献   

3.
申文斌 《测绘学报》2012,41(5):670-675
确定全球大地水准面最常用的方法是斯托克司方法。然而,除了人们熟知的缺陷之外,斯托克司方法还存在人们没有意识到的理论困难：当大地水准面位于参考椭球（WGS84椭球）内部时,在大地水准面上及其与参考椭球面界定的区域中扰动位没有定义,当然在这部分区域也不调和。为了解决这一困难,可以选取一个包含在大地水准面内部的由四个基本参数唯一确定其外部正常重力位的参考椭球（简称内部椭球）,其中心与 WGS84 椭球的中心重合,其中的两个基本参数,旋转角速度和地心引力常数,与 WGS84 椭球面的相同,另外两个参数,半长轴和扁率,如此选取,使得内部椭球产生的新的正常重力位在 WGS84 椭球面上与大地水准面上的重力位 相等。这样,传统的斯托克司方法中存在的理论困难不复存在。  相似文献   

4.
The formulas for the determination of the coefficients of the spherical harmonic expansion of the disturbing potential of the earth are defined for data given on a sphere. In order to determine the spherical harmonic coefficients, the gravity anomalies have to be analytically downward continued from the earth's surface to a sphere—at least to the ellipsoid. The goal of this paper is to continue the gravity anomalies from the earth's surface downward to the ellipsoid using recent elevation models. The basic method for the downward continuation is the gradient solution (theg 1 term). The terrain correction has also been computed because of the role it can play as a correction term when calculating harmonic coefficients from surface gravity data. Theg 1 term and the terrain correction were expanded into the spherical harmonics up to180 th order. The corrections (theg 1 term and the terrain correction) have the order of about 2% of theRMS value of degree variance of the disturbing potential per degree. The influences of theg 1 term and the terrain correction on the geoid take the order of 1 meter (RMS value of corrections of the geoid undulation) and on the deflections of the vertical is of the order 0.1″ (RMS value of correction of the deflections of the vertical).  相似文献   

5.
The Meissl scheme for the geodetic ellipsoid   总被引：2，自引：1，他引：1  
We present a variant of the Meissl scheme to relate surface spherical harmonic coefficients of the disturbing potential of the Earth’s gravity field on the surface of the geodetic ellipsoid to surface spherical harmonic coefficients of its first- and second-order normal derivatives on the same or any other ellipsoid. It extends the original (spherical) Meissl scheme, which only holds for harmonic coefficients computed from geodetic data on a sphere. In our scheme, a vector of solid spherical harmonic coefficients of one quantity is transformed into spherical harmonic coefficients of another quantity by pre-multiplication with a transformation matrix. This matrix is diagonal for transformations between spheres, but block-diagonal for transformations involving the ellipsoid. The computation of the transformation matrix involves an inversion if the original coefficients are defined on the ellipsoid. This inversion can be performed accurately and efficiently (i.e., without regularisation) for transformation among different gravity field quantities on the same ellipsoid, due to diagonal dominance of the matrices. However, transformations from the ellipsoid to another surface can only be performed accurately and efficiently for coefficients up to degree and order 520 due to numerical instabilities in the inversion.  相似文献   

6.
Using the spherical harmonic representations of the earth's disturbing potential and its functionals, we derive the inverse Vening Meinesz formula, which converts deflection of the vertical to gravity anomaly using the gradient of the H function. The deflection-geoid formula is also derived that converts deflection to geoidal undulation using the gradient of the C function. The two formulae are implemented by the 1D FFT and the 2D FFT methods. The innermost zone effect is derived. The inverse Vening Meinesz formula is employed to compute gravity anomalies and geoidal undulations over the South China Sea using deflections from Seasat, Geosat, ERS-1 and TOPEX//POSEIDON satellite altimetry. The 1D FFT yields the best result of 9.9-mgal rms difference with the shipborne gravity anomalies. Using the simulated deflections from EGM96, the deflection-geoid formula yields a 4-cm rms difference with the EGM96-generated geoid. The predicted gravity anomalies and geoidal undulations can be used to study the tectonic structure and the ocean circulations of the South China Sea. Received: 7 April 1997 / Accepted: 7 January 1998  相似文献   

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

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.
本文给出了Hotine 函数法的椭球面积分解,以应用于计算精确的大地水准面起伏。计算表明,当积分半径为20°时,我国近海的椭球改正只有10cm,远比stokes公式的椭球改正要小。  相似文献   

10.
程芦颖 《测绘学报》2013,42(2):203-210
基于物理大地测量边值问题的解,利用一阶边界算子定义,推导重力异常Δg、单层密度μ、大地水准面高N,垂线偏差ε、扰动重力δg等扰动场元的解。利用球谐函数的正交特性,通过对核函数的算子运算,可以得到上述扰动场元的有关逆变换公式。相对经典物理大地测量公式应用的边界面条件,笔者将含有因子r的对应扰动场元反演关系的公式称为广义积分公式。针对常用的重力异常Δg、大地水准面高N,垂线偏差ε、扰动重力δg计算,重点分析它们之间的变换关系,给出利用某个选定扰动场元计算其他扰动场元的广义积分公式。同时,通过对积分边界面的讨论,分析经典公式与广义积分公式的差异和联系。最后,给出所有外部扰动场元与核函数映射的关系表。  相似文献   

11.
J. Li 《Journal of Geodesy》2002,76(4):226-231
 A formula for computing the gravity disturbance and gravity anomaly from the second radial derivative of the disturbing potential is derived in detail using the basic differential equation with spherical approximation in physical geodesy and the modified Poisson integral formula. The derived integral in the space domain, expressed by a spherical geometric quantity, is then converted to a convolution form in the local planar rectangular coordinate system tangent to the geoid at the computing point, and the corresponding spectral formulae of 1-D FFT and 2-D FFT are presented for numerical computation. Received: 27 December 2000 / Accepted: 3 September 2001  相似文献   

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

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

14.
The vertical component obtained from the Global Positioning System (GPS) observations is from the ellipsoid (a mathematical surface), and therefore needs to be converted to the orthometric height, which is from the geoid (represented by the mean sea level). The common practice is to use existing bench marks (around the four corners of a project area and interpolate for the rest of the area), but in many areas bench marks may not be available, in which case an existing geoid undulation is used. Present available global geoid undulation values are not generally as detailed as needed, and in many areas they are not known better than ±1 to ±5 m, because of many limitations. This article explains the difficulties encountered in obtaining precise geoid undulation with some example computations, and proposes a technique of applying corrections to the best available global geoid undulations using detailed free-air gravity anomalies (within a 2° × 2° area) to get relative centimeter accuracy. Several test computations have been performed to decide the optimal block sizes and the effective spherical distances to compute the regional and the local effects of gravity anomalies on geoid undulations by using the Stokes integral. In one test computation a 2° × 2° area was subdivided into smaller surface elements. A difference of 37.34 ± 1.6 cm in geoid undulation was obtained over the same 2° × 2° area when 1° × 1° block sizes were replaced by a combination of 5' × 5' and 1' × 1' subdivision integration elements (block sizes).  相似文献   

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

16.
利用最新的全球引力位模型-EGM2008对经典的重力异常阶方差模型进行了分析比较,分析表明,经典的阶方差模型由于限于当时的观测条件,已经不能准确地描述扰动场元在各个频段的频谱分布。在Moritz阶方差模型基础上,利用EGM2008位模型获得的2160阶阶方差重新构建了新的分段重力异常阶方差模型-TSD模型,该模型与EGM2008位模型计算的阶方差比较其标准差和均值分别为0.25mgal2 、0.0 。利用TSD模型计算了不同频段内大地水准面高、重力异常、扰动重力、垂线偏差四个重力场扰动场元的频谱特征,计算结果表明：扰动场元频谱分布较之传统分析结果有较大的变化,其中重力异常、扰动重力及垂线偏差在中、低频部分的能量有明显的增加而高频及甚高频部分的比重有明显的减少。  相似文献   

17.
从经典边值问题理论及球谐函数理论出发,在空域推导获得了由大地水准面高以及垂线偏差计算扰动重力的解析计算公式,为利用卫星测高数据反演海洋扰动重力提供了理论基础。针对全球海洋区域和局部海洋区域的扰动重力反演,在前人已有工作基础上,提出了改进的基于一维FFT的精确快速算法,保证了计算结果与原解析方法完全一致,且计算速度提高约20倍。该算法在提高计算效率的同时避免了由于引入FFT而产生的混叠、边缘效应问题,而且对观测数据的序列长度没有硬性要求,使得应用更加灵活。利用EGM2008地球重力场模型分别生成了2.5'分辨率大地水准面高数据和垂线偏差数据,按照本文提出的改进方法(采用全球积分计算)分别反演获得了全球及局部海洋区域的扰动重力。经比较分析,由大地水准面和垂线偏差分别反演获得的扰动重力其差异在0.8×10-5 m/s2以内,这说明两种反演方法是基本一致的,但在数据包含系统误差的情况下,由垂线偏差反演扰动重力具有一定优势。  相似文献   

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

19.
This paper is devoted to the spherical and spheroidal harmonic expansion of the gravitational potential of the topographic masses in the most rigorous way. Such an expansion can be used to compute gravimetric topographic effects for geodetic and geophysical applications. It can also be used to augment a global gravity model to a much higher resolution of the gravitational potential of the topography. A formulation for a spherical harmonic expansion is developed without the spherical approximation. Then, formulas for the spheroidal harmonic expansion are derived. For the latter, Legendre’s functions of the first and second kinds with imaginary variable are expanded in Laurent series. They are then scaled into two real power series of the second eccentricity of the reference ellipsoid. Using these series, formulas for computing the spheroidal harmonic coefficients are reduced to surface harmonic analysis. Two numerical examples are presented. The first is a spherical harmonic expansion to degree and order 2700 by taking advantage of existing software. It demonstrates that rigorous spherical harmonic expansion is possible, but the computed potential on the geoid shows noticeable error pattern at Polar Regions due to the downward continuation from the bounding sphere to the geoid. The second numerical example is the spheroidal expansion to degree and order 180 for the exterior space. The power series of the second eccentricity of the reference ellipsoid is truncated at the eighth order leading to omission errors of 25 nm (RMS) for land areas, with extreme values around 0.5 mm to geoid height. The results show that the ellipsoidal correction is 1.65 m (RMS) over land areas, with maximum value of 13.19 m in the Andes. It shows also that the correction resembles the topography closely, implying that the ellipsoidal correction is rich in all frequencies of the gravity field and not only long wavelength as it is commonly assumed.  相似文献   

20.
Spherical cap harmonic analysis is the appropriate analytical technique for modelling Laplacian potential and the corresponding field components over a spherical cap. This paper describes the use of this method by means of a least-squares approach for local gravity field representation. Formulations for the geoid undulation and the components ξ, η of the deflection of the vertical are derived, together with some warnings in the application of the technique. Although most of the formulations have been given by another paper, these were confusing or even incorrect, mainly because of an improper application of the spherical cap harmonic analysis. Received: 16 January 1996 / Accepted: 17 March 1997  相似文献   

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

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