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

高斯投影坐标正反算编程报告 编程思想 进行高斯投影坐标正反算的编程需要牵涉到大量的公式,为了使程序条理更清楚,各块的数据复用性更强,这个地方采取了结构化的编程思想。 程序由四大块组成。 GeodｅsyＨoｍｅworｋ、cpp文件用于存放main()函数,是整个程序的入口、通过结构化的编程尽力使main()函数变得简单、 ＭyFunｃｔioｎ、h和MyFuｎctｉon、cｐp用于存放计算过程中进行角度弧度换算时所要用到的一些自定的转换函数。 Zheｎgsuan、h和Zheｎgsuan、cpｐ用于存放Zｈｅｎｇsｕａｎ类,在Zheｎgsuａｎ类中声明了高斯投影坐标正算所要用到的所有变量,在类的构造函数中进行成员变量的初始化及正算计算、通过ｇeｔ函数获得相应的正算结果、 Ｆaｎsuan。ｈ和Fansuａｎ、cpp用于存放Fansｕａn类,类似于Zhengsuan类,Faｎｓuan类中声明了高斯投影坐标反算所要用到的所有变量,在类的构造函数中进行成员变量的初始化及反算计算。通过geｔ函数获得相应的反算结果。 计算模型 高斯投影正算公式 高斯投影反算公式 　 　 程序框图 开始 开始 输入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 结束 正算 反算 计算结果 附录:程序代码 ///／／主函数入口 GｅｏdｅsｙHomework、cpp #inclｕde　”MyFunction、h ＃iｎclｕde Zhengｓuan、ｈ ＃inｃludｅ Faｎｓuan、h #include 〈iostrｅam usｉng nａmｅｓpace sｔd; ｖoｉd　fanｓuａn(); ｖoid zｈenｇsuaｎ(); void main(){ zｈｅngsuan(); fａｎsuaｎ(); pｒintf(”/n　oｖer!”); } voｉd zｈeｎgsuan()｛ ?dｏｕbｌｅ　ｍyＢ,mｙL; coｕt＜〈【正算】〈enｄl; ?cout请输入大地纬度B＂＜endｌ; myB=angｌeＴoDegrｅｅ(); coｕｔ＜＜”请输入大地经度L＜＜enｄl; ?mｙL=ａnｇｌeToDegｒee(); Zhenｇsuan myZhengｓuaｎ1(mｙB,myL); ?ｐriｎｔｆ(”Raｄiａn B=％ｆ L=%f \ｎ＂,ｍyＺheｎgsuan１、getrB(),myＺｈenｇsuan1。getrL()); ?myＺheｎｇsｕａn1。pｒｉntLｏcａｔｉon(); } voiｄ　faｎsuａn(){ double ｍyＸ,ｍyY; ?cｏut【反算】〈endl; ?cout＜”请输入国家统一坐标　X　Y、例如 3３78627、18１9４517ｅndl; cin〉＞myX＞＞mｙＹ; ?Faｎsuan　myFａnsuan１(ｍｙＸ,myＹ); ?ｍyFansｕan１。printLoｃａtiｏn(); ｝ ///自定功能函数库 MｙFunction、ｈ #defｉne PＩ 3、1415926 #include iｏsｔrｅam〉 using　namespａce std; doublｅ angｌeTｏDegｒee(ｉnｔ ｄu,int feｎ,ｆｌoaｔ　mｉao); ｄouｂｌe angleToＤegrｅe();//将度分秒换算为度 doubｌｅ　deｇreeTｏRadｉaｎ(dｏuble degrｅe); doｕｂlｅ dｅgｒeｅTｏRadｉan();／/将角度换算为弧度 MyFuncｔion、ｃpp #ｉnclude ＂MyFunctｉoｎ、h dｏuble ａnｇlｅTｏＤegree(iｎt du,int fen,ｆlｏａｔ ｍiao){ double　resuｌt=0; ?reｓｕlt=miao／3600、0＋fen/60。0＋du; reｔuｒn ｒｅsult; } ｄoubｌe aｎgleToDｅgｒeｅ()｛ ?inｔ du,ｆeｎ; ?fｌoat ｍiaｏ; ?dｏｕble　reｓｕｌt; ?ｃｏｕt＂请输入度　分 秒。例如:３0 ２0 0０〈eｎdｌ; ?cin＞＞dufeｎ〉〉ｍiao; reｓｕｌt＝aｎgleTｏDegree(du,fen,miaｏ); ?rｅturn　ｒｅsult; } double ｄegrｅｅToRadian(doｕbｌe deｇreｅ)｛ ｄｏuble result＝0; ｒesult＝ｄｅgree/5７。295７７９５1308２