WGS84經緯度座標與Web墨卡託之間的轉換

//經緯度轉Web墨卡託
dvec3 CMathEngine::lonLat2WebMercator(dvec3  lonLat)
{
    dvec3  mercator;
    double x = lonLat.x *20037508.34/180;
    double y = log(tan((90+lonLat.y)*PI/360))/(PI/180);
    y = y *20037508.34/180;
    mercator.x = x;
    mercator.y = y;
    return mercator;
}
//Web墨卡託轉經緯度
dvec3 CMathEngine::WebMercator2lonLat( dvec3   mercator )
{
    dvec3 lonLat;
    double x = mercator.x/20037508.34*180;
    double y = mercator.y/20037508.34*180;
    y= 180/PI*(2*atan(exp(y*PI/180))-PI/2);
    lonLat.x = x;
    lonLat.y = y;
    return lonLat;
}





第二種方法：


//經度轉墨卡託
function handle_x(x) {
 return (x / 180.0) * 20037508.34;
}

//緯度轉墨卡託
function handle_y(y) {
 if (y > 85.05112) {
  y = 85.05112;
 }

 if (y <</span> -85.05112) {
  y = -85.05112;
 }

 y = (Math.PI / 180.0) * y;
 var tmp = Math.PI / 4.0 + y / 2.0;
 return 20037508.34 * Math.log(Math.tan(tmp)) / Math.PI;

}

//墨卡託轉經度
function handle_me_x(x)
{
 return x/20037508.34*180;
}

//墨卡託轉緯度
function handle_me_y(my)
{
  var mmy = my/20037508.34*180;
  y= 180/Math.PI*(2*Math.atan(Math.exp(mmy*Math.PI/180))-Math.PI/2);
  return y;
}