Merge pull request #35 from EOPF-Sample-Service/konstntokas-xxx-bugfix_sen2

Bugfix in cubegen of sen2

Author: konstntokas

CHANGES.md

@@ -1,4 +1,8 @@
-## Changes in 0.1.1 (under development)
+## Changes in 0.1.1
+
+* Fixed a bug in Sentinel-2 cube generation where, during mosaicking of adjacent tiles 
+  from the same solar day, data from one spectral band would overwrite all other bands 
+  in the final cube.
 
 
 ## Changes in 0.1.0

environment.yml

@@ -10,7 +10,7 @@ dependencies:
   - pystac
   - pystac-client
   - xarray
-  - xarray-eopf
+  - xarray-eopf >=0.1.1
   - xcube >=1.11.0
   # Development Dependencies - Tools
   - black

pyproject.toml

@@ -47,7 +47,7 @@ dependencies = [
   "pystac",
   "pystac-client",
   "xarray",
-  "xarray-eopf",
+  "xarray-eopf>=0.1.1",
   "xcube>=1.11.0"
 ]
 

xcube_eopf/prodhandlers/sentinel2.py

@@ -709,14 +709,14 @@ def _create_empty_dataset(
     )
     x = np.arange(x_start + half_res, x_end, spatial_res)
 
-    empty_data = da.full(
-        (grouped_items.sizes["time"], len(y), len(x)),
-        np.nan,
-        chunks=(1, _TILE_SIZE, _TILE_SIZE),
-    )
+    chunks = (1, _TILE_SIZE, _TILE_SIZE)
+    shape = (grouped_items.sizes["time"], len(y), len(x))
     ds = xr.Dataset(
         {
-            key: (("time", "y", "x"), empty_data.astype(var))
+            key: (
+                ("time", "y", "x"),
+                da.full(shape, np.nan, dtype=var.dtype, chunks=chunks),
+            )
             for (key, var) in sample_ds.data_vars.items()
         },
         coords={

xcube_eopf/version.py

@@ -2,4 +2,4 @@
 #  Permissions are hereby granted under the terms of the Apache 2.0 License:
 #  https://opensource.org/license/apache-2-0.
 
-version = "0.1.1.dev0"
+version = "0.1.1"