测地主题正反解解算

由测地坐标系中大地线的微分方程式推导出其微分关系式，得出在地球椭球面上基于测地坐标进行正反解的算法和公式，它与大地主题解算公式相比，更为简捷明了。由实际计算数据表明，对于100km以下的距离解算，它亦能达到相当高的精度。因此测地坐标的点位表述不仅可用于DEM和GIS三维可视化，也可用于三维GIS建模以及空间度量和分析。     

引力场梯度张量的非奇异公式推导

传统的引力场梯度计算公式在两极附近存在奇异性,需要换用其他的非奇异计算公式。从奇异性产生的原因入手,并结合勒让德函数的有关性质,推导了一组新的计算公式。实际计算验证了该公式的正确性和有效性。     

General non-iterative solution of the inverse and direct geodetic problems

Summary Improved practical and theoretical formulas are presented for the calculation of geodetic distances, azimuths, and positions on a spheroid. The formulas are designed for use with either electronic computers or desk calculators. For the latter, the formulas lend themselves to the construction of useful interpolation tables. The report includes convenient computation forms and auxiliary equations which assure a high degree of accuracy for any geodetic line, no matter how short or long (up to half or fully around the earth) and regardless of its orientation or location. Numerical examples illustrate the complete calculation procedure.     

Numerical Evaluation of the Robinson Projection

The Robinson projection is one of the most preferred projections for world reference maps in atlas cartography. The projection is constructed from Robinson's look-up table since there are no analytical formulas. This deficiency has led to a number of requests for the plotting formulas to which cartographers have responded by deriving analytical equations using different interpolation algorithms applied to Robinson's table values. The Robinson projection was examined with regard to its deformations calculated by four different algorithms, including the multiquadratic method. The numerical evaluations were then used to compare the algorithms. Solutions have been presented including some criticisms about this projection. The latitudes along which the scale is true and on which the maximum angular distortion equals zero have been determined.     

Analytical computation of gravity effects for polyhedral bodies

On the basis of recent analytical results we derive new formulas for computing the gravity effects of polyhedral bodies which are expressed solely as function of the coordinates of the vertices of the relevant faces. We thus prove that such formulas exhibit no singularity whenever the position of the observation point is not aligned with an edge of a face. In the opposite case, the contribution of the edge to the potential to its first-order derivative and to the diagonal entries of the second-order derivative is deemed to be zero on the basis of some claims which still require a rigorous mathematical proof. In contrast with a common statement in the literature, it is proved that only the off-diagonal entries of the second-order derivative of the potential do exhibit a noneliminable singularity when the observation point is aligned with an edge of a face. The analytical provisions on the range of validity of the derived formulas have been fully confirmed by the Matlab $^{\textregistered }$ program which has been coded and thoroughly tested by computing the gravity effects induced by real asteroids at arbitrarily placed observation points.     

基于DOA和TOA的机载单站无源定位

研究可对地面固定目标进行定位跟踪的机载单站无源定位方法,建立二维等效定位模型,利用DOA和TOA信息建立状态方程及测量方程,经过扩展卡尔曼滤波算法进行定位及跟踪。通过计算机仿真结果验证,该方法稳定、有效,且可达到相当的精度。     

基于左右正射影像的三维量测方法的研究

介绍了利用左右正射影像上高层建筑物存留的投影差,形成立体,在正射影像立体观察的环境下,推导出一套对地物进行三维量测的公式。实验表明,该公式严密准确,不仅可以使用户直观地检查DEM的正确与否,而且可以实时地获取地物的三维坐标,以大大减少用户对于向量的需求。     

基于逆Vening-Meinesz公式的测高重力中央区效应精密计算

为提高利用逆Vening-Meinesz公式反演测高重力中央区效应的精度,视中央区为矩形域,将垂线偏差分量表示成双二次多项式插值形式,引入非奇异变换,推导出了重力异常的计算公式。以低纬度区域2′×2′的垂线偏差实际数据为背景场进行了计算,结果表明,当中央区包含4个网格时,传统公式与推导出的重力异常计算公式误差的最大值大于1 mGal。推导出的公式可为高精度测高重力中央区效应的计算提供理论依据。     

太阳光压模型中的偏导数

用矢量微分推导出光压摄动加速度对卫星位置矢量的偏导数公式，并用数值微分方法检验了公式的正确性。     

大地主题问题的非迭代新解

首先对大地主题问题解算中的符号体系进行改进与统一,使之更符合大地测量的使用习惯;在此基础上讨论非迭代的大地主题问题的正反解,借助计算机代数系统Mathematica推导更准确的正解公式,得到非迭代的反解公式,与以往的解算相比该公式具有更简明的数学表达形式。该过程丰富和完善了大地主题问题的解算体系。最后进行正反解的数据验证。     

经纬仪工业测量系统直接坐标传递后的点位坐标及精度分析

对经纬仪工业测量系统应用于亚毫米级的大型工件的测量 ,提出了一种快速坐标转换的方法。利用定向参数作为坐标转换参数 ,导出了在统一坐标系下点位坐标的计算公式 ,并对转换后的坐标进行了精度分析。     

