文档章节

关于百度坐标系 (BD-09)与火星坐标系 (GCJ-02)以及WGS84坐标之间的互相转换

冰泪_
 冰泪_
发布于 2017/04/11 11:11
字数 781
阅读 4.5K
收藏 0

「深度学习福利」大神带你进阶工程师,立即查看>>>

 

由于自己工作需要 坐标系的转换 所以查资料翻译了坐标系转换的 oc 版本  

//
//  coordinatesChangeHelper.h
//
//  Created by 冰泪 on 2017/4/10.
//  Copyright © 2017年 冰泪. All rights reserved.
//

#import <Foundation/Foundation.h>
#import <MapKit/MapKit.h>
//#import "coordinatesChange.h"

@interface BLCoordinatesChangeHelper : NSObject
@property (nonatomic,assign) float  x_PI;
@property (nonatomic,assign) float  PI;
@property (nonatomic,assign) float  a;
@property (nonatomic,assign) float  ee;
+ (instancetype)shared;
/**
 * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换
 * 即 百度 转 谷歌、高德
 */
- (CLLocationCoordinate2D)bl_bd09togcj02:(double) lng andLat:(double) lat ;
/**
 * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换
 * 即谷歌、高德 转 百度
 */
- (CLLocationCoordinate2D)bl_gcj02tobd09:(double)lng andLat:(double) lat ;

/**
 * WGS84转GCj02
 */
- (CLLocationCoordinate2D)bl_wgs84togcj02:(double)lng andLat:(double) lat;

/**
 * GCJ02 转换为 WGS84
 */
- (CLLocationCoordinate2D)bl_gcj02towgs84:(double)lng andLat:(double)lat ;

@end

 

 

//
//  coordinatesChangeHelper.m
//
//  Created by 冰泪 on 2017/4/10.
//  Copyright © 2017年 冰泪. All rights reserved.
//

#import "BLCoordinatesChangeHelper.h"

@implementation BLCoordinatesChangeHelper
/**
 * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换
 * 即 百度 转 谷歌、高德
 */
- (CLLocationCoordinate2D)bl_bd09togcj02:(double) lng andLat:(double) lat {
    double x_pi = 3.14159265358979324 * 3000.0 / 180.0;
    double x = lng - 0.0065;
    double y = lat - 0.006;
    double z = sqrt(x * x + y * y) - 0.00002 * sin(y * x_pi);
    double theta = atan2(y, x) - 0.000003 * cos(x * x_pi);
    double bl_lng = z * cos(theta);
    double bl_lat = z * sin(theta);
    CLLocationCoordinate2D coor;
    coor.longitude = bl_lng;
    coor.latitude = bl_lat;
    return coor;
}
/**
 * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换
 * 即谷歌、高德 转 百度
 */
- (CLLocationCoordinate2D)bl_gcj02tobd09:(double)lng andLat:(double) lat {
    double z = sqrt(lng * lng + lat * lat) + 0.00002 * sin(lat * _x_PI);
    double theta = atan2(lat, lng) + 0.000003 * cos(lng * _x_PI);
    double bl_lng = z * cos(theta) + 0.0065;
    double bl_lat = z * sin(theta) + 0.006;
    CLLocationCoordinate2D coor;
    coor.longitude = bl_lng;
    coor.latitude = bl_lat;
    return coor;
}

/**
 * WGS84转GCj02
 */
- (CLLocationCoordinate2D)bl_wgs84togcj02:(double)lng andLat:(double) lat {
    if ([self out_of_china:lng andLat:lat ]) {
        CLLocationCoordinate2D coor;
        coor.longitude = lng;
        coor.latitude = lat;
        return coor;
    }else {
        double dlat = [self transformlat:lng - 105.0 andLat:lat - 35.0];
        double dlng = [self transformlng:lng - 105.0 andLat:lat - 35.0];
        double radlat = lat / 180.0 * _PI;
        double magic = sin(radlat);
        magic = 1 - _ee * magic * magic;
        double sqrtmagic = sqrt(magic);
        dlat = (dlat * 180.0) / ((_a * (1 - _ee)) / (magic * sqrtmagic) * _PI);
        dlng = (dlng * 180.0) / (_a / sqrtmagic * cos(radlat) * _PI);
        double bl_lat = lat + dlat;
        double bl_lng = lng + dlng;
        
        CLLocationCoordinate2D coor;
        coor.longitude = bl_lng;
        coor.latitude = bl_lat;
        return coor;
    }
}

/**
 * GCJ02 转换为 WGS84
 */
- (CLLocationCoordinate2D)bl_gcj02towgs84:(double)lng andLat:(double)lat {
   if ([self out_of_china:lng andLat:lat ]) {
        CLLocationCoordinate2D coor;
        coor.longitude = lng;
        coor.latitude = lat;
        return coor;
    }
    else {
        double dlat = [self transformlat:lng - 105.0 andLat:lat - 35.0];
        double dlng = [self transformlng:lng - 105.0 andLat:lat - 35.0];
        double radlat = lat / 180.0 * _PI;
        double magic = sin(radlat);
        magic = 1 - _ee * magic * magic;
        double sqrtmagic = sqrt(magic);
        dlat = (dlat * 180.0) / ((_a * (1 - _ee)) / (magic * sqrtmagic) * _PI);
        dlng = (dlng * 180.0) / (_a / sqrtmagic * cos(radlat) * _PI);
        double bl_lat = lat + dlat;
        double bl_lng = lng + dlng;
        CLLocationCoordinate2D coor;
        coor.longitude = bl_lng;
        coor.latitude = bl_lat;
        return coor;
    }
}

