更新：我在http://www.mullie.eu/geographic-searches/上写了一篇关于这个特定主题的更详细的博文
--
使用 Google Maps API 遍历所有可用城镇以获取其纬度和经度。将这些保存在某处（数据库）。 - 请注意，Google 不会接听大量电话，因此请限制您的电话。
那么，在获取城镇的时候，可以使用类似下面的代码来获取一定范围内的城市：
public static function getNearby($lat, $lng, $type = 'cities', $limit = 50, $distance = 50, $unit = 'km')
{
// radius of earth; @note: the earth is not perfectly spherical, but this is considered the 'mean radius'
if ($unit == 'km') $radius = 6371.009; // in kilometers
elseif ($unit == 'mi') $radius = 3958.761; // in miles
// latitude boundaries
$maxLat = (float) $lat + rad2deg($distance / $radius);
$minLat = (float) $lat - rad2deg($distance / $radius);
// longitude boundaries (longitude gets smaller when latitude increases)
$maxLng = (float) $lng + rad2deg($distance / $radius / cos(deg2rad((float) $lat)));
$minLng = (float) $lng - rad2deg($distance / $radius / cos(deg2rad((float) $lat)));
// get results ordered by distance (approx)
$nearby = (array) FrontendDB::getDB()->retrieve('SELECT *
FROM table
WHERE lat > ? AND lat < ? AND lng > ? AND lng < ?
ORDER BY ABS(lat - ?) + ABS(lng - ?) ASC
LIMIT ?;',
array($minLat, $maxLat, $minLng, $maxLng, (float) $lat, (float) $lng, (int) $limit));
return $nearby;
}
以上代码注意事项：
- 使用了自己的数据库包装器，因此转换为 mysql_query、PDO、...
- 这并不准确。我们无法在 DB 中进行精确的球面计算，因此我们采用了经纬度上下限。这基本上意味着一个位置比您的距离稍远（例如，在远东北，就在实际半径之外（实际上几乎是一个圆），但仍在最大纬度和经度之内（因为我们比较它在数据库中的平方范围内）。这只会在您的半径范围内提供粗略但绝对 100% 准确的城市选择。
我将尝试说明这一点：
_________________
| / \ |
| Y / \ |
| / \ |
|( X )|
| \ / |
| \ / |
|______\_/______|
上面的圆圈（有点）是你想要在其中找到位置的实际半径，基于位置 X。这很难直接从你的数据库中完成，所以我们实际上从数据库中获取的是周围的正方形.如您所见，位置（如 Y）可能位于这些最大和最小边界内，尽管它们实际上并不在请求的半径范围内。不过这些可以稍后通过 PHP 过滤掉。
要解决最后一个问题，您可以循环所有结果并计算您的根位置和找到的接近匹配项之间的确切距离，以计算它们是否真的在您的半径范围内。为此，您可以使用以下代码：
public static function getDistance($lat1, $lng1, $lat2, $lng2, $unit = 'km')
{
// radius of earth; @note: the earth is not perfectly spherical, but this is considered the 'mean radius'
if ($unit == 'km') $radius = 6371.009; // in kilometers
elseif ($unit == 'mi') $radius = 3958.761; // in miles
// convert degrees to radians
$lat1 = deg2rad((float) $lat1);
$lng1 = deg2rad((float) $lng1);
$lat2 = deg2rad((float) $lat2);
$lng2 = deg2rad((float) $lng2);
// great circle distance formula
return $radius * acos(sin($lat1) * sin($lat2) + cos($lat1) * cos($lat2) * cos($lng1 - $lng2));
}
这将计算位置 X 和位置 Y 之间的（准）精确距离，然后您可以准确地过滤掉那些足够接近以通过粗略的 db-fetch 的城市，但不仅仅是靠近到实际在您的范围内界限。