大旋转角坐标变换非迭代公式

针对求解7参数的过程中,经典的线性化最小二乘法因需线性化、迭代及初值以及存在算法耗时出现不收敛现象的问题,该文对无须迭代的7参数坐标变换公式进行了研究。为避免各类参数间的相关性,采用消去法并按照依次求解旋转参数、比例系数和平移参数的顺序解得坐标变换参数。先利用最小二乘法求解旋转参数,然后通过构建目标函数的方式求解比例系数与平移参数,最终得到无须线性化、无须迭代、无须初值的,可用于大旋转角的7参数坐标变换公式。与线性化最小二乘方法进行相比,该方法具有相当的精度及更高的运算效率,可在一定程度上丰富坐标变换理论。     

Changes of spheroid and of datum

Summary A datum change between two geodetic systems with points in common may be derived in three stages; slight adjustments of coordinates to make the networks of common points geometrically similar in the two systems; a scale factor to make them geometrically congruent; finally, an orthogonal transformation to swing them into coincidence. The geometrical concept is developed of a “datum screw”, not arbitrarily chosen as is the “origin” or “datum point” of a geodetic survey, but intrinsic to the geometry. The conditions under which it degenerates to a simple “datum shift” are discussed. Differential and other formulae for changes of spheroid and of datum are given, together with a set of tables of coefficients.     

Nonlinear analysis of the three-dimensional datum transformation [conformal group ℂ7(3)]

 The weighted Procrustes algorithm is presented as a very effective tool for solving the three-dimensional datum transformation problem. In particular, the weighted Procrustes algorithm does not require any initial datum parameters for linearization or any iteration procedure. As a closed-form algorithm it only requires the values of Cartesian coordinates in both systems of reference. Where there is some prior information about the variance–covariance matrix of the two sets of Cartesian coordinates, also called pseudo-observations, the weighted Procrustes algorithm is able to incorporate such a quality property of the input data by means of a proper choice of weight matrix. Such a choice is based on a properly designed criterion matrix which is discussed in detail. Thanks to the weighted Procrustes algorithm, the problem of incorporating the stochasticity measures of both systems of coordinates involved in the seven parameter datum transformation problem [conformal group ℂ₇(3)] which is free of linearization and any iterative procedure can be considered to be solved. Illustrative examples are given. Received: 7 January 2002 / Accepted: 9 September 2002 Correspondence to: E. W. Grafarend     

Transformation of coordinates between two horizontal geodetic datums

The following topics are discussed in this paper: the geocentric coordinate system and its different realizations used in geodetic practice; the definition of a horizontal geodetic datum (reference ellipsoid) and its positioning and orientation with respect to the geocentric coordinate system; positions on a horizontal datum and errors inherent in the process of positioning; and distortions of geodetic networks referred to a horizontal datum. The problem of determining transformation parameters between a horizontal datum and the geocentric coordinate system from known positions is then analysed. It is often found necessary to transform positions from one horizontal datum to another. These transformations are normally accomplished through the geocentric coordinate system and they include the transformation parameters of the two datums as well as the representation of the respective network distortions. Problems encountered in putting these transformations together are pointed out.     

综合不同时期的GPS数据求解转换参数的方法

提出了综合利用不同时期的GPS数据解算三维坐标转换参数的一种方法 ,并结合上海地区 1999年和 2 0 0 2年施测的两个GPS控制网数据 ,验证了给出的模型和方法。结果表明 ,利用本方法是合理而有效的     

已知平面与高程点不同时的严密三维基准变换

由于常规控制网的平面与高程网分开布设,当求解常规控制网与GNSS三维控制网的基准变换参数时,需要分别利用其平面坐标和高程的已知成果。传统控制网与三维GNSS控制网都存在误差,文中在解算变换参数时,同时对两套坐标引入改正数,并顾及公共点与转换点观测误差的相关性,利用公共点的改正数推估转换点的改正数,实现两套坐标的基准变换。模拟数据的分析结果表明,文中方法能够在常规平面与高程控制网分开布设时,实现与GNSS三维控制网的基准变换。     

