PyGeoFetch 1.8.0__tar.gz → 1.9.1__tar.gz

{pygeofetch-1.8.0 → pygeofetch-1.9.1}/CHANGELOG.md

@@ -1,9 +1,9 @@
 # Changelog
 
 All notable changes to PyGeoFetch are documented here.
-Format follows [Keep a Changelog](https://keepachangelog.com/en/1.8.0/).
+Format follows [Keep a Changelog](https://keepachangelog.com/en/1.9.1/).
 
-## [1.8.0] — 2026-07-12
+## [1.9.1] — 2026-07-12
 
 ### Added
 - Federated search across 22+ satellite data providers

{pygeofetch-1.8.0/pygeofetch.egg-info → pygeofetch-1.9.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pygeofetch
-Version: 1.8.0
+Version: 1.9.1
 Summary: Universal satellite data pipeline + geospatial processing platform
 Author-email: Samuel Appiah Kubi <appiahkubis14@gmail.com>
 License: MIT License
@@ -94,8 +94,8 @@ Requires-Dist: plotly>=5.14.0; extra == "viz"
 Requires-Dist: whitebox>=2.3.0; extra == "viz"
 Requires-Dist: whiteboxgui>=2.3.0; extra == "viz"
 Requires-Dist: geemap>=0.28.0; extra == "viz"
-Requires-Dist: osmnx>=1.8.0; extra == "viz"
-Requires-Dist: contextily>=1.8.0; extra == "viz"
+Requires-Dist: osmnx>=1.9.1; extra == "viz"
+Requires-Dist: contextily>=1.9.1; extra == "viz"
 Requires-Dist: xyzservices>=2023.10.0; extra == "viz"
 Requires-Dist: mapclassify>=2.5.0; extra == "viz"
 Requires-Dist: geoplot>=0.5.0; extra == "viz"
@@ -467,6 +467,70 @@ pygeofetch proc-pipeline template change_detection
 
 ---
 
+
+### 🌐 InSAR Processing (`pygeofetch.insar`)
+
+A full Interferometric SAR chain for Sentinel-1 — search through SBAS time series inversion, pure Python, no external InSAR software required beyond `snaphu-py`'s bundled unwrapper.
+
+```bash
+pip install "pygeofetch[insar]"        # native SBAS inversion
+pip install "pygeofetch[insar-full]"   # + MintPy passthrough for advanced corrections
+```
+
+| Component | Description |
+|---|---|
+| `InterferogramGenerator` | Coregistration (real orbit-based, or shape-based fallback) + Enhanced Spectral Diversity (ESD) refinement for TOPS burst continuity + topographic phase removal |
+| `PhaseUnwrapper` | SNAPHU (Chen & Zebker 2001) via the official `snaphu-py` bindings — the same algorithm used by ASF, ISCE2/3, GAMMA, and SNAP |
+| `AtmosphericCorrector` | Elevation-correlated regression or ERA5 reanalysis-based tropospheric delay correction |
+| `SBASTimeSeries` | Small BAseline Subset displacement/velocity inversion (Berardino et al. 2002) |
+| `DataValidator` | SLC, coherence, and SBAS-network sanity checks, wired in at every real pipeline entry point |
+
+Real orbit-based coregistration computes genuine per-pixel offsets from actual satellite orbit state vectors and acquisition timing (not a shape-matching guess) — supply a DEM, both SAFE archives, and both `.EOF` orbit files, and it's used automatically; otherwise falls back cleanly to shape-based resampling.
+
+Every stage supports automatic visualization (`auto_visualize=True`) and optional GPU acceleration (`use_gpu=True`, CuPy) for coherence estimation and SBAS inversion on large scenes.
+
+See [`pygeofetch/insar/README.md`](pygeofetch/insar/README.md) for the full processing chain, verification methodology, and current limitations.
+
+
+```python
+# InSAR — search to SBAS time series
+from pygeofetch.insar import (
+    InterferogramGenerator, PhaseUnwrapper, AtmosphericCorrector,
+    SBASTimeSeries, DataValidator,
+)
+from pygeofetch.insar.timeseries import InterferogramPair
+from pygeofetch.core.orbits import fetch_orbit_file
+
+# Coregistration + interferogram formation (real orbit-based when a DEM,
+# both SAFE archives, and both orbit files are supplied; falls back to
+# shape-based resampling otherwise)
+gen = InterferogramGenerator(coherence_window=5, esd_enabled=True, use_gpu=False)
+result = gen.process_pair(
+    reference="slc_ref.tif", secondary="slc_sec.tif", dem="dem.tif",
+    reference_safe_zip="ref.SAFE.zip", secondary_safe_zip="sec.SAFE.zip",
+    reference_orbit_file=fetch_orbit_file("S1A_..._ref"),
+    secondary_orbit_file=fetch_orbit_file("S1A_..._sec"),
+)
+result.save("./output", auto_visualize=True)
+
+# Atmospheric correction + phase unwrapping
+atm = AtmosphericCorrector(method="elevation")
+corrected_phase = atm.correct(result.interferogram, dem="dem.tif")
+
+unwrapper = PhaseUnwrapper(cost_mode="defo", init_method="mcf")
+unwrapped, conncomp = unwrapper.unwrap(corrected_phase, result.coherence)
+
+# SBAS time series inversion across a network of interferograms
+pairs = [
+    InterferogramPair("2026-01-01", "2026-01-13", unwrapped, result.coherence),
+    # ... additional pairs forming a connected network
+]
+sbas = SBASTimeSeries(wavelength_m=0.05546576, reference_date="2026-01-01")
+ts_result = sbas.invert(pairs, reference_pixel=(10, 15))  # known-stable pixel
+ts_result.save("./timeseries", auto_visualize=True)
+print(f"Mean velocity: {ts_result.velocity.mean()*1000:.1f} mm/year")
+```
+
 ## 🖥️ Complete CLI Reference
 
 ```
@@ -653,3 +717,8 @@ pytest tests/unit/ -v
 pygeofetch is free and open source software, licensed under the [MIT License](LICENSE).
 
 © 2026 Samuel Appiah Kubi. Part of the **PyGeoVision** platform — [pygeofetch](https://github.com/appiahkubis14/pygeofetch) (data + processing)  complete Earth observation pipeline.
+
+
+
+
+

{pygeofetch-1.8.0 → pygeofetch-1.9.1}/README.md

@@ -331,6 +331,70 @@ pygeofetch proc-pipeline template change_detection
 
 ---
 
+
+### 🌐 InSAR Processing (`pygeofetch.insar`)
+
+A full Interferometric SAR chain for Sentinel-1 — search through SBAS time series inversion, pure Python, no external InSAR software required beyond `snaphu-py`'s bundled unwrapper.
+
+```bash
+pip install "pygeofetch[insar]"        # native SBAS inversion
+pip install "pygeofetch[insar-full]"   # + MintPy passthrough for advanced corrections
+```
+
+| Component | Description |
+|---|---|
+| `InterferogramGenerator` | Coregistration (real orbit-based, or shape-based fallback) + Enhanced Spectral Diversity (ESD) refinement for TOPS burst continuity + topographic phase removal |
+| `PhaseUnwrapper` | SNAPHU (Chen & Zebker 2001) via the official `snaphu-py` bindings — the same algorithm used by ASF, ISCE2/3, GAMMA, and SNAP |
+| `AtmosphericCorrector` | Elevation-correlated regression or ERA5 reanalysis-based tropospheric delay correction |
+| `SBASTimeSeries` | Small BAseline Subset displacement/velocity inversion (Berardino et al. 2002) |
+| `DataValidator` | SLC, coherence, and SBAS-network sanity checks, wired in at every real pipeline entry point |
+
+Real orbit-based coregistration computes genuine per-pixel offsets from actual satellite orbit state vectors and acquisition timing (not a shape-matching guess) — supply a DEM, both SAFE archives, and both `.EOF` orbit files, and it's used automatically; otherwise falls back cleanly to shape-based resampling.
+
+Every stage supports automatic visualization (`auto_visualize=True`) and optional GPU acceleration (`use_gpu=True`, CuPy) for coherence estimation and SBAS inversion on large scenes.
+
+See [`pygeofetch/insar/README.md`](pygeofetch/insar/README.md) for the full processing chain, verification methodology, and current limitations.
+
+
+```python
+# InSAR — search to SBAS time series
+from pygeofetch.insar import (
+    InterferogramGenerator, PhaseUnwrapper, AtmosphericCorrector,
+    SBASTimeSeries, DataValidator,
+)
+from pygeofetch.insar.timeseries import InterferogramPair
+from pygeofetch.core.orbits import fetch_orbit_file
+
+# Coregistration + interferogram formation (real orbit-based when a DEM,
+# both SAFE archives, and both orbit files are supplied; falls back to
+# shape-based resampling otherwise)
+gen = InterferogramGenerator(coherence_window=5, esd_enabled=True, use_gpu=False)
+result = gen.process_pair(
+    reference="slc_ref.tif", secondary="slc_sec.tif", dem="dem.tif",
+    reference_safe_zip="ref.SAFE.zip", secondary_safe_zip="sec.SAFE.zip",
+    reference_orbit_file=fetch_orbit_file("S1A_..._ref"),
+    secondary_orbit_file=fetch_orbit_file("S1A_..._sec"),
+)
+result.save("./output", auto_visualize=True)
+
+# Atmospheric correction + phase unwrapping
+atm = AtmosphericCorrector(method="elevation")
+corrected_phase = atm.correct(result.interferogram, dem="dem.tif")
+
+unwrapper = PhaseUnwrapper(cost_mode="defo", init_method="mcf")
+unwrapped, conncomp = unwrapper.unwrap(corrected_phase, result.coherence)
+
+# SBAS time series inversion across a network of interferograms
+pairs = [
+    InterferogramPair("2026-01-01", "2026-01-13", unwrapped, result.coherence),
+    # ... additional pairs forming a connected network
+]
+sbas = SBASTimeSeries(wavelength_m=0.05546576, reference_date="2026-01-01")
+ts_result = sbas.invert(pairs, reference_pixel=(10, 15))  # known-stable pixel
+ts_result.save("./timeseries", auto_visualize=True)
+print(f"Mean velocity: {ts_result.velocity.mean()*1000:.1f} mm/year")
+```
+
 ## 🖥️ Complete CLI Reference
 
 ```
@@ -517,3 +581,8 @@ pytest tests/unit/ -v
 pygeofetch is free and open source software, licensed under the [MIT License](LICENSE).
 
 © 2026 Samuel Appiah Kubi. Part of the **PyGeoVision** platform — [pygeofetch](https://github.com/appiahkubis14/pygeofetch) (data + processing)  complete Earth observation pipeline.
+
+
+
+
+

{pygeofetch-1.8.0 → pygeofetch-1.9.1}/pygeofetch/__init__.py

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

{pygeofetch-1.8.0 → pygeofetch-1.9.1}/pygeofetch/core/engine.py

@@ -176,7 +176,7 @@ class PyGeoFetch:
         if _product.upper() == "SLC" and _effective:
             _effective = self._route_slc_providers(_effective)
         results = self.searcher.search(query, providers=_effective, use_cache=use_cache)
-        self._warn_if_outdated_constellation(results)
+        self._warn_if_outdated_constellation(results, query)
         return results
 
     def search_and_save(
@@ -217,19 +217,43 @@ class PyGeoFetch:
     }
     _GRD_ONLY = {"planetary_computer", "aws_earth", "element84"}
 
-    def _warn_if_outdated_constellation(self, results: list) -> None:
-        """Log a warning if results contain only decommissioned satellites after S1A decomm."""
-        from datetime import date
+    def _warn_if_outdated_constellation(self, results: list, query) -> None:
+        """Log a warning if results contain only decommissioned satellites
+        for a search whose OWN requested date range extends past S1A's
+        real decommissioning date — not just because today happens to be
+        past it. A search for genuine historical data (e.g. November
+        2024, well before S1A's actual 1 July 2026 end of mission) must
+        never trigger this, regardless of what today's date is; S1A
+        results from before decommissioning are completely legitimate."""
+        from datetime import date, datetime
 
         S1A_DECOM = date(2026, 7, 1)
-        if date.today() >= S1A_DECOM:
+
+        query_end = getattr(query, "end_date", None)
+        if query_end is None:
+            # No end date specified -- fall back to whether the search
+            # was actually run after decommissioning (the old behaviour,
+            # but only as a last resort when there's no query date to
+            # check against).
+            query_end_date = date.today()
+        elif isinstance(query_end, str):
+            query_end_date = datetime.fromisoformat(query_end.replace("Z", "")).date()
+        elif isinstance(query_end, datetime):
+            query_end_date = query_end.date()
+        elif isinstance(query_end, date):
+            query_end_date = query_end
+        else:
+            query_end_date = date.today()
+
+        if query_end_date >= S1A_DECOM:
             platforms = {getattr(r, "satellite", "") or "" for r in results}
             normalised = {_normalise_satellite_name(p) for p in platforms}
             if any("S1" in p for p in normalised):  # only check if S1 results
                 if not any(p in {"S1C", "S1D"} for p in normalised):
                     logger.warning(
-                        "Search returned only S1A/S1B results after Sentinel-1A "
-                        "decommissioning date. This may indicate a provider collection "
+                        "Search returned only S1A/S1B results for a date range "
+                        "extending past Sentinel-1A's real decommissioning date "
+                        "(1 July 2026). This may indicate a provider collection "
                         "routing issue. Check platform filter settings."
                     )
 
@@ -460,4 +484,4 @@ class PyGeoFetch:
 
     def __repr__(self) -> str:
         authed = [item["provider"] for item in self.auth.list()]
-        return f"PyGeoFetch(authenticated={authed}, cache_entries={len(self.searcher.cache._cache)})"
+        return f"PyGeoFetch(authenticated={authed}, cache_entries={len(self.searcher.cache._cache)})"

{pygeofetch-1.8.0 → pygeofetch-1.9.1}/pygeofetch/meta.yaml

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

{pygeofetch-1.8.0 → pygeofetch-1.9.1/pygeofetch.egg-info}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: pygeofetch
-Version: 1.8.0
+Version: 1.9.1
 Summary: Universal satellite data pipeline + geospatial processing platform
 Author-email: Samuel Appiah Kubi <appiahkubis14@gmail.com>
 License: MIT License
@@ -94,8 +94,8 @@ Requires-Dist: plotly>=5.14.0; extra == "viz"
 Requires-Dist: whitebox>=2.3.0; extra == "viz"
 Requires-Dist: whiteboxgui>=2.3.0; extra == "viz"
 Requires-Dist: geemap>=0.28.0; extra == "viz"
-Requires-Dist: osmnx>=1.8.0; extra == "viz"
-Requires-Dist: contextily>=1.8.0; extra == "viz"
+Requires-Dist: osmnx>=1.9.1; extra == "viz"
+Requires-Dist: contextily>=1.9.1; extra == "viz"
 Requires-Dist: xyzservices>=2023.10.0; extra == "viz"
 Requires-Dist: mapclassify>=2.5.0; extra == "viz"
 Requires-Dist: geoplot>=0.5.0; extra == "viz"
@@ -467,6 +467,70 @@ pygeofetch proc-pipeline template change_detection
 
 ---
 
+
+### 🌐 InSAR Processing (`pygeofetch.insar`)
+
+A full Interferometric SAR chain for Sentinel-1 — search through SBAS time series inversion, pure Python, no external InSAR software required beyond `snaphu-py`'s bundled unwrapper.
+
+```bash
+pip install "pygeofetch[insar]"        # native SBAS inversion
+pip install "pygeofetch[insar-full]"   # + MintPy passthrough for advanced corrections
+```
+
+| Component | Description |
+|---|---|
+| `InterferogramGenerator` | Coregistration (real orbit-based, or shape-based fallback) + Enhanced Spectral Diversity (ESD) refinement for TOPS burst continuity + topographic phase removal |
+| `PhaseUnwrapper` | SNAPHU (Chen & Zebker 2001) via the official `snaphu-py` bindings — the same algorithm used by ASF, ISCE2/3, GAMMA, and SNAP |
+| `AtmosphericCorrector` | Elevation-correlated regression or ERA5 reanalysis-based tropospheric delay correction |
+| `SBASTimeSeries` | Small BAseline Subset displacement/velocity inversion (Berardino et al. 2002) |
+| `DataValidator` | SLC, coherence, and SBAS-network sanity checks, wired in at every real pipeline entry point |
+
+Real orbit-based coregistration computes genuine per-pixel offsets from actual satellite orbit state vectors and acquisition timing (not a shape-matching guess) — supply a DEM, both SAFE archives, and both `.EOF` orbit files, and it's used automatically; otherwise falls back cleanly to shape-based resampling.
+
+Every stage supports automatic visualization (`auto_visualize=True`) and optional GPU acceleration (`use_gpu=True`, CuPy) for coherence estimation and SBAS inversion on large scenes.
+
+See [`pygeofetch/insar/README.md`](pygeofetch/insar/README.md) for the full processing chain, verification methodology, and current limitations.
+
+
+```python
+# InSAR — search to SBAS time series
+from pygeofetch.insar import (
+    InterferogramGenerator, PhaseUnwrapper, AtmosphericCorrector,
+    SBASTimeSeries, DataValidator,
+)
+from pygeofetch.insar.timeseries import InterferogramPair
+from pygeofetch.core.orbits import fetch_orbit_file
+
+# Coregistration + interferogram formation (real orbit-based when a DEM,
+# both SAFE archives, and both orbit files are supplied; falls back to
+# shape-based resampling otherwise)
+gen = InterferogramGenerator(coherence_window=5, esd_enabled=True, use_gpu=False)
+result = gen.process_pair(
+    reference="slc_ref.tif", secondary="slc_sec.tif", dem="dem.tif",
+    reference_safe_zip="ref.SAFE.zip", secondary_safe_zip="sec.SAFE.zip",
+    reference_orbit_file=fetch_orbit_file("S1A_..._ref"),
+    secondary_orbit_file=fetch_orbit_file("S1A_..._sec"),
+)
+result.save("./output", auto_visualize=True)
+
+# Atmospheric correction + phase unwrapping
+atm = AtmosphericCorrector(method="elevation")
+corrected_phase = atm.correct(result.interferogram, dem="dem.tif")
+
+unwrapper = PhaseUnwrapper(cost_mode="defo", init_method="mcf")
+unwrapped, conncomp = unwrapper.unwrap(corrected_phase, result.coherence)
+
+# SBAS time series inversion across a network of interferograms
+pairs = [
+    InterferogramPair("2026-01-01", "2026-01-13", unwrapped, result.coherence),
+    # ... additional pairs forming a connected network
+]
+sbas = SBASTimeSeries(wavelength_m=0.05546576, reference_date="2026-01-01")
+ts_result = sbas.invert(pairs, reference_pixel=(10, 15))  # known-stable pixel
+ts_result.save("./timeseries", auto_visualize=True)
+print(f"Mean velocity: {ts_result.velocity.mean()*1000:.1f} mm/year")
+```
+
 ## 🖥️ Complete CLI Reference
 
 ```
@@ -653,3 +717,8 @@ pytest tests/unit/ -v
 pygeofetch is free and open source software, licensed under the [MIT License](LICENSE).
 
 © 2026 Samuel Appiah Kubi. Part of the **PyGeoVision** platform — [pygeofetch](https://github.com/appiahkubis14/pygeofetch) (data + processing)  complete Earth observation pipeline.
+
+
+
+
+

{pygeofetch-1.8.0 → pygeofetch-1.9.1}/pygeofetch.egg-info/requires.txt

@@ -90,8 +90,8 @@ plotly>=5.14.0
 whitebox>=2.3.0
 whiteboxgui>=2.3.0
 geemap>=0.28.0
-osmnx>=1.8.0
-contextily>=1.8.0
+osmnx>=1.9.1
+contextily>=1.9.1
 xyzservices>=2023.10.0
 mapclassify>=2.5.0
 geoplot>=0.5.0

{pygeofetch-1.8.0 → pygeofetch-1.9.1}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "pygeofetch"
-version = "1.8.0"
+version = "1.9.1"
 description = "Universal satellite data pipeline + geospatial processing platform"
 readme = "README.md"
 license = { file = "LICENSE" }
@@ -97,8 +97,8 @@ viz = [
     "whitebox>=2.3.0",
     "whiteboxgui>=2.3.0",
     "geemap>=0.28.0",
-    "osmnx>=1.8.0",
-    "contextily>=1.8.0",
+    "osmnx>=1.9.1",
+    "contextily>=1.9.1",
     "xyzservices>=2023.10.0",
     "mapclassify>=2.5.0",
     "geoplot>=0.5.0",
@@ -143,7 +143,7 @@ dev = [
 insar-full = [
     "snaphu>=0.1",
     "scipy>=1.12",
-    "mintpy>=1.5",
+    "mintpy>=1.5",      # <- almost certainly this one
     "pyaps3>=0.3",
 ]