Draw OSM geometries over each other projected to 1km units

draw-projected

#!/usr/bin/env python
#
# Draw OSM geometries over each other projected to 1km units
#
# Install requirements:
#     pip install osmnx matplotlib
#
# Usage:
#     draw-projected "Cyprus island" "RU-MOW"
#
# Copyright (c) 2024 Andrew Grigorev <[email protected]>
#
# This software is licensed under the MIT License - see the MIT License for details:
# https://opensource.org/licenses/MIT

import sys

import matplotlib.pyplot as plt
import osmnx as ox


def norm(gdf):
    gdf = gdf.to_crs(gdf.estimate_utm_crs())
    centroid = gdf.geometry.unary_union.centroid
    gdf = gdf.translate(-centroid.x, -centroid.y)
    return gdf.scale(.001, .001, origin=(0, 0))


ax = plt.axes()
for i in sys.argv[1:]:
    gdf = ox.geocode_to_gdf(i)
    norm(gdf).plot(ax=ax, facecolor="none")
plt.show()