空间直角坐标计算大地坐标的抛物线逼近法

采用抛物线逼近法求解大地纬度和大地高,先计算空间点在椭球面上的子午面坐标,然后求解点的大地纬度和大地高。     

以地心纬度为变量的常用纬度正反解符号表达式

借助具有强大符号运算功能的计算机代数系统Mathematica,推导了地图投影学中等距离纬度、等角纬度、等面积纬度与地心纬度之间的正反解直接展开式,并将式中的系数统一表示成关于椭圆偏心率e和椭球第三扁率n的幂级数形式。与以往反解方法不同的是,采用符号迭代法进行以地心纬度为变量的等距离纬度、等角纬度、等面积纬度的反解,并使用最大差异值作为衡量精度的标准。算例分析表明,以地心纬度为变量的常用纬度展开式在结构和形式上与以大地纬度为变量的辅助纬度保持一致,基于第三扁率n的幂级数表达式具有更紧凑的形式和更好的收敛性,其直接展开式的精度分别优于(1×10^-8)″和(1×10^-10)″,可以满足大地测量和地图投影精密计算的需要。     

辅助纬度与大地纬度间的无穷展开

利用无穷级数理论和拉格朗日反演定理,详细推导了大地测量和制图学中常用的辅助纬度与大地纬度间的无穷展开,主要表现为参考椭球第一偏心率的幂级数形式。通过建立一系列严格的系数递推公式,得到了等量纬度反解展开式和等角纬度反解展开式;同时,推导了古德曼函数的泰勒展开式,进而得到了等角纬度正解展开式;利用级数除法公式,得到了等距离纬度正解展开式系数的行列式表示。通过比较本文方法与计算机代数系统Mathematica直接推导求得的辅助纬度正反解展开式e^0～e⁴⁰阶系数和相应的程序用时,表明本文算法是正确的、快速的。以CGCS2000参考椭球为例,对辅助纬度正反解进行了算例分析,也进一步验证了本文公式的正确性。     

等面积纬度函数和等量纬度变换的直接解算公式

为实现等面积投影和等角投影间的直接变换,借助计算机代数系统Mathematica,推导出了等面积纬度函数和等量纬度变换的直接解算公式,并将式中系数统一表示为椭球第一偏心率的幂级数形式,可解决不同参考椭球下的变换问题。算例分析表明,本导出公式的计算误差分别小于10 m2和10-4(″),可供实际使用。     

GPS网WGS-84平差坐标向地方独立坐标的转换

论述了采用椭球膨胀法确定区域性椭球面的方法,并给出了有关的大地经纬度计算公式,结合具体工程实例,建立了地方独立坐标系,并对相关结果与数据进行了比较和分析。     

等分分层组合投影研究

针对中国地理格网(1°、10°等多级格网系统)的分割方法,设计了一种适合该格网系统的新型地图投影——分层组合投影。从微分几何的观点出发,把地球椭球按等纬度分割成若干层圆台,分别建立每个圆台的投影模型,即可得到一种地图投影。这种投影还可根据格网间隔的不同进行细分,从而发展成为一种适合多分辨率格网模型的动态地图投影。通过对该投影进行变形计算表明,该投影可以保持等角,而且面积和长度变形都很小,特别是在高纬度地区,与Mercator投影相比变形明显减小。     

常用纬度与归化纬度差异极值分析表达式

从常用纬度与归化纬度的定义出发,借助于计算机代数分析,对常用纬度与归化纬度差值进行了系统的分析,推导出了常用纬度与归化纬度的差异极值分析表达式,并将式中系数展开为椭球第一偏心率e和第三扁率n的幂级数形式.研究结果表明,常用纬度与归化纬度差值极值点均位于45°附近,基于n展开的分析表达式比基于e展开的分析表达式形式上更加...     

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.     

The solution of the Korn–Lichtenstein equations of conformal mapping: the direct generation of ellipsoidal Gauß–Krüger conformal coordinates or the Transverse Mercator Projection

