文档介绍：高斯投影正算与反算的理论方法与实现代码
高斯投影是正形投影的一种，同一坐标系中的高斯投影换带计算公式是根据正形投影原理推导出的两个高斯坐标系间的显函数式。在同一大地坐标系中（例如19５４北京坐标系或1980西安坐标系)，如果两个高斯坐标系只是主子午线的经度不同，那么显函数式前的系数可以根据坐标系使用的椭球元素和主子午线经度唯一确定。但如果两个高斯坐标系除了主子午线的经度不同以外，还存在其他线性系，则将线性变换公式代入换带计算的显函数式中，仍然可以得到严密的坐标变换公式。此时显函数式前的系数等价于使用两个坐标系主子午线的经度和线性变换参数联合求解得到的,可以唯一确定。
//6度带宽 54北京坐标系
//高斯投影由大地坐标(Uｎit:Ｍｅｔres）反算经纬度（Uｎit：DD)
ｖoid GausｓＰrojＩnvCal（ｄｏuble X， double Y， ｄｏuble ＊lonｇｉtｕdｅ， douｂle ＊ｌａtituｄe)
｛
int ＰrojNo; int ZoneWide； ////带宽
doｕｂle lｏｎgitｕde１，lａｔitude1， loｎｇiｔｕｄe0,latiｔude0， Ｘ0,Y0， xｖal,ｙvaｌ；
dｏuｂle e1，ｅ２,f，a, eｅ, ＮＮ, T,C, M， D，R,u,fai, iPＩ；ﻫiPI = ； //／/3。14/; ﻫa ＝ 637824５．０； ｆ = 1．0/２98。3； ／/5４年北京坐标系参数
////a=； f＝１/2９8。257； //80年西安坐标系参数 ﻫＺonｅWiｄe ＝ 6； ／／//6度带宽
PrｏjNo = (iｎｔ）（X／
100００00L) ; //查找带号
ｌonｇiｔude0 = （ＰｒojＮo-1) * ZoneWide ＋ ZoｎeWidｅ / 2； ﻫｌoｎｇｉtude0 ＝ lｏnｇｉtudｅ0 ＊ iPI ; //中央经线ﻫX0 = ＰｒｏjＮo＊１000００0L+５０0000Ｌ; ﻫY０ = 0; ﻫxval = Ｘ—X0； yval ＝ Y—Y0； //带内大地坐标ﻫｅ２ ＝ 2＊f－f*f；
e1 ＝ (—sqrt(1－e2)）/（+sｑｒt(1—e２))；ﻫeｅ = e2／（1-e２);
M = ｙval；
ｕ = Ｍ/（ａ*（1—e2/4－3*e2＊e2/64－5＊e2*e2*e2/2５6));
ｆaｉ ＝ u+（3*e1/2—2７＊e１＊ｅ1＊e1／３２)＊ｓin（2＊u)+(21*e1＊ｅ１/1６—55*e1＊ｅ1*e1*e１／3２）*sｉｎ（ﻫ４＊u)＋（１５１*e1＊ｅ１*e1/96)＊ｓin（６＊u)+（１097*e1*e1*e1＊e1／512)＊ｓin(8＊u）；
Ｃ = eｅ＊cos（ｆai）*coｓ(fai）；
Ｔ ＝ tan(fai）*tａｎ(ｆai);
ＮN = a／sqrt(—e２*sｉｎ(fａi)＊sin（fａi))；
R = a＊（1—e2）/sｑrt(（１-e2＊sin(faｉ）＊ｓin(fai）)＊(1-e2*sin(fａi）*sin(fai)）＊（１—ｅ2＊sｉｎ
(ｆai）＊sｉｎ（fai）