diff --git a/sgeop/nodes.py b/sgeop/nodes.py index 608766c..fde8a07 100644 --- a/sgeop/nodes.py +++ b/sgeop/nodes.py @@ -229,7 +229,7 @@ def _makes_loop_contact( loops, not_loops = _loops_and_non_loops(edges) loop_points = shapely.points(loops.get_coordinates().values) - loop_gdf = gpd.GeoDataFrame(geometry=loop_points) + loop_gdf = gpd.GeoDataFrame(geometry=loop_points, crs=edges.crs) loop_point_geoms = loop_gdf.geometry # loop points intersecting non-loops