PyGeoFetch 2.6.3.0__tar.gz → 2.6.5.0__tar.gz

{pygeofetch-2.6.3.0 → pygeofetch-2.6.5.0}/CHANGELOG.md

@@ -20,7 +20,7 @@ date, and marked accordingly rather than filled in with invented
 bullet points.
 ```
 
-## [Unreleased] — v2.6.3.0.4
+## [Unreleased] — v2.6.5.0.4
 
 ```{note}
 Not yet published to PyPI. This entry documents work completed and
@@ -130,7 +130,7 @@ audit and remain open for a follow-up pass.
 
 ---
 
-## [2.6.3.0.4] — 2026-08-25 (latest)
+## [2.6.5.0.4] — 2026-08-25 (latest)
 
 Confirmed via direct fetch of the current PyPI listing.
 
@@ -149,7 +149,7 @@ Confirmed via direct fetch of the current PyPI listing.
 - `SearchQuery.sort_ascending: bool` — earlier examples showed a `sort_order: "asc"/"desc"` string field that doesn't exist on the model (the CLI's `--sort-order` flag does convert to this internally, so CLI docs were unaffected)
 - `DownloadOptions.post_process` requires a list of `PostProcessAction` objects — passing plain strings raises a real `pydantic.ValidationError`
 
-## [2.6.3.0] — 2026-08-24
+## [2.6.5.0] — 2026-08-24
 Real, verified release. Detailed notes not independently recovered from PyPI's release-history page (no distinct per-version notes exposed for this release specifically).
 
 ## [2.6.1] — 2026-08-22
