wgs84转gcj02坐标php类库实现方法

admin

类库代码：

1. <?php
   
2. 
   
3. namespace addons\eshop\library;
   
4. 
   
5. class Zbzh
   
6. {
   
7.     /**
   
8.      * WGS84 转 GCJ02 (国测局坐标，适用于高德、腾讯地图)
   
9.      * @param float $lat  WGS84 纬度
   
10.      * @param float $lng  WGS84 经度
    
11.      * @return array  ['lat'=>GCJ纬度, 'lng'=>GCJ经度]
    
12.      */
    
13.     public function wgs84_to_gcj02($lat, $lng)
    
14.     {
    
15.         // 常量定义
    
16.         $pi = 3.1415926535897932384626;
    
17.         $a = 6378245.0;
    
18.         $ee = 0.00669342162296594323;
    
19.    
    
20.         // 判断是否在国内，不在国内不做偏移
    
21.         if (!$this->is_in_china($lat, $lng)) {
    
22.             return ['lat' => $lat, 'lng' => $lng];
    
23.         }
    
24.    
    
25.         $dlat = $this->transform_lat($lng - 105.0, $lat - 35.0);
    
26.         $dlng = $this->transform_lng($lng - 105.0, $lat - 35.0);
    
27.         $radlat = $lat / 180.0 * $pi;
    
28.         $magic = sin($radlat);
    
29.         $magic = 1 - $ee * $magic * $magic;
    
30.         $sqrtmagic = sqrt($magic);
    
31.    
    
32.         $dlat = ($dlat * 180.0) / (($a * (1 - $ee)) / ($magic * $sqrtmagic) * $pi);
    
33.         $dlng = ($dlng * 180.0) / ($a / $sqrtmagic * cos($radlat) * $pi);
    
34.    
    
35.         $gcj_lat = $lat + $dlat;
    
36.         $gcj_lng = $lng + $dlng;
    
37.    
    
38.         return [
    
39.             'lat' => round($gcj_lat, 6),
    
40.             'lng' => round($gcj_lng, 6)
    
41.         ];
    
42.     }
    
43.    
    
44.     // 纬度转换
    
45.     public function transform_lat($x, $y)
    
46.     {
    
47.         $pi = 3.1415926535897932384626;
    
48.         $ret = -100.0 + 2.0 * $x + 3.0 * $y + 0.2 * $y * $y + 0.1 * $x * $y + 0.2 * sqrt(abs($x));
    
49.         $ret += (20.0 * sin(6.0 * $x * $pi) + 20.0 * sin(2.0 * $x * $pi)) * 2.0 / 3.0;
    
50.         $ret += (20.0 * sin($y * $pi) + 40.0 * sin($y / 3.0 * $pi)) * 2.0 / 3.0;
    
51.         $ret += (160.0 * sin($y / 12.0 * $pi) + 320.0 * sin($y * $pi / 30.0)) * 2.0 / 3.0;
    
52.         return $ret;
    
53.     }
    
54.    
    
55.     // 经度转换
    
56.     public function transform_lng($x, $y)
    
57.     {
    
58.         $pi = 3.1415926535897932384626;
    
59.         $ret = 300.0 + $x + 2.0 * $y + 0.1 * $x * $x + 0.1 * $x * $y + 0.1 * sqrt(abs($x));
    
60.         $ret += (20.0 * sin(6.0 * $x * $pi) + 20.0 * sin(2.0 * $x * $pi)) * 2.0 / 3.0;
    
61.         $ret += (20.0 * sin($x * $pi) + 40.0 * sin($x / 3.0 * $pi)) * 2.0 / 3.0;
    
62.         $ret += (150.0 * sin($x / 12.0 * $pi) + 300.0 * sin($x / 30.0 * $pi)) * 2.0 / 3.0;
    
63.         return $ret;
    
64.     }
    
65.    
    
66.     // 判断是否在中国境内
    
67.     public function is_in_china($lat, $lng)
    
68.     {
    
69.         return ($lng > 73.666666 && $lng < 135.055555 && $lat > 3.866666 && $lat < 53.555555);
    
70.     }
    
71. }

复制代码

方法调用：
在需要的页面use类库，比如：
use addons\eshop\library\Zbzh;

调用代码：

1. // 实例化坐标类
   
2.                     $zbzh = new Zbzh();
   
3.                     $latlng = $zbzh->wgs84_to_gcj02($qxlogs->lat, $qxlogs->lng);
   
4.                     $at = $latlng['lat'];
   
5.                     $lng = $latlng['lng'];

复制代码

就可以了。