- (double)transformlat:(double)lng andLat:(double)lat {
    double ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + 0.1 * lng * lat + 0.2 * sqrt(fabs(lng));
    ret += (20.0 * sin(6.0 * lng * _PI) + 20.0 * sin(2.0 * lng * _PI)) * 2.0 / 3.0;
    ret += (20.0 * sin(lat * _PI) + 40.0 * sin(lat / 3.0 * _PI)) * 2.0 / 3.0;
    ret += (160.0 * sin(lat / 12.0 * _PI) + 320 * sin(lat * _PI / 30.0)) * 2.0 / 3.0;
    return ret;
}

- (double)transformlng:(double)lng andLat:(double)lat{
    double ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * sqrt(fabs(lng));
    ret += (20.0 * sin(6.0 * lng * _PI) + 20.0 * sin(2.0 * lng * _PI)) * 2.0 / 3.0;
    ret += (20.0 * sin(lng * _PI) + 40.0 * sin(lng / 3.0 * _PI)) * 2.0 / 3.0;
    ret += (150.0 * sin(lng / 12.0 * _PI) + 300.0 * sin(lng / 30.0 * _PI)) * 2.0 / 3.0;
    return ret;
}

/**
 * 判断是否在国内,不在国内则不做偏移
 */
- (BOOL)out_of_china:(double)lng andLat:(double)lat {
    return (lng < 72.004 || lng > 137.8347) || ((lat < 0.8293 || lat > 55.8271) || false);
}
-(instancetype)init{
    _x_PI = 3.14159265358979324 * 3000.0 / 180.0;
    _PI = 3.1415926535897932384626;
    _a = 6378245.0;
    _ee = 0.00669342162296594323;
    return self;
}
+ (instancetype)shared {
    static BLCoordinatesChangeHelper *manager = nil;
    static dispatch_once_t onceToken;
    dispatch_once(&onceToken, ^{
        manager = [[BLCoordinatesChangeHelper alloc] init];
    });
    return manager;
}

@end

 

转载请注明出处谢谢:https://my.oschina.net/iceTear/blog/876670

冰泪_
粉丝 16
博文 91
码字总数 64172
作品 0
青浦
程序员
私信 提问
加载中
请先登录后再评论。
Flappy Bird(安卓版)逆向分析(一)

更改每过一关的增长分数 反编译的步骤就不介绍了,我们直接来看反编译得到的文件夹 方法1:在smali目录下,我们看到org/andengine/,可以知晓游戏是由andengine引擎开发的。打开/res/raw/at...

enimey
2014/03/04
6K
18
【opencv】图形的绘制

1.矩形图像的绘制: 原函数:void cvRectangle(CvArr* img, CvPoint pt1, CvPoint pt2, CvScalar color, int thickness=1, int line_type=8,int shift=0) img就是需要绘制的图像 pt1 and pt......

其实我是兔子
2014/10/08
1.2K
1
beego API开发以及自动化文档

beego API开发以及自动化文档 beego1.3版本已经在上个星期发布了,但是还是有很多人不了解如何来进行开发,也是在一步一步的测试中开发,期间QQ群里面很多人都问我如何开发,我的业余时间实在...

astaxie
2014/06/25
2.7W
22
Javascript图元绘制库--ternlight

基于HTML CANVAS API的Javascript库,提供在HTML页面上绘制图元——如流程图的能力。 目前已支持简单的矩形图元和图元间的连线(直线、直角连线两种),拖拽图元等能力。 该javascript librar...

fancimage1
2013/02/07
6.3K
1
StreetPass

StreetPass,最初是想模拟任天堂NDS掌机中的StreetPass开发的(但未完成,大体框架已可行)。适用于记录每天在街上偶遇的纸妹Or帅锅的信息。 原理是想采用Wifi hot技术,由此可以在搜索到附近...

口米巴
2013/03/18
1.6K
0

没有更多内容

加载失败,请刷新页面

加载更多

jQuery Ajax调用后如何管理重定向请求 - How to manage a redirect request after a jQuery Ajax call

问题: I'm using $.post() to call a servlet using Ajax and then using the resulting HTML fragment to replace a div element in the user's current page. 我使用$.post()使用Ajax调用......

javail
46分钟前
15
0
没有指定分支的“git push”的默认行为 - Default behavior of “git push” without a branch specified

问题: I use the following command to push to my remote branch: 我使用以下命令推送到我的远程分支: git push origin sandbox If I say 如果我说 git push origin does that push ch......

技术盛宴
今天
21
0
为什么在允许某些Unicode字符的注释中执行Java代码?

问题: The following code produces the output "Hello World!" 以下代码生成输出“Hello World!” (no really, try it). (不,真的,试试吧)。 public static void main(String... args......

富含淀粉
今天
18
0
什么是按位移位(位移)运算符以及它们如何工作? - What are bitwise shift (bit-shift) operators and how do they work?

问题: I've been attempting to learn C in my spare time, and other languages (C#, Java, etc.) have the same concept (and often the same operators) ... 我一直在尝试在业余时间学习......

法国红酒甜
今天
32
0
OSChina 周二乱弹 —— 卧槽 李荣浩的契约兽啊

Osc乱弹歌单(2020)请戳(这里) 【今日歌曲】 @薛定谔的兄弟 :分享洛神有语创建的歌单「我喜欢的音乐」: 《红色的回忆》- 痛仰乐队 手机党少年们想听歌,请使劲儿戳(这里) 动弹, 又好多...

小小编辑
今天
89
1

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部