fix(GridIntersect): combine list of geometries using unary_union (#1923)

modflowpy · Aug 25, 2023 · 22205f4 · 22205f4
1 parent 4e00489
commit 22205f4
Showing 1 changed file with 5 additions and 4 deletions.
diff --git a/flopy/utils/gridintersect.py b/flopy/utils/gridintersect.py
@@ -18,6 +18,10 @@
     if Version(shapely.__version__) < Version("1.8"):
         warnings.warn("GridIntersect requires shapely>=1.8.")
         shapely = None
+    if SHAPELY_GE_20:
+        from shapely import unary_union
+    else:
+        from shapely.ops import unary_union
 else:
     SHAPELY_GE_20 = False
 
@@ -1487,10 +1491,7 @@ def _intersect_linestring_structured(
                     tempverts.append(vertices[i])
                     ishp = ixshapes[i]
                     if isinstance(ishp, list):
-                        if len(ishp) > 1:
-                            ishp = shapely_geo.MultiLineString(ishp)
-                        else:
-                            ishp = ishp[0]
+                        ishp = unary_union(ishp)
                     tempshapes.append(ishp)
             nodelist = tempnodes
             lengths = templengths