YiluPHP
这家伙很懒，什么都没有留下...
获赞
32

经验 国际坐标WGS-84、百度地图（BD-09）坐标、高德（GCJ-02）坐标的相互转换（PHP版）

浏览数 186521

简单介绍一下常见的地理经纬度坐标系：
WGS-84坐标系：是国际标准坐标系，GPS坐标（Google Earth使用、或者GPS模块）。

GCJ-02坐标系：火星坐标系，GCJ-02是由中国国家测绘局制订的地理信息系统的坐标系统。它是一种对经纬度数据的加密算法，即加入随机的偏差。国内出版的各种地图系统（包括电子形式），出于国家安全考虑，必须至少采用GCJ-02对地理位置进行首次加密。国内腾讯地图和高德地图都采用了GCJ-02火星坐标系。

BD-09坐标系：百度坐标系，也就是在GCJ-02的基础上，进行了BD-09二次加密措施生成的一种坐标系。

GCJ-02(火星，高德) 坐标和BD-09(百度) 坐标相互转换的函数PHP版本：

//火星坐标(GCJ02坐标，高德，谷歌，腾讯坐标)到百度坐标BD-09

//百度坐标BD-09到火星坐标GCJ02(高德，谷歌，腾讯坐标)

// WGS-84(GPS坐标，谷歌地球坐标) 到 GCJ-02(火星坐标) 的转换


class GpsTransUtil
{
private static $x_pi = 3.14159265358979324*3000.0/180.0;
private static $pi = 3.14159265358979324;
private static $a = 6378137.0; //6378245.0
private static $ee = 0.00669342162296594323;

//火星坐标(GCJ02坐标，高德，谷歌，腾讯坐标)到百度坐标BD-09
public static function gcjTObd($lng, $lat){
$x = $lng;
$y = $lat;
$z = sqrt($x * $x + $y * $y) + 0.00002 * sin($y * self::$x_pi);
$theta = atan2($y, $x) + 0.000003 * cos($x * self::$x_pi);
$to_lng = $z * cos($theta) + 0.0065;
$to_lat = $z * sin($theta) + 0.006;

return ['lng' => $to_lng, 'lat' => $to_lat];
}

//百度坐标BD-09到火星坐标GCJ02(高德，谷歌，腾讯坐标)
public static function bdTOgcj($lng, $lat){
$x = $lat - 0.0065;
$y = $lng - 0.006;
$z = sqrt($x * $x + $y * $y) - 0.00002 * sin($y * self::$x_pi);
$theta = atan2($y, $x) - 0.000003 * cos($x * self::$x_pi);
$to_lng = $z * cos($theta);
$to_lat = $z * sin($theta);
return ['lng' => $to_lng, 'lat' => $to_lat];
}

// WGS-84(GPS坐标，谷歌地球坐标) 到 GCJ-02(火星坐标) 的转换
public static function wgsTOgcj($lng, $lat){
//double wgLat, double wgLon, out double mgLat, out double mgLon
$wgLat = $lat;
$wgLon = $lng;
if (self::outOfChina($wgLat, $wgLon)){
return ['lng' => $wgLon, 'lat' => $wgLat];
}

$dLat = self::transformLat($wgLon - 105.0, $wgLat - 35.0);
$dLon = self::transformLon($wgLon - 105.0, $wgLat - 35.0);
$radLat = $wgLat / 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);
$dLon = ($dLon * 180.0) / (self::$a / $sqrtMagic * cos($radLat) * self::$pi);

$mgLat = $wgLat + $dLat;
$mgLon = $wgLon + $dLon;
return ['lng' => $mgLon, 'lat' => $mgLat];
}

// GCJ-02 到 WGS-84 的转换
public static function gcjTOwgs($lng, $lat)
{
$to = self::wgsTOgcj($lng, $lat);
$g_lat = $to['lat'];
$g_lon = $to['lng'];
$d_lat = $g_lat - $lat;
$d_lon = $g_lon - $lng;
return ['lng' => $lng - $d_lon, 'lat' => $lat - $d_lat];
}

//把地球经纬度转换为lat和lng
public static function earthCoordinate($t1, $t2) {
$tt1 = $t1->x + $t1->y/60.0 + $t1->z/3600.0;
$tt2 = $t2->x + $t2->y/60.0 + $t2->z/3600.0;
return ['lng' => $tt2, 'lat' => $tt1];
}

private static function outOfChina($lat,$lon)
{
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;
}
}
我来说说