Спасибо, что поделился!
В случае использования для других, которые прокомментировали вопрос об отсутствующем коде на третьем шаге, где загружается производная сеть OSMnx, ниже приведен один подход (с благодарностью коллеге, от которого это было адаптировано), который, я думаю, будет работать, как только у вас будет Сетевая переменная OSMnx 'G' определена:
# преобразовать граф в кадры геоданных с узлами и ребрами
gdf_nodes = ox.graph_to_gdfs(G, узлы=Истина, края=Ложь)
gdf_edges = ox.graph_to_gdfs(G, узлы=Ложь, края=Истина)
# загрузка сети
network=pandana.network.Network(gdf_nodes["x"], gdf_nodes["y"], gdf_edges["u"], gdf_edges["v"],gdf_edges[["длина"]])
# установить граничное расстояние для анализа в метрах; обратите внимание, что если доступ не найден, расстояние доступа будет записано как это значение
отсечка = 500
network.precompute (отсечка + 1)
# теперь вы можете перейти к шагу 4, например
network.set_pois(category='all', maxdist=cutoff, maxitems=1, x_col=destinations['long'], y_col=destinations['lat'])
result = network.nearest_pois(distance=cutoff, категория='все', num_pois=1)
результат[‘доступ’] = результат ‹ отсечка
