注意:虽然我使用包含荷兰邮政编码的邮政编码数据库,但这个问题与国家/地区无关。
我有一个数据库,其中包含荷兰的每个邮政编码及其 x 和 y 坐标(纬度/经度)。
例如,我有邮政编码:$baseZipCode = 1044;
,坐标如下:
x coordinate = 4,808855
y coordinate = 52,406332
现在,我想从 $baseZipCode
中查找带有 $range
的所有其他邮政编码。
例如:
SELECT
zipcode
FROM
zipcodes
WHERE
????? // Need help here
问题是地球并不完全是圆的。我找到了很多关于从 a 到 b
计算的教程,但这不是我需要的。
有人知道吗?
更新 感谢 Captaintokyo 我发现了这个:
想要查找距另一个邮政编码或点一定英里/公里半径内的所有邮政编码和相应距离吗?这个问题需要经纬度坐标来解决。对地址进行地理编码可提供地址的纬度/经度坐标。
首先,您需要一个包含所有邮政编码及其相应纬度和经度坐标的数据库:
CREATE TABLE `zipcodes` (
`zipcode` varchar(5) NOT NULL DEFAULT '',
`city` varchar(100) NOT NULL DEFAULT '',
`state` char(2) NOT NULL DEFAULT '',
`latitude` varchar(20) NOT NULL DEFAULT '',
`longitude` varchar(20) NOT NULL DEFAULT '',
KEY `zipcode` (`zipcode`),
KEY `state` (`state`)
)
因此,一旦您拥有了数据库,您就希望找到以中心点为中心的特定英里半径内的所有邮政编码。如果中心点是另一个邮政编码,只需在数据库中查询该邮政编码的纬度和经度坐标即可。那么代码如下:
// ITITIAL POINT
$coords = array('latitude' => "32.8", 'longitude' => "-117.17");
//RADIUS
$radius = 30;
// SQL FOR KILOMETERS
$sql = "SELECT zipcode, ( 6371 * acos( cos( radians( {$coords['latitude']} ) ) * cos( radians( latitude ) ) * cos( radians( longitude ) - radians( {$coords['longitude']} ) ) + sin( radians( {$coords['latitude']} ) ) * sin( radians( latitude ) ) ) ) AS distance FROM zipcodes HAVING distance <= {$radius} ORDER BY distance";
// SQL FOR MILES
$sql = "SELECT zipcode, ( 3959 * acos( cos( radians( {$coords['latitude']} ) ) * cos( radians( latitude ) ) * cos( radians( longitude ) - radians( {$coords['longitude']} ) ) + sin( radians( {$coords['latitude']} ) ) * sin( radians( latitude ) ) ) ) AS distance FROM zipcodes HAVING distance <= {$radius} ORDER BY distance";
// OUTPUT THE ZIPCODES AND DISTANCES
$query = mysql_query($sql);
while($row = mysql_fetch_assoc($query)){
echo "{$row['zipcode']} ({$row['distance']})<br>\n";
}
(雅虎和 Google 都提供免费地理编码服务。)
最佳答案
您必须使用名为 Haversine formula 的东西:
$sql = "
SELECT zipcode
FROM zipcodes
WHERE ".mysqlHaversine($lat, $lon, $distance)."
";
公式:
function mysqlHaversine($lat = 0, $lon = 0, $distance = 0)
{
if($distance > 0)
{
return ('
((6372.797 * (2 *
ATAN2(
SQRT(
SIN(('.($lat*1).' * (PI()/180)-latitude*(PI()/180))/2) *
SIN(('.($lat*1).' * (PI()/180)-latitude*(PI()/180))/2) +
COS(latitude * (PI()/180)) *
COS('.($lat*1).' * (PI()/180)) *
SIN(('.($lon*1).' * (PI()/180)-longitude*(PI()/180))/2) *
SIN(('.($lon*1).' * (PI()/180)-longitude*(PI()/180))/2)
),
SQRT(1-(
SIN(('.($lat*1).' * (PI()/180)-latitude*(PI()/180))/2) *
SIN(('.($lat*1).' * (PI()/180)-latitude*(PI()/180))/2) +
COS(latitude * (PI()/180)) *
COS('.($lat*1).' * (PI()/180)) *
SIN(('.($lon*1).' * (PI()/180)-longitude*(PI()/180))/2) *
SIN(('.($lon*1).' * (PI()/180)-longitude*(PI()/180))/2)
))
)
)) <= '.($distance/1000). ')');
}
return '';
}
通常,我不会在不了解代码工作方式的情况下使用代码,但我必须承认这个功能有点超出我的理解......
关于php - mySQL 选择 y 范围内 x 公里/英里内的邮政编码,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/4011944/