高斯投影正反算的代码

最新推荐文章于 2025-06-02 13:03:38 发布

weixin_34362875

最新推荐文章于 2025-06-02 13:03:38 发布

阅读量3.7k

点赞数 2

CC 4.0 BY-SA版权

文章标签： matlab c/c++ c#

原文链接：http://www.cnblogs.com/imeiba/p/5696967.html

这是一个关于高斯投影计算的代码实现，包括从经纬度到大地坐标的正算和从大地坐标到经纬度的反算过程。代码支持6度带宽的54年北京坐标系，涉及的关键参数有NN（卯酉圈曲率半径）、M（子午线弧长）、fai（底点纬度）和R（底点曲率半径）。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

/高斯投影正、反算
//////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(latitude1);
C=ee*cos(latitude1)*cos(lat