MySQL大圆距离(Haversine公式)

Question

MySQL大圆距离(Haversine公式)

我有一个可以工作的PHP脚本，它获取经度和纬度值，然后将它们输入到MySQL查询中。我只想让它成为MySQL。下面是我当前的PHP代码：

if ($distance != "Any" && $customer_zip != "") { //get the great circle distance

    //get the origin zip code info
    $zip_sql = "SELECT * FROM zip_code WHERE zip_code = '$customer_zip'";
    $result = mysql_query($zip_sql);
    $row = mysql_fetch_array($result);
    $origin_lat = $row['lat'];
    $origin_lon = $row['lon'];

    //get the range
    $lat_range = $distance/69.172;
    $lon_range = abs($distance/(cos($details[0]) * 69.172));
    $min_lat = number_format($origin_lat - $lat_range, "4", ".", "");
    $max_lat = number_format($origin_lat + $lat_range, "4", ".", "");
    $min_lon = number_format($origin_lon - $lon_range, "4", ".", "");
    $max_lon = number_format($origin_lon + $lon_range, "4", ".", "");
    $sql .= "lat BETWEEN '$min_lat' AND '$max_lat' AND lon BETWEEN '$min_lon' AND '$max_lon' AND ";
    }

有谁知道如何把这个变成完整的MySQL吗？我浏览过一些互联网，但它上的大部分文献都令人困惑。

Answer 1

从…GoogleCodeFAQ-使用PHP、MySQL和GoogleMaps创建一个商店定位器:这是SQL语句，它将找到距离37，-122坐标半径25英里内最接近的20个位置。它根据该行的纬度/经度和目标纬度/经度计算距离，然后只请求距离值小于25的行，按距离对整个查询进行排序，并将其限制为20个结果。要搜索公里而不是英里，请将3959替换为6371。SELECT&nbsp;id,&nbsp;(&nbsp;3959&nbsp;*&nbsp;acos(&nbsp;cos(&nbsp;radians(37)&nbsp;)&nbsp;*&nbsp;cos(&nbsp;radians(&nbsp;lat&nbsp;)&nbsp;)&nbsp;*&nbsp;cos(&nbsp;radians(&nbsp;lng&nbsp;)&nbsp;-&nbsp;radians(-122)&nbsp;)&nbsp;+&nbsp;sin(&nbsp;radians(37)&nbsp;)&nbsp;*&nbsp; sin(radians(lat))&nbsp;)&nbsp;)&nbsp;AS&nbsp;distance&nbsp; FROM&nbsp;markers&nbsp; HAVING&nbsp;distance&nbsp;<&nbsp;25&nbsp;ORDER&nbsp;BY&nbsp;distance&nbsp; LIMIT&nbsp;0&nbsp;,&nbsp;20;

Answer 2

$greatCircleDistance = acos( cos($latitude0) * cos($latitude1) * cos($longitude0 - $longitude1) + sin($latitude0) * sin($latitude1));纬度和经度为弧度。所以SELECT&nbsp; &nbsp;&nbsp;acos(&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;cos(radians(&nbsp;$latitude0&nbsp;)) &nbsp;&nbsp;&nbsp;&nbsp;*&nbsp;cos(radians(&nbsp;$latitude1&nbsp;)) &nbsp;&nbsp;&nbsp;&nbsp;*&nbsp;cos(radians(&nbsp;$longitude0&nbsp;)&nbsp;-&nbsp;radians(&nbsp;$longitude1&nbsp;)) &nbsp;&nbsp;&nbsp;&nbsp;+&nbsp;sin(radians(&nbsp;$latitude0&nbsp;))&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;*&nbsp;sin(radians(&nbsp;$latitude1&nbsp;)) &nbsp;&nbsp;)&nbsp;AS&nbsp;greatCircleDistance&nbsp; &nbsp;FROM&nbsp;yourTable;是您的sql查询。若要获得以公里或英里为单位的结果，请将结果乘以地球的平均半径(3959迈尔斯6371公里或3440(海里)在你的例子中，你正在计算的是一个边框。如果将坐标数据放在空间启用MySQL列，你可以用MySQL内置功能查询数据。SELECT&nbsp; &nbsp;&nbsp;id FROM&nbsp;spatialEnabledTable WHERE&nbsp; &nbsp;&nbsp;MBRWithin(ogc_point,&nbsp;GeomFromText('Polygon((0&nbsp;0,0&nbsp;3,3&nbsp;3,3&nbsp;0,0&nbsp;0))'))