大地坐標(biāo)系和空間直角坐標(biāo)系轉(zhuǎn)換

1、#include#include#include#include#define MAX 100#define PI 3.1415926/結(jié)構(gòu)體存放 B, L, Hstruct BLH(char namel20;char name220;double B;char name320;double L;char name420;double H;lllllllllllllllllllllllllllllll結(jié)構(gòu)體存放 X,Y,Zstruct XYZ (char namel20; char name220; double X;char name320;double Y;char name420;dou

2、ble Z;lllllllllllllllllllllllllllllllll文件輸入 B,L,Hdouble FILEIN(BLH &b,char *fname) (ifstream in(fname,ios:nocreate); /建立文件流，并與輸入文件名建立關(guān)聯(lián) if(!in) coutfname error: file does not exist! endl;cout請手動(dòng)輸入數(shù)據(jù) nMendl;return 0;l;2;inb.B;3;inb.L;4;inb.H;in.close();return 1;lllllll

3、lllllllllllllllllllllllllllllllllllllllll文件輸入 X,YZdouble FILEIN1(XYZ &b,char *fname) (ifstream in(fname,ios:nocreate); /建立文件流，并與輸入文件名建立關(guān)聯(lián) if(!in) coutfname error: file does not exist! endl;cout請手動(dòng)輸入數(shù)據(jù)nHendl;return 0;l;2;inb.X;3;inb.Y;4;inb.Z;in.close();return 1;llllll

4、lllllllllllllllllllllllllllllllllllllllllllllllllldouble d_h(double angle)角度化弧度(double a,b;angle=modf(angle,&a);angle=modf(angle*100.0,&b);return (a+b/60.0+angle/36.0)*(PI+3.0E-16)/180.0;llllllllllllllllllllllllllllllllllllllllllllllllllllllllllldouble h_d(double angle)弧度化角度(double a,b,c;angle=modf(

5、angle*180.0/(PI-3.0E-16),&a);angle=modf(angle*60.0,&b);angle=modf(angle*60.0,&c);return a+b*0.01+c*0.0001+angle*0.0001;lllllllllllllllllllllllllllllllllllllllllll鍵盤輸入 B,L,Hjianpan(BLH &a)(cout請輸入大地坐標(biāo)B endl;cina.B;coutiW輸入大地坐標(biāo)L Hendl;cina.L;coutiW輸入大地坐標(biāo)H endl;cina.H;jianpanl(XYZ &a)(coutn請輸入空間直角坐標(biāo)X e

6、ndl;cina.X;coutn請輸入空間直角坐標(biāo)Y endl;cina.Y;coutn請輸入空間直角坐標(biāo)Z endl;cina.Z;llllllllllllllllllllllllllllllllllllllllllllllllllllll 標(biāo)系double B_X(BLH &a)(double aaMAXMAX;double al,bl,N,El,X,YZ;int kznzm;char outfilename20/fname20;cout請選擇：endl;coutHl：文件輸入endl;cout,2：鍵盤輸入endl;cinm;if(m=l)(coutniW 輸入文件名 Hendl;cin

7、fname;if(FILEIN(a,fname);else(if(jianpan(a);else鍵盤輸入xxz大地坐標(biāo)系換算到空間直角坐選擇數(shù)據(jù)輸入方式cout”輸入有誤，程序中止，請重新運(yùn)行Hendl;return false;else if(m=2)jianpan(a);elsecout “輸入代碼有誤，程序中止，請重新運(yùn)行” endl;return false;if(a.B90)cout“輸入B有誤,程序中止，請重新運(yùn)行Hendl;克拉索夫斯基橢球體Hendl;1975年國際橢球體endl;WGS-84 橢球體endl;2000中國大地坐標(biāo)系Mendl;選擇橢球體return false

8、;if(a.L360)cout“輸入L有誤，程序中止, return false;if(a.H0)cout“輸入H有誤，程序中止, return false;aa 0 0=6378245.0000000000;aal0=6378140.0000000000;aa20=6378137.0000000000;aa30=6378137.000;cout請選擇橢球體endl;coutH0:coutl:coutn2:coutH3:cink;請重新運(yùn)行” endl;請重新運(yùn)行Hendl;aa0l=6356863.0187730473;aall=6356755.2881575287;aa2l=6356752

9、.3142;aa3l=6356752.314;屏幕輸出結(jié)文件輸出結(jié)if(k!=0&k!=l&k!=2&k!=3)cout“輸入代碼有誤，程序中止，請重新運(yùn)行”endl;return false;couta.Bendl;al=aak0;bl=aakl;El=(al*al-bl*bl)/(al*al);a.B=d_h(a.B);a.L=d_h(a.L);N=al/sqrt(l-El*sin(a.B)*sin(a.B);X=(N+a.H)*cos(a.B)*cos(a.L);Y=(N+a.H)*cos(a.B)*sin(a.L);Z=(N*(l-El)+a.H)*sin(a.B);coutH轉(zhuǎn)換后的

10、空間直角坐標(biāo)是：” endl;果coutHX nXendl;coutY Yendl;coutHZ nZendl;cout要保存數(shù)據(jù)嗎？ endl;coutl:保存endl;coutH2 不保存endl;cinn;if(n=l)(coutn請輸入保存的文件名endl;果cinoutfilename;ofstream out(outfilename); out.precision(lO);if(!out) coutHcan not open save file!Nendl; outH空間直角坐標(biāo)endl; outX,endl; outXendl; out,Y,endl; outYendl; out

11、Z,endl; outZendl; out.close(); else if(n=2); else(cout輸入代碼有誤，程序中止，請重新運(yùn)行endl; return false; cout 謝謝您使用龔曉鵬編的程序，所有任務(wù)均巳完成，歡迎下次使用，祝生活愉快 endl;IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIII空間直角坐標(biāo)系換算到大地坐標(biāo)系 double X_B(XYZ &a) double al,bl,N,E1,B,L,H,BO,aaMAXMAX; double tanBMAXl; int i=0,n,k,m;cha

12、r outfilename20zfname20;選擇數(shù)據(jù)輸入方式 coutiW選擇:endl; coutHl：文件輸入endl; cout2：鍵盤輸入endl; cinm; if(m=l) (coutniW 輸入文件名 Hendl; cinfname; if(FILEINl(a,fname); else (if(jianpanl(a); else(cout輸入有誤，程序中止，請重新運(yùn)行endl; return false; else if(m=2)jianpanl(a);elsecout“輸入代碼有誤，程序中止，請重新運(yùn)行Hendl; return false;aa 0 0=6378245.0

13、000000000;aal0=6378140.0000000000;aa20=6378137.0000000000;aa30=6378137.000;aa0l=6356863.0187730473;aall=6356755.2881575287;aa2l=6356752.3142;aa3l=6356752.314;cout0.000000000000001);coutiendl;N=al/sqrt(l-El*sin(B)*sin(B);H=sqrt(a.X*a.X+a.Y*a.Y)/cos(B)-N;B=h_d(B);L=h_d(L);if(a.X0&a.Y0)L=180+L;else if(

14、a.X0)L=180+L;else if(a.X0&a.Y0)L=360+L;cout“轉(zhuǎn)換后的大地坐標(biāo)是：Nendl;屏幕輸 出結(jié)果coutHB Bendl;coutHL Lendl; coutH NHendl;cout要保存數(shù)據(jù)嗎？ endl;coutl:保存endl;cout2 不保存Mendl;cinn;if(n=l)cout“請輸入保存的文件名”endl;將結(jié)果保存到文件cinoutfilename;ofstream out(outfilename);out.precision(15);if(!out) coutMcan not open save file!Nendl;outH空間直角坐標(biāo)endl;out,B,endl;outBendl;outL,endl;outLendl;outH,endl;outHendl;out.close();else if(n=2);else(cout輸入代碼有誤，程序中止，請重新運(yùn)行,endl;return fa