经纬度墨卡托互相转换-WGS84 and Mercator transfer each other

转自barry.z-
http://bbs.esrichina-bj.cn/ESRI/thread-78245-1-1.html


//经纬度转墨卡托
-(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 ;
}
//墨卡托转经纬度
-(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;
}

[本日志由 tangc 于 2011-01-26 03:37 PM 编辑]
文章来自: 转载
引用通告: 查看所有引用 | 我要引用此文章
Tags: 经纬度 墨卡托 转换 WGS84 Mercator
相关日志:
评论: 0 | 引用: 0 | 查看次数: 368
发表评论
昵 称:
密 码: 游客发言不需要密码.
内 容:
验证码: 验证码
选 项:
虽然发表评论不用注册,但是为了保护您的发言权,建议您注册帐号.
字数限制 1000 字 | UBB代码 开启 | [img]标签 关闭