-
归化纬度是大地测量与地图投影变换中常用的一种纬度,在地图投影理论分析中,常引入归化纬度来简化问题。例如,基于归化纬度的子午线弧长微分表达式更简单,与第二类椭圆积分的联系更加密切,精度也更高[1-3];归化纬度表示的子午线弧长与大地主题解算具有理论上的一致性[4-5];归化纬度也是贝塞尔大地主题解算方法中重要的中间变量[6],在椭球大地测量的许多领域,归化纬度都可以使复杂的问题得以简化[7-9]。除归化纬度外,在地图投影变换中,也常会用到大地纬度、地心纬度、等距离纬度、等面积纬度、等角纬度5种纬度[10-13]。为了分析6种纬度间的关系,文献[14]利用计算机代数系统Mathematica[15-18]分析了辅助纬度与大地纬度以及辅助纬度间的差值问题;文献[19]提出了针对含有小参数的复杂方程求解的迭代解法,分析了常用纬度与地心纬度的差异极值问题。
从目前看,常用纬度与归化纬度差异一般采用数值形式表示,没有对应的分析表达式,不适用于分析不同参数的椭球。本文着重研究常用纬度与归化纬度的差异极值问题,利用计算机代数分析方法,根据归化纬度与其他常用纬度间函数关系,采用针对一般小参数方程的迭代解法,分析常用纬度与归化纬度间的差异极值问题,并以克拉索夫斯基椭球等4种常用椭球模型为例进行算例分析,分析常用纬度与归化纬度之间的差异。
-
如图 1所示,以子午面椭圆中心为原点
,长半轴 方向为 轴,短半轴 方向为y轴,建立平面直角坐标系。设地球椭球面上任意一点为 ,作以原点为中心、半径为 的辅助圆,延长 点的纵坐标线与圆交于 ,连接 ,则 为 点的归化纬度, 点法线方向与赤道平面的夹角为 点的大地纬度 , 为 点的地心纬度。大地纬度 、地心纬度 和归化纬度 之间的数学表达式为: ]]>t a n ϕ = 1 - e 2 t a n u = 1 - e 2 t a n B 式中,e为第一偏心率。
根据式(1),大地纬度
及地心纬度 以归化纬度 为变量的表达式分别为: B = a r c t a n ( t a n u 1 - e 2 ) ]]>ϕ = a r c t a n ( 1 - e 2 t a n u ) ]]>等距离纬度
的定义为[20]: ψ = X R 1 = ]]>∫ 0 B 1 - e 2 1 - e 2 s i n 2 B 3 / 2 d B 1 - e 2 K 0 式中,X为子午线弧长;
;R1为常数,等于子午线弧长X积分至π/2时对应的数值。 根据子午线弧长与归化纬度
之间的关系[1],等距离纬度 与归化纬度 之间的表达式为: ψ = X R 1 = ]]>∫ 0 u 1 - e 2 c o s 2 u d u 1 - e 2 K 0 同理,根据等面积纬度的定义[21],可得等面积纬度
与归化纬度 之间的表达式为: s i n ϑ = F R 2 = ]]>∫ 0 u c o s u 1 - e 2 c o s 2 u d u ∫ 0 π / 2 c o s u 1 - e 2 c o s 2 u d u 式中,等面积纬度函数
可积;R2为常数,等于等面积纬度函数F积分至π/2时对应的数值。基于归化纬度的等面积纬度函数表达式为: ]]>F = 1 - e 2 2 e a r t a n h e s i n u 1 - e 2 c o s 2 u + 1 2 s i n u 1 - e 2 c o s 2 u ( 7 ) 借助中间变量等量纬度
[22],等角纬度 与归化纬度 之间的关系可表示为: φ = a r c s i n ( t a n h q ) ]]> ]]>q = ∫ 0 u 1 - e 2 c o s 2 u c o s u d u = a r t a n h ( s i n u 1 - e 2 c o s 2 u ) - e a r c s i n ( e s i n u 1 - e 2 ) -
大地纬度、地心纬度、等距离纬度与归化纬度差异极值点可根据表达式直接求得,无需将这些纬度表示为归化纬度的级数展开式。等面积纬度、等角纬度与归化纬度差值表达式复杂,对表达式直接求导后无法计算得到导数为0的点,故需级数展开。
以等面积纬度为例,由式(6)可知等面积纬度与归化纬度的关系,将等面积纬度基于椭球偏心率
幂级数展开,截去高阶项,取至 ,对展开式逐项积分,化三角函数幂形式为倍角形式,整理可得等面积纬度关于归化纬度的正解表达式为: ]]>ϑ = u + α 2 s i n 2 u + α 4 s i n 4 u + α 6 s i n 6 u + α 8 s i n 8 u + α 10 s i n 10 u ( 10 ) 其中,
]]>α 2 = - 1 12 e 2 - 17 360 e 4 - 121 4 032 e 6 - 6 203 302 400 e 8 - 3 535 339 239 500 800 e 10 α 4 = - 7 1 440 e 4 - 17 3 360 e 6 - 15 997 3 628 800 e 8 - 441 499 119 750 400 e 10 α 6 = - 83 181 440 e 6 - 311 453 600 e 8 - 72 437 95 800 320 e 10 α 8 = - 797 14 515 200 e 8 - 12 917 119 750 400 e 10 α 10 = - 3 673 479 001 600 e 10 同理可得由归化纬度直接求解等角纬度的表达式为:
φ = u + β 2 s i n 2 u + β 4 s i n 4 u + ]]>β 6 s i n 6 u + β 8 s i n 8 u + β 10 s i n 10 u ]]>其中,
]]>β 2 = 1 4 e 2 - 1 12 e 4 - 7 192 e 6 - 113 5 760 e 8 - 7 576 e 10 β 4 = 1 96 e 4 + 1 240 e 6 + 1 13 440 e 8 + 29 80 640 e 10 β 6 = - 1 960 e 6 - 13 13 440 e 8 - 83 107 520 e 10 β 8 = 17 322 560 e 8 + 1 32 256 e 10 β 10 = - 1 107 520 e 10 -
为得到常用纬度与归化纬度之间差异极值的分析表达式,将常用纬度与归化纬度的差值表达式对归化纬度进行求导,求出当
时的极值点 。由地球椭球的对称性可知,极值点在南北半球呈对称分布,可求得 时的极值点。 由式(2)知大地纬度与归化纬度的关系式,解得大地纬度与归化纬度的差值关于归化纬度的导数为0时对应的归化纬度值为:
d ( B - u ) d u = 1 - e 2 s e c 2 u 1 - e 2 + t a n 2 u - 1 = 0 ]]>解得
时大地纬度与归化纬度差异极值点与对应极值为: ]]>u m = a r c t a n 1 - e 2 4 B - u m a x = a r c t a n ( 1 / 1 - e 2 4 - 1 - e 2 4 2 ) 此时二阶导数小于0,该极值点为极大值点,将差异极值点与对应极值展开为
的幂级数形式,取至 ,结果见表 1。 表 1 常用纬度与归化纬度差异极值基于椭球第一偏心率e展开的分析表达式
Table 1. Analytical Expressions of Difference Extrema Between Commonly Used and Reduced Latitudes Expressed in a Power Series of Eccentricity e
差值 差异极值点 对应差异极值 同理可求得
时地心纬度与归化纬度差异极值点与对应极值为: ]]>u m = a r c t a n ( 1 / 1 - e 2 4 ) ϕ - u m i n = a r c t a n ( 1 - e 2 4 - 1 / 1 - e 2 4 2 ) 此时二阶导数大于0,该极值点为极小值点,将差异极值点与对应极值展开为
的幂级数形式,取至 ,结果见表 1。由式(13)和式(14)知,大地纬度、地心纬度与归化纬度差异极值互为相反数。 结合式(6)知,等距离纬度与归化纬度差值
关于归化纬度 的导数为: d ( ψ - u ) d u = 1 - e 2 c o s 2 u 1 - e 2 K 0 - 1 = 0 ]]>解得
时等距离纬度与归化纬度差异极值点与对应极值为: ]]>u m = a r c c o s 1 e 1 - 1 - e 2 2 K 0 2 ψ - u m i n = ∫ 0 u m 1 - e 2 c o s 2 u d u 1 - e 2 K 0 - u m 此时二阶导数大于0,该极值点为极小值点,将差异极值点与对应极值展开为
的幂级数形式,取至 ,结果见表 1。 -
将等面积纬度与归化纬度差值
关于归化纬度 的导数展开为 的幂级数形式,取至 ,可以表示为: ]]>d ( ϑ - u ) d u = A 2 ( u ) e 2 + A 4 ( u ) e 4 + A 6 ( u ) e 6 + A 8 ( u ) e 8 + A 10 ( u ) e 10 ( 17 ) 其中,
]]>A 2 = - 1 6 c o s 2 u A 4 = - 17 180 c o s 2 u - 7 360 c o s 4 u A 6 = - 121 2 016 c o s 2 u - 17 840 c o s 4 u - 83 30 240 c o s 6 u A 8 = - 6 203 151 200 c o s 2 u - 15 997 907 200 c o s 4 u - 311 75 600 c o s 6 u - 797 1 814 400 c o s 8 u A 10 = - 3 535 339 119 750 400 c o s 2 u - 441 499 29 937 600 c o s 4 u - 72 437 15 966 720 c o s 6 u - 12 917 14 968 800 c o s 8 u - 3 673 47 900 160 c o s 10 u 由式(17)可知,由于
,故导数为0,即: ∑ k = 1 5 A 2 k ( u ) e 2 k - 2 = 0 ]]>为求解式(18),根据文献[19]中的迭代方法,先求
的反函数为: u = 1 2 a r c c o s ( - 6 t ) ]]>解得式(18)的近似解 ,则极值点在 附近,令迭代初值 建立迭代式: u l + 1 = 1 2 a r c c o s ( 6 ∑ k = 2 5 A 2 k ( u l ) e 2 k - 2 ) ]]>将第
次迭代结果 代入式(20)计算第 次迭代结果,每一次迭代所得结果均展开至 ,截去高阶项。当迭代至第5次时,展开式系数已无变化,至此迭代完成,得到迭代结果为: u m = π 4 - 7 120 e 2 - 697 25 200 e 4 - 296 263 18 144 000 e 6 - 16 876 243 1 552 320 000 ]]>e 8 此时二阶导数大于0,该极值点为极小值点,计算迭代得到的
与 之间的差异 ,展开至 ,截去高阶项,得到 时极值点对应的差异极小值为: ϑ - u m i n = - 1 12 e 2 - 17 360 e 4 - 7 807 259 200 e 6 - 375 329 18 144 000 ]]>e 8 同理求得
时,差异值 的极值为: d ( φ - u ) d u = ∑ k = 1 5 ]]>B 2 k ( u ) e 2 k 令迭代初值
建立迭代式: u l + 1 = 1 2 a r c c o s ( 2 ∑ k = 2 5 B 2 k ( u l ) e 2 k - 2 ) ]]>得到相应差异极小值点与对应的差异极小值,基于椭球偏心率e展开至
,结果见表 1。 传统的椭球大地测量和地图投影数学问题多基于偏心率e进行展开,表达式较为繁琐冗长,近年来椭球第三扁率n经常被用于大地测量计算中,它与第一偏心率e的关系为:
n = a - b a + b = ]]>1 - 1 - e 2 1 + 1 - e 2 由式(23)可知,椭球第三扁率n在数值上比偏心率e更小,收敛速度也更快,可在简化表达式的同时提高精度。将
时,差异值 、 的极值点和对应的差异极值表达式基于椭球第三扁率 展开至 ,结果见表 2。 表 2 常用纬度与归化纬度差异极值基于椭球第三扁率n展开的分析表达式
Table 2. Analytical Expressions of Difference Extrema Between Commonly Used and Reduced Latitudes Expressed in a Power Series of the Third Flattening n
差值 差异极值点 对应差异极值 由表 1、表 2可知,基于椭球第三扁率n展开的分析表达式在表示差异值
、 、 的极值时,只有奇次项,没有二次项与四次项,比基于偏心率 展开的分析表达式形式上更加简洁。差异值 、 、 的极值点位于 右侧, 、 的极值点位于 左侧, 、 的极值点关于 左右对称,不仅如此, 和 互为相反数且绝对值最大。 因为e和n的数值很小、收敛速度快,
约为0.006 694, 约为 , 约为0.001 679, 约为 ,故差异极值的数值大小和正负主要由首项系数,即 项或 项系数决定,此项系数越大,对应差异极值越大,此项系数的正负直接决定差异极值的正负。当差异极值基于椭球第三扁率n展开时,差异极值 、 、 系数的绝对值为1, 系数的绝对值为1/2, 系数的绝对值为1/3,综合考虑差异极值的正负,则可以近似认为这5种常用纬度与归化纬度的差异极值满足 。 -
§2.1、§2.2中分别推导了以椭球第一偏心率e与第三扁率n展开的常用纬度与归化纬度差异极值的分析表达式,便于理论分析,适用于不同椭球参数。为直观体现常用纬度与归化纬度之间差异极值的数值大小情况,以克拉索夫斯基椭球、IUGG1975椭球、WGS84椭球和CGCS2000椭球4种参考椭球为例,计算对应参考椭球下差异极值的数值,4种参考椭球的e与n的数值如表 3所示。
表 3 4种常用参考椭球参数
Table 3. Four Commonly Used Referenced Ellipsoids
参数 克拉索夫斯基椭球 0.081 813 334 0 0.001 678 980 657 IUGGl975椭球 0.081 819 221 4 0.001 679 244 312 WGS84椭球 0.081 819 190 8 0.001 679 220 386 CGCS2000椭球 0.081 819 191 0 0.001 679 220 394 将4种椭球参数分别代入常用纬度与归化纬度的差异极值表达式,得到对应参考椭球模型下的常用纬度与归化纬度差异极值的数值,结果见表 4。其中,WGS84椭球和CGCS2000椭球下常用纬度与归化纬度差异极值点与差异极值在0.01″的精度下完全一致。
表 4 常用纬度与归化纬度差异极值
Table 4. Difference Extrema Between Commonly Used and Reduced Latitudes
参考椭球 差值 差异极值点 差异极值 克拉索夫斯基椭球 IUGGl975椭球 WGS84椭球/CGCS2000椭球 以CGCS2000参考椭球为例,在
范围内,绘制常用纬度与归化纬度间差异示意曲线如图 2所示。 由图 2与表 4可知,在
范围内,各常用纬度与归化纬度差异的绝对值先增大后减小,且在归化纬度为 附近出现极值。由文献[14]中的计算结果可知,常用纬度与大地纬度的差异极值绝对值最小为 ,与其相比,常用纬度与归化纬度的差值明显较小,绝对值最大仅为 ,相当于常用纬度与大地纬度的差异极值绝对值最小值,即归化纬度与常用纬度之间差值的绝对值比大地纬度与常用纬度差值绝对值小,归化纬度在数值上与其他常用纬度更为接近。 -
本文对5种常用纬度与归化纬度差异极值进行了研究,将差异极值点与差异极值分别表示为椭球第一偏心率
和第三扁率n的幂级数形式,并进行了算例分析,得出如下结论。 1)极值点位置方面。在归化纬度从
°到 °变化的过程中,各常用纬度与归化纬度差值的绝对值先增大后减小,在归化纬度为45 附近出现极值。 2)差异极值的大小。差异值
、 的极值互为相反数且绝对值最大, 的极值最小。可近似认为5种常用纬度与归化纬度差异极值满足3( - )min≈2( - )min≈( - )min 。 3)基于椭球第三扁率n展开的分析表达式在表示差异值
、 、 的极值时,没有偶次项,比基于椭球第一偏心率e展开的分析在表达式形式上更加简洁。 本文推导的常用纬度与归化纬度差异极值分析表达式,更加直观地显示了常用纬度与归化纬度的差异,适用于分析不同参数的椭球。
Analytical Expressions of Differences Between Commonly Used Latitudes and Reduced Latitude
-
摘要: 从常用纬度与归化纬度的定义出发,借助于计算机代数分析,对常用纬度与归化纬度差值进行了系统的分析,推导出了常用纬度与归化纬度的差异极值分析表达式,并将式中系数展开为椭球第一偏心率
和第三扁率n的幂级数形式。研究结果表明,常用纬度与归化纬度差值极值点均位于45 附近,基于n展开的分析表达式比基于 展开的分析表达式形式上更加简洁,该分析表达式适用于任何椭球参数,具有通用性,一定程度上丰富了椭球几何学常用纬度分析理论。 Abstract:Objectives The reduced latitude can simplify many complex problems in ellipsoid geodesy. In addition to the reduced latitude, geodetic, geocentric, rectifying, authalic, and conformal latitudes are commonly used in map projection transformation. The differences between commonly used latitudes and reduced latitude are generally expressed in numerical form, and there are no corresponding analytical expressions. Analytical expressions of differences between commonly used and reduced latitudes can visually display the differences between latitudes. Methods Starting from the definition of commonly used and reduced latitudes, with the aid of computer algebra analysis methods, the function relations between commonly used latitudes and reduced latitude are derived, then the differences between commonly used and reduced latitudes can be obtained by finding the extreme points. Analytical expressions of differences between them are derived and their coefficients are expressed in a power series of the eccentricity e, and the third flattening n. Results Theoretical and numerical analysis results show that the analytical expressions in terms of n are more concise than those in terms of e , the absolute value of the differences between commonly used latitudes and reduced latitude first increase and then decrease, and the extreme points are all near , the reduced latitude is closer to commonly used latitudes in numerical values. Conclusions Expressions of differences between commonly used latitudes and reduced latitude are more intuitively and applicable to any referenced ellipsoid, to some extent, they could enrich the analytical theory of the various commonly used latitudes. -
Key words:
- reduced latitude /
- commonly used latitudes /
- difference extrema /
- analytical expressions
-
表 1 常用纬度与归化纬度差异极值基于椭球第一偏心率e展开的分析表达式
Table 1. Analytical Expressions of Difference Extrema Between Commonly Used and Reduced Latitudes Expressed in a Power Series of Eccentricity e
差值 差异极值点 对应差异极值 表 2 常用纬度与归化纬度差异极值基于椭球第三扁率n展开的分析表达式
Table 2. Analytical Expressions of Difference Extrema Between Commonly Used and Reduced Latitudes Expressed in a Power Series of the Third Flattening n
差值 差异极值点 对应差异极值 表 3 4种常用参考椭球参数
Table 3. Four Commonly Used Referenced Ellipsoids
参数 克拉索夫斯基椭球 0.081 813 334 0 0.001 678 980 657 IUGGl975椭球 0.081 819 221 4 0.001 679 244 312 WGS84椭球 0.081 819 190 8 0.001 679 220 386 CGCS2000椭球 0.081 819 191 0 0.001 679 220 394 表 4 常用纬度与归化纬度差异极值
Table 4. Difference Extrema Between Commonly Used and Reduced Latitudes
参考椭球 差值 差异极值点 差异极值 克拉索夫斯基椭球 IUGGl975椭球 WGS84椭球/CGCS2000椭球 -
[1] 过家春, 李厚朴, 庄云玲, 等. 依不同纬度变量的子午线弧长正反解公式的级数展开[J]. 测绘学报, 2016, 45(5): 560-565 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201605008.htm Guo Jiachun, Li Houpu, Zhuang Yunling, et al. Series Expansion for Direct and Inverse Solutions of Meridian in Terms of Different Latitude Variables[J]. Acta Geodaetica et Cartographica Sinica, 2016, 45(5): 560-565 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201605008.htm [2] 边少锋, 李厚朴, 李忠美. 地图投影计算机代数分析研究进展[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 [3] 叶彤, 李厚朴, 钟业勋, 等. 3种辅助纬度关于归化纬度的正反解表达式[J]. 测绘科学技术学报, 2020, 37(1): 106-110 https://www.cnki.com.cn/Article/CJFDTOTAL-JFJC202001021.htm Ye Tong, Li Houpu, Zhong Yexun, et al. The Forward and Inverse Expressions of the Three Auxiliary Latitudes About the Reduced Latitude[J]. Journal of Geomatics Science and Technology, 2020, 37 (1): 106-110 https://www.cnki.com.cn/Article/CJFDTOTAL-JFJC202001021.htm [4] Bermejo-Solera M, Otero J. Simple and Highly Accurate Formulas for the Computation of Transverse Mercator Coordinates from Longitude and Isometric Latitude[J]. Journal of Geodesy, 2009, 83(1): 1-12 doi: 10.1007/s00190-008-0224-y [5] 史国友, 周晓明, 贾传荧. 贝塞尔大地主题正解的改进算法[J]. 大连海事大学学报, 2008, 34(1): 15-19 doi: 10.3969/j.issn.1671-7031.2008.01.005 Shi Guoyou, Zhou Xiaoming, Jia Chuanying. Improved Algorithm for Direct Solution of Bessel's Geodetic Problem[J]. Journal of Dalian Maritime University, 2008, 34(1): 15-19 doi: 10.3969/j.issn.1671-7031.2008.01.005 [6] 杨启和. 地图投影变换原理与方法[M]. 北京: 解放军出版社, 1989 Yang Qihe. The Principle and Method of Map Projection Transformation[M]. Beijing: Publishing House of Liberation Army, 1989 [7] 边少锋, 李厚朴. 大地测量计算机代数分析[M]. 北京: 科学出版社, 2018 Bian Shaofeng, Li Houpu. Computer Algebra Analysis on Geodesy[M]. Beijing: Science Press, 2018 [8] 熊介. 椭球大地测量学[M]. 北京: 解放军出版社, 1988 Xiong Jie. Ellipsoidal Geodesy[M]. Beijing: Publishing House of Liberation Army, 1988 [9] 任留成. 空间地图投影原理[M]. 北京: 测绘出版社, 2013 Ren Liucheng. The Principle of Space Map Projection[M]. Beijing: Surveying and Mapping Publishing House, 2013 [10] 吕晓华, 李少梅. 地图投影原理与方法[M]. 北京: 测绘出版社, 2016 Lü Xiaohua, Li Shaomei. Principles and Methods of Map Projection[M]. Beijing: Surveying and Mapping Publishing House, 2016 [11] 边少锋, 李厚朴. 高斯投影的复变函数表示[M]. 北京: 科学出版社, 2021 Bian Shaofeng, Li Houpu. The Expressions of Gauss Projection by Complex Numbers[M]. Beijing: Science Press, 2004 [12] 李厚朴, 边少锋. 高斯投影与墨卡托投影解析变换的复变函数表达式[J]. 武汉大学学报·信息科学版, 2009, 34(3): 277-279 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200903006.htm Li Houpu, Bian Shaofeng. Expressions for Analytical Transformation Between Gauss and Mercator Projections by Complex Numbers[J]. Geomatics and Information Science of Wuhan University, 2009, 34(3): 277-279 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200903006.htm [13] 李厚朴, 李海波, 唐庆辉. 椭球情形下等角和等面积正圆柱投影间的直接变换[J]. 海洋技术学报, 2019, 38(5): 15-20 https://www.cnki.com.cn/Article/CJFDTOTAL-HYJS201905003.htm Li Houpu, Li Haibo, Tang Qinghui. The Direct Transformations Between Conformal and Equal-Area Normal Cylindrical Projections Based on Earth Ellipsoid[J]. Journal of Ocean Technology, 2019, 38(5): 15-20 https://www.cnki.com.cn/Article/CJFDTOTAL-HYJS201905003.htm [14] 李忠美, 李厚朴, 边少锋. 常用纬度差异极值符号表达式[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 [15] 李厚朴, 边少锋, 刘敏. 地图投影中三种纬度间变换直接展开式[J]. 武汉大学学报·信息科学版, 2013, 38(2): 217-220 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH201302022.htm Li Houpu, Bian Shaofeng, Liu Min. Direct Expansions of Transformations Between Three Kinds of Latitudes Used in Map Projection[J]. Geomatics and Information Science of Wuhan University, 2013, 38(2): 217-220 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH201302022.htm [16] 郑彤, 边少锋. 子午线弧长反问题新解[J]. 武汉大学学报·信息科学版, 2007, 32 (3): 255-258 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200703016.htm Zheng Tong, Bian Shaofeng. New Solutions on Inverse Problem of Meridian Arc[J]. Geomatics and Information Science of Wuhan University, 2007, 32 (3): 255-258 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200703016.htm [17] 李厚朴, 边少锋. 辅助纬度反解公式的Hermite插值法新解[J]. 武汉大学学报·信息科学版, 2008, 33 (6): 623-626 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200806020.htm Li Houpu, Bian Shaofeng. Derivation of Inverse Expansions for Auxiliary Latitudes by Hermite Interpolation Method[J]. Geomatics and Information Science of Wuhan University, 2008, 33(6): 623-626 https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200806020.htm [18] 宗敬文, 李厚朴, 边少锋, 等. 常用地球半径差异符号表达式[J]. 测绘学报, 2019, 48(2): 238-244 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201902014.htm Zong Jingwen, Li Houpu, Bian Shaofeng, et al. Symbolic Expressions of Differences Between Earth Radius[J]. Acta Geodaetica et Cartographica Sinica, 2019, 48(2): 238-244 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201902014.htm [19] 陈成, 金立新, 李忠美. 常用纬度与地心纬度差异分析[J]. 海军工程大学学报, 2018, 30 (2): 72-77 https://www.cnki.com.cn/Article/CJFDTOTAL-HJGX201802014.htm Chen Cheng, Jin Lixin, Li Zhongmei. Differences Between Regular and Geocentric Latitudes[J]. Journal of Naval University of Engineering, 2018, 30 (2): 72-77 https://www.cnki.com.cn/Article/CJFDTOTAL-HJGX201802014.htm [20] 李松林, 李厚朴, 边少锋, 等. 等距离纬度与等量纬度间的直接变换[J]. 测绘科学技术学报, 2018, 35(6): 643-647 https://www.cnki.com.cn/Article/CJFDTOTAL-JFJC201806018.htm Li Songlin, Li Houpu, Bian Shaofeng, et al. Direct Transformation Between Rectifying Latitude and Isometric Latitude[J]. Journal of Geomatics Science and Technology, 2018, 35(6): 643-647 https://www.cnki.com.cn/Article/CJFDTOTAL-JFJC201806018.htm [21] 李松林, 李厚朴, 边少锋, 等. 等面积纬度和等量纬度间的直接变换[J]. 海洋测绘, 2019, 39(3): 26-30 https://www.cnki.com.cn/Article/CJFDTOTAL-HYCH201903006.htm Li Songlin, Li Houpu, Bian Shaofeng, et al. Direct Expansions of Transformation between Authalic Latitude and Isometric Latitude[J]. Hydrographic Surveying and Charting, 2019, 39(3): 26-30 https://www.cnki.com.cn/Article/CJFDTOTAL-HYCH201903006.htm [22] 孔祥元, 郭际明, 刘宗泉. 大地测量学基础[M]. 武汉: 武汉大学出版社, 2001 Kong Xiangyuan, Guo Jiming, Liu Zongquan. Geodetic Foundation[M]. Wuhan: Wuhan University Press, 2001 -