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加载到数据库中（数据库>>吐>>将Shapefile导入PostgreSQL）

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

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

要获取USA States shapefile数据，请转到“ 老虎普查数据”页面。

如果要在全部加载数据后使用Python，则可以使用psycopg2连接到数据库。

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