大地测量学 总结

大地测量技术的任务: 高精度数据采集、量测

测量学与大地的不同：测量学的范围小，观测基准是铅垂线和平面，并且认为铅垂线是平行的，计算基准是垂线和平面。

大地范围是整个地球，观测基准是铅垂线，但铅垂线是是不平行的。计算基准是参考椭球的法线和大地水准面。

大地测量学：研究地球形状及行星几何和物理形态特征及其变化规律的基础科学（物理

几何、卫星、空间）

1.几何：确定地球的形状和大小及确定地面点的几何位置。

2.物理：用物理的方法（重力测量）确定地球形状及其外部重力场。

3.空间：以人造卫星及其它空间探测器为代表的空间大地测量学的理论、技术和方法。

发展的四个阶段：地球圆球、地球椭球、大地水准面、现代大地测量

地球公转：开普勒三大定律 黄道：绕太阳的椭圆轨道

地球自转：地轴：瞬时旋转轴 地轴方向相对于空间的变化（岁差和章动）

岁差：地轴在空间绕黄极发生缓慢的旋转的现象。 周期为26000年。

章动：由于月球引力的影响，导致地轴在岁差的基础上叠加了周期为18.6年的短周期运动，

极移：地轴相对于地球本体内部结构的相对位置变化（国际协议原点CIO）

时间系统：描述卫星或天文现象相应的时间（时空合一）。

组成：一维时间坐标轴＋时间原点为＋时间度量单位

IAG国际大地测量协会IAU国际天文联合会，IUGG地球物理联合会IPMS国际极移服务

BIH国际时间局；CIO国际协议原点；

椭球定位：确定椭球中心的位置 （局部定位：参考椭球 地心定位：总地球椭球）

椭球定向：确定椭球旋转轴的方向。⑴ 椭球短轴平行于地球自转轴

⑵ 大地起始子午面平行于天文起始子午线