@@ -261,9 +261,9 @@ recovered.
 
 ```{note}
 The previous version of this changelog had an entry labeled
-**"2.6.3.0.4 — 2026-07-12"** describing 22+ providers, Sentinel-1C/1D
+**"2.6.5.0.4 — 2026-07-12"** describing 22+ providers, Sentinel-1C/1D
 support, SLC search with provider routing, precise orbit file
-management, 17 spectral indices, and various bug fixes. **2.6.3.0.4 was
+management, 17 spectral indices, and various bug fixes. **2.6.5.0.4 was
 only ever released once, on 2026-08-25** — it was never released on
 2026-07-12. The date matches this real release, **1.0.8**, exactly.
 That content is reproduced below under the version it actually
@@ -406,8 +406,8 @@ The version of this file uploaded for this update had **5 entries
 covering 5 dates**. PyPI's real release history has **54 releases**.
 Specific issues found and fixed:
 
-1. **Two different entries were both labeled `2.6.3.0.4`**, dated
-   2026-07-29 and 2026-07-12. **2.6.3.0.4 was only ever released once**,
+1. **Two different entries were both labeled `2.6.5.0.4`**, dated
+   2026-07-29 and 2026-07-12. **2.6.5.0.4 was only ever released once**,
    on 2026-08-25 (confirmed directly from the live PyPI listing).
 2. **The 2026-07-12 entry's date exactly matches a real release**
    — version **1.0.8** — and its content (22+ providers, Sentinel-1C/D,

{pygeofetch-2.6.3.0/pygeofetch.egg-info → pygeofetch-2.6.5.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pygeofetch
-Version: 2.6.3.0
+Version: 2.6.5.0
 Summary: Universal satellite data pipeline + geospatial processing platform
 Author-email: Samuel Appiah Kubi <appiahkubis14@gmail.com>
 License: MIT License
@@ -26,7 +26,7 @@ License: MIT License
         SOFTWARE.
         
 Project-URL: Homepage, https://www.eocoreint.com/
-Project-URL: Documentation, https://appiahkubis14.github.io/pygeofetch-docs/
+Project-URL: Documentation, https://pygeofetch.readthedocs.io/en/latest/ 
 Project-URL: Repository, https://github.com/EOCoreINT/pygeofetch
 Project-URL: Bug Tracker, https://github.com/EOCoreINT/pygeofetch/issues
 Keywords: satellite,remote-sensing,GIS,earth-observation,geospatial,SAR,Sentinel,Landsat
@@ -43,6 +43,8 @@ Requires-Python: >=3.9
 Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: httpx>=0.27
+Requires-Dist: pydantic
+Requires-Dist: rich
 Requires-Dist: click>=8.1
 Requires-Dist: pydantic>=2.5
 Requires-Dist: pydantic-settings>=2.1
@@ -369,7 +371,7 @@ If you use PyGeoFetch in your research or operational work, please cite it using
   title        = {PyGeoFetch: A Unified Python Framework for Multi-Provider Satellite Data Acquisition, Pre-download Quality Control, and Geospatial Processing},
   year         = {2026},
   publisher    = {Zenodo},
-  version      = {2.6.3.0.1},
+  version      = {2.6.5.0.1},
   doi          = {10.5281/zenodo.22087230},
   url          = {https://doi.org/10.5281/zenodo.22087230}
 }

{pygeofetch-2.6.3.0 → pygeofetch-2.6.5.0}/README.md

@@ -229,7 +229,7 @@ If you use PyGeoFetch in your research or operational work, please cite it using
   title        = {PyGeoFetch: A Unified Python Framework for Multi-Provider Satellite Data Acquisition, Pre-download Quality Control, and Geospatial Processing},
   year         = {2026},
   publisher    = {Zenodo},
-  version      = {2.6.3.0.1},
+  version      = {2.6.5.0.1},
   doi          = {10.5281/zenodo.22087230},
   url          = {https://doi.org/10.5281/zenodo.22087230}
 }
@@ -263,4 +263,4 @@ pip install -e ".[test]" && pytest tests/ -v
 
 ## License
 
-MIT — see [LICENSE](LICENSE). © 2026 Samuel Appiah Kubi.
+MIT — see [LICENSE](LICENSE). © 2026 Samuel Appiah Kubi.

{pygeofetch-2.6.3.0 → pygeofetch-2.6.5.0}/pygeofetch/__init__.py

@@ -30,13 +30,16 @@ CLI::
     pygeofetch doctor
 """
 
-__version__ = "2.6.3.0.4"
+__version__ = "2.6.5.0.4"
 __author__ = "PyGeoFetch Contributors"
 __license__ = "MIT"
 
 from pygeofetch.core.engine import PyGeoFetch  # noqa: F401
+from pygeofetch.core.logging import (
+    quiet_logs,  # noqa: F401  (notebook-friendly logging)
+)
 
-__all__ = ["PyGeoFetch", "__version__"]
+__all__ = ["PyGeoFetch", "quiet_logs", "__version__"]
 
 
 # Optional modules — lazy imported to avoid hard dependency errors

{pygeofetch-2.6.3.0 → pygeofetch-2.6.5.0}/pygeofetch/cli/classify_commands.py

@@ -45,7 +45,14 @@ def classify() -> None:
 
 @classify.command("kmeans")
 @click.argument("inputs", nargs=-1, type=click.Path(exists=True))
-@click.option("--clusters", "-k", default=5, show_default=True, type=int, help="Number of clusters.")
+@click.option(
+    "--clusters",
+    "-k",
+    default=5,
+    show_default=True,
+    type=int,
+    help="Number of clusters.",
+)
 @click.option("--output", "-o", default=None)
 @click.option("--random-state", default=42, show_default=True, type=int)
 # @clip_options
@@ -59,21 +66,62 @@ def kmeans_cmd(inputs, clusters, output, random_state, bbox, geometry, geometry_
         sys.exit(1)
     e = _engine()
     r = e.classify.kmeans(
-        inputs=list(inputs), n_clusters=clusters, output=output, random_state=random_state,
-        bbox=bbox, geometry=geometry, geometry_crs=geometry_crs,
+        inputs=list(inputs),
+        n_clusters=clusters,
+        output=output,
+        random_state=random_state,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
     )
     _pr(r, f"K-means ({clusters} clusters)")
 
 
 @classify.command("train")
 @click.argument("inputs", nargs=-1, type=click.Path(exists=True))
-@click.option("--training-data", required=True, type=click.Path(exists=True), help="Vector file (GeoJSON/Shapefile) with labeled training geometries.")
-@click.option("--class-field", required=True, help="Attribute column holding the real class label.")
-@click.option("--output-model", "-o", required=True, help="Where to save the trained model.")
-@click.option("--model-type", default="random_forest", type=click.Choice(["random_forest", "svm"]), show_default=True)
-@click.option("--n-estimators", default=100, show_default=True, type=int, help="Random Forest tree count.")
-@click.option("--test-size", default=0.2, show_default=True, type=float, help="Fraction of samples held out for real accuracy reporting.")
-def train_cmd(inputs, training_data, class_field, output_model, model_type, n_estimators, test_size):
+@click.option(
+    "--training-data",
+    required=True,
+    type=click.Path(exists=True),
+    help="Vector file (GeoJSON/Shapefile) with labeled training geometries.",
+)
+@click.option(
+    "--class-field",
+    required=True,
+    help="Attribute column holding the real class label.",
+)
+@click.option(
+    "--output-model", "-o", required=True, help="Where to save the trained model."
+)
+@click.option(
+    "--model-type",
+    default="random_forest",
+    type=click.Choice(["random_forest", "svm"]),
+    show_default=True,
+)
+@click.option(
+    "--n-estimators",
+    default=100,
+    show_default=True,
+    type=int,
+    help="Random Forest tree count.",
+)
+@click.option(
+    "--test-size",
+    default=0.2,
+    show_default=True,
+    type=float,
+    help="Fraction of samples held out for real accuracy reporting.",
+)
+def train_cmd(
+    inputs,
+    training_data,
+    class_field,
+    output_model,
+    model_type,
+    n_estimators,
+    test_size,
+):
     """Train a supervised classifier from labeled vector training samples.
 
     Example: pygeofetch classify train B02.tif B03.tif B04.tif B08.tif
@@ -84,8 +132,12 @@ def train_cmd(inputs, training_data, class_field, output_model, model_type, n_es
         sys.exit(1)
     e = _engine()
     r = e.classify.train(
-        inputs=list(inputs), training_data=training_data, class_field=class_field,
-        output_model=output_model, model_type=model_type, n_estimators=n_estimators,
+        inputs=list(inputs),
+        training_data=training_data,
+        class_field=class_field,
+        output_model=output_model,
+        model_type=model_type,
+        n_estimators=n_estimators,
         test_size=test_size,
     )
     _pr(r, f"train ({model_type})")
@@ -93,7 +145,13 @@ def train_cmd(inputs, training_data, class_field, output_model, model_type, n_es
 
 @classify.command("predict")
 @click.argument("inputs", nargs=-1, type=click.Path(exists=True))
-@click.option("--model", "-m", required=True, type=click.Path(exists=True), help="A model trained via 'classify train'.")
+@click.option(
+    "--model",
+    "-m",
+    required=True,
+    type=click.Path(exists=True),
+    help="A model trained via 'classify train'.",
+)
 @click.option("--output", "-o", default=None)
 # @clip_options
 def predict_cmd(inputs, model, output, bbox, geometry, geometry_crs):
@@ -106,7 +164,11 @@ def predict_cmd(inputs, model, output, bbox, geometry, geometry_crs):
         sys.exit(1)
     e = _engine()
     r = e.classify.predict(
-        inputs=list(inputs), model=model, output=output,
-        bbox=bbox, geometry=geometry, geometry_crs=geometry_crs,
+        inputs=list(inputs),
+        model=model,
+        output=output,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
     )
     _pr(r, "predict")

{pygeofetch-2.6.3.0 → pygeofetch-2.6.5.0}/pygeofetch/cli/index_commands.py

@@ -22,9 +22,13 @@ def _parse_bbox(ctx, param, value):
     try:
         parts = [float(x.strip()) for x in value.split(",")]
     except ValueError:
-        raise click.BadParameter("must be 4 comma-separated numbers: minx,miny,maxx,maxy") from None
+        raise click.BadParameter(
+            "must be 4 comma-separated numbers: minx,miny,maxx,maxy"
+        ) from None
     if len(parts) != 4:
-        raise click.BadParameter("must be exactly 4 comma-separated numbers: minx,miny,maxx,maxy")
+        raise click.BadParameter(
+            "must be exactly 4 comma-separated numbers: minx,miny,maxx,maxy"
+        )
     return tuple(parts)
 
 
@@ -33,15 +37,21 @@ def clip_options(f):
     index command whose underlying method supports post-computation
     clipping."""
     f = click.option(
-        "--bbox", callback=_parse_bbox, default=None,
+        "--bbox",
+        callback=_parse_bbox,
+        default=None,
         help="Clip the result to minx,miny,maxx,maxy (comma-separated).",
     )(f)
     f = click.option(
-        "--geometry", default=None, type=click.Path(exists=True),
+        "--geometry",
+        default=None,
+        type=click.Path(exists=True),
         help="Clip the result to a GeoJSON file's geometry instead of a bbox.",
     )(f)
     f = click.option(
-        "--geometry-crs", default="EPSG:4326", show_default=True,
+        "--geometry-crs",
+        default="EPSG:4326",
+        show_default=True,
         help="CRS of --bbox/--geometry coordinates.",
     )(f)
     return f
@@ -82,10 +92,17 @@ def index() -> None:
 )
 @click.option("--output", "-o", default=None)
 @click.option(
-    "--chunked", is_flag=True,
+    "--chunked",
+    is_flag=True,
     help="Process in memory-safe tiles for large rasters (e.g. drone orthomosaics) that exceed available RAM or the 4GB GeoTIFF limit.",
 )
-@click.option("--tile-size", default=1024, show_default=True, type=int, help="Tile edge length in pixels when --chunked is set.")
+@click.option(
+    "--tile-size",
+    default=1024,
+    show_default=True,
+    type=int,
+    help="Tile edge length in pixels when --chunked is set.",
+)
 @clip_options
 def ndvi_cmd(red, nir, output, chunked, tile_size, bbox, geometry, geometry_crs):
     """NDVI — Normalized Difference Vegetation Index. Range: -1 to +1. Vegetation > 0.3.
@@ -95,7 +112,16 @@ def ndvi_cmd(red, nir, output, chunked, tile_size, bbox, geometry, geometry_crs)
     real, embedded band descriptions.
     """
     e = _engine()
-    r = e.indices.ndvi(red=red, nir=nir, output=output, chunked=chunked, tile_size=tile_size, bbox=bbox, geometry=geometry, geometry_crs=geometry_crs)
+    r = e.indices.ndvi(
+        red=red,
+        nir=nir,
+        output=output,
+        chunked=chunked,
+        tile_size=tile_size,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
+    )
     _pr(r, "NDVI")
 
 
@@ -103,15 +129,20 @@ def ndvi_cmd(red, nir, output, chunked, tile_size, bbox, geometry, geometry_crs)
 @click.option(
     "--rededge", required=True, type=click.Path(exists=True), help="Red edge band"
 )
-@click.option(
-    "--nir", required=True, type=click.Path(exists=True), help="NIR band"
-)
+@click.option("--nir", required=True, type=click.Path(exists=True), help="NIR band")
 @click.option("--output", "-o", default=None)
 @click.option(
-    "--chunked", is_flag=True,
+    "--chunked",
+    is_flag=True,
     help="Process in memory-safe tiles for large rasters that exceed available RAM or the 4GB GeoTIFF limit.",
 )
-@click.option("--tile-size", default=1024, show_default=True, type=int, help="Tile edge length in pixels when --chunked is set.")
+@click.option(
+    "--tile-size",
+    default=1024,
+    show_default=True,
+    type=int,
+    help="Tile edge length in pixels when --chunked is set.",
+)
 @clip_options
 def ndre_cmd(rededge, nir, output, chunked, tile_size, bbox, geometry, geometry_crs):
     """NDRE — Normalized Difference Red Edge Index. Range: -1 to +1.
@@ -122,7 +153,16 @@ def ndre_cmd(rededge, nir, output, chunked, tile_size, bbox, geometry, geometry_
     both bands from one merged multi-band image.
     """
     e = _engine()
-    r = e.indices.ndre(rededge=rededge, nir=nir, output=output, chunked=chunked, tile_size=tile_size, bbox=bbox, geometry=geometry, geometry_crs=geometry_crs)
+    r = e.indices.ndre(
+        rededge=rededge,
+        nir=nir,
+        output=output,
+        chunked=chunked,
+        tile_size=tile_size,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
+    )
     _pr(r, "NDRE")
 
 
@@ -131,7 +171,9 @@ def ndre_cmd(rededge, nir, output, chunked, tile_size, bbox, geometry, geometry_
 @click.option("--red", required=True, type=click.Path(exists=True))
 @click.option("--nir", required=True, type=click.Path(exists=True))
 @click.option("--output", "-o", default=None)
-@click.option("--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters.")
+@click.option(
+    "--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters."
+)
 @click.option("--tile-size", default=1024, show_default=True, type=int)
 @clip_options
 def evi_cmd(blue, red, nir, output, chunked, tile_size, bbox, geometry, geometry_crs):
@@ -141,7 +183,17 @@ def evi_cmd(blue, red, nir, output, chunked, tile_size, bbox, geometry, geometry
     three bands from one merged multi-band image.
     """
     e = _engine()
-    r = e.indices.evi(blue=blue, red=red, nir=nir, output=output, chunked=chunked, tile_size=tile_size, bbox=bbox, geometry=geometry, geometry_crs=geometry_crs)
+    r = e.indices.evi(
+        blue=blue,
+        red=red,
+        nir=nir,
+        output=output,
+        chunked=chunked,
+        tile_size=tile_size,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
+    )
     _pr(r, "EVI")
 
 
@@ -150,17 +202,31 @@ def evi_cmd(blue, red, nir, output, chunked, tile_size, bbox, geometry, geometry
 @click.option("--nir", required=True, type=click.Path(exists=True))
 @click.option("--soil-l", default=0.5, show_default=True, type=float)
 @click.option("--output", "-o", default=None)
-@click.option("--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters.")
+@click.option(
+    "--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters."
+)
 @click.option("--tile-size", default=1024, show_default=True, type=int)
 @clip_options
-def savi_cmd(red, nir, soil_l, output, chunked, tile_size, bbox, geometry, geometry_crs):
+def savi_cmd(
+    red, nir, soil_l, output, chunked, tile_size, bbox, geometry, geometry_crs
+):
     """SAVI — Soil Adjusted Vegetation Index. Reduces soil brightness influence.
 
     Pass the SAME path to --red and --nir to auto-detect both bands
     from one merged multi-band image.
     """
     e = _engine()
-    r = e.indices.savi(red=red, nir=nir, L=soil_l, output=output, chunked=chunked, tile_size=tile_size, bbox=bbox, geometry=geometry, geometry_crs=geometry_crs)
+    r = e.indices.savi(
+        red=red,
+        nir=nir,
+        L=soil_l,
+        output=output,
+        chunked=chunked,
+        tile_size=tile_size,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
+    )
     _pr(r, "SAVI")
 
 
@@ -168,7 +234,9 @@ def savi_cmd(red, nir, soil_l, output, chunked, tile_size, bbox, geometry, geome
 @click.option("--green", required=True, type=click.Path(exists=True))
 @click.option("--nir", required=True, type=click.Path(exists=True))
 @click.option("--output", "-o", default=None)
-@click.option("--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters.")
+@click.option(
+    "--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters."
+)
 @click.option("--tile-size", default=1024, show_default=True, type=int)
 @clip_options
 def ndwi_cmd(green, nir, output, chunked, tile_size, bbox, geometry, geometry_crs):
@@ -178,7 +246,16 @@ def ndwi_cmd(green, nir, output, chunked, tile_size, bbox, geometry, geometry_cr
     from one merged multi-band image.
     """
     e = _engine()
-    r = e.indices.ndwi(green=green, nir=nir, output=output, chunked=chunked, tile_size=tile_size, bbox=bbox, geometry=geometry, geometry_crs=geometry_crs)
+    r = e.indices.ndwi(
+        green=green,
+        nir=nir,
+        output=output,
+        chunked=chunked,
+        tile_size=tile_size,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
+    )
     _pr(r, "NDWI")
 
 
@@ -186,7 +263,9 @@ def ndwi_cmd(green, nir, output, chunked, tile_size, bbox, geometry, geometry_cr
 @click.option("--green", required=True, type=click.Path(exists=True))
 @click.option("--swir1", required=True, type=click.Path(exists=True))
 @click.option("--output", "-o", default=None)
-@click.option("--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters.")
+@click.option(
+    "--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters."
+)
 @click.option("--tile-size", default=1024, show_default=True, type=int)
 @clip_options
 def mndwi_cmd(green, swir1, output, chunked, tile_size, bbox, geometry, geometry_crs):
@@ -196,7 +275,16 @@ def mndwi_cmd(green, swir1, output, chunked, tile_size, bbox, geometry, geometry
     bands from one merged multi-band image.
     """
     e = _engine()
-    r = e.indices.mndwi(green=green, swir1=swir1, output=output, chunked=chunked, tile_size=tile_size, bbox=bbox, geometry=geometry, geometry_crs=geometry_crs)
+    r = e.indices.mndwi(
+        green=green,
+        swir1=swir1,
+        output=output,
+        chunked=chunked,
+        tile_size=tile_size,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
+    )
     _pr(r, "MNDWI")
 
 
@@ -204,7 +292,9 @@ def mndwi_cmd(green, swir1, output, chunked, tile_size, bbox, geometry, geometry
 @click.option("--nir", required=True, type=click.Path(exists=True))
 @click.option("--swir1", required=True, type=click.Path(exists=True))
 @click.option("--output", "-o", default=None)
-@click.option("--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters.")
+@click.option(
+    "--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters."
+)
 @click.option("--tile-size", default=1024, show_default=True, type=int)
 @clip_options
 def ndbi_cmd(nir, swir1, output, chunked, tile_size, bbox, geometry, geometry_crs):
@@ -214,7 +304,16 @@ def ndbi_cmd(nir, swir1, output, chunked, tile_size, bbox, geometry, geometry_cr
     from one merged multi-band image.
     """
     e = _engine()
-    r = e.indices.ndbi(nir=nir, swir1=swir1, output=output, chunked=chunked, tile_size=tile_size, bbox=bbox, geometry=geometry, geometry_crs=geometry_crs)
+    r = e.indices.ndbi(
+        nir=nir,
+        swir1=swir1,
+        output=output,
+        chunked=chunked,
+        tile_size=tile_size,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
+    )
     _pr(r, "NDBI")
 
 
@@ -222,7 +321,9 @@ def ndbi_cmd(nir, swir1, output, chunked, tile_size, bbox, geometry, geometry_cr
 @click.option("--green", required=True, type=click.Path(exists=True))
 @click.option("--swir1", required=True, type=click.Path(exists=True))
 @click.option("--output", "-o", default=None)
-@click.option("--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters.")
+@click.option(
+    "--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters."
+)
 @click.option("--tile-size", default=1024, show_default=True, type=int)
 @clip_options
 def ndsi_cmd(green, swir1, output, chunked, tile_size, bbox, geometry, geometry_crs):
@@ -232,7 +333,16 @@ def ndsi_cmd(green, swir1, output, chunked, tile_size, bbox, geometry, geometry_
     bands from one merged multi-band image.
     """
     e = _engine()
-    r = e.indices.ndsi(green=green, swir1=swir1, output=output, chunked=chunked, tile_size=tile_size, bbox=bbox, geometry=geometry, geometry_crs=geometry_crs)
+    r = e.indices.ndsi(
+        green=green,
+        swir1=swir1,
+        output=output,
+        chunked=chunked,
+        tile_size=tile_size,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
+    )
     _pr(r, "NDSI")
 
 
@@ -240,7 +350,9 @@ def ndsi_cmd(green, swir1, output, chunked, tile_size, bbox, geometry, geometry_
 @click.option("--nir", required=True, type=click.Path(exists=True))
 @click.option("--swir1", required=True, type=click.Path(exists=True))
 @click.option("--output", "-o", default=None)
-@click.option("--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters.")
+@click.option(
+    "--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters."
+)
 @click.option("--tile-size", default=1024, show_default=True, type=int)
 @clip_options
 def ndmi_cmd(nir, swir1, output, chunked, tile_size, bbox, geometry, geometry_crs):
@@ -250,7 +362,16 @@ def ndmi_cmd(nir, swir1, output, chunked, tile_size, bbox, geometry, geometry_cr
     from one merged multi-band image.
     """
     e = _engine()
-    r = e.indices.ndmi(nir=nir, swir1=swir1, output=output, chunked=chunked, tile_size=tile_size, bbox=bbox, geometry=geometry, geometry_crs=geometry_crs)
+    r = e.indices.ndmi(
+        nir=nir,
+        swir1=swir1,
+        output=output,
+        chunked=chunked,
+        tile_size=tile_size,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
+    )
     _pr(r, "NDMI")
 
 
@@ -258,7 +379,9 @@ def ndmi_cmd(nir, swir1, output, chunked, tile_size, bbox, geometry, geometry_cr
 @click.option("--nir", required=True, type=click.Path(exists=True))
 @click.option("--swir2", required=True, type=click.Path(exists=True))
 @click.option("--output", "-o", default=None)
-@click.option("--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters.")
+@click.option(
+    "--chunked", is_flag=True, help="Process in memory-safe tiles for large rasters."
+)
 @click.option("--tile-size", default=1024, show_default=True, type=int)
 @clip_options
 def nbr_cmd(nir, swir2, output, chunked, tile_size, bbox, geometry, geometry_crs):
@@ -268,7 +391,16 @@ def nbr_cmd(nir, swir2, output, chunked, tile_size, bbox, geometry, geometry_crs
     from one merged multi-band image.
     """
     e = _engine()
-    r = e.indices.nbr(nir=nir, swir2=swir2, output=output, chunked=chunked, tile_size=tile_size, bbox=bbox, geometry=geometry, geometry_crs=geometry_crs)
+    r = e.indices.nbr(
+        nir=nir,
+        swir2=swir2,
+        output=output,
+        chunked=chunked,
+        tile_size=tile_size,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
+    )
     _pr(r, "NBR")
 
 
@@ -388,15 +520,33 @@ def albedo_cmd(inputs, sensor, output):
 )
 @click.option("--output", "-o", default=None)
 @click.option(
-    "--chunked", is_flag=True,
+    "--chunked",
+    is_flag=True,
     help="Process in memory-safe tiles for large rasters (e.g. drone orthomosaics) that exceed available RAM or the 4GB GeoTIFF limit.",
 )
-@click.option("--tile-size", default=1024, show_default=True, type=int, help="Tile edge length in pixels when --chunked is set.")
+@click.option(
+    "--tile-size",
+    default=1024,
+    show_default=True,
+    type=int,
+    help="Tile edge length in pixels when --chunked is set.",
+)
 @clip_options
-def band_math_cmd(inputs, expr, output, chunked, tile_size, bbox, geometry, geometry_crs):
+def band_math_cmd(
+    inputs, expr, output, chunked, tile_size, bbox, geometry, geometry_crs
+):
     """Arbitrary band arithmetic. Example: --expr '(B[1]-B[0])/(B[1]+B[0])'."""
     e = _engine()
-    r = e.indices.band_math(inputs=list(inputs), expression=expr, output=output, chunked=chunked, tile_size=tile_size, bbox=bbox, geometry=geometry, geometry_crs=geometry_crs)
+    r = e.indices.band_math(
+        inputs=list(inputs),
+        expression=expr,
+        output=output,
+        chunked=chunked,
+        tile_size=tile_size,
+        bbox=bbox,
+        geometry=geometry,
+        geometry_crs=geometry_crs,
+    )
     _pr(r, f"band-math ({expr[:40]})")
 
 
@@ -407,4 +557,4 @@ def stack_cmd(inputs, output):
     """Stack single-band rasters into a multi-band GeoTIFF."""
     e = _engine()
     r = e.indices.stack(inputs=list(inputs), output=output)
-    _pr(r, f"stack ({len(inputs)} bands)")
+    _pr(r, f"stack ({len(inputs)} bands)")