On symmetrical three-dimensional datum conversion

A 3-D similarity transformation is frequently used to convert GPS-WGS84-based coordinates to those in a local datum using a set of control points with coordinate values in both systems. In this application, the Gauss-Markov (GM) model is often employed to represent the problem, and a least-squares approach is used to compute the parameters within the mathematical model. However, the Gauss–Markov model considers the source coordinates in the data matrix (A) as fixed or error-free; this is an imprecise assumption since these coordinates are also measured quantities and include random errors. The errors-in-variables (EIV) model assumes that all the variables in the mathematical model are contaminated by random errors. This model may be solved using the relatively new total least-squares (TLS) estimation technique, introduced in 1980 by Golub and Van Loan. In this paper, the similarity transformation problem is analyzed with respect to the EIV model, and a novel algorithm is described to obtain the transformation parameters. It is proved that even with the EIV model, a closed form Procrustes approach can be employed to obtain the rotation matrix and translation parameters. The transformation scale may be calculated by solving the proper quadratic equation. A numerical example and a practical case study are presented to test this new algorithm and compare the EIV and the GM models.     

大角度三维基准转换的解析封闭解

传统大地测量应用中的基准转换往往涉及小角度旋转,可只考虑旋转角的一阶量采用线性化方法求解。现代空间测量技术成果应用的基准转换涉及大角度旋转,通过将旋转矩阵所有元素作为未知数并利用旋转矩阵正交条件采用附约束条件平差法迭代求解。本文以空间三维基准转换为例,采用多元模型的矩阵形式将多点坐标组成矩阵处理,并利用旋转矩阵的正交条件导出了大角度三维基准转换的解析分步解。同时引入两套公共点坐标误差对传统三维基准转换模型扩展,导出了同时顾及两套公共点坐标误差的大角度三维基准转换模型的解析解。试验表明:给出的大角度三维基准转换解析解能在实现与传统迭代解等效转换结果的同时,有效避免复杂耗时的迭代计算,提高计算效果。     

三维基准转换的约束加权混合整体最小二乘的迭代解法

基于传统的三维基准转换方法及经典最小二乘（LS）理论,引入混合整体最小二乘（LS-TLS）的原理,在此基础上,根据测量平差原理建立附有约束条件的加权模型,并通过Euler-Lagrange逼近法推导出含有初始值的迭代算法,解决了模型与实际测量条件更加匹配的问题;给出了该模型下三维基准转换的实例,并分析了不同条件下权因素的影响。     

Robust estimation of geodetic datum transformation

The robust estimation of geodetic datum transformation is discussed. The basic principle of robust estimation is introduced. The error influence functions of the robust estimators, together with those of least-squares estimators, are given. Particular attention is given to the robust initial estimates of the transformation parameters, which should have a high breakdown point in order to provide reliable residuals for the following estimation. The median method is applied to solve for robust initial estimates of transformation parameters since it has the highest breakdown point. A smooth weight function is then used to improve the efficiency of the parameter estimates in successive iterative computations. A numerical example is given on a datum transformation between a global positioning system network and the corresponding geodetic network in China. The results show that when the coordinates are contaminated by outliers, the proposed method can still give reasonable results. Received: 25 September 1997 / Accepted: 1 March 1999     

两种大旋转角三维基准转换算法比较研究

分析了两种大旋转角条件下的三维基准转换算法,比较了两种算法的优缺点和使用范围,讨论了选定基准转换所需的公共点及其精度和数量。     

顾及起始数据误差影响时单位权中误差的估值

本文推导了在附有条件的间接平差和附有未知数的条件平差两种方法中,顾及起始数据误差影响时,计算单位权中误差估值(?)的公式,并以两个实例说明,平差时应该考虑起始数据误差对(?)的影响。     

