The adjustment of intersecting chains of triangulation using a medium sized computer

Summary A method of adjusting intersecting chains of triangulation by direct elimination, using variation of coordinates, is given which is suitable for a computer of modest storage capacity (about 4000–8000 words). The network is considered to be made up of “nodes” (where two or more chains intersect) and “links” (which connect two nodes). The normal equations for the points in each link are computed, the coefficients of the unknown coordinates of the associated nodes being treated as extra right hand sides. These equations are then solved to express the coordinates of any link station in terms of the associated nodes. The process is repeated for all links. The normal equations for the nodes are then set up and substitution made for all “link” terms. The method is particularly powerful for the adjustment of a few, long, interconnected chains since the reduced normal equation coefficients are then banded about the diagonal, the semi band width being a little greater than the largest number of unknowns in the nodes along a single tie chain. For a computer of 7000 words capacity it is possible to solve for a band width of over 100 unknowns i.e. more than 50 stations.     

基于史赖伯规则的机载InSAR区域网平差

采用区域网平差的方法进行合成孔径雷达干涉测量（interferometric synthetic aperture radar,InSAR）干涉参数定标,是减小大区域、多个干涉像对重叠处反演高程差异的有效方法。为了提高InSAR区域网平差的答解效能,设计一种利用史赖伯规则的机载InSAR区域网平差干涉参数定标方法。该方法利用史赖伯规则,由连接点误差方程式构建等效误差方程式,可有效减小法方程系数矩阵大小,降低对计算机配置的要求,提高平差的答解效能。采用我国机载双天线InSAR数据进行干涉参数定标试验,改化前后两种方案,干涉参数定标结果没有明显变化,改化后的平差耗时均值减小了,验证了利用史赖伯规则的机载InSAR区域网平差干涉参数定标的正确性和有效性。     

The orthometric height and the holonomity problem

When height networks are being adjusted, many geodesists advocate the approach where the adjustment should be done by using geopotential numbers rather than the orthometric or normal heights used in practice. This is based on a conviction that neither orthometric nor normal heights can be used for the adjustment because these height systems are not holonomic, meaning–among other things–that height increments (orthometric or normal) when summed around a closed loop do not sum up to zero. If this was the case, then the two height systems could not be used in the adjustment; the non-zero loop closure would violate the basic, usually unspoken, assumption behind the adjustment, namely that the model claiming that height differences are observable is correct. In this paper, we prove in several different ways that orthometric and normal heights are theoretically just as holonomic as the geopotential numbers are, when they are obtained from levelled height differences using actual gravity values. This disposes of the argument that geopotential numbers should be used in the adjustment. Both orthometric and normal heights are equally qualified to be used in the adjustment directly.     

利用二维DLT及光束法平差进行数字摄像机标定

利用共线方程和二维DLT之间的对应关系导出了由二维DLT的8个参数表达的主纵线方程式，讨论了主点初值和求解方法，给出了单张像片摄像机参数分解不惟一性及临界序列的证明，详细推导了利用二维DLT参数分解摄像机外方位元素初值的实用算法，论述了利用光束法平差进行摄像机标定的数学模型，实际图像数据实验取得了很好的结果，验证了本文所提出的摄像机标定算法的可行性。     

我国天文大地网与1980年大地坐标系的建立及应用

全面论述了我国天文大地网整体平差和１９８０年大地坐标系的建立及其应用，包括：第一，我国天文大地网的结构和参数整体平差的规模。第二，参考椭球的选择和定位，长度基准；垂线偏差和高程异常的计算；在椭球面上进行平差的问题；平差中拉普拉斯方位角的处理；各级大地网整体统一平差的问题。第三，用于我国天文大地网整体平差中的条件联系平差法的基本原理和它的几个特殊问题的处理方法；一种适用于电子计算机计算的大规模法方程     

Total least squares adjustment in partial errors-in-variables models: algorithm and statistical analysis

The weighted total least squares (TLS) method has been developed to deal with observation equations, which are functions of both unknown parameters of interest and other measured data contaminated with random errors. Such an observation model is well known as an errors-in-variables (EIV) model and almost always solved as a nonlinear equality-constrained adjustment problem. We reformulate it as a nonlinear adjustment model without constraints and further extend it to a partial EIV model, in which not all the elements of the design matrix are random. As a result, the total number of unknowns in the normal equations has been significantly reduced. We derive a set of formulae for algorithmic implementation to numerically estimate the unknown model parameters. Since little statistical results about the TLS estimator in the case of finite samples are available, we investigate the statistical consequences of nonlinearity on the nonlinear TLS estimate, including the first order approximation of accuracy, nonlinear confidence region and bias of the nonlinear TLS estimate, and use the bias-corrected residuals to estimate the variance of unit weight.     

