墨卡托坐标与经纬度相互转换的方法

地质BOSS

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

麦卡托投影法

Web Mercator

墨卡托坐标与经纬度相互转换的方法

经纬度转Wev墨卡托

版本1

1. dvec3 CMathEngine::lonLat2WebMercator(dvec3  lonLat)  
   
2. {  
   
3.     dvec3  mercator;  
   
4.     double x = lonLat.x *20037508.34/180;  
   
5.     double y = log(tan((90+lonLat.y)*PI/360))/(PI/180);  
   
6.     y = y *20037508.34/180;  
   
7.     mercator.x = x;  
   
8.     mercator.y = y;  
   
9.     return mercator ;  
   
10. }

复制代码

版本2

1. (CGPoint )lonLat2Mercator:(CGPoint ) lonLat{
   
2.     CGPoint mercator;
   
3.     double x = lonLat.x *20037508.34/180;
   
4.     double  y = log(tan((90+lonLat.y)*M_PI/360))/(M_PI/180);
   
5.     y = y *20037508.34/180; mercator.x = x; mercator.
   
6.     y = y; return mercator ;
   
7. }

复制代码

版本3

1. function lonlat2mercator(lonlat)
   
2. {   
   
3.     var mercator={x:0,y:0};
   
4.     var x = lonlat.x *20037508.34/180;  
   
5.     var y = Math.log(Math.tan((90+lonlat.y)*Math.PI/360))/(Math.PI/180);
   
6.     y = y *20037508.34/180;
   
7.     mercator.x = x;
   
8.     mercator.y = y;
   
9.     return mercator ;
   
10. }

复制代码

Web墨卡托转经纬度

版本1

1. dvec3 CMathEngine::WebMercator2lonLat( dvec3   mercator )  
   
2. {  
   
3.     dvec3 lonLat;  
   
4.     double x = mercator.x/20037508.34*180;  
   
5.     double y = mercator.y/20037508.34*180;  
   
6.     y= 180/PI*(2*atan(exp(y*PI/180))-PI/2);  
   
7.     lonLat.x = x;  
   
8.     lonLat.y = y;  
   
9.     return lonLat;  
   
10. }  

复制代码

版本2

1. (CGPoint )Mercator2lonLat:(CGPoint ) mercator
   
2. {
   
3.     CGPoint lonLat;
   
4.     double x = mercator.x/20037508.34*180;
   
5.     double y = mercator.y/20037508.34*180;
   
6.     y= 180/M_PI*(2*atan(exp(y*M_PI/180))-M_PI/2);
   
7.     lonLat.x = x; lonLat.y = y;
   
8.     return lonLat;
   
9. }

复制代码

版本3

1. function mercator2lonlat(mercator){
   
2.     var lonlat={x:0,y:0};   
   
3.     var x = mercator.x/20037508.34*180;
   
4.     var y = mercator.y/20037508.34*180;
   
5.     y= 180/Math.PI*(2*Math.atan(Math.exp(y*Math.PI/180))-Math.PI/2);
   
6.     lonlat.x = x;   lonlat.y = y;   return lonlat;
   
7. }

复制代码

参考：
http://www.thinkgis.cn/topic/54199646bcea5603e8016e1d

qinccckencn

謝謝分享好好學習一下

tsnc

感谢分享！

路鹏123

感谢分享:mg

wsbd2016

回复是对楼主最基本的尊重！！！

sskcdn

谢谢楼主分享