update(Raster): added min and max resampling for hydrologic resamplin…

…g and streamflow representation
modflowpy · Nov 10, 2021 · d6da642 · d6da642
1 parent 3e770ec
commit d6da642
Showing 1 changed file with 12 additions and 4 deletions.
diff --git a/flopy/utils/rasters.py b/flopy/utils/rasters.py
@@ -399,8 +399,8 @@ def resample_to_grid(
                 method=method,
             )
 
-        elif method in ("median", "mean"):
-            # these methods are slow and could use a speed u
+        elif method in ("median", "mean", "min", "max"):
+            # these methods are slow and could use a speed up
             ncpl = modelgrid.ncpl
             data_shape = modelgrid.xcellcenters.shape
             if isinstance(ncpl, (list, np.ndarray)):
@@ -455,8 +455,12 @@ def resample_to_grid(
 
                     if method == "median":
                         val = np.nanmedian(rstr_data)
-                    else:
+                    elif method == "mean":
                         val = np.nanmean(rstr_data)
+                    elif method == "max":
+                        val = np.nanmax(rstr_data)
+                    else:
+                        val = np.nanmin(rstr_data)
 
                     data[node] = val
         else:
@@ -533,8 +537,12 @@ def __threaded_resampling(
 
         if method == "median":
             val = np.nanmedian(rstr_data)
-        else:
+        elif method == "mean":
             val = np.nanmean(rstr_data)
+        elif method == "max":
+            val = np.nanmax(rstr_data)
+        else:
+            val = np.nanmin(rstr_data)
         q.put((node, val))
         container.release()