The differential equations which generate a general conformal mapping of a two-dimensional Riemann manifold found by Korn and Lichtenstein are reviewed. The Korn–Lichtenstein equations subject to the integrability conditions of type vectorial Laplace–Beltrami equations are solved for the geometry of an ellipsoid of revolution (International Reference Ellipsoid), specifically in the function space of bivariate polynomials in terms of surface normal ellipsoidal longitude and ellipsoidal latitude. The related coefficient constraints are collected in two corollaries. We present the constraints to the general solution of the Korn–Lichtenstein equations which directly generates Gau?–Krüger conformal coordinates as well as the Universal Transverse Mercator Projection (UTM) avoiding any intermediate isometric coordinate representation. Namely, the equidistant mapping of a meridian of reference generates the constraints in question. Finally, the detailed computation of the solution is given in terms of bivariate polynomials up to degree five with coefficients listed in closed form. Received: 3 June 1997 / Accepted: 17 November 1997     

Taylor expansion of the theoretical gravity formula

A recurrence relation has been derived to obtain the derivatives required for the Taylor expansion of the theoretical gravity formula in powers of latitude. The computations using the relationship derived can provide easily all required derivatives for any reference ellipsoid.     

Global optimization of the Gauss conformal mappings of an ellipsoid to a sphere

The Gauss conformal mappings (GCMs) of an oblate ellipsoid of revolution to a sphere are those that transform the meridians into meridians, and the parallels into parallels of the sphere. The infinitesimal-scale function associated with these mappings depends on the geodetic latitude and contains three parameters, including the radius of the sphere. Gauss derived these constants by imposing local optimum conditions on certain parallel. We deal with the problem of finding the constants to minimize the Chebyshev or maximum norm of the logarithm of the infinitesimal-scale function on a given ellipsoidal segment (the region contained between two parallels). We show how to solve this minimax problem using the intrinsic function fminsearch of Matlab. For a particular ellipsoidal segment, we get the solution and show the alternation property characteristic of best Chebyshev approximations. For a pair of points relatively close in the ellipsoid at different latitudes, the best minimax GCM on the segment defined by these points is used to approximate the geodesic distance between them by the spherical distance between their projections on the corresponding sphere. This approach, combined with the best locally GCM if the points are on the same parallel, is illustrated by applying it to some case studies but specially to a 10° × 10° region contained between portions of two parallels and two meridians. In this case, the maximum absolute error of this spherical approximation is equal to 2.9 mm occurring at a distance about 1,360 km. This error decreases up to 0.94 mm on an 8° × 8° region of this type. So, the spherical approximation to the solution of the inverse geodesic problem by best GCM can be acceptable in many practical geodetic activities.     

GNSS相应坐标系之间的关系及其影响

不同的GNSS采用的坐标系定义几乎相近,但参考椭球及其坐标实现不同,这将影响多GNSS融合导航定位效果。根据各GNSS坐标系所采用参考椭球的基本常数,计算比较了不同坐标系参考椭球参数的差异;导出了相应的正常重力公式,比较了这些正常重力公式确定的正常重力值差异;最后分别从坐标系统的定义与实现两个方面分析了其对定位结果的影响。结果表明:1)GPS(BDS)与Galileo和GLONASS所使用的参考椭球引起正常重力差约为0.15和0.30 mgal;2)GPS与BDS,Galileo及GLONASS所使用参考椭球引起纬度分量最大差异约为0.1 mm,3 cm和3 cm,高程分量约为0.1 mm,0.5 m和1 m;3)各GNSS所使用坐标框架间转换参数引起的坐标变化达到厘米级。     

Models for Fitting Gravimetric Geoids and GPS Results

The aim of this investigation is to study how to use a gravimetric(quasi) geoid for levelling by GPS data in an optimal way.The advent of precise geodetic GPS has made the use of a technique possible,which might be called GPS- gravimetric geoid determination.In this approach,GPS heights above the reference ellipsoid are determined for points whose levelled (orthometric) height H is above sea level people have already surveyed;for these points,we thus have the values of the geoid undulation N.These values are then used to constrain the geoid undulations N‘ obtained from the gravimetric solution.     

等角圆锥投影基准纬度非迭代算法

