Это пример, чтобы показать, как отобразить точки с долготой и широтой в окрестности.Первый шаг - найти файл геойсона с координатами границ окрестностей.
Считайте файл геойсона:
`
import json
from shapely.geometry import shape, Point
with open('new_dcneighorhoodboundaries.json') as f:
js = json.load(f)
`
Используйте следующую функцию для итерации по каждой строке:
def find_neighborhood(js,point):
# check each polygon to see if it contains the point
i=0
for feature in js['features']:
polygon = shape(feature['geometry'])
if polygon.contains(point):
return(feature['properties']['id'])
Используйте функцию внутри цикла for для определения окрестности каждой записи / строки:
X=[]
for i in range(0,len(df)):
pt=Point(df.iloc[i][21],df.iloc[i][20])
a=find_neighborhood(js=js,point=pt)
X.append(a)
Имейте в виду, что вы можетенеобходимо внести изменения в приведенный выше код в зависимости от структуры вашего файла Geojson.Вы можете использовать следующий код, чтобы заглянуть в файл геоджона:
import json
from pprint import pprint
with open('new_dcneighorhoodboundaries.json') as f:
data = json.load(f)
pprint(data)
PS: хотя описанный выше метод предназначен для окрестности, его можно экстраполировать на любойформа (город / штат и т. д.)