您好,登錄后才能下訂單哦!
這篇文章主要介紹“php gps如何轉gcj-02”,在日常操作中,相信很多人在php gps如何轉gcj-02問題上存在疑惑,小編查閱了各式資料,整理出簡單好用的操作方法,希望對大家解答”php gps如何轉gcj-02”的疑惑有所幫助!接下來,請跟著小編一起來學習吧!
php gps轉gcj-02的方法:1、創建一個php示例文件;2、通過“public static function wgs84ToGcj02 (float $lng, float $lat):array {...}”方法將WGS84轉為GCJ02即可。
php 經緯度坐標轉換 WGS84、火星坐標 (GCJ-02)、百度坐標 (BD-09)
項目有gps上報的功能, 由于前端插件問題導致大量gps定位數據轉換百度坐標(BD-09)時產生極大偏移, 故需要后端做經緯度坐標轉換, 看到一篇java的相關技術帖, 拿來做了修改
Ps: 坐標轉換稍有偏差,在可接受范圍內
php代碼:
<?php
namespace App\Tool;
/**
* Class GpsUtils
* GCJ-02 -- 由國測局制定的GCJ-02 標準,高德地圖,騰訊地圖,谷歌地圖中國大陸板塊均采用此標準
* @package App\Tool
*/
class GpsUtils {
const x_pi = 3.14159265358979324 * 3000.0 / 180.0;
// π
const pi = 3.1415926535897932384626;
// 長半軸
const a = 6378245.0;
// 扁率
const ee = 0.00669342162296594323;
/**
* 百度坐標系(BD-09)轉WGS坐標
*
* @param float $lng 百度坐標緯度
* @param float $lat 百度坐標經度
* @return array WGS84坐標數組
*/
public static function bd09ToWgs84 (float $lng, float $lat): array {
$gcj = self::bd09ToGcj02($lng, $lat);
return self::gcj02ToWgs84($gcj[0], $gcj[1]);
}
/**
* WGS坐標轉百度坐標系(BD-09)
*
* @param float $lng WGS84坐標系的經度
* @param float $lat WGS84坐標系的緯度
* @return array 百度坐標數組
*/
public static function wgs84ToBd09 (float $lng, float $lat): array {
$gcj = self::wgs84ToGcj02($lng, $lat);
return self::gcj02ToBd09($gcj[0], $gcj[1]);
}
/**
* 火星坐標系(GCJ-02)轉百度坐標系(BD-09)
*
* @param float $lng 火星坐標經度
* @param float $lat 火星坐標緯度
* @return array 百度坐標數組
* @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];
}
/**
* 百度坐標系(BD-09)轉火星坐標系(GCJ-02)
*
* @param float $bd_lon 百度坐標緯度
* @param float $bd_lat 百度坐標經度
* @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轉GCJ02(火星坐標系)
*
* @param float $lng WGS84坐標系的經度
* @param float $lat WGS84坐標系的緯度
* @return array 火星坐標數組
*/
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(火星坐標系)轉GPS84
* @param float $lng 火星坐標系的經度
* @param float $lat 火星坐標系緯度
* @return array WGS84坐標數組
*/
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];
}
/**
* 緯度轉換
* @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;
}
/**
* 經度轉換
* @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;
}
}
相關拓展:
GCJ-02是由中國國家測繪局(G表示Guojia國家,C表示Cehui測繪,J表示Ju局)制訂的地理信息系統的坐標系統。
它是一種對經緯度數據的加密算法,即加入隨機的偏差。
國內出版的各種地圖系統(包括電子形式),必須至少采用GCJ-02對地理位置進行首次加密。
到此,關于“php gps如何轉gcj-02”的學習就結束了,希望能夠解決大家的疑惑。理論與實踐的搭配能更好的幫助大家學習,快去試試吧!若想繼續學習更多相關知識,請繼續關注億速云網站,小編會繼續努力為大家帶來更多實用的文章!
免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。