经纬坐标转换、判断、计算两点之间距离
WGS84转GCj02(北斗转高德),GCJ02 转换为 WGS84 (高德转北斗),百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换,GCJ-02 转换为 BD-09 （火星坐标系 转百度即谷歌、高德 转 百度）,根据经纬坐标计算两点之间距离,判断坐标点是否在某一个区域中
小编
1、在common文件夹创建CoordinateWGBConvert类
2、CoordinateWGBConvert类内容
<?php
namespace app\common;
class CoordinateWGBConvert
{
const x_PI = 52.35987755982988;
const PI = 3.1415926535897932384626;
// const a = 6378245.0;
const a = 6371393.0;//半径米
const ee = 0.00669342162296594323;
/**
* 根据经纬坐标计算两点之间距离
* @lat1:起始纬度
* @lng1:起始经度
* @lat2:终点纬度
* @lng2:终点经度
* 返回：
* @distance：距离
* @company：单位
* */
public function getDistance($lat1, $lng1, $lat2, $lng2)
{
$radLat1 = deg2rad(floatval($lat1));// deg2rad()函数将角度转换为弧度
$radLat2 = deg2rad(floatval($lat2));
$radLng1 = deg2rad(floatval($lng1));
$radLng2 = deg2rad(floatval($lng2));
$a = $radLat1 - $radLat2;
$b = $radLng1 - $radLng2;
// $s = 2 * asin(sqrt(pow(sin($a / 2), 2) + cos($radLat1) * cos($radLat2) * pow(sin($b / 2), 2))) * 6378.137;
$s = 2 * asin(sqrt(pow(sin($a / 2), 2) + cos($radLat1) * cos($radLat2) * pow(sin($b / 2), 2))) * (self::a / 1000);
if ($s * 1000 < 1000) {
return ["distance" => round($s * 1000, 2), "company" => "m"];
} else {
return ["distance" => round($s, 2), "company" => "km"];
}
}
/**
* 判断坐标点是否在某一个区域中
* @param $aLon 经度或者纬度
* @param $aLat 纬度或者经度
* @param $pointList 区域地点多边形点的顺序需根据顺时针或逆时针，不能乱
* @return bool|int true:范围之内，false:范围之外
*/
function is_ptin_poly($aLon, $aLat, $pointList = array())
{
$iSum = 0;
$iCount = count($pointList);
if ($iCount < 3) {
return false;
}
foreach ($pointList as $key => $row) {
$pLon1 = $row[0];
$pLat1 = $row[1];
if ($key === $iCount - 1) {
$pLon2 = $pointList[0][0];
$pLat2 = $pointList[0][1];
} else {
$pLon2 = $pointList[$key + 1][0];
$pLat2 = $pointList[$key + 1][1];
}
if ((($aLat >= $pLat1) && ($aLat < $pLat2)) || (($aLat >= $pLat2) && ($aLat < $pLat1))) {
if (abs($pLat1 - $pLat2) > 0) {
$pLon = $pLon1 - (($pLon1 - $pLon2) * ($pLat1 - $aLat)) / ($pLat1 - $pLat2);
if ($pLon < $aLon) {
$iSum += 1;
}
}
}
}
if ($iSum % 2 != 0) {
return true;
} else {
return false;
}
}
/**
* WGS84转GCj02(北斗转高德)
* @param lng
* @param lat
* @returns {*[]}
*/
public function wgs84togcj02($lng, $lat) {
if ($this->out_of_china($lng, $lat)) {
return array($lng, $lat);
} else {
$dlat = $this->transformlat($lng - 105.0, $lat - 35.0);
$dlng = $this->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 ["lng"=>$mglng, "lat"=>$mglat];
}
}
/**
* GCJ02 转换为 WGS84 (高德转北斗)
* @param lng
* @param lat
* @return array(lng, lat);
*/
public function gcj02towgs84($lng, $lat) {
if ($this->out_of_china($lng, $lat)) {
return array($lng, $lat);
} else {
$dlat = $this->transformlat($lng - 105.0, $lat - 35.0);
$dlng = $this->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 ["lng"=>$lng * 2 - $mglng, "lat"=>$lat * 2 - $mglat];
}
}
/**
* 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换
* 即 百度 转 谷歌、高德
* @param bd_lon
* @param bd_lat
* @returns
*/
public function bd09togcj02 ($bd_lon, $bd_lat) {
$x_pi = 3.14159265358979324 * 3000.0 / 180.0;
$x = $bd_lon - 0.0065;
$y = $bd_lat - 0.006;
$z = sqrt($x * $x + $y * $y) - 0.00002 * sin($y * $x_pi);
$theta = atan2($y, $x) - 0.000003 * cos($x * $x_pi);
$gg_lng = $z * cos($theta);
$gg_lat = $z * sin($theta);
return ["lng"=>$gg_lng, "lat"=>$gg_lat];
}
/**
* GCJ-02 转换为 BD-09 （火星坐标系 转百度即谷歌、高德 转 百度）
* @param $lng
* @param $lat
* @returns array(bd_lng, bd_lat)
*/
public function gcj02tobd09($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 ["lng"=>$bd_lng, "lat"=>$bd_lat];
}
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;
}
private 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 * 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;
}
private function rad($param)
{
return $param * self::PI / 180.0;
}
/**
* 判断是否在国内，不在国内则不做偏移
* @param $lng
* @param $lat
* @returns {boolean}
*/
private function out_of_china($lng, $lat) {
return ($lng < 72.004 || $lng > 137.8347) || (($lat < 0.8293 || $lat > 55.8271) || false);
}
}