高斯投影坐標正反算編程報告參考模板

1、高斯投影坐標正反算編程報告1. 編程思想進行高斯投影坐標正反算的編程需要牽涉到大量的公式，為了使程序條理更清楚，各塊的數據復用性更強，這里采取了結構化的編程思想。程序由四大塊組成。GeodesyHomework.cpp文件用于存放main()函數，是整個程序的入口。通過結構化的編程盡力使main()函數變得簡單。MyFunction.h和MyFunction.cpp用于存放計算過程中進行角度弧度換算時所要用到的一些自定的轉換函數。Zhengsuan.h和Zhengsuan.cpp用于存放Zhengsuan類，在Zhengsuan類中聲明了高斯投影坐標正算所要用到的所有變量，在類的構造函數中進行

2、成員變量的初始化及正算計算。通過get函數獲得相應的正算結果。Fansuan.h和Fansuan.cpp用于存放Fansuan類，類似于Zhengsuan類，F(xiàn)ansuan類中聲明了高斯投影坐標反算所要用到的所有變量，在類的構造函數中進行成員變量的初始化及反算計算。通過get函數獲得相應的反算結果。2. 計算模型 高斯投影正算公式 高斯投影反算公式 3. 程序框圖1 / 12開始輸入B，L求定帶號N，中央緯度L0，緯度差l按照實用公式計算x,y換算為國家統(tǒng)一坐標X，Y輸出X,Y輸入國家統(tǒng)一坐標X,Y由Y取定帶號N，并換算出x，y求出中央經線L0按照實用公式計算B，lL=L0+l求出大地經度L輸

3、出B，L結束正算反算4. 計算結果5. 附錄：程序代碼/主函數入口GeodesyHomework.cpp#include "MyFunction.h"#include "Zhengsuan.h"#include "Fansuan.h"#include <iostream>using namespace std;void fansuan();void zhengsuan();void main()zhengsuan();fansuan();printf("/n over!");void zhengsuan

4、()double myB,myL;cout<<"【正算】"<<endl;cout<<"請輸入大地緯度B"<<endl;myB=angleToDegree();cout<<"請輸入大地經度L"<<endl;myL=angleToDegree();Zhengsuan myZhengsuan1(myB,myL);printf("Radian B=%f L=%f n",myZhengsuan1.getrB(),myZhengsuan1.getrL();

5、myZhengsuan1.printLocation();void fansuan()double myX,myY;cout<<"【反算】"<<endl;cout<<"請輸入國家統(tǒng)一坐標 X Y。例如 3378627.1819 20243953.4517"<<endl;cin>>myX>>myY;Fansuan myFansuan1(myX,myY);myFansuan1.printLocation();/自定功能函數庫MyFunction.h#define PI 3.1415926

6、#include <iostream>using namespace std;double angleToDegree(int du,int fen,float miao);double angleToDegree();/將度分秒換算為度double degreeToRadian(double degree);double degreeToRadian();/將角度換算為弧度MyFunction.cpp#include "MyFunction.h"double angleToDegree(int du,int fen,float miao)double resu

7、lt=0;result=miao/3600.0+fen/60.0+du;return result;double angleToDegree()int du,fen;float miao;double result;cout<<"請輸入度 分 秒。例如：30 20 00"<<endl;cin>>du>>fen>>miao;result=angleToDegree(du,fen,miao);return result;double degreeToRadian(double degree)double result=

8、0;result=degree/57.295779513082321;return result;double degreeToRadian()double result,degree;degree=angleToDegree();result=degreeToRadian(degree);return result;/正算類Zhengsuan.h/ Zhengsuan.h: interface for the Zhengsuan class./#if !defined(AFX_ZHENGSUAN_H_2655EA28_E810_44A3_8F14_56421A7B4466_INCLUDED_

9、)#define AFX_ZHENGSUAN_H_2655EA28_E810_44A3_8F14_56421A7B4466_INCLUDED_#if _MSC_VER > 1000#pragma once#endif / _MSC_VER > 1000#define rouSecond 206264.806247096355#include "MyFunction.h"#include <iostream>#include <math.h>using namespace std;class Zhengsuan public:Zhengsu

10、an();Zhengsuan(double fB,double fL);double getX();double getY();double getrB();double getrL();void printLocation();virtual Zhengsuan();private:double x;double y;/大地坐標double X;double Y;/國家統(tǒng)一坐標double B;double rB;int Bsecond;double L;double rL;/輸入的大地緯度B，大地經度L，rB,rL為對應弧度表示值，Bsecond為換算成秒數值int n;/帶號ndoubl

