经纬度与网络墨卡托(Web Mercator)投影坐标的转换

在GIS中，常会遇到不同投影坐标系之间的转换，下面以JavaScript对转换的方法进行实现：

一、 经纬度转网络墨卡托
方法一和方法二转换的结果有极小的差别，但都是正确的。

 //方法一  
latLng2WebMercator: function(lng, lat){
  //[114.32894, 30.585748]  
  var earthRad = 6378137.0;  
  var x = lng * Math.PI / 180 * earthRad;  
  var a = lat * Math.PI / 180;  
  var y = earthRad / 2 * Math.log((1.0 + Math.sin(a)) / (1.0 - Math.sin(a)));  
  return [x, y]; //[12727039.383734727, 3579066.6894065146]  
}  

    //方法二  
latLng2WebMercator2 : function(lng, lat){  
    var x = lng *20037508.34/180;  
    var y = Math.log(Math.tan((90+lat)*Math.PI/360))/(Math.PI/180);  
    y = y *20037508.34/180;  
    return [x, y];  
}  

二、网络墨卡托转经纬度

webMercator2LngLat: function(x, y) {//[12727039.383734727, 3579066.6894065146]  
    var lng = x / 20037508.34 * 180;  
    var lat = y / 20037508.34 * 180;      
    lat = 180 / Math.PI * (2 * Math.atan(Math.exp(lat * Math.PI / 180)) - Math.PI / 2);  
    return [lng, lat]; //[114.32894001591471, 30.58574800385281]  
}  

以上转换方法的结果经实测都是正确的，只需根据需要用其他语言实现即可。

参考文档：

http://blog.csdn.net/shaxiaozilove/article/details/54908569

http://malagis.com/mercator-longitude-latitude.html

https://www.cnblogs.com/yhlx125/p/3863748.html

https://www.cnblogs.com/DHUtoBUAA/p/6706642.html

http://blog.csdn.net/ryfdizuo/article/details/20793669(重要)

https://en.wikipedia.org/wiki/Scale_(map)

http://blog.csdn.net/williham/article/details/8496153

http://blog.csdn.net/witch_soya/article/details/7597598

http://blog.csdn.net/cleverysm/article/details/2121633（墨卡托投影实现）