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

高斯投影坐标正反算编程报告 编程思想 进行高斯投影坐标正反算的编程需要涉及到大批的公式，为了使程序条理更清楚， 各块的数据复用性更强，这里采纳了却构化的编程思想。 程序由四大块构成。 文件用于寄存main( )函数，是整个程序的进口。经过结构化的编程全力使main( ) 函数变得简单。 和用于寄存计算过程中进行角度弧度换算时所要用到的一些自定的变换函数。 和用于寄存Zhengsuan类，在Zhengsuan类中声了然高斯投影坐标正算所要用到的 全部变量，在类的结构函数中进行成员变量的初始化及正算计算。经过get函数获取相 应的正算结果。 和用于寄存Fansuan类，近似于Zhengsuan类，Fansuan类中声了然高斯投影坐标 反算所要用到的全部变量，在类的结构函数中进行成员变量的初始化及反算计算。经过 get函数获取相应的反算结果。 计算模型 高斯投影正算公式 x X N2sinBcosBl 2 N4 simBcos3B(5 t2 92 44)l4 2 24 N 6sinBcos5B(61 58t2 t4)l 6 720 y NcosBl N 3 cos3 B(1 t2 2)l 3 6 N 5cos5B(5 18t2 t4 14 2 58 2t2)l 5 120 高斯投影反算公式 BBf tf y 2 tf 2 2 9 2 2 y 4 2Mf 53tf f f tf Nf 24MfN3f tf Nf5y61 90tf2 45t4f y6 720Mf y y3 1 2 2 l 6N 2tf f NfcosBf 3fcosBf y5 5 28t2f 24tf4 62f 8 2ft2f 120N5f cosBf 程序框图 开始 输入B，L 求定带号N，中央纬度L0，纬度差l 正 依据适用公式计算x,y 换算为国家一致坐标X，Y 输出X,Y 输入国家一致坐标X,Y 由Y取定带号N，并换算出x，y 求出中央经线L0 反 依据适用公式计算B，l L=L0+l求出大地经度L 输出B，L 结束 计算结果 附录：程序代码 517endl; cinmyXmyY; FansuanmyFansuan1(myX,myY);( ); } // ////////////////////////////////////////////////////////////////////// #include ////////////////////////////////////////////////////////////////////// Construction/Destruction ////////////////////////////////////////////////////////////////////// Fansuan::Fansuan( ) { } Fansuan::Fansuan(doubleX,doubleY){ this-X=X; this-Y=Y;//初始化x,y N=(int)(Y/1000000);//拿出带号 L0=6*N-3;//初始化该带号的中央经线经度 x=X; y=Y-1000000*N-500000; beta=x/;//初始化beta，弧度单位 betaSecond=beta*rouSecond;//初始化beta，秒单位 betaDegree=betaSecond/3600;//初始化beta,整度数单位 Bf=beta+(+(293622+(2350+22*cos(beta)*cos(beta))*cos(beta)*cos(beta))*co s(beta)*cos(beta))*(1e-10)*sin(beta)*cos(beta);//初始化Bf,弧度单位 BfSecond=Bf*rouSecond;//初始化Bf，秒单位 BfDegree=BfSecond/3600;//初始化Bf，整度数单位 Nf=Z=y/(Nf*cos(Bf)); b2=+*cos(Bf)*cos(Bf))*sin(Bf)*cos(Bf); b4=++*cos(Bf)*cos(Bf))*cos(Bf)*cos(Bf); b3=b5=Bsecond=BfSecond-(1-*Z*Z)*Z*Z)*Z*Z*b2*rouSecond;//计算大地经度 B，单位为秒 B=Bsecond/3600;//用整度数表示B lsecond=(1-(b3-b5*Z*Z)*Z*Z)*Z*rouSecond;//计算经度差l，单位为秒 l=lsecond/3600;//用整度数表示l L=L0+l;//计算大地经度L } doubleFansuan::getB( ){ returnB; } doubleFansuan::getL( ){