高斯反算三种实用算法的比较研究.pdfVIP

第1期 矿 山 测 量 No．1 2008年3月 MINE SURVEYING MaL．20 高斯反算三种实用算法的比较研究 王彦伟’，李海洋 ，裴 亮 ，黄吉来 ，张恒壕 (1．葫芦岛市勘查测绘院，辽宁葫芦岛 125000；2．辽阳市国土资源勘查研究院，辽宁辽阳 ll1000； 3．辽宁工程技术大学测绘与地理科学学院，辽宁阜新 123000) 摘要：设计实现了三种简单的高斯反算算法，即一种直接算法和两种迭代算法。试验表明，三种算法 结果相同，均达到了0．00001秒的精度，可以代替此等精度比较复杂的反算算法，．具有实用价值。 关键词：高斯反算；迭代；直接算法 中图分类号：P209 文献标识码：B 文章编号：1001—358X(2oo8)Ol一0048—03 随着国民经济和城市建设的发展，各种工程建 算中极容易出现各种计算错误，从而影响了计算结 设对坐标的需要越来越多。。我国常用的有 1954年 果。而用计算机，只要是算法设计和程序正确，只需 北京坐标系、1980年国家大地坐标系等，他们都是高 几秒钟便可以计算完毕，即缩短了计算周期，提高了 斯平面直角坐标系。但是以前的坐标都是大地坐 效率，又可以避免各种计算错误，从而提高了计算结 标，因此，工程建设中就需要将他们两者进行转换。 果的可靠性。 这是一个旧问题，同时也是一个薪问题，如有研究用 论文中的高斯反算公式采用的是现在普遍使用 cASIOfx 4800P型可编程计算器进行经纬度、3。带和 的一般公式，在高斯反算中设计了3种方法来计算 6。带坐标转换应用程序实现了使用率较高的普遍编 L，即两种迭代和一种直接算法。 程计算器进行一般精度要求的高斯投影面坐标转换 1 算法设计 的计算…，也有研究高斯反算的 ，还有高斯反算用 于地图投影 的等等。 经过推导高斯反算公式 ： 一 以前，坐标系之间互相转化是通过人工计算来 实现的；人工计算不仅计算周期长，效率低，而且计 = 一 。+南4 叫+ 器 … 45 z 一 丽y3 (1+2 + )+ 5 ‘5+28 +24 +6 +8 2。2) 式中凡脚注有‘ ’的，表明这些函数符号都是以 的。 垂足纬度 ，代入求得的。而垂足纬度 ，可以根据子 在计算中设计了三种算法，它们是两种迭代和 午线弧长公式，由 = 很快反算求出。这样，根据式 一种直接算法 。以54坐标系为例介绍这三种算 (3)就可以计算大地纬度B和经差f，进而求出大地 法。 纬度 。式(3)中B、z均以弧度为单位表示，如果以度 第一种迭代首先令 ，=0；然后进入循环体，进 为单位，等号右端应乘以 。 行循环。 仃 … ．． 先计算 ： 当z3．5。时，按式(3)进行坐标换算的精度为 0．0001”，这与公式(I)的换算精度0．001in是相适应 =c×(B／(1一÷ +4 5e4：一丽175e6：+ )+((一 3 e：2+ 45e4：一丽175e6：+ 翳 )×c。s( )+( e4：一 6 e ：+ e：8)×c。s( ) +(一 35e6： + e8)×c。s( ) + e8：×c。s( ) ×sin( )))