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

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

2、过得到函数获得相应的正算结果。fansuan.h和fansuan.cpp用于存放fansuan类，类似于正义类fansuan类中声明了高斯投影坐标反算所要用到的所有变量，在类的构造函数中进行成员变量的初始化及反算计算。通过得到函数获得相应的反算结果。2.计算模型高斯投影正算公式高斯投影反算公式3.程序框图开始输入b，l求定带号n、中央纬度l0，纬度差l按照实用公式计算x，y换算为国家统一坐标x，y输出x，y输入国家统一坐标x，y由y取定带号n、并换算出x，y求出中央经线劳拉西泮按照实用公式计算b，ll=l0 l求出大地经度l输出b，l结束正算反算4.计算结果5.附录：程序代码/主函数入口geo

3、desyhomework.cpp#包含 myfunction.h #包括正算h#包括fansuan.h #包括使用命名空间标准；虚空扇(suan)；无效正算()；void main()正义()；范suan()；printf(/n)结束！”)；void zhengsuan()双myb，myl【cout】正算】我的圣诞节；范suan范苏安1(myx，myy)；我的安素1号。打印位置()；/自定功能函数库myfunction.h#定义pi 3.1415926#包括使用命名空间标准；双角度(整数，整数，浮点苗)；双角度度()；/将度分秒换算为度双度数弧度（双度数)；双度数弧度()；/将角度换算为弧度m

4、yfunction.cpp#包含 myfunction.h 双角度度(整数，整数，浮点苗)双重结果=0；结果=苗/3600.0分/60.0度；返回结果；双角度度()国际分；浮苗；双重结果；cout 请输入度分秒。例如：30 20 00 dufenmiao结果=角度度（度、分、苗)；返回结果；双度数弧度（双度数)双重结果=0；结果=度/57.295779513082321；返回结果；双度数弧度()双重结果，程度；度=角度度()；结果=度数弧度（度)；返回结果；/正算类zhengsuan.h/正义类的正义。h:接口。/#如果！已定义(afx _ zonsuan _ h _ 2655 ea28 _

5、e810 _ 44a 3 _ 8f 14 _ 56421 a7 b 4466 _ included _)#定义afx _ zonsuan _ h _ 2655 ea28 _ e810 _ 44a 3 _ 8f 14 _ 56421 a7 b 4466 _ included u#如果_msc_ver 1000#pragma一次#endif /_msc_ver 1000#定义rousecond 20060。888888888686#包含 myfunction.h #包括#包括使用命名空间标准；正算班public:正义()；正算（双fb，双fl)；双getx()；双gety()；双getrb()；双

6、getrl()；无效打印位置()；虚拟正算()；private:双x；双y ./大地坐标双十；双y。/国家统一坐标双乙；双rb。b秒内；双倍升；双rl。/输入的大地纬度b，大地经度l，rb，rl为对应弧度表示值b秒为换算成秒数值int n；/带号n双l0。/中央经线纬度劳拉西泮双ldot/纬度差l-l0int ldotsecond/换算成秒的纬度差双倍l .双氮；双a0。双a3 .双a4；双a5；双a6；/七个计算参数；#endif /！已定义(afx _ zonsuan _ h _ 2655 ea28 _ e810 _ 44a 3 _ 8f 14 _ 56421 a7 b 4466 _ in

7、cluded _)zhengsuan.cpp/正算类的实现。/#包括正算h/建设/破坏/正算：正算()正算：正算（双fb，双fl)b=fb .rb=度数弧度(瑞士法郎)；l=flrl=度数弧度(f1)；b秒=b * 3600/初始化大地经度我，大地纬度第二，按弧度的大地纬度元素铷的符号n=(整数)(1/6 1)；/初始化带号nl0=6 * n-3；/中央经线经度，角度单位ldot=l-l0；/经度差ldot秒=ldot * 3600l=(ldot)* 3600/rouscond；/计算参数ln=6399698.902-(21562.267-(108.973-0.612 * cos(rb)* c

8、os(rb)* cos(rb)* cos(rb)* cos(rb)* cos(rb)；/计算参数na0=32140.404-(135.3302-(0.7092-0.004 * cos(rb)* cos(rb)* cos(rb)* cos(rb)* cos(rb)* cos(rb)；/计算参数a0a4=(0.25 0.00252 * cos(rb)* cos(rb)* cos(rb)* cos(rb)-0.04166；/计算参数a4a6=(0.166 * cos(rb)* cos(rb)-0.084)* cos(rb)* cos(rb)；/计算参数a6a3=(0.333333333 0.0011

9、23 * cos(rb)* cos(rb)* cos(rb)* cos(rb)-0.1666667；/计算参数a3a5=0.0083-(0.1667-(0.1968 0.004 * cos(rb)* cos(rb)* cos(rb)* cos(rb)* cos(rb)* cos(rb)；/计算参数a5x=6367558.4969 * bsecond/rosecond-(a0-(0.5(a4 a6 * l * l)* l * l * n)* sin(rb)* cos(rb)；/正算xy=(1(a3 a5 * l * l)* l * l)* l * n * cos(rb)；/正算yx=x .y=n

10、*1000000 y 500000 ./国家统一坐标正算： 正算()双正数：getx()返回x。双正数：gety()返回y。无效字符：打印位置()printf()正算得国家统一坐标为：x=%8.8f y=%8.8f n ，x，y)；双正数：getrb()返回rb .双正数：getrl()返回rl。/反算类fansuan.h/fansuan。h : fansuan类的接口。/#如果！已定义(afx _ fansaun _ h _ 5 b5 e 4668 _ ee81 _ 463 f _ 8d 65 _ ffe 2 eace 745 _ included _)#定义afx _ fansuan _ h _ 5 b5 e 4668 _ ee81 _ 463 f _ 8d 65 _ ffe 2 eace 745 _ included u#如果_msc_ver 1000#pragma一次#endif /_msc_ver 1000#定义rousecond 20060。888888888686#包括#包含 my