经纬度BL换算到高斯平面直角坐标XY（高斯投影正算）的源码及算法 收藏
新一篇: C#的6种常用集合类大比拼 | 旧一篇: Silverlight中使用图片及常见问题
一、
经纬度BL换算到高斯平面直角坐标XY（高斯投影正算）的源码及算法
// GaussBL2xy.cpp : Defines the entry point for the console application. //
#include \#include \#include \#include
using namespace std;
void main(int argc, char* argv[]) {
double MyL0 = 108; //中央子午线
double MyB = 33.44556666; //33 du 44 fen 55.6666 miao 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 InputMyX = 3736714.783; //如果是独立计算，应该给出中央子午线L0 double InputMyY = 627497.303; MyPrj.Setxy(InputMyX, InputMyY);
MyPrj.GetBL(&MyPrj.B, &MyPrj.L); //把计算出的BL的弧度值换算为dms形式 double OutputMyB; 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; int Sign; double Rad; if(Dms > = 0) Sign = 1; else
Sign = -1; 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 Rad2Dms(double Rad) {
double Degree, Miniute; double Second; int Sign; double Dms; if(Rad > = 0) Sign = 1; else
Sign = -1;
Rad = fabs(Rad * 180.0 / PI); Degree = floor(Rad);
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; }
/////////////////////////////////////////////////// // Definition of PrjPoint
/////////////////////////////////////////////////// bool PrjPoint::BL2xy() {
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);