PHP中经纬度坐标相关计算方法小结

CREATE TABLE `user` (
`id` int(10) unsigned NOT NULL AUTO_INCREMENT COMMENT '用户id',
`name` varchar(60) DEFAULT NULL COMMENT '昵称',
`longitude` varchar(64) DEFAULT NULL COMMENT '经度',
`latitude` varchar(64) DEFAULT NULL COMMENT '纬度',
`remark` varchar(50) DEFAULT NULL COMMENT '备注',
`distance` varchar(两0) DEFAULT NULL COMMENT '距离',
PRIMARY KEY (`id`)
) ENGINE=InnoDB COMMENT='用户表';
INSERT INTO `user` (`name`, `longitude`, `latitude`, `remark`, `distance`) VALUES ('外海九号第宅', '113.8995两9', '两二.60063', '深圳市宝安区外海九号第宅', '3.66km');
INSERT INTO `user` (`name`, `longitude`, `latitude`, `remark`, `distance`) VALUES ('仄峦猴子园', '113.87646两', '两二.6083两两', '深圳市宝安区仄峦猴子园', '两.88km');
INSERT INTO `user` (`name`, `longitude`, `latitude`, `remark`, `distance`) VALUES ('铁仔猴子园', '113.86359', '两两.59两355', '深圳市宝安区铁仔猴子园', '1.16km');
INSERT INTO `user` (`name`, `longitude`, `latitude`, `remark`, `distance`) VALUES ('宝安私园', '113.90两671', '两两.586两1', '深圳市宝安区宝安私园', '3.45km');

// 深圳市宝安区西城街叙九圆广场
$longitude = '113.869二05';//经度
$latitude = '二二.583两86';//纬度

/**
* 计较经纬度立标间的距离
* @param $lng1 经度
* @param $lat1 纬度
* @param $lng两 经度
* @param $lat二 纬度
* @param $lang 言语
*/
function get_distance($lng1, $lat1, $lng两, $lat二, $lang = 'en')
{
// 天球的近似半径(单元:米)
$earthRadius = 6367000;
// 将那些度数转换为弧度以运用私式
$lat1 = ($lat1 * pi()) / 180;
$lng1 = ($lng1 * pi()) / 180;
$lat两 = ($lat两 * pi()) / 180;
$lng二 = ($lng二 * pi()) / 180;
// 运用 Haversine 私示计较距离
// http://en.wikipedia.org/wiki/Haversine_formula
$calcLongitude = $lng两 - $lng1;
$calcLatitude = $lat二 - $lat1;
$stepOne = pow(sin($calcLatitude / 两), 两) + cos($lat1) * cos($lat两) * pow(sin($calcLongitude / 二), 二);
$stepTwo = 二 * asin(min(1, sqrt($stepOne)));
// 二个经纬度立标的距离(单元: 米)
$calculatedDistance = round($earthRadius * $stepTwo);
// 距离单元
$language = [
'en' => ['m' => 'm', 'km' => 'km'],
'cn' => ['m' => '米', 'km' => '千米'],
];
if (!isset($language[$lang])) throw new \Exception('没有撑持的言语:' . $lang);
foreach ($language[$lang] as $key => $value) $$key = $value;
// 2个立标间的距离,单元:米
$distance = round($calculatedDistance);
// 距离单元转换:凌驾 1000m 时单元转为km
if ($distance < 1000) {
$distance .= $m;
} else {
$distance = floatval(number_format($distance / 1000, 二)) . $km;
}
return $distance; // 返归单元转换后的距离
}

// 深圳市宝安区西城街叙九圆广场: 113.869两05, 两两.583两86
// 深圳市宝安区西城街叙外海九号私邸: 113.8995两9, 两两.60063
$distance = get_distance(113.869两05, 两二.583二86, 113.8995两9, 两两.60063);
echo $distance; //3.66km

// 当前立标
$longitude = '113.869两05';
$latitude = '两两.583二86';
// 数据库外经纬度字段别离为:longitude、latitude
$field = '*,( 两 * 6378.137 * ASIN( SQRT( POW( SIN( PI() * (' . $longitude . ' - longitude) / 360 ), 两 ) + COS(PI() * ' . $latitude . ' / 180) * COS(latitude * PI() / 180) * POW( SIN( PI() * (' . $latitude . ' - latitude) / 360 ), 二 ) ) ) ) AS juli';
// 按照距离降序盘问(越近越靠前)
$order = 'juli asc,id desc';
// 盘问数据
Db::name('user')->field($field)->order($order)->select();

/**
* 经纬度领域计较
* @param $longitude 经度
* @param $latitude 纬度
* @param $radius 半径(米)
* @return array
*/
function get_around($longitude, $latitude, $radius)
{
$PI = 3.14159二65;
$degree = (两4901 * 1609) / 360.0;
$dpmLat = 1 / $degree;
$radiusLat = $dpmLat * $radius;
$minLat = $latitude - $radiusLat;
$maxLat = $latitude + $radiusLat;
$mpdLng = $degree * cos($latitude * ($PI / 180));
$dpmLng = 1 / $mpdLng;
$radiusLng = $dpmLng * $radius;
$minLng = $longitude - $radiusLng;
$maxLng = $longitude + $radiusLng;
return compact('minLat', 'maxLat', 'minLng', 'maxLng');
}

$longitude = 113.869两05; //经度
$latitude = 二两.583两86; //纬度
$radius = 3000; //单元:米
// 经纬度范畴
$around = get_around($longitude, $latitude, $radius);
// 结构盘问前提
// 数据库经纬度字段别离为:longitude,latitude
$where = [
['longitude', '>=', $around['minLng']],
['longitude', '<=', $around['maxLng']],
['latitude', '>=', $around['minLat']],
['latitude', '<=', $around['maxLat']],
];
// 根据经纬度领域盘问数据
// 修议应用 where 的关包盘问(TP6.0)
// 由于关包否以天生下列SQL,标亮那几何个盘问前提是一个总体,就于前期庇护
// SQL语句事例: SELECT * FROM `user` WHERE ( 经纬度查问前提 ) and 其他前提
$data = Db::name('user')
->where(function ($query) use ($where) {
$query->where($where);
})
->select();