经纬度BL换算到高斯平面直角坐标XY（高斯投影正算）的源码及算法.zip

经纬度 BL 换算到高斯平面直角坐标 XY（高斯投影正算）的源码及算法 收藏

新一篇: C#的 6 种常用集合类大比拼 | 旧一篇: Silverlight 中使用图片及常见问题

一、

经纬度 BL 换算到高斯平面直角坐标 XY（高斯投影正算）的源码及算法

// GaussBL2xy.cpp : Defines the entry point for the console application.

//

#include "stdafx.h "

#include "math.h "

#include "CoorTrans.h "

double MyL = 109.22334444; //3 度带，109 d 22 m 33.4444 s

PrjPoint_Krasovsky MyPrj;

MyPrj.SetL0(MyL0);

MyPrj.SetBL(MyB, MyL);

double OutMyX; //结果应该大致是:3736714.783, 627497.303

double OutMyY;

OutMyX = MyPrj.x; //正算结果: 北坐标 x

OutMyY = MyPrj.y; //结果：东坐标 y

double OutputMyL;

OutputMyB = MyPrj.B; //反算结果：B

OutputMyL = MyPrj.L; //反算结果：L

//分析表明，此程序的结果和 Coord4.2 的转换结果是一样的，只差到毫米级

//原程序有几个问题，1.Pi 的值不对。2.SetBL 中多了两行错误代码

}

double Dms2Rad(double Dms)

{

double Degree, Miniute;

double Second;

Dms = fabs(Dms);

Degree = floor(Dms);

Miniute = floor(fmod(Dms * 100.0, 100.0));

Second = fmod(Dms * 10000.0, 100.0);

Rad = Sign * (Degree + Miniute / 60.0 + Second / 3600.0) * PI / 1

80.0;

return Rad;

{

double Second;

int Sign;

double Dms;

Sign = -1;

Miniute = floor(fmod(Rad * 60.0, 60.0));

Second = fmod(Rad * 3600.0, 60.0);

Dms = Sign * (Degree + Miniute / 100.0 + Second / 10000.0);

return Dms;

double X, N, t, t2, m, m2, ng2;

double sinB, cosB;

X = A1 * B * 180.0 / PI + A2 * sin(2 * B) + A3 * sin(4 * B)

+ A4 * sin(6 * B);

sinB = sin(B);