经纬度坐标与高斯坐标的转换代码精编版.docVIP

……………………………………………………………最新资料推荐………………………………………………… PAGE PAGE 1 经纬度坐标与高斯坐标的转换代码 /功能说明： 将绝对高斯坐标(y,x)转换成绝对的地理坐标(wd,jd)。 / // double y; 输入参数: 高斯坐标的横坐标，以米为单位 // double x; 输入参数: 高斯坐标的纵坐标，以米为单位 // short DH; 输入参数: 带号，表示上述高斯坐标是哪个带的 // double *L; 输出参数: 指向经度坐标的指针，其中经度坐标以秒为单位 // double *B; 输出参数: 指向纬度坐标的指针，其中纬度坐标以秒为单位 void GaussToGeo(double y, double x, short DH, double *L, double *B, double LP) { double l0; // 经差 double tf; // tf = tg(Bf0),注意要将Bf转换成以弧度为单位 double nf ; // n = y * sqrt( 1 + etf ** 2) / c, 其中etf = e**2 * cos(Bf0) ** 2 double t_l0; // l0，经差，以度为单位 double t_B0; // B0，纬度，以度为单位 double Bf0; // Bf0 double etf; // etf,其中etf = e**2 * cos(Bf0) ** 2 double X_3 ; double PI=3.14159265358979; double b_e2=0.0067385254147; double b_c=6399698 X_3 = x / 1000000.00 - 3 ; // 以兆米（1000000）为单位 // 对于克拉索夫斯基椭球，计算Bf0 Bf0 = 27.11115372595 + 9.02468257083 * X_3 - 0.00579740442 * pow(X_3,2) - 0.00043532572 * pow(X_3,3) + 0.00004857285 * pow(X_3,4) + 0.00000215727 * pow(X_3,5) - 0.00000019399 * pow(X_3,6) ; tf = tan(Bf0*PI/180); // tf = tg(Bf),注意这里将Bf转换成以弧度为单位 etf = b_e2 * pow(cos(Bf0*PI/180),2); // etf = e**2 * cos(Bf) ** 2 nf = y * sqrt( 1 + etf ) / b_c; // n = y * sqrt( 1 + etf ** 2) / c // 计算纬度，注意这里计算出来的结果是以度为单位的 t_B0 = Bf0 - (1.0+etf) * tf / PI * (90.0 * pow(nf,2) - 7.5 * (5.0 + 3 * pow(tf,2) + etf - 9 * etf * pow(tf,2)) * pow(nf,4) + 0.25 * (61 + 90 * pow(tf,2) + 45 * pow(tf,4)) * pow(nf,6)) ; // 计算经差，注意这里计算出来的结果是以度为单位的 t_l0 = (180 * nf - 30 * ( 1 + 2 * pow(tf,2) + etf ) * pow(nf,3) + 1.5 * (5 + 28 * pow(tf,2) + 24 * pow(tf,4)) * pow(nf,5)) / ( PI * cos(Bf0*PI/180) ) ; l0 = (t_l0 * 3600.0); // 将经差转成秒 if (LP == -1000) { *L = (double)((DH * 6 - 3) * 3600.0 + l0); // 根据带号计算出以秒为单位的绝对经度，返回指针 } else { *L = LP * 3600.0 + l0; // 根据带号计算出以秒为单位的绝对经度，返回指针 }