【PHP版】火星坐标系 (GCJ02) 与百度坐标系 (BD09ll)转换算法

首先感谢java版作者@宋宋宋伟，java版我是看http://blog.csdn.net/coolypf/article/details/8569813

然后根据java代码修改成了php代码。

    <?php

    /**
     * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法 将 GCJ-02 坐标转换成 BD-09 坐标
     * 
     * @param gg_lat
     * @param gg_lon
     * @return
     */
    function gcj02_To_Bd09($gg_lon, $gg_lat) {
        $x = $gg_lon;
        $y = $gg_lat;
        $z = Math.sqrt($x * $x + $y * $y) + 0.00002 * Math.sin($y * pi());
        $theta = Math.atan2($y, $x) + 0.000003 * Math.cos($x * pi());
        $bd_lon = $z * Math.cos($theta) + 0.0065;
        $bd_lat = $z * Math.sin($theta) + 0.006;
        return array($bd_lon, $bd_lat);
    }

    /**
     * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法   将 BD-09 坐标转换成GCJ-02 坐标 
     * 
     * @param bd_lon
     * @param bd_lat
     * @return
     */
    function bd09_To_Gcj02($bd_lon, $bd_lat) {
        $x = $bd_lon - 0.0065;
        $y = $bd_lat - 0.006;
        $z = sqrt($x * $x + $y * $y) - 0.00002 * sin($y * pi());
        $theta = atan2($y, $x) - 0.000003 * cos($x * pi());
        $gg_lon = $z * cos($theta);
        $gg_lat = $z * sin($theta);
        return array($gg_lon, $gg_lat);
    }

    function coordinate_switch($a,$b){//百度转腾讯坐标转换
      $x = (double)$b - 0.0065;
      $y = (double)$a - 0.006;
      $x_pi = 3.14159265358979324;
      $z = sqrt($x * $x+$y * $y) - 0.00002 * sin($y * $x_pi);
      $theta = atan2($y,$x) - 0.000003 * cos($x*$x_pi);
      $gb = number_format($z * cos($theta),15);
      $ga = number_format($z * sin($theta),15);
      return ['Latitude'=>$ga,'Longitude'=>$gb];
    }
     
    function coordinate_switchf($a,$b){//腾讯转百度坐标转换
      $x = (double)$b ;
      $y = (double)$a;
      $x_pi = 3.14159265358979324;
      $z = sqrt($x * $x+$y * $y) + 0.00002 * sin($y * $x_pi);
      $theta = atan2($y,$x) + 0.000003 * cos($x*$x_pi);
      $gb = number_format($z * cos($theta) + 0.0065,6);
      $ga = number_format($z * sin($theta) + 0.006,6);
     
      return ['Latitude'=>$ga,'Longitude'=>$gb];
     
    }

    // // 113.139278,23.112388
    $bd_lon = 23.112388;
    $bd_lat = 113.139278;
    // print_r(bd09_To_Gcj02($bd_lon, $bd_lat));
    print_r(coordinate_switch($bd_lon, $bd_lat));
    // // 23.106200，113.132840
    $gg_lon = 23.106200;
    $gg_lat = 113.132840;
    // print_r(gcj02_To_Bd09($gg_lon, $gg_lat));
    print_r(coordinate_switchf($gg_lon, $gg_lat));