我之前没有接触过地图，还以为只有一套标准。最近也是稍有接触，查看一下资料，原来还有这么七弯八拐的事。
是国际标准坐标系，GPS坐标（Google Earth使用、或者GPS模块）。
火星坐标系，GCJ-02是由中国国家测绘局制订的地理信息系统的坐标系统。它是一种对经纬度数据的加密算法，即加入随机的偏差。国内出版的各种地图系统（包括电子形式），出于国家安全考虑，必须至少采用GCJ-02对地理位置进行首次加密。国内腾讯地图和高德地图都采用了GCJ-02火星坐标系。
百度坐标系，也就是在GCJ-02的基础上，进行了BD-09二次加密措施生成的一种坐标系。
转换源码
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
| <?php
// 命名空间
define('X_PI',3.14159265358979324 * 3000.0 / 180.0);
class Coordinate {
public $x = 0;
public $y = 0;
/**
* Coordinate constructor.
* @param $lon float 经度
* @param $lat float 纬度
*/
public function __construct ($lon,$lat) {
$this->x = $lon;
$this->y = $lat;
}
}
class CoordinateTool {
private static $pi = 3.14159265358979324; // 圆周率
private static $a = 6378245.0; // WGS 长轴半径
private static $ee = 0.00669342162296594323; // WGS 偏心率的平方
/**
* 将火星坐标系GCJ-02 坐标 转换成百度坐标系 BD-09 坐标
* @param $gc_loc 火星坐标点(Class Coordinate)
* @return $bg_loc Coordinate对象，百度地图经纬度坐标
*/
public static function gcj_bd($gc_loc) {
$x_pi = X_PI;
$x = $gc_loc->x;
$y = $gc_loc->y;
$z = sqrt($x * $x + $y * $y) + 0.00002 * sin($y * $x_pi);
$theta = atan2($y, $x) + 0.000003 * cos($x * $x_pi);
$bd_x = $z * cos($theta) + 0.0065;
$bd_y = $z * sin($theta) + 0.006;
$bg_loc = new Coordinate($bd_x,$bd_y);
return $bg_loc;
}
/**
* 将百度坐标系 BD-09 坐标 转换成 火星坐标系GCJ-02 坐标
* @param $bd_loc 火星坐标点(Class Coordinate)
* @return $bg_loc Coordinate对象，火星坐标系经纬度坐标
*/
public static function bd_gcj($bd_loc) {
$x_pi = X_PI;
$x = $bd_loc->x - 0.0065;
$y = $bd_loc->y - 0.006;
$z = sqrt($x * $x + $y * $y) - 0.00002 * sin($y * $x_pi);
$theta = atan2($y, $x) - 0.000003 * cos($x * $x_pi);
$gc_x = $z * cos($theta);
$gc_y = $z * sin($theta);
$gc_loc = new Coordinate($gc_x,$gc_y);
return $gc_loc;
}
/**
* 将国际通用坐标系WGS84坐标 转换成 火星坐标系GCJ-02 坐标
* @param $wgs_loc WGS84坐标点(Class Coordinate)
* @return $bg_loc Coordinate对象，火星坐标系经纬度坐标
*/
public static function wgs_gcj($wgs_loc) {
$wgs_lon = $wgs_loc->x;
$wgs_lat = $wgs_loc->y;
if (self::outOfChina($wgs_lon,$wgs_lat)){
return (new Coordinate($wgs_lon,$wgs_lat));
}
$x_pi = X_PI;
$dLat = self::transformLat($wgs_lon - 105.0, $wgs_lat - 35.0);
$dLon = self::transformLon($wgs_lon - 105.0, $wgs_lat - 35.0);
$radLat = $wgs_lat / 180.0 * $x_pi;
$magic = sin($radLat);
$magic = 1 - self::$ee * $magic * $magic;
$sqrtMagic = sqrt($magic);
$dLat = ($dLat * 180.0) / ((self::$a * (1 - self::$ee)) / ($magic * $sqrtMagic) * self::$pi);
$dLon = ($dLon * 180.0) / (self::$a / $sqrtMagic * cos($radLat) * self::$pi);
$mgLat = $wgs_lat + $dLat;
$mgLon = $wgs_lon + $dLon;
$gcj_loc = new Coordinate($mgLon,$mgLat);
return $gcj_loc;
}
/**
* 将火星坐标系GCJ-02坐标 转换成 国际通用坐标系WGS84坐标
* @param $gcj_loc GCJ-02坐标点(Class Coordinate)
* @return $wgs_loc Coordinate对象，国际通用坐标系WGS84坐标
*/
public static function gcj_wgs($gcj_loc) {
$to = self::wgs_gcj($gcj_loc);
$lon = $gcj_loc->x;
$lat = $gcj_loc->y;
$gcl_lon = $to->x;
$gcl_lat = $to->y;
$d_lat = $gcl_lat - $lat;
$d_lon = $gcl_lon - $lon;
return new Coordinate($lon - $d_lon,$lat - $d_lat);
}
private static function outOfChina($lon,$lat) {
if ($lon < 72.004 || $lon > 137.8347)
return true;
if ($lat < 0.8293 || $lat > 55.8271)
return true;
return false;
}
private static function transformLat($x,$y) {
$ret = -100.0 + 2.0 * $x + 3.0 * $y + 0.2 * $y * $y + 0.1 * $x * $y + 0.2 * sqrt(abs($x));
$ret += (20.0 * sin(6.0 * $x * self::$pi) + 20.0 * sin(2.0 * $x * self::$pi)) * 2.0 / 3.0;
$ret += (20.0 * sin($y * self::$pi) + 40.0 * sin($y / 3.0 * self::$pi)) * 2.0 / 3.0;
$ret += (160.0 * sin($y / 12.0 * self::$pi) + 320 * sin($y * self::$pi / 30.0)) * 2.0 / 3.0;
return $ret;
}
private static function transformLon($x, $y) {
$ret = 300.0 + $x + 2.0 * $y + 0.1 * $x * $x + 0.1 * $x * $y + 0.1 * sqrt(abs($x));
$ret += (20.0 * sin(6.0 * $x * self::$pi) + 20.0 * sin(2.0 * $x * self::$pi)) * 2.0 / 3.0;
$ret += (20.0 * sin($x * self::$pi) + 40.0 * sin($x / 3.0 * self::$pi)) * 2.0 / 3.0;
$ret += (150.0 * sin($x / 12.0 * self::$pi) + 300.0 * sin($x / 30.0 * self::$pi)) * 2.0 / 3.0;
return $ret;
}
}
|
例子
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
| <?php
include 'CoordinateTool.php';
$coordinate = new Coordinate("104.07007156866", "30.54792635074");
$coord = CoordinateTool::wgs_gcj($coordinate);
$coord2 = CoordinateTool::gcj_wgs($coord);
$coord3 = CoordinateTool::gcj_bd($coordinate);
$coord4 = CoordinateTool::bd_gcj($coord3);
echo "<pre>";
var_dump($coord);
var_dump($coord2);
var_dump($coord3);
var_dump($coord4);
?>
|