Fix for avg filtering non-nan masked raster

noaa-ocs-modeling · Aug 24, 2023 · ab591ea · ab591ea
1 parent 17321c4
commit ab591ea
Showing 1 changed file with 6 additions and 1 deletion.
diff --git a/ocsmesh/raster.py b/ocsmesh/raster.py
@@ -965,8 +965,13 @@ def average_filter(
                     mask_below = ma.getdata(outband) < drop_below
                     outband[mask_below] = np.nan
 
+                # Since the nbmean checks for isnan, other fill values
+                # will be treated as normal numbers
                 outband_new = generic_filter(
-                    outband, LowLevelCallable(nbmean.ctypes), size=size)
+                    outband.filled(np.nan),
+                    LowLevelCallable(nbmean.ctypes),
+                    size=size
+                )
 
                 # Mask raster based on result of filter?
                 if drop_above is not None: