留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

地球椭球向径和平均曲率半径的积分表达式

宗敬文 李厚朴 钟业勋

宗敬文, 李厚朴, 钟业勋. 地球椭球向径和平均曲率半径的积分表达式[J]. 武汉大学学报 ● 信息科学版, 2022, 47(7): 1063-1070. doi: 10.13203/j.whugis20200059
引用本文: 宗敬文, 李厚朴, 钟业勋. 地球椭球向径和平均曲率半径的积分表达式[J]. 武汉大学学报 ● 信息科学版, 2022, 47(7): 1063-1070. doi: 10.13203/j.whugis20200059
ZONG Jingwen, LI Houpu, ZHONG Yexun. Integral Expressions of Earth Ellipsoid Radius Vector and Mean Radius of Curvature[J]. Geomatics and Information Science of Wuhan University, 2022, 47(7): 1063-1070. doi: 10.13203/j.whugis20200059
Citation: ZONG Jingwen, LI Houpu, ZHONG Yexun. Integral Expressions of Earth Ellipsoid Radius Vector and Mean Radius of Curvature[J]. Geomatics and Information Science of Wuhan University, 2022, 47(7): 1063-1070. doi: 10.13203/j.whugis20200059

地球椭球向径和平均曲率半径的积分表达式

doi: 10.13203/j.whugis20200059
基金项目: 

国家自然科学基金 41671459

国家自然科学基金 41871376

国家自然科学基金 41771487

湖北省杰出青年科学基金 2019CFA086

详细信息
    作者简介:

    宗敬文,博士生,主要从事几何大地测量学和物理大地测量的理论与方法研究。zjw19950613@163.com

    通讯作者: 李厚朴,博士,副教授。lihoupu1985@126.com
  • 中图分类号: P282

Integral Expressions of Earth Ellipsoid Radius Vector and Mean Radius of Curvature

Funds: 

The National Natural Science Foundation of China 41671459

The National Natural Science Foundation of China 41871376

The National Natural Science Foundation of China 41771487

the Hubei Province Science Fund for Distinguished Young Scholars 2019CFA086