二维坐标转换公式推导：平面极坐标公式sincos11ryrx 旋转公式)sin()cos(22ryrx

公式展开写成矩阵形式1122cossinsincosyxyx平移缩放bayxyx1122cossinsincos差计算转化为线性模型进行平sin,cosdc

112112cydxbydycxax

为尺度参数是旋转参数是平移参数和，，ba)1(

反算：byaxcddcyxˆˆˆˆˆˆ12211

三维坐标转换数学模型的解法：

1、用两个公共点计算出两个坐标的距离之比就是尺度参数的近似值。

2、消去平衡移参数，计算R。4个公共点可列如下12个方程。

3、计算平移参数（用第1个点的三个方程）

4、计算旋转参数

5、线性化列误差方程

6、平差计算和精度评定。

重力场：地球上不同位置的点的重力所形成的空间

引力2rmMfF .离心力2mP

重力：地球上任一质点处重力是地球引力与地球自转产生的离心力的合力

引力位：rMfV物理意义：单位质点从无穷远处移动到该点引力所做的功

离心力位：)(2222yxQ 22mP

二阶导数：022222222zyxQ（布阿桑算子）

重力位：QVW )(2222yxrdmfW

二阶算子QV, 外部点 ：QVW 0222222zVyVxVV拉普拉斯方程

引力位是调和函数” 离心位不是调和函数“重力位不是调和函数 正常重力位就是正常椭球上任意一点的重力位。可以由公式精确计算出

扰动位：正常重力位与地球重力位差异。所以地球重力位＝正常重力位＋扰动位

任意一点的扰动位的确定是确定重力位之关键

莫洛金斯基：莫洛金斯基于1954年提出，通过地面上观测的重力值精确求定地面点的扰动位，但其边界不是大地水位面而是地球表面，即直接利用地面上的重力观测值求定地球形状和外部重力场，而不是通过大地水准面求解。

大地水准面高度又称大地水准面差距，似大地水准面高度又称高程异常

正高定义：地面上任一点沿垂线方向到大地水准面的距离。CBCBBdHHH正

⑴不因水准路线而异，是惟一的⑵不能精确测定，正高无法精确求得。

正高不可能得到：从大地水准面沿重力方向到地面点的各占点重力无法测得。

正高改正：用正常重力值代替实际重力 得到的正高称为近似正高 1100OBBBmBmBmBOdhrrrdhrrr是正常重力，r

重力异常：由于地球表面起伏及其内部质量密度分布不均匀，因此正常重力值实测重力值g并不相同，它们的差数g

正常高高差：＝常常ABHHABdh 正常位水准面不平行引起的高差改正重力异常

似大地水准面：由地面点向下量取正常高所得的点形成的连续曲面；

正高：大地水准面 正常高：似大地水准面

垂线偏差：地面上一点的重力向量和相应椭球面上的法线向量之间的夹角定义为该点的垂线偏差

分类：绝对垂线偏差：总地球椭球 相对垂线偏差：参考椭球

测量方法：1、天文大地测量方法：⑴精度高⑵作业量大，只适用于少数的天文大地点上。

2、重力测量方法：重力异常垂线偏差 ⑴假设大地水准面外没有扰动物质，全球重力异常已知; ⑵没有独立应用

3、天文重力方法：天文大地测量方法＋重力测量 说明：地形平坦，基线短精度要求较低

4、还有平面法

大地水准面差距：大地水准面上一点P到正常椭球面的距离， 称为该点的大地水准面差距

扰动位:大地水准面上一点p的实际重力位与正常重力位

之差，称之为该点的扰动位。即：T=W－U

方法：1.地球重力场模型法2.斯托克斯方法3.卫星无线电测高法4.GPS高程拟合法

确定地球形状：天文大地测量方法 重力测量方法 空间大地测量方法

天文大地测量方法，重力测量方法适用于陆地；空间大地测量方法适用于全球最佳是综合

旋转椭球的形状和大小5个基本几何参数：长短半轴 扁率aba

第一偏心率abae22 第二偏心率bbae22'

引入符号：bac2 Bttan Be222cos'

辅助量：BeW22sin1BeV22cos'1 推导2'1eWV 21eVW

222222222222222221cos'11sin''11sin'111sin)1('1sin1eVBeeBeeeBeeeBeeBeW 211eWV2'1eW

推导222e

•

2222222222222222222)(2222)(eeabaabaaabaabaababaabaa>b,V>W,e’>e,c>a21e大小2'1e小大

WVbac1

1. 大地坐标系以椭球面和法线为基础建立的大地纬度B过P点的椭球法线与赤道面的夹角

2. 空间直角坐标系XYZO3. 子午面直角坐标系),,(yxLP4. 地心纬度坐标系),,(LP和归化纬度坐标系),(uLP5. 大地极坐标系P（S，A）

1. 子午面直角坐标系同大地坐标系的关系

BBdxdyocot)90tan(12222byax

BexyyxeyxabByxabdxdytan)1()1(cot222222 1tan)1(2222222bBexax VBbBeWayWBaBeBaxBaBexsinsin)1(cossin1coscos)sin1(22222222

BeNyWaNBNxsin)1(cos2

)1(sin2eNPQBPQy又 2NeQn

2. 空间直角坐标系同子午面直角坐标系

yZLxYLxXsincos

3. 空间直角坐标系同大地坐标系的关系

BeNZLBNYLBNXsin)1(sincoscoscos2点在椭球面上 不在需加高程

法截面：过椭球面上任意一点可作一条垂直于椭球面的法线，包含这条法线的平面叫法截面。 法截线（弧）：法截面与椭球面的交线叫法截线。

1.子午圈曲率半径M

32)1(sin1sinWeaMBdBdxMBdxdSdBdSMBeW22sin1

23VNMVcM bac2

BeVWaN22cos'1

B M 说明

B=0o

0o

B=90o M0=a(1-e2)

a(1-e2)

M90=c M0

M↗ B↗

M90=c

c：极点处（两极）子午圈的曲率半径。

2.卯酉圈曲率半径

卯酉圈：过椭球面上一点的法线，可作无限个法截面。其中一个与该点子午面相垂直的法截面同椭球面相截形成的闭合的圈

VcNWaNWBarxBNrcoscos

3.主曲率半径的计算 子午圈曲率半径M 卯酉圈曲率半径N

4、任意法截线的曲率半径

尤拉公式：NAMARA22sincos1

ARABeN222coscos'1

⑴ 法截线的方位角以子午圈的北方向为基准

⑵ RA与B有关，与A有关。

当A=0o（或180o）时，RA=M（最小值）

当A=90o（或270o）时，RA=N（最大值）

当A：0o→90o，RA：M→N

A：90o→180o，RA：N→M

5.平均曲率半径：过椭球面上一点的一切法截弧（从0到2π），当其数目趋于无穷时，它们的曲率半径的算术平均值的极限，用R表示。MNR

6.M，N，R的关系：N>R>M cMRNooo909090

子午线弧长计算公式与对应椭球对应，不同的椭球有不同的公式。

每秒纬度变化，引起子午线弧长的变化约30米

意义：大地测量正算和反算。高斯反算

⑵子午线上纬度为 ， 间的弧长：12XXX

2.由子午弧长求大地纬度：⑴ 迭代解法⑵直接解法3.平行圈弧长公式

大地线的定义和性质：

定义：椭球面上两点间的最短程曲线叫大地线

⑴ 大地线是一条空间曲线；

⑵ 大地线惟一，位于相对法截线之间。

3.大地线的微分方程和克莱劳方程：描述p到p1时，dS与 dA、dL、dB之间的关系

dSMAdBAdSdBMcoscos dSBNAdLAdSdLBNcossinsincosdSBNAdAtansin B N 说明

B=0o

0o

B=90o N0=a

a

N90=c 卯酉圈即赤道

N↗B↗

卯酉圈子午圈

1B2B