11、e L0;/中央經線緯度L0double LDot;/緯度差L-L0int LDotSecond;/換算成秒的緯度差double l;double N;double a0;double a3;double a4;double a5;double a6;/七個計算參數;#endif / !defined(AFX_ZHENGSUAN_H_2655EA28_E810_44A3_8F14_56421A7B4466_INCLUDED_)Zhengsuan.cpp/ Zhengsuan.cpp: implementation of the Zhengsuan class./#include "Z

12、hengsuan.h"/ Construction/Destruction/Zhengsuan:Zhengsuan()Zhengsuan:Zhengsuan(double fB,double fL)B=fB;rB=degreeToRadian(fB);L=fL;rL=degreeToRadian(fL);Bsecond=B*3600;/初始化大地經度L，大地緯度B，Bsecond,按弧度的大地緯度rBn=(int)(L/6+1);/初始化帶號nL0=6*n-3;/中央經線經度，角度單位LDot=L-L0;/經度差LDotSecond=LDot*3600;l=(LDot)*3600/r

13、ouSecond;/計算參數lN=6399698.902-(21562.267-(108.973-0.612*cos(rB)*cos(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)*co

14、s(rB)*cos(rB);/計算參數a6a3=(0.3333333+0.001123*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/rouSecond-(a0-(0.5+(a4+a6*l*l)*l*l)*l*l*N)*sin(rB)*cos(rB);/正算xy=(1+(a3+a5*l*l)*l*l)*l*N*cos(rB);

15、/正算yX=x;Y=n*1000000+y+500000;/國家統(tǒng)一坐標Zhengsuan:Zhengsuan()double Zhengsuan:getX()return X;double Zhengsuan:getY()return Y;void Zhengsuan:printLocation()printf("正算得國家統(tǒng)一坐標為： X= %8.8f Y=%8.8f n",X,Y);double Zhengsuan:getrB()return rB;double Zhengsuan:getrL()return rL;/反算類Fansuan.h/ Fansuan.h:

16、interface for the Fansuan class./#if !defined(AFX_FANSUAN_H_5B5E4668_EE81_463F_8D65_FFE2AEACE745_INCLUDED_)#define AFX_FANSUAN_H_5B5E4668_EE81_463F_8D65_FFE2AEACE745_INCLUDED_#if _MSC_VER > 1000#pragma once#endif / _MSC_VER > 1000#define rouSecond 206264.806247096355#include <math.h>#inc

17、lude "MyFunction.h"#include <iostream>using namespace std;class Fansuan public:Fansuan();Fansuan(double X,double Y);double getB();double getL();void printLocation();virtual Fansuan();private:double x;double y;/高斯投影坐標double X;double Y;int N;/國家統(tǒng)一坐標,N為帶號double B,Bsecond;double L;/最后反算得

18、到B、Ldouble L0;/中央經線經度double l,lsecond;/L=L0+l，L0=6*N-3double Bf,BfSecond,BfDegree;double beta,betaSecond,betaDegree;double Z;double Nf;double b2;double b3;double b4;double b5;/計算的8個參數;#endif / !defined(AFX_FANSUAN_H_5B5E4668_EE81_463F_8D65_FFE2AEACE745_INCLUDED_)Fansuan.cpp/ Fansuan.cpp: implementat

19、ion of the Fansuan class./#include "Fansuan.h"/ Construction/Destruction/Fansuan:Fansuan()Fansuan:Fansuan(double X,double Y)this->X=X;this->Y=Y;/初始化x,yN=(int)(Y/1000000);/取出帶號L0=6*N-3;/初始化該帶號的中央經線經度x=X;y=Y-1000000*N-500000;beta=x/6367558.4969;/初始化beta，弧度單位betaSecond=beta*rouSecond;/初始化beta，秒單位betaDegree=betaSecond/3600;/初始化beta,整度數單位Bf=beta+(50221746+(293622+(2350+22*cos(beta)*cos(beta)*cos(beta)*cos(beta)*cos(beta)*cos(beta)*(1e-10)*sin(beta)*cos(beta);/初始化Bf,弧度單位BfSecond=Bf*rouSecond;/初始化Bf，秒單位BfDegree=BfSecond/3600;/初始化Bf，整度數單位Nf=6399698.902-(215