Merge branch '269-bugfixes-various' of github.com:rosepearson/GeoFabr…

…ics into 269-bugfixes-various
rosepearson · Jan 9, 2025 · f150bc8 · f150bc8
2 parents b4c1d64 + b0f19d6
commit f150bc8
Showing 1 changed file with 8 additions and 4 deletions.
diff --git a/src/geofabrics/dem.py b/src/geofabrics/dem.py
@@ -894,8 +894,9 @@ def interpolate_ocean_chunked(
                 if chunk_region_to_tile.area.sum() == 0:
                     self.logger.debug(f"\t\tReturning empty tile as out of RIO")
                     delayed_chunked_x.append(
-                        dask.array.empty(
+                        dask.array.full(
                             shape=(len(dim_y), len(dim_x)),
+                            fill_value=numpy.nan,
                             dtype=raster_options["raster_type"],
                         )
                     )
@@ -1229,8 +1230,9 @@ def add_points_within_polygon_chunked(
                 if chunk_region_to_tile.area.sum() == 0:
                     self.logger.debug(f"\t\tReturning empty tile as out of RIO")
                     delayed_chunked_x.append(
-                        dask.array.empty(
+                        dask.array.full(
                             shape=(len(dim_y), len(dim_x)),
+                            fill_value=numpy.nan,
                             dtype=raster_options["raster_type"],
                         )
                     )
@@ -1457,8 +1459,9 @@ def add_points_within_polygon_nearest_chunked(
                 if chunk_region_to_tile.area.sum() == 0:
                     self.logger.debug(f"\t\tReturning empty tile as out of RIO")
                     delayed_chunked_x.append(
-                        dask.array.empty(
+                        dask.array.full(
                             shape=(len(dim_y), len(dim_x)),
+                            fill_value=numpy.nan,
                             dtype=raster_options["raster_type"],
                         )
                     )
@@ -2040,8 +2043,9 @@ def _add_tiled_lidar_chunked(
                             f"\t\tReturning empty tile as no LiDAR or out of ROI"
                         )
                         delayed_chunked_x.append(
-                            dask.array.empty(
+                            dask.array.full(
                                 shape=(len(dim_y), len(dim_x)),
+                                fill_value=numpy.nan,
                                 dtype=raster_options["raster_type"],
                             )
                         )