坐标转换之计算公式
一、参心大地坐标与参心空间直角坐标转换
1名词解释:
A:参心空间直角坐标系: a) 以参心0为坐标原点;
b) Z轴与参考椭球的短轴(旋转轴)相重合; c) X轴与起始子午面和赤道的交线重合;
d) Y轴在赤道面上与X轴垂直,构成右手直角坐标系0-XYZ; e) 地面点P的点位用(X,Y,Z)表示; B:参心大地坐标系:
a) 以参考椭球的中心为坐标原点,椭球的短轴与参考椭球旋转轴重合; b) 大地纬度B:以过地面点的椭球法线与椭球赤道面的夹角为大地纬度B; c) 大地经度L:以过地面点的椭球子午面与起始子午面之间的夹角为大地经度L; d) 大地高H:地面点沿椭球法线至椭球面的距离为大地高H; e) 地面点的点位用(B,L,H)表示。 2 参心大地坐标转换为参心空间直角坐标:
X(NH)*cosB*cosLY(NH)*cosB*sinL Z[N*(1e2)H]*sinB公式中,N为椭球面卯酉圈的曲率半径,e为椭球的第一偏心率,a、b 椭球的长短半径,f 椭球扁率,W为第一辅助系数
ea2b2 或 ea22*f1 fW(1e*sin2BaNW
3 参心空间直角坐标转换参心大地坐标
YLarctan()XBarctan(HZ*(NH)(XY)*N*(1e)H222)
X2Y2NcosB二 高斯投影及高斯直角坐标系
1、高斯投影概述
高斯-克吕格投影的条件:1. 是正形投影;2. 子午线不变形
高斯投影的性质:1. 投影后角度不变;2. 长度比与点位有关,与方向无关; 3. 离子午线越远变形越大
为控制投影后的长度变形,采用分带投影的方法。常用3度带或6度带分带,城市或工程控制网坐标可采用不按3度带子午线的任意带。
2、高斯投影正算公式:
xXNNsinBcosBl2sinBcos3B(5t29244)l4224N sinBcos5B(6158t2t4)l6720
NyNcosBlcos3B(1t22)l36Ncos5B(518t2t414258t22)l51203、高斯投影反算公式:
1ylcosBfNf12y1(12t2ff)6Nf2124222y528tf24tf6f8ftfN 120f4
1153t2BBf292t2ffff122Mf4124y 6190tf45tf360NftfyyNfyNf2
表1 BJ54与WGS84基准参数
参考椭球体 长半轴 短半轴 扁率 BJ54基准参数 Krasovsky_1940 6378245 6356863.0188 298.3 WGS84基准参数 WGS 84 6378137 6356752.3142 298.257224
NN232244xXsinBcosBlsimBcosB(5t94)l22244N5246sinBcosB(6158tt)l7206
yNN3223cosBlcosB(1t)l63
N5242225cosB(518tt1458t)l1205(8-42)
8.3.2高斯投影坐标反算公式
BBftfy2tf353t2f22429tyfff2MfNf24MfNftf720M5y6190t246f45tfyfNflyNcosBy36N312t22fffffcosBfy5120N5528t24222f24tf6f8ftffcosBf
(8-56)
1.高斯投影正算公式:
111xXNtm2(5t29244)m4(6158t2t4)m6
24720211yNm(1t22)m3(518t2t4142582t2)m5
6120180m60(13t224)m312(2t2)m5
式中,x,y分别为高斯平面纵坐标与横坐标,为子午线收敛角,单位为度。X为子午线弧长,对于克氏椭球:
X111134.8611B(32005.7799sinB133.9238sin3B0.6976sin5B0.0039sin7B)cosB
对于“IAG 75”椭球:
X111134.0047B(32009.8575sinB133.9602sin3B0.6976sin5B0.0039sin7B)cosB
其余符号为:
ttgB,2e'2cos2B,Nc12,mcosB180l,lLL0
e'a2a2b2,称作第二偏心率;c,称作极曲率半径。L0为
bb2子午线经度。
对于克氏椭球:
e'20.0067385254147,c6399698.90178271对于“IAG 75”椭球:
e'20.0067395018195,c6399596.65198801算出的横坐标y应加上500公里,再在前冠以带号,才是常见的横坐标形式。
2.高斯投影反算公式:
BBf12f2224246 tf90n27.5(53t29t)n0.25(6190t45t)nffffffl12345 180n30(12t2t2ff)n1.5(528f24tf)ncosBftf180n60(1t2f32452 f)n12(25tf3tf)n 式中,Bf为底点纬度,以度为单位。n公式,只是以底点纬度代替大地纬度。
y12fc,其余符号同正算