水准网闭合差自动解算的新方法

介绍一种在水准测量中独立闭合环或附合路线自动搜索及闭合差自动解算的新方法。该方法将间接平差函数模型转化为条件平差函数模型,利用条件方程式的所有信息,搜索出独立闭合环或附合路线并快速解算出闭合差,与给定的限差进行比较,较好地检验并剔除粗差。最后经过编程实现该方法,用实例验证了该方法的正确性和用之进行观测值质量检核的有效性。     

Local crossover analysis of exactly repeating satellite altimeter data

The crossover adjustment plays a central role in the processing of satellite altimeter measurements. The usual procedure is to form sea surface height differences at crossover points, solve for the radial orbit error (with due attention to the singular nature of the estimation problem) and then to construct altimetric sea-level maps using the mean sea surface heights at the crossover points. Our approach is very different, to make direct use of measurements at crossover points without differencing and to estimate simultaneously orbit parameters, mean sea surface height and sea surface height variability in a single, unified adjustment. The technique is suited for repeat data over an area small enough that adjoining passes may be considered to be parallel and to permit the solution of a set of linear equations of dimension equal to the number of crossover points. The size of the numerical problem is almost independent of the number of repeat cycles of the altimeter mission. Explicit recognition is given to the rank defect of the least-squares estimation problem; we show that, for an orbit model with r parameters, the rank defect of the local crossover problem is exactly r ². The defect may be overcome by choosing an appropriate set of constraints – either giving a best fit of mean sea surface heights to a reference surface, or minimising orbit parameters, or a minimum norm solution in which both mean sea surface heights and orbit parameters are minimised. There is no need to choose a reference pass, all passes are treated equally and data gaps are easily accommodated. Numerical results are presented for the south-western Indian Ocean, based on the first 2 years of altimeter data from the Geosat Exact Repeat Mission. Received: 31 May 1996 / Accepted: 19 April 1997     

THE HEIGHT ADJUSTMENT OF A SINGLE STRIP

This note represents an extension of some earlier work on the plan adjustment of a fully controlled strip (Thompson, 1967). A similar method was developed for the height adjustment. Unlike other adjustment methods using polynomials, both plan and height adjustment formulae are very simple, so that adjustment can be rapidly carried out by hand computation.     

薄板样条与大区域高程异常插值

精确确定高程异常是将大地高转换成正常高的基础,也是制约GPS应用于高程测量的关键。目前诸多高程异常的内插或拟合方法都有其局限性与适用条件,尤其对于大区域范围内的高程异常插值没有一个好的办法。本文利用薄板样条在二维曲面插值中的优越性,提出了大区域高程异常插值的薄板样条法,简要介绍了薄板样条理论,给出了薄板样条的实用计算方法,并用模拟算例与样条法进行对比,表明了该方法的可行性与优越性。     

Vertical velocities of European coastal sites derived from continuous GPS observations

Vertical velocities of 30 European permanent Global Positioning System (GPS) stations at or close to tide gauge sites are estimated from more than 3 years of continuous observations. The results of two different solution strategies are presented and compared. The first approach accumulates the daily free network normal equations, the second introduces all sets of daily ellipsoidal height estimates and their covariance matrix into a subsequent common least squares adjustment. In both solutions, mean station heights at a reference epoch, linear vertical velocities, height discontinuities and short period height offsets are estimated. The second approach solves in addition for periodic annual signals and for site-specific pressure loading coefficients. The vertical velocities range from +8 mm/year in the center of the Fennoscandian uplift area to –4 mm/year at a few subsiding locations. Apart from these extrema, most of the sites experience only very small vertical motions. The standard deviations from the second approach providing more realistic error estimates are well below 0.15 mm/year. Some specific data problems are discussed.     

基于Matlab的GPS网间接平差教学设计与实现

GPS网间接平差是测绘类本科专业误差理论与测量平差基础教学的重要组成部分。本文针对GPS网间接平差存在的观测数据多、计算量大、在实际课堂教学中很难讲解清楚等问题,提出基于Matlab的GPS网间接平差教学设计与实现方法。给出GPS网间接平差课堂教学的内容及难点,阐述基于Matlab的GPS网间接平差程序设计流程及关键步骤的实现方法。结合3个具体的GPS网平差实例,展示程序运行的结果。结果表明,基于Matlab进行GPS网间接平差教学既有助于加深学生对GPS网间接平差理论知识的理解,培养学生理论联系实际的能力,又可以实现学生所学计算机高级编程语言与专业知识的有机结合,提高学生的编程能力和学习兴趣。     

SAR影像几何校正

