From 59037ef3cf7286eec921d12dfa3bdb0bca6bce41 Mon Sep 17 00:00:00 2001 From: jlarsen Date: Wed, 7 Aug 2024 13:27:53 -0700 Subject: [PATCH] update(resample_to_grid): filter raster nan values from scipy resampling routines --- flopy/utils/rasters.py | 7 +++++++ 1 file changed, 7 insertions(+) diff --git a/flopy/utils/rasters.py b/flopy/utils/rasters.py index fe9155121..4fabf0482 100644 --- a/flopy/utils/rasters.py +++ b/flopy/utils/rasters.py @@ -537,6 +537,13 @@ 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),