如果我的回答不够清晰，请在评论中告诉我或者需要更多的信息。也许已经有PHP中解决我所需的问题的方法了。
我正在寻找一个函数，可以从经度或纬度值中加减距离。
原因：我有一个包含所有纬度和经度的数据库，并想形成一个查询来提取X公里（或英里）内的所有城市。我的查询可能会像这样：
Select * From Cities Where (Longitude > X1 and Longitude < X2) And (Latitude > Y1 and Latitude < Y2)
Where X1 = Longitude - (distance)
Where X2 = Longitude + (distance)
Where Y1 = Latitude - (distance)
Where Y2 = Latitude + (distance)
我正在使用PHP和MySql数据库。
也欢迎任何建议! :)
这是一个MySQL查询语句，可以完全满足您的要求。请注意，像这样的东西通常是近似值，因为地球不是完美的球体，也没有考虑到山脉、丘陵、山谷等因素。我们在AcademicHomes.com上使用PHP和MySQL代码，返回距离$latitude、$longitude $radius英里范围内的记录。
$res = mysql_query("SELECT
*
FROM
your_table
WHERE
(
(69.1 * (latitude - " . $latitude . ")) *
(69.1 * (latitude - " . $latitude . "))
) + (
(69.1 * (longitude - " . $longitude . ") * COS(" . $latitude . " / 57.3)) *
(69.1 * (longitude - " . $longitude . ") * COS(" . $latitude . " / 57.3))
) < " . pow($radius, 2) . "
ORDER BY
(
(69.1 * (latitude - " . $latitude . ")) *
(69.1 * (latitude - " . $latitude . "))
) + (
(69.1 * (longitude - " . $longitude . ") * COS(" . $latitude . " / 57.3)) *
(69.1 * (longitude - " . $longitude . ") * COS(" . $latitude . " / 57.3))
) ASC");
的一度宽度的公式：
老实说，这个问题背后的复杂性是如此之大，以至于您最好使用Google Maps等服务来获取数据。 具体而言，地球不是一个完美的球体，两度之间的距离会随着您离赤道越近/越远而变化。
请参见http://en.wikipedia.org/wiki/Geographic_coordinate_system，了解我所说的内容，并查看Google Maps API。
我尝试使用上面的代码，但当点之间的距离在20-30英里范围内时，答案误差太大了，但我可以接受几英里的误差。与我的一个地图伙伴交谈后，我们想出了这个解决方案。代码是使用Python编写的，但您可以很容易地将其翻译成其他语言。为了避免不断转换为弧度，我重做了我的数据库，将纬度/经度点从度转换为弧度。好处在于大部分计算只需要进行一次。
ra = 3963.1906 # radius @ equator in miles, change to km if you want distance in km
rb = 3949.90275 # radius @ poles in miles, change to km if you want distance in km
ra2 = ra * ra
rb2 = rb * rb
phi = self.lat
big_ol_constant = (math.pow(ra2*math.cos(phi), 2) + pow(rb2*math.sin(phi), 2))/ (pow(ra*math.cos(phi), 2) + pow(rb*math.sin(phi), 2))
sqlWhere = "%(distance)g > sqrt((power(lat - %(lat)g,2) + power(lng-%(lng)g,2)) * %(big_ol_constant)g)" % {
'big_ol_constant': big_ol_constant, 'lat': self.lat, 'lng': self.lng, 'distance': distance}
# This is the Django portion of it, where the ORM kicks in. sqlWhere is what you would put after the WHERE part of your SQL Query.
qs = ZipData.objects.extra(where=[sqlWhere]);
当距离较小时，似乎非常准确，而在距离增加到200英里时（当然到那时，你会遇到“直线距离”与“铺设道路”的问题）。
这是我上面提到的ZipData模型。
class ZipData(models.Model):
zipcode = ZipCodeField(null=False, blank=False, verbose_name="ZipCode", primary_key=True)
city = models.CharField(max_length=32, null=False, blank=False)
state = models.CharField(max_length=2)
lat = models.FloatField(null=False, blank=False)
lng = models.FloatField(null=False, blank=False)
根据您包含的城市数量，您可以预先计算列表。我们在这里为内部应用程序执行此操作，其中+100m的不准确度对于我们的设置来说太多了。它通过具有位置1、位置2、距离的两个关键表来工作。然后，我们可以非常快速地从位置1返回距离x位置。
由于计算可以离线完成，因此不会影响系统的运行。用户还可以获得更快的结果。
使用以下URL中的设置，我构建了下面的查询。（请注意，我正在使用codeIgnitor查询数据库）
http://howto-use-mysql-spatial-ext.blogspot.com/2007/11/using-circular-area-selection.html
function getRadius($point="POINT(-29.8368 30.9096)", $radius=2)
{
$km = 0.009;
$center = "GeomFromText('$point')";
$radius = $radius*$km;
$bbox = "CONCAT('POLYGON((',
X($center) - $radius, ' ', Y($center) - $radius, ',',
X($center) + $radius, ' ', Y($center) - $radius, ',',
X($center) + $radius, ' ', Y($center) + $radius, ',',
X($center) - $radius, ' ', Y($center) + $radius, ',',
X($center) - $radius, ' ', Y($center) - $radius, '
))')";
$query = $this->db->query("
SELECT id, AsText(latLng) AS latLng, (SQRT(POW( ABS( X(latLng) - X({$center})), 2) + POW( ABS(Y(latLng) - Y({$center})), 2 )))/0.009 AS distance
FROM crime_listing
WHERE Intersects( latLng, GeomFromText($bbox) )
AND SQRT(POW( ABS( X(latLng) - X({$center})), 2) + POW( ABS(Y(latLng) - Y({$center})), 2 )) < $radius
ORDER BY distance
");
if($query->num_rows()>0){
return($query->result());
}else{
return false;
}
}
ALTER FUNCTION [dbo].[DistanceBetween] (@Lat1 as real,
@Long1 as real, @Lat2 as real, @Long2 as real)
RETURNS real
AS
BEGIN
DECLARE @dLat1InRad as float(53);
SET @dLat1InRad = @Lat1 * (PI()/180.0);
DECLARE @dLong1InRad as float(53);
SET @dLong1InRad = @Long1 * (PI()/180.0);
DECLARE @dLat2InRad as float(53);
SET @dLat2InRad = @Lat2 * (PI()/180.0);
DECLARE @dLong2InRad as float(53);
SET @dLong2InRad = @Long2 * (PI()/180.0);
DECLARE @dLongitude as float(53);
SET @dLongitude = @dLong2InRad - @dLong1InRad;
DECLARE @dLatitude as float(53);
SET @dLatitude = @dLat2InRad - @dLat1InRad;
/* Intermediate result a. */
DECLARE @a as float(53);
SET @a = SQUARE (SIN (@dLatitude / 2.0)) + COS (@dLat1InRad)
* COS (@dLat2InRad)
* SQUARE(SIN (@dLongitude / 2.0));
/* Intermediate result c (great circle distance in Radians). */
DECLARE @c as real;
SET @c = 2.0 * ATN2 (SQRT (@a), SQRT (1.0 - @a));
DECLARE @kEarthRadius as real;
/* SET kEarthRadius = 3956.0 miles */
SET @kEarthRadius = 6376.5; /* kms */
DECLARE @dDistance as real;
SET @dDistance = @kEarthRadius * @c;
return (@dDistance);
END我猜这可能会有帮助。
你可以使用勾股定理来计算两对纬度/经度点之间的距离。
如果你有两个位置（Alpha和Beta），你可以使用以下公式来计算它们之间的距离：
SQRT( POW(Alpha_lat - Beta_lat,2) + POW(Alpha_lon - Beta_lon,2) )
有很多（糟糕的选择）
使用数学公式计算距离（将X1-X2和Y1-Y2视为向量）。
预先创建一个查找表，包含所有组合并保存距离。
考虑使用MySQL的GIS特定扩展。这里是我找到的一篇文章。