v0.0.9

CHANGELOG.md

@@ -4,9 +4,13 @@ All notable changes to this project will be documented in this file.
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
-## [0.0.9] - Unreleased
+## [0.0.9] - 2024-02-29
 
+### Fixed
 
+-  `_typer` has a better args management.
+-  `available_indices` returns only indices that can be computed.
+
 ## [0.0.8] - 2024-02-28
 
 ### Added

earthdaily/__init__.py

@@ -1,4 +1,8 @@
 from . import earthdatastore, datasets
 from .accessor import EarthDailyAccessorDataArray, EarthDailyAccessorDataset
+import warnings
+
+# to hide warnings from rioxarray or nano seconds conversion
+warnings.filterwarnings("ignore")
 
 __version__ = "0.0.9"

earthdaily/accessor/__init__.py

@@ -28,12 +28,10 @@ def force(*args, **kwargs):
             idx = 1
             for key, val in func.__annotations__.items():
                 is_kwargs = key in kwargs.keys()
-                if (
-                    val not in _SUPPORTED_DTYPE
-                    or kwargs.get(key, None) is None
-                    and is_kwargs
-                    or len(args) == 1
-                ):
+                if not is_kwargs and idx>= len(args):
+                    continue
+                input_value = kwargs.get(key, None) if is_kwargs else args[idx]
+                if type(input_value) == val:
                     continue
                 if raise_mistype and (
                     val != type(kwargs.get(key))
@@ -256,20 +254,23 @@ def _auto_mapper(self):
                 params[_BAND_MAPPING[v]] = self._obj[v]
         return params
 
-    def available_index(self, details=False):
+    def available_indices(self, details=False):
         mapper = list(self._auto_mapper().keys())
         indices = spyndex.indices
         available_indices = []
         for k, v in indices.items():
             needed_bands = v.bands
+            missing_bands = False
             for needed_band in needed_bands:
                 if needed_band not in mapper:
+                    missing_bands = True
                     break
+            if missing_bands is False:
                 available_indices.append(spyndex.indices[k] if details else k)
         return available_indices
 
     @_typer()
-    def add_index(self, index: list, **kwargs):
+    def add_indices(self, index: list, **kwargs):
         """
         Uses spyndex to compute and add index.
 

examples/field_evolution.py

@@ -44,9 +44,7 @@
 # Plots cube with SCL with at least 50% of clear data
 # ----------------------------------------------------
 pivot_cube = pivot_cube.load()
-pivot_cube.to_array(dim="band").plot.imshow(
-    vmin=0, vmax=0.4, col="time", col_wrap=3
-)
+pivot_cube.ed.plot_rgb(col_wrap=3)
 
 plt.title("Clear cover percent with SCL")
 plt.title("Pivot evolution with SCL masks")

examples/multisensors_cube.py

@@ -43,7 +43,7 @@
 )
 
 # Add the NDVI
-datacube = datacube.ed.add_index('NDVI')
+datacube = datacube.ed.add_indices(['NDVI'])
 
 # Load in memory
 datacube = datacube.load()
@@ -52,18 +52,15 @@
 # See the evolution in RGB
 # -------------------------------------------
 
-datacube[["red", "green", "blue"]].to_array(dim="band").plot.imshow(
-    col="time", col_wrap=3, vmax=0.2
-)
+datacube.ed.plot_rgb(col_wrap=3)
 plt.show()
 
+
 ##############################################################################
 # See the NDVI evolution
 # -------------------------------------------
 
-datacube["NDVI"].plot.imshow(
-    col="time", col_wrap=3, vmin=0, vmax=0.8, cmap="RdYlGn"
-)
+datacube["NDVI"].ed.plot_index(col_wrap=3, vmin=0, vmax=0.8, cmap="Greens")
 plt.show()
 
 ##############################################################################