给定坐标，如何获取 10 英里半径内的所有邮政编码？

Question

我有一个位置（纬度和经度）。如何获取部分或全部位于我所在位置 10 英里半径范围内的邮政编码列表？

该解决方案可以是调用众所周知的 Web 服务（google 地图、bing 地图等）或本地数据库解决方案（客户端有 sql server 2005）或算法。

我看到了一些类似的问题，但其中的所有答案几乎都与使用 SQL Server 2008 地理功能有关，而我无法使用该功能。

Answer 1

从包含邮政编码及其相应纬度和经度坐标的邮政编码数据库开始：

http://www.zipcodedownload.com/Products/Product/Z5Commercial/Standard/Overview/

要获取纬度和经度之间的距离，您需要一个好的距离公式。该网站有几个变体：

http: //www.meridianworlddata.com/distance-calculation/

“大圆距离”公式有点极端。根据我的经验，这个效果很好：

sqrt(x * x + y * y)

where x = 69.1 * (lat2 - lat1)
and y = 69.1 * (lon2 - lon1) * cos(lat1/57.3)

然后你的 SQL 查询将如下所示：

select zd.ZipCode
from ZipData zd
where 
    sqrt(
        square(69.1 * (zd.Latitude - @Latitude)) +
        square(69.1 * (zd.Longitude - @Longitude) * cos(@Latitude/57.3))
    ) < @Distance

祝你好运！

Answer 2

首先，您需要一个包含所有邮政编码及其相应纬度和经度的数据库。在澳大利亚，只有几千个（而且信息很容易获得），但我认为这在美国可能是一项更困难的任务。

其次，如果您知道自己所在的位置，并且知道要查找的半径，则可以查找该半径内的所有邮政编码。用 PHP 编写的简单内容如下：（抱歉，它不是用 C# 编写的）

function distanceFromTo($latitude1,$longitude1,$latitude2,$longitude2,$km){
  $latitude1  = deg2rad($latitude1);
  $longitude1 = deg2rad($longitude1);
  $latitude2  = deg2rad($latitude2);
  $longitude2 = deg2rad($longitude2);
  $delta_latitude  = $latitude2  - $latitude1;
  $delta_longitude = $longitude2 - $longitude1;
  $temp = pow(sin($delta_latitude/2.0),2) + cos($latitude1) * cos($latitude2) * pow(sin($delta_longitude/2.0),2);
  $earth_radius = 3956;
  $distance = $earth_radius * 2 * atan2(sqrt($temp),sqrt(1-$temp));
  if ($km)
    $distance = $distance * 1.609344;
  return $distance;
}

Answer 3

大多数搜索都使用质心。为了使用 10 英里内的部分邮政编码，您必须购买邮政编码多边形数据库 (*)。然后实施一种算法，检查顶点在 10 英里半径内的邮政编码。为了正确完成，您应该使用半正矢公式进行距离测量。通过一些巧妙的数据结构，您可以显着减少搜索空间。同样，通过存储并最初与 zipcoe 范围（北、西、东、南）进行比较，可以大大加快搜索速度。

(*) 注意：从技术上讲，邮政编码不是多边形！我知道我们都这样认为它们，但实际上它们是数据点（街道地址）的集合，这就是 USPS 真正使用它们的方式。这意味着邮政编码可以包含其他邮政编码；邮政编码可以由多个“多边形”组成；并且邮政编码可以与其他邮政编码重叠。大多数情况应该不是问题，但您必须处理可以定义为多个多边形的邮政编码。