使用经纬度使用 MySQL 查找最近的 10 个城市？答案

【问题标题】：Find closest 10 cities with MySQL using latitude and longitude?使用经纬度使用 MySQL 查找最近的 10 个城市？
【发布时间】：2012-11-21 21:30:36
【问题描述】：

我已经在网上搜索了很多以找到解决我的问题的方法，但我仍然无法弄清楚。

我有一个非常简单的数据库，其中包含 id、城市名称、纬度、经度和 city_info。当有人进入一个城市页面时，想显示附近的 10 个城市。

如何用 MySQL 计算并用 PHP 返回？

我在这个网站上看到了很多建议，但是这些都不起作用。

我尝试了什么但没有成功。我没有得到任何结果。

<?php
$slatitude = 43.2141341;
$slongitude = 64.4368684;
$miles = 200;

//connect

$query = "SELECT *, 
( 3959 * acos( cos( radians('$slatitude') ) * 
cos( radians( latitude ) ) * 
cos( radians( longitude ) - 
radians('$slongitude') ) + 
sin( radians('$slatitude') ) * 
sin( radians( latitude ) ) ) ) 
AS distance FROM cities HAVING distance < '$miles' ORDER BY distance ASC LIMIT 0, 10";

$query = mysql_query($query);
$numrows = mysql_num_rows($query);
if ($numrows > 0){

while ($row = mysql_fetch_assoc($query)){
$id = $row['id'];
$cityname = $row['cityname'];
$latitude = $row['latitude'];
$longitude = $row['longitude'];

echo "$cityname<br />";

}}

