diff --git a/autotest/t065_test_gridintersect.py b/autotest/t065_test_gridintersect.py
index f37a7f2466..768955b9e7 100644
--- 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
index 3aa6a64559..6d4e856a66 100644
--- 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