从现代测绘技术发展谈测绘继续教育

20世纪90年代中后期,以“3S”技术和“4D”产品的普及推广为标志,我国的传统测绘产业开始了向现代空间地理信息产业的过渡,本文从时空信息获取、加工、管理和服务4个方面对地球空间信息未来的技术发展作了简要的叙述,这对测绘行业本身、对其他行业、对整个社会发展都具有重大意义。继续教育是社会发展加快、竞争激烈的现代性产物。在高科技环境下,测绘继续教育对提高在职测绘专业技术人员的素质,适应国内外测绘行业日益激烈的竞争起着关键性的作用。     

球内Dirichlet问题解及其应用

本文基于球内调和函数的Ｄｉｒｉｃｈｌｅｔ问题的球谐解式,推导了球内调和空间的Ｐｏｉｓｓｏｎ积分,将其应用于航空重力测量数据的向下延拓时,积分边界面是空中面,边界值是空中重力异常或纯重力异常,推求地面重力异常可直接积分计算,而勿需像球外Ｐｏｉｓｓｏｎ积分那样迭代求解积分方程。     

基于改进泰勒级数法的总强度磁异常稳健向下延拓

位场向下延拓隶属于经典不适定问题,观测数据中细微的误差在向下延拓过程中都被严重放大,甚至会掩盖真实信息。如何精确地求解总强度磁异常（Bm）在垂直方向的各阶导数,是利用泰勒级数实现稳健向下延拓的关键。为此,本文首先分析了调和函数的相关性质,从理论上证明了Bm为准调和函数的结论,在精确计算各阶垂向导数基础上,提出利用改进泰勒级数实现磁场稳健向下延拓。为降低边界效应对向下延拓计算结果的影响,提出采用半余弦函数对磁场在4个方向上进行平滑扩边处理。通过球体与长方体仿真试验以及航空、船载实测磁场数据对提出方法进行了验证。结论表明,提出的技术方法可实现磁场稳健向下延拓,当观测数据无噪声时,计算结果精度要明显优于现行的FFT法、常规泰勒级数法以及积分-迭代法;当观测数据含有噪声时,本文方法和积分-迭代法计算结果精度相当。     

重力数据的多尺度边缘分析及其应用

In order to improve the processing and interpretation of gravity data, multiscale edge theory in image processing is introduced into the study of gravity field. Multiscale edges of gravity anomaly are analyzed based on a special wavelet. It shows that the multiscale edges are the extrema points of the horizontal gravity gradient at different heights, which are related to the sharp discontinuities of underground sources. The applications of multiscale edge in downward continuation and gravity inversion are discussed. The simulated examples show that the multiscale edges can be applied to stabilize the downward continuation operator when the continuation height is low. The multiscale edges also have a convenient application to infer the geometry of the underground source. Moreover, the gravity inversion algorithm based on the multiscale edges has a good antinoise property. Supported by the National Natural Science Foundation of China(No.40704003), the National 973 Program of China(No.2007CB714405), the Open Research Fund from Key Laboratory of Geospace Environment and Geodesy(No.04-01-08).     

Precise geoid determination over Sweden using the Stokes–Helmert method and improved topographic corrections

 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     

卫星重力梯度向下延拓的谱方法

本文提出在平面近似下解算卫星重力梯度向下延拓问题的谱方法，并采用模拟数据进行了试算，结果表明该方法是有效的。这为利用卫星重力梯度数据精化局部重力场提供了可供参考的方法。     

A comparison of methods for the inversion of airborne gravity data

Four integral-based methods for the inversion of gravity disturbances, derived from airborne gravity measurements, into the disturbing potential on the Bjerhammar sphere and the Earths surface are investigated and compared with least-squares (LS) collocation. The performance of the methods is numerically investigated using noise-free and noisy observations, which have been generated using a synthetic gravity field model. It is found that advanced interpolation of gravity disturbances at the nodes of higher-order numerical integration formulas significantly improves the performance of the integral-based methods. This is preferable to the commonly used one-point composed Newton–Cotes integration formulas, which intrinsically imply a piecewise constant interpolation over a patch centered at the observation point. It is shown that the investigated methods behave similarly for noise-free observations, but differently for noisy observations. The best results in terms of root-mean-square (RMS) height-anomaly errors are obtained when the gravity disturbances are first downward continued (inverse Poisson integral) and then transformed into potential values (Hotine integral). The latter has a strong smoothing effect, which damps high-frequency errors inherent in the downward-continued gravity disturbances. An integral method based on the single-layer representation of the disturbing potential shows a similar performance. This representation has the advantage that it can be used directly on surfaces with non-spherical geometry, whereas classical integral-based methods require an additional step if gravity field functionals have to be computed on non-spherical geometries. It is shown that defining the single-layer density on the Bjerhammar sphere gives results with the same quality as obtained when using the Earths topography as support for the single-layer density. A comparison of the four integral-based methods with LS collocation shows that the latter method performs slightly better in terms of RMS height-anomaly errors.     

Geoid determination using one-step integration

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.     

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.     

航空重力向下延拓分组修正的正则化解法

在航空重力向下延拓中,针对病态性对解算结果各个部分影响的不同进行分组修正,提出了分组修正的正则化解法。利用信噪比评估病态性影响,得到将参数分组的方式;依分组修正思想构造正则化矩阵;通过极小化均方误差选取正则化参数。基于EGM2008重力场模型仿真一组航空重力数据,验证了该方法对航空重力数据向下延拓过程的有效性,并与另外3种方法作比较。结果表明,新方法具有更高的精度。