任意椭圆轨道下星载SAR多普勒参数的精确计算

国内外已有文献在推导星载SAR多普勒参数计算公式的过程中,均未考虑椭圆轨道下卫星航迹角不为零的影响,这使得已有公式只能近似地使用在近圆轨道上。如果用于大偏心率轨道,误差较大,将严重降低多普勒参数估计过程的收敛速度。为此,本文将航迹角的影响因素考虑在内,推导出可在任意偏心率轨道下使用的计算公式。针对小偏心率和大偏心率轨道,分别设计了一组仿真试验并进行了对比,结果表明:在大偏心率轨道下,航迹角对多普勒参数确有重要影响。     

Hyperboloidal coordinates: transformations and applications in special constructions

The vector-based algorithms for biaxial and triaxial ellipsoidal coordinates presented by Feltens (J Geod 82:493–504, 2008; 83:129–137, 2009) have been extended to hyperboloids of one sheet. For the backward transformation from Cartesian to hyperboloidal coordinates, of two iterative process candidates one was identified to be well suited. It turned out that a careful selection of the center of curvature is essential for the establishment of a stable and reliable iteration process. In addition, for zero hyperboloidal heights a closed solution is presented. The hyperboloid algorithms are again based on simple formulae and have been successfully tested for various theoretical hyperboloids. The paper concludes with a practical application example on a cooling tower construction.     

The parallel sum of matrices and its application in geodetic adjustments

Summary . It is well known that for the comparison and combination of geodetic networks their heterogeneous datum definitions are well to be considered. Various algorithms have been developed for this purpose. As an alternative concept to deal with hybrid datum problems, the operator parallel sum of matrices is introduced in this paper. To begin with, a definition is given and some basic properties are explained. To demonstrate the usefulness of the operator, two practical applications are given. The first deals with the estimation of parameters describing the deformation of two networks which are to be compared to each other. The second one treats the estimation of parameters representing the heterogeneous datum definitions of two networks which are to be merged into a hybrid network. It will be shown that – regardless of the datum definitions of the preadjusted individual networks – the parallel sum of matrices can be used to simplify the algorithms applied for the estimation of those parameters. Received 15 December 1995; Accepted 24 September 1996     

永久性潮汐与大地测量基准

社会、经济和科技的发展,对大地测量基准的精度、历元、框架等方面的要求愈来愈高,因此必须顾及具有同等精度量级的地球动力现象如潮汐、极移、地壳运动、冰后期反弹等对大地测量基准的影响。本文仅讨论永久性潮汐对重力基准、高程基准、大地水准面、垂线偏差等方面的影响。对大地测量基准中使用的无潮汐、平均湖汐和零潮汐等概念和相关计算进行了推导和讨论。     

适用于大旋转角的三维基准转换的一种简便模型

基于大角度的空间直角坐标转换，提出了以方向余弦为参数、适用于任何角度旋转的空间直角坐标转换的简明计算方法。     

地心坐标反解大地坐标的非线性方程组牛顿迭代法分析

目前的一些地心坐标向大地坐标的转换模型在计算速度,稳定性或精度方面存在一定的局限性。本文探讨了非线性方程组数值迭代的求解方法及其在MATLAB 7.0中的编程实现,并将结果与基于一元三次方程求解的严密方法做了比较分析,结果证明该方法是可以在实际中参考应用的。     

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

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

Error analysis of the 3D similarity coordinate transformation

The 3D similarity coordinate transformation with the Gauss–Helmert error model is investigated. The first-order error analysis of an analytical least-squares solution to this problem is developed in detail. While additive errors are assumed in the translation and scale estimates, a 3 × 1 multiplicative error vector is defined to effectively parameterize the rotation matrix estimation error. The propagation of the errors in the coordinate measurements to the errors in the estimated transformation parameters is derived step-by-step, and the formulae for calculating the variance–covariance matrix of the estimated parameters are presented.