I cannot find any real answer so I'd like to ask how to crop shapefile with roads in Python.
The goal is to leave roads (nodes) in black circle only, then output that to dataframe in geopandas (possibly save later as shapefile).
Assuming, that I have some place's coordinates in lat + lon and a very big file of roads, how to clip this file by radius, let's say, of 1km?
Theoretically I can create shapefile of the circle basing on coordinates and then crop roads by this shapefile using some external library like ogr2ogr, but it seems to be very slow and I don't want to use subprocess.
Is it possible to do that directly in geopandas? (I can load all roads to pd dataframe first)
Is there any already implemented option to crop roads but leave all these nodes that are even partially in black circle?