SAR影像的几何校正是阻碍其应用的瓶颈问题，目前所用的方法都有其局限性。介绍描述SAR构像几何的坐标系统及距离－多普勒方程，同时，介绍卫星轨道及描述卫星轨道的参数。为提高几何校正的精度，提出一种由几个离散点的值内插整个弧段上卫星状态向量的轨道模拟算法。在此基础上，给出利用DEM进行SAR影像几何校正的方法。该方法以时间参数为自变量，迭代解算距离－多普勒方程，计算DEM中点每一点对应的影像坐标，最后进行了实验与讨论。     

THE USE OF AIRBORNE AUXILIARY DATA IN THE RIGOROUS LEAST SQUARES ADJUSTMENT OF A BLOCK OF AERIAL TRIANGULATION

The problem of incorporating airborne auxiliary data into an aerial triangulation block adjustment is considered. The conclusion is reached that this problem can be most readily solved by the use of a fully rigorous adjustment method. The first order adjustment equations for all types of airborne auxiliary control are developed and the method of incorporating these equations into the block solution is demonstrated. The general computing problems of the aerial triangulation block adjustment are briefly discussed and a method outlined for the incorporation of tie strips, at any angle across the block, into the adjustment without unduly increasing the storage requirements for the solution. The problem of correctly weighting the different types of available data is discussed. The paper concludes with a demonstration of the improved accuracy possible if precise aerial triangulation data are combined with apparently crude absolute airborne data.     

等价条件平差模型的方差-协方差分量最小二乘估计方法

提出等价条件闭合差的方差-协方差分量最小二乘估计方法，简称LSV-ECM法。首先，利用等价条件平差模型建立了基于等价条件闭合差二次型的方差-协方差分量估计方程，由矩阵半拉直算子将其变换为线性Gauss-Markov形式，进而通过最小二乘准则导出了具有模型通用性、形式简洁性且满足无偏性和最优性的方差-协方差分量估计公式。其次，证明了LSV-ECM方法与残差型VCE方法的等价性，并在此基础上通过计算复杂度定量分析了所提方法的计算高效性。最后，通过边角网平差和中国区域GNSS站坐标时序建模及其结果分析，验证了所提新方法的正确性和计算高效性。     

Improvement of parameter accuracy by choice and quality of observation

In a least squares adjustment (a minimum variance solution) using the technique of variation of coordinates (observation equations), a key result is the co-variance (dispersion) matrix of parameters. Assuming that standard errors of observations are used in the formation of the normal equations, rather than relative weights, this dispersion matrix gives the estimates of standard errors for the parameters solved for in the adjustment. A method will be presented which allows the designer of the observing plan to alter this dispersion matrix, which may not meet user requirements, so that it will meet user requirements and, from its inverse, solve mathematically for the selection and quality (accuracy) of the observations required to form this altered dispersion matrix of parameters.     

方差分量估计的MINQE理论及其在摄影测量中的应用

本文深入研究了异方差分量模型的MINQE(U,I)准则的算法实现的可能性,所导出的约化算法具有压缩内存,提高运算速度等明显优点,因而具有实用意义。文中介绍了实现约化算法的计算机程序框图,并提供了数据例子。作者将方差分量估计的MINQE理论引入解析摄影测量平差实践,分别采用模拟数据和新滩滑坡监测数据研究了自检校光束法平差中随机模型的确定,结果表明,MINQE理论用于自检校的验后权估计是可行的。最后,作者就平差随机模型的扩展,以及MINQE理论的应用提出了进一步研究的几点设想。     

整体大地测量的确定性参数模型

本文将整体大地测量各类观测方程中作为随机信号的重力场参数（ΔN，ξ，η）均表达为重力异常Δg的泛函，又引入Δg的B样条函数并以样条系数C_(11)作待定参数，提出一种整体大地测量的确定性参数模型，给出了相应系数矩阵的表达式和算法，最后对该模型的求解问题作了说明。     

简便实用的GPS网平差模型

采用数值微商求偏导数的方法，从而使ＧＰＳ网平差的计算模型非常简洁。选用平面坐标和大地高作为平差参数，可以方便地简化为平面网或高程网，加入地面归心数据和常规观测数据也变得十分容易。对旋转角大的地方独立坐标系，提出了处理办法。经对实测数据的处理，证明了模型的正确性     

An Implementation of the Robinson Map Projection Based on Cubic Splines

The Robinson world map projection has been in existence since 1963. Mapping equations for the projection are based on table interpolation. Cubic splines with the choice of appropriate boundary conditions are shown to possess favorable characteristics as interpolative functions. The evaluation of cubic splines for given latitudes provides a basis for efficient forward mapping equations. Inverse mapping equations are based on an efficient inversion of the cubic spline functions utilizing Newton's method. Derivatives of the cubic spline functions lead to formulas for scale factors along projected meridians and parallels and for areal scale factors. Details for a computer implementation of the mapping equations and the computation of scale factors are given.