子午线收敛角的计算公式与计算精度分析
子午线收敛角的计算公式与计算精度分析
应红立
(中国矿业大学环境与测绘学院)
摘要:介绍了子午线收敛角的概念和计算方法,采用大地坐标五次项展开式计算值近似地作摘要:
为子午线收敛角的真值,分析了随纬度纬度、经差的变化规律,分析了在经差3.5°,纬度5°-85°时近似公式γ=l *sinB 计算值与真值的差值,给出了纬度与差值三次四项式的拟合公式,同时采用高斯平面直角坐标五次展开式与其简化公式γ=ky进行分析比较,发现在精度要求不高的情况下,简化公式完全符合精度要求。
关键词:子午线收敛角经差高斯投影坐标方位角计算公式分析
0引言
在大地控制测量以及高精度工程测量中常常来进行大地方位角和平面坐标方位角两类方位角之间的转换,如利用陀螺经纬仪所测的方位角为大地方位角,当我们需要获得坐标方位角时,就势必要先计算子午线收敛角。
子午线收敛角的计算精度会直接影响坐标方位角的精度,目前,关于子午线收敛角计算公式在许多大地测量书籍中[]都已经给出,但由于精度要求不同以及公式推导的方法不同,使得计算公式不尽相同,但它们的计算精度如何有关文献[2-4]中涉及较少。本文描述子午线收敛角的概念,分析各种计算方法及其精度。
系:
A=α+γ(1)
γ即为子午线收敛角,它是某点的真子午线方向与该点的坐标北方向之间的夹角。γ有“+,-”之分,点位于中央子午线的东侧,取正号;反之,取负号。由公式(1)
图1.子午线收敛角概念图下式是真子午线与坐标纵线的角度关变形,子午线投影后向中央子午线弯曲,并向两极收敛。
如图1所示,地球椭球体面上一点处子午线收敛角是该点真子午线与位于此点所在的投影带的中央子午线之间的夹角,即在高斯平面上的真子午线与坐标纵线的夹角,
通常用γ表示。
1子午线收敛角的概念
地球的形状接近于旋转椭球面,椭球面是一个不可直接展开的曲面,大地坐标经过高斯投影得到平面直角坐标后,必然有投影
可知,只要求得α和γ,一切问题就迎刃而解,其中坐标方位角的计算方法比较简单,在此就不详加叙述了,接下来介绍几种γ的计算方法。
2子午线收敛角的计算公式
2.1用大地坐标计算
2.1.1文献[2]中给出的展开五次项公式
γ=l *sinB[1+l 2
cos 2B (1+3η2
3
4
(2)+2η4) +l
15
cos 4B (2−t 2)]
式中,l 为该点距离中央子午线的经差
l =L −L 0,B 为该点的纬度,η为垂线偏差u 的卯酉圈分量(子午圈分量为ξ)
,
t =M
N
tan A
,A 为大地方位角。2.1.2文献[3]中给出的近似公式γ=l *sinB (3)
式中,l 为经差,即某点与中央子午线的经度差,B 为点位所处的纬度,γ值为子午线收敛角。
2.2用高斯平面坐标计算
2.2.1.文献[2]中给出的展开五次项公式
γ=y N t y 3(1+t 22
f −3t f f −ηf −f 3N f
y 2(4)
2η424f ) +15N 5t f
(2+5t f +3t f ) f
公式(4)中,N f 、t f 、ηf 均为垂足纬度(底点纬度)B f 的函数,将B f 的数值代入相应公式后可求得。
2.2.2文献[4]中给出的公式γ=ky(5)
式中γ=ky———以分为单位;y ———点的横坐标,以km 为单位;k ———系数,以纵坐标x (以km 计)
为引数由表一中查取。
3子午线收敛角的计算精度分析
3.1公式分析[2]由公式(2)分析可知
(1)γ为l 的奇函数,而且l 愈大,γ也愈大;
(2)γ有正负,点位于中央子午线以东时,γ为正;在西时,为负;
(3)当l 不变时,γ随纬度增加而增大。3.2计算精度分析
3.2.1采用大地坐标时精度分析:此处选择克拉索夫斯基椭球参数计算:当B 取5°—85°时由于实际工作涉及到换带计算,所以取l =3.5°来计算子午线收敛角差值△;当l 取3.5°时,γ有不同值,如图2所示:
B (°)
γ(2)(″)
51099.547281098.162341.[1**********].163089652.159865.00322
102190.662282187.967002.[1**********]5.5675210321.315644.25188
153264.981483261.119913.[1**********]5.3299610911.919983.40998
204314.27784309.453744.[1**********]2.0234811419.478012.54547
255330.52725324.990015.[1**********]1.8529611840.126931.72603
306305.972046299.999905.[1**********]1.6802812170.665341.01494
357233.181567227.062996.[1**********]9.0430412408.577640.46540
408105.109128099.123775.[1**********]2.1718412552.053170.11867
458915
.144768909.545325.59944
γ(3)(″)
△(″)
B (°)
γ(2)(″)
γ(3)(″)
△(″)
图2.B 与△的数值表
由上表数据分析可得如下结论:(1)用近似公式得到的子午线收敛角与真值之间的差值△随着纬度的不同而不同;差值△随着纬度的增加先增加再减少,在B 取35°时取得最大值,6.11857″;
(2)可以看出在同一子午线上,差值△在中纬(20°-55°)之间误差较大;
(3)在中低纬度(B=5°-80°)时,公式(3)的计算精度能达到秒级。
绘制纬度B 和差值△的函数图象,如图3:
图3.B 与△的函数图像为便于实际应用,同时考虑我国所处的经纬度,我们选择四次三项式来拟合经差为3.5°时的纬度B 和子午线收敛角差值△间的关系[1],得到拟合多项式公式如下(y 代表B ;x 代表子午线收敛角差值△):
由于实际测量工作中,大多数情况下得到的只是高斯平面直角坐标,如果再化为大地坐标计算子午线收敛角,过程较为复杂,
3.2.2采用高斯平面直角坐标时精度分析:
实际上公式(4)是由公式(2)计算得来的,把公式(2)中的l 换成高斯平面直角坐标,把B 换成B f 推导可得公式(4)
子午线收敛角系数K 值表
x/kmK
1000.
00852000.
01703000.
02554000.
03415000.
04266000.
05127000.
05988000.
06849000.
077110000.
085811000.
094512000.
103313000.
112114000.
121015000.
1300
△[***********][1**********]0
x/kmK 16000.
139017000.
148118000.
157319000.
166620000.
175921000.
185422000.
194923000.
204624000.
214325000.
224226000.
234227000.
244428000.
254729000.
265130000.
2753
△[***********][***********]7
x/kmK 31000.
286532000.
297533000.
308634000.
320035000.
331636000.
343437000.
355438000.
367739000.
380240000.
393141000.
406242000.
419643000.
433444000.
447545000.
4619
△[***********][***********]141144149
x/kmK 46000.
476847000.
492148000.
507849000.
524050000.
540751000.
557952000.
575753000.
594154000.
613155000.
632856000.
653357000.
674558000.
696759000.
719760000.
7437
△[***********][***********]230240248
所以通常采用公式(5)来进行简化计算,关于公式(5)有一张子午线收敛角系数K 值表,如下表,实际工作中可用它来进行简化计算[4]。
本文分别取X 值如下(单位:km )
X 100500
[***********][***********]55006000
当经差为3.5°时,在赤道上Y 值最大为390km ,所以本文取Y 值为390km 来计算。
X (km )
γ(4)(″)γ(5)(″)
通过编程计算,公式(4)和(5)计算结果如下:15003041.9723042-0.[1**********]0.36212652.38
-2.018
20004116.3724116.060.[1**********]7.20814807.52
-0.312
25005246.6785246.28
30006452
.5856442.02
100199.086198.90.[1**********].5957759.44
-0.845
5001000
997.2007.378058996.842007.
720.[1**********].2579198.54
-2.283
-0.[1**********]7.56610808.46
-0.894
△(″)X (km )
γ(4)(″)γ(5)(″)
△(″)
0.39810.[1**********]1.36817402.58
-1.212
分析其子午线收敛角差值△可知:(1)在Y 值一定的情况下,子午线收敛角差值△随X 值的变化而变化,但变化呈现无规律状态,在X=3000,4000,5000km ,误差达到秒级,X 取其他值时误
差达到0.1″级;
(2)利用高斯平面直角坐标计算出的子午线收敛角与真值之间的差值较小,计算精度可达到0.1″,准确度可达到1″。X 值与△函数图像如图4:
图4.X 值与△函数图像
因为子午线收敛角差值△呈现无规律状态,所以不便于拟合其函数式,此处就不再做拟合工作。
着经差和纬度的不同而变化,对比分析计算数据可得到如下结论[1-2]:
(1)同一平行圈上,经差越大,差值愈大;(2)在同一子午线,点位处于中纬度(25°-35°)时较差△较大;
(3)在中低纬度(B=5°-80°)时,公式(3)
4结束语
通过以上计算可知,利用公式(2)和公式(3)计算的子午线收敛角存在差值,并随
的计算精度能达到秒级。
(4)采用公式(5)计算时,在Y 值一定的情况下,子午线收敛角差值△随X 值的呈现无规律变化,在X=3000,4000,5000km ,误差达到秒级,X 取其他值时误差达到0.1″级。
(5)在实际工作中,大多数情况下都是利用高斯平面坐标计算子午线收敛角,公式(5)不仅计算方便而且具有足够的计算精度,所以公式(5)具有较大的实际意义和作用,应在实际工作中大力推广。
参考文献:
[1]邬熙娟,江国焰,高俊强.子午线收敛角计算公式及计算精度分析[J].现代测绘,2005,28(6):22-25.
[2]张华海等,等.应用大地测量学基础[M].徐州:中国矿业大学出版社,2011.8.
[3]高井祥,等.数字测图原理与方法[M].徐州:中国矿业大学出版社,2010.8.
[4]张国良,等.矿山测量学[M].徐州:中国矿业大学出版社,2008.10.
.