PyGeoFetch 2.5.0__py3-none-any.whl → 2.6.1__py3-none-any.whl

pygeofetch/__init__.py

@@ -30,7 +30,7 @@ CLI::
     pygeofetch doctor
 """
 
-__version__ = "2.5.0"
+__version__ = "2.6.1"
 __author__ = "PyGeoFetch Contributors"
 __license__ = "MIT"
 

pygeofetch/insar/preflight.py

@@ -610,11 +610,14 @@ class PreflightGate:
         advisory," never "risk confirmed absent."
         """
         import tempfile
+        import shutil
+        import importlib
         try:
-            from pygeofetch.providers.copernicus_nodes import fetch_annotation_zip
-            from pygeofetch.core.orbits import fetch_orbit_file
-            from pygeofetch.insar.stack_selection import select_burst_synchronized_dates
-            from pygeofetch.insar.geolocation import geodetic_to_ecef
+            for _mod_name in (
+                "pygeofetch.providers.copernicus_nodes", "pygeofetch.core.orbits",
+                "pygeofetch.insar.stack_selection", "pygeofetch.insar.geolocation",
+            ):
+                importlib.import_module(_mod_name)
         except ImportError as exc:
             logger.info("Real pre-download burst check unavailable (%s).", exc)
             return None
@@ -625,6 +628,42 @@ class PreflightGate:
             logger.info("Real pre-download burst check: couldn't create a scratch dir (%s).", exc)
             return None
 
+        # REAL BUG FOUND AND FIXED: work_dir was never removed on any exit
+        # path -- not on success, not on the early "too few usable dates"
+        # return, not on the exception fallback. Every real call to this
+        # method (i.e. every real preflight run with attempt_real_burst_
+        # check=True, which is the default) left behind a full scratch
+        # directory of real annotation mini-zips and real, full orbit
+        # files (~4MB+ each) that was never cleaned up. Re-running the
+        # same preflight cell repeatedly while iterating on a notebook --
+        # completely normal, everyday usage -- accumulates one orphaned
+        # directory per run, unbounded, on whatever drive the OS temp
+        # folder lives on (on Windows, typically the same small system
+        # drive as everything else). A real, concrete, plausible path to
+        # exhausting that drive, which is a well-known way to make a
+        # Jupyter kernel become unresponsive or fail to save cleanly --
+        # matching a real report of a notebook going inactive with blank
+        # cell outputs immediately after a preflight run. The real body
+        # of this method is now factored into _run_real_burst_family_check
+        # below, wrapped here so work_dir is removed on every real exit
+        # path, success or failure, not just the ones tested before.
+        try:
+            return self._run_real_burst_family_check(selected, work_dir)
+        finally:
+            shutil.rmtree(work_dir, ignore_errors=True)
+
+    def _run_real_burst_family_check(self, selected: List[Any], work_dir: Path):
+        """
+        The real body of _try_real_burst_family_check(), factored out so
+        that method's own try/finally can guarantee work_dir cleanup on
+        every real exit path here -- including the early returns --
+        without duplicating the rmtree call at each one.
+        """
+        from pygeofetch.providers.copernicus_nodes import fetch_annotation_zip
+        from pygeofetch.core.orbits import fetch_orbit_file
+        from pygeofetch.insar.stack_selection import select_burst_synchronized_dates
+        from pygeofetch.insar.geolocation import geodetic_to_ecef
+
         safe_zips: Dict[str, Any] = {}
         orbit_files: Dict[str, Any] = {}
 

pygeofetch/meta.yaml

@@ -1,4 +1,4 @@
-{% set version = "2.5.0" %}
+{% set version = "2.6.1" %}
 
 package:
   name: pygeofetch

pygeofetch-2.6.1.dist-info/METADATA

