php gps轉(zhuǎn)gcj-02的方法:1、創(chuàng)建一個(gè)php示例文件;2、通過“public static function wgs84ToGcj02 (float $lng, float $lat):array {…}”方法將WGS84轉(zhuǎn)為GCJ02即可。
本教程操作環(huán)境:windows10系統(tǒng)、PHP8.1版、DELL G3電腦
php gps怎么轉(zhuǎn)gcj-02?
php 經(jīng)緯度坐標(biāo)轉(zhuǎn)換 WGS84、火星坐標(biāo) (GCJ-02)、百度坐標(biāo) (BD-09)
項(xiàng)目有g(shù)ps上報(bào)的功能, 由于前端插件問題導(dǎo)致大量gps定位數(shù)據(jù)轉(zhuǎn)換百度坐標(biāo)(BD-09)時(shí)產(chǎn)生極大偏移, 故需要后端做經(jīng)緯度坐標(biāo)轉(zhuǎn)換, 看到一篇java的相關(guān)技術(shù)帖, 拿來做了修改
Ps: 坐標(biāo)轉(zhuǎn)換稍有偏差,在可接受范圍內(nèi)
百度經(jīng)緯度糾偏api: http://api.map.baidu.com/ag/coord/convert
php代碼:
<?php namespace AppTool; /** * Class GpsUtils * GCJ-02 -- 由國測局制定的GCJ-02 標(biāo)準(zhǔn),高德地圖,騰訊地圖,谷歌地圖中國大陸板塊均采用此標(biāo)準(zhǔn) * @package AppTool */ class GpsUtils { const x_pi = 3.14159265358979324 * 3000.0 / 180.0; // π const pi = 3.1415926535897932384626; // 長半軸 const a = 6378245.0; // 扁率 const ee = 0.00669342162296594323; /** * 百度坐標(biāo)系(BD-09)轉(zhuǎn)WGS坐標(biāo) * * @param float $lng 百度坐標(biāo)緯度 * @param float $lat 百度坐標(biāo)經(jīng)度 * @return array WGS84坐標(biāo)數(shù)組 */ public static function bd09ToWgs84 (float $lng, float $lat): array { $gcj = self::bd09ToGcj02($lng, $lat); return self::gcj02ToWgs84($gcj[0], $gcj[1]); } /** * WGS坐標(biāo)轉(zhuǎn)百度坐標(biāo)系(BD-09) * * @param float $lng WGS84坐標(biāo)系的經(jīng)度 * @param float $lat WGS84坐標(biāo)系的緯度 * @return array 百度坐標(biāo)數(shù)組 */ public static function wgs84ToBd09 (float $lng, float $lat): array { $gcj = self::wgs84ToGcj02($lng, $lat); return self::gcj02ToBd09($gcj[0], $gcj[1]); } /** * 火星坐標(biāo)系(GCJ-02)轉(zhuǎn)百度坐標(biāo)系(BD-09) * * @param float $lng 火星坐標(biāo)經(jīng)度 * @param float $lat 火星坐標(biāo)緯度 * @return array 百度坐標(biāo)數(shù)組 * @see 谷歌、高德——>百度 */ public static function gcj02ToBd09 (float $lng, float $lat): array { $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 [$bd_lng, $bd_lat]; } /** * 百度坐標(biāo)系(BD-09)轉(zhuǎn)火星坐標(biāo)系(GCJ-02) * * @param float $bd_lon 百度坐標(biāo)緯度 * @param float $bd_lat 百度坐標(biāo)經(jīng)度 * @return array * @see 百度——>谷歌、高德 */ public static function bd09ToGcj02 (float $bd_lon, float $bd_lat): array { $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 [$gg_lng, $gg_lat]; } /** * WGS84轉(zhuǎn)GCJ02(火星坐標(biāo)系) * * @param float $lng WGS84坐標(biāo)系的經(jīng)度 * @param float $lat WGS84坐標(biāo)系的緯度 * @return array 火星坐標(biāo)數(shù)組 */ public static function wgs84ToGcj02 (float $lng, float $lat): array { $d_lat = self::transformlat($lng - 105.0, $lat - 35.0); $d_lng = self::transformlng($lng - 105.0, $lat - 35.0); $rad_lat = $lat / 180.0 * self::pi; $magic = sin($rad_lat); $magic = 1 - self::ee * $magic * $magic; $sqrt_magic = sqrt($magic); $d_lat = ($d_lat * 180.0) / ((self::a * (1 - self::ee)) / ($magic * $sqrt_magic) * self::pi); $d_lng = ($d_lng * 180.0) / (self::a / $sqrt_magic * cos($rad_lat) * self::pi); $mg_lat = $lat + $d_lat; $mg_lng = $lng + $d_lng; return [$mg_lng, $mg_lat]; } /** * GCJ02(火星坐標(biāo)系)轉(zhuǎn)GPS84 * @param float $lng 火星坐標(biāo)系的經(jīng)度 * @param float $lat 火星坐標(biāo)系緯度 * @return array WGS84坐標(biāo)數(shù)組 */ public static function gcj02ToWgs84 (float $lng, float $lat): array { $d_lat = self::transformlat($lng - 105.0, $lat - 35.0); $d_lng = self::transformlng($lng - 105.0, $lat - 35.0); $rad_lat = $lat / 180.0 * self::pi; $magic = sin($rad_lat); $magic = 1 - self::ee * $magic * $magic; $sqrt_magic = sqrt($magic); $d_lat = ($d_lat * 180.0) / ((self::a * (1 - self::ee)) / ($magic * $sqrt_magic) * self::pi); $d_lng = ($d_lng * 180.0) / (self::a / $sqrt_magic * cos($rad_lat) * self::pi); $mg_lat = $lat + $d_lat; $mg_lng = $lng + $d_lng; return [$lng * 2 - $mg_lng, $lat * 2 - $mg_lat]; } /** * 緯度轉(zhuǎn)換 * @param float $lng * @param float $lat * @return float|int */ public static function transFormLat (float $lng, float $lat): float { $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; } /** * 經(jīng)度轉(zhuǎn)換 * @param float $lng * @param float $lat * @return float */ public static function transFormLng (float $lng, float $lat): float { $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; } }
相關(guān)拓展:
GCJ-02是由中國國家測繪局(G表示Guojia國家,C表示Cehui測繪,J表示Ju局)制訂的地理信息系統(tǒng)的坐標(biāo)系統(tǒng)。
它是一種對(duì)經(jīng)緯度數(shù)據(jù)的加密算法,即加入隨機(jī)的偏差。
國內(nèi)出版的各種地圖系統(tǒng)(包括電子形式),必須至少采用GCJ-02對(duì)地理位置進(jìn)行首次加密。
推薦學(xué)習(xí):《PHP視頻教程》