set dataframe kwarg in gridintersect to suppress warnings

Author: dbrakenhoff

.docs/Notebooks/groundwater2023_watershed_example.py

@@ -95,7 +95,7 @@ def densify_geometry(line, step, keep_internal_nodes=True):
 # function to set the active and inactive model area
 def set_idomain(grid, boundary):
     ix = GridIntersect(grid, rtree=True)
-    result = ix.intersect(Polygon(boundary))
+    result = ix.intersect(Polygon(boundary), geo_dataframe=False)
     idx = list(result.cellids)
     idx = np.array(idx, dtype=int)
     nr = idx.shape[0]
@@ -244,7 +244,7 @@ def set_idomain(grid, boundary):
 ixs = flopy.utils.GridIntersect(struct_grid)
 cellids = []
 for sg in sgs:
-    v = ixs.intersect(LineString(sg), sort_by_cellid=True)
+    v = ixs.intersect(LineString(sg), sort_by_cellid=True, geo_dataframe=False)
     cellids += v["cellids"].tolist()
 intersection_sg = np.zeros(struct_grid.shape[1:])
 for loc in cellids:
@@ -318,7 +318,7 @@ def set_idomain(grid, boundary):
 ixs = flopy.utils.GridIntersect(struct_vrc_grid)
 cellids = []
 for sg in sgs:
-    v = ixs.intersect(LineString(sg), sort_by_cellid=True)
+    v = ixs.intersect(LineString(sg), sort_by_cellid=True, geo_dataframe=False)
     cellids += v["cellids"].tolist()
 intersection_sg_vrc = np.zeros(struct_vrc_grid.shape[1:])
 for loc in cellids:
@@ -422,7 +422,7 @@ def set_idomain(grid, boundary):
 ixs = flopy.utils.GridIntersect(struct_gridp)
 cellids = []
 for sg in sgs:
-    v = ixs.intersect(LineString(sg), sort_by_cellid=True)
+    v = ixs.intersect(LineString(sg), sort_by_cellid=True, geo_dataframe=False)
     cellids += v["cellids"].tolist()
 intersection_ngp = np.zeros(struct_gridp.shape[1:])
 for loc in cellids:
@@ -432,7 +432,7 @@ def set_idomain(grid, boundary):
 ixs = flopy.utils.GridIntersect(struct_gridc)
 cellids = []
 for sg in sgs:
-    v = ixs.intersect(LineString(sg), sort_by_cellid=True)
+    v = ixs.intersect(LineString(sg), sort_by_cellid=True, geo_dataframe=False)
     cellids += v["cellids"].tolist()
 intersection_ngc = np.zeros(struct_gridc.shape[1:])
 for loc in cellids:
@@ -518,7 +518,7 @@ def set_idomain(grid, boundary):
 ixs = flopy.utils.GridIntersect(quadtree_grid)
 cellids = []
 for sg in sgs:
-    v = ixs.intersect(LineString(sg), sort_by_cellid=True)
+    v = ixs.intersect(LineString(sg), sort_by_cellid=True, geo_dataframe=False)
     cellids += v["cellids"].tolist()
 intersection_qg = np.zeros(quadtree_grid.shape[1:])
 for loc in cellids:
@@ -590,6 +590,7 @@ def set_idomain(grid, boundary):
         LineString(sg),
         return_all_intersections=True,
         sort_by_cellid=True,
+        geo_dataframe=False,
     )
     cellids += v["cellids"].tolist()
 intersection_tg = np.zeros(triangular_grid.shape[1:])
@@ -639,6 +640,7 @@ def set_idomain(grid, boundary):
         LineString(sg),
         return_all_intersections=True,
         sort_by_cellid=True,
+        geo_dataframe=False,
     )
     cellids += v["cellids"].tolist()
 intersection_vg = np.zeros(voronoi_grid.shape[1:])

.docs/Notebooks/mf6_parallel_model_splitting_example.py

@@ -355,7 +355,7 @@ def string2geom(geostring, conversion=None):
 # +
 # calculate and set idomain
 ix = flopy.utils.GridIntersect(modelgrid, rtree=True)
-result = ix.intersect(Polygon(boundary_polygon))
+result = ix.intersect(Polygon(boundary_polygon), geo_dataframe=False)
 idxs = tuple(zip(*result.cellids))
 idomain = np.zeros((nrow, ncol), dtype=int)
 idomain[idxs] = 1
@@ -372,7 +372,7 @@ def string2geom(geostring, conversion=None):
 ixs = flopy.utils.GridIntersect(modelgrid)
 cellids = []
 for seg in segs:
-    v = ixs.intersect(LineString(seg), sort_by_cellid=True)
+    v = ixs.intersect(LineString(seg), sort_by_cellid=True, geo_dataframe=False)
     cellids += v["cellids"].tolist()
 intersection_rg = np.zeros(modelgrid.shape[1:])
 for loc in cellids:
@@ -405,7 +405,7 @@ def string2geom(geostring, conversion=None):
 drn_cellids = []
 drn_lengths = []
 for seg in segs:
-    v = ixs.intersect(LineString(seg), sort_by_cellid=True)
+    v = ixs.intersect(LineString(seg), sort_by_cellid=True, geo_dataframe=False)
     drn_cellids += v["cellids"].tolist()
     drn_lengths += v["lengths"].tolist()
 # -

autotest/test_model_splitter.py

@@ -1014,6 +1014,7 @@ def string2geom(geostring, conversion=None):
             boundary,
         ],
         shapetype="Polygon",
+        geo_dataframe=False,
     )
     r, c = list(zip(*list(result.cellids)))
     idomain = np.zeros(modelgrid.shape, dtype=int)
@@ -1033,7 +1034,12 @@ def string2geom(geostring, conversion=None):
     lengths = []
     for sg in stream_segs:
         sg = string2geom(sg)
-        v = ixs.intersect(sg, shapetype="LineString", sort_by_cellid=True)
+        v = ixs.intersect(
+            sg,
+            shapetype="LineString",
+            sort_by_cellid=True,
+            geo_dataframe=False,
+        )
         cellids += v["cellids"].tolist()
         lengths += v["lengths"].tolist()