PHP 不同地圖坐标系經緯度轉換 GCj02 WGS84 BD-09

項目中用到了 北鬥車輛導航，app用的高德插件，因為高德(GCj02坐标系)與北鬥(WGS84 坐标系)采用的坐标系不同，就在網上找了轉換方法，用PHP寫了一遍
           

代碼裡面還有百度的 BD-09 坐标，關于各種地圖坐标的坐标系大家可以自己查資料了解。
           

誤差測試
           

1、北鬥轉高德


 源資料 ：117.093457,39.077387 
 高德轉換 117.099876,39.078529
 本地轉換 117.099884,39.078533


2、高德轉北鬥，再轉高德


 源資料 ：				117.095563,39.077410  高德擷取的高新區發展六道坐标
 本地轉換成北鬥坐标：   117.089142,39.076267
 本地轉換成高德坐标：	117.095565,39.077412


  源資料 ：				117.087199,39.090175	117大廈
  本地轉換成北鬥坐标：  117.080784,39.089027
  本地轉換成高德坐标：	117.087198,39.090174
           

3、天津 117大廈坐标經過本地轉成北鬥，再轉回高德後的誤差
           

下面是代碼：
           

備注：因為項目中隻用到了北鬥與高德之前的經緯度轉換，是以對于其他的轉換函數未做測試，如果要用代碼的話請自己測試誤差大小
           

<?php
class Map
{
    const x_PI  = 52.35987755982988;
    const PI  = 3.1415926535897932384626;
    const a = 6378245.0;
    const ee = 0.00669342162296594323;


    /**
     * WGS84轉GCj02(北鬥轉高德)
     * @param lng
     * @param lat
     * @returns {*[]}
     */
    public function wgs84togcj02($lng,  $lat) {
        if ($this->out_of_china($lng, $lat)) {
            return array($lng, $lat);
        } else {
            $dlat = $this->transformlat($lng - 105.0, $lat - 35.0);
            $dlng = $this->transformlng($lng - 105.0, $lat - 35.0);
            $radlat = $lat / 180.0 * self::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);
            $dlng = ($dlng * 180.0) / (self::a / $sqrtmagic * cos($radlat) * self::PI);
            $mglat = $lat + $dlat;
            $mglng = $lng + $dlng;
            return array($mglng, $mglat);
        }
    }
    /**
     * GCJ02 轉換為 WGS84 (高德轉北鬥)
     * @param lng
     * @param lat
     * @return array(lng, lat);
     */
    public function gcj02towgs84($lng, $lat) {
        if ($this->out_of_china($lng, $lat)) {
            return array($lng, $lat);
        } else {
            $dlat = $this->transformlat($lng - 105.0, $lat - 35.0);
            $dlng = $this->transformlng($lng - 105.0, $lat - 35.0);
            $radlat = $lat / 180.0 * self::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);
            $dlng = ($dlng * 180.0) / (self::a / $sqrtmagic * cos($radlat) * self::PI);
            $mglat = $lat + $dlat;
            $mglng = $lng + $dlng;
            return array($lng * 2 - $mglng, $lat * 2 - $mglat);
        }
    }


        /**
    　　* 百度坐标系 (BD-09) 與 火星坐标系 (GCJ-02)的轉換
    　　* 即 百度 轉 谷歌、高德
    　　* @param bd_lon
    　　* @param bd_lat
    　　* @returns
　　*/
        public function bd09togcj02 ($bd_lon, $bd_lat) {
            $x_pi = 3.14159265358979324 * 3000.0 / 180.0;
            $x = $bd_lon - 0.0065;
            $y = $bd_lat - 0.006;
            $z = sqrt($x * $x + $y * $y) - 0.00002 * sin($y * $x_pi);
            $theta = atan2($y, $x) - 0.000003 * cos($x * $x_pi);
            $gg_lng = $z * cos(theta);
            $gg_lat = $z * sin(theta);
            return array($gg_lng, $gg_lat);
        }

    /**
    * GCJ-02 轉換為 BD-09  （火星坐标系 轉百度即谷歌、高德 轉 百度）
    * @param $lng
    * @param $lat
    * @returns array(bd_lng, bd_lat)
    */
    public function gcj02tobd09($lng, $lat) { 
        $z = sqrt($lng * $lng + $lat * $lat) + 0.00002 * Math.sin($lat * x_PI);
        $theta = Math.atan2($lat, $lng) + 0.000003 * Math.cos($lng * x_PI);
        $bd_lng = $z * cos($theta) + 0.0065;
        $bd_lat = z * sin($theta) + 0.006;
        return array($bd_lng, $bd_lat);
    }


    private function transformlat($lng, $lat) {
        $ret = -100.0 + 2.0 * $lng + 3.0 * $lat + 0.2 * $lat * $lat + 0.1 * $lng * $lat + 0.2 * sqrt(abs($lng));
        $ret += (20.0 * sin(6.0 * $lng * self::PI) + 20.0 * sin(2.0 * $lng * self::PI)) * 2.0 / 3.0;
        $ret += (20.0 * sin($lat * self::PI) + 40.0 * sin($lat / 3.0 * self::PI)) * 2.0 / 3.0;
        $ret += (160.0 * sin($lat / 12.0 * self::PI) + 320 * sin($lat * self::PI / 30.0)) * 2.0 / 3.0;
        return $ret;
    }
    private function transformlng($lng, $lat) {
        $ret = 300.0 + $lng + 2.0 * $lat + 0.1 * $lng * $lng + 0.1 * $lng * $lat + 0.1 * sqrt(abs($lng));
        $ret += (20.0 * sin(6.0 * $lng * self::PI) + 20.0 * sin(2.0 * $lng * self::PI)) * 2.0 / 3.0;
        $ret += (20.0 * sin($lng * self::PI) + 40.0 * sin($lng / 3.0 * self::PI)) * 2.0 / 3.0;
        $ret += (150.0 * sin($lng / 12.0 * self::PI) + 300.0 * sin($lng / 30.0 * self::PI)) * 2.0 / 3.0;
        return $ret;
    }


    private function rad($param)
    {
      return  $param * self::PI / 180.0;
    }
    /**
    * 判斷是否在國内，不在國内則不做偏移
    * @param $lng
    * @param $lat
    * @returns {boolean}
    */
    private function out_of_china($lng, $lat) { 
        return ($lng < 72.004 || $lng > 137.8347) || (($lat < 0.8293 || $lat > 55.8271) || false);
    }

}

$test = new Map();
var_dump($test->wgs84togcj02(117.087199,39.090175));

?>