WGS84,GCJ-02,BD-09坐标系间的经纬度坐标转换，PHP版

//高德GCJ02    天地图WGS84
define('xPI',(3.14159265358979324 * 3000.0) / 180.0);
define('PI',3.1415926535897932384626);
define('a',6378245.0);
define('ee',0.00669342162296594323);

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 * 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;
}

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 * 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;
}
/**
 * 判断是否在国内，不在国内则不做偏移
 * @param $lng
 * @param $lat
 * @returns {boolean}
 */
function outOfChina($lng, $lat) {
    return $lng < 72.004 || $lng > 137.8347 || $lat < 0.8293 || $lat > 55.8271 || false;
}

//高德GCJ02    天地图WGS84

/**
 * GCJ02 转换为 WGS84
 * @param $lng
 * @param $lat
 * @returns {*[]}
 */
function gcj02towgs84($lng, $lat) {
    if (outOfChina($lng, $lat)) {
        return [$lng, $lat];
    } else {
        $dlat = transformlat($lng - 105.0, $lat - 35.0);
        $dlng = transformlng($lng - 105.0, $lat - 35.0);
        $radlat = ($lat / 180.0) * PI;
        $magic = sin($radlat);
        $magic = 1 - ee * $magic * $magic;
        $sqrtmagic = sqrt($magic);
        $dlat = ($dlat * 180.0) / (((a * (1 - ee)) / ($magic * $sqrtmagic)) * PI);
        $dlng = ($dlng * 180.0) / ((a / $sqrtmagic) * cos($radlat) * PI);
        $mglat = $lat + $dlat;
        $mglng = $lng + $dlng;
        return [$lng * 2 - $mglng, $lat * 2 - $mglat];
    }
}

/**
 * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换
 * 即 百度 转 谷歌、高德
 * @param $bdlng
 * @param $bdLat
 * @returns {*[]}
 */
function bd09togcj02($bdlng, $bdLat) {
    $x = $bdlng - 0.0065;
    $y = $bdLat - 0.006;
    $z = sqrt($x * $x + $y * $y) - 0.00002 * sin($y * xPI);
    $theta = atan2($y, $x) - 0.000003 * cos($x * xPI);
    $ggLng = $z * cos($theta);
    $ggLat = $z * sin($theta);
    return [$ggLng, $ggLat];
}

/**
 * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换
 * 即谷歌、高德 转 百度
 * @param $lng
 * @param $lat
 * @returns {*[]}
 */
function gcj02tobd09($lng, $lat) {
    $z = sqrt($lng * $lng + $lat * $lat) + 0.00002 * sin($lat * xPI);
    $theta = atan2($lat, $lng) + 0.000003 * cos($lng * xPI);
    $bdLng = $z * cos($theta) + 0.0065;
    $bdLat = $z * sin($theta) + 0.006;
    return [$bdLng, $bdLat];
}

/**
 * WGS84转GCj02
 * @param $lng
 * @param $lat
 * @returns {*[]}
 */
function wgs84togcj02($lng, $lat) {
    if (outOfChina($lng, $lat)) {
        return [$lng, $lat];
    } else {
        $dlat = transformlat($lng - 105.0, $lat - 35.0);
        $dlng = transformlng($lng - 105.0, $lat - 35.0);
        $radlat = ($lat / 180.0) * PI;
        $magic = sin($radlat);
        $magic = 1 - ee * $magic * $magic;
        $sqrtmagic = sqrt($magic);
        $dlat = ($dlat * 180.0) / (((a * (1 - ee)) / ($magic * $sqrtmagic)) * PI);
        $dlng = ($dlng * 180.0) / ((a / $sqrtmagic) * cos($radlat) * PI);
        $mglat = $lat + $dlat;
        $mglng = $lng + $dlng;
        return [$mglng, $mglat];
    }
}