Python - 返回给定纬度和经度的美国州名

Question

有谁知道Python中的离线方法可以在给定经纬度对的情况下返回美国州？我需要运行超过 700k 个纬度/经度对，因此 Google 的 API 不适合。

纬度/经度值存储在数据框中，但此处可以视为列表。

示例:

lat_vals = [32]
lon_vals = [-109]

for lat, lon in zip(lat_vals, lon_vals):
    get_state(lat, lon)

>>> NM

Answer 1

由于您拥有大量数据，您可能必须考虑使用 PostGIS对您的数据进行一些分析。 Python 没有任何好的工具来进行空间分析，即使您尝试将其与 OGR GDAL 一起使用也是如此。

最好的办法是建立一个基本的 PostgreSQL 和 PostGIS 数据库，使用 Quantum GIS将 shapefile 加载到数据库中(数据库 >> Spit >> 将 Shapefile 导入 PostgreSQL)

然后您可以运行这些基本的 SQL 命令来提取数据。

SELECT states.name, points.name
FROM states, points
WHERE ST_Contains(states.geom, points.geom)

要获取美国各州 shapefile 数据，请转至 Tiger Census Data页。

如果数据全部加载完毕后想使用Python，可以使用psycopg2连接到您的数据库。

我希望这能有所帮助，而不是一个简单的解决方案。