高斯投影坐标正反算编程报告.doc

《高斯投影坐标正反算编程报告.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、高斯投影坐标正反算编程报告1.编程思想进行高斯投影坐标正反算的编程需要牵涉到大量的公式，为了使程序条理更清楚，各块的数据复用性更强，这里采取了结构化的编程思想。程序由四大块组成。GeodesyHomework.cpp文件用于存放main()函数，是整个程序的入口。通过结构化的编程尽力使main()函数变得简单。MyFunction.h和MyFunction.cpp用于存放计算过程中进行角度弧度换算时所要用到的一些自定的转换函数。Zhengsuan.h和Zhengsuan.cpp用于存放Zhengsuan类，在Zhengsuan类中声明

2、了高斯投影坐标正算所要用到的所有变量，在类的构造函数中进行成员变量的初始化及正算计算。通过get函数获得相应的正算结果。Fansuan.h和Fansuan.cpp用于存放Fansuan类，类似于Zhengsuan类，Fansuan类中声明了高斯投影坐标反算所要用到的所有变量，在类的构造函数中进行成员变量的初始化及反算计算。通过get函数获得相应的反算结果。2.计算模型高斯投影正算公式高斯投影反算公式3.程序框图11开始输入B，L求定带号N，中央纬度L0，纬度差l按照实用公式计算x,y换算为国家统一坐标X，Y输出X,Y输入国家统一坐标X

3、,Y由Y取定带号N，并换算出x，y求出中央经线L0按照实用公式计算B，lL=L0+l求出大地经度L输出B，L结束正算反算1.计算结果111.附录：程序代码/////主函数入口GeodesyHomework.cpp#include"MyFunction.h"#include"Zhengsuan.h"#include"Fansuan.h"#includeusingnamespacestd;voidfansuan();voidzhengsuan();voidmain(){zhengsuan();fansuan();pri

4、ntf("/nover!");}voidzhengsuan(){doublemyB,myL;cout<<"【正算】"<

5、ocation();11}voidfansuan(){doublemyX,myY;cout<<"【反算】"<>myX>>myY;FansuanmyFansuan1(myX,myY);myFansuan1.printLocation();}///自定功能函数库MyFunction.h#definePI3.1415926#includeusingnamespacestd;doubleangle

6、ToDegree(intdu,intfen,floatmiao);doubleangleToDegree();//将度分秒换算为度doubledegreeToRadian(doubledegree);doubledegreeToRadian();//将角度换算为弧度MyFunction.cpp#include"MyFunction.h"doubleangleToDegree(intdu,intfen,floatmiao){doubleresult=0;result=miao/3600.0+fen/60.0+du;returnresul

7、t;}doubleangleToDegree(){intdu,fen;floatmiao;doubleresult;cout<<"请输入度分秒。例如：302000"<>du>>fen>>miao;result=angleToDegree(du,fen,miao);returnresult;}doubledegreeToRadian(doubledegree){doubleresult=0;result=degree/57.295779513082321;11returnresult;}doubledegreeToR

8、adian(){doubleresult,degree;degree=angleToDegree();result=degreeToRadian(degree);returnresult;}///正算类Zhengsuan