一组高精度椭球面电_子计_算实用公式.docVIP

一组高精度椭球面电子计算实用公式顾旦生(中国测绘科学研究院 100039)我国天文大地网整体平差于 1982 年 5 月完成。 经整体平差的我国天文大地网的精度得到极大提高一组高精度椭球面电子计算实用公式顾旦生(中国测绘科学研究院 100039)我国天文大地网整体平差于 1982 年 5 月完成。 经整体平差的我国天文大地网的精度得到极大提高, 能 完全满足我国国家级各种比例尺测图控制以及其它国民经济建设和国防建设等有关方面的需要。近年来, 由 于全国各地普遍开展大面积的大比例尺测图和各种大型工程建设, 因此, 整体平差成果已被许多部门重视和 应用。本文笔者长期从事超大规模大地网平差的研究, 是国家测绘局系统全国天文大地网整体平差的技术负 责, 也是整体平差方法和整套计算公式的研究和提供 者, 因而常被询问或索取有关整体平差计算公式。该组公式至今仍不失其实用价值, 故现将其发表供有关同志参考。我国天文大地网整体平差是在中国 1980 大地坐 标系 (建议简称 80 坐标系) 的参考椭球面上进行的。在椭球面上进行生产性特大规模大地网 ( 近 5 万个大地 点) 平差, 在国内外是首次。 椭球面大地测量计算公式一般采用级数形式, 这就涉及到公式的精度问题。我国 天文大地网是一个大自由网, 当采用坐标间接平差法进行整体平差时, 法方程的阶数很高, 近 20 万阶, 法方 程的条件数很大, 解算结果受计算凑整误差和公式误差的影响十分敏感, 特别是观测误差方程系数和常数项误差的影响最为重要。另外, 由于对极其复杂和庞大 的整体平差工作要求全部在电子计算机上自动完成, 因此, 在整体平差时, 除采用足够字长的计算机外, 必 须采用一组顾及所有特殊情况的、适用于电子计算机 自动运算的特高精度实用计算公式。 这组公式是笔者 于 1977 年在上海从事本项研究时提出的, 后根据生产 部门的意见, 于 1981 年笔者对观测误差方程系数和常 数项的计算公式又作了扩充, 进一步提高了公式的精 度。这组公式的一部分是取自有关文献, 但为了适应电 子计算自动判别和运算的需要, 在形式和算法等方面 作了必要的改变和补充; 另一部分公式是笔者当时新 推导而得。由于整体平差是采用以坐标间接平差为基础的条 件联系平差法1 、2 , 平差时除了输入必不可少的各种观测值和表示大地网结构的信息外, 还输入了全网各大地点的近似坐标值。鉴于这一情况, 为了尽量减少输 入数据和信息量, 笔者提供的下面各种有关计算公式的系数, 全部采用以大地点近似坐标为参数的中纬度 形式。这组公式的推导思想和过程十分繁杂, 本文只能给出其结果性的最后实用计算公式。 这组公式经过大 规模整体平差的使用和检验, 已充分说明其正确、可靠和有效。80 坐标系3 、54 坐标系3 和 8054 坐标系3 采用的椭 球参数长半径 a 和扁率 Α如下:80 坐标系:a = 63 78140 mΑ= 1?298. 25754 坐标系和 8054 坐标系:a = 6 378 245 mΑ= 1?298. 3一、大地坐标 B 、L 和高斯平面坐标 x、y 之 间的换算公式符号和常数80 坐标系:a = 6 378 140 mΑ= 1?298. 257b= 6 356 755. 288 157 528 6 m2e = 0. 006 694 384 999 587 952(1)2e′= 0. 006 739 501 819 472 927c= 6 399 596. 651 988 010 4 mΘ″= 206 264. 806 247 1Π= 3. 141 592 653 589 8t= tgB2 2 2Γ = e′co s BV = 1+ Γ22 1(2)1+ e′2V 0 =2N = c?V (m )M = c?V 3 (m )R = c?V 2 (m )c0 = 6 367 452. 132 79 mc1 = 32 009. 857 5 mc2 = 133. 960 2 mc3 = 0. 697 6 mK 0 = 0. 157 048 687 473×10- 6 m - 1K 1 = 0. 005 052 505 593K 2 = 0. 000 029 847 335K 3 = 0. c0 = 6 367 452. 132 79 mc1 = 32 009. 857 5 mc2 = 133. 960 2 mc3 = 0. 697 6 mK 0 = 0. 157 048 687 473×10- 6 m - 1K 1 = 0. 005 052 505 593K 2 = 0. 000 029 847 335K 3 = 0.