osmnx

import geopandas as gpd
import folium
import osmnx as ox
 
 
 
 
def visualize_polygons(geometry):
    lats, lons = get_lat_lon(geometry)
    m = folium.Map(location=[sum(lats) / len(lats), sum(lons) / len(lons)], zoom_start=12, tiles='cartodbpositron')
    overlay = gpd.GeoSeries(geometry).to_json()
    folium.GeoJson(overlay, name='boundary').add_to(m)
 
# Добавим ПЧ5 на карту
    folium.Marker(location=[55.740441, 37.966298],
                  popup='<b>улица Керамическая, д. 1а</b>',
                  tooltip='ПЧ-5',
                  icon=folium.Icon(color="green")).add_to(m)  # подкрасили метку на карте
 
    m.save('Geo.html')
 
 
# выводим центроиды полигонов
def get_lat_lon(geometry):
    lon = geometry.apply(lambda x: x.x if x.type == 'Point' else x.centroid.x)
    lat = geometry.apply(lambda x: x.y if x.type == 'Point' else x.centroid.y)
    return lat, lon
 
 
# выгрузим границы Балашихи из OSM
cities = ['Балашиха']
polygon_krd = ox.geometries_from_place(cities, {'boundary': 'administrative'}).reset_index()
polygon_krd = polygon_krd[(polygon_krd['name'] =='городской округ Балашиха')]#Кучинский лес
 
 
 
# посмотрим что получилось
visualize_polygons(polygon_krd['geometry'])