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

HYPERLINK /cuijie/archive/2008/01/06/2027845.aspx 经纬度BL换算到高斯平面直角坐标XY（高斯投影正算）的源码及算法 HYPERLINK JavaScript:d=document;t=d.selection?(d.selection.type!=None?d.selection.createRange().text:):(d.getSelection?d.getSelection():);void(saveit=window.open(/storeit.aspx?t=+escape(d.title)+u=+escape(d.location.href)+c=+escape(t),saveit,scrollbars=no,width=590,height=300,left=75,top=20,status=no,resizable=yes));saveit.focus(); \o 收藏到我的网摘中，并分享给我的朋友 收藏 HYPERLINK /cuijie/archive/2008/01/06/2027850.aspx 新一篇:?C#的6种常用集合类大比拼?|? HYPERLINK /cuijie/archive/2008/01/03/2013838.aspx 旧一篇:?Silverlight中使用图片及常见问题 一、经纬度BL换算到高斯平面直角坐标XY（高斯投影正算）的源码及算法 //???GaussBL2xy.cpp???:???Defines???the???entry???point???for???the???console???application.?//?#include???stdafx.h??#include???math.h??#include???CoorTrans.h??#include???iostream?using???namespace???std;?void???main(int???argc,???char*???argv[])?{???double???MyL0???=???108;???//中央子午线???double???MyB???=???33???//33???du???44???fen???55.6666???miao???double???MyL???=???109???//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???=??