update(resample_to_grid): filter raster nan values from scipy resampl…

…ing routines (#2285)

Address #2273 (comment)

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),