A solution to the downward continuation effect on the geoid determined by Stokes' formula

 The analytical continuation of the surface gravity anomaly to sea level is a necessary correction in the application of Stokes' formula for geoid estimation. This process is frequently performed by the inversion of Poisson's integral formula for a sphere. Unfortunately, this integral equation corresponds to an improperly posed problem, and the solution is both numerically unstable, unless it is well smoothed, and tedious to compute. A solution that avoids the intermediate step of downward continuation of the gravity anomaly is presented. Instead the effect on the geoid as provided by Stokes' formula is studied directly. The practical solution is partly presented in terms of a truncated Taylor series and partly as a truncated series of spherical harmonics. Some simple numerical estimates show that the solution mostly meets the requests of a 1-cm geoid model, but the truncation error of the far zone must be studied more precisely for high altitudes of the computation point. In addition, it should be emphasized that the derived solution is more computer efficient than the detour by Poisson's integral. Received: 6 February 2002 / Accepted: 18 November 2002 Acknowledgements. Jonas ?gren carried out the numerical calculations and gave some critical and constructive remarks on a draft version of the paper. This support is cordially acknowledged. Also, the thorough work performed by one unknown reviewer is very much appreciated.     

非线性二维转换模型的稳健总体最小二乘算法

针对应用线性最小二乘估计准则求解非线性平面转换模型参数时,通过定义间接参数将模型线性化的方法不能直接求解转换模型参数的问题,该文在非线性平面转换模型的基础上,建立线性模型,实现平面坐标的转换。为解决控制点已知坐标与观测坐标中均含有误差对转换参数求解的影响,对应用稳健总体最小二乘求解线性模型参数的算法进行讨论。最后,通过算例比较稳健总体最小二乘算法与最小二乘算法在抗差性方面的优势。结果表明,稳健总体最小二乘算法更适用于应用线性模型求解未知控制点的转换坐标。     

Truncated geoid and gravity inversion for one point-mass anomaly

The truncated geoid, defined by the truncated Stokes' integral transform, an integral convolution of gravity anomalies with the Stokes' function on a spherical cap, is often used as a mathematical tool in geoid computations via Stokes' integral to overcome computational difficulties, particularly the need to integrate over the entire boundary spheroid. The objective of this paper is to demonstrate that the truncated geoid does, besides having mathematical applications, have physical interpretation, and thus may be used in gravity inversion. A very simple model of one point-mass anomaly is chosen and a method for inverting its synthetic gravity field with the use of the truncated geoid is presented. The method of inverting the synthetic field generated by one point-mass anomaly has become fundamental for the authors' inversion studies for sets of point-mass anomalies, which are published in a separate paper. More general applications are currently under investigation. Since an inversion technique for physically meaningful mass distributions based on the truncated geoid has not yet been developed, this work is not related to any of the existing gravity inversion techniques. The inversion for one point mass is based on the onset of the so-called dimple event, which occurs in the sequence of surfaces (or profiles) of the first derivative of the truncated geoid with respect to the truncation parameter (radius of the integration cap), its only free parameter. Computing the truncated geoid at various values of the truncation parameter may be understood as spatial filtering of surface gravity data, a type of weighted spherical windowing method. Studying the change of the truncated geoid represented by its first derivative may be understood as a data enhancement method. The instant of the dimple onset is practically independent of the mass of the point anomaly and linearly dependent on its depth. Received: 26 September 1996 /Accepted: 28 September 1998     

The effect of topography on the determination of the geoid using analytical downward continuation

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.     

基于模糊综合评判的土木工程测量实习成绩评定方法研究

本文针对目前工程测量实习课成绩评定中存在的问题,以模糊数学理论为依据,提出一种科学合理的评价标准,它利用模糊综合评判方法结合实际实习的指标体系,计算学生的实习成绩。该方法也对其他学科的实习成绩的评定具有一定的参考价值。     

An Adaptive Water Extraction Method from Remote Sensing Image Based on NDWI

Water is one of the most common and important objects on the earth, and its extraction is of great significance to many related researches in remote sensing domain. However, water always appears diversely, which makes its extraction not so simple. Many former methods are developed to extract water, which mainly based on a single model and only use spectral information, but the results are not so satisfying. An adaptive extraction method based on normalized difference water index (NDWI) is proposed here to extract water completely and accurately from remote sensing image. This study first compute NDWI to enhance water??s spectral information, and then it is redefined so as to use the modified histogram auto-segmentation method to initially separate water from background; next, after segmentation, water pixels can be searched out and are taken as seed points to proceed region growing to get the local area of water; last, the edge of the local area is searched by a window template, and iterative classification within it is employed to precisely extract water??s precise partition. Experiments are carried out here on an ETM+ image of a paralic area to extract water. Through comparison with other commonly used methods, it shows that the performance of the proposed method is superior to the others.     

Methods for computing deflections of the vertical by modifying Vening-Meinesz' function

Summary The application of combined data (satellite and terrestrial data) to the practical computation of height anomalies or the deflections of the vertical was originally suggested by (Molodensky et al. 1962). This idea usually leads to the modification of Stokes' or Vening-Meinesz' functions in the integration procedure. In the recent decade there were various suggestions in this regard especially for the computation of height anomalies. For example, a considerable mathematical insight into the modification of Stokes' function and the truncation of its integral has been provided by (Meissl 1971, Houtze et al. 1979, Rapp 1980, Jekeli 1980). Five different methods for computing deflections of the vertical by modifying Vening-Meinesz' function are studied and compared with each other. The corresponding formulae, the values of the coefficients in each method and the estimations of their corresponding potential coefficient error and truncation error are given in this article. This paper was written at the Institut f. Angewandte Geod?sie, Technische Universit?t Graz, Austria.     

A method of evaluating the truncation error coefficients for geoidal height

Neglecting distant zones in the computation of geoidal height using Stokes' formula gives rise to some truncation error. This truncation error is expressible as a weighted summation of the zonal harmonic components of the gravity anomaly. Making use of the well-known properties of Legendre polynomials, a compact method of computing these theoretical coefficients has been developed in this paper.     

APPROXIMATE METHODS OF SOLVING NORMAL EQUATIONS

Abstract

My attention has recently been drawn to an article by G. H. Menzies (E.5.R., vi, 46, 474) on this subject. The present note is intended to point out improvements in the Anér method, which he favours, and to refute some of his criticisms of the Relaxation method. References are to page and table numbers in Menzies's paper.     

Some modifications of Stokes' formula that account for truncation and potential coefficient errors

 Stokes' formula from 1849 is still the basis for the gravimetric determination of the geoid. The modification of the formula, originating with Molodensky, aims at reducing the truncation error outside a spherical cap of integration. This goal is still prevalent among various modifications. In contrast to these approaches, some least-squares types of modification that aim at reducing the truncation error, as well as the error stemming from the potential coefficients, are demonstrated. The least-squares estimators are provided in the two cases that (1) Stokes' kernel is a priori modified (e.g. according to Molodensky's approach) and (2) Stokes' kernel is optimally modified to minimize the global mean square error. Meissl-type modifications are also studied. In addition, the use of a higher than second-degree reference field versus the original (Pizzetti-type) reference field is discussed, and it is concluded that the former choice of reference field implies increased computer labour to achieve the same result as with the original reference field. Received: 14 December 1998 / Accepted: 4 October 1999