From 15f1b94a45487e42aa36afdd2abb2b111c2197a6 Mon Sep 17 00:00:00 2001 From: Joshua Larsen Date: Wed, 9 Oct 2024 15:19:50 -0700 Subject: [PATCH] fix(resample_to_grid): fix unintended extrapolation (#2331) * remove nan masking from standard scipy resampling methods * allow extrapolate_edges kwarg to be applied to nearest neighbor resampling --- flopy/utils/rasters.py | 9 +-------- 1 file changed, 1 insertion(+), 8 deletions(-) diff --git a/flopy/utils/rasters.py b/flopy/utils/rasters.py index 4fabf0482..39697f97c 100644 --- a/flopy/utils/rasters.py +++ b/flopy/utils/rasters.py @@ -537,13 +537,6 @@ def resample_to_grid( arr = self.get_array(band, masked=True) arr = arr.flatten() - # filter out nan values from the original dataset - if np.isnan(np.sum(arr)): - idx = np.isfinite(arr) - rxc = rxc[idx] - ryc = ryc[idx] - arr = arr[idx] - # step 3: use griddata interpolation to snap to grid data = griddata( (rxc, ryc), @@ -574,7 +567,7 @@ def resample_to_grid( else: raise TypeError(f"{method} method not supported") - if extrapolate_edges and method != "nearest": + if extrapolate_edges: xc = modelgrid.xcellcenters yc = modelgrid.ycellcenters