Only set fill value as nan for floating point arrays (#113)

carbonplan · Mar 13, 2024 · 8e56c02 · 8e56c02
1 parent 0c1c734
commit 8e56c02
Showing 1 changed file with 10 additions and 2 deletions.
diff --git a/ndpyramid/core.py b/ndpyramid/core.py
@@ -8,7 +8,12 @@
 import xarray as xr
 
 from .common import Projection
-from .utils import add_metadata_and_zarr_encoding, get_levels, get_version, multiscales_template
+from .utils import (
+    add_metadata_and_zarr_encoding,
+    get_levels,
+    get_version,
+    multiscales_template,
+)
 
 
 def pyramid_coarsen(
@@ -132,7 +137,10 @@ def pyramid_reproject(
         dst_transform = projection_model.transform(dim=dim)
 
         def reproject(da, var):
-            da.encoding['_FillValue'] = np.nan
+            # Set float FillValue that rasterio will use to initialize array
+            if da.encoding.get('_FillValue') is None and np.issubdtype(da.dtype, np.floating):
+                da.encoding['_FillValue'] = np.nan
+
             da = da.rio.reproject(
                 projection_model._crs,
                 resampling=Resampling[resampling_dict[var]],