fixed esri_land bug (time mean instead of sel=0)

Author: lukegre

cryogrid_pytools/data/from_planetary_computer.py

@@ -222,12 +222,16 @@ def get_esri_land_cover(bbox_WSEN: tuple, res: int = 30, epsg: int = 32643):
         datetime="2020-01-01/2020-12-31",
     )
 
-    da = _stackstac.stack(
-        items,
-        epsg=epsg,
-        resolution=res,
-        bounds_latlon=bbox_WSEN,
-    ).isel(time=0, band=0)
+    da = (
+        _stackstac.stack(
+            items,
+            epsg=epsg,
+            resolution=res,
+            bounds_latlon=bbox_WSEN,
+        )
+        .mean("time")
+        .isel(band=0)
+    )
 
     da = da.rename("land_use_and_cover").assign_attrs(
         long_name="Land use and land cover (LULC) from ESA Sentinel-2",

pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "hatchling.build"
 
 [project]
 name = "cryogrid_pytools"
-version = "0.3.13"
+version = "0.3.14"
 description = "Tools to read in CryoGrid MATLAB data to Python and create forcing data"
 readme = "README.md"
 license = "MIT"