【200905】WGS84墨卡托投影转换C#源码

 

代码：

WGS84 Mercator project转换的算法（C#）

    N0 = 6378137.0 / Math.Sqrt( 1-Math.Pow(0.081819190843,2)*Math.Pow(Math.Sin（RefLat*Math.PI/180),2) );
      q1 = Math.Log( Math.Tan( (180.0/4.0+dLat/2.0)*Math.PI/180.0 ) );
      q2 = 0.081819190843/2 * Math.Log( (1+0.081819190843*Math.Sin(dLat*Math.PI/180.0) ) / 
       (1-0.081819190843*Math.Sin(dLat*Math.PI/180.0) )  );
      q = q1 - q2 ;
      x = N0 * Math.Cos(RefLat*Math.PI/180.0) * ((dLong-CentralMeridian)/57.29577951) ;
      y = N0 * Math.Cos(RefLat*Math.PI/180.0) * q ;


式中
输入参数：dLat----纬度

               dLong----经度

               RefLat----基准纬度

              CentralMeridian----中央子午线

输出：x----横坐标,y----纵坐标

其余参数为中间过程参数