高斯克吕格与经纬度坐标值转换代码.doc

《高斯克吕格与经纬度坐标值转换代码.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、高斯克吕格与经纬度坐标值转换代码[转帖]'WritenbyRodgerYuan952006'参考文献'v00.1'用于在经纬度坐标和高斯克吕格坐标之间的转换。'高斯克吕格为一种投影，根据椭球体和基准面不同又有所区分，常用的北京54和西安80即'采用这种投影方式，投影后的坐标为平面坐标系，单位为米'现在参数的坐标系采用测绘坐标系，x为纵坐标，y为横坐标'返回参数为自定义类型，双精度点'调用转换函数前需要调用初始化过程进行初始化'---------------------------------------------------------------

2、----------------'PublicSubinit(ByValTuoqiuCanshuAsCanshu,ByValDaihaoAsInteger)'说明:用于初始化转换参数'TuoqiuCanshu  枚举类型，提供北京54、西安80和WGS84三个椭球参数'Daihao  整型  为高斯克吕格投影六度分带带号，取值为1～60'-------------------------------------------------------------------------------'PublicSubinitEx(ByValdEAsDou

3、ble,ByValdNAsDouble,ByValk0AsDouble)'说明:用于进一步初始化转换参数（暂不提供）'dE  东偏移'dE  北偏移'k0  比例因子'-------------------------------------------------------------------------------'PublicFunctionJWgetGK(ByValWAsDouble,ByValJAsDouble)AsPointD'*****************************************************

4、*******************'基本变量定义DimaAsDouble'椭球体长半轴DimbAsDouble'椭球体短半周DimfAsDouble'扁率DimeAsDouble'第一偏心率DimeqAsDouble'第二偏心率DimdhAsInteger'带号DimFEAsDouble'东偏移DimFNAsDouble'北偏移DimL0AsDouble'中央经度Dimk0AsDouble'比例因子ConstPIAsDouble=3.14159265358979PublicEnumCanshu  Beijing54=0  Xian80=1  WG

5、S84=2EndEnumPublicTypePointD  XAsDouble  YAsDoubleEndTypePublicSubinit(ByValTuoqiuCanshuAsCanshu,ByValDaihaoAsInteger)SelectCaseTuoqiuCanshu'Krassovsky（北京54采用）63782456356863.0188'IAG75（西安80采用）63781406356755.2882'WGS8463781376356752.3142Case0:'北京五四  a=6378245  b=6356863.0188Case

6、1:'西安八零  a=6378140  b=6356755.2882Case2:'WGS84  a=6378137  b=6356752.3142EndSelectf=(a-b)/ae=Sqr(1-(b/a)^2)eq=Sqr((a/b)^2-1)IfDaihao<1OrDaihao>60ThenExitSubdh=DaihaoL0=(6*dh-3)*PI/180k0=1FE=500000+dh*1000000FN=0EndSubPublicSubinitEx(ByValdEAsDouble,ByValdNAsDouble,ByValdk0AsDou

7、ble)EndSubPublicFunctionJWgetGK(ByValWAsDouble,ByValJAsDouble)AsPointD'给出经纬度坐标，转换为高克投影坐标DimBYAsDoubleDimLXAsDoubleDimTCAsDoubleDimCCAsDoubleDimACAsDoubleDimMCAsDoubleDimNCAsDoubleDimrxAsDoubleDimryAsDoubleDimresultPAsPointDBY=W*PI/180LX=J*PI/180TC=Math.Tan(BY)^2CC=eq^2*Cos(BY)^

8、2AC=(LX-L0)*Cos(BY)MC=a*((1-e^2/4-3*e^4/64-5*e^6/256)*