class Gps { const x_PI = 52.35987755982988; // π * 3000.0 / 180.0 这里π = 3.14159265358979324 const PI = 3.1415926535897932384626; // π const a = 6378245.0; // WGS长轴半径 const ee = 0.00669342162296594323; // WGS偏心率的平方 /** * WGS84转GCj02(北斗转高德) * @param int|string $lng 经度 * @param int|string $lat 纬度 * @return array */ static public function beiDouToAMap($lng = 0, $lat = 0) { if (self::outOfChina($lng, $lat)) { return array($lng, $lat); } else { $dlat = self::transFormLat($lng - 105.0, $lat - 35.0); $dlng = self::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 [$mglng, $mglat, 'lng' => $mglng, 'lat' => $mglat]; } } /** * GCJ02 转换为 WGS84 (高德转北斗) * @param int|string $lng 经度 * @param int|string $lat 纬度 * @return array */ static public function aMapToBeiDou($lng = 0, $lat = 0) { if (self::outOfChina($lng, $lat)) { return array($lng, $lat); } else { $dlat = self::transFormLat($lng - 105.0, $lat - 35.0); $dlng = self::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 string $bd_lon 百度经度 * @param string $bd_lat 百度纬度 * @return array */ static public function baiDuToAMap($bd_lon = '', $bd_lat = '') { $x = $bd_lon - 0.0065; $y = $bd_lat - 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); $gg_lng = $z * cos($theta); $gg_lat = $z * sin($theta); return array($gg_lng, $gg_lat); } /** * GCJ-02 转换为 BD-09 (火星坐标系转百度 即谷歌、高德转百度) * @param string $lng 经度 * @param string $lat 纬度 * @return array */ static public function aMapToBaiDu($lng = '', $lat = '') { $z = sqrt($lng * $lng + $lat * $lat) + 0.00002 * sin($lat * self::x_PI); $theta = atan2($lat, $lng) + 0.000003 * cos($lng * self::x_PI); $bd_lng = $z * cos($theta) + 0.0065; $bd_lat = $z * sin($theta) + 0.006; return array($bd_lng, $bd_lat); } static 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; } static private function transFormLng($lng = 0, $lat = 0) { $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; } /** * 判断是否在国内,不在国内不做偏移 * @param int $lng 经度 * @param int $lat 纬度 * @return bool */ static private function outOfChina($lng = 0, $lat = 0) { return ($lng < 72.004 || $lng > 137.8347) || (($lat < 0.8293 || $lat > 55.8271) || false); } }
留个备份。
评论
评论功能已经关闭!