不同国家坐标系下子午线收敛角计算公式及其精度分析

介绍了子午线收敛角的概念和计算公式,并用C＋＋语言对高斯坐标正反算和子午线收敛角的计算公式进行了编译,提高了计算效率。选取徐州地区14个点进行试算,发现1954国家大地坐标系和1980国家大地坐标系在徐州地区可以互用收敛角计算公式,其误差对于高精度工程测量可以忽略不计。但是,通过运用误差传播定律对子午线收敛角公式分析发现,对于投影带边缘、纬度较高的地区计算子午线收敛角时,要区分使用。     

一类球谐函数与三角函数乘积积分的计算

本文根据球谐函数的跨次递推公式和三角函数的性质,详细推导了在重力梯度调和分析中出现的一类球谐函数积分的跨次递推公式和递推初始值的计算公式。数值试验表明,球谐函数跨次递推算法具有快速、稳定的优点。该类积分的跨次递推实现,为卫星重力梯度调和分析奠定了算法基础。     

双等距离投影及其应用

双等距离投影是指自图上两定点出发至图上任意点间的距离都与实地等长的投影。本文先从球面极坐标出发，讨论了双等距离投影公式的建立，继而给出了在地理坐标下的正轴和斜轴投影计算方法。本文证明了该投影的性质属任意投影，其计算方法是可行的，有实用价值。双等距离投影在测量上可用于圆圆系统定位，在军事上可用于敌情监测等。     

附加约束条件的序贯平差

从观测量分组的参数平差出发,详细推导了附加约束条件的序贯平差及其精度估计,讨论了附加约束条件的抗差序贯平差的计算过程.算例表明,推导的公式正确有效,简单实用,附加约束条件的序贯平差抗差估计能有效地抵制粗差的影响.     

Gravity recovery using COSMIC GPS data: application of orbital perturbation theory

 COSMIC is a joint Taiwan–US mission to study the atmosphere using the Global Positioning System (GPS) occultation technique. Improved formulas are developed for the radial, along-track, and cross-track perturbations, which are more accurate than the commonly used order-zero formulas. The formulas are used to simulate gravity recovery using the geodetic GPS data of COSMIC in the operational phase. Results show that the EGM96 model can be improved up to degree 26 using 1 year of COSMIC data. TOPEX/POSEIDON altimeter data are used to derive a temporal gravity variation. COSMIC cannot reproduce this gravity variation perfectly because of data noise and orbital configuration, but the recovered field clearly shows the gravity signature due to mass movement in an El Ni?o. Received: 3 March 2000 / Accepted: 10 November 2000     

BAYES方差分量估计在虚拟权平差中的应用

本文把近似Bayes方基分量估计公式应用于矿区扩建网的虚拟权平差中，解决了具体计算中的一些实际问题，最后，通过一矿区扩建网平差实例验证了该方法的可行性．     

On the evaluation of the gravity effects of polyhedral bodies and a consistent treatment of related singularities

We show that the singularities which can affect the computation of the gravity effects (potential, gravity and tensor gradient fields) can be systematically addressed by invoking distribution theory and suitable formulas of differential calculus. Thus, differently from previous contributions on the subject, the use of a-posteriori corrections of the formulas derived in absence of singularities can be ruled out. The general approach presented in the paper is further specialized to the case of polyhedral bodies and detailed for a rectangular prism having a constant mass density. With reference to this last case, we derive novel expressions for the related gravitational field, as well as for its first and second derivative, at an observation point coincident with a prism vertex and show that they turn out to be more compact than the ones reported in the specialized literature.     

Spherical integral formulas for upward/downward continuation of gravitational gradients onto gravitational gradients

New integral formulas for upward/downward continuation of gravitational gradients onto gravitational gradients are derived in this article. They provide more options for continuation of gravitational gradient combinations and extend available mathematical apparatus formulated for this purpose up to now. The starting point represents the analytical solution of the spherical gradiometric boundary value problem in the spatial domain. Applying corresponding differential operators on the analytical solution of the spherical gradiometric boundary value problem, a total of 18 integral formulas are provided. Spatial and spectral forms of isotropic kernels are given and their behaviour for parameters of a GOCE-like satellite is investigated. Correctness of the new integral formulas and the isotropic kernels is tested in a closed-loop simulation. The derived integral formulas and the isotropic kernels form a theoretical basis for validation purposes and geophysical applications of satellite gradiometric data as provided currently by the GOCE mission. They also extend the well-known Meissl scheme.     

等角投影变换的常系数公式及其在高斯—克吕格投影换带中的应用

讨论了等角投影变换的常系数一般公式及应用模型,墨卡托投影和高斯－克吕格投影问题的正解变换及其在高斯－克吕格投影换带中的应用,常系数计算公式优于传统的变系数计算公式,是基于计算机的等角投影变换的最佳模型,它在计算机制图,地理数据库,ＧＩＳ等领域中有着广泛的应用。