高斯正算C语言程序,一个老师给的高斯投影正、反算c++源码(最新整理)-CSDN博客

本文提供了高斯投影的C/C++源码，包括从经纬度到大地坐标的正算和从大地坐标到经纬度的反算。代码适用于6度带宽的54年北京坐标系和80年西安坐标系，详细解释了关键参数和计算步骤。

摘要由CSDN通过智能技术生成

《一个老师给的高斯投影正、反算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=