PHP 经纬度坐标相关计算方法
-
1. 前言
-
2. 计算经纬度坐标间的距离
-
3. 根据经纬度坐标距离排序
-
4. 经纬度范围查询
1. 前言
想要测试本文提供的几个功能函数，可以使用下面这个数据表结构及其数据
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(20) DEFAULT NULL COMMENT '距离', PRIMARY KEY (`id`)) ENGINE=InnoDB COMMENT='用户表';INSERT INTO `user` (`name`, `longitude`, `latitude`, `remark`, `distance`) VALUES ('中海九号公馆', '113.899529', '22.60063', '深圳市宝安区中海九号公馆', '3.66km');INSERT INTO `user` (`name`, `longitude`, `latitude`, `remark`, `distance`) VALUES ('平峦山公园', '113.876462', '22.608322', '深圳市宝安区平峦山公园', '2.88km');INSERT INTO `user` (`name`, `longitude`, `latitude`, `remark`, `distance`) VALUES ('铁仔山公园', '113.86359', '22.592355', '深圳市宝安区铁仔山公园', '1.16km');INSERT INTO `user` (`name`, `longitude`, `latitude`, `remark`, `distance`) VALUES ('宝安公园', '113.902671', '22.58621', '深圳市宝安区宝安公园', '3.45km');
本文内容测试各个功能函数时，使用的当前位置坐标均为:
// 深圳市宝安区西乡街道九方广场$longitude = '113.869205';//经度$latitude = '22.583286';//纬度
2. 计算经纬度坐标间的距离
计算经纬度坐标间的距离 功能函数 (前四个参数为两组经纬度坐标)
/** * 计算经纬度坐标间的距离 * @param $lng1 经度 * @param $lat1 纬度 * @param $lng2 经度 * @param $lat2 纬度 * @param $lang 语言 */function get_distance($lng1, $lat1, $lng2, $lat2, $lang = 'en'){ // 地球的近似半径(单位:米) $earthRadius = 6367000; // 将这些度数转换为弧度以使用公式 $lat1 = ($lat1 * pi()) / 180; $lng1 = ($lng1 * pi()) / 180; $lat2 = ($lat2 * pi()) / 180; $lng2 = ($lng2 * pi()) / 180; // 使用 Haversine 公示计算距离 // http://en.wikipedia.org/wiki/Haversine_formula $calcLongitude = $lng2 - $lng1; $calcLatitude = $lat2 - $lat1; $stepOne = pow(sin($calcLatitude / 2), 2) + cos($lat1) * cos($lat2) * pow(sin($calcLongitude / 2), 2); $stepTwo = 2 * 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; // 两个坐标间的距离,单位:米 $distance = round($calculatedDistance); // 距离单位转换:超出 1000m 时单位转为km if ($distance < 1000) { $distance .= $m; } else { $distance = floatval(number_format($distance / 1000, 2)) . $km; } return $distance; // 返回单位转换后的距离}
使用示例:
我在 九方广场，手机上的高德地图导航至 中海九号公馆 显示的距离为 3.6公里，计算结果还是很准确的
// 深圳市宝安区西乡街道九方广场: 113.869205, 22.583286// 深圳市宝安区西乡街道中海九号公馆: 113.899529, 22.60063$distance = get_distance(113.869205, 22.583286, 113.899529, 22.60063);echo $distance; //3.66km
3. 根据经纬度坐标距离排序
项目中经常有距离显示数据的场景，根据距离排序，越近越靠前显示；比如: 店铺地址、房源信息等。代码示例:
// 当前坐标$longitude = '113.869205';$latitude = '22.583286';// 数据库中经纬度字段分别为:longitude、latitude$field = '*,( 2 * 6378.137 * ASIN( SQRT( POW( SIN( PI() * (' . $longitude . ' - longitude) / 360 ), 2 ) + COS(PI() * ' . $latitude . ' / 180) * COS(latitude * PI() / 180) * POW( SIN( PI() * (' . $latitude . ' - latitude) / 360 ), 2 ) ) ) ) AS juli';// 根据距离升序查询(越近越靠前)$order = 'juli asc,id desc';// 查询数据Db::name('user')->field($field)->order($order)->select();
4. 经纬度范围查询
经纬度范围计算 功能函数
/** * 经纬度范围计算 * @param $longitude 经度 * @param $latitude 纬度 * @param $radius 半径(米) * @return array */function get_around($longitude, $latitude, $radius){ $PI = 3.14159265; $degree = (24901 * 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');}
使用示例
查询 3 公里内的数据。首先，根据当前位置获取 3 公里内的经纬度范围，然后带上查询条件查询数据库即可
$longitude = 113.869205; //经度$latitude = 22.583286; //纬度$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();
