使用 Open Street Maps 获取特定邮政编码的街道

Question

我想编写一个将 Countrycode 和 Postcode 作为输入的代码，输出是使用一些使用 GSM 的 api 在给定邮政编码中的街道。

我的策略如下：

1. 我需要获取地区的关系 ID。例如，1991416 是奥地利维也纳第三区的关系 ID。它由 nominatim api 提供：http://nominatim.openstreetmap.org/details.php?place_id=158947085

2. 把id放到这个api url中：http://polygons.openstreetmap.fr/get_wkt.py?id=1991416&params=0

3. 下载多边形后，我可以将收集到的多边形放在立交桥 API 上的此查询中

   ( 方法 （poly：“多边形数据”） ["高速公路"~"^(小学|中学|大学|住宅)$"] [“名称”];

   ); 出几何；

这给了我搜索区的街道。我对这个解决方案的两个问题是 1. 这需要相当长的时间，因为在资源和资源上每个请求询问三个不同的 API 并不容易 2. 我不知道如何从第一步自动收集关系 ID。当我输入像 http://nominatim.openstreetmap.org/search?format=json&country=austria&postalcode=1030 这样的 Nominatim 查询时，我只是得到了该地区的各个点，而不是搜索地区的关系 id 以获得所需的多边形。

所以我的问题是，是否有人可以告诉我如何获取relation_Id 以执行上述工作流程，或者是否有其他更好的方法来解决这个问题。

感谢您的帮助！

最好的问候 丹尼尔

Answer 1

假设您定义了一些相关标签以匹配相关关系，您可以大大简化您的方法，只需调用一次 Overpass API。特别是，您根本不必使用poly，即无需将关系转换为纬度/经度对列表。现在可以使用区域的概念来查询由方式或关系定义的多边形中的某些对象。请查看documentation for more details on areas。

为了获得关系 1991416 的匹配区域，我使用了 postal_code=1030 和 boundary=administrative 作为过滤条件。然后，您可以使用该区域在该特定多边形中搜索方式：

//uncomment the following line, if you need csv output
//[out:csv(::id, ::type, name)];

//adjust area to your needs, filter critera are the same as for relations
area[postal_code=1030][boundary=administrative]->.a;

// Alternative:  {{geocodeArea:name}} -> see
// http://wiki.openstreetmap.org/wiki/Overpass_turbo/Extended_Overpass_Queries

way(area.a)["highway"~"^(primary|secondary|tertiary|residential)$"]["name"];
(._;>;);out meta;

// just for checking if we're looking at the right area
rel(pivot.a);out geom;

在立交桥上试用链接：http://overpass-turbo.eu/s/6uN

注意：并非所有方式/关系都有相应的区域，即适用某些区域生成规则（参见上面的 wiki 页面）。但是，对于您的特定用例，您应该没问题。