椭球面上的常用坐标系及其相互关系
§6.2 椭球面上的常用坐标系及其相互关系
6.2.1大地坐标系 P点的子午面NPS与起始子午面NGS所构成的二面角
叫做P点的大地经度,由起始子午面起算,向东为正,L,
o
叫东经(0°~180°),向西为负,叫西经(0~180°)。P点的法线Pn与赤道面的夹角B,叫做P点的大地纬度。由赤道面起算,向北为正,叫北纬(0°~90°);向南为负,叫南纬(0°~90°)。
大地坐标系是用大地经度L、大地纬度B和大地高H表示地面点位的。过地面点P的子午面与起始子午面间的夹角叫P点的大地经度。由起始子午面起算,向东为正,叫东经(0°~180°),向西为负,叫西经(0°~-180°)。过P点的椭球法线与赤道面的夹角叫P点的大地纬度。由赤道面起算,向北为正,叫北纬(0°~90°),向南为负,叫南纬(0°~-90°)。从地面点P沿椭球法线到椭球面的距离叫大地高。大地坐标坐标系中,P点的位置用L,B表示。如果点不在椭球面上,表示点的位置除L,B外,还要附加另一参数——大地高H,它同正常高H正常及正高H正有如下关系 HH正常(高程异常)
HH正N(大地水准面差距)
6.2.2空间直角坐标系
以椭球体中心O为原点,起始子午面与赤道面交线为X轴,在赤道面上与X轴正交的方向为Y轴,椭球体的旋转轴为Z轴,构成右手坐标系O-XYZ,在该坐标系中,P点的位置用X,Y,Z表示。
地球空间直角坐标系的坐标原点位于地球质心(地心坐标系)或参考椭球中心(参心坐标系),z轴指向地球北极,x轴指向起始子午面与地球赤道的交点,y轴垂直于XOZ面并构成右手坐标系。
6.2.3子午面直角坐标系
设P点的大地经度为L,在过P点的子午面上,以子午圈椭圆中心为原点,建立x,y平面直角坐标系。在该坐标系中,P点的位置用L,x,y表示。
6.2.4大地极坐标系
M为椭球体面上任意一点,MN为过M点的子午线,S为连结MP的大地线长,A为大地线在M点的方位角。以M为极点,MN为极轴,S为极半径,A为极角,这样就构成大地极坐标系。在该坐标系中P点的位置用S,A表示。
椭球面上点的极坐标(S,A)与大地坐标(L,B)可以互相换算,这种换算叫做大地主题解算。
6.2.5各坐标系间的关系
椭球面上的点位可在各种坐标系中表示,由于所用坐标系不同,表现出来的坐标值也不同。
1.子午面直角坐标系同大地坐标系的关系
过P点作法线Pn,它与x轴之夹角为B,过P点作子午圈的切线TP,它与x轴的夹角为(90°+B)。子午面直角坐标x,y同大地纬度B的关系式如下:
x
acosBe2sin2B
acosB
W
y
a(1e2)sinBe2sin2B
absinB(1e2)sinB WV
2.空间直角坐标系同子午面直角坐标系的关系
空间直角坐标系中的P2P相当于子午平面直角坐标系中的y,前者的OP2相当于后者的x,并且二者的经度L相同。
XxcosL
YxsinL
Zy
3.空间直角坐标系同大地坐标系的关系
同一地面点在地球空间直角坐标系中的坐标和在大地坐标系中的坐标可用如下两组公式转换。
xNHcosBcosL
yNHcosBsinL
zN1e2HsinB
y
L
x
2
zNesinB
B 22
xy
z2
HN1e
sinB
式中:e——子午椭圆第一偏心率,可由长短半径按式e2a2b2/a2算得。
N——法线长度,可由式Na/
e2sin2B算得。
§6.3 几种主要的椭球公式
过椭球面上任意一点可作一条垂直于椭球面的法线,包含这条法线的平面叫做法截面,法截面同椭球面交线叫法截线(或法截弧)。包含椭球面一点的法线,可作无数多个法截面,相应有无数多个法截线。椭球面上的法截线曲率半径不同于球面上的法截线曲率半径都等于圆球的半径,而是不同方向的法截弧的曲率半径都不相同。
6.3.1子午圈曲率半径
子午椭圆的一部分上取一微分弧长DKds,相应地有坐标增量dx,点n是微分弧dS的曲率中心,于是线段Dn及Kn便是子午圈曲率半径M。
任意平面曲线的曲率半径的定义公式为:
M
dS dB
子午圈曲率半径公式为:
MM
a(1e2)W3
cN 或 M32VV
M与纬度B有关.它随B的增大而增大,变化规律如下表所示:
6.3.2卯酉圈曲率半径
过椭球面上一点的法线,可作无限个法截面,其中一个与该点子午面相垂直的法截面同椭球面相截形成的闭合的圈称为卯酉圈。在图中PEE即为过P点的卯酉圈。卯酉圈的曲率半径用N表示。
为了推导N的表达计算式,过P点作以O为中心的平行圈PHK的切线PT,该切线位于垂直于子午面的平行圈平面内。因卯酉圈也垂直于子午面,故PT也是卯酉圈在P点处的切线。即
所以PT是平行圈PHK及卯酉圈PT垂直于Pn。
PEE在P点处的公切线。
卯酉圈曲率半径可用下列两式表示:
NN
a Wc V
6.3.3 任意法截弧的曲率半径
子午法截弧是南北方向,其方位角为0°或180°。卯酉法截弧是东西方向,其方位角为90°或270°。现在来讨论方位角为A的任意法截弧的曲率半径RA的计算公式。 任意方向A的法截弧的曲率半径的计算公式如下:
RA
NN
12cos2A1e2cos2Bcos2A
(7-87)
6.3.4 平均曲率半径
在实际际工程应用中,根据测量工作的精度要求,在一定范围内,把椭球面当成具有适当半径的球面。取过地面某点的所有方向RA的平均值来作为这个球体的半径是合适的。这个球面的半径——平均曲率半径R:
RMN
或
R
bcNa222
VWWV
(1e2)
因此,椭球面上任意一点的平均曲率半径R等于该点子午圈曲率半径M和卯酉圈曲率
半径N的几何平均值。 6.3.5 子午线弧长计算公式
子午椭圆的一半,它的端点与极点相重合;而赤道又把子午线分成对称的两部分。 如图所示,取子午线上某微分弧PPdx,令P点纬度为B,P点纬度为BdB,P点的子午圈曲率半径为M,于是有: dxMdB
从赤道开始到任意纬度B的平行圈之间的弧长可由下列积分求出:
XMdB
0B
式中M可用下式表达:
M
a0a2cos2Ba4cos4Ba6cos6Ba8cos8B
a0m0a2
其中: a4
a6a8
m23535
m4m6m82816128m2m4157m6m8
223216
m437m6m8 81632m6m83216m8
128
经积分,进行整理后得子午线弧长计算式:
aaaa
Xa0B2sin2B4sin4B6sin6B8sin8B
2468
为求子午线上两个纬度B1及B2间的弧长,只需按上式分别算出相应的X1及X2,而后取差:XX2X1,该X即为所求的弧长。
克拉索夫斯基椭球子午线弧长计算公式:
X111134.861B16036.480sin2B16.828sin4B0.022sin6B
X111134.861B32005.780sinBcosB133.929sin3BcosB0.697sin5BcosB
1975年国际椭球子午线弧长计算公式:
X111133.005B16038.528sin2B16.833sin4B0.022sin6B
X111133.005B32009.858sinBcosB133.960sin3BcosB0.698sin5BcosB
6.3.6 底点纬度计算
在高斯投影反算时,已知高斯平面直角坐标(X,Y)反求其大地坐标(L,B)。首先X当作中央子午线上弧长,反求其纬度,此时的纬度称为底点纬度或垂直纬度。计算底点纬度的公式可以采用迭代解法和直接解法。 (1)迭代法
在克拉索夫斯基椭球上计算时,迭代开始时设
B1.8611 fX/111134
以后每次迭代按下式计算:
Bif1(XF(Bif))/111134.8611
F(Bif)16036.4803sin2Bif16.8281sin4Bif0.0220sin6Bif
重复迭代直至Bif1Bif为止。
在1975年国际椭球上计算时,也有类似公式。
(2)直接解法 1975年国际椭球:
X/6367452.133
BfB{50228976[293697(238322cos2)cos2]cos2}1010sincos
克拉索夫斯基椭球: X/6367588.4969
Bf{50221746[293622(235022cos2)cos2]cos2}
6.3.7 大地线
椭球面上两点间的最短程曲线叫做大地线。在微分几何中,大地线(又称测地线)另有这样的定义:“大地线上每点的密切面(无限接近的三个点构成的平面)都包含该点的曲面法线”,亦即“大地线上各点的主法线与该点的曲面法线重合”。因 假如在椭球模型表面A,B两点之间,画出相对法截线如图所示,然后在A,B两点上各插定一个大头针,并紧贴着椭球面在大头针中间拉紧一条细橡皮筋,并设橡皮筋和椭球面之间没有摩擦力,则橡皮筋形成一条曲线,恰好位于相对法截线之间,这就是一条大地线。由于橡皮筋处于拉力之下,所以它实际上是两点间的最短线。
在椭球面上进行测量计算时,应当以两点间的大地线为依据。在地面上测得的方向、距离等,应当归算成相应大地线的方向、距离。