?>`

【问题讨论】：

您的数据库中*肯定*有 200 英里范围内的城市吗？
如果没有HAVING 子句，你会得到结果吗？您收到错误消息吗？
我的数据库中肯定有 200 英里范围内的城市。我根本没有收到任何错误。
对于短距离（200 英里以下很短），我建议使用“等角投影”而不是“余弦定律”。

标签： php mysql geocoding latitude-longitude

【解决方案1】：

你不能使用有，因为你没有按任何东西分组。您需要做的是在 where 中重复您在 select 中所做的操作。

$query = "SELECT *, 
( 3959 * acos( cos( radians('$slatitude') ) * 
cos( radians( latitude ) ) * 
cos( radians( longitude ) - 
radians('$slongitude') ) + 
sin( radians('$slatitude') ) * 
sin( radians( latitude ) ) ) ) 
AS distance FROM cities WHERE ( 3959 * acos( cos( radians('$slatitude') ) * 
cos( radians( latitude ) ) * 
cos( radians( longitude ) - 
radians('$slongitude') ) + 
sin( radians('$slatitude') ) * 
sin( radians( latitude ) ) ) ) < '$miles' ORDER BY distance ASC LIMIT 0, 10";

或者你可以这样做：

$query = "
SELECT * FROM (
  select *, 
  ( 3959 * acos( cos( radians('$slatitude') ) * 
  cos( radians( latitude ) ) * 
  cos( radians( longitude ) - 
  radians('$slongitude') ) + 
  sin( radians('$slatitude') ) * 
  sin( radians( latitude ) ) )) as distance from cities
) WHERE distance < '$miles' ORDER BY distance ASC LIMIT 0, 10";

【讨论】：

您需要将第二个查询中的AS distance 移动到子查询内部。
我相信查询会按照最初使用HAVING 子句编写的方式工作。别名表达式可以在 MySQL 中不带 GROUP BY 的 HAVING 子句中使用，但它不是标准 SQL。
感谢您的快速回复。 @Ian Overton，好的，我用您的代码替换了我的代码。我仍然没有得到任何结果...
你能用 SELECT * FROM ( select *, ( 3959 * acos( cos( radians('$slatitude') ) * cos( radians( latitude ) ) * cos( radians(经度）-弧度（'$slongitude'））+ sin（弧度（'$slatitude'））* sin（弧度（纬度））））作为与城市的距离）？
感觉你的sql有错误。您是否能够直接通过数据库回显生成的 sql 并手动运行 sql？请尝试我的两个查询，而不需要 where、order by 和 limit。不过它们应该可以工作。

【解决方案2】：

为了提高查询速度，您可以首先限制您将使用以下子选择之类的计算结果集（请注意，我还使用了不同的方法来计算距离 - 它适用于我，ymmv）。在我的测试中，使用公差为 1 的 where 子句比不使用它的查询快 100 倍以上。

...
$tol = 1; // limits the search to lat/long within 1 from the given values
$query_args = array($lat,$long,$lat-$tol,$lat+$tol,$long-$tol,$long+$tol);
$query = "
  SELECT *,latitude, longitude, 
    SQRT( POW( 69.1 * ( latitude - %s) , 2 ) 
        + POW( 69.1 * ( %s - longitude ) 
          * COS( latitude / 57.3 ) , 2 ) ) 
    AS distance FROM zipcodes
      WHERE latitude > %d AND latitude < %d 
      AND longitude > %d AND longitude < %d 
  ORDER BY distance ASC limit 10
";
...

【讨论】：

【解决方案3】：

这是我用来从给定的纬度/经度集中获取最近的邮政编码的查询。表格很简单：

id、zip、lat、lng

此查询假设您的中心点使用“$slatitude: and :$slongitude”。这将返回距离变量“$miles”（即 3、2.4、1000）内的所有匹配项。在我的脚本中，结果被放入一个数组中，并根据最接近的 X 个结果进行分类。

$getzip = mysql_query("SELECT zip,((ACOS(SIN(".$slatitude." * PI() / 180) * SIN(lat * PI() / 180) + COS(".$slatitude." * PI() / 180) * COS(lat * PI() / 180) * COS((".$slongitude." - lng) * PI() / 180)) * 180 / PI()) * 60 * 1.1515) AS distance FROM zips HAVING distance<='$miles' ORDER BY 'distance' ASC", $link2) or die (mysql_error($link2));

希望对你有帮助！

【讨论】：

【解决方案4】：

最简单的解决方案是：

获取城市的坐标并

select ...
order by abs(`Latt`-reqLatt) * abs(`Long`-reqLong)

它并不完美，但你可以从前 50 名中挑选出你想要的。

【讨论】：

一般来说这可能会起作用，但它不使用投影，因此在半径边缘附近的城市和距离相似的城市排序时会出现一些错误。 “半径”也是“正方形”。
这也没有提供限制 200 英里的方法。

【解决方案5】：

我已经为此使用了 PDO，因为 mysql_ 函数正在被弃用

<?php
 //Connect to database
$dbh = new PDO("mysql:host=$host;dbname=$database", $username, $password);//Change to suit 
$dbh->setAttribute(PDO::ATTR_ERRMODE, PDO::ERRMODE_EXCEPTION);
try {
    // Prepare statement
    $stmt = $dbh->prepare("SELECT  *, ( 3959 * acos( cos( radians(?) ) * cos( radians( lat ) ) * cos( radians( lng ) - radians(?) ) + sin( radians(?) ) * sin( radians( lat ) ) ) ) AS distance FROM cities HAVING distance < ? ORDER BY distance LIMIT 0 , 10");
    // Assign parameters
    $stmt->bindParam(1,$slatitude);
    $stmt->bindParam(2,$slongitude);
    $stmt->bindParam(3,$slatitude);
    $stmt->bindParam(4,$miles);
    //Execute query
    $stmt->setFetchMode(PDO::FETCH_ASSOC);
    $stmt->execute();
    if ($stmt->rowCount()>0) { 
    // Iterate through the rows
    while($row = $stmt->fetch()) {
        $id = $row['id'];
        $cityname = $row['cityname'];
        $latitude = $row['latitude'];
        $longitude = $row['longitude'];
        echo "$cityname<br />";]);
        }
    }
   else{
        echo "No Records";
       }   
}

catch(PDOException $e) {
    echo "I'm sorry I'm afraid you can't do that.". $e->getMessage() ;// Remove or modify after testing 
    file_put_contents('PDOErrors.txt',date('[Y-m-d H:i:s]').", mapSelect.php, ". $e->getMessage()."\r\n", FILE_APPEND);  
 }
//Close the connection
$dbh = null; 
?>

这个DEMO 使用这个查询

【讨论】：