fix(Raster): resample_to_grid failure no data masking failure with in…

…t dtype (#1328) * update(test_raster_sampling_methods): update for mean resampling
modflowpy · Jan 21, 2022 · dced106 · dced106
1 parent 16b84e8
commit dced106
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/autotest/t065_test_gridintersect.py b/autotest/t065_test_gridintersect.py
@@ -1398,7 +1398,7 @@ def test_raster_sampling_methods():
     methods = {
         "min": 2088.52343,
         "max": 2103.54882,
-        "mean": 2097.05053,
+        "mean": 2097.05035,
         "median": 2097.36254,
         "nearest": 2097.81079,
         "linear": 2097.81079,

diff --git a/flopy/utils/rasters.py b/flopy/utils/rasters.py
@@ -454,7 +454,7 @@ def resample_to_grid(
             else:
                 for node in range(ncpl):
                     verts = modelgrid.get_cell_vertices(node)
-                    rstr_data = self.sample_polygon(verts, band)
+                    rstr_data = self.sample_polygon(verts, band).astype(float)
                     msk = np.in1d(rstr_data, self.nodatavals)
                     rstr_data[msk] = np.nan
 
@@ -539,7 +539,7 @@ def __threaded_resampling(
         """
         container.acquire()
         verts = modelgrid.get_cell_vertices(node)
-        rstr_data = self.sample_polygon(verts, band)
+        rstr_data = self.sample_polygon(verts, band).astype(float)
         msk = np.in1d(rstr_data, self.nodatavals)
         rstr_data[msk] = np.nan