《一个老师给的高斯投影正、反算c++源码(最新整理)》由会员分享，可在线阅读，更多相关《一个老师给的高斯投影正、反算c++源码(最新整理)(4页珍藏版)》请在人人文库网上搜索。

1、高斯投影正、反算/6 度带宽54 年北京坐标系/高斯投影由经纬度(Unit:DD)反算大地坐标(含带号，Unit:Metres)void GaussProjCal(double longitude, double latitude, double *X, doubl e *Y)int ProjNo=0;int ZoneWide;/带宽doublelongitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval; doublea,f, e2,ee, NN, T,C,A, M, iPI;iPI = 0.0174532925199433;/3.。

2、1415926535898/180.0;ZoneWide = 6;/6 度带宽a=6378245.0;f=1.0/298.3;/54 年北京坐标系参数/a=6378140.0;f=1/298.257;/80 年西安坐标系参数ProjNo = (int)(longitude / ZoneWide) ; longitude0 = ProjNo * ZoneWide + ZoneWide / 2; longitude0 = longitude0 * iPI ;latitude0=0;longitude1 = longitude * iPI ; /经度转换为弧度latitude1 = latitude。

3、 * iPI ; /纬度转换为弧度e2=2*f-f*f;ee=e2*(1.0-e2);NN=a/sqrt(1.0-e2*sin(latitude1)*sin(latitude1); T=tan(latitude1)*tan(latitude1); C=ee*cos(latitude1)*cos(latitude1); A=(longitude1-longitude0)*cos(latitude1);M=a*(1-e2/4-3*e2*e2/64-5*e2*e2*e2/256)*latitude1-(3*e2/8+3*e2*e2/32+45*e2*e2*e2/1024)*sin(2*latitud。

4、e1)+(15*e2*e2/256+45*e2*e2*e2/1024)*sin(4*latitude1)-(35*e2*e2*e2/3072)* sin(6*latitude1);xval = NN*(A+(1-T+C)*A*A*A/6+(5-18*T+T*T+72*C-58*ee)*A*A*A*A*A/120); yval = M+NN*tan(latitude1)*(A*A/2+(5-T+9*C+4*C*C)*A*A*A*A/24+(61-58*T+T*T+600*C-330*ee)*A*A*A*A*A*A/720); X0 = 1000000L*(ProjNo+1)+500000L;Y0。

5、 = 0;xval = xval+X0;yval = yval+Y0;*X = xval;*Y = yval;/高斯投影由大地坐标(Unit:Metres)反算经纬度(Unit:DD)void GaussProjInvCal(double X, double Y, double *longitude, double *l atitude)int ProjNo;int ZoneWide;/带宽double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval; double e1,e2,f,a, ee, NN, T,C, M, 。

6、D,R,u,fai, iPI;iPI = 0.0174532925199433;/3.1415926535898/180.0;a = 6378245.0;f = 1.0/298.3;/54 年北京坐标系参数/a=6378140.0;f=1/298.257;/80 年西安坐标系参数ZoneWide = 6;/6 度带宽ProjNo = (int)(X/1000000L) ;/查找带号longitude0 = (ProjNo-1) * ZoneWide + ZoneWide / 2; longitude0 = longitude0 * iPI ;/中央经线X0 = ProjNo*1000000L+。

7、500000L; Y0 = 0;xval = X-X0;yval = Y-Y0;/带内大地坐标e2 = 2*f-f*f;e1 = (1.0-sqrt(1-e2)/(1.0+sqrt(1-e2); ee = e2/(1-e2);M = yval;u = M/(a*(1-e2/4-3*e2*e2/64-5*e2*e2*e2/256);fai = u+(3*e1/2-27*e1*e1*e1/32)*sin(2*u)+(21*e1*e1/16-55*e1*e1*e1*e1/32)*sin(4*u)+(151*e1*e1*e1/96)*sin(6*u)+(1097*e1*e1*e1*e1/512)*si。

8、n(8*u); C = ee*cos(fai)*cos(fai);T = tan(fai)*tan(fai);NN = a/sqrt(1.0-e2*sin(fai)*sin(fai);R = a*(1-e2)/sqrt(1-e2*sin(fai)*sin(fai)*(1-e2*sin(fai)*sin(fai)* (1-e2*sin(fai)*sin(fai);D = xval/NN;/计算经度(Longitude) 纬度(Latitude)longitude1 = longitude0+(D-(1+2*T+C)*D*D*D/6+(5-2*C+28*T-3*C*C+8*ee+ 24*T*T)*。

9、D*D*D*D*D/120)/cos(fai);latitude1 = fai -(NN*tan(fai)/R)*(D*D/2-(5+3*T+10*C-4*C*C-9*ee)*D*D* D*D/24+(61+90*T+298*C+45*T*T-256*ee-3*C*C)*D*D*D*D*D*D/720);/转换为度 DD*longitude = longitude1 / iPI;*latitude = latitude1 / iPI;如果有需要程序的，可以直接跟我联系，呵呵附：高斯正反算参数pi=0.0174532925 0.0174532925199433 /长半轴 a=6378245.0;。

10、扁率 f=1.0/298.3;/54 年北京坐标系参数长半轴 a=6378140.0;扁率 f=1/298.257;/80 年西安坐标系参数长半轴 a=6378137m；扁率 f=1:298.257223563。/WGS-84 坐标系“”“”At the end, Xiao Bian gives you a passage. Minand once said, people who learn to learn are very happy people. In every wonderful life, learning is an eternal theme. As a professio。

11、nal clerical and teaching position, I understand the importance of continuous learning, life is diligent, nothing can be gained, only continuous learning can achieve better self. Only by constantly learning and mastering the latest relevant knowledge, can employees from all walks of life keep up with the pace of enterprise development and innovate to meet the needs of the market. This document is also edited by my studio professionals, there may be errors in the document, if there are errors, please correct, thank you。