为简化传统正轴等角圆锥投影求解基准纬度时繁琐的迭代算法,引入平均纬度和平均纬差的概念,借助计算机代数系统Mathematica,在平均纬差处级数展开,导出了基于球体模型的正轴等角圆锥投影求解基准纬度的非迭代算法。以全国和不同纬差的省区为例,将其与传统椭球迭代算法进行对比分析。结果表明,推导的基于球体模型的非迭代公式计算基准纬度B0、B1、B2的相对误差最大值为2.011%,长度变形的相对误差小于1×10-6,基本可满足全国以及各省区地图制图的精度要求,从而验证了所研究算法的精确性与实用性。     

Transforming ellipsoidal heights and geoid undulations between different geodetic reference frames

Transforming height information that refers to an ellipsoidal Earth reference model, such as the geometric heights determined from GPS measurements or the geoid undulations obtained by a gravimetric geoid solution, from one geodetic reference frame (GRF) to another is an important task whose proper implementation is crucial for many geodetic, surveying and mapping applications. This paper presents the required methodology to deal with the above problem when we are given the Helmert transformation parameters that link the underlying Cartesian coordinate systems to which an Earth reference ellipsoid is attached. The main emphasis is on the effect of GRF spatial scale differences in coordinate transformations involving reference ellipsoids, for the particular case of heights. Since every three-dimensional Cartesian coordinate system ‘gauges’ an attached ellipsoid according to its own accessible scale, there will exist a supplementary contribution from the scale variation between the involved GRFs on the relative size of their attached reference ellipsoids. Neglecting such a scale-induced indirect effect corrupts the values for the curvilinear geodetic coordinates obtained from a similarity transformation model, and meter-level apparent offsets can be introduced in the transformed heights. The paper explains the above issues in detail and presents the necessary mathematical framework for their treatment. An erratum to this article can be found at     

采用双重投影法的椭球面日晷投影

研究了从地球椭球面描写到海图平面的采用双重投影法的日晷投影方程式,并阐明了椭球面日晷投影的 3个重要特性及其应用。     

基于VB6.0的高斯投影坐标转换的实现

介绍了高斯投影坐标转换的方法,包含坐标的正算和反算。具体来说就是经纬度坐标（B,L）转换为本椭球系的平面直角坐标（x,y）,以及平面直角坐标（x,y）转换为相应椭球系的经纬度坐标（B,L）。本文还介绍了转换软件的开发过程、功能及其转换精度的验证。     

基于椭球不确定性的平差模型与算法

测量平差模型中的参数通常存在一些不确定的附加信息或先验信息,充分利用它们可以对部分参数进行约束,从而保证参数解的唯一性和稳定性。本文利用椭球集合描述不确定性,建立了一个新的带有椭球不确定性的平差模型。以两个椭球交集的外接椭球的特征矩阵的迹最小平差准则,分析了不确定度的传播规律,给出了带有椭球不确定性的平差方法。最后,通过算例验证了算法的有效性,说明了平差解与带权混合估计的关系。     

不同变形性质正轴圆柱投影和正轴圆锥投影间的直接变换

为避免不同变形性质正轴圆柱投影和正轴圆锥投影间传统间接变换繁琐的计算过程,利用子午线弧长、等量纬度和等面积纬度函数间变换的直接展开式,建立了相应投影坐标间的直接变换模型,无需计算大地纬度即可完成变换。本文导出公式均为含参考椭球第一偏心率的符号形式,可解决两类投影在不同参考椭球下的变换问题。算例分析表明与传统间接变换模型相比,本文建立的直接变换模型提高了计算效率和计算精度,可供实际使用。     

正常重力场的确定以及相关的一个理论问题

给定参考椭球体的四个基本参数,椭球体产生的外部正常重力位场可利用传统方法确定。本文基于WGS84参考椭球的四个基本参数和压缩恢复法,确定了一个定义在半径为6000km内部球的外部正则调和的虚拟正常场。在mm级精度水平下,所确定的虚拟正常场在椭球体外部与真实正常场一致。基于虚拟正常场,可解决物理大地测量学中存在的一个无定义问题。