python - 将纬度/经度分配给邮政编码

Question

我在 R 中有一个数据框纬度和经度点。使用 R data science toolkit ，我可以将这些点分配给国家/州/县/市/选区/邻里，这很有用。

我想将这些点分配给 R 或 Python 中的 5 位邮政编码(甚至 9 位邮政编码!)。有没有简单的方法可以做到这一点？

编辑:我找到了一个 file that contains zip code boundaries .不幸的是，它是一个 .lpk 文件。如何将其导入 R？
编辑 2:我找到了一个 shape file ，这将更易于使用。

Answer 1

1. 在某处找到一个邮政编码 shapefile(.shp 格式，有时称为 ArcGIS 形状文件)。

2. 使用 maptools 包的 readShapePoly 将其加载到 R 中 命令。

3. 创建一个 SpatialPointsDataFrame 来保存点。

4. 确保您的预测是正确的。

5. 使用 sp 包的 overlay 命令将点叠加到多边形中。

您可能会发现 taRifx 包中的 cleanLatLon 很有用。