/高斯投影正、反算
//6度带宽 54年北京坐标系
//高斯投影由经纬度(UnitD)反算大地坐标(含带号,Unit:Metres)
void GaussProjCal(double longitude, double latitude, double *X, double *Y)
{
int ProjNo=0; int ZoneWide; 带宽
double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval;
double a,f, e2,ee, NN, T,C,A, M, iPI;
iPI = 0.0174532925199433; 3.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 * 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(lat
高斯投影正反算的代码
最新推荐文章于 2024-01-25 14:16:50 发布
这是一个关于高斯投影计算的代码实现,包括从经纬度到大地坐标的正算和从大地坐标到经纬度的反算过程。代码支持6度带宽的54年北京坐标系,涉及的关键参数有NN(卯酉圈曲率半径)、M(子午线弧长)、fai(底点纬度)和R(底点曲率半径)。
摘要由CSDN通过智能技术生成