-
地球向径和平均曲率半径是测量和地球科学计算中常用的基本参数,根据地球科学、空间科学的要求和某些需求,常用到地球向径、平均曲率半径、平均球半径、等距离球半径、等面积球半径和等体积球半径[1-7]。由于地球是一个旋转椭球,向径与曲率半径是背离的。向径
最大 , 为椭球长半轴)时,曲率半径最小;向径 最小 , 为椭球短半轴)时,曲率半径最大。只有将地球看作是球体时,向径与半径才相等。 随着空间技术和计算机技术在大地测量及地图学[8]中的应用和发展,研究地球向径和常用地球半径之间的关系具有更加重要的实用价值。文献[9]提出使用不同的球体半径表示球体特性,将不同球体半径代入克拉索夫斯基椭圆参数求得数值解,再代入到球体面积公式,求得不同半径下大圆航线长度和角度变形大小,发现使用等角球半径计算大圆航线可以满足航海使用精度要求。文献[10]将子午圈曲率半径和卯酉圈曲率半径展开为大地纬度的幂级数,并借助辅助函数
导出了不同形式的平均曲率半径展开式。文献[11]利用地球向径和计算机代数系统对大地测量学和地理信息系统中的问题进行了研究。 从目前来看,前人对这一领域做了很多卓有成效的工作,但主要是针对其中一种或几种地球半径单独进行计算和应用,多是对等角纬度、等面积纬度及等距离纬度的计算与分析[12-17],很少有文献对向径进行推导和研究。为丰富对这一问题的理论研究,引入向径积分平均值和半径积分平均值的概念,并利用计算机代数系统(计算机代数系统可以将基本的数学公式展开成幂级数形式,推导出的公式相对于人工计算具有更高的精确度[18-21])将向径和半径分别进行椭圆曲线积分和椭球曲面积分,得到向径曲线积分平均值、向径曲面积分平均值、曲率半径曲线积分平均值和曲率半径曲面积分平均值。
-
地球向径
是地心O到椭球面上任意一点P的距离,如图 1所示。图 1中,X为横轴,Y为纵轴。 关于地心纬度的表达式为: 式中,
为地心纬度,地心纬度为参考椭球上一观测点和椭球中心的连线与赤道平面之间的夹角; 为椭球偏心率; 为地球短半轴。地心纬度 与大地纬度 的关系表达式为: 略去推导过程,由式(1)和式(2)得到地球向径
与大地纬度 的关系表达式为: 将式(3)展开成偏心率
的幂级数形式: 式中,
为地球长半轴,系数 ~ 为: 平均曲率半径
是描述椭球面曲率的几何量,即椭球面上某一点所有方向法截线曲率半径的算术平均值,常用于制作地球表面上局部地区地图[22],如图 2所示。 的表达式为: 式中,
为卯酉圈曲率半径; 为子午圈曲率半径。略去推导过程,得到平均曲率半径与大地纬度 的关系表达式为: 将式(7)展开成偏心率
的幂级数形式: 式(8)的系数
~ 为: 从图 1中可以看出,地球向径
是地心到椭球面上任意一点的距离。由于曲率半径是随着大地纬度 不断变化的,曲率半径不过地球球心,因此所构成的曲率圆也不同。于是,随着点位不断的变化,曲率圆中心形成了图 2所示的渐屈线。也就是说,曲率半径的中心是在不断变化的,并不是从地心到地球表面的一段线段,不能准确地表示地球半径的准确位置。从这一点来看,地球向径可以更精准地描绘地心到地球上某一点的距离,更具有代表性。 图 3为地球向径、平均曲率半径变化示意图。从图 3中可以看出,地球向径随着大地纬度的增大而减小,平均曲率半径随着大地纬度的增大而增大,地球向径与平均曲率半径在一定意义上是背离的。向径最大时,曲率半径最小;向径最小时,曲率半径最大。
-
基于地球向径和平均曲率半径关于大地纬度
的表达式,本文提出一个新的概念,即地球向径的椭圆曲线积分平均值和平均曲率半径的椭圆曲线积分平均值。以推导地球向径的椭圆曲线积分平均值为例:首先对地球向径关于大地纬度 的表达式进行椭圆曲线积分,再将经过椭圆曲线积分后的地球向径表达式取平均值,得到平面椭圆上的地球向径曲线积分平均值符号表达式。 将式(3)展开成幂级数形式,并在
范围内积分后取平均得到地球向径曲线积分平均值级数解表达式: $$ \begin{aligned} \bar{\rho}=& \frac{2}{\pi} \int \rho \mathrm{d} B=\frac{2}{\pi}\left[1-\frac{1}{2} \sin ^{2} B e^{2}+\left(\frac{\sin ^{2} B}{2}-\right.\right.\\ &\left.\frac{5 \sin ^{4} B}{8}\right) e^{4}+\left(\frac{3 \sin ^{4} B}{4}-\frac{13 \sin ^{6} B}{16}\right) e^{6}+\\ & \frac{1}{128}\left(-16 \sin ^{4} B+152 \sin ^{6} B-\right.\\ &\left.\left.141 \sin ^{8} B\right) e^{8}\right]\left.\right|_{0} ^{\pi / 2}=\\ & a\left(1-\frac{e^{2}}{4}+\frac{e^{4}}{64}+\frac{7 e^{6}}{256}+\frac{377 e^{8}}{16384}\right) \end{aligned} $$ (10) 式中,
为地球向径曲线积分平均值。 同理,对平均曲率半径关于大地纬度
的表达式进行椭圆曲线积分,再将经过椭圆曲线积分后的平均曲率半径表达式取平均,得到平面椭圆上的平均曲率半径曲线积分平均值符号表达式。 将式(7)在
范围内进行积分计算,得到平均曲率半径在平面椭圆上的曲线积分平均值: $$ \begin{aligned} \bar{R}_{a}=& \frac{2}{\pi} \int_{0}^{\pi / 2} R_{a} \mathrm{~d} B=\\ &\left.a\left[\frac{2}{\pi} \arctan \left(\sqrt{1-e^{2}} \tan B\right)\right]\right|_{0} ^{\pi / 2}=a \end{aligned} $$ (11) 式中,
为平均曲率半径曲线积分平均值。 从式(11)中可以看出,平均曲率半径曲线积分平均值为
,要大于一般意义上的地球半径曲线平均值(即地球半径平均值为 ),与通常所认为的平均曲率半径在平面椭圆上的平均值应该在椭圆长半轴 和短半轴 之间不同,证明存在系统误差。也就是说,平均曲率半径代表地球半径平均值是不精确的,这与平均曲率半径不在地心有关,然而地球向径是从地心到椭球面上任意一点的距离,将式(10)和式(11)进行对比可以看出,地球向径曲线积分平均值在某种程度上更能代表地球半径平均值进行测量和地球科学计算。 -
基于上述结果,本文提出另一个新的概念,即地球向径在地球椭球面上的曲面积分平均值和平均曲率半径在地球椭球面上的曲面积分平均值。将地球向径关于大地纬度
的表达式进行椭球面曲面积分后再取平均,得到地球向径在地球椭球面上的曲面积分平均值符号表达式。将式(3)代入到椭球面曲面积分公式中可得: 式中,
表示地球向径在椭球面上的曲面积分平均值; 表示球面积分 ; 为椭球面积的分析解。 由于此积分公式是椭圆积分,得不到分析解,因此将式(12)展开成幂级数形式再进行积分,得到地球向径曲面积分平均值级数解:
$$ \begin{aligned} {\overline{\overline \rho } }=& a\left(1-\frac{1}{6} e^{2}-\frac{17}{360} e^{4}-\frac{83}{3024} e^{6}+\right.\\ &\left.\frac{974573}{1814400} e^{8}\right) \end{aligned} $$ (13) 同理,将平均曲率半径关于大地纬度
的表达式进行椭球面曲面积分后再取平均,得到平均曲率半径在地球椭球面上的曲面积分平均值符号表达式。 将式(7)代入到椭球面曲面积分公式中得:
式中,
表示平均曲率半径在椭球曲面上的曲面积分平均值。 为便于将平均曲率半径和地球向径的椭球面曲面积分平均值进行比较,将式(7)代入到椭球面曲面积分公式中先展开成
的幂级数形式再进行积分,得到平均曲率半径椭球面曲面积分平均值的级数解: $$ \begin{aligned} {\overline{\overline R } }_{a}=& a\left(1-\frac{1}{6} e^{2}+\frac{31}{360} e^{4}+\frac{1697}{15120} e^{6}+\right.\\ &\left.\frac{1214861}{1814400} e^{8}\right) \end{aligned} $$ (15) 从式(13)和式(15)中可以看出,在地球椭球面上,地球向径和平均曲率半径的曲面积分平均值的级数解在前两项是一致的,从第3项开始产生了差距。因此,将式(13)和(15)作差得到:
由式(16)可以看出,在地球椭球面上,地球向径和平均曲率半径的曲面积分平均值差异很小,差异首项为
。由于曲率半径的中心是在不断变化的,并不能准确地表示地球半径的准确位置,从这一点来看,地球向径更加精准。 -
就地球总体而言,平均球半径(取地球椭球三轴半轴长的算术平均值,用于简单决定球半径) [23]、等面积球半径(保持球体表面积等于地球椭球面相应全面积而决定的球半径,主要用于等面积投影) [24]、等距离球半径(使球面经线总长等于地球椭球面经线总长而决定的球半径,主要用于等距离投影中)[25]和等体积球半径(使地球球体的体积等于地球椭球体的体积来决定的球半径)的表达式与大地纬度
无关,只取决于地球参考椭球模型的参数 和 ,其表达式分别为: 式中,
表示平均球半径; 表示等面积球半径; 表示等距离球半径; 表示等体积球半径。将4种常用球体半径分别与地球向径在椭球面上的曲面积分平均值和平均曲率半径在椭球面上的曲面积分平均值作差, , ,推导出其差异符号表达式列于表 1、表 2。 表 1 4种常用球体半径与向径平均值差异符号表达式
Table 1. Different Symbolical Expressions Between Four Global Radii and Integral Mean Values ofEllipsoid Radius Vector
差异符号表达式 表 2 4种常用球体半径分别与曲率半径平均值差异符号表达式
Table 2. Different Symbolical Expressions Between Four Global Radii and Mean Radius of Curvature
差异符号表达式 对比表 1、表 2可知,4种常用球体半径与向径曲面积分平均值之差比4种常用球体半径与平均曲率半径曲面积分平均值之差要小。平均球半径与向径曲面积分平均值之差首项系数绝对值是平均球半径与平均曲率半径曲面积分平均值之差的
;等面积球半径与向径曲面积分平均值之差比等面积球半径与平均曲率半径曲面积分平均值之差多了一个数量级 ;等距离球半径与向径曲面积分平均值之差第二项系数绝对值是等距离球半径与平均曲率半径曲面积分平均值之差的 ;等体积球半径与向径曲面积分平均值之差第二项系数绝对值是等体积球半径与平均曲率半径曲面积分平均值之差的 。 本文以CGCS2000参考椭球为例,对地球向径和平均曲率半径曲面积分平均值与4种球体半径差异问题进行数值分析,将地球椭球长半轴a=6 378 137 m、e=0.081 819 19代入表 1、表 2中得到表 3。
表 3 地球向径和平均曲率半径曲面积分平均值与4种球体半径数值差异/m
Table 3. Value of Different Between Integral Mean Values of Radius Vector and Mean Radius of Curvature/m
i 差异 1 2 3 4 1.59 0.003 -3 556.36 -6.38 -36.75 -38.34 -3 594.71 -44.73 观察表 3,通过数值比较看出4种常用球体半径与向径曲面积分平均值之差比4种常用球体半径与平均曲率半径曲面积分平均值之差要小,符合表 1、表 2结论。为了更直观地表现出4种常用球体半径与向径曲面积分平均值之差的差异问题,以及4种常用球体半径与平均曲率半径曲面积分平均值之间的差异问题,将椭球体的偏心率进行变化,将椭球偏心率的范围确定为
,以便更好地放大它们之间的差异关系并分析讨论,差异曲线见图 4~7。 由图 4~7可以看出,椭球偏心率
在(0,0.8)范围内变化时,4种常用球体半径与向径曲面积分平均值之差比4种常用球体半径与平均曲率半径曲面积分平均值之差小,差值更加稳定,反映出地球向径曲面积分平均值更接近4种常用球体半径。 将地球向径在椭球面上的曲面积分平均值、平均曲率半径在地球椭球面上的曲面积分平均值、平均球半径、等面积球半径、等距离球半径和等体积球半径随偏心率
在(0,0.8)的变化曲线绘制如图 8所示。由图 8可知,在椭球偏心率 在(0,0.8)范围内变化时,地球向径曲面积分平均值与4种常用球体半径之间更为接近,特别是在偏心率 逐渐变大的过程中体现得更为突出。因此,地球向径积分平均值比平均曲率半径积分平均值能更加准确地代表地球半径平均值进行测量和地球科学计算。 -
本文从地球向径和平均曲率半径的基本概念入手,引入了两个新的概念:地球向径与平均曲率半径的曲线积分平均值和曲面积分平均值,通过计算机代数系统对它们进行推导计算,并进行对比分析,为大地测量基准问题提供理论参考,得到如下结论:
1)地球是一个旋转椭球体,地球向径与平均曲率半径在一定意义上是背离的:向径最大
时,平均曲率半径最小;向径最小 时,平均曲率半径最大;只有将地球看作是球体时,向径与平均曲率半径相等。这表明平均曲率半径不能精准地代表地球半径。 2)地球向径是从地心到椭球面上任意一点的距离,地球向径曲线积分平均值介于椭球长半轴
与短半轴 之间,相比较于平均曲率半径曲线积分平均值,可以更加准确地代表地球半径的平均值进行测量和地球科学计算。 3)地球向径曲面积分平均值与4种常用球体半径之差比平均曲率半径曲面积分平均值与4种常用球体半径之差小,特别是偏心率
在(0,0.8)范围内变化时显得尤为明显。从这一点来看,地球向径积分平均值能更加准确地代表地球半径平均值进行测量和地球科学计算。 4)将地球向径与常用地球半径的差异表示为符号形式,并统一展开为偏心率
的幂级数形式。该表达式易于比较分析,一定程度上丰富了测量及地图学数学分析理论。
-
摘要: 引入地球向径积分平均值和地球平均曲率半径积分平均值的概念,借助计算机代数系统推导出了两者的符号表达式,并将它们表示为偏心率e的幂级数形式。将地球向径积分平均值和地球平均曲率半径积分平均值分别与平均球半径、等面积球半径、等距离球半径、等体积球半径这4种常用球体半径进行比较,研究表明地球向径积分平均值与4种常用球体半径间的差异更小。由于地球是一个旋转椭球体,向径与曲率半径是背离的,向径最大时,曲率半径最小, 向径最小时,曲率半径最大,传统思维所认为的曲率半径并不能准确地代表地球半径平均值,因此在一定程度上,地球向径的积分平均值更能代表地球半径平均值。这些研究结果可为地球科学、空间科学、导航定位提供基础理论依据。
-
关键词:
- 大地测量 /
- 地球向径积分表达式 /
- 平均曲率半径积分表达式 /
- 计算机代数系统 /
- 差异符号表达式
Abstract:Objectives The Earth's radial diameter and mean radius of curvature are the basic parameters commonly used in measurement and Earth science calculation.According to the requirements of Earth scien-ce and space science and some requirements, the Earth's radial diameter, mean radius of curvature, mean radius of sphere, equal distance radius of sphere, equal area radius of sphere and equal volume radius of sphere are commonly used. The concept of integral mean values of ellipsoid radius vector and mean radius of curvature is introduced.With the application and development of space technology and computer technology in geodesy and cartography, it is of more important practical value to study the relationship between the Earth's radial diameter and the common Earth radius. Methods The symbolic expressions of them are deduced by computer algebraic system and are expressed as the power series of eccentricity.We use the method that comparing the radial vector integral mean and the radius integral mean with the four common sphere radius respectively. Results The results shows that the difference between the radial integral mean and the four common sphere radii is smaller. Since the Earth is a rotating ellipsoid, the radial and the radius of curvature deviate. When the radial vector is the largest, the radius of curvature is the smallest. When the radial is the smallest, the radius of curvature is the largest. Conclusions The radius of curvature considered by traditional thinking cannot accurately represent the radius of the earth average values, to a certain extent, which means that integral mean value of ellipsoid radius vector is more representative of the average of the Earth's radius. These research results can provide theoretical basis for Earth science, space science and navigation and positioning. -
表 1 4种常用球体半径与向径平均值差异符号表达式
Table 1. Different Symbolical Expressions Between Four Global Radii and Integral Mean Values ofEllipsoid Radius Vector
差异符号表达式 表 2 4种常用球体半径分别与曲率半径平均值差异符号表达式
Table 2. Different Symbolical Expressions Between Four Global Radii and Mean Radius of Curvature
差异符号表达式 表 3 地球向径和平均曲率半径曲面积分平均值与4种球体半径数值差异/m
Table 3. Value of Different Between Integral Mean Values of Radius Vector and Mean Radius of Curvature/m
i 差异 1 2 3 4 1.59 0.003 -3 556.36 -6.38 -36.75 -38.34 -3 594.71 -44.73 -
[1] 王家耀, 成毅. 论地图学的属性和地图的价值[J]. 测绘学报, 2015, 44(3): 237-241 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201503002.htm Wang Jiayao, Cheng Yi. Discussions on the Attributes of Cartography and the Value of Map[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(3): 237-241 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201503002.htm [2] Skopeliti A, Tsoulos L. Choosing a Suitable Projection for Navigation in the Arctic[J]. Marine Geodesy, 2013, 36(2): 234-259 doi: 10.1080/01490419.2013.781087 [3] 吕志平, 张建军, 乔书波. 大地测量学基础[M]. 北京: 解放军出版社, 2005 Zhiping Lü, Zhang Jianjun, Qiao Shubo. Basis of Geodesy [M]. Beijing: Publishing House of Liberation Army, 2005 [4] 华棠. 海图数学基础[M]. 北京: 海潮出版社, 1985 Hua Tang. The Mathematical Basis of Chart [M]. Beijing: Haichao Press, 1985 [5] 王瑞, 李厚朴. 基于地球椭球模型的符号形式的航迹计算法[J]. 测绘学报, 2010, 39(2): 151-155 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201002011.htm 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 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201002011.htm [6] 李厚朴, 边少锋, 钟斌. 地理坐标系计算机代数精密分析理论[M]. 北京: 国防工业出版社, 2015 Li Houpu, Bian Shaofeng, Zhong Bin. Geographic Coordinate System Computer Algebraic Precision Analysis Theory[M]. Beijing: National Defense Industry Press, 2015 [7] 冯黎刚, 金立新, 边少锋, 等. 大椭圆椭球参数模型在长线工程中的应用[J]. 武汉大学学报∙信息科学版, 2020, 45(2): 219-225 doi: 10.13203/j.whugis20180133 Feng Ligang, Jin Lixin, Bian Shaofeng, et al. Application of Large Ellipsoid Parameter Model in Long Linear Engineering[J]. Geomatics and Information Science of Wuhan University, 2020, 45(2): 219-225 doi: 10.13203/j.whugis20180133 [8] 李志林, 刘启亮, 高培超. 地图信息论: 从狭义到广义的发展回顾[J]. 测绘学报, 2016, 45(7): 757-767 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201607002.htm Li Zhilin, Liu Qiliang, Gao Peichao. Entropy-Based Cartographic Communication Models: Evolution from Special to General Cartographic Information Theory[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(7): 757-767 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201607002.htm [9] 丁佳波. 拟大地线航法地理坐标的解算公式[M]. 天津: 海军出版社, 2003 Ding Jiabo. The Formula Calculating the Geographic Coordinates of the Approximate Geodetic Line[M]. Tianjin: Publishing House of Navy, 2003 [10] 熊介. 椭球面曲率半径展为纬差的幂级数[J]. 解放军测绘学院学报, 1985, 2(2): 15-21 https://www.cnki.com.cn/Article/CJFDTOTAL-JFJC198502003.htm 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(2): 15-21 https://www.cnki.com.cn/Article/CJFDTOTAL-JFJC198502003.htm [11] Awange J L, Grafarend E W, Paláncz B, et al. Algebraic Geodesy and Geoinformatics[M]. Berlin, Heidelberg: Springer, 2010 [12] Adams O S. Latitude Developments Connected with Geodesy and Cartography, with Tables Including a Table for Lambert Equal-Area Meridional Projection[J]. The American Mathematical Monthly, 1921, 28(11/12): 454 [13] Snyder J P. Map Projections: A Working Manual[M]. Washington DC: Government Printing Office, 1987 [14] Yang Qihe, Snyder J P, Tobler W R. Map Projection Transformation: Principles and Application[M]. London: Taylor and Francis, 2000 [15] 吴忠性, 杨启和. 数学制图学原理[M]. 北京: 测绘出版社, 1989 Wu Zhongxing, Yang Qihe. Principle of Mathematical Cartography[M]. Beijing: Surveying and Mapping Press, 1989 [16] 李忠美, 李厚朴, 边少锋. 常用纬度差异极值符号表达式[J]. 测绘学报, 2014, 43(2): 214-220 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201402018.htm 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 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201402018.htm [17] 边少锋, 纪兵. 等距离纬度等量纬度和等面积纬度展开式[J]. 测绘学报, 2007, 36(2): 218-223 doi: 10.3321/j.issn:1001-1595.2007.02.018 Bian Shaofeng, Ji Bing. The Expansions of Rectifying Latitude, Conformal Latitude and Authalic Latitude[J]. Acta Geodaetica et Cartographica Sinica, 2007, 36(2): 218-223 doi: 10.3321/j.issn:1001-1595.2007.02.018 [18] 边少锋, 李厚朴, 李忠美. 地图投影计算机代数分析研究进展[J]. 测绘学报, 2017, 46(10): 1557-1569 doi: 10.11947/j.AGCS.2017.20170396 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 [19] Bian S F, Li H P. Mathematical Analysis in Cartography by Means of Computer Algebra System[M]//Carlos Bateira. Cartography: A Tool for Spatial Analysis. Croatia: InTech, 2012 [20] 边少锋, 陈良友, 纪兵. Mathematica软件在大地测量中的应用[J]. 测绘科学与工程, 2004, 24(2): 25-28 https://www.cnki.com.cn/Article/CJFDTOTAL-JZSX201306108.htm Bian Shaofeng, Chen Liangyou, Ji Bing. Application of Mathematica in Geodesy[J]. Geomatics Science and Engineering, 2004, 24(2): 25-28 https://www.cnki.com.cn/Article/CJFDTOTAL-JZSX201306108.htm [21] 边少锋, 李厚朴. 大地测量计算机代数分析[M]. 北京: 科学出版社, 2018 Bian Shaofeng, Li Houpu. Computer Algebra Analysis on Geodesy[M]. Beijing: Science Press, 2018 [22] 陈武军. 地球曲率半径取值精度对控制测量精度的影响分析[J]. 地矿测绘, 2012, 28(4): 27-28 https://www.cnki.com.cn/Article/CJFDTOTAL-DKCH201204010.htm 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 https://www.cnki.com.cn/Article/CJFDTOTAL-DKCH201204010.htm [23] 陈丽华, 汪孔政. 关于参考椭球平均半径的探讨[J]. 测绘通报, 2000(10): 15-17 https://www.cnki.com.cn/Article/CJFDTOTAL-CHTB200010005.htm Chen Lihua, Wang Kongzheng. On Mean Radius of Reference Ellipsoid[J]. Bulletin of Surveying and Mapping, 2000(10): 15-17 https://www.cnki.com.cn/Article/CJFDTOTAL-CHTB200010005.htm [24] 李厚朴, 边少锋, 陈良友. 等面积纬度函数和等量纬度变换的直接解算公式[J]. 武汉大学学报∙信息科学版, 2011, 36(7): 843-846 http://ch.whu.edu.cn/article/id/588 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 http://ch.whu.edu.cn/article/id/588 [25] 陈成, 金立新, 李厚朴, 等. 等距离球面高斯投影[J]. 测绘通报, 2017(10): 1-6 https://www.cnki.com.cn/Article/CJFDTOTAL-CHTB201710001.htm 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 https://www.cnki.com.cn/Article/CJFDTOTAL-CHTB201710001.htm -