最近用到了墨卡托坐标转换经纬度的需求，原理什么的不想多说，需要的可以查看这里：

麦卡托投影法

Web Mercator

经纬度转Wev墨卡托

版本1

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 ;

}

版本2

(CGPoint )lonLat2Mercator:(CGPoint )

lonLat{

CGPoint mercator;

double x = lonLat.x *20037508.34/180;

double  y =

log(tan((90+lonLat.y)*M_PI/360))/(M_PI/180);

y

= y *20037508.34/180; mercator.x = x; mercator.

y

= y; return mercator ;

}

版本3

function lonlat2mercator(lonlat)

{

var mercator={x:0,y:0};

var x = lonlat.x *20037508.34/180;

var y = Math.log(Math.tan((90+lonlat.y)*Math.PI/360))/(Math.PI/180);

y

= y *20037508.34/180;

mercator.x = x;

mercator.y = y;

return mercator ;

}

Web墨卡托转经纬度

版本1

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;

}

版本2

(CGPoint )Mercator2lonLat:(CGPoint )

mercator

{

CGPoint lonLat;

double x = mercator.x/20037508.34*180;

double y = mercator.y/20037508.34*180;

y= 180/M_PI*(2*atan(exp(y*M_PI/180))-M_PI/2);

lonLat.x = x; lonLat.y = y;

return lonLat;

}

版本3

function mercator2lonlat(mercator){

var lonlat={x:0,y:0};

var x = mercator.x/20037508.34*180;

var y = mercator.y/20037508.34*180;

y= 180/Math.PI*(2*Math.atan(Math.exp(y*Math.PI/180))-Math.PI/2);

lonlat.x = x;   lonlat.y = y;   return lonlat;

}