More Information
    Author Bio:

    ZONG Jingwen, PhD candidate, specializes in the theories and methods of geometric geodesy and physical geodesy. E-mail: zjw19950613@163.com

    Corresponding author: LI Houpu, PhD, associate professor. E-mail: lihoupu1985@126.com
  • 摘要: 引入地球向径积分平均值和地球平均曲率半径积分平均值的概念,借助计算机代数系统推导出了两者的符号表达式,并将它们表示为偏心率e的幂级数形式。将地球向径积分平均值和地球平均曲率半径积分平均值分别与平均球半径、等面积球半径、等距离球半径、等体积球半径这4种常用球体半径进行比较,研究表明地球向径积分平均值与4种常用球体半径间的差异更小。由于地球是一个旋转椭球体,向径与曲率半径是背离的,向径最大时,曲率半径最小, 向径最小时,曲率半径最大,传统思维所认为的曲率半径并不能准确地代表地球半径平均值,因此在一定程度上,地球向径的积分平均值更能代表地球半径平均值。这些研究结果可为地球科学、空间科学、导航定位提供基础理论依据。
  • 图  1  地球向径示意图

    Figure  1.  Digram of Earth Radius Vector

    图  2  平均曲率半径示意图

    Figure  2.  Diagram of Average Curvature Radius

    图  3  地球向径、平均曲率半径变化示意图

    Figure  3.  Diagram of Changes About Earth Radial Direction and Average Radius of Curvature

    图  4  Δ1D1线

    Figure  4.  Change Graph of Δ1,D1

    图  5  Δ2D2线

    Figure  5.  Change Graph of Δ2,D2

    图  6  Δ3D3线

    Figure  6.  Change Graph of Δ3,D3

    图  7  Δ4D4线

    Figure  7.  Change Graph of Δ4D4

    图  8  变化曲线图

    Figure  8.  Change Graph of R¯¯a,ρ¯¯,Re,RF,RS,RV

    表  1  4种常用球体半径与向径平均值差异符号表达式

    Table  1.   Different Symbolical Expressions Between Four Global Radii and Integral Mean Values ofEllipsoid Radius Vector

    Δi 差异符号表达式
    Δ1 a(1180e4+5756e6-499 099907 200e8)
    Δ2 a(1189e6-2 2274 050e8)
    Δ3 a(-112e2+12 880e4+38348 384e6-127 225 969232 243 200e8)
    Δ4 a(-145e4-171 134e6-192 959340 200 e8)
    下载: 导出CSV

    表  2  4种常用球体半径分别与曲率半径平均值差异符号表达式

    Table  2.   Different Symbolical Expressions Between Four Global Radii and Mean Radius of Curvature

    Di 差异符号表达式
    D1 a(-23180e4-5033 780e6-619 243907 200e8)
    D2 a(-215e4-127945e6-38 68756 700e8)
    D3 a(-112e2-3832 880e4-31 877241 920e6-157 982 833232 243 200e8)
    D4 a(-745e4-8775 670e6-238 013340 200e8)
    下载: 导出CSV

    表  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
    Δi 1.59 0.003 -3 556.36 -6.38
    Di -36.75 -38.34 -3 594.71 -44.73
    下载: 导出CSV
  • [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
  • [1] 叶彤, 李厚朴, 钟业勋, 金立新.  常用纬度与归化纬度差异极值分析表达式 . 武汉大学学报 ● 信息科学版, 2022, 47(3): 473-480. doi: 10.13203/j.whugis20200004
    [2] 李金岭, 黄飞, 孙中苗.  Ka波段天测与测地VLBI观测的利弊分析 . 武汉大学学报 ● 信息科学版, 2017, 42(2): 257-262. doi: 10.13203/j.whugis20140723
    [3] 王磊, 郭际明, 申丽丽, 高奋生.  顾及椭球面不平行的椭球膨胀法在高程投影面变换中的应用 . 武汉大学学报 ● 信息科学版, 2013, 38(6): 725-728.
    [4] 李厚朴, 边少锋, 刘敏.  地图投影中三种纬度间变换直接展开式 . 武汉大学学报 ● 信息科学版, 2013, 38(2): 217-220.
    [5] 李厚朴, 边少锋, 陈良友.  等面积纬度函数和等量纬度变换的直接解算公式 . 武汉大学学报 ● 信息科学版, 2011, 36(7): 843-846.
    [6] 李世祥, 樊红, 王玉莉.  基因表达式程序设计的卫星遥感影像恢复 . 武汉大学学报 ● 信息科学版, 2010, 35(7): 877-.
    [7] 朱智勤, 李征航, 刘万科.  相位中心改正模式的转变对GPS数据处理的影响 . 武汉大学学报 ● 信息科学版, 2009, 34(11): 1301-1304.
    [8] 李厚朴, 边少锋.  高斯投影与墨卡托投影解析变换的复变函数表达式 . 武汉大学学报 ● 信息科学版, 2009, 34(3): 277-279.
    [9] 李厚朴, 边少锋.  辅助纬度反解公式的Hermite插值法新解 . 武汉大学学报 ● 信息科学版, 2008, 33(6): 623-626.
    [10] 杨元喜, 曾安敏.  大地测量数据融合模式及其分析 . 武汉大学学报 ● 信息科学版, 2008, 33(8): 771-774.
    [11] 刘经南, 陈庭, 许才军, 易长荣.  球坐标下数值流形方法有限单元覆盖的单元矩阵表达式 . 武汉大学学报 ● 信息科学版, 2005, 30(12): 1035-1038.
    [12] 陈俊勇.  地球重力场、重力、重力梯度在三维直角坐标系中的表达式 . 武汉大学学报 ● 信息科学版, 2004, 29(5): 377-379. doi: 10.13203/j.whugis2004.05.001
    [13] 宁津生, 汪海洪, 罗志才.  小波分析在大地测量中的应用及其进展 . 武汉大学学报 ● 信息科学版, 2004, 29(8): 659-663.
    [14] 王永尚.  国家大地测量数据库网络化建设 . 武汉大学学报 ● 信息科学版, 2003, 28(S1): 41-43,47.
    [15] 许才军.  用大地测量资料检验青藏高原隆升的两种构造演化模式 . 武汉大学学报 ● 信息科学版, 1995, 20(2): 163-167.
    [16] 许才军, 陶本藻, 晁定波.  论大地测量研究青藏高原地学问题的方法 . 武汉大学学报 ● 信息科学版, 1995, 20(1): 57-61.
    [17] 金标仁, 申文斌.  相对论大地测量 . 武汉大学学报 ● 信息科学版, 1994, 19(3): 227-231.
    [18] 张德涵.  关于青藏高原地壳的隆升与增厚 . 武汉大学学报 ● 信息科学版, 1993, 18(3): 34-40.
    [19] 钟业勋, 胡毓钜.  地图的集合模型(表达式)及比较应用初探 . 武汉大学学报 ● 信息科学版, 1990, 15(1): 58-65.
    [20] 单杰.  联合平差中大地测量观测值的数学模型 . 武汉大学学报 ● 信息科学版, 1988, 13(4): 105-111.
  • 加载中
图(8) / 表(3)
计量
  • 文章访问数:  108
  • HTML全文浏览量:  49
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-02-27
  • 刊出日期:  2022-07-05

地球椭球向径和平均曲率半径的积分表达式

doi: 10.13203/j.whugis20200059
    基金项目:

    国家自然科学基金 41671459

    国家自然科学基金 41871376

    国家自然科学基金 41771487

    湖北省杰出青年科学基金 2019CFA086

    作者简介:

    宗敬文,博士生,主要从事几何大地测量学和物理大地测量的理论与方法研究。zjw19950613@163.com

    通讯作者: 李厚朴,博士,副教授。lihoupu1985@126.com
  • 中图分类号: P282

摘要: 引入地球向径积分平均值和地球平均曲率半径积分平均值的概念,借助计算机代数系统推导出了两者的符号表达式,并将它们表示为偏心率e的幂级数形式。将地球向径积分平均值和地球平均曲率半径积分平均值分别与平均球半径、等面积球半径、等距离球半径、等体积球半径这4种常用球体半径进行比较,研究表明地球向径积分平均值与4种常用球体半径间的差异更小。由于地球是一个旋转椭球体,向径与曲率半径是背离的,向径最大时,曲率半径最小, 向径最小时,曲率半径最大,传统思维所认为的曲率半径并不能准确地代表地球半径平均值,因此在一定程度上,地球向径的积分平均值更能代表地球半径平均值。这些研究结果可为地球科学、空间科学、导航定位提供基础理论依据。

English Abstract

宗敬文, 李厚朴, 钟业勋. 地球椭球向径和平均曲率半径的积分表达式[J]. 武汉大学学报 ● 信息科学版, 2022, 47(7): 1063-1070. doi: 10.13203/j.whugis20200059
引用本文: 宗敬文, 李厚朴, 钟业勋. 地球椭球向径和平均曲率半径的积分表达式[J]. 武汉大学学报 ● 信息科学版, 2022, 47(7): 1063-1070. doi: 10.13203/j.whugis20200059
ZONG Jingwen, LI Houpu, ZHONG Yexun. Integral Expressions of Earth Ellipsoid Radius Vector and Mean Radius of Curvature[J]. Geomatics and Information Science of Wuhan University, 2022, 47(7): 1063-1070. doi: 10.13203/j.whugis20200059
Citation: ZONG Jingwen, LI Houpu, ZHONG Yexun. Integral Expressions of Earth Ellipsoid Radius Vector and Mean Radius of Curvature[J]. Geomatics and Information Science of Wuhan University, 2022, 47(7): 1063-1070. doi: 10.13203/j.whugis20200059
  • 地球向径和平均曲率半径是测量和地球科学计算中常用的基本参数,根据地球科学、空间科学的要求和某些需求,常用到地球向径、平均曲率半径、平均球半径、等距离球半径、等面积球半径和等体积球半径[1-7]。由于地球是一个旋转椭球,向径与曲率半径是背离的。向径ρ最大(ρ=aa为椭球长半轴)时,曲率半径最小;向径ρ最小(ρ=bb为椭球短半轴)时,曲率半径最大。只有将地球看作是球体时,向径与半径才相等。

    随着空间技术和计算机技术在大地测量及地图学[8]中的应用和发展,研究地球向径和常用地球半径之间的关系具有更加重要的实用价值。文献[9]提出使用不同的球体半径表示球体特性,将不同球体半径代入克拉索夫斯基椭圆参数求得数值解,再代入到球体面积公式,求得不同半径下大圆航线长度和角度变形大小,发现使用等角球半径计算大圆航线可以满足航海使用精度要求。文献[10]将子午圈曲率半径和卯酉圈曲率半径展开为大地纬度的幂级数,并借助辅助函数V导出了不同形式的平均曲率半径展开式。文献[11]利用地球向径和计算机代数系统对大地测量学和地理信息系统中的问题进行了研究。

    从目前来看,前人对这一领域做了很多卓有成效的工作,但主要是针对其中一种或几种地球半径单独进行计算和应用,多是对等角纬度、等面积纬度及等距离纬度的计算与分析[12-17],很少有文献对向径进行推导和研究。为丰富对这一问题的理论研究,引入向径积分平均值和半径积分平均值的概念,并利用计算机代数系统(计算机代数系统可以将基本的数学公式展开成幂级数形式,推导出的公式相对于人工计算具有更高的精确度[18-21])将向径和半径分别进行椭圆曲线积分和椭球曲面积分,得到向径曲线积分平均值、向径曲面积分平均值、曲率半径曲线积分平均值和曲率半径曲面积分平均值。

    • 地球向径ρ是地心O到椭球面上任意一点P的距离,如图 1所示。图 1中,X为横轴,Y为纵轴。ρ关于地心纬度的表达式为:

      ρ=b1-e2cos2ϕ ]]>

      图  1  地球向径示意图

      Figure 1.  Digram of Earth Radius Vector

      式中,ϕ为地心纬度,地心纬度为参考椭球上一观测点和椭球中心的连线与赤道平面之间的夹角;e为椭球偏心率;b为地球短半轴。地心纬度ϕ与大地纬度B的关系表达式为:

      tanϕ=(1-e2)tanB ]]>

      略去推导过程,由式(1)和式(2)得到地球向径ρ与大地纬度B的关系表达式为:

      ρ=1-e2(2-e2)(1-cos2B)1-e2(1-cos2B) ]]>

      将式(3)展开成偏心率e的幂级数形式:

      ρ=a(a0+a1cos2B+a2cos4B+a3cos6B+a4cos8B) ]]>

      式中,a为地球长半轴,系数a0~a4为:

      a0=1-12e2-18e4-116e6-5128e8a1=12e2+34e4-516e6+3532e8a2=-58e4-2716e6-20364e8a3=1316e6+10332e8a4=-141128e8 ]]>

      平均曲率半径Ra是描述椭球面曲率的几何量,即椭球面上某一点所有方向法截线曲率半径的算术平均值,常用于制作地球表面上局部地区地图[22],如图 2所示。

      图  2  平均曲率半径示意图

      Figure 2.  Diagram of Average Curvature Radius

      Ra的表达式为:

      Ra=MNN=a(1-e2sin2B)1/2M=a(1-e2)(1-e2sin2B)3/2 ]]>

      式中,N为卯酉圈曲率半径;M为子午圈曲率半径。略去推导过程,得到平均曲率半径与大地纬度B的关系表达式为:

      Ra=a(1-e2)1-e2(1-cos2B)/2 ]]>

      将式(7)展开成偏心率e的幂级数形式:

      Ra=ab1+b2cos2B+b3cos22B+b4cos32B+b5cos42B ]]>

      式(8)的系数b1~b5为:

      b1=1-18e4-18e6-13128e8b2=-12e2-14e4-116e6+132e8b3=14e4-14e6+532e8b4=-18e6-316e8b5=116e8 ]]>

      图 1中可以看出,地球向径ρ是地心到椭球面上任意一点的距离。由于曲率半径是随着大地纬度B不断变化的,曲率半径不过地球球心,因此所构成的曲率圆也不同。于是,随着点位不断的变化,曲率圆中心形成了图 2所示的渐屈线。也就是说,曲率半径的中心是在不断变化的,并不是从地心到地球表面的一段线段,不能准确地表示地球半径的准确位置。从这一点来看,地球向径可以更精准地描绘地心到地球上某一点的距离,更具有代表性。

      图 3为地球向径、平均曲率半径变化示意图。从图 3中可以看出,地球向径随着大地纬度的增大而减小,平均曲率半径随着大地纬度的增大而增大,地球向径与平均曲率半径在一定意义上是背离的。向径最大时,曲率半径最小;向径最小时,曲率半径最大。

      图  3  地球向径、平均曲率半径变化示意图

      Figure 3.  Diagram of Changes About Earth Radial Direction and Average Radius of Curvature

    • 基于地球向径和平均曲率半径关于大地纬度B的表达式,本文提出一个新的概念,即地球向径的椭圆曲线积分平均值和平均曲率半径的椭圆曲线积分平均值。以推导地球向径的椭圆曲线积分平均值为例:首先对地球向径关于大地纬度B的表达式进行椭圆曲线积分,再将经过椭圆曲线积分后的地球向径表达式取平均值,得到平面椭圆上的地球向径曲线积分平均值符号表达式。

      将式(3)展开成幂级数形式,并在[0,π/2]范围内积分后取平均得到地球向径曲线积分平均值级数解表达式:

      $$ \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)

      式中,ρ¯为地球向径曲线积分平均值。

      同理,对平均曲率半径关于大地纬度B的表达式进行椭圆曲线积分,再将经过椭圆曲线积分后的平均曲率半径表达式取平均,得到平面椭圆上的平均曲率半径曲线积分平均值符号表达式。

      将式(7)在[0,π/2]范围内进行积分计算,得到平均曲率半径在平面椭圆上的曲线积分平均值:

      $$ \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)

      式中,R¯a为平均曲率半径曲线积分平均值。

      从式(11)中可以看出,平均曲率半径曲线积分平均值为a,要大于一般意义上的地球半径曲线平均值(即地球半径平均值为(a+b)/2),与通常所认为的平均曲率半径在平面椭圆上的平均值应该在椭圆长半轴a和短半轴b之间不同,证明存在系统误差。也就是说,平均曲率半径代表地球半径平均值是不精确的,这与平均曲率半径不在地心有关,然而地球向径是从地心到椭球面上任意一点的距离,将式(10)和式(11)进行对比可以看出,地球向径曲线积分平均值在某种程度上更能代表地球半径平均值进行测量和地球科学计算。

    • 基于上述结果,本文提出另一个新的概念,即地球向径在地球椭球面上的曲面积分平均值和平均曲率半径在地球椭球面上的曲面积分平均值。将地球向径关于大地纬度B的表达式进行椭球面曲面积分后再取平均,得到地球向径在地球椭球面上的曲面积分平均值符号表达式。将式(3)代入到椭球面曲面积分公式中可得:

      ρ¯¯=ρdσdσ=02π-π2-π2a1-e2(2-e2)sin2B1-e2sin2Bb2cosB(1-e2sin2B)2dBdL2πb2(arctanh(e)e+11-e2) ]]>

      式中,ρ¯¯表示地球向径在椭球面上的曲面积分平均值;b2cosB(1-e2sin2B)2dBdL表示球面积分dσdσ=2πb2(arctanh(e)e+11-e2)为椭球面积的分析解。

      由于此积分公式是椭圆积分,得不到分析解,因此将式(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)

      同理,将平均曲率半径关于大地纬度B的表达式进行椭球面曲面积分后再取平均,得到平均曲率半径在地球椭球面上的曲面积分平均值符号表达式。

      将式(7)代入到椭球面曲面积分公式中得:

      R¯¯a=Radσdσ=02π-π/2π/2a1-e2(1-e2sin2B)b2cosB(1-e2sin2B)2dBdL2πb2(arctanh(e)e+11-e2)= a(1-e2)(3arctanh(e)e+5-3e2(1-e2)2)8π(arctanh(e)e+11-e2) ]]>

      式中,R¯¯a表示平均曲率半径在椭球曲面上的曲面积分平均值。

      为便于将平均曲率半径和地球向径的椭球面曲面积分平均值进行比较,将式(7)代入到椭球面曲面积分公式中先展开成e的幂级数形式再进行积分,得到平均曲率半径椭球面曲面积分平均值的级数解:

      $$ \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)作差得到:

      R¯¯a-ρ¯¯=a(215e4+44315e6+2 50318 900e8) ]]>

      由式(16)可以看出,在地球椭球面上,地球向径和平均曲率半径的曲面积分平均值差异很小,差异首项为215e4。由于曲率半径的中心是在不断变化的,并不能准确地表示地球半径的准确位置,从这一点来看,地球向径更加精准。

    • 就地球总体而言,平均球半径(取地球椭球三轴半轴长的算术平均值,用于简单决定球半径) [23]、等面积球半径(保持球体表面积等于地球椭球面相应全面积而决定的球半径,主要用于等面积投影) [24]、等距离球半径(使球面经线总长等于地球椭球面经线总长而决定的球半径,主要用于等距离投影中)[25]和等体积球半径(使地球球体的体积等于地球椭球体的体积来决定的球半径)的表达式与大地纬度B无关,只取决于地球参考椭球模型的参数ae,其表达式分别为:

      Re=a(2+1-e2)/3RF=a12+1-e24eln1+e1-eRS=a(1-14e2-364e4-5256e6-17516 384e8)RV=a(1-e2)1/6 ]]>

      式中,Re表示平均球半径;RF表示等面积球半径;RS表示等距离球半径;RV表示等体积球半径。将4种常用球体半径分别与地球向径在椭球面上的曲面积分平均值和平均曲率半径在椭球面上的曲面积分平均值作差,Δ1=a(Re-ρ¯¯),Δ2=a (RF- ρ¯¯), Δ3 = a (RS - ρ¯¯), Δ4 = a (RV - ρ¯¯)D1=a(Re-R¯¯a), D2=a(RF-R¯¯a), D3=a(RS-R¯¯a),D4=a(RV-R¯¯a),推导出其差异符号表达式列于表 1表 2

      表 1  4种常用球体半径与向径平均值差异符号表达式

      Table 1.  Different Symbolical Expressions Between Four Global Radii and Integral Mean Values ofEllipsoid Radius Vector

      Δi 差异符号表达式
      Δ1 a(1180e4+5756e6-499 099907 200e8)
      Δ2 a(1189e6-2 2274 050e8)
      Δ3 a(-112e2+12 880e4+38348 384e6-127 225 969232 243 200e8)
      Δ4 a(-145e4-171 134e6-192 959340 200 e8)

      表 2  4种常用球体半径分别与曲率半径平均值差异符号表达式

      Table 2.  Different Symbolical Expressions Between Four Global Radii and Mean Radius of Curvature

      Di 差异符号表达式
      D1 a(-23180e4-5033 780e6-619 243907 200e8)
      D2 a(-215e4-127945e6-38 68756 700e8)
      D3 a(-112e2-3832 880e4-31 877241 920e6-157 982 833232 243 200e8)
      D4 a(-745e4-8775 670e6-238 013340 200e8)

      对比表 1表 2可知,4种常用球体半径与向径曲面积分平均值之差比4种常用球体半径与平均曲率半径曲面积分平均值之差要小。平均球半径与向径曲面积分平均值之差首项系数绝对值是平均球半径与平均曲率半径曲面积分平均值之差的1/23;等面积球半径与向径曲面积分平均值之差比等面积球半径与平均曲率半径曲面积分平均值之差多了一个数量级e4;等距离球半径与向径曲面积分平均值之差第二项系数绝对值是等距离球半径与平均曲率半径曲面积分平均值之差的1/383;等体积球半径与向径曲面积分平均值之差第二项系数绝对值是等体积球半径与平均曲率半径曲面积分平均值之差的1/7

      本文以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
      Δi 1.59 0.003 -3 556.36 -6.38
      Di -36.75 -38.34 -3 594.71 -44.73

      观察表 3,通过数值比较看出4种常用球体半径与向径曲面积分平均值之差比4种常用球体半径与平均曲率半径曲面积分平均值之差要小,符合表 1表 2结论。为了更直观地表现出4种常用球体半径与向径曲面积分平均值之差的差异问题,以及4种常用球体半径与平均曲率半径曲面积分平均值之间的差异问题,将椭球体的偏心率进行变化,将椭球偏心率的范围确定为e(0, 0.8),以便更好地放大它们之间的差异关系并分析讨论,差异曲线见图 4~7

      图  4  Δ1D1线

      Figure 4.  Change Graph of Δ1,D1

      图  5  Δ2D2线

      Figure 5.  Change Graph of Δ2,D2

      图  6  Δ3D3线

      Figure 6.  Change Graph of Δ3,D3

      图  7  Δ4D4线

      Figure 7.  Change Graph of Δ4D4

      图 4~7可以看出,椭球偏心率e在(0,0.8)范围内变化时,4种常用球体半径与向径曲面积分平均值之差比4种常用球体半径与平均曲率半径曲面积分平均值之差小,差值更加稳定,反映出地球向径曲面积分平均值更接近4种常用球体半径。

      将地球向径在椭球面上的曲面积分平均值、平均曲率半径在地球椭球面上的曲面积分平均值、平均球半径、等面积球半径、等距离球半径和等体积球半径随偏心率e在(0,0.8)的变化曲线绘制如图 8所示。由图 8可知,在椭球偏心率e在(0,0.8)范围内变化时,地球向径曲面积分平均值与4种常用球体半径之间更为接近,特别是在偏心率e逐渐变大的过程中体现得更为突出。因此,地球向径积分平均值比平均曲率半径积分平均值能更加准确地代表地球半径平均值进行测量和地球科学计算。

      图  8  变化曲线图

      Figure 8.  Change Graph of R¯¯a,ρ¯¯,Re,RF,RS,RV

    • 本文从地球向径和平均曲率半径的基本概念入手,引入了两个新的概念:地球向径与平均曲率半径的曲线积分平均值和曲面积分平均值,通过计算机代数系统对它们进行推导计算,并进行对比分析,为大地测量基准问题提供理论参考,得到如下结论:

      1)地球是一个旋转椭球体,地球向径与平均曲率半径在一定意义上是背离的:向径最大(ρ=a)时,平均曲率半径最小;向径最小(ρ=b)时,平均曲率半径最大;只有将地球看作是球体时,向径与平均曲率半径相等。这表明平均曲率半径不能精准地代表地球半径。

      2)地球向径是从地心到椭球面上任意一点的距离,地球向径曲线积分平均值介于椭球长半轴a与短半轴b之间,相比较于平均曲率半径曲线积分平均值,可以更加准确地代表地球半径的平均值进行测量和地球科学计算。

      3)地球向径曲面积分平均值与4种常用球体半径之差比平均曲率半径曲面积分平均值与4种常用球体半径之差小,特别是偏心率e在(0,0.8)范围内变化时显得尤为明显。从这一点来看,地球向径积分平均值能更加准确地代表地球半径平均值进行测量和地球科学计算。

      4)将地球向径与常用地球半径的差异表示为符号形式,并统一展开为偏心率e的幂级数形式。该表达式易于比较分析,一定程度上丰富了测量及地图学数学分析理论。

参考文献 (25)

目录

    /

    返回文章
    返回