子午线弧长的一种新算法应用.pdfVIP

第 23 卷第 3 期 海 　洋 　测 　绘 Vol 23 ,N o3 2003 年 5 月 HYDRO GRA PHIC SURV EYIN G AND CHARTIN G May ,2003 子午线弧长的一种新算法应用 严 伯 铎 ( 山东省城乡建设勘察院 ,山东 济南 　250031) A New Algorithm f or Arc L ength of Meridan YAN B oduo 　　摘要 : 介绍一种椭球子午线弧长的计算方法 ,并以我国常用的 1954 年北京坐标系椭球为例进行了实例计算 。 该法与椭圆积分计算方法相比较 ,公式简单明了 ,计算方法简便 ,并能满足精度要求 。 关键词 : 子午线弧长 ;计算方法 ;精度分析 ( ) 中图分类号 : P226 　　　文献标识码 : B 　　　文章编号 : 2003 α ( ) 2 β B ( ) 2 1 　引　言 则 : = A a 1 - e ; 　 = - a 1 - e ; 2 γ C ( ) 2 δ D ( ) 2 ( ) 在处理大地测量 、天文测量 、航天航空技术以及 = - a 1 - e ; 　 = - a 1 - e 4 2 2 地理信息系统环境中的某些几何问题时 ,椭球的子 午线弧长的计算是经常要进行的一项工作 。关于子 在这些计算式中都是以椭球的第一偏心率 e 作 午线弧长的计算式 ,通常都是把以椭圆积分的表达 为椭球的形状参数来表述的。为了保证子午线弧长 ( ) ( ) 式 ,依二项式展开成级数形式后 ,再以逐项积分的方 的计算达到必要的精度, 通常对 3 或 4 式的各系 6 8 ( ) ( ) 数的计算都要取位到 e 甚至 e 以上 。而对 1 式 法 ,在保证满足必要的精度 通常为 0001m 要求下 [ 1 ,2 ] ( ) 给出实用的计算式 即: