2. 32022部队, 湖北 武汉 430033
2. 32022 Troops, Wuhan 430033, China
地球半径是测量和地球科学计算中最常用的基本参数,根据地球科学、空间科学、导航定位的要求和某些需求,常用到平均曲率半径、平均球半径、等距离球半径、等面积球半径和等体积球半径这5种地球半径[1-4]。随着空间技术和计算机技术在大地测量及地图学中的应用和发展,研究常用地球半径间的关系及它们之间的差异问题具有更加重要的实用价值。对于这一问题,国内外许多学者进行了研究,取得了显著成果。文献[5—6]提出使用不同的球体半径表示球体特性,将不同球体半径代入克拉索夫斯基椭圆参数求得数值解,再代入球体面积公式,求得不同半径下大圆航线长度和角度变形大小,发现使用等角球半径计算大圆航线可以满足航海使用精度要求。文献[7]将子午圈曲率半径和卯酉圈曲率半径展为大地纬度的幂级数,并且借助辅助函数V导出了不同形式的平均曲率半径展开式。文献[8]中利用等距离球半径和等面积球半径,推导出常用纬度函数应用于求解地图投影变换,提出了变系数线性插值方法,解决了测量和地图学中所要求的计算精度问题。文献[9]利用球体半径和计算机代数系统对大地测量学和地理信息系统中的问题进行了研究。
从目前来看,前人对这一领域作了很多卓有成效的工作,但主要是针对其中一种或几种地球半径进行计算和应用,多是对等角纬度、等面积纬度及等距离纬度的计算与分析[10-15],很少有文献将这几种常用的地球半径间的差异在符号和数值上进行系统的比较。计算机代数系统可以将基本的数学公式展开成幂级数形式,推导出的公式相对于人工计算具有更高的精确度[16-19]。为丰富对这一问题的研究,使人们对这几种常用地球半径形成较直观的认识,本文着重研究了常用地球半径间的差异问题,推导出了常用地球半径间差异符号表达式,最后以CGCS2000参考椭球[20]为例,对常用地球半径间的差异进行了数值分析和比对,发现平均球半径和等体积球半径在等面积地图投影中可以对等面积球半径进行替换。
1 常用地球半径定义平均曲率半径Ra常用于制作地球表面上局部地区地图[21],也可以取制图区域中心点,记椭球偏心率为e,则它关于大地纬度B的表达式
式中,N为卯酉圈曲率半径;M为子午圈曲率半径,将其展开成偏心率e的幂级数形式
式(2)中的系数计算如下
就地球总体而言,平均球半径(取地球椭球三轴半轴长的算术平均值,用于简单决定球半径)[22]、等面积球半径(保持球体表面积等于地球椭球面相应全面积而决定的球半径,主要用于等面积投影)[23]、等距离球半径(使球面经线总长等于地球椭球面经线总长而决定的球半径,主要用于等距离投影中)[24]和等体积球半径(使地球球体的体积等于地球椭球体的体积来决定的球半径)的表达式与大地纬度B无关,只取决于地球椭球模型的参数a、e,其表达式为
式中, Re表示平均球半径;RF表示等面积球半径;RS表示等距离球半径;RV表示等体积半径。
2 常用地球半径间的差异符号表达式由于5种常用地球半径具有一定差异,在实际应用中会用到它们的差值表达式。为了解常用地球半径间的差异情况,特别是差异最值点及其对应差异最值,可基于它们之间的差值表达式推导出其差异符号表达式。
2.1 平均曲率半径与4种常用球体半径间的差异符号表达式通过对平均曲率半径与4种常用球体半径的表达式进行分析,即可推算出其差异符号表达式。以平均球半径为例
将式(5)展开成幂级数的形式
式(6)中系数计算公式为
从上式中可以看出,在确定地球椭球体参考模型之后可以确定常数值e,此时两者差值只与大地纬度B有关,由于大地纬度B的取值范围为
当B取0时,存在最小值, 将其最小值进一步展开成偏心率e的幂级数形式
类似的,将平均曲率半径与等面积地球半径、等距离地球半径和等体积地球半径分别进行作差,分析结果同平均曲率半径与平均球半径差值分析结果:当B取
当平均曲率半径与平均球半径相等时,有如下公式
对式(10)进行整理,略去推导过程,可得
式(11)中,B0e为平均曲率半径与平均球半径相等时的大地纬度,将其展开为偏心率e的幂级数形式
故采用与平均球半径类似的方法,推导出它们与平均曲率半径相等时大地纬度B的符号表达式,最后结果列于表 3。
2.2 4种常用球体半径间的差异符号表达式
为系统地比较各常用球体半径之间的差异,除对平均曲率半径与4种常用球体半径间的差异进行分析以外,需要对其余4种常用球体半径间的差异进行分析。以下分别对平均球半径、等面积球半径、等距离球半径和等体积球半径间的差异进行分析。就地球总体而言,它们的表达式与大地纬度B无关,只取决于地球椭球模型参数a、e,与推导平均曲率半径和4种常用球体半径间差异值符号表达式类似,通过作差之后级数展开的形式进行分析(以平均球半径和等体积球半径间差异值符号表达式为例),其差异可以表示为
将式(13)展开成幂级数的形式
相对于平均曲率半径与其余4种常用球体半径间差异,它们相互之间的差异就显得相对简单。类似的,可以推导出平均球半径与等面积球半径、平均球半径与等距离球半径、等距离球半径与等面积球半径、等距离球半径与等体积球半径和等面积球半径与等体积球半径间的差异符号表达式,最后结果列于表 4。
由表 4可以看出,4种常用球体半径间的差异只与e有关,它们之间差异值随着e的增大而增大,且等距离球半径和等体积球半径之间的差异最小,平均球半径和等距离球半径之间的差异最大。
3 算例分析为使人们对各地球半径间的差异在数值上有一个直观的认识,下面以CGCS 2000参考椭球(a=6 378 137, e=0.081 819 191 042 8)为例,对5种常用地球半径间的差异进行数值比较与分析。
3.1 平均曲率半径与4种常用球体半径间的差异比较为了解平均曲率半径与平均球半径、等面积球半径、等距离球半径和等体积球半径间的差异,可以绘制大地纬度B∈[0°, 90°],平均曲率半径与各常用球体半径的差异曲线图,如图 1所示;表 5为大地纬度B∈[0°, 90°],每隔15°对应的平均曲率半径与各常用球体半径间的差异。
m | |||||||
B | 0° | 15° | 30° | 45° | 60° | 75° | 90° |
Ra-Re | -14 256.5 | -11 404.6 | -3 599.9 | 7 092.3 | 17 820.5 | 25 697.0 | 28 584.9 |
Ra-RF | -14 254.9 | -11 403.0 | -3 598.4 | 7 093.9 | 17 822.1 | 25 698.6 | 28 586.4 |
Ra-RS | -10 696.8 | -7 844.9 | -4 036.8 | 10 651.9 | 21 380.1 | 29 256.6 | 32 144.5 |
Ra-RV | -14 248.5 | -11 396.6 | -3 592.0 | 7 100.2 | 17 828.5 | 25 705.0 | 28 592.8 |
如图 1、表 5所示,平均曲率半径与其余4种常用球体半径之间差异值随着大地纬度B∈[0°, 90°]增大而增大,平均曲率半径与平均球半径、等面积球半径和等体积球半径差异值在所选取的7个点上差异很小,在图 1上表现为3条曲线几乎重合成一条曲线。为了便于观察3条曲线之间的差异,取大地纬度B∈[0°, 1°]范围内,画出3条曲线,如图 2所示。
由图 2可以看出,在1°的变化范围内,平均曲率半径与平均球半径、等面积球半径和等体积球半径间差异的变化图,3条曲线的差异值在10 m以内,可以看出平均球半径、等面积球半径和等体积球半径间相差较小,因此,在实际应用中可以根据需要选择合适的半径以避免公式计算和推导过程中的复杂性。
根据已推导出的平均曲率半径与其余4种常用球体半径间的差异符号表达式,可以计算出在大地纬度B∈[0°, 90°]范围内存在一个点使得两个半径相等,即差异值为0。
由图 1和表 6可以看出,差异曲线与坐标轴的交点即平均曲率半径与平均球半径、等面积球半径、等距离球半径和等体积球半径相等的点,具体数值总结于表 6。其中,等距离球半径和平均曲率半径相等点的大地纬度在30°3′45.04″,与其余3个相等点大地纬度相差5°左右,与图 1反映的差异曲线相符。
由于平均曲率半径是随着大地纬度B逐点变化的,是微观上的,而4种常用球体半径是宏观意义上地球球体的平均值,它们之间定义不同,差异逐点地变化,运用的场合也不同,因此宏观无法代替微观上地平均曲率半径,但在一定程度上,平均曲率半径与4种常用球体半径在某一个特殊点上相等,此时它们之间可以互相替代。
3.2 4种常用球体半径间的差异比较为全面分析常用地球半径间的差异问题,现以CGCS2000参考椭球(a=6 378 137, e =0.081 819 191 0)为例,对4种常用球体半径间的差异进行数值比较与分析,最后结果列于表 7中。
由表 7可以看出,平均球半径和等面积球半径之间的差异值的绝对值最小,平均球半径和等距离球半径之间的差异值绝对值最大。由于4种球体半径是不同意义上定义的球体半径平均值,因此,它们之间差异很大。等距离球半径与其余3种地球半径差异较大,因此在进行等距离投影时不宜使用其余地球半径进行替换;反观平均球半径和等体积球半径与等面积球半径之间,在CGCS2000参考系下数值相差精度在8 m以内,将平均球半径、等面积球半径和等体积球半径分别代入等面积投影长度比公式中,可以看出,椭球面在球面上的等面积投影的长度变形在赤道处最大,其值分别为±0.111 7%、±0.111 8%和±0.111 9%,由此可见,这3种地球半径对等面积投影的影响很小。由于等面积球半径公式较复杂,对于复杂公式计算推导存在一定的困难,因此,在实际应用过程中可以用平均球半径和等体积球半径来替代等面积球半径。
4 结论本文对5种常用地球半径间的差异进行了研究,推导出了常用地球半径间的差异符号表达式,并以CGCS2000椭球为例对它们进行了数值分析和比对,得出如下结论。
(1) 平均曲率半径与4种常用球体半径差异值随大地纬度的增大而增大,差异绝对值随大地纬度的增大先减小后增大,在某一个特定点处差异值为0。其中,当B取
(2) 平均曲率半径是微观上的平均值,是随着大地纬度B逐点变化的;4种常用球体半径是宏观意义上的地球球体上的平均值,与大地纬度B无关。平均曲率半径与4种常用球体半径除在特殊点等价外,其余点无法进行相互代替应用。
(3) 4种常用球体半径间,平均球半径和等面积半径间的差异绝对值最小,绝对值最小符号表达式的首项为
(4) 将常用地球半径间差异值表示为符号形式,并统一展开为偏心率e的幂级数形式,该表达式易于比较分析,一定程度上丰富了测量及地图学数学分析理论。
[1] |
吕志平, 张建军, 乔书波.
大地测量学基础[M]. 北京: 解放军出版社, 2005.
LÜ Zhiping, ZHANG Jianjun, QIAO Shubo. Basis of geodesy[M]. Beijing: Publishing House of Liberation Army, 2005. |
[2] |
华棠.
海图数学基础[M]. 北京: 海潮出版社, 1985.
HUA Tang. The mathematical basis of chart[M]. Beijing: Haichao Press, 1985. |
[3] |
王瑞, 李厚朴.
基于地球椭球模型的符号形式的航迹计算法[J]. 测绘学报, 2010, 39(2): 151–155.
WANG Rui, LI Houpu. Symbolic expressions of the track calculation methods based on the ellipsoidal model[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(2): 151–155. |
[4] |
熊介.
椭球大地测量学[M]. 北京: 解放军出版社, 1988: 18-23.
XIONG Jie. Ellipsoidal geodesy[M]. Beijing: Publishing House of Liberation Army, 1988: 18-23. |
[5] |
丁佳波.
拟大地线航法的球体半径[J]. 中国航海, 1988(1): 10–14.
DING Jiabo. The spherical radius of the approximate geodetic line nautical method[J]. Navigation of China, 1988(1): 10–14. |
[6] |
丁佳波.
拟大地线航法地理坐标的解算公式[J]. 天津航海, 2003(1): 4–6.
DING Jiabo. Solving formula for the geographical coordinate of the geodesic line sailing[J]. Tianjin Navigation Aditorial Department, 2003(1): 4–6. DOI:10.3969/j.issn.1005-9660.2003.01.002 |
[7] |
熊介.
椭球面曲率半径展为纬差的幂级数[J]. 解放军测绘学院学报, 1985(2): 15–21.
XIONG Jie. Expansions of radius of curvature of the ellipsoid in power series of latitude difference[J]. Journal of the PLA Institute of Surveying and Mapping, 1985(2): 15–21. |
[8] |
杨启和, 杨晓梅.
测量和地图学中应用的三种纬度函数及其反解变换的线性插值方法[J]. 测绘学报, 1997, 26(1): 92–94.
YANG Qihe, YANG Xiaomei. Three types of latitudinal function and linear interpolation algorithms used in surveying and cartography[J]. Acta Geodaetica et Cartographica Sinica, 1997, 26(1): 92–94. |
[9] | AWANGE J L, GRAFAREND E W. Solving algebraic computational problems in geodesy and geoinformatics:the answer to modern challenges[M]. Berlin: Springer, 2005. |
[10] | ADAMS O S. Latitude developments connected with geodesy and cartography with tables, including a table for lambert equal-area meridional projection[M]. Washington, DC: Government Printing Office, 1921. |
[11] | SNYDER J P. Map projections-a working manual[M]. Washington, DC: Government Printing Office, 1987. |
[12] |
杨启和.
地图投影变换原理与方法[M]. 北京: 解放军出版社, 1989.
YANG Qihe. The theory and method of map projection[M]. Beijing: PLA Press, 1989. |
[13] | YANG Qihe, SNYDER J P, TOBLER W R. Map projection transformation:principles and applications[M]. London: Taylor & Francis, 2000. |
[14] |
吴忠性, 杨启和.
数学制图学原理[M]. 北京: 测绘出版社, 1989.
WU Zhongxing, YANG Qihe. Principle of mathematical cartography[M]. Beijing: Surveying and Mapping Press, 1989. |
[15] |
李忠美, 李厚朴, 边少锋.
常用纬度差异极值符号表达式[J]. 测绘学报, 2014, 43(2): 214–220.
LI Zhongmei, LI Houpu, BIAN Shaofeng. Symbolic expressions of difference extrema between regular latitudes[J]. Acta Geodaetica et Cartographica Sinica, 2014, 43(2): 214–220. DOI:10.13485/j.cnki.11-2089.2014.0031 |
[16] |
边少锋, 李厚朴, 李忠美.
地图投影计算机代数分析研究进展[J]. 测绘学报, 2017, 46(10): 1557–1569.
BIAN Shaofeng, LI Houpu, LI Zhongmei. Research progress in mathematical analysis of map projection by computer algebra[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(10): 1557–1569. DOI:10.11947/j.AGCS.2017.20170396 |
[17] |
李厚朴, 边少锋, 钟斌.
地理坐标系计算机代数精密分析理论[M]. 北京: 国防工业出版社, 2015.
LI Houpu, BIAN Shaofeng, ZHONG Bin. Geographic coordinate system computer algebraic precision analysis theory[M]. Beijing: National Defense Industry Press, 2015. |
[18] |
边少锋, 许江宁.
计算机代数系统与大地测量数学分析[M]. 北京: 国防工业出版社, 2004.
BIAN Shaofeng, XU Jiangning. Computer algebra system and mathematical analysis in geodesy[M]. Beijing: National Defend Industry Press, 2004. |
[19] |
边少锋, 陈良友, 纪兵.
Mathematica软件在大地测量中的应用[J]. 测绘科学与工程, 2004, 24(2): 25–28.
BIAN Shaofeng, CHEN Liangyou, JI Bing. Applications of mathematica system in geodesy[J]. Geomatic Science and Engineering, 2004, 24(2): 25–28. |
[20] |
陈俊勇.
中国现代大地基准——中国大地坐标系统2000(CGCS 2000)及其框架[J]. 测绘学报, 2008, 37(3): 269–271.
CHEN Junyong. Chinese modern geodetic datum-Chinese Geodetic Coordinate System 2000(CGCS 2000) and its frame[J]. Acta Geodaetica et Cartographica Sinica, 2008, 37(3): 269–271. DOI:10.3321/j.issn:1001-1595.2008.03.001 |
[21] |
陈武军.
地球曲率半径取值精度对控制测量精度的影响分析[J]. 地矿测绘, 2012, 28(4): 27–28, 31.
CHEN Wujun. Analysis of the influence earth curvature radius value precision to control survey precision[J]. Surveying and Mapping of Geology and Mineral Resources, 2012, 28(4): 27–28, 31. DOI:10.3969/j.issn.1007-9394.2012.04.009 |
[22] |
陈丽华, 汪孔政.
关于参考椭球平均半径的探讨[J]. 测绘通报, 2000(10): 15–17.
CHEN Lihua, WANG Kongzheng. On mean radius of reference Ellipsoid[J]. Bulletin of Surveying and Mapping, 2000(10): 15–17. DOI:10.3969/j.issn.0494-0911.2000.10.006 |
[23] |
李厚朴, 边少锋, 陈良友.
等面积纬度函数和等量纬度变换的直接解算公式[J]. 武汉大学学报(信息科学版), 2011, 36(7): 843–846.
LI Houpu, BIAN Shaofeng, CHEN Liangyou. The direct calculating formulae for transformations between authalic latitude function and isometric latitude[J]. Geomatics and Information Science of Wuhan University, 2011, 36(7): 843–846. |
[24] |
陈成, 金立新, 李厚朴, 等.
等距离球面高斯投影[J]. 测绘通报, 2017(10): 1–6.
CHEN Cheng, JIN Lixin, LI Houpu, et al. A kind of gauss projection based on the equidistant sphere[J]. Bulletin of Surveying and Mapping, 2017(10): 1–6. DOI:10.13474/j.cnki.11-2246.2017.0306 |