@@ -0,0 +1,407 @@
+Metadata-Version: 2.4
+Name: pygeofetch
+Version: 2.6.1
+Summary: Universal satellite data pipeline + geospatial processing platform
+Author-email: Samuel Appiah Kubi <appiahkubis14@gmail.com>
+License: MIT License
+        
+        Copyright (c) 2026 Samuel Appiah Kubi
+        
+        Permission is hereby granted, free of charge, to any person obtaining a copy
+        of this software and associated documentation files (the "Software"), to deal
+        in the Software without restriction, including without limitation the rights
+        to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+        copies of the Software, and to permit persons to whom the Software is
+        furnished to do so, subject to the following conditions:
+        
+        The above copyright notice and this permission notice shall be included in all
+        copies or substantial portions of the Software.
+        
+        THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+        IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+        FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+        AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+        LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+        OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+        SOFTWARE.
+        
+Project-URL: Homepage, https://www.eocoreint.com/
+Project-URL: Documentation, https://appiahkubis14.github.io/pygeofetch-docs/
+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
+Classifier: Development Status :: 5 - Production/Stable
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3.9
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Classifier: Topic :: Scientific/Engineering :: GIS
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: httpx>=0.27
+Requires-Dist: click>=8.1
+Requires-Dist: pydantic>=2.5
+Requires-Dist: pydantic-settings>=2.1
+Requires-Dist: PyYAML>=6.0
+Requires-Dist: cryptography>=42.0
+Requires-Dist: keyring>=24.3
+Requires-Dist: tenacity>=8.2
+Requires-Dist: python-dateutil>=2.9
+Requires-Dist: requests>=2.31
+Requires-Dist: anyio>=4.2
+Requires-Dist: numpy>=1.26
+Requires-Dist: tqdm>=4.66
+Provides-Extra: geo
+Requires-Dist: rasterio>=1.3; extra == "geo"
+Requires-Dist: geopandas>=0.14; extra == "geo"
+Requires-Dist: pyarrow>=15.0; extra == "geo"
+Requires-Dist: shapely>=2.0; extra == "geo"
+Requires-Dist: pyproj>=3.6; extra == "geo"
+Provides-Extra: processor
+Requires-Dist: rasterio>=1.3; extra == "processor"
+Requires-Dist: rioxarray>=0.15; extra == "processor"
+Requires-Dist: xarray>=2024.1; extra == "processor"
+Requires-Dist: eoreader>=0.24; extra == "processor"
+Requires-Dist: spyndex>=0.5; extra == "processor"
+Requires-Dist: scipy>=1.12; extra == "processor"
+Requires-Dist: scikit-image>=0.23; extra == "processor"
+Requires-Dist: geopandas>=0.14; extra == "processor"
+Provides-Extra: providers
+Requires-Dist: eodag>=3.0; extra == "providers"
+Provides-Extra: insar
+Requires-Dist: sarxarray>=0.1; extra == "insar"
+Requires-Dist: scipy>=1.12; extra == "insar"
+Provides-Extra: ost
+Requires-Dist: opensartoolkit>=3.0; extra == "ost"
+Provides-Extra: viz
+Requires-Dist: leafmap>=0.30; extra == "viz"
+Requires-Dist: folium>=0.16; extra == "viz"
+Requires-Dist: matplotlib>=3.7; extra == "viz"
+Requires-Dist: localtileserver>=1.0.0; extra == "viz"
+Requires-Dist: ipysheet>=0.6.0; extra == "viz"
+Requires-Dist: ipyleaflet>=0.17.0; extra == "viz"
+Requires-Dist: ipywidgets>=8.0.0; extra == "viz"
+Requires-Dist: ipyevents>=2.0.2; extra == "viz"
+Requires-Dist: ipyfilechooser>=0.6.0; extra == "viz"
+Requires-Dist: sidecar>=0.6.0; extra == "viz"
+Requires-Dist: bokeh>=3.0.0; extra == "viz"
+Requires-Dist: seaborn>=0.12.0; extra == "viz"
+Requires-Dist: plotly>=5.14.0; extra == "viz"
+Requires-Dist: whitebox>=2.4.0; extra == "viz"
+Requires-Dist: whiteboxgui>=2.4.0; extra == "viz"
+Requires-Dist: geemap>=0.28.0; extra == "viz"
+Requires-Dist: osmnx>=1.9.4; extra == "viz"
+Requires-Dist: contextily>=1.9.4; extra == "viz"
+Requires-Dist: xyzservices>=2023.10.0; extra == "viz"
+Requires-Dist: mapclassify>=2.6.0; extra == "viz"
+Requires-Dist: geoplot>=0.5.0; extra == "viz"
+Requires-Dist: mplleaflet>=0.0.5; extra == "viz"
+Requires-Dist: branca>=0.6.0; extra == "viz"
+Requires-Dist: plotly-geo>=1.0.0; extra == "viz"
+Requires-Dist: streamlit>=1.28.0; extra == "viz"
+Requires-Dist: pyvista; extra == "viz"
+Provides-Extra: viz-3d
+Requires-Dist: keplergl>=0.3.0; extra == "viz-3d"
+Provides-Extra: cloud
+Requires-Dist: boto3>=1.34; extra == "cloud"
+Requires-Dist: pystac>=1.9; extra == "cloud"
+Requires-Dist: planetary-computer>=1.0; extra == "cloud"
+Provides-Extra: notebook
+Requires-Dist: ipython>=8.0; extra == "notebook"
+Requires-Dist: ipywidgets>=8.0; extra == "notebook"
+Requires-Dist: jupyter>=1.0; extra == "notebook"
+Provides-Extra: schedule
+Requires-Dist: croniter>=2.0; extra == "schedule"
+Provides-Extra: dev
+Requires-Dist: pytest>=8.0; extra == "dev"
+Requires-Dist: pytest-cov>=4.1; extra == "dev"
+Requires-Dist: pytest-mock>=3.12; extra == "dev"
+Requires-Dist: ruff>=0.3; extra == "dev"
+Requires-Dist: mypy>=1.8; extra == "dev"
+Requires-Dist: build>=1.0; extra == "dev"
+Requires-Dist: twine>=5.0; extra == "dev"
+Provides-Extra: insar-full
+Requires-Dist: snaphu>=0.1; extra == "insar-full"
+Requires-Dist: scipy>=1.12; extra == "insar-full"
+Provides-Extra: full
+Requires-Dist: pygeofetch[cloud,geo,notebook,processor,providers,sar,schedule,viz]; extra == "full"
+Provides-Extra: all
+Requires-Dist: pygeofetch[full]; extra == "all"
+Dynamic: license-file
+
+<div align="center">
+
+<img src="https://raw.githubusercontent.com/EOCoreINT/pygeofetch/main/icon/pygeofetch_logo.png" alt="PyGeoFetch Logo" width="350">
+
+[![PyPI version](https://badge.fury.io/py/pygeofetch.svg)](https://pypi.org/project/pygeofetch/)
+[![Python Versions](https://img.shields.io/pypi/pyversions/pygeofetch.svg)](https://pypi.org/project/pygeofetch/)
+[![License: MIT](https://img.shields.io/badge/License-MIT-yellow.svg)](https://opensource.org/licenses/MIT)
+[![Tests](https://github.com/appiahkubis14/pygeofetch/actions/workflows/tests.yml/badge.svg)](https://github.com/appiahkubis14/pygeofetch/actions)
+[![Coverage](https://codecov.io/gh/appiahkubis14/pygeofetch/branch/main/graph/badge.svg)](https://codecov.io/gh/appiahkubis14/pygeofetch)
+
+**A unified satellite data and geospatial processing platform — one CLI, one Python API, 24 providers, and a fully verified InSAR chain.**
+
+[Quick Start](#-quick-start) · [Mexico City Case Study](#-case-study-mapping-mexico-city-land-subsidence-with-insar) · [Case Study Repo ↗](https://github.com/EOCoreINT/mexico-subsidence-project) · [Documentation](https://appiahkubis14.github.io/pygeofetch-docs/) · [Notebooks](#-notebooks)
+
+</div>
+
+---
+
+## Overview
+
+pygeofetch is a production-ready framework for acquiring and processing Earth observation data. It provides authenticated, unified access to **24 satellite repositories** — Sentinel, Landsat, Planet, Maxar, Airbus, Copernicus, USGS, NASA, JAXA, and more — through a single CLI and Python API, and layers a complete geospatial processing engine on top.
+
+Where pygeofetch goes further than most data-access libraries is in its **InSAR chain**: burst-aware coregistration, real flat-earth and topographic phase removal, ERA5 tropospheric and ionospheric correction, phase unwrapping, and SBAS time-series inversion — each stage independently verified against synthetic ground truth or real, published deformation studies, not assumed correct from theory alone. The [Mexico City case study](#-case-study-mapping-mexico-city-land-subsidence-with-insar) below walks through that chain end to end on real Sentinel-1 data.
+
+**Core capabilities:**
+
+| | |
+|---|---|
+| 🔐 **Authenticated access** | 24 providers, credentials stored via system keyring (Keychain / Credential Manager / Secret Service) |
+| 🔍 **Federated search** | One query across all providers → STAC 1.0 GeoJSON / GeoParquet / CSV, with real footprint geometry where the provider supplies it |
+| 📥 **Resilient downloads** | Parallel, resumable, checksum-verified, band-selective, atomic writes |
+| ⚙️ **Preprocessing** | Atmospheric correction, cloud masking, reprojection, resampling, pan-sharpening, mosaicking |
+| 📊 **17 spectral indices** | NDVI, EVI, SAVI, NDWI, MNDWI, NDBI, LST, Albedo, dNBR, GLCM texture, and more |
+| 🌐 **Verified InSAR chain** | Coregistration → phase removal → atmospheric/ionospheric correction → unwrapping → SBAS inversion |
+| 📋 **YAML pipelines** | Chainable, schedulable, repeatable workflows with full run history |
+
+---
+
+<p align="center">
+  <img src="https://raw.githubusercontent.com/EOCoreINT/pygeofetch/main/docs/images/trend_map.png" width="48%" />
+  <img src="https://raw.githubusercontent.com/EOCoreINT/pygeofetch/main/docs/images/trend_classification.png" width="48%" />
+</p>
+<p align="center"><em>NDVI trend (2018–2024) and severity classification for the Obuasi Municipal District, Ghana — computed end-to-end with pygeofetch from boundary-clipped USGS Landsat data.</em></p>
+
+Across 2018–2024, 32.0% of the district shows measurable vegetation decline (9.2% strong, 22.8% moderate), against 57.5% stable and 10.5% increasing. The decline is spatially concentrated west of −1.70° longitude, consistent with the district's documented small-scale mining activity — though an NDVI trend alone can't isolate cause from correlated signals like logging or agricultural clearing without ground verification. The eastern two-thirds of the district is overwhelmingly stable, which is itself informative: it suggests a localized driver rather than a district-wide seasonal artifact.
+
+---
+
+## Why pygeofetch
+
+Satellite data access is fragmented — every provider has its own auth scheme, query API, and file format, and most tools that unify search stop well short of processing. pygeofetch is the only package that combines broad provider coverage with a full processing engine and a genuinely verified InSAR chain in one place.
+
+| Feature | pygeofetch | EODAG | pystac-client | satpy | sentinelsat |
+|---|---|---|---|---|---|
+| Providers | **24** | 10+ | STAC only | Limited | Sentinel only |
+| Processing engine | ✅ Full | ❌ | ❌ | Partial | ❌ |
+| Spectral indices | ✅ 17+ | ❌ | ❌ | ❌ | ❌ |
+| Full InSAR chain | ✅ Verified | ❌ | ❌ | ❌ | ❌ |
+| Real footprint geometry | ✅ | ❌ | ✅ (STAC only) | ❌ | ❌ |
+| YAML pipelines + scheduling | ✅ | ❌ | ❌ | ❌ | ❌ |
+| Commercial providers | ✅ Planet/Maxar | ❌ | ❌ | ❌ | ❌ |
+
+---
+
+## Providers
+
+**11 open-access providers**, no login required — including `planetary_computer`, `aws_earth`, `element84`, `noaa_big_data`, `esa_scihub`, `jaxa_earth`, and `geoserver_generic` for any OGC WMS/WFS/WCS endpoint.
+
+**13 authenticated providers** — USGS, Copernicus CDSE, NASA Earthdata (+ Cloud), Alaska SAR Facility, OpenTopography, Planet Labs, Sentinel Hub, Maxar GBDX, Airbus OneAtlas, Google Earth Engine, TerraBotics, and Earth Explorer.
+
+Every provider has been directly verified to return at minimum a correct bounding box; several are further confirmed to return real, precise footprint geometry against a live API response (marked in `pygeofetch providers info PROVIDER`). Run `pygeofetch providers list` for the full, current table.
+
+---
+
+## Installation
+
+```bash
+pip install pygeofetch                 # core — free providers work immediately
+pip install "pygeofetch[geo]"          # + rasterio, geopandas, shapely
+pip install "pygeofetch[insar]"        # + native SBAS inversion, phase unwrapping
+pip install "pygeofetch[insar-full]"   # + MintPy passthrough for advanced corrections
+pip install "pygeofetch[all]"          # everything
+```
+
+Requires Python 3.9+. Run `pygeofetch doctor` to verify your install and provider connectivity.
+
+---
+
+## ⚡ Quick Start
+
+```python
+from pygeofetch import PyGeoFetch
+from pygeofetch.models.search_query import SearchQuery, BoundingBox
+from pygeofetch.models.download_task import DownloadOptions
+
+client = PyGeoFetch()
+client.add_credentials("copernicus", username="you@example.com", password="...")
+
+results = client.search(
+    SearchQuery(bbox=BoundingBox.from_string("-74.1,40.6,-73.7,40.9"),
+                start_date="2024-01-01", cloud_cover_max=20),
+    providers=["copernicus", "planetary_computer", "aws_earth"],
+)
+
+downloads = client.download(results[:5], destination="./data/",
+                             options=DownloadOptions(parallel=4, bands=["B02", "B03", "B04"]))
+
+ndvi = client.indices.ndvi(red="B04.tif", nir="B08.tif")
+```
+
+The equivalent is available as a CLI (`pygeofetch search run ...` / `pygeofetch download run ...`) and as chainable, schedulable YAML pipelines — see the [full CLI reference](#cli-reference) and [notebooks](#-notebooks).
+
+---
+
+## 🌐 Case Study: Mapping Mexico City Land Subsidence with InSAR
+
+**Full project, notebook, and data:** [github.com/EOCoreINT/mexico-subsidence-project](https://github.com/EOCoreINT/mexico-subsidence-project)
+
+Mexico City sits on a drained lakebed and is subsiding at some of the fastest rates on Earth as its aquifer is depleted. It's a demanding, real-world test for an InSAR chain: a huge, largely incoherent urban scene, a fragmented acquisition history, and a deformation signal large enough to make phase unwrapping genuinely hard. What follows is a condensed walkthrough of the [full project above](https://github.com/EOCoreINT/mexico-subsidence-project) — the pygeofetch InSAR chain run end to end against real Sentinel-1 data over Iztapalapa, no synthetic shortcuts. For the complete methodology, every intermediate number, and an unusually candid validation discussion, read the project README itself.
+
+**1. Authenticated federated search.** `client.search()` against Copernicus Data Space for SLC scenes over the city (bbox `-99.183, 19.278, -99.003, 19.438`, July 2016 – September 2017) returns a stack of Sentinel-1A/1B candidates in seconds, filterable straight down to a single consistent track.
+
+<p align="center"><img src="https://raw.githubusercontent.com/EOCoreINT/mexico-subsidence-project/main/analysis/data-search.png" width="85%" /></p>
+
+Results carry real footprint geometry and drop straight onto a map via `MapViewer`, with hover info for scene ID, date, satellite, and provider — useful for spotting swath overlap before committing to a download.
+
+<p align="center">
+  <img src="https://raw.githubusercontent.com/EOCoreINT/mexico-subsidence-project/main/analysis/search_1.png" width="48%" />
+  <img src="https://raw.githubusercontent.com/EOCoreINT/mexico-subsidence-project/main/analysis/search.png" width="48%" />
+</p>
+
+**2. DEM acquisition.** A matching search against OpenTopography returns seven DEM products for the AOI; pygeofetch downloads and clips the selected one (SRTM 30 m) automatically for topographic phase removal.
+
+<p align="center"><img src="https://raw.githubusercontent.com/EOCoreINT/mexico-subsidence-project/main/analysis/dem_acquisition.png" width="85%" /></p>
+
+**3. Preflight validation.** Before anything downloads, the `PreflightGate` checks search truncation, AOI coverage, temporal network connectivity, and burst-timing family compatibility from lightweight annotation XMLs — catching acquisitions that would waste compute before they cost bandwidth. This run avoided roughly 360 GB of downloads by excluding 48 of 115 candidate scenes at this stage alone.
+
+<p align="center"><img src="https://raw.githubusercontent.com/EOCoreINT/mexico-subsidence-project/main/analysis/preflight-log.png" width="85%" /></p>
+
+Real orbit state vectors then pre-filter the candidate pair list against Sentinel-1's own mission thresholds, rejecting pairs for excess burst-sync error, temporal baseline, or spatial baseline before any interferogram is attempted.
+
+<p align="center"><img src="https://raw.githubusercontent.com/EOCoreINT/mexico-subsidence-project/main/analysis/pre-filter-for-interferogram.png" width="85%" /></p>
+
+**4. Interferogram formation.** Of 2,211 theoretically possible pairs across 67 downloaded scenes, real orbit-based coregistration, per-burst Enhanced Spectral Diversity, deburst, and flat-earth/topographic phase removal reduce that to 78 valid pairs — the dense urban fringe pattern below is raw wrapped phase, before any correction.
+
+<p align="center"><img src="https://raw.githubusercontent.com/EOCoreINT/mexico-subsidence-project/main/analysis/wrapped_phase.png" width="85%" /></p>
+
+**5. Coherence.** Amplitude and coherence are estimated together for every pair; the urban core holds noticeably higher coherence (0.614) than the scene-wide average (0.508) — exactly the contrast expected between stable built structures and vegetated or agricultural surrounds.
+
+<p align="center"><img src="https://raw.githubusercontent.com/EOCoreINT/mexico-subsidence-project/main/analysis/comparism.png" width="85%" /></p>
+
+**6. Phase unwrapping.** SNAPHU (Chen & Zebker, 2001) resolves each wrapped interferogram into continuous phase, both as a raw array and reprojected onto the real city footprint. Reliable coverage is intentionally sparse (~0.1–0.2%) — dense vegetation cover means most of the scene decorrelates between 18–24 day passes, so the pipeline marks those pixels NaN rather than interpolating across them.
+
+<p align="center">
+  <img src="https://raw.githubusercontent.com/EOCoreINT/mexico-subsidence-project/main/analysis/unwrapped_phase.png" width="48%" />
+  <img src="https://raw.githubusercontent.com/EOCoreINT/mexico-subsidence-project/main/analysis/unwrap.png" width="48%" />
+</p>
+
+**7. SBAS time-series inversion.** The full 58-date descending stack, once low-coherence and unreliable-reference-pixel pairs are excluded, fractures into 17 disconnected network "islands" — a genuine finding about this AOI's data quality, not a bug. The largest connected island (10 dates, 11 pairs) carries the final inversion, referenced near Cerro de la Estrella.
+
+<p align="center"><img src="https://raw.githubusercontent.com/EOCoreINT/mexico-subsidence-project/main/analysis/1.png" width="85%" /></p>
+
+The resulting velocity field, uncertainty map, reliability mask, and displacement time series recover a 2nd–98th-percentile vertical velocity of **−35.5 to −2.9 cm/year** where reliable pixels exist.
+
+<p align="center"><img src="https://raw.githubusercontent.com/EOCoreINT/mexico-subsidence-project/main/analysis/island1_sbas_full_analysis.png" width="90%" /></p>
+
+That result is **directionally consistent** — same sign, same order of magnitude — with the peer-reviewed peak of −39.1 cm/year reported for Iztapalapa by Cigna & Tapete (2021, *Remote Sensing of Environment*), a study built from 300+ scenes across six years. It is explicitly *not* a validated replication of that study: this run used a fraction of the data volume, and the project's own README is deliberately careful not to overclaim the comparison — worth reading in full if you want the honest version of what this number does and doesn't prove.
+
+The full, runnable notebook — real search and download throughout, no synthetic substitutions — is [`mexico_city_full_confirmed (2).ipynb`](<https://github.com/EOCoreINT/mexico-subsidence-project/blob/main/mexico_city_full_confirmed%20(2).ipynb>) in the [mexico-subsidence-project](https://github.com/EOCoreINT/mexico-subsidence-project) repo; see also `piton_fournaise_full.ipynb` and `amatrice_full.ipynb` in [`notebooks/`](notebooks/) for the same chain applied to the 2021 Piton de la Fournaise eruption and the 2016 Amatrice earthquake.
+
+---
+
+## Processing Engine
+
+**Preprocessing** (`client.preprocess`): `atmos` (DOS1/DOS2, Sen2Cor, FLAASH, 6S, iCOR) · `cloud_mask` · `cloud_fill` · `topo_correct` · `clip` · `reproject` · `resample` · `pansharpen` · `tile` · `mosaic` · `composite`
+
+**Spectral indices** (`client.indices`): `ndvi` · `evi` · `savi` · `ndwi` · `mndwi` · `ndbi` · `ndsi` · `ndmi` · `nbr` / `dnbr` · `tct` · `pca` · `texture` (GLCM) · `lst` · `albedo` · `band_math`
+
+**Post-processing** (`client.post`): `vectorize` → `smooth` → `regularize` → `zonal_stats` → `buffer` → `centroids` → `compress` → `cog`
+
+**SAR** (`client.sar`): `despeckle` (Lee, Enhanced Lee, Frost, Gamma MAP) · `calibrate` (σ0/γ0/β0) · `flood_map` · `coherence`
+
+**InSAR** (`pygeofetch.insar`): `SLCExtractor` · `InterferogramGenerator` · `AtmosphericCorrector` (elevation-correlated or real ERA5) · `IonosphericCorrector` (real IONEX/CDDIS) · `PhaseUnwrapper` (SNAPHU) · `SBASTimeSeries` · `DataValidator` · `InSARProject` (high-level, search-to-interferogram wrapper)
+
+Full method signatures and options are in the [documentation](https://appiahkubis14.github.io/pygeofetch-docs/) and [`pygeofetch/insar/README.md`](pygeofetch/insar/README.md).
+
+---
+
+## YAML Pipelines
+
+```yaml
+name: weekly-sentinel2-ndvi
+schedule: "0 6 * * 1"
+steps:
+  - search: {providers: [copernicus, aws_earth], date_range: last_7_days, cloud_cover: "0-10"}
+  - filter: {expression: "data.cloud_cover < 5"}
+  - download: {parallel: 4, output: ./raw/, bands: [B04, B08]}
+  - ndvi: {red: B04.tif, nir: B08.tif}
+  - vectorize: {threshold: 0.3, format: geojson}
+  - cog: {compress: deflate}
+```
+
+```bash
+pygeofetch proc-pipeline run weekly-sentinel2.yaml --input scene.tif
+pygeofetch pipeline schedule weekly-sentinel2.yaml --name ndvi-weekly
+```
+
+Six built-in templates ship out of the box: `ndvi` · `change_detection` · `flood_map` · `urban_mapping` · `sar_analysis` · `land_cover`.
+
+---
+
+## CLI Reference
+
+Every capability above is also exposed as a CLI command, grouped by area:
+
+| Group | Example | Purpose |
+|---|---|---|
+| `auth` | `pygeofetch auth add copernicus --username U --password P` | Manage per-provider credentials |
+| `providers` | `pygeofetch providers list --capabilities sar` | Discover and inspect providers |
+| `search` | `pygeofetch search run --bbox "..." --providers aws_earth,copernicus` | Federated search, 7 output formats |
+| `download` | `pygeofetch download run --from-search results.geojson --parallel 4` | Parallel, resumable downloads |
+| `preprocess` | `pygeofetch preprocess clip scene.tif --bbox "..."` | Atmospheric correction, clipping, reprojection, etc. |
+| `index` | `pygeofetch index ndvi --red B04.tif --nir B08.tif` | Spectral index computation |
+| `post` | `pygeofetch post cog scene.tif --compress deflate` | Vectorization, zonal stats, COG conversion |
+| `sar` | `pygeofetch sar coherence slc1.tif slc2.tif` | Despeckling, calibration, coherence, flood mapping |
+| `proc-pipeline` / `pipeline` | `pygeofetch pipeline schedule weekly.yaml` | YAML pipeline execution and scheduling |
+| `cache` / `config` | `pygeofetch cache prune --max-size 1GB` | Cache and configuration management |
+
+Run `pygeofetch --help` or any subcommand with `--help` for the complete, current option list, or see the [CLI reference docs](https://appiahkubis14.github.io/pygeofetch-docs/).
+
+---
+
+<!-- ## 📚 Notebooks
+
+| Notebook | Topics |
+|---|---|
+| `01_getting_started.ipynb` | Install, doctor, first search, first download |
+| `02_authentication_and_providers.ipynb` | All 24 providers, credentials, capability filters |
+| `03_advanced_search.ipynb` | Federated search, CQL2 filters, output formats, caching |
+| `04_download_and_postprocessing.ipynb` | Band selection, parallel downloads, post-processing |
+| `05_pipelines_and_scheduling.ipynb` | YAML pipelines, scheduling, Python builder API |
+| `06_real_world_workflows.ipynb` | NDVI time series, change detection, multi-sensor fusion |
+| `07_copernicus_and_authenticated_providers.ipynb` | Copernicus, USGS, NASA, Planet, ASF, OpenTopography |
+| `08_cli_complete_reference.ipynb` | Every CLI command with runnable examples |
+| `09_processing_complete.ipynb` | Full processing engine: preprocessing, indices, SAR, pipelines |
+| [`mexico_city_full_confirmed (2).ipynb`](<https://github.com/EOCoreINT/mexico-subsidence-project/blob/main/mexico_city_full_confirmed%20(2).ipynb>) ↗ | The [Mexico City case study](#-case-study-mapping-mexico-city-land-subsidence-with-insar) above, in a separate companion repo — real search/download, full InSAR chain, SBAS inversion |
+| `piton_fournaise_full.ipynb` | Verified InSAR pipeline against the real April 2021 Piton de la Fournaise eruption |
+| `amatrice_full.ipynb` | Verified InSAR pipeline against the real 2016 Amatrice, Italy earthquake |
+| `provider_search_footprint_test.ipynb` | Search and footprint-map validation across all 24 providers |
+
+```bash
+cd notebooks/ && jupyter lab
+``` -->
+
+---
+
+## Documentation
+
+Full documentation: **https://appiahkubis14.github.io/pygeofetch-docs/** — CLI reference, provider auth guides, pipeline configuration, post-processing catalogue, and InSAR verification methodology.
+
+## Contributing
+
+Contributions are welcome — see [CONTRIBUTING.md](CONTRIBUTING.md). Good first issues include extending real footprint geometry to remaining bbox-only providers, wiring stub providers to full API integrations, and adding new post-processing actions.
+
+```bash
+git clone git@github.com:EOCoreINT/pygeofetch.git && cd pygeofetch
+pip install -e ".[dev,all]" && pytest tests/unit/ -v
+```
+
+## License
+
+MIT — see [LICENSE](LICENSE). © 2026 Samuel Appiah Kubi. Part of the **PyGeoVision** platform, alongside [pygeofetch](https://github.com/appiahkubis14/pygeofetch) itself.

{pygeofetch-2.5.0.dist-info → pygeofetch-2.6.1.dist-info}/RECORD

@@ -1,5 +1,5 @@
-pygeofetch/__init__.py,sha256=GWQJPhEHRE912aufB3A3thLi87N5Zj_ZjsPVa08UQUA,1990
-pygeofetch/meta.yaml,sha256=5hj9c3zPaZ-aRyPnb1xpN-2jgrw7L_AB2qW3hfoLURg,1280
+pygeofetch/__init__.py,sha256=ZCRJRFP3DKU1xf4yo6eQHjHKT4nUjWjWohebYE1Yc7M,1990
+pygeofetch/meta.yaml,sha256=zSwzpicrlpHWUze00OmFGFoNRHRGkLaAY7WYMXsancA,1280
 pygeofetch/cli/__init__.py,sha256=QZMIIHEgmjGFieucqq9fZ5uephxP9-jdu1s8Q4qqw8o,30
 pygeofetch/cli/auth_commands.py,sha256=YqbrN_BYcAYCaUIYLpSJ_VIw35tKFUiD9Kc55SWbDWg,7742
 pygeofetch/cli/config_commands.py,sha256=1Pps0HLCLPODMlM6SFV8ru3IKCDpexSTz8uKIXvIIF8,3918
@@ -36,7 +36,7 @@ pygeofetch/insar/geolocation.py,sha256=-vwTYFhvutKjI4c_f-3BGkfV3XnnDUh5r7KwQbLlE
 pygeofetch/insar/gpu.py,sha256=E4Dx9-78s2GyKl2_dJrc6lXIHiy67tlqHlpkYue_418,3933
 pygeofetch/insar/interferogram.py,sha256=6T9k7LAGupWcib2S26fyacNwKi3DBP2jAKbTM8M2Bvw,111471
 pygeofetch/insar/ionosphere.py,sha256=vG60pyUtxyoJqLCMb_KjFEJvkDE1Bl8fqTzTTbrTIX8,30101
-pygeofetch/insar/preflight.py,sha256=s5MtHVNocHhWtQqudpcrhCvNCBUubfuB68Hd3EJlcQ0,39070
+pygeofetch/insar/preflight.py,sha256=EDrSxzlPIrBHtYgcMyaoaBr7CmInvg6FlHRb9k6EJfw,41270
 pygeofetch/insar/provenance.py,sha256=b2Ql40ZVeI9xHwFUvstBl4npK2o4lua39sRcpZI_lAQ,1876
 pygeofetch/insar/splitspectrum_diagnostic.py,sha256=_IlfBs__BW7Yz2RwdmWxfrOspe67jHMgJPyNyzK_rYE,7705
 pygeofetch/insar/stack_selection.py,sha256=FaRBXDNnLwwiXPGeALSW1AaUdpLFOEepVTgEGVfNC4o,37555
@@ -108,9 +108,9 @@ pygeofetch/utils/validators.py,sha256=f-k3mGFP4byjimyDymhPEdsz5I16jZacx1HSnwrARj
 pygeofetch/viz/__init__.py,sha256=xf2C9IL4n3bv3oxk8RYkogq7nkuNyt18V-urAoXWRQY,449
 pygeofetch/viz/map.py,sha256=vYO99NvTumCAVXVtXcEzMx_cpmH-FvSVNWl8X0PNSgE,15494
 pygeofetch/viz/plot.py,sha256=I16rdOWb80gcjYDsDvTMzuyek7QHgUNOQxY2RejBAPs,53675
-pygeofetch-2.5.0.dist-info/licenses/LICENSE,sha256=zavMyhODD6NvZTmmKFEtATrngwY4CU_NEtY7lhjxGnI,1075
-pygeofetch-2.5.0.dist-info/METADATA,sha256=frDzX4t6B0Vc_j-ggBrzZUHzgJS9KEsjQCZ7OHD0YBQ,43254
-pygeofetch-2.5.0.dist-info/WHEEL,sha256=YVMoNqKzERt-wjUZwJ33xBGAwnFl-4cqbYkTtWa4itE,91
-pygeofetch-2.5.0.dist-info/entry_points.txt,sha256=dHiUXLF5xNJ8vgtXZoUQNleEMQUeDPZiKQOE1tW8SDY,55
-pygeofetch-2.5.0.dist-info/top_level.txt,sha256=rOmnZargIfNk99fwVOqt6OJyJvpSTBBAYO5HU37yhHk,11
-pygeofetch-2.5.0.dist-info/RECORD,,
+pygeofetch-2.6.1.dist-info/licenses/LICENSE,sha256=zavMyhODD6NvZTmmKFEtATrngwY4CU_NEtY7lhjxGnI,1075
+pygeofetch-2.6.1.dist-info/METADATA,sha256=Lp15FqyJcP359D5zenSKnO98FmOLT92Rvk2RffGoUtw,25946
+pygeofetch-2.6.1.dist-info/WHEEL,sha256=YVMoNqKzERt-wjUZwJ33xBGAwnFl-4cqbYkTtWa4itE,91
+pygeofetch-2.6.1.dist-info/entry_points.txt,sha256=dHiUXLF5xNJ8vgtXZoUQNleEMQUeDPZiKQOE1tW8SDY,55
+pygeofetch-2.6.1.dist-info/top_level.txt,sha256=rOmnZargIfNk99fwVOqt6OJyJvpSTBBAYO5HU37yhHk,11
+pygeofetch-2.6.1.dist-info/RECORD,,