mysql getdistance,mysql 根据经纬度和距离筛选数据

先根据经纬度和距离计算一个矩形,然后筛选经纬度在矩形范围内的数据。

/**

*

* @param Select $select

* @param string $columnName

* @param string $longitude

* @param string $latitude

* @param int $distance

* metre

* @return string

*/

static function ST_MBRWithinRectangle(Select $select, $columnName,

$longitude, $latitude, $distance)

{

list ($min, $max) = self::rectangle($longitude, $latitude, $distance);

$where = "$columnName!=PointFromText('Point(0 0)') and $columnName!='' and";

$where .= " ST_MBRWithin(`$columnName`,ST_GeomFromText('Polygon(($min[0] $min[1],";

$where .= "$min[0] $max[1],$max[0] $max[1],$max[0] $min[1],$min[0] $min[1]))'))";

$select->where($where);

}

/**

* get rectangle longitude and latitude

*

* @param float $longitude

* @param float $latitude

* @param integer $distance

* meter

* @return array

*/

static function rectangle($longitude, $latitude, $distance)

{

$radius = 6371 * 1000;

// latitude boundaries

$maxlat = $latitude + rad2deg($distance / $radius);

$minlat = $latitude - rad2deg($distance / $radius);

// longitude boundaries (longitude gets smaller when latitude increases)

$maxlng = $longitude +

rad2deg($distance / $radius / cos(deg2rad($latitude)));

$minlng = $longitude -

rad2deg($distance / $radius / cos(deg2rad($latitude)));

return array(

array(

'lng' => $minlng,

'lat' => $minlat

),

array(

'lng' => $maxlng,

'lat' => $maxlat

)

);

}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
计算距离可以使用 Haversine 公式,计算方向角可以使用 Vincenty 公式。下面是 Java 代码示例: ```java import java.lang.Math; public class LocationUtils { /** * 计算两个经纬度之间的距离(单位:米) * @param lat1 第一个点的纬度 * @param lon1 第一个点的经度 * @param lat2 第二个点的纬度 * @param lon2 第二个点的经度 * @return 两个经纬度之间的距离(单位:米) */ public static double getDistance(double lat1, double lon1, double lat2, double lon2) { double R = 6371000; // 地球半径(米) double dLat = Math.toRadians(lat2 - lat1); double dLon = Math.toRadians(lon2 - lon1); double a = Math.sin(dLat / 2) * Math.sin(dLat / 2) + Math.cos(Math.toRadians(lat1)) * Math.cos(Math.toRadians(lat2)) * Math.sin(dLon / 2) * Math.sin(dLon / 2); double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a)); return R * c; } /** * 计算从第一个点到第二个点的方向角(单位:度) * @param lat1 第一个点的纬度 * @param lon1 第一个点的经度 * @param lat2 第二个点的纬度 * @param lon2 第二个点的经度 * @return 从第一个点到第二个点的方向角(单位:度) */ public static double getAzimuth(double lat1, double lon1, double lat2, double lon2) { double dLon = Math.toRadians(lon2 - lon1); double y = Math.sin(dLon) * Math.cos(Math.toRadians(lat2)); double x = Math.cos(Math.toRadians(lat1)) * Math.sin(Math.toRadians(lat2)) - Math.sin(Math.toRadians(lat1)) * Math.cos(Math.toRadians(lat2)) * Math.cos(dLon); double azimuth = Math.toDegrees(Math.atan2(y, x)); return azimuth >= 0 ? azimuth : azimuth + 360; // 范围:0~360度 } } ``` 使用示例: ```java double distance = LocationUtils.getDistance(39.9042, 116.4074, 31.2304, 121.4737); double azimuth = LocationUtils.getAzimuth(39.9042, 116.4074, 31.2304, 121.4737); System.out.println("距离:" + distance + " 米"); System.out.println("方向角:" + azimuth + " 度"); ``` 其中,纬度和经度的单位都是度。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值