python - 在 Python 中使用 geoJSON 指向多边形

Question

我有一个 geoJSON 数据库，其中包含很多多边形(特别是人口普查区)，并且我有很多长纬度点。

我希望存在一个高效的 Python 代码来识别给定坐标所在的人口普查区，但是到目前为止，我的谷歌搜索还没有发现任何信息。

谢谢!

Answer 1

我发现了一个有趣的 article描述如何做你想做的事情。

TL;DR:使用 Shapely

你会在文末找到这段代码:

import json
from shapely.geometry import shape, Point
# depending on your version, use: from shapely.geometry import shape, Point

# load GeoJSON file containing sectors
with open('sectors.json') as f:
    js = json.load(f)

# construct point based on lon/lat returned by geocoder
point = Point(-122.7924463, 45.4519896)

# check each polygon to see if it contains the point
for feature in js['features']:
    polygon = shape(feature['geometry'])
    if polygon.contains(point):
        print 'Found containing polygon:', feature