MySQL 大圆距离（哈弗正弦公式）

php mysql great-circle

2022-08-30 06:26:41

我有一个工作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？我浏览了一下互联网，但大多数关于它的文献都非常混乱。

答案 1

来自 Google Code FAQ - 使用 PHP、MySQL 和 Google Maps 创建商店定位器：

下面是 SQL 语句，它将查找半径为 37，-122 坐标的半径范围内最接近的 20 个位置。它根据该行的纬度/经度和目标纬度/经度计算距离，然后仅要求输入距离值小于 25 的行，按距离对整个查询进行排序，并将其限制为 20 个结果。要按公里而不是英里进行搜索，请将 3959 替换为 6371。

SELECT id, ( 3959 * acos( cos( radians(37) ) * cos( radians( lat ) ) 
* cos( radians( lng ) - radians(-122) ) + sin( radians(37) ) * sin(radians(lat)) ) ) AS distance 
FROM markers 
HAVING distance < 25 
ORDER BY distance 
LIMIT 0 , 20;

答案 2

$greatCircleDistance = acos( cos($latitude0) * cos($latitude1) * cos($longitude0 - $longitude1) + sin($latitude0) * sin($latitude1));

纬度和经度以弧度为单位。

所以

SELECT 
  acos( 
      cos(radians( $latitude0 ))
    * cos(radians( $latitude1 ))
    * cos(radians( $longitude0 ) - radians( $longitude1 ))
    + sin(radians( $latitude0 )) 
    * sin(radians( $latitude1 ))
  ) AS greatCircleDistance 
 FROM yourTable;

是您的 SQL 查询

要获得以公里或英里为单位的结果，请将结果乘以地球的平均半径（英里，公里或海里）395963713440

您在示例中计算的是一个边界框。如果将坐标数据放在启用了空间的 MySQL 列中，则可以使用 MySQL 的内置功能来查询数据。

SELECT 
  id
FROM spatialEnabledTable
WHERE 
  MBRWithin(ogc_point, GeomFromText('Polygon((0 0,0 3,3 3,3 0,0 0))'))