选择半径内的纬度和经度点

Question

我正在使用论坛在给定位置和半径的数据库中选择纬度/经度点，但是我无法正常工作。只有当我使用 6371 公里（世界的大小）的半径时，我才能让查询工作，我认为论坛存在问题，但我似乎无法弄清楚在哪里......

这是我的 sql 查询

   $latitude = $_POST["latitude"]; // latitude of centre of bounding circle in degrees
   $longitude = $_POST["longitude"];; // longitude of centre of bounding circle in degrees
   $rad = $_POST["radius"]; // radius of bounding circle in kilometers
   $R = 6371;// earth's mean radius, km
   $maxLat = $latitude + rad2deg($rad/$R);
   $minLat = $latitude - rad2deg($rad/$R);
    // compensate for degrees longitude getting smaller with increasing latitude
   $maxLon = $longitude + rad2deg($rad/$R/cos(deg2rad($latitude)));
   $minLon = $longitude - rad2deg($rad/$R/cos(deg2rad($latitude)));



   $con = mysqli_connect("localhost","root","Winter#05");

   if (!$con)
   {
   die('Could not connect: ' . mysql_error());
   }
   mysqli_select_db($con, "appdatabase");


   $result = mysqli_query($con,"select id, PostTitle, SubmitDate, PostVotes, ImagePath, comments, latitude, longitude, acos(sin($latitude)*sin(radians(latitude)) + cos($latitude)*cos(radians(latitude))*cos(radians(longitude)-($longitude))) * $R as D 
   from ( select id, PostTitle, SubmitDate, PostVotes, comments, ImagePath, latitude, longitude 
   from posts where latitude > $minLat and latitude < $maxLat and longitude > $minLon and longitude < $maxLon ) as first_cut 
   where acos(sin($latitude)*sin(radians(latitude)) + cos($latitude)*cos(radians(latitude))*cos(radians(longitude) - ($longitude))) * $R < $rad order by D") or die('Errant query:');

   while($row = mysqli_fetch_assoc($result))
   {
        $output[]=$row;
   }


   print(json_encode($output));

   mysqli_close($con);

Answer 1

一个问题是$latitude 和$longitude 似乎正在以度 为单位合并到SQL 文本中。该公式（在 SQL 查询中实现）看起来期望这些位置的值以 弧度 为单位。

---

如果我必须阅读/维护这段代码，如果它是这样编写的，那么理解查询会更容易：

SELECT p.id
     , p.PostTitle
     , p.SubmitDate
     , p.PostVotes
     , p.ImagePath
     , p.comments
     , p.latitude
     , p.longitude
     , ACOS( SIN(       (  $latitude ))
           * SIN(RADIANS( p.latitude ))
           + COS(       (  $latitude ))
           * COS(RADIANS( p.latitude ))
           * COS(RADIANS( p.longitude ) -        ( $longitude ))
           ) * $R
       AS D
  FROM posts p
 WHERE p.latitude  > $minLat
   AND p.latitude  < $maxLat 
   AND p.longitude > $minLon
   AND p.longitude < $maxLon 
HAVING D < $rad
 ORDER BY D

大圆距离公式看起来不错，但我想知道为什么 $latitude 和 $longitude 以弧度而不是度为单位提供。

然后我查看 PHP 代码，你瞧，它们实际上是以度为单位的，我发现缺少到弧度的转换。

返回距离“D”的表达式中缺少的只是$latitude 和$longitude 周围的RADIANS 函数

所以这个：

     , ACOS( SIN(       (  $latitude ))
           * SIN(RADIANS( p.latitude ))
           + COS(       (  $latitude ))
           * COS(RADIANS( p.latitude ))
           * COS(RADIANS( p.longitude ) -        ( $longitude ))
           ) * $R
       AS D

被替换为

     , ACOS( SIN(RADIANS(  $latitude ))
           * SIN(RADIANS( p.latitude ))
           + COS(RADIANS(  $latitude ))
           * COS(RADIANS( p.latitude ))
           * COS(RADIANS( p.longitude ) - RADIANS( $longitude ))
           ) * $R
       AS D

如果我要维护它，我也会将它转换为带有 bind placeholders 的 prepared statement。至少，我会正确地转义任何被合并到 SQL 文本中的潜在不安全值...mysqli_real_escape_string

例如

$sql .= "  + COS(RADIANS( " . mysqli_real_escape_string($con,$latitude) . " ))"

---

我怀疑$minLat、$maxLat、$minLon、$maxLon 的计算。我一定会测试这些。为了测试没有这些的查询，我可以注释掉整个 WHERE 子句，然后计算表中 每个 翻转行的距离 D，然后通过 @987654338 过滤它们@ 子句。

我们确实希望在 WHERE 子句中包含谓词（如果正确完成）将提供一个“边界框”，它将限制我们需要通过大圆计算进行曲折的行数。

消除内联视图应该会提高性能，因为 MySQL 不会有具体化派生表的开销。

Answer 2

我目前正在使用以下工具，并且对（相对）准确性感到满意。

( 
    3958*3.1415926 * sqrt( 
        ( table.latitude - $latitude ) * ( table.latitude - $latitude ) 
        + cos( table.latitude / 57.29578 ) 
        * cos( $latitude / 57.29578 ) 
        * ( table.longitude - $longitude ) 
        * ( table.longitude - $longitude ) 
    ) / 180
) <= $radius