《高斯投影坐标正反算编程报告》由会员分享，可在线阅读，更多相关《高斯投影坐标正反算编程报告（13页珍藏版）》请在金锄头文库上搜索。
1、高斯投影坐标正反算编程报告1 .编程思想进行高斯投影坐标正反算的编程需要牵涉到大量的公式，为了使程序条理更清楚， 各块的数据食用性更强，这里采取了结构化的编程思想。程序由四大块组成。GcodcsyHomork.cpp文件用于存放main。函数，是整个程序的入口。通过结构化 的编程尽力使main。函数变得简单。MyFunction.h和MyFunclion.cpp用于存放计算过程中进行角度弧度换算时所要用到 的一些自定的转换函数。Zhciigsujin.h 和 Zhciigsuaii.cpp 用于存放 Zhciigsium 类，在 Zhcngsuan 类中声明了高 斯投影坐标正算所要用到的所有变
2、量，在类的构造函数中进行成员变量的初始化及正算 计算，通过get函数获得相应的正算结果。Fansuan.h 和 Fansiuui.cpp 用于存放 Fansuan 类，类似于 Zhcngsuan 类，Fansuan 类中 声明了高斯投影坐标反算所要用到的所有变量，在类的构造函数中进行成员变量的初始 化及反算计算。通过get函数获得相应的反算结果。2 .计算模型高斯投影正算公式x = X +、sin BcosB二 +-1r( 5/n?Bcos3 B(5-r2 +97/2 +44)/4+ sinBcos5 8(61-58/2+力/67201y = cosB - r + - -cos3 B(l-r2
3、 +z/2)/w3. P”6pN尸cos 8(5-18产 + +142 -582/)/5高斯投影反算公式y- +一匚(5 + 3片+#-9疝汾4 24M720M,N；yy(61 + 90r； +45，卜6N,cosB, 6Nos8,JJJJy5+1207VjcosBrJJ(5 + 28/； +24/； +6； +8)；)3 .程序框图Ha输出X,Y /输入国家统一坐标X,Y结束4.计算结果5.附录：程序代码/主函数入口GeodesyHomeworkcpp#includc HMvFunction.hH#includc HZhcngsu;in.hH#includc Fansuan.h#includ
4、c using namespace std;void fansuaiiQ;void zhcligsuaiiQ;void main。 zhciigsuruiQ; faiisurinO； printf(n/n ovcr!,r);void zhciigsuaiiQ double myB,myL;coutvv【正算】VVundl;coutvv”请输入大地纬度B,cndl；myB=r)cgrccO；coutve”请输入大地经度LMendl;myL=anglcTr)cgrcc0;Zhcngsuaii myZhciigsuuil (myB,myL);printffRadian B=%f L=%f nM,my
5、Zhengsuanl .getrBQ,myZhengsuanl .getrL0);myZhcngsuanl .printLocationQ;void fansuan()double myXjnyY;coutM【反算】“ vvendl;coutv”请输入国家统一坐标 XY。例如 3378627.1819 20243953.4517”mvX mvY;Fansuan myFansuanl(myX,myY);myFansuaiil.printLocationO；/自定功能函数库MyFuncHon.h#dcfinc Pl 3.1415926#indudc using namespace std;doub
6、le aiiglcT)r)cgrcc(ini du,int fen,float miao);double angkToDugrsO；/将度分秒换算为度double dcgrccToRadiail(double degree);double dugreuToRadianO；将角度将算为弧度MyFuncHon.cpp#indudu NlyFunction.h Jdouble uiglcT)r)cgrcc(ini du,int fen,float miao)double rcsult=O;rcsult=miao/3600.0+fcii/60.0+du;return result;double iil
7、lglcTdufciimiao;rcsultanglcToDcgrccfdu,fcn.miao);return result;)double dcgrccT)Radiail(doublc degree) double rcsult=();rcsult=dcgrcc/57.295779513082321; return result;double dcgrccToRadiallQ !double result,degree；dcgrcc=ailglcT)r)cgrcc();rcsult=dcgrccT 1000#pragma once#undif/_MSC_VER 1000#dcfinc rou
8、Second 206264.806247096355#includc HMvFuiicti()n.hH#includc #includc using namespace std;class Zhungsuanpublic:ZhcngsuaiiQ;Zhcngsuaii(doublc fB,doublc fL);double grtXQ;double grtYO；double grtrBO；double grtrLO；void printLocationQ;virtual ZhengsuanQ;private：double x；double y;大地坐标double X;double Y;/国家统
9、一坐标double B；double rB；int Bsceond;double L;double rL;输入的大地纬度B,大地经度L,B,rL为对应弧度表示值，Bsucend为换算成秒 数值int n；带号 ndouble L0;中央经线纬度L0double LDot;纬度差 L-L0im LDotSucond;换算成秒的纬度差double 1;double N;double aO;double a3;double a4;double a5;double a6;七个计算参数;#clldif/ !dcfincd(AFX_ZHENGSUAN_H_2655EA28_E810_44A3_8F14_5
10、6421A7B4466_lNCLUDED_)ZhengsuaiLcpp/ Zhengsuan.cpp: implementation of the Zhciigsuaii class./#includc Zhungsuan上”/ Construction/Dustruction/Zhciigsu;ui:Zhciigsuan()Zhcngsum：:Zhcngsu;in(doublc fB,double fL)B=fB;rB=degrccToRadian(fB);L=fL;rL=dcgrccT()Radiail(fL);Bscc(md=B*360();初始化大地经度L,大地纬度B, Bsucond
11、,按弧度的大地纬度rBn=(int)(L/6+1);初始化带号 nL0=6*n-3;中央经线经度，角度单位LDt=L-LO;/经度差LD()iScc()nd=LDot3600;l=(LDot)*3600/imiSucond; 计算参数 1N=6399698.902-(21562.267-(108.973-0.612*cos(rB)*cos(rB)*cos(rB)*8s(rB)*cos(rB)*cos (田);计算参数Na0=32140.404(135.3302-(0.7092-0.004*cos(rB)*cos(rB)*cos(rB)*cos(rB)*cos(rB)*cos(rB); 计算参数
12、aO遍=(0.25+0.00252*(:(屈1记)*(:(彩(4)*口彩(市)*：05(田)-0.04166;计算参数 a446=(0.166*(:(夙电3以由)-0.084)*(:(与(电*。双市);计算参数26a3=(O.3333333+O.OO1123*cos(rB)*cos(rB)*8s(rB)*cos(rB)-0.l 666667;/计算参数 a3a5=0.0083-(0.l 667-(0.1968+0.004*cos(rB)*cos(rB)*cos(rB)*cos(rB)*cos(rB)*8s(rB); 计 算参数a5x=6367558.4969*Bsecond/rouSecond-(a0-(0.5+(a4+a6*l*l)*l*l)*l*l*N)*sin(rB)*cos(rB);/ 正算Xy=(l +(a3+a5*l*l)*l*l)*l*N*cos(rB);/IE# yX=x;Y=n*1000000+y+500000;国家统一坐标
